High-order geometric integrators for the variational Gaussian approximation
arXiv:2306.17608 · doi:10.1063/5.0165489
Abstract
Among the single-trajectory Gaussian-based methods for solving the time-dependent Schrödinger equation, the variational Gaussian approximation is the most accurate one. In contrast to Heller's original thawed Gaussian approximation, it is symplectic, conserves energy exactly, and may partially account for tunneling. However, the variational method is also much more expensive. To improve its efficiency, we symmetrically compose the second-order symplectic integrator of Faou and Lubich and obtain geometric integrators that can achieve an arbitrary even order of convergence in the time step. We demonstrate that the high-order integrators can speed up convergence drastically compared to the second-order algorithm and, in contrast to the popular fourth-order Runge-Kutta method, are time-reversible and conserve the norm and the symplectic structure exactly, regardless of the time step. To show that the method is not restricted to low-dimensional systems, we perform most of the analysis on a non-separable twenty-dimensional model of coupled Morse oscillators. We also show that the variational method may capture tunneling and, in general, improves accuracy over the non-variational thawed Gaussian approximation.
arXiv comment17 pages, 11 figures
References in corpus (21)
- Gaussian Wavepacket Dynamics: semiquantal and semiclassical phase space formalism
- Computing quantum dynamics in the semiclassical regime
- Geometry of variational methods: dynamics of closed quantum systems
- Semiclassical "Divide-and-Conquer" Method for Spectroscopic Calculations of High Dimensional Molecular Systems
- On-the-fly ab initio semiclassical dynamics: Identifying degrees of freedom essential for emission spectra of oligothiophenes
- On-the-fly ab initio semiclassical evaluation of absorption spectra of polyatomic molecules beyond the Condon approximation
- Single-Hessian thawed Gaussian approximation: The missing rung on the ladder
- Strong coupling theory for tunneling and vibrational relaxation in driven bistable systems
- Approximating tunneling rates in multi-dimensional field spaces
- Relation of exact Gaussian basis methods to the dephasing representation: Theory and application to time-resolved electronic spectra
- On-the-fly ab initio semiclassical evaluation of third-order response functions for two-dimensional electronic spectroscopy
- Vibrational Investigation of Nucleobases by Means of Divide and Conquer Semiclassical Dynamics
- Symplectic Semiclassical Wave Packet Dynamics
- Efficient geometric integrators for nonadiabatic quantum dynamics. II. The diabatic representation
- Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation
- The Siegel Upper Half Space is a Marsden-Weinstein Quotient: Symplectic Reduction and Gaussian Wave Packets
- Problem-free time-dependent variational principle for open quantum systems
- High-order geometric integrators for representation-free Ehrenfest dynamics
- Calculation of eigenvalues of a strongly chaotic system using Gaussian wavepacket dynamics
- Validity of the time-dependent variational approximation to the Gaussian wavepacket method applied to double-well systems
- Wave packets in QFT: leading order width corrections to decay rates and clock behaviour under Lorentz boosts