Speeding up $N$-body simulations of modified gravity: Chameleon screening models

Abstract

We describe and demonstrate the potential of a new and very efficient method for simulating certain classes of modified gravity theories, such as the widely studied $f(R)$ gravity models. High resolution simulations for such models are currently very slow due to the highly nonlinear partial differential equation that needs to be solved exactly to predict the modified gravitational force. This nonlinearity is partly inherent, but is also exacerbated by the specific numerical algorithm used, which employs a variable redefinition to prevent numerical instabilities. The standard Newton-Gauss-Seidel iterative method used to tackle this problem has a poor convergence rate. Our new method not only avoids this, but also allows the discretised equation to be written in a form that is analytically solvable. We show that this new method greatly improves the performance and efficiency of $f(R)$ simulations. For example, a test simulation with $512^3$ particles in a box of size $512\mathrm{Mpc}/h$ is now 5 times faster than before, while a Millennium-resolution simulation for $f(R)$ gravity is estimated to be more than 20 times faster than with the old method. Our new implementation will be particularly useful for running very high resolution, large-sized simulations which, to date, are only possible for the standard model, and also makes it feasible to run large numbers of lower resolution simulations for covariance analyses. We hope that the method will bring us to a new era for precision cosmological tests of gravity.

arXiv comment14 pages, 3 figures, 1 table; version accepted for publication in JCAP

References in corpus (116)

• Planck 2015 results. XIII. Cosmological parameters
• Seven-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Interpretation
• Dynamics of dark energy
• Nine-Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: Cosmological Parameter Results
• Simulating the joint evolution of quasars, galaxies and their large-scale distribution
• f(R) Theories Of Gravity
• Modified Gravity and Cosmology
• The Confrontation between General Relativity and Experiment
• 4D Gravity on a Brane in 5D Minkowski Space
• The 6dF Galaxy Survey: Baryon Acoustic Oscillations and the Local Hubble Constant
• Cosmological Hydrodynamics with Adaptive Mesh Refinement: a new high resolution code called RAMSES
• The galileon as a local modification of gravity
• Baryon Acoustic Oscillations in the Sloan Digital Sky Survey Data Release 7 Galaxy Sample
• Resummation of Massive Gravity
• Models of f(R) Cosmic Acceleration that Evade Solar-System Tests
• Chameleon Cosmology
• The DESI Experiment Part I: Science,Targeting, and Survey Design
• Covariant Galileon
• A Redetermination of the Hubble Constant with the Hubble Space Telescope from a Differential Distance Ladder
• Cosmological Tests of Modified Gravity
• Symmetron Fields: Screening Long-Range Forces Through Local Symmetry Restoration
• The Supernova Legacy Survey 3-year sample: Type Ia Supernovae photometric distances and cosmological constraints
• Dark Energy vs. Modified Gravity
• The clustering of galaxies in the SDSS-III Baryon Oscillation Spectroscopic Survey: measurements of the growth of structure and expansion rate at z=0.57 from anisotropic clustering
• The DESI Experiment, a whitepaper for Snowmass 2013
• Evading Equivalence Principle Violations, Cosmological and other Experimental Constraints in Scalar Field Theories with a Strong Coupling to Matter
• The Cosmology of f(R) Gravity in the Metric Variational Approach
• f(R) Gravity and Chameleon Theories
• Galileon Cosmology
• Cosmology of a covariant Galileon field
• Non-linear Evolution of f(R) Cosmologies III: Halo Statistics
• ECOSMOG: An Efficient Code for Simulating Modified Gravity
• The DESI Experiment Part II: Instrument Design
• The Dilaton and Modified Gravity
• A Unified Description of Screened Modified Gravity
• Non-linear evolution of f(R) cosmologies II: power spectrum
• N-body Simulations for f(R) Gravity using a Self-adaptive Particle-Mesh Code
• Self-Accelerating Universe in Galileon Cosmology
• Astrophysical Tests of Modified Gravity: Constraints from Distance Indicators in the Nearby Universe
• Testing Gravity using Cosmic Voids
• Modified Gravity-GADGET: A new code for cosmological hydrodynamical simulations of modified gravity models
• Skewness and Kurtosis in Large-Scale Cosmic Fields
• Vainshtein screening in a cosmological background in the most general second-order scalar-tensor theory
• Non-linear evolution of f(R) cosmologies I: methodology
• Effective Field Theory of Cosmic Acceleration: constraining dark energy with CMB data
• Halos and Voids in f(R) Gravity
• Galileon gravity and its relevance to late time cosmic acceleration
• No-Go Theorems for Generalized Chameleon Field Theories
• Modified Gravity N-body Code Comparison Project
• New constraints on $f(R)$ gravity from clusters of galaxies
• Observational constraints on Galileon cosmology
• The nonlinear matter and velocity power spectra in f(R) gravity
• Skewness Induced by Gravity
• Nonlinear structure formation in the Cubic Galileon gravity model
• ISIS: a new N-body cosmological code with scalar fields based on RAMSES. Code presentation and application to the shapes of clusters
• Redshift space distortions in f(R) gravity
• Modified gravity a la Galileon: Late time cosmic acceleration and observational constraints
• A Compendium of Chameleon Constraints
• Linear perturbations in Galileon gravity models
• Testing chameleon gravity with the Coma cluster
• Structure Formation in the Symmetron Model
• Exploring Vainshtein mechanism on adaptively refined meshes
• Testing Gravity Using Dwarf Stars
• Systematic simulations of modified gravity: symmetron and dilaton models
• The Paths of Gravity in Galileon Cosmology
• Constraining chameleon models with cosmology
• Accurate estimators of power spectra in N-body simulations
• Vainshtein mechanism in second-order scalar-tensor theories
• Simulating the quartic Galileon gravity model on adaptively refined meshes
• The XMM Cluster Survey: Testing chameleon gravity using the profiles of clusters
• Testing General Relativity using the Environmental Dependence of Dark Matter Halos
• Revisiting fifth forces in the Galileon model
• Cosmological Structure Formation under MOND: a new numerical solver for Poisson's equation
• Constraining $f(R)$ Gravity Theory Using Weak Lensing Peak Statistics from the Canada-France-Hawaii-Telescope Lensing Survey
• A Convenient Set of Comoving Cosmological Variables and Their Application
• Vainshtein Mechanism In ${\Lambda }_3$ - Theories
• An Extended Excursion Set Approach to Structure Formation in Chameleon Models
• Laboratory Tests of the Galileon
• Modified Gravity Tomography
• Structure Formation by Fifth Force I: N-Body vs. Linear Simulations
• Chameleon f(R) gravity in the virialized cluster
• K-mouflage Cosmology: the Background Evolution
• Characterising Vainshtein Solutions in Massive Gravity
• Halo model and halo properties in Galileon gravity cosmologies
• Cosmic Web and Environmental Dependence of Screening: Vainshtein vs. Chameleon
• Systematic simulations of modified gravity: chameleon models
• Chameleon Halo Modeling in f(R) Gravity
• Galaxy Infall Kinematics as a Test of Modified Gravity
• Astrophysical Tests of Modified Gravity: the Morphology and Kinematics of Dwarf Galaxies
• The Vainshtein Mechanism in the Cosmic Web
• Modeling the nonlinear clustering in modified gravity models I: A fitting formula for matter power spectrum of f(R) gravity
• Cluster abundance in chameleon $f(R)$ gravity I: toward an accurate halo mass function prediction
• Galileons on Trial
• Testing the quasi-static approximation in $f(R)$ gravity simulations
• K-mouflage Cosmology: Formation of Large-Scale Structures
• Galileon Hairs of Dyson Spheres, Vainshtein's Coiffure and Hirsute Bubbles
• Structure Formation by the Fifth Force III: Segregation of Baryons and Dark Matter
• Constraints on f(R) theory and Galileons from the latest data of galaxy redshift surveys
• Experimental constraints on the uncoupled Galileon model from SNLS3 data and other cosmological probes
• Hierarchical clustering in chameleon $f(R)$ gravity
• A Fast Route to Non-Linear Clustering Statistics in Modified Gravity Theories
• Exploring the liminality: properties of haloes and subhaloes in borderline $f(R)$ gravity
• Speeding up N-body simulations of modified gravity: Vainshtein screening models
• Non-Gaussianity in single field models without slow-roll
• Hydrodynamic Effects in the Symmetron and $f(R)$-gravity Models
• Structure Formation by Fifth Force: Power Spectrum from N-Body Simulations
• Simulations of Galileon modified gravity: Clustering statistics in real and redshift space
• Rapid simulation rescaling from standard to modified gravity models
• Simulation tests of galaxy cluster constraints on chameleon gravity
• Domain walls coupled to matter: the symmetron example
• Weak lensing by galaxy troughs with modified gravity
• Modified Gravity Spins Up Galactic Halos
• Zoomed cosmological simulations of Milky Way sized halos in f(R)-gravity
• Constraining f(R) Gravity with Planck Sunyaev-Zel'dovich Clusters
• Local Group velocity versus gravity: the coherence function
• Can background cosmology hold the key for modified gravity tests?

Cited by in corpus (41)

• Dark energy in Horndeski theories after GW170817: A review
• The Santiago-Harvard-Edinburgh-Durham void comparison I: SHEDding light on chameleon gravity tests
• On the road to percent accuracy: nonlinear reaction of the matter power spectrum to dark energy and modified gravity
• COLA with scale-dependent growth: applications to screened modified gravity models
• Finite-time Cosmological Singularities and the Possible Fate of the Universe
• The Novel Probes Project -- Tests of Gravity on Astrophysical Scales
• Testing the theory of gravity with DESI: estimators, predictions and simulation requirements
• The Santiago-Harvard-Edinburgh-Durham void comparison II: unveiling the Vainshtein screening using weak lensing
• Matter Power Spectrum Emulator for f(R) Modified Gravity Cosmologies
• Realistic simulations of galaxy formation in f(R) modified gravity
• Growth of Cosmic Structure
• Large-scale redshift space distortions in modified gravity theories
• FORGE -- the f(R) gravity cosmic emulator project I: Introduction and matter power spectrum emulator
• Marked clustering statistics in $f(R)$ gravity cosmologies
• Fast full N-body simulations of generic modified gravity: conformal coupling models
• An accurate perturbative approach to redshift space clustering of biased tracers in modified gravity
• Fast full $N$-body simulations of generic modified gravity: derivative coupling models
• Simulating galaxy formation in f(R) modified gravity: Matter, halo, and galaxy-statistics
• Redshift space power spectrum beyond Einstein-de Sitter kernels
• The shape dependence of chameleon screening
• Convolution Lagrangian perturbation theory for biased tracers beyond general relativity
• Real- and redshift-space halo clustering in $f(R)$ cosmologies
• The e-MANTIS emulator: fast predictions of the non-linear matter power spectrum in $f(R)$CDM cosmology
• Simulation techniques for modified gravity
• Testing $f(R)$ Gravity With Scale Dependent Cosmic Void Velocity Profiles
• Cosmological Probes of Structure Growth and Tests of Gravity
• Fifth-Force Screening around Extremely Compact Sources
• Probing gravity with redshift-space distortions: effects of tracer bias and sample selection
• Accurate Computation of the Screening of Scalar Fifth Forces in Galaxies
• Robustness of the pairwise kinematic Sunyaev-Zeldovich power spectrum shape as a cosmological gravity probe
• Improved analytical modeling of the non-linear power spectrum in modified gravity cosmologies
• N-body simulations, halo mass functions, and halo density profile in $f(T)$ gravity
• Screening maps of the local Universe I -- Methodology
• On the galaxy 3-point correlation function in Modified Gravity
• Euclid preparation LXXI. Simulations and nonlinearities beyond $\mathsf{\Lambda }$CDM. 3. Constraints on $f(R)$ models from the photometric primary probes
• Euclid preparation. Simulations and nonlinearities beyond $\Lambda$CDM. 1. Numerical methods and validation
• The Galaxy Halo Connection in Modified Gravity Cosmologies: Environment Dependence of Galaxy Luminosity function
• Breaking the ${\tau }_{\mathrm{T}}-f$ degeneracy of the kinetic Sunyaev Zel'dovich cosmology in redshift space
• PySCo: A fast Particle-Mesh $N$-body code for modified gravity simulations in Python
• Bending the web: exploring the impact of modified gravity on the density field and halo properties within the cosmic web
• Universality of the halo mass function in modified gravity cosmologies