ﻻ يوجد ملخص باللغة العربية
We address the generation of initial conditions (ICs) for GRAMSES, a code for nonlinear general relativistic (GR) $N$-body cosmological simulations recently introduced in Ref. [1]. GRAMSES adopts a constant mean curvature slicing with a minimal distortion gauge, where the linear growth rate is scale-dependent, and the standard method for realising initial particle data is not straightforwardly applicable. A new method is introduced, in which the initial positions of particles are generated from the displacement field realised for a matter power spectrum as usual, but the velocity is calculated by finite-differencing the displacement fields around the initial redshift. In this way, all the information required for setting up the initial conditions is drawn from three consecutive input matter power spectra, and additional assumptions such as scale-independence of the linear growth factor and growth rate are not needed. We implement this method in a modified 2LPTic code, and demonstrate that in a Newtonian setting it can reproduce the velocity field given by the default 2LPTic code with subpercent accuracy. We also show that the matter and velocity power spectra of the initial particle data generated for GRAMSES simulations using this method agree very well with the linear-theory predictions in the particular gauge used by GRAMSES. Finally, we discuss corrections to the finite difference calculation of the velocity when radiation is present, as well as additional corrections implemented in GRAMSES to ensure consistency. This method can be applied in ICs generation for GR simulations in generic gauges, and simulations of cosmological models with scale-dependent linear growth rate.
We present GRAMSES, a new pipeline for nonlinear cosmological $N$-body simulations in General Relativity (GR). This code adopts the Arnowitt-Deser-Misner (ADM) formalism of GR, with constant mean curvature and minimum distortion gauge fixings, which
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 back-scaling m
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 perturbations and
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 Zeldovich Approximation to no
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 scales, coin