English

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

Chemical Physics 2024-09-26 v2 Computational Physics Quantum Physics

Abstract

Geometric integrators of the Schr\"{o}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 101010^{-10}. In a companion paper [J. Roulet, S. Choi, and J. Van\'{\i}\v{c}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.

Keywords

Cite

@article{arxiv.1902.04661,
  title  = {Efficient geometric integrators for nonadiabatic quantum dynamics. I. The adiabatic representation},
  author = {Seonghoon Choi and Jiří Vaníček},
  journal= {arXiv preprint arXiv:1902.04661},
  year   = {2024}
}

Comments

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

R2 v1 2026-06-23T07:39:20.861Z