Related papers: A Hybrid N-Body Code Incorporating Algorithmic Reg…
We discuss the implementation of a new regular algorithm for simulation of the gravitational few-body problem. The algorithm uses components from earlier methods, including the chain structure, the logarithmic Hamiltonian, and the…
In this study, an $N$-body simulation code was developed for self-gravitating systems with a limited first-order post-Newtonian approximation. The code was applied to a special case in which the system consists of one massive object and…
This paper studies the formation and evolution of binary supermassive black holes (SMBHs) in rotating galactic nuclei, focusing on the role of stellar dynamics. We present the first N-body simulations that follow the evolution of the SMBHs…
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…
We report on results of fully consistent N-body simulations of globular cluster models with N = 100 000 members containing neutron stars and black holes. Using the improved `algorithmic regularization' method of Hellstrom and Mikkola for…
The "gravitational million-body problem," to model the dynamical evolution of a self-gravitating, collisional N-body system with ~10^6 particles over many relaxation times, remains a major challenge in computational astrophysics.…
Space-based observations are beginning to yield detailed information about the stellar kinematics at the very centers of galaxies, within the sphere of gravitational influence of the black hole. The structure and dynamics of these regions…
This paper presents a fast, economical particle-multiple-mesh N-body code optimized for large-N modelling of collisionless dynamical processes, such as black-hole wandering or bar-halo interactions, occurring within isolated galaxies. The…
We present a clear N-body realization of the growth of a Bahcall-Wolf f ~ E^0.25 (rho ~ 1/r^1.75) density cusp around a massive object ("black hole") at the center of a stellar system. Our N-body algorithm incorporates a novel…
Astrophysical Challenges which demand the solution of the one million (or more) gravitating body problem are briefly discussed for the fields of cosmology, galactic nuclei and globular star clusters. Results from the classical three-body…
We developed a new direct-tree hybrid N-body algorithm for fully self-consistent N-body simulations of star clusters in their parent galaxies. In such simulations, star clusters need high accuracy, while galaxies need a fast scheme because…
We describe efforts over the last six years to implement regularization methods suitable for studying one or more interacting black holes by direct N-body simulations. Three different methods have been adapted to large-N systems: (i)…
We review the recent optimizations of gravitational $N$-body kernels for running them on graphics processing units (GPUs), on single hosts and massive parallel platforms. For each of the two main $N$-body techniques, direct summation and…
The increasing role of general relativity in the dynamics of stellar systems with central massive black holes and in the evolution of hierarchical triple systems inspires a close examination of how post-Newtonian effects are incorporated…
We discuss the performance of direct summation codes used in the simulation of astrophysical stellar systems on highly distributed architectures. These codes compute the gravitational interaction among stars in an exact way and have an…
Special high-accuracy direct force summation N-body algorithms and their relevance for the simulation of the dynamical evolution of star clusters and other gravitating N-body systems in astrophysics are presented, explained and compared…
We follow the sinking of two massive black holes in a spherical stellar system where the black holes become bound under the influence of dynamical friction. Once bound, the binary hardens by three-body encounters with surrounding stars. We…
Motivated by the closest major merger, the Antennae Galaxies (NGC4038/4039), we want to improve our genetic algorithm based modeling code Minga (Theis 1999). The aim is to reveal the major interaction and galaxy parameters, e.g. orbital…
We present a new implementation of the numerical integration of the classical, gravitational, N-body problem based on a high order Hermite's integration scheme with block time steps, with a direct evaluation of the particle-particle forces.…
Direct-summation N-body algorithms compute the gravitational interaction between stars in an exact way and have a computational complexity of O(N^2). Performance can be greatly enhanced via the use of special-purpose accelerator boards like…