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.

In this paper, we study the Starobinsky model of inflation in the context of gravity's rainbow theory. We propose that gravity rainbow functions can be written in the power-law form of the Hubble parameter. We present a detailed derivation of the spectral index of curvature perturbation and the tensor-to-scalar ratio and compare the predictions of our models with Planck 2015 data. We discover that in order to be consistent with Planck data up to $2\sigma$ C.L., the viable values of $N_{k}$ $e$-folds would satisfy $42\lesssim N_{k}\lesssim 87$ and the rainbow parameter $\lambda$ is nicely constrained to be $\lambda \lesssim6.0$.