We construct a fully covariant theory of massive gravity which does not require the introduction of an external reference metric, and overcomes the usual problems of massive gravity theories (fatal ghosts instabilities, acausality and/or vDVZ discontinuity). The equations of motion of the theory are non-local, but respect causality. The starting point is the quadratic action proposed in the context of the degravitation idea. We show that it is possible to extended it to a fully non-linear covariant theory. This theory describes the five degrees of freedom of a massive graviton plus a scalar ghost. However, contrary to generic non-linear extensions of Fierz-Pauli massive gravity, the ghost has the same mass m as the massive graviton, independently of the background, and smoothly goes into a non-radiative degree of freedom for m-> 0. As a consequence, for $m\sim H_0$ the vacuum instability induced by the ghost is irrelevant even over cosmological time-scales. We finally show that an extension of the model degravitates a vacuum energy density of order $M_{Planck}^4$ down to a value of order $M_{Planck}^2 m^2$, which for $m \sim H_0$ is of order of the observed value of the vacuum energy density.

We re-examine the classic problem of the renormalization of zero-point quantum fluctuations in a Friedmann-Robertson-Walker background. We discuss a number of issues that arise when regularizing the theory with a momentum-space cutoff, and show explicitly how introducing non-covariant counter-terms allows to obtain covariant results for the renormalized vacuum energy-momentum tensor. We clarify some confusion in the literature concerning the equation of state of vacuum fluctuations. Further, we point out that the general structure of the effective action becomes richer if the theory contains a scalar field phi with mass m smaller than the Hubble parameter H(t). Such an ultra-light particle cannot be integrated out completely to get the effective action. Apart from the volume term and the Einstein-Hilbert term, that are reabsorbed into renormalizations of the cosmological constant and Newton's constant, the effective action in general also has a term proportional to F(phi)R, for some function F(phi). As a result, vacuum fluctuations of ultra-light scalar fields naturally lead to models where the dark energy density has the form rho_{DE}(t)=rho_X(t)+rho_Z(t), where rho_X is the component that accelerates the Hubble expansion at late times and rho_Z(t) is an extra contribution proportional to H^2(t). We perform a detailed comparison of such models with CMB, SNIa and BAO data.

We examine a cosmological model with a dark energy density of the form $\rho_{DE}(t)=\rho_X(t)+\rho_Z(t)$, where $\rho_X$ is the component that accelerates the Hubble expansion at late times and $\rho_Z(t)$ is an extra contribution proportional to $H^2(t)$. This form of $\rho_Z(t)$ follows from the recent proposal that the contribution of zero-point fluctuations of quantum fields to the total energy density should be computed by subtracting the Minkowski-space result from that computed in the FRW space-time. We discuss theoretical arguments that support this subtraction. By definition, this eliminates the quartic divergence in the vacuum energy density responsible for the cosmological constant problem. We show that the remaining quadratic divergence can be reabsorbed into a redefinition of Newton's constant only under the assumption that the energy-momentum tensor of vacuum fluctuations is conserved in isolation. However, in the presence of an ultra-light scalar field $X$ with $m_X<H_0$, as typical of some dark energy models, the gravity effective action depends both on the gravitational field and on the $X$ field. In this case general covariance only requires the conservation of the total energy-momentum tensor, including both the classical term $T^X_{\mu\nu}$ and the vacuum expectation value of T_{\mu\nu}. If there is an exchange of energy between these two terms, there are potentially observable consequences. We construct an explicit model with an interaction between $\rho_X$ and $\rho_Z$ and we show that the total dark energy density $\rho_{DE}(t)=\rho_X(t)+\rho_Z(t)$ always remains a finite fraction of the critical density at any time, providing a specific model of early dark energy. We discuss the implication of this result for the coincidence problem and we estimate the model parameters by means of a full likelihood analysis using current CMB, SNe Ia and BAO data.