A study of symplectic integrators for planetary system problems: error analysis and comparisons

Abstract

The symplectic Wisdom-Holman map revolutionized long-term integrations of planetary systems. There is freedom in such methods of how to split the Hamiltonian and which coordinate system to employ, and several options have been proposed in the literature. These choices lead to different integration errors, which we study analytically and numerically. The Wisdom-Holman method in Jacobi coordinates and the method of Hernandez, H16, compare favorably and avoid problems of some of the other maps, such as incorrect center-of-mass position or truncation errors even in the one-planet case. We use H16 to compute the evolution of Pluto's orbital elements over 500 million years in a new calculation.

24 pages, 9 figures, 1 table. Matches accepted MNRAS version

References in corpus (7)

• 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
• TTVFast: An efficient and accurate code for transit timing inversion problems
• A Fast and Accurate Universal Kepler Solver without Stumpff Series
• Symplectic fourth-order maps for the collisional N-body problem
• Symplectic integration for the collisional gravitational $N$-body problem
• Fast and reliable symplectic integration for planetary system $N$-body problems

Cited by in corpus (25)

• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• Hybrid Symplectic Integrators for Planetary Dynamics
• Construction of Explicit Symplectic Integrators in General Relativity. I. Schwarzschild Black Holes
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• Realistic collisional water transport during terrestrial planet formation: Self-consistent modeling by an N-body--SPH hybrid code
• TRACE: a code for Time-Reversible Astrophysical Close Encounters
• celmech: A Python package for celestial mechanics
• Are long-term $N$-body simulations reliable?
• The use of Kepler solver in numerical integrations of quasi-Keplerian orbits
• A differentiable N-body code for transit timing and dynamical modeling. I. Algorithm and derivatives
• Stepsize errors in the $N$-body problem: discerning Mercury's true possible long-term orbits
• Should $N$-body integrators be symplectic everywhere in phase space?
• Switching integrators reversibly in the astrophysical $N$-body problem
• orbitN: A symplectic integrator for planetary systems dominated by a central mass -- Insight into long-term solar system chaos
• EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour
• Improving the accuracy of simulated chaotic $N$-body orbits using smoothness
• Adjustment of force-gradient operator in symplectic methods
• WHFast512: A symplectic N-body integrator for planetary systems optimized with AVX512 instructions
• On the influence of equilibrium tides on transit-timing variations of close-in super-Earths. I. Application to single-planet systems and the case of K2-265 b
• Terrestrial Exoplanet Simulator (TES): an error optimal planetary systems integrator that permits close encounters
• An efficient integrator for stellar dynamics in effective gravity fields based on the isochrone potential
• A Criterion for the Onset of Chaos in Compact, Eccentric Multiplanet Systems
• 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
• A close-encounter method for simulating the dynamics of planetesimals