Are long-term $N$-body simulations reliable?

Abstract

$N$-body integrations are used to model a wide range of astrophysical dynamics, but they suffer from errors which make their orbits diverge exponentially in time from the correct orbits. Over long time-scales, their reliability needs to be established. We address this reliability by running a three-body planetary system over about $200$ e-folding times. Using nearby initial conditions, we can construct statistics of the long-term phase-space structure and compare to rough estimates of resonant widths of the system. We compared statistics for a wide range of numerical methods, including a Runge--Kutta method, Wisdom--Holman method, symplectic corrector methods, and a method by Laskar and Robutel. "Improving" an integrator did not increase the phase space accuracy, but simply increasing the number of initial conditions did. In fact, the statistics of a higher order symplectic corrector method were inconsistent with the other methods in one test.

arXiv comment14 pages, 13 figures, 4 tables. Matches accepted MNRAS version

References in corpus (21)

• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• 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
• Dynamical Chaos in the Wisdom-Holman Integrator: Origins and Solutions
• A Criterion for the Onset of Chaos in Systems of Two Eccentric Planets
• High precision Symplectic Integrators for the Solar System
• On the minimal accuracy required for simulating self-gravitating systems by means of direct N-body methods
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• A Fast and Accurate Universal Kepler Solver without Stumpff Series
• A study of symplectic integrators for planetary system problems: error analysis and comparisons
• Symplectic fourth-order maps for the collisional N-body problem
• Symplectic integration for the collisional gravitational $N$-body problem
• Numerical verification of the microscopic time reversibility of Newton's equations of motion: Fighting exponential divergence
• Fast and reliable symplectic integration for planetary system $N$-body problems
• Time-symmetric integration in astrophysics
• Should $N$-body integrators be symplectic everywhere in phase space?
• JANUS: A bit-wise reversible integrator for N-body dynamics
• On the accuracy of symplectic integrators for secularly evolving planetary systems
• Improving the accuracy of simulated chaotic $N$-body orbits using smoothness
• Shadowing unstable orbits of the Sitnikov elliptic 3-body problem

Cited by in corpus (11)

• High-order implicit time integration scheme based on Padé expansions
• Chaos in self-gravitating many-body systems: Lyapunov time dependence of $N$ and the influence of general relativity
• High-order implicit time integration scheme with controllable numerical dissipation based on mixed-order Padé expansions
• Stepsize errors in the $N$-body problem: discerning Mercury's true possible long-term orbits
• Energy-conserving integrator for conservative Hamiltonian systems with ten-dimensional phase space
• Simple physics and integrators accurately reproduce Mercury instability statistics
• Switching integrators reversibly in the astrophysical $N$-body problem
• EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour
• Extreme Gravitational Interactions in the Problem of Three Black Holes in General Relativity
• An efficient integrator for stellar dynamics in effective gravity fields based on the isochrone potential
• Significant improvement in planetary system simulations from statistical averaging