High order symplectic integrators for planetary dynamics and their implementation in REBOUND

Abstract

Direct N-body simulations and symplectic integrators are effective tools to study the long-term evolution of planetary systems. The Wisdom-Holman (WH) integrator in particular has been used extensively in planetary dynamics as it allows for large timesteps at good accuracy. One can extend the WH method to achieve even higher accuracy using several different approaches. In this paper we survey integrators developed by Wisdom et al. (1996), Laskar & Robutel (2001, and Blanes et al. (2013). Since some of these methods are harder to implement and not as readily available to astronomers compared to the standard WH method, they are not used as often. This is somewhat unfortunate given that in typical simulations it is possible to improve the accuracy by up to six orders of magnitude (!) compared to the standard WH method without the need for any additional force evaluations. To change this, we implement a variety of high order symplectic methods in the freely available N-body integrator REBOUND. In this paper we catalogue these methods, discuss their differences, describe their error scalings, and benchmark their speed using our implementations.

arXiv comment10 pages, 2 figures, accepted by MNRAS, code to reproduce figures available at https://github.com/hannorein/ReinTamayoBrown2019

References in corpus (9)

• 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
• New families of symplectic splitting methods for numerical integration in dynamical astronomy
• Dynamical Chaos in the Wisdom-Holman Integrator: Origins and Solutions
• A new paradigm for reproducing and analyzing N-body simulations of planetary systems
• A study of symplectic integrators for planetary system problems: error analysis and comparisons

Cited by in corpus (30)

• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• On the long-term stability of the Solar System in the presence of weak perturbations from stellar flybys
• Rare Event Sampling Improves Mercury Instability Statistics
• Revisiting high-order Taylor methods for astrodynamics and celestial mechanics
• TRACE: a code for Time-Reversible Astrophysical Close Encounters
• The influence of the secular perturbation of an intermediate-mass companion: I. Eccentricity excitation of disk stars at the Galactic center
• The Stability Boundary of the Distant Scattered Disk
• Resonance in the K2-19 system is at odds with its high reported eccentricities
• Are long-term $N$-body simulations reliable?
• Stepsize errors in the $N$-body problem: discerning Mercury's true possible long-term orbits
• Stellar Evolution and Tidal Dissipation in REBOUNDx
• On the accuracy of symplectic integrators for secularly evolving planetary systems
• The influence of the secular perturbation of an intermediate-mass companion: II. Ejection of hypervelocity stars from the Galactic Center
• General relativistic precession and the long-term stability of the solar system
• Worlds Next Door: A Candidate Giant Planet Imaged in the Habitable Zone of $\alpha$ Cen A. I. Observations, Orbital and Physical Properties, and Exozodi Upper Limits
• A repository of vanilla long term integrations of the Solar System
• EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour
• The Yarkovsky effect in REBOUNDx
• DMPP-3: confirmation of short-period S-type planet(s) in a compact eccentric binary star system, and warnings about long-period RV planet detections
• A low-mass sub-Neptune planet transiting the bright active star HD 73344
• An implicit symplectic solver for high-precision long term integrations of the Solar System
• Close-in compact super-Earth systems emerging from resonant chains: slow destabilization by unseen remnants of formation
• Embedded operator splitting methods for perturbed systems
• The possible origin of three Apollo asteroids 3200 Phaethon, 2005UD, and 1999YC
• Faster Algorithmic Quantum and Classical Simulations by Corrected Product Formulas
• Uncovering the Rapidly Evolving Orbits of the Dynamic TOI-201 System
• An efficient integrator for stellar dynamics in effective gravity fields based on the isochrone potential
• Significant improvement in planetary system simulations from statistical averaging
• Galactic tides in the Solar System within a non-axisymmetric Milky Way model adjusted to Gaia data