Significant improvement in planetary system simulations from statistical averaging

Abstract

Symplectic integrators are widely used for the study of planetary dynamics and other $N$-body problems. In a study of the outer Solar system, we demonstrate that individual symplectic integrations can yield biased errors in the semi-major axes and possibly other orbital elements. The bias is resolved by studying an ensemble of initial conditions of the outer Solar system. Such statistical sampling could significantly improve measurement of planetary system properties like their secular frequencies. We also compared the distributions of action-like variables between high and low accuracy integrations; traditional statistical metrics are unable to distinguish the distribution functions.

arXiv comment4 pages, 1 figure. Accepted by RNAAS

References in corpus (8)

• La2010: A new orbital solution for the long term motion of the Earth
• High order symplectic integrators for perturbed Hamiltonian systems
• New families of symplectic splitting methods for numerical integration in dynamical astronomy
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• A study of symplectic integrators for planetary system problems: error analysis and comparisons
• Are long-term $N$-body simulations reliable?
• EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour
• Terrestrial Exoplanet Simulator (TES): an error optimal planetary systems integrator that permits close encounters

Cited by in corpus (1)

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