Stepsize errors in the $N$-body problem: discerning Mercury's true possible long-term orbits

Abstract

Numerical integrations of the Solar System have been carried out for decades. Their results have been used, for example, to determine whether the Solar System is chaotic, whether Mercury's orbit is stable, or to help discern Earth's climate history. We argue that all of the past studies we consider in this work are affected by numerical chaos to different degrees, affecting the possible orbits and instability probability of Mercury, sometimes significantly. We show how to eliminate the effects of numerical chaos by resolving Mercury's pericentre passage. We also show that several higher order symplectic maps do not exhibit significant differences in resolving pericentre passage of Mercury (at fixed time step), making their advantages suspect for calculating long-term orbits. Resolving pericentre passage affects a wide array of orbital numerical studies, like exoplanet studies, studies of the galactic centre, and other $N$-body problems.

arXiv comment6 pages, 6 figures. Matches accepted MNRAS version

References in corpus (17)

• High order symplectic integrators for perturbed Hamiltonian systems
• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• Chaotic diffusion in the Solar System
• Long-term planetary integration with individual time steps
• Dynamical Chaos in the Wisdom-Holman Integrator: Origins and Solutions
• Chaotic Disintegration of the Inner Solar System
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• Dynamic stability of the Solar System: Statistically inconclusive results from ensemble integrations
• A study of symplectic integrators for planetary system problems: error analysis and comparisons
• Symplectic fourth-order maps for the collisional N-body problem
• Rare Event Sampling Improves Mercury Instability Statistics
• Highly stable evolution of Earth's future orbit despite chaotic behavior of the Solar System
• Time-symmetric integration in astrophysics
• Are long-term $N$-body simulations reliable?
• Should $N$-body integrators be symplectic everywhere in phase space?
• A repository of vanilla long term integrations of the Solar System
• Significant improvement in planetary system simulations from statistical averaging

Cited by in corpus (6)

• TRACE: a code for Time-Reversible Astrophysical Close Encounters
• Simple physics and integrators accurately reproduce Mercury instability statistics
• Recurring Activity Discovered on Quasi-Hilda 2009 DQ118
• orbitN: A symplectic integrator for planetary systems dominated by a central mass -- Insight into long-term solar system chaos
• Switching integrators reversibly in the astrophysical $N$-body problem
• WHFast512: A symplectic N-body integrator for planetary systems optimized with AVX512 instructions