A simple prediction of the nonlinear matter power spectrum in Brans-Dicke gravity from linear theory

Abstract

Brans-Dicke (BD), one of the first proposed scalar-tensor theories of gravity, effectively makes the gravitational constant of general relativity (GR) time-dependent. Constraints on the BD parameter $\omega$ serve as a benchmark for testing GR, which is recovered in the limit $\omega \to \infty$. Current small-scale astrophysical constraints $\omega \gtrsim 10^5$ are much tighter than large-scale cosmological constraints $\omega \gtrsim 10^3$, but the two decouple if the true theory of gravity features screening. On the largest cosmological scales, BD approximates the most general second-order scalar-tensor (Horndeski) theory, so constraints here have wider implications. These constraints will improve with upcoming large-scale structure and cosmic microwave background surveys. To constrain BD with weak gravitational lensing, one needs its nonlinear matter power spectrum $P_{\mathrm{BD}}$. By comparing the boost $B=P_{\mathrm{BD}}/P_{\mathrm{GR}}$ from linear theory and nonlinear $N$-body simulations, we show that the nonlinear boost can simply be predicted from linear theory if the BD and GR universes are parameterized in a way that makes their early cosmological evolution and quasilinear power today similar. In particular, they need the same $H_0/\sqrt{G_{\mathrm{eff}}(a=0)}$ and ${\sigma }_8$, where $G_{\mathrm{eff}}$ is the (effective) gravitational strength. Our prediction is $1\%$ accurate for $\omega \ge 100$, $z\le 3$, and $k\le 1h/\mathrm{Mpc}$; and $2\%$ up to $k\le 5h/\mathrm{Mpc}$. It also holds for $G_{\mathrm{BD}}$ that do not match Newton's constant today, so one can study GR with different gravitational constants $G_{\mathrm{GR}}$ by sending $\omega \to \infty$. We provide a code that computes $B$ with the linear Einstein-Boltzmann solver hi_class and multiplies it by the nonlinear $P_{\mathrm{GR}}$ from EuclidEmulator2 to predict $P_{\mathrm{BD}}$.

12 pages, 8 figures. Accepted for publication in A&A. A code that computes the non-linear boost can be found at https://github.com/hersle/jbd

References in corpus (20)

• The Confrontation between General Relativity and Experiment
• Dark Energy after GW170817: dead ends and the road ahead
• The relativistic pulsar-white dwarf binary PSR J1738+0333 II. The most stringent test of scalar-tensor gravity
• HMcode-2020: Improved modelling of non-linear cosmological power spectra with baryonic feedback
• 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
• Cosmological performance of SKA HI galaxy surveys
• On the road to percent accuracy IV: ReACT -- computing the non-linear power spectrum beyond $\Lambda$CDM
• An improved test of the strong equivalence principle with the pulsar in a triple star system
• A comparison of Einstein-Boltzmann solvers for testing General Relativity
• The Mira-Titan Universe IV. High Precision Power Spectrum Emulation
• Testing Gravity on Cosmic Scales: A Case Study of Jordan-Brans-Dicke Theory
• Accelerating Large-Scale-Structure data analyses by emulating Boltzmann solvers and Lagrangian Perturbation Theory
• Cosmological constraints on the gravitational constant
• Enabling matter power spectrum emulation in beyond-$\Lambda$CDM cosmologies with COLA
• The e-MANTIS emulator: fast predictions of the non-linear matter power spectrum in $f(R)$CDM cosmology
• Testing extended Jordan-Brans-Dicke theories with future cosmological observations
• Euclid: Constraining linearly scale-independent modifications of gravity with the spectroscopic and photometric primary probes
• Sesame: A power spectrum emulator pipeline for beyond-$\Lambda$CDM models
• Tracking the validity of the quasi-static and sub-horizon approximations in modified gravity