Fourth-order leapfrog algorithms for numerical time evolution of classical and quantum systems

Abstract

Chau et al. [New J. Phys. 20, 073003 (2018)] presented a new and straight-forward derivation of a fourth-order approximation '$U_7$' of the time-evolution operator and hinted at its potential value as a symplectic integrator. $U_7$ is based on the Suzuki-Trotter split-operator method and leads to an algorithm for numerical time propagation that is superior to established methods. We benchmark the performance of $U_7$ and other algorithms, including a Runge-Kutta method and another recently developed Suzuki-Trotter-based scheme, that are exact up to fourth order in the evolution parameter, against various classical and quantum systems. We find $U_7$ to deliver any given target accuracy with the lowest computational cost, across all systems and algorithms tested here. This study is accompanied by open-source numerical software that we hope will prove valuable in the classroom.

arXiv comment14 pages, 6 figures; for accompanying open-source program, see https://github.com/huehou/Fourth-Order-Leapfrog

References in corpus (7)

• High order symplectic integrators for perturbed Hamiltonian systems
• On the construction of high-order force gradient algorithms for integration of motion in classical and quantum systems
• Numerical integration of variational equations
• Fourth-Order Algorithms for Solving the Imaginary Time Gross-Pitaevskii Equation in a Rotating Anisotropic Trap
• Symplectic fourth-order maps for the collisional N-body problem
• Chaos in the Honeycomb Optical Lattice Unit Cell
• First-principles quantum corrections for carrier correlations in double-layer two-dimensional heterostructures