Multidimensional quantum dynamics with explicitly correlated Gaussian wave packets using Rothe's method
arXiv:2411.05459 · doi:10.1063/5.0247732
Abstract
In a previous publication [J. Chem. Phys., 161, 044105 (2024)], it has been shown that Rothe's method can be used to solve the time-dependent Schrödinger equation (TDSE) for the hydrogen atom in a strong laser field using time-dependent Gaussian wave packets. Here, we generalize these results, showing that Rothe's method can propagate arbitrary numbers of thawed, complex-valued, explicitly correlated Gaussian functions (ECGs) with dense correlation matrices for systems with varying dimensionality. We consider the multidimensional Henon-Heiles potential, and show that the dynamics can be quantitatively reproduced using only 30 Gaussians in 2D, and that accurate spectra can be obtained using 20 Gaussians in 2D and 30 to 40 Gaussians in 3D and 4D. Thus, the relevant multidimensional dynamics can be described at high quality using only a small number of ECGs that give a very compact representation of the wave function. This efficient representation, along with the demonstrated ability of Rothe's method to propagate Gaussian wave packets in strong fields and ECGs in complex potentials, paves the way for accurate molecular dynamics calculations beyond the Born-Oppenheimer approximation in strong fields.
17 pages, 13 figures
References in corpus (15)
- SciPy 1.0--Fundamental Algorithms for Scientific Computing in Python
- Barren plateaus in quantum neural network training landscapes
- Multilayer multi-configuration time-dependent Hartree method: implementation and applications to a Henon-Heiles Hamiltonian and to pyrazine
- Multiconfiguration Time-Dependent Hartree-Fock Treatment of Electronic and Nuclear Dynamics in Diatomic Molecules
- Large-scale quantum-dynamics with matrix product states
- Modeling electron dynamics coupled to continuum states in finite volumes
- A general approach to few-cycle intense laser interactions with complex atoms
- A tensor network view of multilayer multiconfiguration time-dependent Hartree methods
- An optimized absorbing potential for ultrafast, strong-field problems
- Regularized relativistic corrections for polyelectronic and polyatomic systems with explicitly correlated Gaussians
- Various variational approximations of quantum dynamics
- Gaussians for Electronic and Rovibrational Quantum Dynamics
- Time evolution as an optimization problem: The hydrogen atom in strong laser fields in a basis of time-dependent Gaussian wave packets
- No need for a grid: Adaptive fully-flexible gaussians for the time-dependent Schrödinger equation
- Thawed Gaussian wave packet dynamics: a critical assessment of three propagation schemes