Exponential Integrators for Quantum-Classical Molecular Dynamics

We study time integration methods for equations of mixed quantum-classical molecular dynamics in which Newtonian equations of motion and Schrödinger equations are nonlinearly coupled. Such systems exhibit different time scales in the classical and the quantum evolution, and the solutions are typically highly oscillatory. The numerical methods use the exponential of the quantum Hamiltonian whose product with a state vector is approximated using Lanczos' method. This allows time steps that are much larger than the inverse of the highest frequencies.