BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code

Abstract

We present BIFROST, an extended version of the GPU-accelerated hierarchical fourth-order forward symplectic integrator code FROST. BIFROST (BInaries in FROST) can efficiently evolve collisional stellar systems with arbitrary binary fractions up to $f_{\mathrm{bin}}=100\%$ by using secular and regularised integration for binaries, triples, multiple systems or small clusters around black holes within the fourth-order forward integrator framework. Post-Newtonian (PN) terms up to order PN3.5 are included in the equations of motion of compact subsystems with optional three-body and spin-dependent terms. PN1.0 terms for interactions with black holes are computed everywhere in the simulation domain. The code has several merger criteria (gravitational-wave inspirals, tidal disruption events and stellar and compact object collisions) with the addition of relativistic recoil kicks for compact object mergers. We show that for systems with $N$ particles the scaling of the code remains good up to $N_{\mathrm{GPU}}\sim 40\times N/10^6$ GPUs and that the increasing binary fractions up to 100 per cent hardly increase the code running time (less than a factor $\sim 1.5$). We also validate the numerical accuracy of BIFROST by presenting a number of star clusters simulations the most extreme ones including a core collapse and a merger of two intermediate mass black holes with a relativistic recoil kick.

arXiv comment25 pages, 16 figures, accepted for publication in MNRAS

References in corpus (93)

• Observation of Gravitational Waves from a Binary Black Hole Merger
• GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral
• On the variation of the Initial Mass Function
• Embedded Clusters in Molecular Clouds
• Binary interaction dominates the evolution of massive stars
• Evolution of binary stars and the effect of tides on binary populations
• A Survey of Stellar Families: Multiplicity of Solar-Type Stars
• Comprehensive analytic formulae for stellar evolution as a function of mass and metallicity
• GADGET: A code for collisionless and gasdynamical cosmological simulations
• Stellar Multiplicity
• Mind your Ps and Qs: the Interrelation between Period (P) and Mass-ratio (Q) Distributions of Binary Stars
• Young massive star clusters
• Compact Remnant Mass Function: Dependence on the Explosion Mechanism and Metallicity
• The Eccentric Kozai-Lidov Effect and Its Applications
• Coalescing binary systems of compact objects to (post)$^{5/2}-Newtonian order. V. Spin Effects
• IAS15: A fast, adaptive, high-order integrator for gravitational dynamics, accurate to machine precision over a billion orbits
• On The Maximum Mass of Stellar Black Holes
• The Gravitational Million-Body Problem
• The Effects of Binary Evolution on the Dynamics of Core Collapse and Neutron-Star Kicks
• 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
• Secular Dynamics in Hierarchical Three-Body Systems
• The Effect of Pair-Instability Mass Loss on Black Hole Mergers
• Nuclear Star Clusters
• From binaries to multiples II: hierarchical multiplicity of F and G dwarfs
• Toward Complete Statistics of Massive Binary Stars: Penultimate Results from the Cygnus OB2 Radial Velocity Survey
• The DRAGON simulations: globular cluster evolution with a million stars
• NBODY6++GPU: Ready for the gravitational million-body problem
• Third post-Newtonian accurate generalized quasi-Keplerian parametrization for compact binaries in eccentric orbits
• Relativistic Dynamics and Extreme Mass Ratio Inspirals
• Prospects for detection of detached double white dwarf binaries with Gaia, LSST and LISA
• Modeling Dense Star Clusters in the Milky Way and Beyond with the $\text{CMC}$ Cluster Catalog
• From binaries to multiples I: Data on F and G dwarfs within 67 pc of the Sun
• Black holes and core expansion in massive star clusters
• Mass segregation and fractal substructure in young massive clusters: (I) the McLuster code and method calibration
• Relativistic Binaries in Globular Clusters
• Implementing Few-Body Algorithmic Regularization with Post-Newtonian Terms
• STARFORGE: Toward a comprehensive numerical model of star cluster formation and feedback
• Relative Frequencies of Blue Stragglers in Galactic Globular Clusters: Constraints for the Formation Mechanisms
• Third post-Newtonian dynamics of compact binaries: Equations of motion in the center-of-mass frame
• An Updated Look at Binary Characteristics of Massive Stars in the Cygnus OB2 Association
• 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
• SPHGal: Smoothed Particle Hydrodynamics with improved accuracy for Galaxy simulations
• BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences
• Formation of Massive Black Holes in Galactic Nuclei: Runaway Tidal Encounters
• The core binary fractions of star clusters from realistic simulations
• Algorithmic regularization with velocity-dependent forces
• Intermediate Mass Black Hole Formation in compact Young Massive Star Clusters
• Hydrodynamical Neutron-star Kicks in Electron-capture Supernovae and Implications for the CRAB Supernova
• Mergers and ejections of black holes in globular clusters
• Binary stars in the Galactic thick disc
• 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
• On incorporating post-Newtonian effects in N-body dynamics
• Star Clusters with Primordial Binaries: I. Dynamical Evolution of Isolated Models
• Radii of Young Star Clusters in Nearby Galaxies
• Post-Newtonian dynamical modeling of supermassive black holes in galactic-scale simulations
• N-body Integrators with Individual Time Steps from Hierarchical Splitting
• The formation of extremely diffuse galaxy cores by merging supermassive black holes
• VINE -- A numerical code for simulating astrophysical systems using particles I: Description of the physics and the numerical methods
• Stellar Collisions and Blue Straggler Stars in Dense Globular Clusters
• Modeling the remnant mass, spin, and recoil from unequal-mass, precessing black-hole binaries: The Intermediate Mass Ratio Regime
• A slow-down time-transformed symplectic integrator for solving the few-body problem
• Black hole mergers in compact star clusters and massive black hole formation beyond the mass-gap
• Multiple Stellar Evolution: a population synthesis algorithm to model the stellar, binary, and dynamical evolution of multiple-star systems
• The Formation and Early Evolution of Young Massive Clusters
• Modelling of the Effects of Stellar Feedback during Star Cluster Formation Using a Hybrid Gas and N-Body Method
• Formation and evolution of young massive clusters in galaxy mergers: the SMUGGLE view
• Relativistic three-body effects in hierarchical triples
• MSTAR -- a fast parallelised algorithmically regularised integrator with minimum spanning tree coordinates
• Binaries in star clusters and the origin of the field stellar population
• Secular and tidal evolution of circumbinary systems
• SIRIUS project. I. Star formation models for star-by-star simulations of star clusters and galaxy formation
• Signatures of the Many Supermassive Black Hole Mergers in a Cosmologically Forming Massive Early-Type Galaxy
• A Fast and Accurate Universal Kepler Solver without Stumpff Series
• Algebraic and machine learning approach to hierarchical triple-star stability
• Resolving the Complex Evolution of a Supermassive Black Hole Triplet in a Cosmological Simulation
• The Keplerian three-body encounter I. Insights on the origin of the S-stars and the G-objects in the Galactic center
• Structure and rotation of young massive star clusters in a simulated dwarf starburst
• Symplectic fourth-order maps for the collisional N-body problem
• The formation and early evolution of embedded star clusters in spiral galaxies
• FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
• Gravitational Waves from the Inspiral of Supermassive Black Holes in Galactic-scale Simulations
• 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
• SpaceHub: A high-performance gravity integration toolkit for few-body problems in astrophysics
• MYRIAD: A new N-body code for simulations of Star Clusters
• A general-purpose timestep criterion for simulations with gravity
• A connected component-based method for efficiently integrating multiscale $N$-body systems
• Fast Multipole Methods for $N$-body Simulations of Collisional Star Systems
• Pairing mechanisms for binary stars
• The impact of Post-Newtonian effects on massive black hole binary evolution at $\sim 1000R_{sch}$ separations

Cited by in corpus (15)

• Structure formation with primordial black holes: collisional dynamics, binaries, and gravitational waves
• Rapid formation of a very massive star >50000 $M_{\odot }$ and subsequently an IMBH from runaway collisions. Direct N-body and Monte Carlo simulations of dense star clusters
• 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
• Little Red Dots are Tidal Disruption Events in Runaway-Collapsing Clusters
• Novel Conservative Methods for Adaptive Force Softening in Collisionless and Multi-Species N-Body Simulations
• Evolution of a black hole cluster in full general relativity
• Evolution of eccentric stellar disks around supermassive black holes: the complex disk disruption dynamics and the milliparsec stars
• Mergers all the way down: stellar collisions and kinematics of a dense hierarchically forming massive star cluster in a dwarf starburst
• Structure Formation with Warm White Noise: Effects of Finite Number Density and Velocity Dispersion in Particle and Wave Dark Matter
• Dynamical evolution of the open clusters with different star formation efficiencies and orbital parameters
• DAWN. I. Simulating the formation and early evolution of stellar clusters with Phantom N-Body
• Radiation magnetohydrodynamics simulations of Population III star formation during the Epoch of Reionization
• FROST-CLUSTERS -- III. Metallicity-dependent intermediate mass black hole formation by runaway collisions in dense star clusters