A Hybrid N-Body Code Incorporating Algorithmic Regularization and Post-Newtonian Forces
arXiv:0803.2310 · doi:10.1111/j.1365-2966.2008.13557.x
Abstract
We describe a novel N-body code designed for simulations of the central regions of galaxies containing massive black holes. The code incorporates Mikkola's 'algorithmic' chain regularization scheme including post-Newtonian terms up to PN2.5 order. Stars moving beyond the chain are advanced using a fourth-order integrator with forces computed on a GRAPE board. Performance tests confirm that the hybrid code achieves better energy conservation, in less elapsed time, than the standard scheme and that it reproduces the orbits of stars tightly bound to the black hole with high precision. The hybrid code is applied to two sample problems: the effect of finite-N gravitational fluctuations on the orbits of the S-stars; and inspiral of an intermediate-mass black hole into the galactic center.
12 pages, 15 figures, accepted for publication in MNRAS
References in corpus (28)
- Stellar Orbits Around the Galactic Center Black Hole
- The First Measurement of Spectral Lines in a Short-Period Star Bound to the Galaxy's Central Black Hole: A Paradox of Youth
- SINFONI in the Galactic Center: young stars and IR flares in the central light month
- Formation of Galactic Nuclei
- Resonant Relaxation in Stellar Systems
- Testing the general relativistic ''no-hair'' theorems using the galactic center black hole SgrA*
- The structure of the nuclear stellar cluster of the Milky Way
- Ejection of Supermassive Black Holes from Galaxy Cores
- The effect of mass-segregation on gravitational wave sources near massive black holes
- The dynamical evolution of massive black hole binaries - II. Self-consistent N-body integrations
- Implementing Few-Body Algorithmic Regularization with Post-Newtonian Terms
- GRAPE-6: The massively-parallel special-purpose computer for astrophysical particle simulation
- Performance Analysis of Direct N-Body Algorithms on Special-Purpose Supercomputers
- A class of symplectic integrators with adaptive timestep for separable Hamiltonian systems
- Stellar Dynamics at the Galactic Center with an Extremely Large Telescope
- Massive Black Holes in Star Clusters. I. Equal-mass clusters
- Probing Post-Newtonian Gravity near the Galactic Black Hole with Stellar Doppler Measurements
- Algorithmic regularization with velocity-dependent forces
- A Successful Targeted Search for Hypervelocity Stars
- Long-Term Evolution of Massive Black Hole Binaries. III. Binary Evolution in Collisional Nuclei
- Evolution of Galactic Nuclei. I. orbital evolution of IMBH
- Tracing intermediate-mass black holes in the Galactic Centre
- Resonant relaxation near a massive black hole: the dependence on eccentricity
- Reconstruction of Stellar Orbits Close to Sagittarius A*: Possibilities for Testing General Relativity
- Black Hole Binary Dynamics
- Post-Newtonian N-body simulations
- Models of Galaxies with Central Black Holes: Simulation Methods
- Rotational Brownian Motion of a Massive Binary
Cited by in corpus (28)
- The final-parsec problem in the collisionless limit
- N-body simulations of gravitational dynamics
- SAPPORO: A way to turn your graphics cards into a GRAPE-6
- Modeling Dense Star Clusters in the Milky Way and Beyond with the Cluster Monte Carlo Code
- Hypercompact Stellar Systems Around Recoiling Supermassive Black Holes
- Mergers and ejections of black holes in globular clusters
- On the distribution of stellar remnants around massive black holes: slow mass segregation, star cluster inspirals and correlated orbits
- Secular Stellar Dynamics near a Massive Black Hole
- The final-parsec problem in non-spherical galaxies revisited
- Post-Newtonian dynamical modeling of supermassive black holes in galactic-scale simulations
- N-body Integrators with Individual Time Steps from Hierarchical Splitting
- A new Fokker-Planck approach for relaxation-driven evolution of galactic nuclei
- The loss cone problem in axisymmetric nuclei
- A new code for orbit analysis and Schwarzschild modelling of triaxial stellar systems
- MSTAR -- a fast parallelised algorithmically regularised integrator with minimum spanning tree coordinates
- A new Monte Carlo method for dynamical evolution of non-spherical stellar systems
- body dynamics of Intermediate mass-ratio inspirals in globular clusters
- KETJU -- resolving small-scale supermassive black hole dynamics in GADGET-4
- Dynamical evolution of massive black holes in galactic-scale N-body simulations - introducing the regularized tree code "rVINE"
- MYRIAD: A new N-body code for simulations of Star Clusters
- A New Hybrid Technique for Modeling Dense Star Clusters
- GraviDy, a GPU modular, parallel direct-summation body integrator: Dynamics with softening
- Influence of tidal dissipation on outcomes of binary-single encounters between stars and black holes in stellar clusters
- Nuclear star clusters as probes of dark matter halos: the case of the Sagittarius Dwarf Spheroidal Galaxy
- How well do STARLAB and NBODY compare? II: Hardware and accuracy
- Survey of Multiscale and Multiphysics Applications and Communities
- Population Synthesis of Gravitational Wave Sources
- KRIOS: A new basis-expansion -body code for collisional stellar dynamics