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

DOI10.1088/1475-7516/2020/04/056MaRDI QIDQ5068477zbMATH OpenOpenAlexFDO

Authors Cristian Barrera Hinojosa, Baojiu Li

Publication date 6 April 2022

Published in Journal of Cosmology and Astroparticle Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/2001.07968

zbMATH Keywords

cosmological simulationsdark matter simulations

Mathematics Subject Classification ID

Nonlinear elliptic equations (35J60) Yang-Mills and other gauge theories in mechanics of particles and systems (70S15) Many-body theory; quantum Hall effect (81V70) Relativistic gravitational theories other than Einstein's, including asymmetric field theories (83D05) Relativistic cosmology (83F05) Lattice gravity, Regge calculus and other discrete methods in general relativity and gravitational theory (83C27) Propagation of singularities; initial value problems on manifolds (58J47) Flows related to mean curvature (53E10) Dark matter and dark energy (83C56) Mathematical modeling or simulation for problems pertaining to relativity and gravitational theory (83-10)

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.

Recommendations

• Relativistic initial conditions for N-body simulations
• General relativistic cosmological N-body simulations. Part I. Time integration
• A relativistic view on large scale \(N\)-body simulations
• \({N}\)-body methods for relativistic cosmology
• scientific article; zbMATH DE number 769955

Cites work

• A cosmologically motivated reference formulation of numerical relativity
• A new approach to cosmological structure formation with massive neutrinos
• General relativistic cosmological N-body simulations. Part I. Time integration
• General relativistic weak-field limit and Newtonian \(N\)-body simulations
• Improved constrained scheme for the Einstein equations: an approach to the uniqueness issue
• Numerical relativity. Solving Einstein's equations on the computer
• Relativistic initial conditions for N-body simulations
• WAVELESS APPROXIMATION THEORIES OF GRAVITY

Cited in

(6)

• IllinoisGRMHD: an open-source, user-friendly GRMHD code for dynamical spacetimes
• General relativistic cosmological N-body simulations. Part I. Time integration
• GRAMSES: a new route to general relativistic N-body simulations in cosmology. Part I. Methodology and code description
• Relativistic initial conditions for N-body simulations
• Fast full N-body simulations of generic modified gravity: conformal coupling models
• EBWeyl: a code to invariantly characterize numerical spacetimes

Describes a project that uses

Uses Software

• nbodykit
• gevolution
• GRAMSES
• CLASS
• DTFE

This page was built for publication: GRAMSES: a new route to general relativistic \(N\)-body simulations in cosmology. II: Initial conditions

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5068477)