GrGadget: an N-body TreePM relativistic code for cosmological simulations

Abstract

We present the merging of the Particle-Mesh (PM) relativistic Gevolution code with the TreePM Gadget-4 code, with the aim of studying general relativity effects in cosmology. Our code, called GrGadget, is able to track the evolution of metric perturbations in the weak field limit by using Gevolution's implementation of a relativistic PM in the Poisson gauge. To achieve this, starting from Gevolution we have written a C++ library called libgevolution, that allows a code to access and use the same abstractions and resources that Gevolution uses for its PM-only N-body simulations. The code works under the assumption that particle interactions at short distances can be approximated as Newtonian, so that we can combine the forces computed with a Newtonian Tree with those computed with a relativistic PM. The result is a TreePM simulation code that represents metric perturbations at the scales where they are relevant, while resolving non-linear structures. We validate our code by closely matching Gadget-4 forces, computed with the Tree switched off, with those computed with libgevolution in the Newtonian limit. With GrGadget we obtain a matter power spectrum that is compatible with Newtonian Gadget at small scales and contains GR features at large scales that are consistent with results obtained with Gevolution. We demonstrate that, due to the better resolution of the highly non-linear regime, the representation of the relativistic fields sampled on the mesh improves with respect to the PM-only simulations.

arXiv comment15 pages, 14 figures

References in corpus (30)

• Observation of Gravitational Waves from a Binary Black Hole Merger
• The cosmological simulation code GADGET-2
• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• LSST: from Science Drivers to Reference Design and Anticipated Data Products
• Weak Gravitational Lensing
• The DESI Experiment Part I: Science,Targeting, and Survey Design
• Wide-Field InfrarRed Survey Telescope-Astrophysics Focused Telescope Assets WFIRST-AFTA 2015 Report
• Simulating cosmic structure formation with the GADGET-4 code
• Approaches to Understanding Cosmic Acceleration
• ECOSMOG: An Efficient Code for Simulating Modified Gravity
• Accurate Estimators of Correlation Functions in Fourier Space
• Modified Gravity-GADGET: A new code for cosmological hydrodynamical simulations of modified gravity models
• gevolution: a cosmological N-body code based on General Relativity
• Dark Energy Survey Year 1 Results: Multi-Probe Methodology and Simulated Likelihood Analyses
• ISIS: a new N-body cosmological code with scalar fields based on RAMSES. Code presentation and application to the shapes of clusters
• Connection between Newtonian simulations and general relativity
• Testing the theory of gravity with DESI: estimators, predictions and simulation requirements
• Relativistic N-body simulations with massive neutrinos
• Computing General Relativistic effects from Newtonian N-body simulations: Frame dragging in the post-Friedmann approach
• N-body methods for relativistic cosmology
• LIGER: mock relativistic light-cones from Newtonian simulations
• $k$-evolution: a relativistic N-body code for clustering dark energy
• Weak-lensing observables in relativistic N-body simulations
• GRAMSES: a new route to general relativistic $N$-body simulations in cosmology. Part I. Methodology and code description
• $N$-body simulations for parametrised modified gravity
• Galaxy formation in the brane world I: overview and first results
• Clustering dark energy imprints on cosmological observables of the gravitational field
• $\text{fRevolution}$ $-$ Relativistic Cosmological Simulations in $f(R)$ Gravity I: Methodology
• Latfield2: A c++ library for classical lattice field theory
• Vector modes in $\Lambda$CDM: the gravitomagnetic potential in dark matter haloes from relativistic $N$-body simulations

Cited by in corpus (3)

• A novel approach to cosmological non-linearities as an effective fluid
• Structure formation and quasi-spherical collapse from initial curvature perturbations with numerical relativity simulations
• Local environmental dependence on weak-lensing shear statistics