On the accuracy of symplectic integrators for secularly evolving planetary systems

Abstract

Symplectic integrators have made it possible to study the long-term evolution of planetary systems with direct N-body simulations. In this paper we reassess the accuracy of such simulations by running a convergence test on 20Myr integrations of the Solar System using various symplectic integrators. We find that the specific choice of metric for determining a simulation's accuracy is important. Only looking at metrics related to integrals of motions such as the energy error can overestimate the accuracy of a method. As one specific example, we show that symplectic correctors do not improve the accuracy of secular frequencies compared to the standard Wisdom-Holman method without symplectic correctors, despite the fact that the energy error is three orders of magnitudes smaller. We present a framework to trace the origin of this apparent paradox to one term in the shadow Hamiltonian. Specifically, we find a term that leads to negligible contributions to the energy error but introduces non-oscillatory errors that result in artificial periastron precession. This term is the dominant error when determining secular frequencies of the system. We show that higher order symplectic methods such as the Wisdom-Holman method with a modified kernel or the SABAC family of integrators perform significantly better in secularly evolving systems because they remove this specific term.

arXiv comment13 pages, 7 figures, accepted for publication in MNRAS

References in corpus (10)

• REBOUND: An open-source multi-purpose N-body code for collisional dynamics
• La2010: A new orbital solution for the long term motion of the Earth
• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• High order symplectic integrators for perturbed Hamiltonian systems
• Hybrid Symplectic Integrators for Planetary Dynamics
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• Frequency map analysis and quasiperiodic decompositions
• A new paradigm for reproducing and analyzing N-body simulations of planetary systems
• An orbital window into the ancient Sun's mass

Cited by in corpus (9)

• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• Born eccentric: constraints on Jupiter and Saturn's pre-instability orbits
• On the long-term stability of the Solar System in the presence of weak perturbations from stellar flybys
• Are long-term $N$-body simulations reliable?
• Energy-conserving integrator for conservative Hamiltonian systems with ten-dimensional phase space
• General relativistic precession and the long-term stability of the solar system
• Nodal Precession in Closely Spaced Planet Pairs
• EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour
• Embedded operator splitting methods for perturbed systems