Fast numerical method to generate halo catalogs in modified gravity (part I): second-order Lagrangian Perturbation Theory

Abstract

We present and test a new numerical method to determine second-order Lagrangian displacement fields in the context of modified gravity (MG) theories. We start from the extension of Lagrangian Perturbation Theory to a class of MG models that can be described by a parametrized Poisson equation, with the introduction of a scale-dependent function. We exploit fast Fourier transforms to compute the full source term of the differential equation for the second-order Lagrangian displacement field. We compare its mean to the source term computed for specific configurations for which a k-dependent solution can be found numerically. We choose the configuration that best matches the full source term, thus obtaining an approximate factorization of the second-order displacement field as the space term valid for standard gravity times a k-dependent, second-order growth factor $D_2(k,t)$. This approximation is used to compute second order displacements for particles. The method is tested against N-body simulations run with standard and $f(R)$ gravity: we rely on the results of a friends-of-friends code run on the N-body snapshots to assign particles to halos, then compute the halo power spectrum. We find very consistent results for the two gravity theories: second-order LPT (2LPT) allows to recover the halo power spectrum of N-body simulations within $\sim 10\%$ precision to $k\sim 0.2-0.4h{\mathrm{Mpc}}^{-1}$, as well as halo positions, with an error that is a fraction of the inter-particle distance. We show that, when considering the same level of non-linearity in the density field, the performance of 2LPT with MG is the same (within $1\%$) as the one obtained for the standard $\Lambda$CDM model with General Relativity. When implemented in a computer code, this formulation of 2LPT can quickly generate dark matter distributions with $f(R)$ gravity, and can easily be extended to other MG theories.

arXiv comment13 pages, 11 figures, submitted to MNRAS for review

References in corpus (38)

• Observational Evidence from Supernovae for an Accelerating Universe and a Cosmological Constant
• Measurements of Omega and Lambda from 42 High-Redshift Supernovae
• Observation of Gravitational Waves from a Binary Black Hole Merger
• Planck 2015 results. XIII. Cosmological parameters
• First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole
• f(R) theories
• Cosmological parameters from CMB and other data: a Monte-Carlo approach
• Models of f(R) Cosmic Acceleration that Evade Solar-System Tests
• Beyond the Cosmological Standard Model
• Cosmology and Fundamental Physics with the Euclid Satellite
• Everything You Always Wanted To Know About The Cosmological Constant Problem (But Were Afraid To Ask)
• Beyond $\Lambda$CDM: Problems, solutions, and the road ahead
• Testing General Relativity in Cosmology
• Solving Large Scale Structure in Ten Easy Steps with COLA
• PINOCCHIO: pinpointing orbit-crossing collapsed hierarchical objects in a linear density field
• Accurate Estimators of Correlation Functions in Fourier Space
• Cosmology with massive neutrinos II: on the universality of the halo mass function and bias
• Modified Gravity-GADGET: A new code for cosmological hydrodynamical simulations of modified gravity models
• Modelling Baryon Acoustic Oscillations with Perturbation Theory and Stochastic Halo Biasing
• Non-linear Evolution of Matter Power Spectrum in Modified Theory of Gravity
• Cosmic Degeneracies I: Joint N-body Simulations of Modified Gravity and Massive Neutrinos
• ICE-COLA: Towards fast and accurate synthetic galaxy catalogues optimizing a quasi $N$-body method
• The mass-Peak Patch algorithm for fast generation of deep all-sky dark matter halo catalogues and its N-Body validation
• COLA with scale-dependent growth: applications to screened modified gravity models
• Comparing approximate methods for mock catalogues and covariance matrices II: Power spectrum multipoles
• Comparing approximate methods for mock catalogues and covariance matrices I: correlation function
• Fast and accurate mock catalogue generation for low-mass galaxies
• Comparing approximate methods for mock catalogues and covariance matrices III: Bispectrum
• Improving fast generation of halo catalogs with higher-order Lagrangian perturbation theory
• Weak Lensing Light-Cones in Modified Gravity simulations with and without Massive Neutrinos
• HALOGEN: A tool for fast generation of mock halo catalogues
• Efficient simulations of large scale structure in modified gravity cosmologies with comoving Lagrangian acceleration
• EFTCAMB/EFTCosmoMC: massive neutrinos in dark cosmologies
• Lagrangian perturbation theory for modified gravity
• Investigating the degeneracy between modified gravity and massive neutrinos with redshift-space distortions
• Nonlinear evolution of initially biased tracers in modified gravity
• Simulating cosmologies beyond $\Lambda$CDM with PINOCCHIO
• Testing approximate predictions of displacements of cosmological dark matter halos

Cited by in corpus (7)

• A Lagrangian Perturbation Theory in the presence of massive neutrinos
• Fast computation of non-linear power spectrum in cosmologies with massive neutrinos
• $\text{Hi-COLA}$: Fast, approximate simulations of structure formation in Horndeski gravity
• Fast generation of mock galaxy catalogues in modified gravity models with COLA
• Non-linear density-velocity dynamics in $f(R)$ gravity from spherical collapse
• Numerical implementation of the Cubic Galileon model in PINOCCHIO
• An implementation of nDGP gravity in Pinocchio