Testing gravity on cosmic scales: A case study of Jordan-Brans-Dicke theory

We provide an end-to-end exploration of a distinct modified gravitational theory in Jordan-Brans-Dicke (JBD) gravity, from an analytical and numerical description of the background expansion and linear perturbations, to the nonlinear regime captured with a hybrid suite of N-body simulations, to the cosmological constraints from existing probes of the expansion history, the large-scale structure, and the cosmic microwave background (CMB). We have focused on JBD gravity as it both approximates a wider class of Horndeski scalar-tensor theories on cosmological scales and allows us to adequately model the nonlinear corrections to the matter power spectrum. In a combined analysis of the Planck 2018 CMB temperature, polarization, and lensing reconstruction, together with Pantheon supernova distances and the Baryon Oscillation Spectroscopic Survey (BOSS) measurements of baryon acoustic oscillation distances, the Alcock-Paczynski effect, and the growth rate, we constrain the JBD coupling constant to ₁₃>970 (95% confidence level; C. L. ) in agreement with the General Relativistic expectation given by ₁₃. In the unrestricted JBD model, where the effective gravitational constant at present, G₌₀ₓₓ₄ₑ/G, is additionally varied, increased dataset concordance (e. g. , within 1 agreement in S₈=₈{₌/0. 3}) enables us to further include the combined (``32pt'') dataset of cosmic shear, galaxy-galaxy lensing, and overlapping redshift-space galaxy clustering from the Kilo Degree Survey and the 2-degree Field Lensing Survey (KiDS2dFLenS). In analyzing the weak lensing measurements, the nonlinear corrections due to baryons, massive neutrinos, and modified gravity are simultaneously modeled and propagated in the cosmological analysis for the first time. In the joint analysis of all datasets, we constrain ₁₃>1540 (95% C. L. ), G₌₀ₓₓ₄ₑ/G=0. 9970. 029, the sum of neutrino masses, m_1460 and ₁₃>2230 (both at 95% C. L. ) in the restricted and unrestricted JBD models, respectively. We explore the impact of the JBD modeling choices, and show that a more restrictive parametrization of the coupling constant degrades the neutrino mass bound by up to a factor of three. In addition to the improved concordance between KiDS2dFLenS and Planck, the tension in the Hubble constant between Planck and the direct measurement of Riess et al. (2019) is reduced to 3; however, we find no substantial model selection preference for JBD gravity relative to. We further show that a positive shift in the effective gravitational constant suppresses the CMB damping tail, which might complicate future inferences of small-scale physics, given its degeneracy with the primordial helium abundance, the effective number of neutrinos, and the running of the spectral index.