@citekhatri

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

, , and . (2019)cite arxiv:1910.08667Comment: 13 pages, 12 figures, to be submitted to MNRAS, comments welcome.

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. Our statistics are approximately consistent for a wide range of numerical methods, including a Runge--Kutta method, Wisdom--Holman method, symplectic corrector methods, and a method by Laskar & Robutel. "Improving" an integrator did not affect the phase space accuracy, but simply increasing the number of initial conditions did.

Description

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

Links and resources

Tags