Related papers: Star Cluster Simulations: The State of the Art
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 first discuss two-body and chain regularization methods for direct N-body simulations on HARP-2 and GRAPE-6. The former is used for accurate integration of perturbed binaries and hierarchies, whereas the latter deals with strong…
The past few years have seen dramatic improvements in the scope and realism of star cluster simulations. Accurate treatments of stellar evolution, coupled with robust descriptions of all phases of binary evolution, have been incorporated…
Numerical simulations of star cluster formation have advanced greatly during the past decade, covering increasingly massive gas clouds while accounting for more and more complex physics. In this review, I discuss the present state of the…
Numerical integration methods are central to the study of self-gravitating systems, particularly those comprised of many bodies or otherwise beyond the reach of analytical methods. Predictor-corrector schemes, both multi-step methods and…
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 have performed fully self-consistent $N$-body simulations of star clusters near the Galactic center (GC). Such simulations have not been performed because it is difficult to perform fast and accurate simulations of such systems using…
The simulation of rich star clusters presents challenging problems of several kinds, including the design of suitable hardware and software, and numerous theoretical problems in stellar dynamics and stellar physics. Great progress has been…
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…
In this paper we describe a new algorithm for the long-term numerical integration of the two-body problem, in which two particles interact under a Newtonian gravitational potential. Although analytical solutions exist in the unperturbed and…
Spherically symmetric equal mass star clusters containing a large amount of primordial binaries are studied using a hybrid method, consisting of a gas dynamical model for single stars and a Monte Carlo treatment for relaxation of binaries…
We present a new hybrid Smoothed Particle Hydrodynamics (SPH)/N-body method for modelling the collisional stellar dynamics of young clusters in a live gas background. By deriving the equations of motion from Lagrangian mechanics we obtain a…
A symplectic integrator algorithm suitable for hierarchical triple systems is formulated and tested. The positions of the stars are followed in hierarchical Jacobi coordinates, whilst the planets are referenced purely to their primary. The…
The new approach outlined in Paper I (Spurzem \& Giersz 1996) to follow the individual formation and evolution of binaries in an evolving, equal point-mass star cluster is extended for the self-consistent treatment of relaxation and close…
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.…
The recent detection of gravitational waves and electromagnetic counterparts emitted during and after the collision of two neutron stars marks a breakthrough in the field of multi-messenger astronomy. Numerical relativity simulations are…
We describe a star cluster formation model that includes individual star formation from self-gravitating, magnetized gas, coupled to collisional stellar dynamics. The model uses the Astrophysical Multi-purpose Software Environment (AMUSE)…
Symplectic integrators are the tool of choice for many researchers studying dynamical systems because of their good long-term energy conservation properties. For systems with a dominant central mass, symplectic integrators are also highly…
A revision of Stodolkiewicz's Monte-Carlo code is used to simulate evolution of star clusters. The new method treats each superstar as a single star and follows the evolution and motion of all individual stellar objects. The first…
We study the bar instability in collisionless, rotating, anisotropic, stellar systems, using N-body simulations and also the matrix technique for calculation of modes with the perturbed collisionless Boltzmann equation. These methods are…