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.

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

References in corpus (35)

• Dynamics of dark energy
• The Confrontation between General Relativity and Experiment
• Resummation of Massive Gravity
• Models of f(R) Cosmic Acceleration that Evade Solar-System Tests
• The DESI Experiment Part I: Science,Targeting, and Survey Design
• Covariant Galileon
• 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
• Non-linear Evolution of f(R) Cosmologies III: Halo Statistics
• The DESI Experiment Part II: Instrument Design
• Non-linear evolution of f(R) cosmologies II: power spectrum
• N-body Simulations for f(R) Gravity using a Self-adaptive Particle-Mesh Code
• Non-linear evolution of f(R) cosmologies I: methodology
• Effective Field Theory of Cosmic Acceleration: constraining dark energy with CMB data
• Observational constraints on Galileon cosmology
• Modified gravity a la Galileon: Late time cosmic acceleration and observational constraints
• A Compendium of Chameleon Constraints
• Linear perturbations in Galileon gravity models
• Constraining chameleon models with cosmology
• Accurate estimators of power spectra in N-body simulations
• The XMM Cluster Survey: Testing chameleon gravity using the profiles of clusters
• 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
• K-mouflage Cosmology: the Background Evolution
• The Vainshtein Mechanism in the Cosmic Web
• Testing the quasi-static approximation in $f(R)$ gravity simulations
• K-mouflage Cosmology: Formation of Large-Scale Structures
• A Fast Route to Non-Linear Clustering Statistics in Modified Gravity Theories
• Non-Gaussianity in single field models without slow-roll
• Hydrodynamic Effects in the Symmetron and $f(R)$-gravity Models
• Domain walls coupled to matter: the symmetron example
• Constraining f(R) Gravity with Planck Sunyaev-Zel'dovich Clusters
• Can background cosmology hold the key for modified gravity tests?

Cited by in corpus (8)

• Finite-time Cosmological Singularities and the Possible Fate of the Universe
• Realistic simulations of galaxy formation in f(R) modified gravity
• Redshift space power spectrum beyond Einstein-de Sitter kernels
• Convolution Lagrangian perturbation theory for biased tracers beyond general relativity
• Simulation techniques for modified gravity
• On the galaxy 3-point correlation function in Modified Gravity
• The Galaxy Halo Connection in Modified Gravity Cosmologies: Environment Dependence of Galaxy Luminosity function
• Universality of the halo mass function in modified gravity cosmologies