Symplectic integration for the collisional gravitational $N$-body problem

Abstract

We present a new symplectic integrator designed for collisional gravitational $N$-body problems which makes use of Kepler solvers. The integrator is also reversible and conserves 9 integrals of motion of the $N$-body problem to machine precision. The integrator is second order, but the order can easily be increased by the method of \citeauthor{yos90}. We use fixed time step in all tests studied in this paper to ensure preservation of symplecticity. We study small $N$ collisional problems and perform comparisons with typically used integrators. In particular, we find comparable or better performance when compared to the 4th order Hermite method and much better performance than adaptive time step symplectic integrators introduced previously. We find better performance compared to SAKURA, a non-symplectic, non-time-reversible integrator based on a different two-body decomposition of the $N$-body problem. The integrator is a promising tool in collisional gravitational dynamics.

11 pages, 7 Figures. Minor changes made and appendix added. To appear in MNRAS: (September 11, 2015) 452 (2): 1934-1944

References in corpus (3)

• The fine-grained phase-space structure of Cold Dark Matter halos
• On the minimal accuracy required for simulating self-gravitating systems by means of direct N-body methods
• A Keplerian-based Hamiltonian Splitting for Gravitational $N$-body Simulations

Cited by in corpus (17)

• The nature of the TRAPPIST-1 exoplanets
• Magnetic activity in the HARPS M-dwarf sample. The rotation-activity relationship for very low-mass stars through R'HK
• 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
• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• TRACE: a code for Time-Reversible Astrophysical Close Encounters
• Kepler Object of Interest Network I. First results combining ground and space-based observations of Kepler systems with transit timing variations
• Fast and reliable symplectic integration for planetary system $N$-body problems
• Time-symmetric integration in astrophysics
• Are long-term $N$-body simulations reliable?
• A differentiable N-body code for transit timing and dynamical modeling. I. Algorithm and derivatives
• Should $N$-body integrators be symplectic everywhere in phase space?
• High order regularised symplectic integrator for collisional planetary systems
• Capture of interstellar objects I: the capture cross-section
• Switching integrators reversibly in the astrophysical $N$-body problem
• Understanding the evolution of Atira-class asteroid 2019 AQ3, a major step towards the future discovery of the Vatira population