Discrepant measurements of the Universe's expansion rate ($H_0$) may signal physics beyond the standard cosmological model. Here I describe two early modified gravity mechanisms that reconcile the value of $H_0$ by increasing the expansion rate in the era of matter-radiation equality. These mechanisms, based on viable Horndeski theories, require significantly less fine-tuned initial conditions than early dark energy with oscillating scalar fields. In Imperfect Dark Energy at Equality (IDEE), the initial energy density dilutes slower than radiation but faster than matter, naturally peaking around the era of equality. The minimal IDEE model, a cubic Galileon, is too constrained by the cosmic microwave background (Planck) and baryon acoustic oscillations (BAO) to relieve the $H_0$ tension. In Enhanced Early Gravity (EEG), the scalar field value modulates the cosmological strength of gravity. The minimal EEG model, an exponentially coupled cubic Galileon, gives a Planck+BAO value $H_0=68.7 \pm 1.5$ (68\% c.l.), reducing the tension with SH0ES from $4.4\sigma$ to $2.6\sigma$. Additionally, Galileon contributions to cosmic acceleration may reconcile $H_0$ via Late-Universe Phantom Expansion (LUPE). Combining LUPE, EEG and $\Lambda$ reduces the tension between Planck, BAO and SH0ES to $2.5\sigma$. I will also describe additional tests of coupled Galileons based on local gravity tests, primordial element abundances and gravitational waves. While further model building is required to fully resolve the $H_0$ problem and satisfy all available observations, these examples show the wealth of possibilities to solve cosmological tensions beyond Einstein's General Relativity.

Future gravitational-wave observations will enable unprecedented and unique science in extreme gravity and fundamental physics answering questions about the nature of dynamical spacetimes, the nature of dark matter and the nature of compact objects.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

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]

The nature of dark matter (DM) remains unknown despite very precise knowledge of its abundance in the universe. An alternative to new elementary particles postulates DM as made of macroscopic compact halo objects (MACHO) such as black holes formed in the very early universe. Stellar-mass primordial black holes (PBHs) are subject to less robust constraints than other mass ranges and might be connected to gravitational-wave signals detected by the Laser Interferometer Gravitational-Wave Observatory (LIGO). New methods are therefore necessary to constrain the viability of compact objects as a DM candidate. Here we report bounds on the abundance of compact objects from gravitational lensing of type Ia supernovae (SNe). Current SNe datasets constrain compact objects to represent less than 35.2% (Joint Lightcurve Analisis) and 37.2% (Union 2.1) of the total matter content in the universe, at 95% confidence-level. The results are valid for masses larger than $\sim 0.01M_\odot$ (solar-masses), limited by the size SNe relative to the lens Einstein radius. We demonstrate the mass range of the constraints by computing magnification probabilities for realistic SNe sizes and different values of the PBH mass. Our bounds are sensitive to the total abundance of compact objects with $M \lesssim 0.01M_\odot$ and complementary to other observational tests. These results are robust against cosmological parameters, outlier rejection, correlated noise and selection bias. PBHs and other MACHOs are therefore ruled out as the dominant form of DM for objects associated to LIGO gravitational wave detections. These bounds constrain early-universe models that predict stellar-mass PBH production and strengthen the case for lighter forms of DM, including new elementary particles.

We present the public version of hi_class (www.hiclass-code.net), an extension of the Boltzmann code CLASS to a broad ensemble of modifications to general relativity. In particular, hi_class can calculate predictions for models based on Horndeski's theory, which is the most general scalar-tensor theory described by second-order equations of motion and encompasses any perfect-fluid dark energy, quintessence, Brans-Dicke, $f(R)$ and covariant Galileon models. hi_class has been thoroughly tested and can be readily used to understand the impact of alternative theories of gravity on linear structure formation as well as for cosmological parameter extraction.

The next generation of surveys will greatly improve our knowledge of cosmological gravity. In this paper we focus on how Stage IV photometric redshift surveys, including weak lensing and multiple tracers of the matter distribution and radio experiments combined with measurements of the cosmic microwave background will lead to precision constraints on deviations from General Relativity. We use a broad subclass of Horndeski scalar-tensor theories to forecast the accuracy with which we will be able to determine these deviations and their degeneracies with other cosmological parameters. Our analysis includes relativistic effects, does not rely on the quasi-static evolution and makes conservative assumptions about the effect of screening on small scales. We define a figure of merit for cosmological tests of gravity and show how the combination of different types of surveys, probing different length scales and redshifts, can be used to pin down constraints on the gravitational physics to better than a few percent, roughly an order of magnitude better than present probes. Future cosmological experiments will be able to constrain the Brans-Dicke parameter at a level comparable to Solar System and astrophysical tests

We address the impact of consistent modifications of gravity on the largest observable scales, focusing on relativistic effects in galaxy number counts and the cross-correlation between the matter large scale structure (LSS) distribution and the cosmic microwave background (CMB). Our analysis applies to a very broad class of general scalar-tensor theories encoded in the Horndeski Lagrangian and is fully consistent on linear scales, retaining the full dynamics of the scalar field and not assuming quasi-static evolution. As particular examples we consider self-accelerating Covariant Galileons, Brans-Dicke theory and parameterizations based on the effective field theory of dark energy, using the \hiclass\, code to address the impact of these models on relativistic corrections to LSS observables. We find that especially effects which involve integrals along the line of sight (lensing convergence, time delay and the integrated Sachs-Wolfe effect -- ISW) can be considerably modified, and even lead to $\mathcal{O}(1000\%)$ deviations from General Relativity in the case of the ISW effect for Galileon models, for which standard probes such as the growth function only vary by $\mathcal{O}(10\%)$. These effects become dominant when correlating galaxy number counts at different redshifts and can lead to $\sim 50\%$ deviations in the total signal that might be observable by future LSS surveys. Because of their integrated nature, these deep-redshift cross-correlations are sensitive to modifications of gravity even when probing eras much before dark energy domination. We further isolate the ISW effect using the cross-correlation between LSS and CMB temperature anisotropies and use current data to further constrain Horndeski models (abridged).

The scale of Baryon Acoustic Oscillations (BAO) imprinted in the matter power spectrum provides an almost-perfect standard ruler: it only suffers sub-percent deviations from fixed comoving length due to non-linear effects. We study the BAO shift in the large Horndeski class of gravitational theories and compute its magnitude in momentum space using second-order perturbation theory and a peak-background split. The standard prediction is affected by the modified linear growth, as well as by non-linear gravitational effects that alter the mode-coupling kernel. For covariant Galileon models, we find a $14-45\%$ enhancement of the BAO shift with respect to standard gravity and a distinct time evolution depending on the parameters. Despite the larger values, the shift remains well below the forecasted precision of next-generation galaxy surveys. Models that produce significant BAO shift would cause large redshift-space distortions or affect the bispectrum considerably. Our computation therefore validates the use of the BAO scale as a comoving standard ruler for tests of general dark energy models.

The theory of bigravity offers one of the simplest possibilities to describe a massive graviton while having self-accelerating cosmological solutions without a cosmological constant. However, it has been shown recently that bigravity is affected by early-time fast growing modes on the tensor sector. Here we argue that we can only trust the linear analysis up to when perturbations are in the linear regime and use a cut-off to stop the growing of the metric perturbations. This analysis, although more consistent, still leads to growing tensor modes that are unacceptably large for the theory to be compatible with measurements of the cosmic microwave background (CMB), both in temperature and polarization spectra. In order to suppress the growing modes and make the model compatible with CMB spectra, we find it necessary to either fine-tune the initial conditions, modify the theory or set the cut-off for the tensor perturbations of the second metric much lower than unity. Initial conditions such that the growing mode is sufficiently suppresed can be achieved in scenarios in which inflation ends at the GeV scale.

It is shown that a disformally coupled theory in which the gravitational sector has the Einstein-Hilbert form is equivalent to a quartic DBI Galileon Lagrangian, possessing non-linear higher derivative interactions, and hence allowing for the Vainshtein effect. This Einstein Frame description considerably simplifies the dynamical equations and highlights the role of the different terms. The study of highly dense, non-relativistic environments within this description unravels the existence of a disformal screening mechanism, while the study of static vacuum configurations reveals the existence of a Vainshtein radius, at which the asymptotic solution breaks down. Disformal couplings to matter also allow the construction of Dark Energy models, which behave differently than conformally coupled ones and introduce new effects on the growth of Large Scale Structure over cosmological scales, on which the scalar force is not screened. We consider a simple Disformally Coupled Dark Matter model in detail, in which standard model particles follow geodesics of the gravitational metric and only Dark Matter is affected by the disformal scalar field. This particular model is not compatible with observations in the linearly perturbed regime. Nonetheless, disformally coupled theories offer enough freedom to construct realistic cosmological scenarios, which can be distinguished from the standard model through characteristic signatures.

It is shown that extensions to General Relativity, which introduce a strongly coupled scalar field, can be viable if the interaction has a non-conformal form. Such disformal coupling depends upon the gradients of the scalar field. Thus, if the field is locally static and smooth, the coupling becomes invisible in the solar system: this is the disformal screening mechanism. A cosmological model is considered where the disformal coupling triggers the onset of accelerated expansion after a scaling matter era, giving a good fit to a wide range of observational data. Moreover, the interaction leaves signatures in the formation of large-scale structure that can be used to probe such couplings.

New constraints on inhomogeneous Lema\^itre-Tolman-Bondi (LTB) models alternative to Dark Energy are presented, focusing on adiabatic profiles with space-independent Big Bang and baryon fraction. The Baryon Acoustic Oscillation (BAO) scale at early times is computed in terms of the asymptotic value and then projected to different redshifts by following the geodesics of the background metric. Additionally, a model-independent method to constraint the local expansion rate using a prior on supernovae luminosity is presented. Cosmologies described by an adiabatic GBH matter profile with \Omega_{out}=1 and \Omega_{out}\leq 1 are investigated using a Markov Chain Monte Carlo analysis including the latest BAO data from the WiggleZ collaboration and the local expansion rate from the Hubble Space Telescope, together with Union-II type Ia supernovae data and the position and height of the Cosmic Microwave Background acoustic peaks. The addition of BAO data at higher redshifts increases considerably their constraining power and represents a new drawback for this type of models, yielding a value of the local density parameter \Omega_{in}>0.2 which is 3\sigma apart from the value \Omega_{in}<0.15 found using supernovae. The situation does not improve if the asymptotic flatness assumption is dropped, and a Bayesian analysis shows that constrained GBH models are ruled out at high confidence. We emphasize that these are purely geometric probes, that only recently have become sufficiently constraining to independently rule out the whole class of adiabatic LTB models.

It has been recently proposed that gravity might be an entropic force. Although a well defined fundamental description for such a mechanism is still lacking, it is still possible to address the viability of phenomenological models of entropic-inspired modified gravities. I will summarize some recent work directed to using cosmology as a tool to constraint scenarios in which the modifications are aimed to explain the physics behind dark energy and inflation. A phenomenological modification is able to explain cosmic acceleration at the background level and fit observations, but simple inflation models with higher curvature corrections are in conflict with late time matter domination.

It has been recently proposed that the interpretation of gravity as an emergent, entropic force might have nontrivial implications to cosmology. Here two approaches are investigated: in one, the Friedman equation receives entropic contributions from the usually neglected surface terms, and in another, the extra terms are derived from quantum corrections to the entropy formula. UV terms may drive inflation, avoiding a recently derived no-go theorem, though in some cases leading to a graceful exit problem. IR terms can generate dark energy, alleviating the cosmological constant problem. The quantum corrections are bounded by their implications to the BBN, and the surface terms are constrained in addition by their effect upon the behavior of matter. Likelihood analyses are performed to constrain the modifications by the SNeIa, BAO and CMB data. It is found that a monomial correction to the area-entropy formula results in late acceleration in very good agreement with observations, which then turn out to be compatible with positive curvature. The evolution of perturbations is deduced by assuming the Jebsen-Birkhoff theorem. Distinct signatures can then be identified in the large scale structure formation. Furthermore, it is shown that the visible universe satisfies the Bekenstein bound.

Disformal transformations have proven to be very useful to devise models of the dark sector. In the present paper we apply such transformation to a single scalar field theory as a way to drive the field into a slow roll phase. The canonical scalar field Lagrangian, when coupled to a disformal metric, turns out to have relations to bimetric dark matter theories and to describe many specific dark energy models at various limits, thus providing a surprisingly simple parametrisation of a wide variety of models including tachyon, Chaplygin gas, K-essence and dilatonic ghost condensate. We investigate the evolution of the background and linear perturbations in disformal quintessence in order to perform a full comparison of the predictions with the cosmological data. The dynamics of the expansion, in particular the mechanism of the transition to accelerating phase, is described in detail. We then study the effects of disformal quintessence on cosmic microwave background (CMB) anisotropies and large scale structures (LSS). A likelihood analysis using the latest data on wide-ranging SNIa, CMB and LSS observations is performed allowing variations in six cosmological parameters and the two parameters specifying the model. We find that while a large region of parameter space remains compatible with observations, models featuring either too much early dark energy or too slow transition to acceleration are ruled out.