Single-Hessian Gaussian wavepacket dynamics (GWD) significantly reduces the computational burden of Heller's local harmonic GWD, while maintaining comparable accuracy in approximating vibronic spectra. Here, we provide a new, symplectic derivation of the equations of motion of single-Hessian GWD and show that, unlike the local harmonic version, this method conserves the non-canonical symplectic structure on the manifold of Gaussian wavepackets and $-$ for bounded dynamics in smooth potentials $-$ avoids the drift of energy. Our numerical results suggest that, despite being much more efficient than the local harmonic variant, the single-Hessian GWD exhibits the same $\mathcal{O}(\hbar)$ asymptotic error in averages of observables. To further accelerate numerical simulations, we implement high-order time-stepping geometric integrators that are time-reversible and conserve the norm and symplectic structure exactly, regardless of the time step. In addition, we present explicit expressions for the exact evolution of the width of a single-Hessian Gaussian wavepacket in a general potential, as well as for the exact evolution of the whole wavepacket in a global harmonic potential. Using on-the-fly ab initio Gaussian wavepacket dynamics on the first excited-state surface of ammonia, we numerically confirm the conservation of geometric properties by these integrators and demonstrate that high-order integrators can enhance both accuracy and computational efficiency. We also compute the photoelectron spectrum of the difluorocarbene anion and the absorption spectrum of methylamine, and find that, in comparison with experiment, single-Hessian GWD outperforms global harmonic models and matches the accuracy of local harmonic GWD. Finally, we identify which spectral features are sensitive to the choice of reference Hessian.Single-Hessian Gaussian wavepacket dynamics (GWD) significantly reduces the computational burden of Heller's local harmonic GWD, while maintaining comparable accuracy in approximating vibronic spectra. Here, we provide a new, symplectic derivation of the equations of motion of single-Hessian GWD and show that, unlike the local harmonic version, this method conserves the non-canonical symplectic structure on the manifold of Gaussian wavepackets and$-$for bounded dynamics in smooth potentials$-$avoids the drift of energy. Our numerical results suggest that, despite being much more efficient than the local harmonic variant, the single-Hessian GWD exhibits the same $\mathcal{O}(\hbar)$ asymptotic error in averages of observables. To further accelerate numerical simulations, we implement high-order time-stepping geometric integrators that are time-reversible and conserve the norm and symplectic structure exactly, regardless of the time step. In addition, we present explicit expressions for the exact evolution of the width of a single-Hessian Gaussian wavepacket in a general potential, as well as for the exact evolution of the whole wavepacket in a global harmonic potential. Using on-the-fly ab initio Gaussian wavepacket dynamics on the first excited-state surface of ammonia, we numerically confirm the conservation of geometric properties by these integrators and demonstrate that high-order integrators can enhance both accuracy and computational efficiency. We also compute the photoelectron spectrum of the difluorocarbene anion and the absorption spectrum of methylamine, and find that, in comparison with experiment, single-Hessian GWD outperforms global harmonic models and matches the accuracy of local harmonic GWD. Finally, we identify which spectral features are sensitive to the choice of reference Hessian.

physics.chem-phCross-listedarXiv:2605.07440May 11, 2026

Linear dependence and energy conservation in Gaussian wavepacket basis sets

Scott Habershon

We propose a method for dealing with the problem of linear dependence in quantum dynamics simulations employing over-complete Gaussian wavepacket (GWP) basis sets. In particular, by periodically projecting out redundant basis functions using the matching pursuit algorithm whilst simultaneously introducing GWPs which avoid linear dependence with the current basis set, we find that numerical conditioning of the equations-of-motion can be readily controlled. In applications to particle tunnelling in one- and two-dimensional potentials, this method allows us to reproduce the exact quantum-mechanical results with fewer GWP basis functions than similar calculations with non-adaptive basis sets, a result which we trace back to the improved energy conservation of our adaptive approach.We propose a method for dealing with the problem of linear dependence in quantum dynamics simulations employing over-complete Gaussian wavepacket (GWP) basis sets. In particular, by periodically projecting out redundant basis functions using the matching pursuit algorithm whilst simultaneously introducing GWPs which avoid linear dependence with the current basis set, we find that numerical conditioning of the equations-of-motion can be readily controlled. In applications to particle tunnelling in one- and two-dimensional potentials, this method allows us to reproduce the exact quantum-mechanical results with fewer GWP basis functions than similar calculations with non-adaptive basis sets, a result which we trace back to the improved energy conservation of our adaptive approach.

The Journal of Chemical PhysicsJan 4, 201245 citations

Gaussian wave packets

Aleksei Prokhorovich Kiselev, Maria V. Perel

Optics and SpectroscopyMar 1, 199923 citations

-

for bounded dynamics in smooth potentials

-

avoids the drift of energy. Our numerical results suggest that, despite being much more efficient than the local harmonic variant, the single-Hessian GWD exhibits the same

\mathcal{O}(\hbar)

asymptotic error in averages of observables. To further accelerate numerical simulations, we implement high-order time-stepping geometric integrators that are time-reversible and conserve the norm and symplectic structure exactly, regardless of the time step. In addition, we present explicit expressions for the exact evolution of the width of a single-Hessian Gaussian wavepacket in a general potential, as well as for the exact evolution of the whole wavepacket in a global harmonic potential. Using on-the-fly ab initio Gaussian wavepacket dynamics on the first excited-state surface of ammonia, we numerically confirm the conservation of geometric properties by these integrators and demonstrate that high-order integrators can enhance both accuracy and computational efficiency. We also compute the photoelectron spectrum of the difluorocarbene anion and the absorption spectrum of methylamine, and find that, in comparison with experiment, single-Hessian GWD outperforms global harmonic models and matches the accuracy of local harmonic GWD. Finally, we identify which spectral features are sensitive to the choice of reference Hessian.