Efficient fourth order symplectic integrators for near-harmonic separable Hamiltonian systems
Abstract
Efficient fourth order symplectic integrators are proposed for numerical integration of separable Hamiltonian systems H(p,q)=T(p)+V(q). Symmetric splitting coefficients with five to nine stages are obtained by higher order decomposition of the simple harmonic oscillator. The performance of the methods is evaluated for various Hamiltonian systems: Integration errors are compared to those of acclaimed integrators composed by S. Blanes et al. (2013), W. Kahan et al. (1999) and H. Yoshida (1990). Numerical tests indicate that the integrators obtained in this paper perform significantly better than previous integrators for common Hamiltonian systems.
- Publication:
-
arXiv e-prints
- Pub Date:
- January 2015
- DOI:
- arXiv:
- arXiv:1501.04345
- Bibcode:
- 2015arXiv150104345M
- Keywords:
-
- Quantum Physics;
- Physics - Computational Physics
- E-Print:
- 9 pages, 6(9) figures, 4 tables. Utilized software available at http://zymplectic.com [Revised benchmark data and corrected typo/figure layout]