Particle-Particle Particle-Tree: A Direct-Tree Hybrid Scheme for Collisional N-Body Simulations
Abstract
In this paper, we present a new hybrid algorithm for the time integration of collisional N-body systems. In this algorithm, the gravitational force between two particles is divided into short- and long-range terms, using a distance-dependent cutoff function. The long-range interaction is calculated using a tree algorithm, and integrated with a constant-timestep leapfrog integrator. The short-range term is calculated directly and integrated with the high-order Hermite scheme. We can reduce the calculation cost per orbital period from O(N2) to O(NlogN), without significantly increasing the long-term integration error. The results of our test simulations show that close encounters are integrated accurately. Long-term errors of the total energy show random-walk behaviour, because they are dominated by the error caused by the tree approximation.
- Publication:
-
Publications of the Astronomical Society of Japan
- Pub Date:
- August 2011
- DOI:
- arXiv:
- arXiv:1101.5504
- Bibcode:
- 2011PASJ...63..881O
- Keywords:
-
- methods: n-body simulations;
- solar system: formation;
- Astrophysics - Earth and Planetary Astrophysics
- E-Print:
- 22 pages, 15 figures