Abstract: We study the impact of non-standard momentum distributions of cosmic neutrinos on the anisotropy spectrum of the cosmic microwave background and the matter power spectrum of the large scale structure. We show that the neutrino distribution has almost no unique observable imprint, as it is almost entirely degenerate with the effective number of neutrino flavours, N-eff, and the neutrino mass, m(nu). Performing a Markov chain Monte Carlo analysis with current cosmological data, we demonstrate that the neutrino mass bound heavily depends on the assumed momentum distribution of relic neutrinos. The message of this work is simple and has to our knowledge not been pointed out clearly before: cosmology allows that neutrinos have larger masses if their average momentum is larger than that of a perfectly thermal distribution. Here we provide an example in which the mass limits are relaxed by a factor of two.

Abstract: We show that the reheating temperature of a matter-domination era in the early universe can be pushed down to the neutrino decoupling temperature at around 2 MeV if the reheating takes place through non-hadronic decays of the dominant matter and neutrino-antineutrino asymmetries are still large enough, vertical bar L vertical bar greater than or similar to O(10(-2)) (depending on the neutrino flavor) at the end of reheating.

Abstract: We revisit the topic of invisible neutrino decay in the precision cosmological context, via a first-principles approach to understanding the cosmic microwave background and large-scale structure phenomenology of such a non-standard physics scenario. Assuming an effective Lagrangian in which a heavier standard-model neutrino nu(H) couples to a lighter one nu(l) and a massless scalar particle phi via a Yukawa interaction, we derive from first principles the complete set of Boltzmann equations, at both the spatially homogeneous and the firstorder inhomogeneous levels, for the phase space densities of nu(H), nu(l), and phi in the presence of the relevant decay and inverse decay processes. With this set of equations in hand, we perform a critical survey of recent works on cosmological invisible neutrino decay in both limits of decay while nu(H) is ultra-relativistic and non-relativistic. Our two main findings are: (i) in the non-relativistic limit, the effective equations of motion used to describe perturbations in the neutrino-scalar system in the existing literature formally violate momentum conservation and gauge invariance, and (ii) in the ultra-relativistic limit, exponential damping of the anisotropic stress does not occur at the commonly-used rate Gamma(T) = (1/tau(0))( m(nu H)/E-nu H)(3), but at a rate similar to (1/ tau(0))(m(nu H)/E-nu H)(5). Both results are model-independent. The impact of the former finding on the cosmology of invisible neutrino decay is likely small. The latter, however, implies a significant revision of the cosmological limit on the neutrino lifetime tau(0) from tau(old)(0) greater than or similar to 1.2 x 10(9) s (m(nu H)/50 meV)(3) to tau(0) greater than or similar to (4 x 10(5) -> 4 x 10(6)) s (m(nu H)/50 meV)(5).

Abstract: We present in this work a new calculation of the standard-model benchmark value for the effective number of neutrinos, N-eff(SM), that quantifies the cosmological neutrinoto-photon energy densities. The calculation takes into account neutrino flavour oscillations, finite-temperature effects in the quantum electrodynamics plasma to O(e(3)), where e is the elementary electric charge, and a full evaluation of the neutrino-neutrino collision integral. We provide furthermore a detailed assessment of the uncertainties in the benchmark N(eff)(SM )value, through testing the value's dependence on (i) optional approximate modelling of the weak collision integrals, (ii) measurement errors in the physical parameters of the weak sector, and (iii) numerical convergence, particularly in relation to momentum discretisation. Our new, recommended standard-model benchmark is N-eff(SM) 3.0440 +/- 0.0002, where the nominal uncertainty is attributed predominantly to errors incurred in the numerical solution procedure (vertical bar delta N-eff vertical bar similar to 10(-4)), augmented by measurement errors in the solar mixing angle sin(2) theta(12) (vertical bar delta N-eff vertical bar similar to 10(-4)).

Abstract: We investigate the non-linear evolution of the relic cosmic neutrino background by running large box-size, high resolution N-body simulations which incorporate cold dark matter (CDM) and neutrinos as independent particle species. Our set of simulations explore the properties of neutrinos in a reference Lambda CDM model with total neutrino masses between 0.05-0.60 eV in cold dark matter haloes of mass 10(11) – 10(15) h(-1) M-circle dot, over a redshift range z = 0 – 2. We compute the halo mass function and show that it is reasonably well fitted by the Sheth-Tormen formula, once the neutrino contribution to the total matter is removed. More importantly, we focus on the CDM and neutrino properties of the density and peculiar velocity fields in the cosmological volume, inside and in the outskirts of virialized haloes. The dynamical state of the neutrino particles depends strongly on their momentum: whereas neutrinos in the low velocity tail behave similarly to CDM particles, neutrinos in the high velocity tail are not affected by the clustering of the underlying CDM component. We find that the neutrino (linear) unperturbed momentum distribution is modified and mass and redshift dependent deviations from the expected Fermi-Dirac distribution are in place both in the cosmological volume and inside haloes. The neutrino density profiles around virialized haloes have been carefully investigated and a simple fitting formula is provided. The neutrino profile, unlike the cold dark matter one, is found to be cored with core size and central density that depend on the neutrino mass, redshift and mass of the halo, for halos of masses larger than similar to 10(13.5) h(-1) M-circle dot. For lower masses the neutrino profile is best fitted by a simple power-law relation in the range probed by the simulations. The results we obtain are numerically converged in terms of neutrino profiles at the 10% level for scales above similar to 200 h(-1) kpc at z = 0, and are stable with respect to box-size and starting redshift of the simulation. Our findings are particularly important in view of upcoming large-scale structure surveys, like Euclid, that are expected to probe the non-linear regime at the percent level with lensing and clustering observations.