Track papers across labs, topics, and trends.
Search papers, labs, and topics across Lattice.
This paper introduces thawed Gaussian Ehrenfest dynamics (TGED), a novel framework that unifies Ehrenfest dynamics and thawed Gaussian wavepacket dynamics to effectively model both electronic nonadiabaticity and nuclear quantum effects. By applying the time-dependent variational principle to a Hartree product of electronic and Gaussian nuclear wavepackets, the authors derive a fully variational formulation that can be adapted through various effective locally quadratic potentials. The study not only establishes TGED's connection to existing methods but also provides geometric integrators and identifies conditions for exact approximations, enhancing the computational efficiency of molecular dynamics simulations.
A unified framework that captures both electronic nonadiabaticity and nuclear quantum effects could revolutionize molecular dynamics simulations.
Ehrenfest dynamics is a widely used mixed quantum--classical approach for nonadiabatic molecular dynamics, whereas thawed Gaussian wavepacket dynamics provides an efficient semiclassical description of adiabatic nuclear quantum dynamics. Here we describe thawed Gaussian Ehrenfest dynamics (TGED), which unifies and generalizes these two methods to capture both electronic nonadiabaticity and nuclear quantum effects within a single framework. The fully variational formulation of TGED is derived by applying the time-dependent variational principle to a Hartree product of electronic and Gaussian nuclear wavepackets. Replacing the effective locally quadratic molecular potential obtained from this variational treatment by alternative effective locally quadratic potentials yields an infinite family of TGED methods, of which we present several members. We analyze the limiting cases of the general formalism and show, in particular, that it reduces to conventional Ehrenfest dynamics in the classical limit for the nuclei and to thawed Gaussian wavepacket dynamics in the absence of electronic coupling. Finally, we present explicit geometric integrators for the entire family of methods and identify the conditions under which the different approximations become exact.