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.

arXiv comment18 pages, 7 figures. Accepted for publication in MNRAS

References in corpus (42)

• On the variation of the Initial Mass Function
• The cosmological simulation code GADGET-2
• Evolution of binary stars and the effect of tides on binary populations
• Comprehensive analytic formulae for stellar evolution as a function of mass and metallicity
• WHFast: A fast and unbiased implementation of a symplectic Wisdom-Holman integrator for long term gravitational simulations
• Simulating cosmic structure formation with the GADGET-4 code
• The DRAGON simulations: globular cluster evolution with a million stars
• Accelerating NBODY6 with Graphics Processing Units
• NBODY6++GPU: Ready for the gravitational million-body problem
• Mass segregation and fractal substructure in young massive clusters: (I) the McLuster code and method calibration
• Implementing Few-Body Algorithmic Regularization with Post-Newtonian Terms
• A class of symplectic integrators with adaptive timestep for separable Hamiltonian systems
• PeTar: a high-performance N-body code for modeling massive collisional stellar systems
• N-body simulations of gravitational dynamics
• Long-term planetary integration with individual time steps
• SAPPORO: A way to turn your graphics cards into a GRAPE-6
• 6th and 8th Order Hermite Integrator for N-body Simulations
• Algorithmic regularization with velocity-dependent forces
• Secular dynamics of hierarchical multiple systems composed of nested binaries, with an arbitrary number of bodies and arbitrary hierarchical structure. First applications to multiplanet and multistar systems
• Post-Newtonian dynamical modeling of supermassive black holes in galactic-scale simulations
• N-body Integrators with Individual Time Steps from Hierarchical Splitting
• 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
• High order symplectic integrators for planetary dynamics and their implementation in REBOUND
• A Time-Symmetric Block Time-Step Algorithm for N-Body Simulations
• A Fast and Accurate Universal Kepler Solver without Stumpff Series
• Gravitational softening as a smoothing operation
• Symplectic fourth-order maps for the collisional N-body problem
• Variational Integrators for the Gravitational N-Body Problem
• Symplectic integration for the collisional gravitational $N$-body problem
• Implementing Tidal and Gravitational Wave Energy Losses in Few-body Codes: A Fast and Easy Drag Force Model
• 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
• Fast Multipole Methods for $N$-body Simulations of Collisional Star Systems
• EnckeHH: an integrator for gravitational dynamics with a dominant mass that achieves optimal error behaviour
• Embedded operator splitting methods for perturbed systems

Cited by in corpus (19)

• 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
• A rapid channel for the collisional formation and gravitational wave driven mergers of supermassive black hole seeds at high redshift
• Efficiency of black hole formation via collisions in stellar systems: An analysis of data from simulations and observations
• Formation and evolution of binary black holes in $N$-body simulations of star clusters with up to two million stars
• Fast Multipole Methods for $N$-body Simulations of Collisional Star Systems
• Evolution of eccentric stellar disks around supermassive black holes: the complex disk disruption dynamics and the milliparsec stars
• Evolution of Massive Black Hole Binaries in Collisionally Relaxed Nuclear Star Clusters -- Impact of Mass Segregation
• Population Synthesis of Gravitational Wave Sources
• DAWN. I. Simulating the formation and early evolution of stellar clusters with Phantom N-Body
• KRIOS: A new basis-expansion $N$-body code for collisional stellar dynamics
• FROST-CLUSTERS -- III. Metallicity-dependent intermediate mass black hole formation by runaway collisions in dense star clusters