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.

Bimetric theory describes gravitational interactions in the presence of an extra spin-2 field. Previous work has suggested that its cosmological solutions are generically plagued by instabilities. We show that by taking the Planck mass for the second metric, $M_f$, to be small, these instabilities can be pushed back to unobservably early times. In this limit, the theory approaches general relativity with an effective cosmological constant which is, remarkably, determined by the spin-2 interaction scale. This provides a late-time expansion history which is extremely close to $\Lambda$CDM, but with a technically-natural value for the cosmological constant. We find $M_f$ should be no larger than the electroweak scale in order for cosmological perturbations to be stable by big-bang nucleosynthesis. We further show that in this limit the helicity-0 mode is no longer strongly-coupled at low energy scales.

Bimetric gravity theories allow for many different types of cosmological solutions, but not all of them are theoretically allowed. In this work we discuss the conditions to satisfy the Higuchi bound and to avoid gradient instabilities in the scalar sector at the linear level. We find that in expanding universes the ratio of the scale factors of the reference and observable metric has to increase at all times. This automatically implies a ghost-free helicity-2 and helicity-0 sector and enforces a phantom dark energy. Furthermore, the condition for the absence of gradient instabilities in the scalar sector will be analyzed. Finally, we discuss whether cosmological solutions can exist, including exotic evolutions like bouncing cosmologies, in which both the Higuchi ghost and scalar instabilities are absent at all times.

There is a no-go theorem forbidding flat and closed FLRW solutions in massive gravity on a flat reference metric, while open solutions are unstable. Recently it was shown that this no-go theorem can be overcome if at least some matter couples to a hybrid metric composed of both the dynamical and the fixed reference metric. We show that this is not compatible with the standard description of cosmological sources in terms of effective perfect fluids, and the predictions of the theory become sensitive either to the detailed field-theoretical modelling of the matter content or to the presence of additional dark degrees of freedom. This is a serious practical complication. Furthermore, we demonstrate that viable cosmological background evolution with a perfect fluid appears to require the presence of fields with highly contrived properties. This could be improved if the equivalence principle is broken by coupling only some of the fields to the composite metric, but viable self-accelerating solutions due only to the massive graviton are difficult to obtain.

Nonlinear, ghost-free massive gravity has two tensor fields; when both are dynamical, the mass of the graviton can lead to cosmic acceleration that agrees with background data, even in the absence of a cosmological constant. Here the question of the stability of linear perturbations in this bimetric theory is examined. Instabilities are presented for several classes of models, and simple criteria for the cosmological stability of massive bigravity are derived. In this way, we identify a particular self-accelerating bigravity model, infinite-branch bigravity (IBB), which exhibits both viable background evolution and stable linear perturbations. We discuss the modified gravity parameters for IBB, which do not reduce to the standard $\Lambda$CDM result at early times, and compute the combined likelihood from measured growth data and type Ia supernovae. IBB predicts a present matter density $\Omega_{m0}=0.18$ and an equation of state $w(z)=-0.79+0.21z/(1+z)$. The growth rate of structure is well-approximated at late times by $f(z)\approx\Omega_{m}^{0.47}[1+0.21z/(1+z)]$. The implications of the linear instability for other bigravity models are discussed: the instability does not necessarily rule these models out, but rather presents interesting questions about how to extract observables from them when linear perturbation theory does not hold.

We discuss in detail a particularly simple example of a bimetric massive gravity model which seems to offer an alternative to the standard cosmological model at background level. For small redshifts, its equation of state is $w(z)\approx-1.22_{-0.02}^{+0.02}-0.64_{-0.04}^{+0.05}z/(1+z)$. Just like $\Lambda$CDM, it depends on a single parameter, has an analytical background expansion law and fits the expansion cosmological data well. However, confirming previous results, we find that the model is unstable at early times at small scales and speculate over possible ways to cure the instability. In the regime in which the model is stable, we find that it fits the linear perturbation observations well and has a growth index of approximately $\gamma=0.47$.

We find the general conditions for viable cosmological solution at the background level in bigravity models. Furthermore, we constrain the parameters by comparing to the Union 2.1 supernovae catalog and identify, in some cases analytically, the best fit parameter or the degeneracy curve among pairs of parameters. We point out that a bimetric model with a single free parameter predicts a simple relation between the equation of state and the density parameter, fits well the supernovae data and is a valid and testable alternative to $\Lambda$CDM. Additionally, we identify the conditions for a phantom behavior and show that viable bimetric cosmologies cannot cross the phantom divide.