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.

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 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.