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.