Symplectic integration for the collisional gravitational
Name
Symplectic integration.pdf
Size
434.8 KB
Format
Adobe PDF
Checksum (MD5)
170baf564db618ad40c343a98fd3019b
Author(s)
Hernandez, David Michael
Bertschinger, Edmund
Date Issued
July 2015
Journal
Monthly Notices of the Royal Astronomical Society
Publisher
Oxford University Press
Citation
Hernandez, David M., and Edmund Bertschinger. “Symplectic Integration for the Collisional Gravitational N -Body Problem.” Monthly Notices of the Royal Astronomical Society 452.2 (2015): 1934–1944.
Version
Author's final manuscript
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 nine 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 Yoshida. 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 fourth-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.
MIT Department
Massachusetts Institute of Technology. Department of Physics
MIT Kavli Institute for Astrophysics and Space Research
Terms of Use
Creative Commons Attribution-Noncommercial-Share Alike
Persistent DSpace Link
DOI of Published Version
https://doi.org/10.1093/mnras/stv1439
