We point out that dark-energy perturbations may become unstable in the presence of a gravitational wave of sufficiently large amplitude. We study this effect for the cubic Horndeski operator (braiding), proportional to $\alpha_{\rm B}$. The scalar that describes dark-energy fluctuations features ghost and/or gradient instabilities for gravitational-wave amplitudes that are produced by typical binary systems. Taking into account the populations of binary systems, we conclude that the instability is triggered in the whole Universe for $|\alpha_{\rm B} |\gtrsim 10^{-2}$, i.e. when the modification of gravity is sizeable. The instability is triggered by massive black-hole binaries down to frequencies corresponding to $10^{10}$ km: the instability is thus robust, unless new physics enters on even longer wavelengths. The fate of the instability and the subsequent time-evolution of the system depend on the UV completion, so that the theory may end up in a state very different from the original one. The same kind of instability is present in beyond-Horndeski theories for $|\alpha_{\rm H}| \gtrsim 10^{-20}$. In conclusion, the only dark-energy theories with sizeable cosmological effects that avoid these problems are $k$-essence models, with a possible conformal coupling with matter.

We study the decay of gravitational waves into dark energy fluctuations $\pi$, taking into account the large occupation numbers. We describe dark energy using the effective field theory approach, in the context of generalized scalar-tensor theories. When the $m_3^3$ (cubic Horndeski) and $\tilde m_4^2$ (beyond Horndeski) operators are present, the gravitational wave acts as a classical background for $\pi$ and modifies its dynamics. In particular, $\pi$ fluctuations are described by a Mathieu equation and feature instability bands that grow exponentially. Focusing on the regime of small gravitational-wave amplitude, corresponding to narrow resonance, we calculate analytically the produced $\pi$, its energy and the change of the gravitational-wave signal. The resonance is affected by $\pi$ self-interactions in a way that we cannot describe analytically. This effect is very relevant for the operator $m_3^3$ and it limits the instability. In the case of the $\tilde m_4^2$ operator self-interactions can be neglected, at least in some regimes. The modification of the gravitational-wave signal is observable for $3 \times 10^{-20} \lesssim \alpha_{\rm H} \lesssim 10^{-17}$ with a LIGO/Virgo-like interferometer and for $10^{-16} \lesssim \alpha_{\rm H} \lesssim 10^{-10}$ with a LISA-like one.

We study the decay of gravitational waves into dark energy fluctuations $\pi$, through the processes $\gamma \to \pi\pi$ and $\gamma \to \gamma \pi$, made possible by the spontaneous breaking of Lorentz invariance. Within the EFT of Dark Energy (or Horndeski/beyond Horndeski theories) the first process is large for the operator $\frac12 \tilde m_4^2(t) \, \delta g^{00}\, \left( {}^{(3)}\! R + \delta K_\mu^\nu \delta K^\mu_\nu -\delta K^2 \right)$, so that the recent observations force $\tilde m_4 =0$ (or equivalently $\alpha_{\rm H}=0$). This constraint, together with the requirement that gravitational waves travel at the speed of light, rules out all quartic and quintic GLPV theories. Additionally, we study how the same couplings affect the propagation of gravitons at loop order. The operator proportional to $\tilde m_4^2$ generates a calculable, non-Lorentz invariant higher-derivative correction to the graviton propagation. The modification of the dispersion relation provides a bound on $\tilde m_4^2$ comparable to the one of the decay. Conversely, operators up to cubic Horndeski do not generate sizeable higher-derivative corrections.

We consider the scale-invariant inflationary model studied in [1]. The Lagrangian includes all the scale-invariant operators that can be built with combinations of $R$, $R^2$ and one scalar field. The equations of motion show that the symmetry is spontaneously broken after an arbitrarily long inflationary period and a fundamental mass scale is generated. Upon symmetry breaking, and in the Jordan frame, both Hubble function and the scalar field undergo damped oscillations that can eventually amplify Standard Model fields and reheat the Universe. In the present work, we study in detail inflation and the reheating mechanism of this model in the Einstein frame and we compare some of the results with the latest observational data.