Improving the accuracy of simulated chaotic $N$-body orbits using smoothness

Abstract

Symplectic integrators are a foundation to the study of dynamical $N$-body phenomena, at scales ranging from from planetary to cosmological. These integrators preserve the Poincaré invariants of Hamiltonian dynamics. The $N$-body Hamiltonian has another, perhaps overlooked, symmetry: it is smooth, or, in other words, it has infinite differentiability class order (DCO) for particle separations greater than $0$. Popular symplectic integrators, such as hybrid methods or block adaptive stepping methods do not come from smooth Hamiltonians and it is perhaps unclear whether they should. We investigate the importance of this symmetry by considering hybrid integrators, whose DCO can be tuned easily. Hybrid methods are smooth, except at a finite number of phase space points. We study chaotic planetary orbits in a test considered by Wisdom. We find that increasing smoothness, at negligible extra computational cost in particular tests, improves the Jacobi constant error of the orbits by about $5$ orders of magnitude in long-term simulations. The results from this work suggest that smoothness of the $N$-body Hamiltonian is a property worth preserving in simulations.

arXiv comment9 pages, 8 figures. Matches accepted MNRAS version

References in corpus (10)

• The cosmological simulation code GADGET-2
• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• Hybrid Symplectic Integrators for Planetary Dynamics
• 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
• Variational Integrators for the Gravitational N-Body Problem
• 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?

Cited by in corpus (5)

• REBOUNDx: A Library for Adding Conservative and Dissipative Forces to Otherwise Symplectic N-body Integrations
• TRACE: a code for Time-Reversible Astrophysical Close Encounters
• Are long-term $N$-body simulations reliable?
• Switching integrators reversibly in the astrophysical $N$-body problem
• EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour