Related papers: A New Hybrid Technique for Modeling Dense Star Clu…
The conservation of energy, linear momentum and angular momentum are important drivers for our physical understanding of the evolution of the Universe. These quantities are also conserved in Newton's laws of motion under gravity…
Using a Monte Carlo code, we construct a dynamic evolutionary model of the Galactic globular cluster M22 (NGC6656). The initial conditions are chosen so that, after about 12Gyr of stellar and dynamical evolution, the model is an approximate…
Approximate N-body methods, such as FastPM and COLA, have been successful in modelling halo and galaxy clustering statistics, but their low resolution on small scales is a limitation for applications that require high precision. Full N-body…
Dense stellar systems such as globular clusters, galactic nuclei and nuclear star clusters are ideal loci to study stellar dynamics due to the very high densities reached, usually a million times higher than in the solar neighborhood; they…
We recently introduced a new method for simulating collisional gravitational N-body systems with approximately linear time scaling with $N$, based on the Multi-Particle Collision (MPC) scheme, previously applied in Plasma Physics. We…
Several self-consistent models have been proposed, aiming at describing the phase space distribution of stars in globular clusters. This study explores the ability of the recently proposed LIMEPY models (Gieles & Zocchi) to reproduce the…
Globular clusters (GCs), the oldest stellar systems observed in the Milky Way, have for long been considered single stellar populations. As such, they provided an ideal laboratory to understand stellar dynamics and primordial star formation…
We present an algorithm for the simulation of the exact real-time dynamics of classical many-body systems with discrete energy levels. In the same spirit of kinetic Monte Carlo methods, a stochastic solution of the master equation is found,…
We present a general scheme for constructing Monte Carlo realizations of equilibrium, collisionless galaxy models with known distribution function (DF) f_0. Our method uses importance sampling to find the sampling DF f_s that minimizes the…
We simulate the assembly of a massive rich cluster and the formation of its constituent galaxies in a flat, low-density universe. Our most accurate model follows the collapse, the star-formation history and the orbital motion of all…
Gravitational waves and electromagnetic signals from merging neutron star binaries provide valuable information about the the properties of dense matter, the formation of heavy elements, and high-energy astrophysics. To fully leverage…
(Abridged) Dynamical friction can be used to distinguish Newtonian gravity and modified Newtonian dynamics (MOND) because it works differently in these frameworks. This concept, however, has yet to be explored very much with MOND. Previous…
We test the modified Newtonian dynamics (MOND) theory with the velocity-dispersion profiles of Galactic globular clusters populating the outermost region of the Milky Way halo, where the Galactic acceleration is lower than the…
The stability of a galaxy model is most easily assessed through N-body simulation. Particle-mesh codes have been widely used for this purpose, since they enable the largest numbers of particles to be employed. We show that the functional…
We describe Monte Carlo models for the dynamical evolution of the nearby globular cluster M4. The code includes treatments of two-body relaxation, three- and four-body interactions involving primordial binaries and those formed dynamically,…
Challenging stellar dynamical problems, such as the study of gravothermal oscillations in star clusters, have in the past initiated the very successful building of GRAPE special purpose computers. It is discussed, that present day tasks…
We describe source code level parallelization for the {\tt kira} direct gravitational $N$-body integrator, the workhorse of the {\tt starlab} production environment for simulating dense stellar systems. The parallelization strategy, called…
Advanced algorithms are necessary to obtain faster-than-real-time dynamic simulations in a number of different physical problems that are characterized by widely disparate time scales. Recent advanced dynamic Monte Carlo algorithms that…
Using our recently improved Monte Carlo evolution code, we study the evolution of the binary fraction in globular clusters. In agreement with previous N-body simulations, we find generally that the hard binary fraction in the core tends to…
We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied $f(R)$ gravity models. High resolution simulations for such models are…