The late time acceleration of the Universe can be characterized in terms of an extra, time dependent, component of the universe -- dark energy. The simplest proposal for dark energy is a scalar-tensor theory -- quintessence -- which consists of a scalar field, $\phi$, whose dynamics is solely dictated by its potential, $V(\phi)$. Such a theory can be uniquely characterized by the equation of state of the scalar field energy momentum-tensor. We find the time dependence of the equation of state for a broad family of potentials and, using this information, we propose an analytic prior distribution for the most commonly used parametrization. We show that this analytic prior can be used to accurately predict the distribution of observables for the next generation of cosmological surveys. Including the theoretical priors in the comparison with observations considerably improves the constraints on the equation of state.

Cosmological datasets have great potential to elucidate the nature of dark energy and test gravity on the largest scales available to observation. Theoretical predictions can be computed with hi_class (www.hiclass-code.net), an accurate, fast and flexible code for linear cosmology, incorporating a wide range of dark energy theories and modifications to general relativity. We introduce three new functionalities into hi_class: (1) Support for models based on covariant Lagrangians, including a constraint-preserving integration scheme for the background evolution and a series of worked-out examples: Galileon, nKGB, quintessence (monomial, tracker) and Brans-Dicke. (2) Consistent initial conditions for the scalar-field perturbations in the deep radiation era, identifying the conditions under which modified-gravity isocurvature perturbations may grow faster than adiabatic modes leading to a loss of predictivity. (3) An automated quasistatic approximation scheme allowing order-of-magnitude improvement in computing performance without sacrificing accuracy for wide classes of models. These enhancements bring the treatment of dark energy and modified gravity models to the level of detail comparable to software tools restricted to standard $\Lambda$CDM cosmologies. The hi_class code is publicly available (https://github.com/miguelzuma/hi_class_public), ready to explore current data and prepare for next-generation experiments.

Since the very beginning of astronomy the location of objects on the sky has been a fundamental observational quantity that has been taken for granted. While precise two dimensional positional information is easy to obtain for observations in the electromagnetic spectrum, the positional accuracy of current and near future gravitational wave detectors is limited to between tens and hundreds of square degrees, which makes it extremely challenging to identify the host galaxies of gravitational wave events or to confidently detect any electromagnetic counterparts. Gravitational wave observations provide information on source properties and distances that is complementary to the information in any associated electromagnetic emission and that is very hard to obtain in any other way. Observing systems with multiple messengers thus has scientific potential much greater than the sum of its parts. A gravitational wave detector with higher angular resolution would significantly increase the prospects for finding the hosts of gravitational wave sources and triggering a multi-messenger follow-up campaign. An observatory with arcminute precision or better could be realised within the Voyage 2050 programme by creating a large baseline interferometer array in space and would have transformative scientific potential. Precise positional information of standard sirens would enable precision measurements of cosmological parameters and offer new insights on structure formation; a high angular resolution gravitational wave observatory would allow the detection of a stochastic background and resolution of the anisotropies within it; it would also allow the study of accretion processes around black holes; and it would have tremendous potential for tests of modified gravity and the discovery of physics beyond the Standard Model.

The $\alpha$-attractor inflationary models are nowadays favored by CMB Planck observations. Their similarity with canonical quintessence models motivates the exploration of a common framework that explains both inflation and dark energy. We study the expected constraints that next-generation cosmological experiments will be able to impose for the dark energy $\alpha$-attractor model. We systematically account for the constraining power of SNIa from WFIRST, BAO from DESI and WFIRST, galaxy clustering and shear from LSST and Stage-4 CMB experiments. We assume a tensor-to-scalar ratio, $10^{-3} < r < 10^{-2}$, which permits to explore the wide regime sufficiently close, but distinct, to a cosmological constant, without need of fine tunning the initial value of the field. We find that the combination S4CMB + LSST + SNIa will achieve the best results, improving the FoM by almost an order of magnitude; respect to the S4CMB + BAO + SNIa case. We find this is also true for the FoM of the $w_0 - w_a$ parameters. Therefore, future surveys will be uniquely able to probe models connecting early and late cosmic acceleration.

Gravitational waves (GWs) provide a new tool to probe the nature of dark energy (DE) and the fundamental properties of gravity. We review the different ways in which GWs can be used to test gravity and models for late-time cosmic acceleration. Lagrangian-based gravitational theories beyond general relativity (GR) are classified into those breaking fundamental assumptions, containing additional fields and massive graviton(s). In addition to Lagrangian based theories we present the effective theory of DE and the $\mu$-$\Sigma$ parametrization as general descriptions of cosmological gravity. Multi-messenger GW detections can be used to measure the cosmological expansion (standard sirens), providing an independent test of the DE equation of state and measuring the Hubble parameter. Several key tests of gravity involve the cosmological propagation of GWs, including anomalous GW speed, massive graviton excitations, Lorentz violating dispersion relation, modified GW luminosity distance and additional polarizations, which may also induce GW oscillations. We summarize present constraints and their impact on DE models, including those arising from the binary neutron star merger GW170817. Upgrades of LIGO-Virgo detectors to design sensitivity and the next generation facilities such as LISA or Einstein Telescope will significantly improve these constraints in the next two decades.

Cosmic shear is one of the primary probes to test gravity with current and future surveys. There are two main techniques to analyse a cosmic shear survey; a tomographic method, where correlations between the lensing signal in different redshift bins are used to recover redshift information, and a 3D approach, where the full redshift information is carried through the entire analysis. Here we compare the two methods, by forecasting cosmological constraints for future surveys like Euclid. We extend the 3D formalism for the first time to theories beyond the standard model, belonging to the Horndeski class. This includes the majority of universally coupled extensions to $\Lambda$CDM with one scalar degree of freedom in addition to the metric, still in agreement with current observations. Given a fixed background, the evolution of linear perturbations in Horndeski gravity is described by a set of four functions of time only. We model their time evolution assuming proportionality to the dark energy density fraction and place Fisher matrix constraints on the proportionality coefficients. We find that a 3D analysis can constrain Horndeski theories better than a tomographic one, in particular with a decrease in the errors of the order of 20$\%$. This paper shows for the first time a quantitative comparison on an equal footing between Fisher matrix forecasts for both a fully 3D and a tomographic analysis of cosmic shear surveys. The increased sensitivity of the 3D formalism comes from its ability to retain information on the source redshifts along the entire analysis.

No Slip Gravity is a simple modified gravity theory with only one free function and the interesting characteristic that -- unlike many modified gravity theories -- it suppresses growth. This allows it to fit current redshift space distortion and $\sigma_8$ mass fluctuation amplitude data better than $\Lambda$CDM in general relativity, while retaining a $\Lambda$CDM background expansion and hence distances. Since it has no gravitational slip it alters equally CMB lensing and matter density growth, and in addition affects the CMB gravitational wave B-mode polarization power spectrum. We investigate and compute the effects of No Slip Gravity for CMB lensing and B-modes, and present a simple analytic approximation. Using a Monte Carlo analysis, we place constraints on the theory from current CMB data.

The possibility of linking inflation and late cosmic accelerated expansion using the $\alpha$-attractor models has received increasing attention due to their physical motivation. In the early universe, $\alpha$-attractors provide an inflationary mechanism compatible with Planck satellite CMB observations and predictive for future gravitational wave CMB modes. Additionally $\alpha$-attractors can be written as quintessence models with a potential that connects a power law regime with a plateau or uplifted exponential, allowing a late cosmic accelerated expansion that can mimic behavior near a cosmological constant. In this paper we study a generalized dark energy $\alpha$-attractor model. We thoroughly investigate its phenomenology, including the role of all model parameters and the possibility of large-scale tachyonic instability clustering. We verify the relation that $1+w\sim 1/\alpha$ (while the gravitational wave power $r\sim\alpha$) so these models predict that a signature should appear in either the primordial B-modes or in late time deviation from a cosmological constant. We constrain the model parameters with current datasets, including the cosmic microwave background (Planck 2015 angular power spectrum, polarization and lensing), baryon acoustic oscillations (BOSS DR12) and supernovae (Pantheon compressed). Our results show that expansion histories close to a cosmological constant exist in large regions of the parameter space, not requiring a fine-tuning of the parameters or initial conditions.

We compare Einstein-Boltzmann solvers that include modifications to General Relativity and find that, for a wide range of models and parameters, they agree to a high level of precision. We look at three general purpose codes that primarily model general scalar-tensor theories, three codes that model Jordan-Brans-Dicke (JBD) gravity, a code that models f(R) gravity, a code that models covariant Galileons, a code that models Ho\v{r}ava-Lifschitz gravity and two codes that model non-local models of gravity. Comparing predictions of the angular power spectrum of the cosmic microwave background and the power spectrum of dark matter for a suite of different models, we find agreement at the sub-percent level. This means that this suite of Einstein-Boltzmann solvers is now sufficiently accurate for precision constraints on cosmological and gravitational parameters.

Multi-messenger gravitational wave (GW) astronomy has commenced with the detection of the binary neutron star merger GW170817 and its associated electromagnetic counterparts. The almost coincident observation of both signals places an exquisite bound on the GW speed $|c_g/c-1|\leq5\cdot10^{-16}$. We use this result to probe the nature of dark energy (DE), showing that a large class of scalar-tensor theories and DE models are highly disfavored. As an example we consider the covariant Galileon, a cosmologically viable, well motivated gravity theory which predicts a variable GW speed at low redshift. Our results eliminate any late-universe application of these models, as well as their Horndeski and most of their beyond Horndeski generalizations. Three alternatives (and their combinations) emerge as the only possible scalar-tensor DE models: 1) restricting Horndeski's action to its simplest terms, 2) applying a conformal transformation which preserves the causal structure and 3) compensating the different terms that modify the GW speed (to be robust, the compensation has to be independent on the background on which GWs propagate). Our conclusions extend to any other gravity theory predicting varying $c_g$ such as Einstein-Aether, Ho\v{r}ava gravity, Generalized Proca, TeVeS and other MOND-like gravities.

Cosmological models with Galileon gravity are an alternative to the standard $\Lambda {\rm CDM}$ paradigm with testable predictions at the level of its self-accelerating solutions for the expansion history, as well as large-scale structure formation. Here, we place constraints on the full parameter space of these models using data from the cosmic microwave background (CMB) (including lensing), baryonic acoustic oscillations (BAO) and the Integrated Sachs-Wolfe (ISW) effect. We pay special attention to the ISW effect for which we use the cross-spectra, $C_\ell^{\rm T g}$, of CMB temperature maps and foreground galaxies from the WISE survey. The sign of $C_\ell^{\rm T g}$ is set by the time evolution of the lensing potential in the redshift range of the galaxy sample: it is positive if the potential decays (like in $\Lambda {\rm CDM}$), negative if it deepens. We constrain three subsets of Galileon gravity separately known as the Cubic, Quartic and Quintic Galileons. The cubic Galileon model predicts a negative $C_\ell^{\rm T g}$ and exhibits a $7.8\sigma$ tension with the data, which effectively rules it out. For the quartic and quintic models the ISW data also rule out a significant region of the parameter space but permit regions where the goodness-of-fit is comparable to $\Lambda {\rm CDM}$. The data prefers a non zero sum of the neutrino masses ($\sum m_\nu\approx 0.5$eV) with $ \sim \! 5\sigma$ significance in these models. The best-fitting models have values of $H_0$ consistent with local determinations, thereby avoiding the tension that exists in $\Lambda {\rm CDM}$. We also identify and discuss a $\sim \! 2\sigma$ tension that Galileon gravity exhibits with recent BAO measurements. Our analysis shows overall that Galileon cosmologies cannot be ruled out by current data but future lensing, BAO and ISW data hold strong potential to do so.

The direct detection of gravitational waves (GWs) is an invaluable new tool to probe gravity and the nature of cosmic acceleration. A large class of scalar-tensor theories predict that GWs propagate with velocity different than the speed of light, a difference that can be $\mathcal{O}(1)$ for many models of dark energy. We determine the conditions behind the anomalous GW speed, namely that the scalar field spontaneously breaks Lorentz invariance and couples to the metric perturbations via the Weyl tensor. If these conditions are realized in nature, the delay between GW and electromagnetic (EM) signals from distant events will run beyond human timescales, making it impossible to measure the speed of GWs using neutron star mergers or other violent events. We present a robust strategy to exclude or confirm an anomalous speed of GWs using eclipsing binary systems, whose EM phase can be exquisitely determined. he white dwarf binary J0651+2844 is a known example of such system that can be used to probe deviations in the GW speed as small as $c_g/c-1\gtrsim 2\cdot 10^{-12}$ when LISA comes online. This test will either eliminate many contender models for cosmic acceleration or wreck a fundamental pillar of general relativity.

We study the role of field redefinitions in general scalar-tensor theories. In particular, we first focus on the class of field redefinitions linear in the spin-2 field and involving derivatives of the spin-0 mode, generically known as disformal transformations. We start by defining the action of a disformal transformation in the tangent space. Then, we take advantage of the great economy of means of the language of differential forms to compute the full transformation of Horndeski's theory under general disformal transformations. We obtain that Horndeski's action maps onto itself modulo a reduced set of non-Horndeski Lagrangians. These new Lagrangians are found to be invariant under disformal transformation that depend only in the first derivatives of the scalar. Moreover, these combinations of Lagrangians precisely appear when expressing in our basis the constraints of the recently proposed Extended Scalar-Tensor (EST) theories. These results allow us to classify the different orbits of scalar-tensor theories invariant under particular disformal transformations, namely the special disformal, kinetic disformal and disformal Horndeski orbits. In addition, we consider generalizations of this framework. We find that there are possible well-defined extended disformal transformations that have not been considered in the literature. However, they generically cannot link Horndeski theory with EST theories. Finally, we study further generalizations in which extra fields with different spin are included. These field redefinitions can be used to connect different gravity theories such as multi-scalar-tensor theories, generalized Proca theories and bi-gravity. We discuss how the formalism of differential forms could be useful for future developments in these lines.

Many theories of modified gravity with higher order derivatives are usually ignored because of serious problems that appear due to an additional ghost degree of freedom. Most dangerously, it causes an immediate decay of the vacuum. However, breaking Lorentz invariance can cure such abominable behavior. By analyzing a model that describes a massive graviton together with a remaining Boulware-Deser ghost mode we show that even ghostly theories of modified gravity can yield models that are viable at both classical and quantum levels and, therefore, they should not generally be ruled out. Furthermore, we identify the most dangerous quantum scattering process that has the main impact on the decay time and find differences to simple theories that only describe an ordinary scalar field and a ghost. Additionally, constraints on the parameters of the theory including some upper bounds on the Lorentz-breaking cutoff scale are presented. In particular, for a simple theory of massive gravity we find that a breaking of Lorentz invariance is allowed to happen even at scales above the Planck mass. Finally, we discuss the relevance to other theories of modified gravity.

We use a description based on differential forms to systematically explore the space of scalar-tensor theories of gravity. Within this formalism, we propose a basis for the scalar sector at the lowest order in derivatives of the field and in any number of dimensions. This minimal basis is used to construct a finite and closed set of Lagrangians describing general scalar-tensor theories invariant under Local Lorentz Transformations in a pseudo-Riemannian manifold, which contains ten physically distinct elements in four spacetime dimensions. Subsequently, we compute their corresponding equations of motion and find which combinations are at most second order in derivatives in four as well as arbitrary number of dimensions. By studying the possible exact forms (total derivatives) and algebraic relations between the basis components, we discover that there are only four Lagrangian combinations producing second order equations, which can be associated with Horndeski's theory. In this process, we identify a new second order Lagrangian, named kinetic Gauss-Bonnet, that was not previously considered in the literature. However, we show that its dynamics is already contained in Horndeski's theory. Finally, we provide a full classification of the relations between different second order theories. This allows us to clarify, for instance, the connection between different covariantizations of Galileons theory. In conclusion, our formulation affords great computational simplicity with a systematic structure. As a first step we focus on theories with second order equations of motion. However, this new formalism aims to facilitate advances towards unveiling the most general scalar-tensor theories.

Kinetic mixing between the metric and scalar degrees of freedom is an essential ingredient in contemporary scalar-tensor theories. This often makes hard to understand their physical content, especially when derivative mixing is present, as it is the case for Horndeski action. In this work we develop a method that allows to write a Ricci curvature-free scalar field equation and discuss some of the advantages of such rephrasing in the study of stability issues in the presence of matter, the existence of an Einstein frame and the generalization of the disformal screening mechanism. For quartic Horndeski theories, such procedure leaves, in general, a residual coupling to curvature, given by the Weyl tensor. This gives rise to a binary classification of scalar-tensor theories into stirred theories, for which the curvature can be substituted for, and shaken theories for which a residual coupling to curvature remains. Quite remarkably, we have found that generalized DBI Galileons belong to the first class. Finally, we discuss kinetic mixing in quintic theories for which non-linear mixing terms appears and in the recently proposed theories beyond Horndeski which display a novel form of kinetic mixing, in which the field equation is sourced by derivatives of the energy-momentum tensor.

We study the structure of scalar-tensor theories of gravity based on derivative couplings between the scalar and the matter degrees of freedom introduced through an effective metric. Such interactions are classified by their tensor structure into conformal (scalar), disformal (vector) and extended disformal (traceless tensor), as well as by the derivative order of the scalar field. Relations limited to first derivatives of the field ensure second order equations of motion in the Einstein frame and hence the absence of Ostrogradski ghost degrees of freedom. The existence of a mapping to the Jordan frame is not trivial in the general case, and can be addressed using the Jacobian of the frame transformation through its eigenvalues and eigentensors. These objects also appear in the study of different aspects of such theories, including the metric and field redefinition transformation of the path integral in the quantum mechanical description. Although sane in the Einstein frame, generic disformally coupled theories are described by higher order equations of motion in the Jordan frame. This apparent contradiction is solved by the use of a hidden constraint: the contraction of the metric equations with a Jacobian eigentensor provides a constraint relation for the higher field derivatives, which allows one to express the dynamical equations in a second order form. This signals a loophole in Horndeski's theorem and allows one to enlarge the set of scalar-tensor theories which are Ostrogradski-stable. The transformed Gauss-Bonnet terms are also discussed for the simplest conformal and disformal relations.