Related papers: A Hybrid N-Body Code Incorporating Algorithmic Reg…
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…
We discuss the effects of the black holes' spin-multipole structure in the orbital dynamics of binary black holes according to general relativity, focusing on the leading-post-Newtonian-order couplings at each order in an expansion in the…
We present analysis of the evolution of dark matter halos in dense environments of groups and clusters in dissipationless cosmological simulations. The premature destruction of halos in such environments, known as the overmerging, reduces…
Direct $N$-body simulations of star clusters are accurate but expensive, largely due to the numerous $\mathcal{O} (N^2)$ pairwise force calculations. To solve the post-million-body problem, it will be necessary to use approximate force…
Inspiralling and coalescing binary black holes are promising sources of gravitational radiation. The orbital motion and gravitational-wave emission of such system can be modelled using a variety of approximation schemes and numerical…
We revisit a method to incorporate the Vainshtein screening mechanism in N-body simulations proposed by R. Scoccimarro in~\cite{Scoccimarro:2009eu}. We further extend this method to cover a subset of Horndeski theories that evade the bound…
Numerical effects are known to plague adaptive mesh refinement (AMR) codes when treating massive particles, e.g. representing massive black holes (MBHs). In an evolving background, they can experience strong, spurious perturbations and then…
The numerical simulations of massive collisional stellar systems, such as globular clusters (GCs), are very time-consuming. Until now, only a few realistic million-body simulations of GCs with a small fraction of binaries (5%) have been…
We revisit the accuracy of the post-Newtonian (PN) approximation and its region of validity for quasi-circular orbits of a point particle in the Kerr spacetime, by using an analytically known highest post-Newtonian order gravitational…
Nuclear star clusters surrounding supermassive black holes (SMBHs) in galactic nuclei contain large numbers of stars, black holes (BHs) and neutron stars (NSs), a fraction of which are likely to form binaries. These binaries were suggested…
For simulations that deal only with dark matter or stellar systems, the conventional N-body technique is fast, memory efficient, and relatively simple to implement. However when including the effects of gas physics, mesh codes are at a…
General relativistic precession limits the ability of gravitational encounters to increase the eccentricity $e$ of orbits near a supermassive black hole (SBH). This "Schwarzschild barrier" (SB) has been shown to play an important role in…
We investigate a hybrid numerical algorithm aimed at the large-scale cosmological N-body simulation for the on-going and the future high precious sky surveys. It makes use of a truncated Fast Multiple Method (FMM) for short-range gravity,…
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…
Constraining neutrino mass through cosmological observations relies on precise simulations to calibrate their effects on large scale structure, while these simulations must overcome computational challenges like dealing with large velocity…
We present a Fortran 95 code for simulating the evolution of astrophysical systems using particles to represent the underlying fluid flow. The code is designed to be versatile, flexible and extensible, with modular options that can be…
Most galaxies have supermassive black holes (SMBH) at their centres, surrounded by stars with binary systems also present in this environment. We use two schemes - post-Newtonian (PN) and a scalar perturbation to a background metric to…
Early-type galaxies exhibit a wealth of photometric and dynamical structures. These signatures are fossil records of their formation and evolution processes. In order to examine these structures in detail, we build models aimed at…
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…
Massive globular clusters lose stars via internal and external processes. Internal processes include mainly two-body relaxation, while external processes include interactions with the Galactic tidal field. We perform a suite of N-body…