Related papers: Accelerating $N$-body simulation of self-gravitati…
The Kernel Polynomial Method (KPM) is one of the fast diagonalization methods used for simulations of quantum systems in research fields of condensed matter physics and chemistry. The algorithm has a difficulty to be parallelized on a…
The development of numerical N -body simulations have allowed to study formation process and evolution of galaxies at different scales. This paper presents the fundamental concepts of N-body systems applied to the cosmological evolution of…
A wide variety of outstanding problems in astrophysics involve the motion of a large number of particles ($N\gtrsim 10^{6}$) under the force of gravity. These include the global evolution of globular clusters, tidal disruptions of stars by…
We review what has been learned recently using N-body simulations about the evolution of globular clusters. While simulations of star clusters have become more realistic, and now include the evolution of single and binary stars, the…
We develop a formalism for General Relativistic N-body simulations in the weak field regime, suitable for cosmological applications. The problem is kept tractable by retaining the metric perturbations to first order, the first derivatives…
The increasing complexity and scale of cosmological N-body simulations, driven by astronomical surveys like Euclid, call for a paradigm shift towards more sustainable and energy-efficient high-performance computing (HPC). The rising energy…
The initial conditions for Newtonian $N$-body simulations are usually generated by applying the Zel'dovich approximation to the initial displacements of the particles using an initial power spectrum of density fluctuations generated by an…
We report a novel application of graphics processing units (GPUs) for the purpose of accelerating the search pipelines for gravitational waves from coalescing binaries of compact objects. A speed-up of 16 fold has been achieved compared…
This is the first in a series of papers devoted to fully general-relativistic $N$-body simulations applied to late-time cosmology. The purpose of this paper is to present the combination of a numerical relativity scheme, discretization…
A major goal of computational astrophysics is to simulate the Milky Way Galaxy with sufficient resolution down to individual stars. However, the scaling fails due to some small-scale, short-timescale phenomena, such as supernova explosions.…
N-body simulations are essential tools in physical cosmology to understand the large-scale structure (LSS) formation of the Universe. Large-scale simulations with high resolution are important for exploring the substructure of universe and…
A number of globular clusters appear to have undergone core collapse, in the sense that their predicted collapse time is much shorter than their current age. Simulations using gas models and Fokker-Planck approximation have shown that the…
We discuss the implementation of a new regular algorithm for simulation of the gravitational few-body problem. The algorithm uses components from earlier methods, including the chain structure, the logarithmic Hamiltonian, and the…
We present a simplified and fast method for simulating minor mergers between galaxy clusters. Instead of following the evolution of the dark matter halos directly by the N-body method, we employ a rigid potential approximation for both…
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 present MG-GLAM, a code developed for the very fast production of full $N$-body cosmological simulations in modified gravity (MG) models. We describe the implementation, numerical tests and first results of a large suite of cosmological…
$N$-body simulations study the dynamics of $N$ particles under the influence of mutual long-distant forces such as gravity. In practice, $N$-body codes will violate Newton's third law if they use either an approximate Poisson solver or…
High-resolution cosmological N-body simulations are excellent tools for modelling the formation and clustering of dark matter haloes. These simulations suggest complex physical theories of halo formation governed by a set of effective…
N-body simulations are the most powerful method to study the non-linear evolution of large-scale structure. However, they require large amounts of computational resources, making unfeasible their direct adoption in scenarios that require…
The role of gravity is crucial in astrophysics. It determines the evolution of any system, over an enormous range of time and space scales. Astronomical stellar systems as composed by N interacting bodies represent examples of…