University of Leicester
Browse
MNRAS-2017-Dehnen-1201-17.pdf (2.22 MB)

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

Download (2.22 MB)
journal contribution
posted on 2017-01-04, 15:32 authored by Walter Dehnen, D. M. Hernandez
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.

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 Astronomy

Version

  • VoR (Version of Record)

Published in

Monthly Notices of the Royal Astronomical Society (February 11

Publisher

Oxford University Press (OUP), Royal Astronomical Society

issn

0035-8711

eissn

1365-2966

Acceptance date

2016-10-23

Available date

2017-01-04

Publisher version

http://mnras.oxfordjournals.org/content/465/1/1201

Language

en