FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator

Abstract

We present a novel hierarchical formulation of the fourth-order forward symplectic integrator and its numerical implementation in the GPU-accelerated direct-summation N-body code FROST. The new integrator is especially suitable for simulations with a large dynamical range due to its hierarchical nature. The strictly positive integrator sub-steps in a fourth-order symplectic integrator are made possible by computing an additional gradient term in addition to the Newtonian accelerations. All force calculations and kick operations are synchronous so the integration algorithm is manifestly momentum-conserving. We also employ a time-step symmetrisation procedure to approximately restore the time-reversibility with adaptive individual time-steps. We demonstrate in a series of binary, few-body and million-body simulations that FROST conserves energy to a level of $|\Delta E/E|\sim 10^{-10}$ while errors in linear and angular momentum are practically negligible. For typical star cluster simulations, we find that FROST scales well up to $N_{\mathrm{GPU}}^{\max }\sim 4\times N/10^5$ GPUs, making direct summation N-body simulations beyond $N=10^6$ particles possible on systems with several hundred and more GPUs. Due to the nature of hierarchical integration the inclusion of a Kepler solver or a regularised integrator with post-Newtonian corrections for close encounters and binaries in the code is straightforward.

18 pages, 7 figures. Accepted for publication in MNRAS

References in corpus (17)

• Implementing Few-Body Algorithmic Regularization with Post-Newtonian Terms
• PeTar: a high-performance N-body code for modeling massive collisional stellar systems
• N-body simulations of gravitational dynamics
• 6th and 8th Order Hermite Integrator for N-body Simulations
• Algorithmic regularization with velocity-dependent forces
• A slow-down time-transformed symplectic integrator for solving the few-body problem
• Multiple Stellar Evolution: a population synthesis algorithm to model the stellar, binary, and dynamical evolution of multiple-star systems
• GANDALF - Graphical Astrophysics code for N-body Dynamics And Lagrangian Fluids
• The fate of planetesimal discs in young open clusters: implications for 1I/'Oumuamua, the Kuiper belt, the Oort cloud and more
• MSTAR -- a fast parallelised algorithmically regularised integrator with minimum spanning tree coordinates
• A Time-Symmetric Block Time-Step Algorithm for N-Body Simulations
• The physics of symplectic integrators: perihelion advances and symplectic corrector algorithms
• Time-symmetric integration in astrophysics
• Towards time symmetric N-body integration
• MYRIAD: A new N-body code for simulations of Star Clusters
• A connected component-based method for efficiently integrating multiscale $N$-body systems
• Embedded operator splitting methods for perturbed systems

Cited by in corpus (18)

• Large-scale dark matter simulations
• Modeling Dense Star Clusters in the Milky Way and Beyond with the Cluster Monte Carlo Code
• Preparing the next gravitational million-body simulations: Evolution of single and binary stars in Nbody6++GPU, MOCCA and McLuster
• The challenge of simulating the star cluster population of dwarf galaxies with resolved interstellar medium
• The growth of intermediate mass black holes through tidal captures and tidal disruption events
• Code Comparison in Galaxy Scale Simulations with Resolved Supernova Feedback: Lagrangian vs. Eulerian Methods
• KETJU -- resolving small-scale supermassive black hole dynamics in GADGET-4
• BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code
• Structure formation with primordial black holes: collisional dynamics, binaries, and gravitational waves
• Efficiency of black hole formation via collisions in stellar systems: An analysis of data from simulations and observations
• A rapid channel for the collisional formation and gravitational wave driven mergers of supermassive black hole seeds at high redshift
• Formation and evolution of binary black holes in $N$-body simulations of star clusters with up to two million stars
• Evolution of Massive Black Hole Binaries in Collisionally Relaxed Nuclear Star Clusters -- Impact of Mass Segregation
• Evolution of eccentric stellar disks around supermassive black holes: the complex disk disruption dynamics and the milliparsec stars
• DAWN. I. Simulating the formation and early evolution of stellar clusters with Phantom N-Body
• Population Synthesis of Gravitational Wave Sources
• FROST-CLUSTERS -- III. Metallicity-dependent intermediate mass black hole formation by runaway collisions in dense star clusters
• KRIOS: A new basis-expansion $N$-body code for collisional stellar dynamics