Related papers: Introducing a new multi-particle collision method …
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…
Core collapse is a prominent evolutionary stage of self-gravitating systems. In an idealised collisionless approximation, the region around the cluster core evolves in a self-similar way prior to the core collapse. Thus, its radial density…
We summarize a series of numerical experiments of collisional dynamics in dense stellar systems such as globular clusters (GCs) and in weakly collisional plasmas using a novel simulation technique, the so-called Multi-particle collision…
Star clusters with multi-mass components dynamically evolve faster than those modeled with equal-mass components. Using a series of direct $N$-body simulations, we investigate the dynamical evolution of star clusters with mass functions,…
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 study the early dynamical evolution of young, dense star clusters using Monte Carlo simulations for systems with up to N~10^7 stars. Rapid mass segregation of massive main-sequence stars and the development of the Spitzer instability can…
In this thesis we study several aspects of dynamical evolution of stellar clusters. The results of more than 200 simulations of single-mass star clusters with different initial total mass, half-mass radius and galactocentric distance, are…
We describe results from large numbers of $N$-body simulations containing from $250$ to $1000$ stars each. The distribution of stellar masses is a power law, and the systems are isolated. While the collapse of the core exhibits the expected…
We have recently written a new code to simulate the long term evolution of spherical clusters of stars. It is based on the pioneering Monte Carlo scheme proposed by Henon in the 70's. Our code has been devised in the specific goal to treat…
We present a new parallel code for computing the dynamical evolution of collisional N-body systems with up to N~10^7 particles. Our code is based on the the Henon Monte Carlo method for solving the Fokker-Planck equation, and makes…
We present a set of 148 independent $N$-body simulations of globular clusters (GCs) computed using the code $\texttt{CMC}$ ($\texttt{Cluster Monte Carlo}$). At an age of $\sim10-13\,$Gyr, the resulting models cover nearly the full range of…
We present a novel approach to the riddle of star cluster multiple populations. Stars form from molecular cores. But not all cores form stars. Following their initial compression, such 'failed' cores re-expand, rather than collapsing. We…
Recent simulations and observations suggest that star clusters form via the assembling of smaller sub-clusters. Because of their short relaxation time, sub-clusters experience core collapse much earlier than virialized solo-clusters, which…
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 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 present a new parallel supercomputer implementation of the Monte-Carlo method for simulating the dynamical evolution of globular star clusters. Our method is based on a modified version of Henon's Monte-Carlo algorithm for solving the…
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.…
We report on the results of a direct N-body simulation of a star cluster that started with N = 200 000, comprising 195 000 single stars and 5 000 primordial binaries. The code used for the simulation includes stellar evolution, binary…
We introduce version two of the fast star cluster evolution code Evolve Me A Cluster of StarS (EMACSS). The first version (Alexander & Gieles) assumed that cluster evolution is balanced for the majority of the life-cycle, meaning that the…
This paper, and its companion, investigate the evolution of dense stellar systems due to the influence of two-body gravitational encounters, physical collisions and stellar evolution. Our goal is the simulation of the densest centers of…