A molecule resembles a dancing pair: electrons whirl quickly, while heavy nuclei move more slowly. But sometimes the partners lose the rhythm — an electron jumps between states, and the dance changes. This is called a nonadiabatic transition. Usually, scientists describe such scenes either with a mixed quantum-classical method (Ehrenfest dynamics) or with a semiclassical approach using wave packets. The former is good for electronic jumps but crude for nuclear quantum effects. The latter is accurate for nuclei but blind to electronic transitions.
The new method — TGED (Ehrenfest dynamics with thawed Gaussian packets)—is like teaching dancers to listen to each other. It unifies both descriptions in a single equation: the electronic part remains quantum, while the nuclear part is represented by a wave packet that can change shape. The authors derived the method from a variational principle, where the total energy of the system is minimized for the correct motion. This yields a family of variants: you can choose how accurately to account for the local curvature of the molecular potential.
In the classical limit (when nuclei are very heavy), TGED reduces to ordinary Ehrenfest dynamics. If electronic transitions are absent, the method reduces to thawed Gaussian wave packet dynamics. The authors also constructed geometric integrators — numerical schemes that preserve important physical properties (such as energy) during computer simulations. This is like a choreographer who ensures the dance does not fall apart at every step.
as written by the authors arXiv:2607.10847
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.
Ehrenfest dynamics is named after Paul Ehrenfest, who in 1927 proposed a theorem linking quantum and classical equations of motion. His work became a bridge between two worlds of physics.