High order symplectic integrators for perturbed Hamiltonian systems

Abstract

We present a class of symplectic integrators adapted for the integration of perturbed Hamiltonian systems of the form $H=A+\hat{I}\mu B$. We give a constructive proof that for all integer $p$, there exists an integrator with positive steps with a remainder of order $O(\ddot{I}{„}^p\hat{I}\mu +\ddot{I}{„}^2\hat{I}{\mu }^2)$, where $\ddot{I}„$ is the stepsize of the integrator. The analytical expressions of the leading terms of the remainders are given at all orders. In many cases, a corrector step can be performed such that the remainder becomes $O(\ddot{I}{„}^p\hat{I}\mu +\ddot{I}{„}^4\hat{I}{\mu }^2)$. The performances of these integrators are compared for the simple pendulum and the planetary 3-Body problem of Sun-Jupiter-Saturn.

24 pages, 6 figurres