On the single-Hessian Gaussian wavepacket dynamics

Abstract

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}(\hslash )$ 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.

References in corpus (28)

• Gaussian Wavepacket Dynamics: semiquantal and semiclassical phase space formalism
• Computing quantum dynamics in the semiclassical regime
• 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
• On-the-fly ab initio semiclassical evaluation of vibronic spectra at finite temperature
• Semiclassical Approach to Photophysics Beyond Kasha's Rule and Vibronic Spectroscopy Beyond the Condon Approximation. The Case of Azulene
• 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 time-resolved electronic spectra
• Finite-temperature, anharmonicity, and Duschinsky effects on the two-dimensional electronic spectra from ab initio thermo-field Gaussian wavepacket dynamics
• On-the-fly ab initio semiclassical evaluation of third-order response functions for two-dimensional electronic spectroscopy
• On-the-fly ab initio three thawed Gaussians approximation: A semiclassical approach to Herzberg-Teller spectra
• 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
• Family of Gaussian wavepacket dynamics methods from the perspective of a nonlinear Schrödinger equation
• The Siegel Upper Half Space is a Marsden-Weinstein Quotient: Symplectic Reduction and Gaussian Wave Packets
• How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation?
• Internal Conversion Rates from the Extended Thawed Gaussian Approximation: Theory and Validation
• Symmetry and Conservation Laws in Semiclassical Wave Packet Dynamics
• High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics
• High-order geometric integrators for the variational Gaussian approximation
• Isotope effects in the electronic spectra of ammonia from ab initio semiclassical dynamics
• Prediction of Fluorescence Quantum Yields using the Extended Thawed Gaussian Approximation
• Thawed Gaussian wavepacket dynamics with $\Delta$-machine learned potentials
• On Hagedorn wavepackets associated with different Gaussians
• On-the-Fly Ab Initio Hagedorn Wavepacket Dynamics: Single Vibronic Level Fluorescence Spectra of Difluorocarbene
• Hamiltonian Dynamics of Semiclassical Gaussian Wave Packets in Electromagnetic Potentials