We study perturbation theory for large-scale structure in the most general scalar-tensor theories propagating a single scalar degree of freedom, which include Horndeski theories and beyond. We model the parameter space using the effective field theory of dark energy. For Horndeski theories, the gravitational field and fluid equations are invariant under a combination of time-dependent transformations of the coordinates and fields. This symmetry allows one to construct a physical adiabatic mode which fixes the perturbation-theory kernels in the squeezed limit and ensures that the well-known consistency relations for large-scale structure, originally derived in general relativity, hold in modified gravity as well. For theories beyond Horndeski, instead, one generally cannot construct such an adiabatic mode. Because of this, the perturbation-theory kernels are modified in the squeezed limit and the consistency relations for large-scale structure do not hold. We show, however, that the modification of the squeezed limit depends only on the linear theory. We investigate the observational consequences of this violation by computing the matter bispectrum. In the squeezed limit, the largest effect is expected when considering the cross-correlation between different tracers. Moreover, the individual contributions to the 1-loop matter power spectrum do not cancel in the infrared limit of the momentum integral, modifying the power spectrum on non-linear scales.

We study the screening mechanism in the most general scalar-tensor theories that leave gravitational waves unaffected and are thus compatible with recent LIGO/Virgo observations. Using the effective field theory of dark energy approach, we consider the general action for perturbations beyond linear order, focussing on the quasi-static limit. When restricting to the subclass of theories that satisfy the gravitational wave constraints, the fully nonlinear effective Lagrangian contains only three independent parameters. One of these, $\beta_1$, is uniquely present in degenerate higher-order theories. We compute the two gravitational potentials for a spherically symmetric matter source and we find that for $\beta_1 \ge 0$ they decrease as the inverse of the distance, as in standard gravity, while the case $\beta_1 < 0$ is ruled out. For $\beta_1 > 0$, the two potentials differ and their gravitational constants are not the same on the inside and outside of the body. Generically, the bound on anomalous light bending in the Solar System constrains $\beta_1 \lesssim 10^{-5}$. Standard gravity can be recovered outside the body by tuning the parameters of the model, in which case $\beta_1 \lesssim 10^{-2}$ from the Hulse-Taylor pulsar.

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 develop a new technique for finding black hole solutions in modified gravity that have "stealth" hair, i.e., hair whose only gravitational effect is to tune the cosmological constant. We consider scalar-tensor theories in which gravitational waves propagate at the speed of light, and show that Einstein metrics can be painted with stealth hair provided there exists a family of geodesics always normal to spacelike surfaces. We also present a novel scalar-dressed rotating black hole that has finite scalar field at both the black hole and cosmological event horizons.

In the context of Degenerate Higher-Order Scalar-Tensor (DHOST) theories, we study cosmological solutions and their stability properties. In particular, we explicitly illustrate the crucial role of degeneracy by showing how the higher order homogeneous equations in the physical frame (where matter is minimally coupled) can be recast in a system of equations that do not involve higher order derivatives. We study the fixed points of the dynamics, finding the conditions for having a de Sitter attractor at late times. Then we consider the coupling to matter field (described for convenience by a k-essence Lagrangian) and find the conditions to avoid gradient and ghost instabilities at linear order in cosmological perturbations, extending previous work. Finally, we apply these results to a simple subclass of DHOST theories, showing that de Sitter attractor conditions, no ghost and no gradient instabilities conditions (both in the self-accelerating era and in the matter dominated era) can be compatible.

The recent simultaneous detection of gravitational waves and a gamma ray burst from a neutron star merger significantly shrank the space of viable scalar-tensor theories by demanding that the speed of gravity is equal to that of light. The survived theories belong to the class of degenerate higher order scalar-tensor theories. We study whether these theories are suitable as dark energy candidates. We find scaling solutions in the matter dominated universe that lead to de Sitter solutions at late times without the cosmological constant, realising self-acceleration. We evaluate quasi-static perturbations around self-accelerating solutions and show that the stringent constraints coming from astrophysical objects and gravitational waves can be satisfied, leaving interesting possibilities to test these theories by cosmological observations.

We consider theories describing the dynamics of a four-dimensional metric, whose Lagrangian is diffeomorphism invariant and depends at most on second derivatives of the metric. Imposing degeneracy conditions we find a set of Lagrangians that, apart form the Einstein-Hilbert one, are either trivial or contain more than two degrees of freedom. Among the partially degenerate theories, we recover Chern-Simons gravity, endowed with constraints whose structure suggests the presence of instabilities. Then, we enlarge the class of parity violating theories of gravity by introducing new "chiral scalar-tensor theories". Although they all raise the same concern as Chern-Simons gravity, they can nevertheless make sense as low energy effective field theories or, by restricting them to the unitary gauge (where the scalar field is uniform), as Lorentz breaking theories with a parity violating sector.

The almost simultaneous detection of gravitational waves and a short gamma-ray burst from a neutron star merger has put a tight constraint on the difference between the speed of gravity and light. In the four-dimensional scalar-tensor theory with second order equations of motion, the Horndeski theory, this translates into a significant reduction of the viable parameter space of the theory. Recently, extensions of Horndeski theory, which are free from Ostrogradsky ghosts despite the presence of higher order derivatives in the equations of motion, have been identified and classified exploiting the degeneracy criterium. In these new theories, the fifth force mediated by the scalar field must be suppressed in order to evade the stringent Solar System constraints. We study the Vainshtein mechanism in the most general degenerate higher order scalar-tensor theory in which light and gravity propagate at the same speed. We find that the Vainshtein mechanism generally works outside a matter source but it is broken inside matter, similarly to beyond Horndeski theories. This leaves interesting possibilities to test these theories that are compatible with gravitational wave observations using astrophysical objects.

We provide a full analysis of ghost free higher derivative field theories with coupled degrees of freedom. Assuming the absence of gauge symmetries, we derive the degeneracy conditions in order to evade the Ostrogradsky ghosts, and analyze which (non)trivial classes of solutions this allows for. It is shown explicitly how Lorentz invariance avoids the propagation of "half" degrees of freedom. Moreover, for a large class of theories, we construct the field redefinitions and/or (extended) contact transformations that put the theory in a manifestly first order form. Finally, we identify which class of theories cannot be brought to first order form by such transformations.

We present all scalar-tensor Lagrangians that are cubic in second derivatives of a scalar field, and that are degenerate, hence avoiding Ostrogradsky instabilities. Thanks to the existence of constraints, they propagate no more than three degrees of freedom, despite having higher order equations of motion. We also determine the viable combinations of previously identified quadratic degenerate Lagrangians and the newly established cubic ones. Finally, we study whether the new theories are connected to known scalar-tensor theories such as Horndeski and beyond Horndeski, through conformal or disformal transformations.

We study new consistent scalar-tensor theories of gravity recently introduced by Langlois and Noui with potentially interesting cosmological applications. We derive the conditions for the existence of a primary constraint that prevents the propagation of an additional dangerous mode associated with higher order equations of motion. We then classify the most general, consistent scalar-tensor theories that are at most quadratic in the second derivatives of the scalar field. In addition, we investigate the possible connection between these theories and (beyond) Horndeski through conformal and disformal transformations. Finally, we point out that these theories can be associated with new operators in the effective field theory of dark energy, which might open up new possibilities to test dark energy models in future surveys.

Determining the most general, consistent scalar tensor theory of gravity is important for building models of inflation and dark energy. In this work we investigate the number of degrees of freedom present in the theory of beyond Horndeski. We discuss how to construct the theory from the extrinsic curvature of the constant scalar field hypersurface, and find a simple expression for the action which guarantees the existence of the primary constraint necessary to avoid the Ostrogradsky instability. Our analysis is completely gauge-invariant. However we confirm that, mixing together beyond Horndeski with a different order of Horndeski, obstructs the construction of this primary constraint. Instead, when the mixing is between actions of the same order, the theory can be mapped to Horndeski through a generalised disformal transformation. This mapping however is impossible with beyond Horndeski alone, since we find that the theory is invariant under such a transformation. The picture that emerges is that beyond Horndeski is a healthy but isolated theory: combined with Horndeski, it either becomes Horndeski, or likely propagates a ghost.

The basic building block for Lorentz invariant and ghost free massive gravity is the square root of the combination $g^{-1}\eta\,$, where $g^{-1}$ is the inverse of the physical metric and $\eta$ is a reference metric. Since the square root of a matrix is not uniquely defined, it is possible to have physically inequivalent potentials corresponding to different branches. We show that around Minkowski background the only perturbatively well defined branch is the potential proposed by de Rham, Gabadadze and Tolley. On the other hand, if Lorentz symmetry is broken spontaneously, other potentials exist with a standard perturbative expansion. We show this explicitly building new Lorentz invariant, ghost-free massive gravity potentials for theories that in the background preserve rotational invariance, but break Lorentz boosts.

We study cosmology in the bigravity formulation of the dRGT model where matter couples to both metrics. At linear order in perturbation theory two mass scales emerge: an hard one from the dRGT potential, and an environmental dependent one from the coupling of bigravity with matter. At early time, the dynamics is dictated by the second mass scale which is of order of the Hubble scale. The set of gauge invariant perturbations that couples to matter follow closely the same behaviour as in GR. The remaining perturbations show no issue in the scalar sector, while problems arise in the tensor and vector sectors. During radiation domination, a tensor mode grows power-like at super-horizon scales. More dangerously, the only propagating vector mode features an exponential instability on sub-horizon scales. We discuss the consequences of such instabilities and speculate on possible ways to deal with them.

Cosmological perturbations of FRW solutions in ghost free massive bigravity, including also a second matter sector, are studied in detail. At early time, we find that sub horizon exponential instabilities are unavoidable and they lead to a premature departure from the perturbative regime of cosmological perturbations.

We study cosmological perturbations for a ghost free massive gravity theory formulated with a dynamical extra metric that is needed to massive deform GR. In this formulation FRW background solutions fall in two branches. In the dynamics of perturbations around the first branch solutions, no extra degree of freedom with respect to GR ispresent at linearized level, likewise what is found in the Stuckelberg formulation of massive gravity where the extra metric isflat and non dynamical. In the first branch, perturbations are probably strongly coupled. On the contrary, for perturbations around the second branch solutions all expected degrees of freedom propagate. While tensor and vector perturbations of the physical metric that couples with matter follow closely the ones of GR, scalars develop an exponential Jeans-like instability on sub-horizon scales. On the other hand, around a de Sitter background there is no instability. We argue that one could get rid of the instabilities by introducing a mirror dark matter sector minimally coupled to only the second metric.

We study FRW homogeneous cosmological solutions in the recently found ghost free massive gravity. In previous works it was shown that when the additional extra metric, needed to generate the mass term, is taken as non-dynamical and flat, no homogeneous FRW cosmology exists. We show that, when the additional metric is a dynamical field, a perfectly safe FRW universe exists. FRW solutions fall in two branches. In the first branch the massive deformation is equivalent to an effectively generated cosmological constant whose scale is determined by the graviton mass. The second branch is quite rich: we have FRW cosmology in the presence of a "gravitational" fluid. The control parameter xi is the ratio of the two conformal factors. When xi is large, generically the cosmological evolution greatly differs from GR at the early time. In the small xi region, the evolution is similar to GR and the universe flows at late time toward an attractor represented by a dS phase.

In GR the static gravitational potential of a self-gravitating body goes as 1/r at large distances and any slower decrease leads to infinity energy. We show that in a class of four-dimensional massive gravity theories there exists spherically symmetric solutions with finite total energy, featuring an asymptotic behavior slower than 1/r and generically of the form $r^\gamma$. This suggests that configurations with nonstandard asymptotics may well turn out to be physical. The effect is due to an extra field coupled only gravitationally, which allows for modifications of the static potential generated by matter, while counterbalancing the apparently infinite energy budget.