The e-MANTIS emulator: fast predictions of the non-linear matter power spectrum in $f(R)$CDM cosmology

Abstract

In order to probe modifications of gravity at cosmological scales, one needs accurate theoretical predictions. N-body simulations are required to explore the non-linear regime of structure formation but are very time consuming. In this work, we release a new public emulator, dubbed e-MANTIS, that performs an accurate and fast interpolation between the predictions of $f(R)$ modified gravity cosmological simulations, run with ECOSMOG. We sample a wide 3D parameter space given by the current background scalar field value $10^{-7}<\left|f_{R_0}\right|<10^{-4}$, matter density $0.24<{\Omega }_{\mathrm{m}}<0.39$, and primordial power spectrum normalisation $0.6<{\sigma }_8<1.0$, with 110 points sampled from a Latin Hypercube. For each model we perform pairs of $f(R)$CDM and $\Lambda$CDM simulations covering an effective volume of ${\left(560h^{-1}\mathrm{Mpc}\right)}^3$ with a mass resolution of $\sim 2\times 10^{10}{h}^{-1}{M}_{\odot }$. We build an emulator for the matter power spectrum boost $B(k)=P_{f(R)}(k)/P_{\Lambda \mathrm{CDM}}(k)$ using a Gaussian Process Regression method. The boost is mostly independent of $h$, $n_s$, and ${\Omega }_b$, which reduces the dimensionality of the relevant cosmological parameter space. Additionally, it is more robust against statistical and systematic errors than the raw power spectrum, thus strongly reducing our computational needs. According to our dedicated study of numerical systematics, the resulting emulator has an estimated maximum error of $3\%$ across the whole cosmological parameter space, for scales $0.03h{\mathrm{Mpc}}^{-1}<k<7h{\mathrm{Mpc}}^{-1}$, and redshifts $0<z<2$, while in most cases the accuracy is better than $1\%$. Such an emulator could be used to constrain $f(R)$ gravity with weak lensing analyses

23 pages, 18 figures. Accepted for publication in MNRAS. Link to the emulator code: https://zenodo.org/doi/10.5281/zenodo.7738362

References in corpus (24)

• Array Programming with NumPy
• Models of f(R) Cosmic Acceleration that Evade Solar-System Tests
• Transients from Initial Conditions in Cosmological Simulations
• Non-linear evolution of f(R) cosmologies I: methodology
• Initial Conditions for Large Cosmological Simulations
• A simple multigrid scheme for solving the Poisson equation with arbitrary domain boundaries
• Euclid preparation: IX. EuclidEmulator2 -- Power spectrum emulation with massive neutrinos and self-consistent dark energy perturbations
• The Mira-Titan Universe II: Matter Power Spectrum Emulation
• Cosmic Calibration: Constraints from the Matter Power Spectrum and the Cosmic Microwave Background
• Testing Screened Modified Gravity
• Emulation of baryonic effects on the matter power spectrum and constraints from galaxy cluster data
• Testing the quasi-static approximation in $f(R)$ gravity simulations
• Matter Power Spectrum Covariance Matrix from the DEUS-PUR ΛCDM simulations: Mass Resolution and non-Gaussian Errors
• Realistic simulations of galaxy formation in f(R) modified gravity
• Matter Power Spectrum Emulator for f(R) Modified Gravity Cosmologies
• Speeding up $N$-body simulations of modified gravity: Chameleon screening models
• The Mira-Titan Universe IV. High Precision Power Spectrum Emulation
• FORGE -- the f(R) gravity cosmic emulator project I: Introduction and matter power spectrum emulator
• Fast full N-body simulations of generic modified gravity: conformal coupling models
• Fast full $N$-body simulations of generic modified gravity: derivative coupling models
• Cosmological tests of modified gravity: constraints on $F(R)$ theories from the galaxy clustering ratio
• Simulation techniques for modified gravity
• The RayGalGroupSims cosmological simulation suite for the study of relativistic effects: an application to lensing-matter clustering statistics
• Fast production of cosmological emulators in modified gravity: the matter power spectrum

Cited by in corpus (16)

• CSST Cosmological Emulator I: Matter Power Spectrum Emulation with one percent accuracy to k = 10 h/Mpc
• Sesame: A power spectrum emulator pipeline for beyond-$\Lambda$CDM models
• Constraining f(R) gravity with cross-correlation of galaxies and cosmic microwave background lensing
• FREmu: Power Spectrum Emulator for $f(R)$ Gravity
• A field-level emulator for modified gravity
• Euclid preparation LXXI. Simulations and nonlinearities beyond $\mathsf{\Lambda }$CDM. 3. Constraints on $f(R)$ models from the photometric primary probes
• Principal Components for Model-Agnostic Modified Gravity with 3x2pt
• Cosmological gravity on all scales IV: 3x2pt Fisher forecasts for pixelised phenomenological modified gravity
• PySCo: A fast Particle-Mesh $N$-body code for modified gravity simulations in Python
• Power Spectrum Emulators from Neural Networks and Tree-Based Methods
• Probing gravity with non-linear clustering in redshift space
• Cosmological gravity on all scales V: MCMC forecasts combining large scale structure and CMB lensing for binned phenomenological modified gravity
• A simple prediction of the nonlinear matter power spectrum in Brans-Dicke gravity from linear theory
• Marked statistics across the cosmic web: Environmental dependent clustering in modified gravity simulations
• Cluster gravitational redshifts: uncertainties and survey requirements
• Toward testing gravity with LSST using $E_G$