Despite the great observational success of the standard cosmological model some discrepancies in the inferred parameter constraints have manifested among a number of cosmological data sets. These include a tension between the expansion rate of our Cosmos as inferred from the cosmic microwave background (CMB) and as found from local measurements, the preference for an enhanced amplitude of CMB lensing, a somewhat low quadrupole moment of the CMB fluctuations as well as a preference for a lower amplitude of matter fluctuations in large-scale structure surveys than inferred from the CMB. We analyse these observational tensions under the addition of spatial curvature and a free CMB background temperature that may deviate from its locally measured value. With inclusion of these parameters, we observe a trend in the parameter constraints from cosmic microwave background and baryon acoustic oscillation data towards an open and hotter universe with larger current expansion rate, standard CMB lensing amplitudes, lower amplitude of matter fluctuations, and lower CMB quadrupole moment, consistently reducing the individual tensions among the cosmological data sets. Combining this data, we find a preference for an open and hotter universe. Finally, we briefly discuss a local void as a possible source for a deviation of the locally measured CMB temperature from its background value and as mimic of negative spatial curvature for CMB photons, which when correctly implemented may further ease parameter tensions.

We bring together two popular formalisms which generically parameterise deviations from General Relativity on astrophysical and cosmological scales, namely the parameterised post-Newtonian (PPN) formalism and the effective field theory (EFT) of dark energy and modified gravity. These separate formalisms are successfully applied to independently perform tests of gravity in their respective regimes of applicability on vastly different length scales. Nonlinear screening mechanisms indeed make it imperative to probe General Relativity across a wide range of scales. For a comprehensive interpretation of the complementary measurements it is important to connect them to effectively constrain the vast gravitational model space. We establish such a connection within the framework of Horndeski scalar-tensor theories restricted to a luminal propagation speed of gravitational waves. This is possible via the reconstruction of the family of linearly degenerate covariant Horndeski actions from the set of EFT functions and the subsequent derivation of the PPN parameters from the reconstructed theory. We outline the required conditions which ensure a reconstructed Horndeski model possesses a screening mechanism that enables significant modifications on cosmological scales while respecting stringent astrophysical bounds. Employing a scaling method, we then perform the general post-Newtonian expansion of the reconstructed models to derive their PPN parameters $\gamma$ and $\beta$ in their screened regimes.

Outside galaxy clusters the competition between the inwards gravitational attraction and the outwards expansion of the Universe leads to a special radius of velocity cancellation, which is called the turn-around radius. Measurements of the turn-around radius hold promises of constraining cosmological parameters, and possibly even properties of gravity. Such a measurement is, however, complicated by the fact that the surroundings of galaxy clusters are not spherical, but instead are a complicated collection of filaments, sheets and voids. In this paper we use the results of numerically simulated universes to quantify realistic error-bars of the measurement of the turn-around radius. We find that for a $\Lambda$CDM cosmology these error-bars are typically of the order of $20\%$. We numerically simulate three different implementations of dark energy models and of a scalar dark sector interaction to address whether the turn-around radius can be used to constrain non-trivial cosmologies, and we find that only rather extreme models can be distinguished from a $\Lambda$CDM universe due to the large error-bars arising from the non-trivial cluster environments.

To effectively exploit large-scale structure surveys, we depend on accurate and reliable predictions of non-linear cosmological structure formation. Tools for efficient and comprehensive computational modelling are therefore essential to perform cosmological parameter inference analyses. We present the public software package ReACT, demonstrating its capability for the fast and accurate calculation of non-linear power spectra from non-standard physics. We showcase ReACT through a series of forecast constraints on the DGP and $f(R)$ gravity models, adopting LSST-like cosmic shear power spectra. Accurate non-linear modelling with ReACT more than doubles LSST's constraining power on the $f(R)$ parameter, in contrast to an analysis that is limited to the quasi-linear regime. We find that ReACT is sufficiently robust for the inference of consistent constraints on theories beyond $\Lambda$CDM for current and ongoing surveys. With further improvement, particularly in terms of the accuracy of the non-linear $\Lambda$CDM power spectrum, ReACT can, in principle, meet the accuracy requirements for future surveys such as Euclid and LSST.

We present $\texttt{MG-evolution}$, an $N$-body code simulating the cosmological structure formation for parametrised modifications of gravity. It is built from the combination of parametrised linear theory with a parametrisation of the deeply nonlinear cosmological regime extrapolated from modified spherical collapse computations that cover the range of known screening mechanisms. We test $\texttt {MG-evolution}$, which runs at the speed of conventional $\Lambda$CDM simulations, against a suit of existing exact model-specific codes, encompassing linearised and chameleon $f(R)$ gravity as well as the normal branch of the Dvali-Gabadadz-Porrati braneworld model, hence covering both large-field value and large-derivative screening effects. We compare the nonlinear power spectra produced by the parametrised and model-specific approaches over the full range of scales set by the box size and resolution of our simulations, $k=(0.05-2.5)$~h/Mpc, and for two redshift slices, $z=0$ and $z=1$. We find sub-percent to one-percent level recovery of all the power spectra generated with the model-specific codes for the full range of scales. $\texttt {MG-evolution}$ can be used for generalised and accurate tests of gravity and dark energy with the increasing wealth of high-precision cosmological survey data becoming available over the next decade.

Positivity bounds provide conditions that a consistent UV-completion exists for a quantum field theory. We examine their application to Horndeski gravity models reconstructed from the effective field theory (EFT) of dark energy. This enables us to assess whether particular phenomenological parameterizations of the EFT functions reconstruct theories that respect or violate the positivity bounds. We find that commonly adopted EFT parametrizations, cast in terms of the dark energy density or power laws of the scale factor, only satisfy the positivity bounds in non-trivial regions of the parameter space. We then examine parameterizations of the inherently stable EFT basis, constructed to avoid gradient and ghost instabilities by default. In stark contrast, in this basis the positivity bounds either only provide constraints in a-priori unrealistic regions of the parameter space or do not provide any constraints on parameter values at all. The application of positivity bounds to common parametrizations of the standard EFT functions can therefore lead to artificial conclusions that the region of viable Horndeski modifications of gravity is highly constrained. Our results provide a strong motivation, in addition to the default avoidance of theoretical instabilities, for instead adopting parametrizations of the inherently stable EFT basis when testing dark energy and modified gravity models with forthcoming cosmological survey data.

Recently the global variation of the Planck mass in the General Relativistic Einstein-Hilbert action was proposed as a self-tuning mechanism of the cosmological constant preventing vacuum energy from freely gravitating. We show that this global mechanism emerges for generic local scalar-tensor theories with additional coupling of the scalar field to the field strength of a three-form gauge field that turns the scalar field constant on the domain of the action. Evaluation of the resulting integral constraint equation over the observable Universe yields a self-consistent framework with General Relativistic field equations and arbitrary radiatively stable residual cosmological constant. We argue that the expectation value for this residual is in good agreement with the magnitude of the observed cosmic acceleration.

A significant tension has become manifest between the current expansion rate of our Universe measured from the cosmic microwave background by the Planck satellite and from local distance probes, which has prompted for interpretations of that as evidence of new physics. Within conventional cosmology a likely source of this discrepancy is identified here as a matter density fluctuation around the cosmic average of the 40 Mpc environment in which the calibration of Supernovae Type Ia separations with Cepheids and nearby absolute distance anchors is performed. Inhomogeneities on this scale easily reach 40% and more. In that context, the discrepant expansion rates serve as evidence of residing in an underdense region of $\delta_{\rm env}\approx-0.5\pm0.1$. The probability for finding this local expansion rate given the Planck data lies at the 95% confidence level. Likewise, a hypothetical equivalent local data set with mean expansion rate equal to that of Planck, while statistically favoured, would not gain strong preference over the actual data in the respective Bayes factor. These results therefore suggest borderline consistency between the local and Planck measurements of the Hubble constant. Generally accounting for the environmental uncertainty, the local measurement may be reinterpreted as a constraint on the cosmological Hubble constant of $H_0=74.7^{+5.8}_{-4.2}$ km/s/Mpc. The current simplified analysis may be augmented with the employment of the full available data sets, an impact study for the immediate $\lesssim10$ Mpc environment of the distance anchors, more prone to inhomogeneities, as well as expansion rates measured by quasar lensing, gravitational waves, currently limited to the same 40 Mpc region, and local galaxy distributions.

Future large-scale structure surveys will measure three-point statistics with high statistical significance. This will offer significant improvements on our understanding of gravity, provided we can model these statistics accurately. We assess the performance of several schemes for theoretical modelling of the matter bispectrum, including halo-model based approaches and fitting formulae. We compare the model predictions against N-body simulations, considering scales up to $k_{\rm max} = 4 h/{\rm Mpc}$, well into non-linear regime of structure formation. Focusing on the equilateral configuration, we conduct this analysis for three theories of gravity: general relativity, $f(R)$ gravity, and the DGP braneworld model. Additionally, we compute the lensing convergence bispectrum for these models. We find that all current modelling prescriptions in modified gravity, in particular for theories with scale-dependent linear growth, fail to attain the accuracy required by the precision of the Stage IV surveys such as \emph{Euclid}. Among these models, we find that a halo-model corrected fitting formula achieves the best overall performance.

We revisit the dark degeneracy that arises from the Einstein equations relating geometry to the total cosmic substratum but not resolving its individual components separately. We establish the explicit conditions for the dark degeneracy in the fluid description of the dark sector. At the background level, this degeneracy can be formally understood in terms of a unified dark sector Equation of State (EoS) that depends both on the dynamical nature of the dark energy (DE) as well as on its interaction with the pressureless dark matter. For linear perturbations, the degeneracy arises for specified DE pressure perturbations (or sound speed, equivalently) and DE anisotropic stress. Specializing to the degeneracy between non-interacting dynamical DE and interacting vacuum DE models, we perform a parameter estimation analysis for a range of dynamical DE parametrizations, where for illustration we explicitly break the degeneracy at the linear level by adopting a luminal sound speed for both scenarios. We conduct this analysis using cosmological background data alone and in combination with Planck measurements of the cosmic microwave background radiation. We find that although the overall phenomenology between the dynamical DE and interacting approaches is similar, there are some intriguing differences. In particular, there is an ambiguity in the strength of constraints on $\Omega_{m0}$ and $\sigma_8$, which are considerably weakened for interacting vacuum DE, indicating that the dark degeneracy can change the significance of tensions in cosmological parameters inferred from different data sets.

Modified gravity theories with an effective Newton constant that varies over cosmological timescales generally predict a different gravitational wave luminosity distance than General Relativity. While this holds for a uniform variation, we show that if locally screened at the source and at the observer as required to pass stringent astrophysical tests of gravity, the General Relativistic distance is restored. In the absence of such a screening, the same effect must modify electromagnetic luminosity distances inferred from supernovae Type Ia, to the extent that the effects can cancel in the comparison. Hence, either the modifications considered employ screening, which leaves no signature in Standard Sirens of a cosmological modification of gravity, or screening does not operate, in which case there can be a signal that is however well below the forseeable sensitivity of the probe when astrophysical bounds are employed. We recover these results both in the Jordan and Einstein frames, paying acute attention to pecularities of each frame such as the notion of redshift or geodesic motions. We emphasise that despite these limitations, Standard Sirens provide valuable independent tests of gravity that differ fundamentally from other probes, a circumstance that is generally important for the wider scope of gravitational modifications and related scenarios. Finally, we use our results to show that the gravitational wave propagation is not affected by dark sector interactions, which restores a dark degeneracy between conformal and disformal couplings that enables observationally viable cosmic self-acceleration to emenate from those.

Standard sirens have been proposed as probes of alternative theories of gravity, such as Horndeski models. Hitherto, all studies have been conducted on a homogeneous-isotropic cosmological background, which is unable to consistently account for realistic distributions of matter, and for inhomogeneities in the Horndeski scalar field. Yet, the latter are essential for screening mechanisms. In this article, we comprehensively analyze the propagation of Horndeski gravitational waves in an arbitrary background spacetime and scalar field. We restrict to the class of theories in which gravitational waves propagate at light speed, and work in the geometric-optics regime. We find that kinetic braiding only produces a nonphysical longitudinal mode, whereas conformal coupling affects the amplitude of the standard transverse modes but not their polarization. We confirm that any observable deviation from general relativity depends on the local value of the effective Planck mass at emission and reception of the wave. This result is interpreted as the conservation of the number of gravitons.

Accurate predictions for the non-linear matter power spectrum are needed to confront theory with observations in current and near future weak lensing and galaxy clustering surveys. We propose a computationally cheap method to create an emulator for modified gravity models by utilizing existing emulators for $\Lambda$CDM. Using a suite of $N$-body simulations we construct a fitting function for the enhancement of both the linear and non-linear matter power spectrum in the commonly studied Hu-Sawicki $f(R)$ gravity model valid for wave-numbers $k \lesssim 5-10\, h\text{Mpc}^{-1}$ and redshifts $z \lesssim 3$. We show that the cosmology dependence of this enhancement is relatively weak so that our fit, using simulations coming from only one cosmology, can be used to get accurate predictions for other cosmological parameters. We also show that the cosmology dependence can, if needed, be included by using linear theory, approximate $N$-body simulations (such as COLA) and semi-analytical tools like the halo model. Our final fit can easily be combined with any emulator or semi-analytical models for the non-linear $\Lambda$CDM power spectrum to accurately, and quickly, produce a non-linear power spectrum for this particular modified gravity model. The method we use can be applied to fairly cheaply construct an emulator for other modified gravity models. As an application of our fitting formula we use it to compute Fisher-forecasts for how well galaxy clustering and weak lensing in a Euclid-like survey will be at constraining modifications of gravity.

Constraints on gravity and cosmology will greatly benefit from performing joint clustering and weak lensing analyses on large-scale structure data sets. Utilising non-linear information coming from small physical scales can greatly enhance these constraints. At the heart of these analyses is the matter power spectrum. Here we employ a simple method, dubbed "Hybrid $P_\ell(k)$", based on the Gaussian streaming model to calculate the quasi non-linear redshift space matter power spectrum multipoles from a non-linear prescription for the matter power spectrum, which remains agnostic about gravity or dark energy. We test this approach against comoving Lagrangian acceleration simulation measurements performed in GR, DGP and $f(R)$ gravity and find that our method performs comparably or better to the dark matter TNS redshift space power spectrum model. When comparing the redshift space multipoles for halos, we find that a linear bias implementation is insufficient to be competitive with the TNS-based model used in the BOSS clustering analysis, which however has additional degrees of freedom. Our approach offers many avenues for improvement in accuracy as well as further unification under the halo model.

An additional variation of the Einstein-Hilbert action with respect to the Planck mass provides a constraint on the average Ricci scalar that prevents vacuum energy from gravitating. Consideration of the evolution of the inhomogeneous matter distribution in the Universe with evaluation of the averaging constraint on disconnected matter cells that ultimately form isolated gravitationally bound structures yields a backreaction effect that self-consistently produces the cosmological constant of the background. A uniform prior on our location in the formation of these isolated structures implies a mean expectation for the present cosmological constant energy density parameter of $\Omega_{\Lambda}=0.704$, giving rise to a late-time acceleration of the cosmic expansion and a coincident current energy density of matter.

With the increasing wealth of high-quality astronomical and cosmological data and the manifold departures from General Relativity in principle conceivable, the development of generalized parametrization frameworks that unify gravitational models and cover a wide range of length scales and a variety of observational probes to enable systematic high-precision tests of gravity has been a stimulus for intensive research. A review is presented here for some of the formalisms devised for this purpose, covering the cosmological large- and small-scale structures, the astronomical static weak-field regime as well as emission and propagation effects for gravitational waves. This includes linear and nonlinear parametrized post-Friedmannian frameworks, effective field theory approaches, the parametrized post-Newtonian expansion, the parametrized post-Einsteinian formalism as well as an inspiral-merger-ringdown waveform model among others. Connections between the different formalisms are highlighted where they have been established and a brief outlook is provided for general steps towards a unified global framework for tests of gravity and dark sector models.

We introduce The Novel Probes Project, an initiative to advance the field of astrophysical tests of the dark sector by creating a forum that connects observers and theorists. This review focuses on tests of gravity and is intended to be of use primarily to observers, but also to theorists with interest in the development of experimental tests. It is twinned with a separate review on tests of dark matter self-interactions (Adhikari et al., in prep.). Our focus is on astrophysical probes of gravity in the weak-field regime, ranging from stars to quasilinear cosmological scales. These are complementary to both strong-field tests and background and linear probes in cosmology. In particular, the nonlinear screening mechanisms that are an integral part of viable modified gravity models lead to characteristic signals specifically on astrophysical scales. The constraining power of these signals is not limited by cosmic variance, but comes with the challenge of building robust theoretical models of the nonlinear dynamics of stars, galaxies, clusters and large scale structure. In this review we lay the groundwork for a thorough exploration of the astrophysical regime with an eye to using the current and next generation of observations for tests of gravity. We begin by setting the scene for how theories beyond General Relativity are expected to behave, focusing primarily on screened fifth forces. We describe the analytic and numerical techniques for exploring the pertinent astrophysical systems, as well as the signatures of modified gravity. With these in hand we present a range of observational tests, and discuss prospects for future measurements and theoretical developments.

We present a general method to compute the nonlinear matter power spectrum for dark energy and modified gravity scenarios with percent-level accuracy. By adopting the halo model and nonlinear perturbation theory, we predict the reaction of a $\Lambda$CDM matter power spectrum to the physics of an extended cosmological parameter space. By comparing our predictions to $N$-body simulations we demonstrate that with no-free parameters we can recover the nonlinear matter power spectrum for a wide range of different $w_0$-$w_a$ dark energy models to better than 1% accuracy out to $k \approx 1 \, h \, {\rm Mpc}^{-1}$. We obtain a similar performance for both DGP and $f(R)$ gravity, with the nonlinear matter power spectrum predicted to better than 3% accuracy over the same range of scales. When including direct measurements of the halo mass function from the simulations, this accuracy improves to 1%. With a single suite of standard $\Lambda$CDM $N$-body simulations, our methodology provides a direct route to constrain a wide range of non-standard extensions to the concordance cosmology in the high signal-to-noise nonlinear regime.

Modifications of General Relativity leave their imprint both on the cosmic expansion history through a non-trivial dark energy equation of state, and on the evolution of cosmological perturbations in the scalar and in the tensor sectors. In particular, the modification in the tensor sector gives rise to a notion of gravitational-wave (GW) luminosity distance, different from the standard electromagnetic luminosity distance, that can be studied with standard sirens at GW detectors such as LISA or third-generation ground based experiments. We discuss the predictions for modified GW propagation from some of the best studied theories of modified gravity, such as Horndeski or the more general degenerate higher order scalar-tensor (DHOST) theories, non-local infrared modifications of gravity, bigravity theories and the corresponding phenomenon of GW oscillation, as well as theories with extra or varying dimensions. We show that modified GW propagation is a completely generic phenomenon in modified gravity. We then use a simple parametrization of the effect in terms of two parameters $(\Xi_0,n)$, that is shown to fit well the results from a large class of models, to study the prospects of observing modified GW propagation using supermassive black hole binaries as standard sirens with LISA. We construct mock source catalogs and perform detailed Markov Chain Monte Carlo studies of the likelihood obtained from LISA standard sirens alone, as well as by combining them with CMB, BAO and SNe data to reduce the degeneracies between cosmological parameters. We find that the combination of LISA with the other cosmological datasets allows one to measure the parameter $\Xi_0$ that characterizes modified GW propagation to the percent level accuracy, sufficient to test several modified gravity theories. [Abridged]

We have previously presented a reconstruction of Horndeski scalar-tensor theories from linear cosmological observables. It includes free nonlinear terms which can be added onto the reconstructed covariant theory without affecting the background and linear dynamics. After discussing the uniqueness of these correction terms, we apply this nonlinear freedom to a range of different applications. First we demonstrate how the correction terms can be configured to endow the reconstructed models with screening mechanisms such as the chameleon, k-mouflage and Vainshtein effects. A further implication is the existence of classes of Horndeski models that are degenerate with standard cosmology to an arbitrary level in the cosmological perturbations. Particularly interesting examples are kinetically self-accelerating models that mimic the dynamics of the cosmological constant to an arbitrary degree in perturbations. Finally, we develop the reconstruction method further to the level of higher-order effective field theory, which under the restriction to a luminal propagation speed of gravitational waves introduces two new free functions per order. These functions determine the corresponding correction terms in the linearly reconstructed action at the same order. Our results enable the connection of linear cosmological constraints on generalised modifications of gravity and dark energy with the nonlinear regime and astrophysical probes for a more global interpretation of the wealth of forthcoming cosmological survey data.

The sequestering mechanism has been proposed as a remedy to the old cosmological constant problem of the non-gravitating vacuum energy in the matter sector. Here it is shown that an extension of this global constraint mechanism arises naturally from an averaging condition for residual cosmological constants produced in different patches of the Universe. The new mechanism naturally yields the small observed value that gives rise to cosmic acceleration at late times, hence, addressing the new cosmological constant problem. Hereby the halo model picture is adopted with all matter content ultimately residing in the last collapsed structures formed. Sequestering in these collapsing patches produces the desired average residual, where a uniform prior on our present location in the collapse predicts $\Omega_{\Lambda}=0.697$. Finally, a fluctuation of the local residual from the cosmological average can naturally give rise to a locally enhanced Hubble constant.

The growing wealth of cosmological observations places increasingly more stringent constraints on dark energy and alternative gravity models. Particularly successful in efficiently probing the vast model space has been the effective field theory of dark energy and modified gravity, providing a unified framework for generalised cosmological predictions. However, the optimal parametrisation of the free time-dependent functions inherent to the formalism is still unresolved. It should respect a multitude of requirements, ranging from simplicity, generality, and representativity of known theories to computational efficiency. But in particular, for theoretical soundness, the parameter space should adhere to strict stability requirements. We have recently proposed an inherently stable effective field theory with physical basis of Planck mass evolution, sound speed of the scalar field fluctuation, kinetic coefficient, and background expansion, which covers Horndeski models with luminal speed of gravity. Here we devise a parametrisation of these basis functions that can straightforwardly be configured to evade theoretical pathologies such as ghost or gradient instabilities or to accommodate further theoretical priors such as (sub)luminal scalar sound speed. The parametrisation is simple yet general, conveniently represents a broad range of known dark energy and gravitational theories, and with a simple additional approximation can be rendered numerically highly efficient. Finally, by operating in our new basis, we show that there are no general limitations from stability requirements on the current values that can be assumed by the phenomenological modification of the Poisson equation and the gravitational slip besides the exclusion of anti-gravity. The inherently stable effective field theory is ready for implementation in parameter estimation analyses employing linear cosmological observations.

Recently we have derived a set of mapping relations that enables the reconstruction of the family of Horndeski scalar-tensor theories which reproduce the background dynamics and linear perturbations of a given set of effective field theory of dark energy coefficients. In this paper we present a number of applications of this reconstruction. We examine the form of the underlying theories behind different phenomenological parameterizations of modified gravity and dark energy used in the literature, as well as examine theories that exhibit weak gravity, linear shielding, and minimal self-acceleration. Finally, we propose a new inherently stable parametrization basis for modified gravity and dark energy models.

Studying the effects of dark energy and modified gravity on cosmological scales has led to a great number of physical models being developed. The effective field theory (EFT) of cosmic acceleration allows an efficient exploration of this large model space, usually carried out on a phenomenological basis. However, constraints on such parametrized EFT coefficients cannot be trivially connected to fundamental covariant theories. In this paper we reconstruct the class of covariant Horndeski scalar-tensor theories that reproduce the same background dynamics and linear perturbations as a given EFT action. One can use this reconstruction to interpret constraints on parametrized EFT coefficients in terms of viable covariant Horndeski theories. We demonstrate this method with a number of well-known models and discuss a range of future applications.

The cosmological coincidences between the matter and radiation energy densities at recombination as well as between the densities of matter and the cosmological constant at present time are well known. We point out that moreover the third intersection between the energy densities of radiation and the cosmological constant coincides with the reionization epoch. To quantify the statistical relevance of this concurrence, we compute the Bayes factor between the concordance cosmology with free Thomson scattering optical depth and a model for which this parameter is inferred from imposing a match between the time of density equality and the epoch of reionization. This is to characterize the potential explanatory gain if one were to find a parameter-free physical connection. We find a very strong preference for such a concurrence on the Jeffreys scale from current cosmological observations. We furthermore discuss the effect of choice of priors, changes in reionization history, and free sum of neutrino masses. We also estimate the impact of adding intermediate polarization data from the Planck High Frequency Instrument and prospects for future 21 cm surveys. In the first case, preference for the correlation remains substantial, whereas future data may give results more decisive in pro or substantial in contra. Finally, we provide a discussion on different interpretations of these findings. In particular, we show how a connection between the star-formation history and the cosmological background dynamics can give rise to this concurrence.

The parameterised post-Newtonian (PPN) formalism has enabled stringent tests of static weak-field gravity in a theory-independent manner. Here we incorporate screening mechanisms of modified gravity theories into the framework by introducing an effective gravitational coupling and defining the PPN parameters as functions of position. To determine these functions we develop a general method for efficiently performing the post-Newtonian expansion in screened regimes. For illustration, we derive all the PPN functions for a cubic galileon and a chameleon model. We also analyse the Shapiro time delay effect for these two models and find no deviations from General Relativity insofar as the signal path and the perturbing mass reside in a screened region of space.

We refine the mass and environment dependent spherical collapse model of chameleon $f(R)$ gravity by calibrating a phenomenological correction inspired by the parameterized post-Friedmann framework against high-resolution $N$-body simulations. We employ our method to predict the corresponding modified halo mass function, and provide fitting formulas to calculate the fractional enhancement of the $f(R)$ halo abundance with respect to that of General Relativity (GR) within a precision of $\lesssim 5\%$ from the results obtained in the simulations. Similar accuracy can be achieved for the full $f(R)$ mass function on the condition that the modeling of the reference GR abundance of halos is accurate at the percent level. We use our fits to forecast constraints on the additional scalar degree of freedom of the theory, finding that upper bounds competitive with current Solar System tests are within reach of cluster number count analyses from ongoing and upcoming surveys at much larger scales. Importantly, the flexibility of our method allows also for this to be applied to other scalar-tensor theories characterized by a mass and environment dependent spherical collapse.

Most of the information on our cosmos stems from either late-time observations or the imprint of early-time inhomogeneities on the cosmic microwave background. We explore to what extent early modifications of gravity, which become significant after recombination but then decay towards the present, can be constrained by current cosmological observations. For the evolution of the gravitational modification, we adopt the decaying mode of a hybrid-metric Palatini $f(\mathcal{R})$ gravity model which is designed to reproduce the standard cosmological background expansion history and due to the decay of the modification is naturally compatible with Solar-System tests. We embed the model in the effective field theory description of Horndeski scalar-tensor gravity with an early-time decoupling of the gravitational modification. Since the quasistatic approximation for the perturbations in the model breaks down at high redshifts, where modifications remain relevant, we introduce a computationally efficient correction to describe the evolution of the scalar field fluctuation in this regime. We compare the decaying early-time modification against geometric probes and recent Planck measurements and find no evidence for such effects in the observations. Current data constrains the scalar field value at $|f_{\mathcal{R}}(z=z_{\rm on})| \lesssim 10^{-2}$ for modifications introduced at redshifts $z_{\rm on}\sim(500-1000)$ with present-day value $|f_{\mathcal{R}0}|\lesssim10^{-8}$. Finally, we comment on constraints that will be achievable with future 21~cm surveys and gravitational wave experiments.

Viable modifications of gravity on cosmological scales predominantly rely on screening mechanisms to recover Einstein's Theory of General Relativity in the Solar System, where it has been well tested. A parametrisation of the effects of such modifications in the spherical collapse model is presented here for the use of modelling the modified nonlinear cosmological structure. The formalism allows an embedding of the different screening mechanisms operating in scalar-tensor theories through large values of the gravitational potential or its first or second derivatives as well as of linear suppression effects or more general transitions between modified and Einstein gravity limits. Each screening or suppression mechanism is parametrised by a time, mass, and environment dependent screening scale, an effective modified gravitational coupling in the fully unscreened limit that can be matched to linear theory, the exponent of a power-law radial profile of the screened coupling, determined by derivatives, symmetries, and potentials in the scalar field equation, and an interpolation rate between the screened and unscreened limits. Along with generalised perturbative methods, the parametrisation may be used to formulate a nonlinear extension to the linear parametrised post-Friedmannian framework to enable generalised tests of gravity with the wealth of observations from the nonlinear cosmological regime.

With the advent of gravitational-wave astronomy marked by the aLIGO GW150914 and GW151226 observations, a measurement of the cosmological speed of gravity will likely soon be realized. We show that a confirmation of equality to the speed of light as indicated by indirect Galactic observations will have important consequences for a very large class of alternative explanations of the late-time accelerated expansion of our Universe. It will break the dark degeneracy of self-accelerated Horndeski scalar-tensor theories in the large-scale structure that currently limits a rigorous discrimination between acceleration from modified gravity and from a cosmological constant or dark energy. Signatures of a self-acceleration must then manifest in the linear, unscreened cosmological structure. We describe the minimal modification required for self-acceleration with standard gravitational-wave speed and show that its maximum likelihood yields a 3-sigma poorer fit to cosmological observations compared to a cosmological constant. Hence, equality between the speeds challenges the concept of cosmic acceleration from a genuine scalar-tensor modification of gravity.

Understanding the reason for the observed accelerated expansion of the Universe represents one of the fundamental open questions in physics. In cosmology, a classification has emerged among physical models for the acceleration, distinguishing between Dark Energy and Modified Gravity. In this review, we give a brief overview of models in both categories as well as their phenomenology and characteristic observable signatures in cosmology. We also introduce a rigorous distinction between Dark Energy and Modified Gravity based on the strong and weak equivalence principles.

The Horndeski action is the most general scalar-tensor theory with at most second-order derivatives in the equations of motion, thus evading Ostrogradsky instabilities and making it of interest when modifying gravity at large scales. To pass local tests of gravity, these modifications predominantly rely on nonlinear screening mechanisms that recover Einstein's Theory of General Relativity in regions of high density. We derive a set of conditions on the four free functions of the Horndeski action that examine whether a specific model embedded in the action possesses an Einstein gravity limit or not. For this purpose, we develop a new and surprisingly simple scaling method that identifies dominant terms in the equations of motion by considering formal limits of the couplings that enter through the new terms in the modified action. This enables us to find regimes where nonlinear terms dominate and Einstein's field equations are recovered to leading order. Together with an efficient approximation of the scalar field profile, one can then further evaluate whether these limits can be attributed to a genuine screening effect. For illustration, we apply the analysis to both a cubic galileon and a chameleon model as well as to Brans-Dicke theory. Finally, we emphasise that the scaling method also provides a natural approach for performing post-Newtonian expansions in screened regimes.

We present an accurate non-linear matter power spectrum prediction scheme for a variety of extensions to the standard cosmological paradigm, which uses the tuned halo model previously developed in Mead (2015b). We consider dark energy models that are both minimally and non-minimally coupled, massive neutrinos and modified gravitational forces with chameleon and Vainshtein screening mechanisms. In all cases we compare halo-model power spectra to measurements from high-resolution simulations. We show that the tuned halo model method can predict the non-linear matter power spectrum measured from simulations of parameterised $w(a)$ dark energy models at the few per cent level for $k<10\,h\mathrm{Mpc}^{-1}$, and we present theoretically motivated extensions to cover non-minimally coupled scalar fields, massive neutrinos and Vainshtein screened modified gravity models that result in few per cent accurate power spectra for $k<10\,h\mathrm{Mpc}^{-1}$. For chameleon screened models we achieve only 10 per cent accuracy for the same range of scales. Finally, we use our halo model to investigate degeneracies between different extensions to the standard cosmological model, finding that the impact of baryonic feedback on the non-linear matter power spectrum can be considered independently of modified gravity or massive neutrino extensions. In contrast, considering the impact of modified gravity and massive neutrinos independently results in biased estimates of power at the level of 5 per cent at scales $k>0.5\,h\mathrm{Mpc}^{-1}$. An updated version of our publicly available HMcode can be found at https://github.com/alexander-mead/HMcode

We identify a scalar-tensor model embedded in the Horndeski action whose cosmological background and linear scalar fluctuations are degenerate with the concordance cosmology. The model admits a self-accelerated background expansion at late times that is stable against perturbations with a sound speed attributed to the new field that is equal to the speed of light. While degenerate in scalar fluctuations, self-acceleration of the model implies a present cosmological tensor mode propagation at < 95% of the speed of light with a damping of the wave amplitude that is > 5% less efficient than in general relativity. We show that these discrepancies are endemic to self-accelerated Horndeski theories with degenerate large-scale structure and are tested with measurements of gravitational waves emitted by events at cosmological distances. Hence, gravitational-wave cosmology breaks the dark degeneracy in observations of the large-scale structure between two fundamentally different explanations of cosmic acceleration - a cosmological constant and a scalar-tensor modification of gravity. The gravitational wave event GW150914 recently detected with the aLIGO instruments and its potential association with a weak short gamma-ray burst observed with the Fermi GBM experiment may have provided this crucial measurement.

Viable modifications of gravity that may produce cosmic acceleration need to be screened in high-density regions such as the Solar System, where general relativity is well tested. Screening mechanisms also prevent strong anomalies in the large-scale structure and limit the constraints that can be inferred on these gravity models from cosmology. We find that by suppressing the contribution of the screened high-density regions in the matter power spectrum, allowing a greater contribution of unscreened low densities, modified gravity models can be more readily discriminated from the concordance cosmology. Moreover, by variation of density thresholds, degeneracies with other effects may be dealt with more adequately. Specializing to chameleon gravity as a worked example for screening in modified gravity, employing N-body simulations of f(R) models and the halo model of chameleon theories, we demonstrate the effectiveness of this method. We find that a percent-level measurement of the clipped power at k < 0.3 h/Mpc can yield constraints on chameleon models that are more stringent than what is inferred from Solar System tests or distance indicators in unscreened dwarf galaxies. Finally, we verify that our method is also applicable to the Vainshtein mechanism.

Linear cosmological perturbations of a large class of modified gravity and dark energy models can be unified in the effective field theory of cosmic acceleration, encompassing Horndeski scalar-tensor theories and beyond. The fully available model space inherent to this formalism cannot be constrained by measurements in the quasistatic small-scale regime alone. To facilitate the analysis of modifications from the concordance model beyond this limit, we introduce a semi-dynamical treatment extrapolated from the evolution of perturbations at a pivot scale of choice. At small scales, and for Horndeski theories, the resulting modifications recover a quasistatic approximation but account for corrections to it near the Hubble scale. For models beyond Horndeski gravity, we find that the velocity field and time derivative of the spatial metric potential can generally not be neglected, even in the small-scale limit. We test the semi-dynamical approximation against the linear perturbations of a range of dark energy and modified gravity models, finding good agreement between the two.

We develop and test an algorithm to rescale a simulated dark-matter particle distribution or halo catalogue from a standard gravity model to that of a modified gravity model. This method is based on that of Angulo & White but with some additional ingredients to account for (i) scale-dependent growth of linear density perturbations and (ii) screening mechanisms that are generic features of viable modified gravity models. We attempt to keep the method as general as possible, so that it may plausibly be applied to a wide range of modified theories, although tests against simulations are restricted to a subclass of $f(R)$ models at this stage. We show that rescaling allows the power spectrum of matter to be reproduced at the $\sim 3$ per cent level in both real and redshift space up to $k=0.1h\,\mathrm{Mpc}^{-1}$ if we change the box size and alter the particle displacement field; this limit can be extended to $k=1h\,\mathrm{Mpc}^{-1}$ if we additionally alter halo internal structure. We simultaneously develop an algorithm that can be applied directly to a halo catalogue, in which case the halo mass function and clustering can be reproduced at the $\sim 5$ per cent level. Finally we investigate the clustering of halo particle distributions, generated from rescaled halo catalogues, and find that a similar accuracy can be reached.

We study the model space generated by the time-dependent operator coefficients in the effective field theory of the cosmological background evolution and perturbations of modified gravity and dark energy models. We identify three classes of modified gravity models that reduce to Newtonian gravity on the small scales of linear theory. These general classes contain enough freedom to simultaneously admit a matching of the concordance model background expansion history. In particular, there exists a large model space that mimics the concordance model on all linear quasistatic subhorizon scales as well as in the background evolution. Such models also exist when restricting the theory space to operators introduced in Horndeski scalar-tensor gravity. We emphasize that whereas the partially shielded scenarios might be of interest to study in connection with tensions between large and small scale data, with conventional cosmological probes, the ability to distinguish the fully shielded scenarios from the concordance model on near-horizon scales will remain limited by cosmic variance. Novel tests of the large-scale structure remedying this deficiency and accounting for the full covariant nature of the alternative gravitational theories, however, might yield further insights on gravity in this regime.

Chameleon fields may modify gravity on cluster scales while recovering general relativity locally. This article reviews signatures of chameleon modifications in the nonlinear cosmological structure, comparing different techniques to model them, summarising the current state of observational constraints, and concluding with an outlook on prospective constraints from future observations and applications of the analytic tools developed in the process to more general scalar-tensor theories. Particular focus is given to the Hu-Sawicki and designer models of f(R) gravity.

We analyse modelling techniques for the large-scale structure formed in scalar-tensor theories of constant Brans-Dicke parameter which match the concordance model background expansion history and produce a chameleon suppression of the gravitational modification in high-density regions. Thereby, we use a mass and environment dependent chameleon spherical collapse model, the Sheth-Tormen halo mass function and linear halo bias, the Navarro-Frenk-White halo density profile, and the halo model. Furthermore, using the spherical collapse model, we extrapolate a chameleon mass-concentration scaling relation from a LCDM prescription calibrated to N-body simulations. We also provide constraints on the model parameters to ensure viability on local scales. We test our description of the halo mass function and nonlinear matter power spectrum against the respective observables extracted from large-volume and high-resolution N-body simulations in the limiting case of f(R) gravity, corresponding to a vanishing Brans-Dicke parameter. We find good agreement between the two; the halo model provides a good qualitative description of the shape of the relative enhancement of the f(R) matter power spectrum with respect to LCDM caused by the extra attractive gravitational force but fails to recover the correct amplitude. Introducing an effective linear power spectrum in the computation of the two-halo term to account for an underestimation of the chameleon suppression at intermediate scales in our approach, we accurately reproduce the measurements from the N-body simulations.

The presence of a scalar field that couples nonminimally and universally to matter can enhance gravitational forces on cosmological scales while restoring general relativity in the Solar neighborhood. In the intermediate regime, kinematically inferred masses experience an additional radial dependence with respect to the underlying distribution of matter, which is caused by the increment of gravitational forces with increasing distance from the Milky Way center. The same effect can influence the internal kinematics of subhalos and cause cuspy matter distributions to appear core-like. Specializing to the chameleon model as a worked example, we demonstrate this effect by tracing the scalar field from the outskirts of the Milky Way halo to its interior, simultaneously fitting observed velocity dispersions of chemo-dynamically discriminated red giant populations in the Fornax and Sculptor dwarf spheroidals. Whereas in standard gravity these observations suggest that the matter distribution of the dwarfs is cored, we find that in the presence of a chameleon field the assumption of a cuspy Navarro-Frenk-White profile becomes perfectly compatible with the data. Importantly, chameleon models also predict the existence of slopes between two stellar subcomponents that in Newtonian gravity would be interpreted as a depletion of matter in the dwarf center. Hence, an observation of such an apparently pathological scenario may serve as a smoking gun for the presence of a chameleon field or a similar modification of gravity, independent of baryonic feedback effects. In general, measuring the dynamic mass profiles of the Milky Way dwarfs provides stronger constraints than those inferred from the screening scale of the Solar System since these are located at greater distances from the halo center.

We investigate the performance of semi-analytical modelling of large-scale structure in Galileon gravity cosmologies using results from N-body simulations. We focus on the Cubic and Quartic Galileon models that provide a reasonable fit to CMB, SNIa and BAO data. We demonstrate that the Sheth-Tormen mass function and linear halo bias can be calibrated to provide a very good fit to our simulation results. We also find that the halo concentration-mass relation is well fitted by a power law. The nonlinear matter power spectrum computed in the halo model approach is found to be inaccurate in the mildly nonlinear regime, but captures reasonably well the effects of the Vainshtein screening mechanism on small scales. In the Cubic model, the screening mechanism hides essentially all of the effects of the fifth force inside haloes. In the case of the Quartic model, the screening mechanism leaves behind residual modifications to gravity, which make the effective gravitational strength time-varying and smaller than the standard value. Compared to normal gravity, this causes a deficiency of massive haloes and leads to a weaker matter clustering on small scales. For both models, we show that there are realistic halo occupation distributions of Luminous Red Galaxies that can match both the observed large-scale clustering amplitude and the number density of these galaxies.

We propose a novel method to test the gravitational interactions in the outskirts of galaxy clusters. When gravity is modified, this is typically accompanied by the introduction of an additional scalar degree of freedom, which mediates an attractive fifth force. The presence of an extra gravitational coupling, however, is tightly constrained by local measurements. In chameleon modifications of gravity, local tests can be evaded by employing a screening mechanism that suppresses the fifth force in dense environments. While the chameleon field may be screened in the interior of the cluster, its outer region can still be affected by the extra force, introducing a deviation between the hydrostatic and lensing mass of the cluster. Thus, the chameleon modification can be tested by combining the gas and lensing measurements of the cluster. We demonstrate the operability of our method with the Coma cluster, for which both a lensing measurement and gas observations from the X-ray surface brightness, the X-ray temperature, and the Sunyaev-Zel'dovich effect are available. Using the joint observational data set, we perform a Markov chain Monte Carlo analysis of the parameter space describing the different profiles in both the Newtonian and chameleon scenarios. We report competitive constraints on the chameleon field amplitude and its coupling strength to matter. In the case of f(R) gravity, corresponding to a specific choice of the coupling, we find an upper bound on the background field amplitude of |f_{R0}|<6*10^{-5}, which is currently the tightest constraint on cosmological scales.

Current constraints on f(R) gravity from the large-scale structure are at the verge of penetrating into a region where the modified forces become nonlinearly suppressed. For a consistent treatment of observables at these scales, we study cluster quantities produced in chameleon and linearized Hu-Sawicki f(R) gravity dark matter N-body simulations. We find that the standard Navarro-Frenk-White halo density profile and the radial power law for the pseudo phase-space density provide equally good fits for f(R) clusters as they do in the Newtonian scenario. We give qualitative arguments for why this should be the case. For practical applications, we derive analytic relations, e.g., for the f(R) scalar field, the gravitational potential, and the velocity dispersion as seen within the virialized clusters. These functions are based on three degrees of freedom fitted to simulations, i.e., the characteristic density, scale, and velocity dispersion. We further analyze predictions for these fitting parameters from the gravitational collapse and the Jeans equation, which are found to agree well with the simulations. Our analytic results can be used to consistently constrain chameleon f(R) gravity with future observations on virialized cluster scales without the necessity of running a large number of simulations.

We explore the signatures of quintessence and modified gravity theories in the relativistic description of galaxy clustering within a parametrized post-Friedmann framework. For this purpose, we develop a calibration method to consistently account for horizon-scale effects in the linear parametrized Post-Friedmann perturbations of minimally and nonminimally coupled scalar-tensor theories and test it against the full model-specific fluctuations. We further study the relativistic effects in galaxy clustering for the normal and self-accelerating branches of the Dvali-Gabadadze-Porrati braneworld model as well as for phenomenological modifications of gravity. We quantify the impact of modified gravity and dark energy models on galaxy clustering by computing the velocity-to-matter density ratio F, the velocity contribution R, and the potential contribution P and give an estimate of their detectability in future galaxy surveys. Our results show that, in general, the relativistic correction contains additional information on gravity and dark energy, which needs to be taken into account in consistent horizon-scale tests of departures from LCDM using the galaxy-density field.

On cosmological scales, observations of the cluster abundance currently place the strongest constraints on f(R) gravity. These constraints lie in the large-field limit, where the modifications of general relativity can correctly be modeled by setting the Compton wavelength of the scalar field to its background value. These bounds are, however, at the verge of penetrating into a regime, where the modifications become nonlinearly suppressed due to the chameleon mechanism and cannot be described by this linearized approximation. For future constraints based on observations subjected to cluster abundance, it is therefore essential to consistently model the chameleon effect. We analyze descriptions of the halo mass function in chameleon f(R) gravity using a mass- and environment-dependent spherical collapse model in combination with excursion set theory and phenomenological fits to N-body simulations in the LCDM and f(R) gravity scenarios. Our halo mass functions consistently incorporate the chameleon suppression and cosmological parameter dependencies, improving upon previous formalisms and providing an important extension to N-body simulations for the application in consistent tests of gravity with observables sensitive to the abundance of clusters.

We use the joint measurement of geometry and growth from anisotropic galaxy clustering in the Baryon Oscillation Spectroscopic Survey Data Release 9 CMASS sample reported by Reid et al. to constrain dark energy properties and possible deviations from the General Relativity. Assuming GR and taking a prior on the linear matter power spectrum at high redshift from the cosmic microwave background (CMB), anisotropic clustering of the CMASS DR9 galaxies alone constrains $\Omega_{\rm m} = 0.308 \pm 0.022$ and $100\Omega_{\rm k} = 5.9 \pm 4.8$ for $w = -1$, or $w = -0.91 \pm 0.12$ for $\Omega_k = 0$. When combined with the full CMB likelihood, the addition of the anisotropic clustering measurements to the spherically-averaged BAO location increases the constraining power on dark energy by a factor of 4 in a flat CDM cosmology with constant dark energy equation of state $w$ (giving $w = -0.87 \pm 0.05$). This impressive gain depends on our measurement of both the growth of structure and Alcock-Paczynski effect, and is not realised when marginalising over the amplitude of redshift space distortions. Combining with both the CMB and Supernovae Type Ia (SNeIa), we find $\Omega_{\rm m} = 0.281 \pm 0.014$ and $1000\Omega_{\rm k}=-9.2\pm5.0$ for $w = -1$, or $w_0 = -1.13 \pm 0.12$ and $w_{\rm a}=0.65 \pm 0.36$ assuming $\Omega_k = 0$. Finally, when a $\Lambda$CDM background expansion is assumed, the combination of our estimate of the growth rate with previous growth measurements provides tight constraints on the parameters describing possible deviations from GR giving $\gamma = 0.64 \pm 0.05$. For one parameter extensions of the flat $\Lambda$CDM model, we find a $\sim 2\sigma$ preference either for $w > -1$ or slower growth than in GR. However, the data is fully consistent with the concordance model, and the evidence for these additional parameters is weaker than $2\sigma$.

We present a new test of gravitational interactions at the r\sim(0.2-20)Mpc scale, around the virial radius of dark matter halos measured through cluster-galaxy lensing of maxBCG clusters from the Sloan Digital Sky Survey (SDSS). We employ predictions from self-consistent simulations of f(R) gravity to find an upper bound on the background field amplitude of f_R0<3.5x10^-3 at the 1D-marginalized 95% confidence level. As a model-independent assessment of the constraining power of cluster profiles measured through weak gravitational lensing, we also constrain the amplitude F_0 of a phenomenological modification based on the profile enhancement induced by f(R) gravity when not including effects from the increased cluster abundance in f(R). In both scenarios, dark-matter-only simulations of the concordance model corresponding to f_R0=0 and F_0=0 are consistent with the lensing measurements, i.e., at the 68% and 95% confidence level, respectively.

The standard model of cosmology LCDM assumes general relativity, flat space, and the presence of a positive cosmological constant. We relax these assumptions allowing spatial curvature, time-dependent effective dark energy equation of state, as well as modifications of the Poisson equation for the lensing potential, and modifications of the growth of linear matter density perturbations in alternate combinations. Using six parameters characterizing these relations, we check LCDM for consistency utilizing cosmic microwave background anisotropies, cross correlations thereof with high-redshift galaxies through the integrated Sachs-Wolfe effect, the Hubble constant, supernovae and baryon acoustic oscillation distances, as well as the relation between weak gravitational lensing and galaxy flows. In all scenarios, we find consistency of the concordance model at the 95% confidence level. However, we emphasize that constraining supplementary background parameters and parametrizations of the growth of large-scale structure separately may lead to a priori exclusion of viable departures from the concordance model.

We study cosmological constraints on metric f(R) gravity models that are designed to reproduce the LCDM expansion history with modifications to gravity described by a supplementary cosmological freedom, the Compton wavelength parameter B_0. We conduct a Markov chain Monte Carlo analysis on the parameter space, utilizing the geometrical constraints from supernovae distances, the baryon acoustic oscillation distances, and the Hubble constant, along with all of the cosmic microwave background data, including the largest scales, its correlation with galaxies, and a probe of the relation between weak gravitational lensing and galaxy flows. The strongest constraints, however, are obtained through the inclusion of data from cluster abundance. Using all of the data, we infer a bound of B_0<0.0011 at the 95% C.L.

Deviations from general relativity, such as could be responsible for the cosmic acceleration, would influence the growth of large scale structure and the deflection of light by that structure. We clarify the relations between several different model independent approaches to deviations from general relativity appearing in the literature, devising a translation table. We examine current constraints on such deviations, using weak gravitational lensing data of the CFHTLS and COSMOS surveys, cosmic microwave background radiation data of WMAP5, and supernova distance data of Union2. Markov Chain Monte Carlo likelihood analysis of the parameters over various redshift ranges yields consistency with general relativity at the 95% confidence level.

Although general relativity underlies modern cosmology, its applicability on cosmological length scales has yet to be stringently tested. Such a test has recently been proposed, using a quantity, EG, that combines measures of large-scale gravitational lensing, galaxy clustering and structure growth rate. The combination is insensitive to 'galaxy bias' (the difference between the clustering of visible galaxies and invisible dark matter) and is thus robust to the uncertainty in this parameter. Modified theories of gravity generally predict values of EG different from the general relativistic prediction because, in these theories, the 'gravitational slip' (the difference between the two potentials that describe perturbations in the gravitational metric) is non-zero, which leads to changes in the growth of structure and the strength of the gravitational lensing effect3. Here we report that EG = 0.39 +/- 0.06 on length scales of tens of megaparsecs, in agreement with the general relativistic prediction of EG $\approx$ 0.4. The measured value excludes a model within the tensor-vector-scalar gravity theory, which modifies both Newtonian and Einstein gravity. However, the relatively large uncertainty still permits models within f(R) theory, which is an extension of general relativity. A fivefold decrease in uncertainty is needed to rule out these models.

We perform a Markov Chain Monte Carlo analysis of the self-accelerating and normal branch of Dvali-Gabadadze-Porrati braneworld gravity. By adopting a parameterized post-Friedmann description of gravity, we utilize all of the cosmic microwave background data, including the largest scales, and its correlation with galaxies in addition to the geometrical constraints from supernovae distances and the Hubble constant. We find that on both branches brane tension or a cosmological constant is required at high significance with no evidence for the unique Dvali-Gabadadze-Porrati modifications. The cross-over scale must therefore be substantially greater than the Hubble scale H_0 r_c > 3 and 3.5 at the 95% CL with and without uncertainties from spatial curvature. With spatial curvature, the limit from the normal branch is substantially assisted by the galaxy cross-correlation which highlights its importance in constraining infrared modifications to gravity.