We use cookies to ensure that we give you the best experience on our website. By continuing to browse this repository, you give consent for essential cookies to be used. You can read more about our Privacy and Cookie Policy.

Durham Research Online
You are in:

GRAMSES : a new route to general relativistic $N$-body simulations in cosmology. Part II. Initial conditions.

Barrera-Hinojosa, Cristian and Li, Baojiu (2020) 'GRAMSES : a new route to general relativistic $N$-body simulations in cosmology. Part II. Initial conditions.', Journal of cosmology and astroparticle physics., 2020 (4). 056.


We address the generation of initial conditions (ICs) for {\sc gramses}, a code for nonlinear general relativistic (GR) N-body cosmological simulations recently introduced in ref. [1]. {\sc 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 {\sc 2LPTic} code, and demonstrate that in a Newtonian setting it can reproduce the velocity field given by the default {\sc 2LPTic} code with subpercent accuracy. We also show that the matter and velocity power spectra of the initial particle data generated for {\sc gramses} simulations using this method agree very well with the linear-theory predictions in the particular gauge used by {\sc gramses}. Finally, we discuss corrections to the finite difference calculation of the velocity when radiation is present, as well as additional corrections implemented in {\sc 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.

Item Type:Article
Full text:(AM) Accepted Manuscript
Download PDF
Publisher Web site:
Publisher statement:The deposited manuscript is available under a CC BY-NC-ND 4.0 licence.
Date accepted:06 April 2020
Date deposited:20 May 2020
Date of first online publication:30 April 2020
Date first made open access:30 April 2021

Save or Share this output

Look up in GoogleScholar