In this paper we show that the Germani-Kehagias model of Higgs inflation (or New Higgs Inflation), where the Higgs boson is kinetically non-minimally coupled to the Einstein tensor is in perfect compatibility with the latest Planck and BICEP2 data. Moreover, we show that the tension between the Planck and BICEP2 data can be relieved within the New Higgs inflation scenario by a negative running of the spectral index. Regarding the unitarity of the model, we argue that it is unitary throughout the evolution of the Universe. Weak couplings in the Higgs-Higgs and Higgs-graviton sectors are provided by a large background dependent cut-off scale during inflation. In the same regime, the W and Z gauge bosons acquire a very large mass, thus decouple. On the other hand, if they are also non-minimally coupled to the Higgs boson, their effective masses can be enormously reduced. In this case, the W and Z bosons are no longer decoupled. After inflation, the New Higgs model is well approximated by a quartic Galileon with a renormalizable potential. We argue that this can unitarily create the right conditions for inflation to eventually start.

We consider general theories of a massive spin-2 particle $h_{\mu\nu}$ on a Minkowski background. A decomposition of $h_{\mu\nu}$ in terms of helicity eigenstates allows us to directly test whether any given theory possesses a consistent description as a massive spin-2 representation of the Poincar\'e group. We demonstrate (i) that any nonlinear theory with an Einsteinian derivative structure either contains ghosts or does not describe a weakly coupled spin-2 and (ii) that there exists a two-parameter family of non-Einsteinian cubic self-interactions which constitute a ghost-free massive spin-2 theory.

The spherical collapse model is often used to follow the evolution of overdensities into the nonlinear regime. We describe the correct approach to be used in coupled dark energy cosmologies, where a fifth force, different from gravity and mediated by the dark energy scalar field, influences the collapse. We reformulate the spherical collapse description by deriving it directly from the set of nonlinear hydrodynamical Navier Stokes equations. By comparing with the corresponding relativistic equations, we show how the fifth force should be taken into account within the spherical collapse picture and clarify the problems arising when an inhomogeneous scalar field is considered within a spherical collapse picture. We then apply our method to the case of coupled quintessence, where the fifth force acts among cold dark matter particles, and to growing neutrino quintessence, where the fifth force acts between neutrinos. Furthermore, we review this method when applied to standard cosmologies and apply our analysis to minimally coupled quintessence and check past results for early dark energy parametrizations.

The interaction between the cosmon and neutrinos may solve the "why now problem" for dark energy cosmologies. Within growing neutrino quintessence it leads to the formation of nonlinear neutrino lumps. For a test of such models by the integrated Sachs-Wolfe effect for the cosmic microwave background (CMB) we estimate the size and time evolution of the gravitational potential induced by these lumps. A population of lumps with size of 100 Mpc or more could lead to observable effects on the CMB anisotropies for low angular momenta. The linear approximation is found to be invalid for the relevant length scales. Quantitative estimates depend strongly on the details of the transition between the linear and nonlinear regimes. In particular, important backreaction effects arise from the nonlinearities of the cosmon interactions. At the present stage the uncertainties of the estimate make it difficult to constrain the parameter space of growing neutrino models. We explicitly discuss scenarios and models that are compatible with the CMB observations.

A quintessence scalar field or cosmon interacting with neutrinos can have important effects on cosmological structure formation. Within growing neutrino models the coupling becomes effective only in recent times, when neutrinos become non-relativistic, stopping the evolution of the cosmon. This can explain why dark energy dominates the universe only in a rather recent epoch by relating the present dark energy density to the small mass of neutrinos. Such models predict the presence of stable neutrino lumps at supercluster scales (~200 Mpc and bigger), caused by an attractive force between neutrinos which is stronger than gravity and mediated by the cosmon. We present a method to follow the initial non-linear formation of neutrino lumps in physical space, by integrating numerically on a 3D grid non-linear evolution equations, until virialization naturally occurs. As a first application, we show results for cosmologies with final large neutrino average mass ~2 eV: in this case, neutrino lumps indeed form and mimic very large cold dark matter structures, with a typical gravitational potential 10^{-5} for a lump size ~10 Mpc, and reaching larger values for lumps of about 200 Mpc. A rough estimate of the cosmological gravitational potential at small k in the non-linear regime, Phi_nu = 10^{-6} (k/k_0)^{-2}, 1.2x10^{-2} h/Mpc < k_0 < 7.8x10^{-2} h/Mpc, turns out to be many orders of magnitude smaller than an extrapolation of the linear evolution of density fluctuations. The size of the neutrino-induced gravitational potential could modify the spectrum of CMB anisotropies for small angular momenta.