Structure-Preserving Integration for Magnetic Gaussian Wave Packet Dynamics
This work addresses the need for accurate long-time simulations in quantum dynamics for applications like materials science or chemistry, though it is incremental as it adapts existing integration techniques to a specific variational framework.
The authors tackled the problem of simulating Gaussian wave packet dynamics under magnetic fields by developing structure-preserving time integration schemes, resulting in methods that conserve key invariants and show near-conservation of the Hamiltonian with rigorous error bounds uniform in the semiclassical parameter.
We develop structure-preserving time integration schemes for Gaussian wave packet dynamics associated with the magnetic Schrödinger equation. The variational Dirac--Frenkel formulation yields a finite-dimensional Hamiltonian system for the wave packet parameters, where the presence of a magnetic vector potential leads to a non-separable structure and a modified symplectic geometry. By introducing kinetic momenta through a minimal substitution, we reformulate the averaged dynamics as a Poisson system that closely parallels the classical equations of charged particle motion. This representation enables the construction of Boris-type integrators adapted to the variational setting. In addition, we propose explicit high-order symplectic schemes based on splitting methods and partitioned Runge--Kutta integrators. The proposed methods conserve the quadratic invariants characterizing the Hagedorn parametrization, preserve linear and angular momentum under symmetry assumptions, and exhibit near-conservation of the averaged Hamiltonian over long time intervals. Rigorous error estimates are derived for both the wave packet parameters and observable quantities, with bounds uniform in the semiclassical parameter. Numerical experiments demonstrate the favorable long-time behavior and structure preservation of the integrators.