Related papers: General relativistic corrections to $N$-body simul…
We discuss the relation between the output of Newtonian N-body simulations on scales that approach or exceed the particle horizon to the description of General Relativity. At leading order, the Zeldovich approximation is correct on large…
The standard numerical tools for studying non-linear collapse of matter are Newtonian $N$-body simulations. Previous work has shown that these simulations are in accordance with General Relativity (GR) up to first order in perturbation…
Initial conditions for (Newtonian) cosmological N-body simulations are usually set by re-scaling the present-day power spectrum obtained from linear (relativistic) Boltzmann codes to the desired initial redshift of the simulation. This…
We use gauge-invariant cosmological perturbation theory to calculate the displacement field that sets the initial conditions for $N$-body simulations. Using first and second-order fully relativistic perturbation theory in the…
We examine the deviation of Cold Dark Matter particle trajectories from the Newtonian result as the size of the region under study becomes comparable to or exceeds the particle horizon. To first order in the gravitational potential, the…
On large-scales, comparable to the horizon, the observable clustering properties of galaxies are affected by various general relativistic effects. To calculate these effects one needs to consistently solve for the metric, densities and…
We describe the irrotational dust component of the universe in terms of a relativistic gradient expansion and transform the resulting synchronous metric to a Newtonian coordinate system. The two metrics are connected via a space-like…
The set-up of the initial conditions in cosmological N-body simulations is usually implemented by rescaling the desired low-redshift linear power spectrum to the required starting redshift consistently with the Newtonian evolution of the…
Initial conditions for cosmological N-body simulations are usually calculated by rescaling the present day linear power spectrum obtained from an Einstein-Boltzmann solver to the initial time employing the scale-independent matter growth…
We present a description for setting initial particle displacements and field values for simulations of arbitrary metric theories of gravity, for perfect and imperfect fluids with arbitrary characteristics. We extend the Zel'dovich…
We present a framework for general relativistic N-body simulations in the regime of weak gravitational fields. In this approach, Einstein's equations are expanded in terms of metric perturbations about a Friedmann-Lema\^itre background,…
We show how standard Newtonian N-body simulations can be interpreted in terms of the weak-field limit of general relativity by employing the recently developed Newtonian motion gauge. Our framework allows the inclusion of radiation…
Cosmological $N$-body simulations are the standard tool to study the emergence of the observed large-scale structure of the Universe. Such simulations usually solve for the gravitational dynamics of matter within the Newtonian…
We show how the Zel'dovich approximation and the second order displacement field of Lagrangian perturbation theory can be obtained from a general relativistic gradient expansion in \Lambda{}CDM cosmology. The displacement field arises as a…
Simulations have become an indispensable tool for accurate modelling of observables measured in galaxy surveys, but can be expensive if very large dynamic range in scale is required. We describe how to combine Lagrangian perturbation theory…
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…
Newtonian N-body simulations have been employed successfully over the past decades for the simulation of the cosmological large-scale structure. Such simulations usually ignore radiation perturbations (photons and massless neutrinos) and…
We study the dynamics of small inhomogeneities in an expanding universe collapsing to form bound structures using full solutions of the Einstein-Vlasov (N-body) equations. We compare these to standard Newtonian N-body solutions using…
We present $N$-body simulations in which either all, or a fraction of, the cold dark matter decays non-relativistically to a relativistic, non-interacting dark radiation component. All effects from radiation and general relativity are…
We present a novel method for including the impact of massive neutrinos in cold dark matter N-body simulations. Our approach is compatible with widely employed Newtonian N-body codes and relies on only three simple modifications. First, we…