Related papers: Towards optimal softening in 3D N-body codes: I. M…
We study the accumulation of errors in cosmological N-body algorithms that are caused by representing the continuous distribution of matter by massive particles, comparing the PPPM and Adaptive Multigrid codes. We use for this a new measure…
Planet-forming discs often contain structures like spiral arms, typically linked to the disc's gravitational forces. In 2D models, an ad hoc softening prescription is commonly used for self-gravity, but this overlooks the vertical…
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…
N-body simulations are essential for understanding the formation and evolution of structure in the Universe. However, the discrete nature of these simulations affects their accuracy when modelling collisionless systems. We introduce a new…
We discuss capability of Smooth Particle Hydrodynamics to represent adequately the dynamics of self-gravitating systems, in particular for what regards the quality of approximation of force fields in the motion equations. When cubic spline…
We perform a 3D reduction of the two-fermion Bethe-Salpeter equation, by series expansion around a positive-energy instantaneous approximation of the Bethe-Salpeter kernel, followed by another series expansion at the 3D level in order to…
We present a novel, fast and precise method for including the effect of light neutrinos in cosmological N-body simulations. The effect of the neutrino component is included by using the linear theory neutrino perturbations in the…
We present a two-dimensional (2D) Particle-Particle-Particle-Mesh (P$^3$M) algorithm with an optimized Green function and adaptive softening length for gravitational lensing studies in N-Body simulations. The analytical form of the…
A method for adaptive model order reduction for nonsmooth discrete element simulation is developed and analysed in numerical experiments. Regions of the granular media that collectively move as rigid bodies are substituted with rigid bodies…
We present a new time-stepping criterion for N-body simulations that is based on the true dynamical time of a particle. This allows us to follow the orbits of particles correctly in all environments since it has better adaptivity than…
We present a novel implementation of a soft sphere, discrete elements code to simulate the dynamics of self-gravitating granular materials. The code is used to study the outcome of sub-sonic collisions between self-gravitating rubble piles…
We present an algorithm for quickly generating multiple realizations of N-body simulations to be used, for example, for cosmological parameter estimation from surveys of large-scale structure. Our algorithm uses a new method to resample the…
We describe a NOVel form of Adaptive softening (NovA) for collisionless $N$-body simulations, implemented in the Ramses adaptive mesh refinement code. We introduce a refinement criterion that the particle distribution within each cell be…
In this work, we investigate the unjamming transition in a three-dimensional granular system composed of frictional spheres, in which the packing fraction is systematically reduced by random particle extractions. Using Discrete Element…
We revisit the r\^{o}le of discreteness and chaos in the dynamics of self-gravitating systems by means of $N$-body simulations with active and frozen potentials, starting from spherically symmetric stationary states and considering the…
We analyse with simple real-space statistics the Virgo consortium's cosmological N-body simulations. Significant clustering rapidly develops well below the initial mean interparticle separation \Lambda_i, where the gravitational force on a…
We introduce the N-body simulation technique to follow structure formation in linear and nonlinear regimes for the extended quintessence models (scalar-tensor theories in which the scalar field has a self-interaction potential and behaves…
I use N-body simulations with two mass species of particles to demonstrate that disk galaxy simulations are subject to collisional relaxation at a higher rate than is widely assumed. Relaxation affects the vertical thickness of the disk…
We investigate the range of applicability of a model for the real-space power spectrum based on N-body dynamics and a (quadratic) Lagrangian bias expansion. This combination uses the highly accurate particle displacements that can be…
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…