Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation

Abstract

Geometric integrators of the Schrödinger equation conserve exactly many invariants of the exact solution. Among these integrators, the split-operator algorithm is explicit and easy to implement, but, unfortunately, is restricted to systems whose Hamiltonian is separable into a kinetic and potential terms. Here, we describe several implicit geometric integrators applicable to both separable and non-separable Hamiltonians, and, in particular, to the nonadiabatic molecular Hamiltonian in the adiabatic representation. These integrators combine the dynamic Fourier method with recursive symmetric composition of the trapezoidal rule or implicit midpoint method, which results in an arbitrary order of accuracy in the time step. Moreover, these integrators are exactly unitary, symplectic, symmetric, time-reversible, and stable, and, in contrast to the split-operator algorithm, conserve energy exactly, regardless of the accuracy of the solution. The order of convergence and conservation of geometric properties are proven analytically and demonstrated numerically on a two-surface NaI model in the adiabatic representation. Although each step of the higher order integrators is more costly, these algorithms become the most efficient ones if higher accuracy is desired; a thousand-fold speedup compared to the second-order trapezoidal rule (the Crank-Nicolson method) was observed for wavefunction convergence error of $10^{-10}$. In a companion paper [J. Roulet, S. Choi, and J. Van\'ıček (2019)], we discuss analogous, arbitrary-order compositions of the split-operator algorithm and apply both types of geometric integrators to a higher-dimensional system in the diabatic representation.

title modified (now part I of a series), old Sec. II.F on split-operator algorithms replaced with new Sec. II.F on diabatic-adiabatic transformation, added new Fig. 5 on convergence of phase, old Figs. 6 and 7 on geometric properties merged into new Fig. 7, Propositions 5 and 6 on convergence of SOD added

References in corpus (1)

• Efficient geometric integrators for nonadiabatic quantum dynamics. II. The diabatic representation

Cited by in corpus (19)

• Grid-based methods for chemistry simulations on a quantum computer
• 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
• Efficient geometric integrators for nonadiabatic quantum dynamics. II. The diabatic representation
• Strong Error Bounds for Trotter & Strang-Splittings and Their Implications for Quantum Chemistry
• Family of Gaussian wavepacket dynamics methods from the perspective of a nonlinear Schrödinger equation
• How important are the residual nonadiabatic couplings for an accurate simulation of nonadiabatic quantum dynamics in a quasidiabatic representation?
• High-order geometric integrators for representation-free Ehrenfest dynamics
• Which form of the molecular Hamiltonian is the most suitable for simulating the nonadiabatic quantum dynamics at a conical intersection?
• Nonadiabatic conical intersection dynamics in the local diabatic representation with Strang splitting and Fourier basis
• High-order geometric integrators for the variational Gaussian approximation
• High-order geometric integrators for the local cubic variational Gaussian wavepacket dynamics
• A time-reversible integrator for the time-dependent Schrödinger equation on an adaptive grid
• Finite-temperature vibronic spectra from the split-operator coherence thermofield dynamics
• Ehrenfest dynamics accelerated with SPEED
• Applicability of the thawed Gaussian wavepacket dynamics to the calculation of vibronic spectra of molecules with double-well potential energy surfaces
• Time-reversible and norm-conserving high-order integrators for the nonlinear time-dependent Schrödinger equation: Application to local control theory
• On the single-Hessian Gaussian wavepacket dynamics
• An implicit split-operator algorithm for the nonlinear time-dependent Schrödinger equation