Journal
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY
Volume 493, Issue 3, Pages 3398-3411Publisher
OXFORD UNIV PRESS
DOI: 10.1093/mnras/staa480
Keywords
methods: numerical; software: simulations; globular clusters: general
Categories
Funding
- JSPS (School of Science, The University of Tokyo)
Ask authors/readers for more resources
An accurate and efficient method dealing with the few-body dynamics is important for simulating collisional N-body systems like star clusters and to follow the formation and evolution of compact binaries. We describe such a method which combines the time-transformed explicit symplectic integrator and the slow-down method. The former conserves the Hamiltonian and the angular momentum for a long-term evolution, while the latter significantly reduces the computational cost for a weakly perturbed binary. In this work, the Hamilton equations of this algorithm are analysed in detail. We mathematically and numerically show that it can correctly reproduce the secular evolution like the orbit averaged method and also well conserve the angular momentum. For a weakly perturbed binary, the method is possible to provide a few orders of magnitude faster performance than the classical algorithm. A publicly available code written in the C++ language, SDAR, is available on GITHUB. It can be used either as a standalone tool or a library to be plugged in other N-body codes. The high precision of the floating point to 62 digits is also supported.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available