MNRAS-2017-Dehnen-1201-17.pdf (2.22 MB)
Symplectic fourth-order maps for the collisional N-body problem
journal contribution
posted on 2017-01-04, 15:32 authored by Walter Dehnen, D. M. HernandezWe 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.
Funding
This work used the DiRAC Complexity system, operated by the University of Leicester IT Services, which forms part of the STFC DiRAC HPC Facility (www.dirac.ac.uk). This equipment is funded by BIS National E-Infrastructure capital grant ST/K000373/1 and STFC DiRAC Operations grant ST/K0003259/1. DiRAC is part of the National E-Infrastructure.
History
Citation
Monthly Notices of the Royal Astronomical Society (February 11, 2017) 465 (1): 1201-1217.Author affiliation
/Organisation/COLLEGE OF SCIENCE AND ENGINEERING/Department of Physics and AstronomyVersion
- VoR (Version of Record)