A connected component-based method for efficiently integrating multiscale $N$-body systems

Abstract

We present a novel method for efficient direct integration of gravitational N-body systems with a large variation in characteristic time scales. The method is based on a recursive and adaptive partitioning of the system based on the connected components of the graph generated by the particle distribution combined with an interaction-specific time step criterion. It uses an explicit and approximately time-symmetric time step criterion, and conserves linear and angular momentum to machine precision. In numerical tests on astrophysically relevant setups, the method compares favourably to both alternative Hamiltonian-splitting integrators as well as recently developed block time step-based GPU-accelerated Hermite codes. Our reference implementation is incorporated in the HUAYNO code, which is freely available as a part of the AMUSE framework.

12 pages, 8 figures, accepted for A&A

References in corpus (15)

• Young massive star clusters
• Multi-physics simulations using a hierarchical interchangeable software interface
• The Astrophysical Multipurpose Software Environment
• The dynamical evolution of fractal star clusters: the survival of substructure
• High Performance Direct Gravitational N-body Simulations on Graphics Processing Units -- II: An implementation in CUDA
• SAPPORO: A way to turn your graphics cards into a GRAPE-6
• 6th and 8th Order Hermite Integrator for N-body Simulations
• High Performance Direct Gravitational N-body Simulations on Graphics Processing Units
• High Performance Direct Gravitational N-body Simulations on Graphics Processing Unit I: An implementation in Cg
• N-body Integrators with Individual Time Steps from Hierarchical Splitting
• A fully parallel, high precision, N-body code running on hybrid computing platforms
• Particle-Particle Particle-Tree: A Direct-Tree Hybrid Scheme for Collisional N-Body Simulations
• On the minimal accuracy required for simulating self-gravitating systems by means of direct N-body methods
• A Keplerian-based Hamiltonian Splitting for Gravitational $N$-body Simulations
• MYRIAD: A new N-body code for simulations of Star Clusters

Cited by in corpus (9)

• The origin of interstellar asteroidal objects like 1I/2017 U1 'Oumuamua
• Survivability of planetary systems in young and dense star clusters
• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code
• The 800pc long tidal tails of the Hyades star cluster: Possible discovery of candidate epicyclic overdensities from an open star cluster
• Oort Cloud Ecology. III. The Sun left the parent star cluster shortly after the giant planets formed
• Venice: a multi-scale operator-splitting algorithm for multi-physics simulations
• Interstellar objects follow the collapse of molecular clouds
• Tracing Milky Way substructure with an RR Lyrae hierarchical clustering forest