Related papers: KRIOS: A new basis-expansion $N$-body code for col…
Stellar systems are broadly divided into collisional and non-collisional. The latter are large-N systems with long relaxation timescales and can be simulated disregarding two-body interactions, while either computationally expensive direct…
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.…
Globular cluster stellar streams probe galaxy-formation processes and can potentially reveal the distribution of dark matter in galaxies. In many theoretical studies, streams are modeled with particle-spray or direct N-body codes. But…
We describe the public release of the Cluster Monte Carlo Code (CMC) a parallel, star-by-star $N$-body code for modeling dense star clusters. CMC treats collisional stellar dynamics using H\'enon's method, where the cumulative effect of…
We describe the astrophysical and numerical basis of N-body simulations, both of collisional stellar systems (dense star clusters and galactic centres) and collisionless stellar dynamics (galaxies and large-scale structure). We explain and…
One of the computational challenges of cluster formation simulations is resolving individual stars and simulating massive clusters with masses of more than $10^4 M_{\odot}$ without gravitational softening. Combining direct $N$-body code…
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 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)…
We study the dynamical evolution of globular clusters using our H\'enon-type Monte Carlo code for stellar dynamics including all relevant physics such as two-body relaxation, single and binary stellar evolution, Galactic tidal stripping,…
We present a hybrid technique of N-body simulation to deal with collisionless stellar systems having an inhomogeneous global structure. We combine a treecode and a self-consistent field code such that each of the codes model a different…
Debris discs are traditionally studied using two distinct types of numerical models: statistical particle-in-a-box codes to study their collisional and size distribution evolution, and dynamical N-body models to study their spatial…
We have developed a new simulation code aimed at studying the stellar dynamics of a galactic central star cluster surrounding a massive black hole. In order to include all the relevant physical ingredients (2-body relaxation, stellar mass…
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…
We present a new, simple, fast algorithm to numerically evolve disks of inelastically colliding particles surrounding a central star. Our algorithm adds negligible computational cost to the fastest existing collisionless N-body codes, and…
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…
In a previous paper we introduced a new method for simulating collisional gravitational $N$-body systems with linear time scaling on $N$, based on the Multi-Particle Collision (MPC) approach. This allows us to simulate globular clusters…
It is well known that numerical errors grow exponentially in $N$-body simulations of gravitational bound stellar systems, but it is not well understood how the accuracy parameters of algorithms affect the physical evolution in simulations.…
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 describe a new hybrid N-body/hydrodynamical code based on the particle-mesh (PM) method and the piecewise-parabolic method (PPM) for use in solving problems related to the evolution of large-scale structure, galaxy clusters, and…
We introduce a GPU-accelerated hybrid hydro/N-body code (Enzo-N) designed to address the challenges of concurrently simulating star clusters and their parent galaxies. This task has been exceedingly challenging, primarily due to the…