MYRIAD: A new N-body code for simulations of Star Clusters
arXiv:1006.3326 · doi:10.1051/0004-6361/200913890
Abstract
We present a new C++ code for collisional N-body simulations of star clusters. The code uses the Hermite fourth-order scheme with block time steps, for advancing the particles in time, while the forces and neighboring particles are computed using the GRAPE-6 board. Special treatment is used for close encounters, binary and multiple sub-systems that either form dynamically or exist in the initial configuration. The structure of the code is modular and allows the appropriate treatment of more physical phenomena, such as stellar and binary evolution, stellar collisions and evolution of close black-hole binaries. Moreover, it can be easily modified so that the part of the code that uses GRAPE-6, could be replaced by another module that uses other accelerating-hardware like the Graphics Processing Units (GPUs). Appropriate choice of the free parameters give a good accuracy and speed for simulations of star clusters up to and beyond core collapse. Simulations of Plummer models consisting of equal-mass stars reached core collapse at t~17 half-mass relaxation times, which compares very well with existing results, while the cumulative relative error in the energy remained below 0.001. Also, comparisons with published results of other codes for the time of core collapse for different initial conditions, show excellent agreement. Simulations of King models with an initial mass-function, similar to those found in the literature, reached core collapse at t~0.17, which is slightly smaller than the expected result from previous works. Finally, the code accuracy becomes comparable and even better than the accuracy of existing codes, when a number of close binary systems is dynamically created in a simulation. This is due to the high accuracy of the method that is used for close binary and multiple sub-systems.
24 pages, 29 figures, accepted for publication to Astronomy & Astrophysics
References in corpus (31)
- Large Merger Recoils and Spin Flips From Generic Black-Hole Binaries
- Maximum gravitational recoil
- Star cluster ecology IVa: Dissection of an open star cluster---photometry
- Supermassive recoil velocities for binary black-hole mergers with antialigned spins
- Ejection of Supermassive Black Holes from Galaxy Cores
- Gravitational Radiation from Post-Newtonian Sources and Inspiralling Compact Binaries
- High Performance Direct Gravitational N-body Simulations on Graphics Processing Units -- II: An implementation in CUDA
- GRAPE-6: The massively-parallel special-purpose computer for astrophysical particle simulation
- Modeling kicks from the merger of generic black-hole binaries
- Remnant Masses, Spins and Recoils from the Merger of Generic Black-Hole Binaries
- Statistics of N-body Simulations I. Equal Masess Before Core Collapse
- Which Globular Clusters contain Intermediate-mass Black Holes?
- SAPPORO: A way to turn your graphics cards into a GRAPE-6
- Binary Black Holes: Spin Dynamics and Gravitational Recoil
- Gravitational Wave Recoil and the Retention of Intermediate Mass Black Holes
- 6th and 8th Order Hermite Integrator for N-body Simulations
- Evidence for primordial mass segregation in globular clusters
- Relativistic Suppression of Black Hole Recoils
- Star Clusters with Primordial Binaries: I. Dynamical Evolution of Isolated Models
- Superkicks in Hyperbolic Encounters of Binary Black Holes
- A new Monte Carlo code for star cluster simulations: I. Relaxation
- Modeling gravitational recoil from precessing highly-spinning unequal-mass black-hole binaries
- High Performance Direct Gravitational N-body Simulations on Graphics Processing Unit I: An implementation in Cg
- High Performance Direct Gravitational N-body Simulations on Graphics Processing Units
- Dynamics of compact objects clusters: A post-Newtonian study
- Star Cluster Simulations: The State of the Art
- A Hybrid N-Body Code Incorporating Algorithmic Regularization and Post-Newtonian Forces
- Gravitational waves from eccentric intermediate-mass black hole binaries
- Star Cluster Ecology: VII The evolution of young dense star clusters containing primordial binaries
- Post-Newtonian simulations of super-massive black hole binaries in galactic nuclei
- How well do STARLAB and NBODY4 compare? I: Simple models
Cited by in corpus (13)
- N-body Integrators with Individual Time Steps from Hierarchical Splitting
- Compact Binary Coalescences: Astrophysical Processes and Lessons Learned
- Investigating the retention of intermediate-mass black holes in star clusters using N-body simulations
- Swarm-NG: a CUDA Library for Parallel n-body Integrations with focus on Simulations of Planetary Systems
- FROST: a momentum-conserving CUDA implementation of a hierarchical fourth-order forward symplectic integrator
- A Performance Comparison of Different Graphics Processing Units Running Direct N-Body Simulations
- BIFROST: simulating compact subsystems in star clusters using a hierarchical fourth-order forward symplectic integrator code
- A hybrid SPH/N-body method for star cluster simulations
- A connected component-based method for efficiently integrating multiscale -body systems
- Sowing the seeds of massive black holes in small galaxies: Young clusters as the building blocks of Ultra-Compact-Dwarf Galaxies
- GraviDy, a GPU modular, parallel direct-summation body integrator: Dynamics with softening
- How well do STARLAB and NBODY compare? II: Hardware and accuracy
- Direct -body code on low-power embedded ARM GPUs