Symplectic fourth-order maps for the collisional N-body problem

Abstract

We study analytically and experimentally certain symplectic and time-reversible N-body integrators which employ a Kepler solver for each pair-wise interaction, including the method of Hernandez & Bertschinger (2015). Owing to the Kepler solver, these methods treat close two-body interactions correctly, while close three-body encounters contribute to the truncation error at second order and above. The second-order errors can be corrected to obtain a fourth-order scheme with little computational overhead. We generalise this map to an integrator which employs a Kepler solver only for selected interactions and yet retains fourth-order accuracy without backward steps. In this case, however, two-body encounters not treated via a Kepler solver contribute to the truncation error.

17 pages, re-submitted to MNRAS

References in corpus (8)

• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• N-body simulations of gravitational dynamics
• On the minimal accuracy required for simulating self-gravitating systems by means of direct N-body methods
• A Time-Symmetric Block Time-Step Algorithm for N-Body Simulations
• A Fast and Accurate Universal Kepler Solver without Stumpff Series
• Symplectic integration for the collisional gravitational $N$-body problem
• Fast and reliable symplectic integration for planetary system $N$-body problems
• A Keplerian-based Hamiltonian Splitting for Gravitational $N$-body Simulations

Cited by in corpus (20)

• The growth of intermediate mass black holes through tidal captures and tidal disruption events
• A study of symplectic integrators for planetary system problems: error analysis and comparisons
• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• TRACE: a code for Time-Reversible Astrophysical Close Encounters
• BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code
• Time-symmetric integration in astrophysics
• Are long-term $N$-body simulations reliable?
• Structure formation with primordial black holes: collisional dynamics, binaries, and gravitational waves
• A differentiable N-body code for transit timing and dynamical modeling. I. Algorithm and derivatives
• Stepsize errors in the $N$-body problem: discerning Mercury's true possible long-term orbits
• A rapid channel for the collisional formation and gravitational wave driven mergers of supermassive black hole seeds at high redshift
• 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
• Evolution of eccentric stellar disks around supermassive black holes: the complex disk disruption dynamics and the milliparsec stars
• Reversible time-step adaptation for the integration of few-body systems
• Improving the accuracy of simulated chaotic $N$-body orbits using smoothness
• Fourth-order leapfrog algorithms for numerical time evolution of classical and quantum systems
• DAWN. I. Simulating the formation and early evolution of stellar clusters with Phantom N-Body
• FROST-CLUSTERS -- III. Metallicity-dependent intermediate mass black hole formation by runaway collisions in dense star clusters