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

Abstract

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.

27 pages, 4 figures. Major revisions of the Newtonian-relativistic correspondence, main results unchanged. Matches accepted version

References in corpus (12)

• Transients from Initial Conditions in Cosmological Simulations
• The Cosmic Linear Anisotropy Solving System (CLASS) I: Overview
• Initial Conditions for Large Cosmological Simulations
• A constrained scheme for Einstein equations based on Dirac gauge and spherical coordinates
• Improved constrained scheme for the Einstein equations: An approach to the uniqueness issue
• Waveless Approximation Theories of Gravity
• Accurate estimators of power spectra in N-body simulations
• Inhomogeneous Cosmology with Numerical Relativity
• Mathematical Issues in a Fully-Constrained Formulation of Einstein Equations
• The Limited Accuracy of Linearized Gravity
• A new approach to cosmological structure formation with massive neutrinos
• General Relativistic Cosmological N-body Simulations I: time integration

Cited by in corpus (12)

• In the Realm of the Hubble tension $-$ a Review of Solutions
• Fast full N-body simulations of generic modified gravity: conformal coupling models
• Higher-order initial conditions with massive neutrinos
• Luminosity distance and anisotropic sky-sampling at low redshifts: a numerical relativity study
• Numerical solutions to Einstein's equations in a shearing-dust Universe: a code comparison
• Cosmology using numerical relativity
• Vector modes in $\Lambda$CDM: the gravitomagnetic potential in dark matter haloes from relativistic $N$-body simulations
• Looking for a twist: probing the cosmological gravitomagnetic effect via weak lensing-kSZ cross correlations
• First investigation of void statistics in numerical relativity simulations
• EBWeyl: a Code to Invariantly Characterize Numerical Spacetimes
• Structure formation and quasi-spherical collapse from initial curvature perturbations with numerical relativity simulations
• Measuring the gravitomagnetic distortion from rotating halos I: methods