High-order geometric integrators for representation-free Ehrenfest dynamics
arXiv:2107.00607 · doi:10.1063/5.0061878
Abstract
Ehrenfest dynamics is a useful approximation for ab initio mixed quantum-classical molecular dynamics that can treat electronically nonadiabatic effects. Although a severe approximation to the exact solution of the molecular time-dependent Schrödinger equation, Ehrenfest dynamics is symplectic, time-reversible, and conserves exactly the total molecular energy as well as the norm of the electronic wavefunction. Here, we surpass apparent complications due to the coupling of classical nuclear and quantum electronic motions and present efficient geometric integrators for "representation-free" Ehrenfest dynamics, which do not rely on a diabatic or adiabatic representation of electronic states and are of arbitrary even orders of accuracy in the time step. These numerical integrators, obtained by symmetrically composing the second-order splitting method and exactly solving the kinetic and potential propagation steps, are norm-conserving, symplectic, and time-reversible regardless of the time step used. Using a nonadiabatic simulation in the region of a conical intersection as an example, we demonstrate that these integrators preserve the geometric properties exactly and, if highly accurate solutions are desired, can be even more efficient than the most popular non-geometric integrators.
Replaced "Ehrenfest dynamics" with "representation-free Ehrenfest dynamics" in the title. Moved most of the supplementary material to the appendix of the main text. Added new references and a new Eq. (54)
References in corpus (8)
- Generalized spin mapping for quantum-classical dynamics
- Exact quantum statistics for electronically nonadiabatic systems using continuous path variables
- Accurate nonadiabatic quantum dynamics on the cheap: making the most of mean field theory with master equations
- Nonadiabatic semiclassical dynamics in the mixed quantum-classical initial value representation
- Ab initio simulations of excited carrier dynamics in carbon nanotubes
- Why do mixed quantum-classical methods describe short-time dynamics through conical intersections so well? Analysis of geometric phase effects
- Mixed quantum-classical dynamics using collective electronic variables: A better alternative to electronic friction theories
- Evaluation of the importance of spin-orbit couplings in the nonadiabatic quantum dynamics with quantum fidelity and with its efficient "on-the-fly" ab initio semiclassical approximation