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.

Abstract: Recent Cosmic Microwave Background (CMB) results from the Planck satellite, combined with previous CMB data and Hubble constant measurements from the Hubble Space Telescope, provide a constraint on the effective number of relativistic degrees of freedom 3.62(-0.48)(+0.50) at 95% CL. New Planck data provide a unique opportunity to place limits on models containing relativistic species at the decoupling epoch. We present here the bounds on sterile neutrino models combining Planck data with galaxy clustering information. Assuming N-eff active plus sterile massive neutrino species, in the case of a Planck+WP+HighL+HST analysis we find m(nu,sterile)(eff) < 0.36 eV and 3.14 < N-eff < 4.15 at 95% CL, while using Planck+WP+HighL data in combination with the full shape of the galaxy power spectrum from the Baryon Oscillation Spectroscopic Survey BOSS Data Relase 9 measurements, we find that 3.30 < N-eff < 4.43 and m(nu,sterile)(eff) < 0.33 eV both at 95% CL with the three active neutrinos having the minimum mass allowed in the normal hierarchy scheme, i.e. Sigma m(nu) similar to 0.06 eV. These values compromise the viability of the (3 + 2) massive sterile neutrino models for the parameter region indicated by global fits of neutrino oscillation data. Within the (3 + 1) massive sterile neutrino scenario, we find m(nu,sterile)(eff) < 0.34 eV at 95% CL. While the existence of one extra sterile massive neutrino state is compatible with current oscillation data, the values for the sterile neutrino mass preferred by oscillation analyses are significantly higher than the current cosmological bound. We review as well the bounds on extended dark sectors with additional light species based on the latest Planck CMB observations.

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 study the decoupling process of neutrinos in the early universe in the presence of three-flavour oscillations. The evolution of the neutrino spectra is found by solving the corresponding momentum-dependent kinetic equations for the neutrino density matrix, including for the first time the proper collision integrals for both diagonal and off-diagonal elements. This improved calculation modifies the evolution of the off-diagonal elements of the neutrino density matrix and changes the deviation from equilibrium of the frozen neutrino spectra. However, it does not vary the contribution of neutrinos to the cosmological energy density in the form of radiation, usually expressed in terms of the effective number of neutrinos, N-eff. We find a value of N-eff = 3.045, in agreement with previous theoretical calculations and consistent with the latest analysis of Planck data. This result does not depend on the ordering of neutrino masses. We also consider the effect of non-standard neutrino-electron interactions (NSI), predicted in many theoretical models where neutrinos acquire mass. For two sets of NSI parameters allowed by present data, we find that Neff can be reduced down to 3.040 or enhanced up to 3.059.

Abstract: Recent constrains on the sum of neutrino masses inferred by analyzing cosmological data, show that detecting a non-zero neutrino mass is within reach of forthcoming cosmological surveys. Such a measurement will imply a direct determination of the absolute neutrino mass scale. Physically, the measurement relies on constraining the shape of the matter power spectrum below the neutrino free streaming scale: massive neutrinos erase power at these scales. However, detection of a lack of small-scale power from cosmological data could also be due to a host of other effects. It is therefore of paramount importance to validate neutrinos as the source of power suppression at small scales. We show that, independent on hierarchy, neutrinos always show a footprint on large, linear scales; the exact location and properties are fully specified by the measured power suppression (an astrophysical measurement) and atmospheric neutrinos mass splitting (a neutrino oscillation experiment measurement). This feature cannot be easily mimicked by systematic uncertainties in the cosmological data analysis or modifications in the cosmological model. Therefore the measurement of such a feature, up to 1% relative change in the power spectrum for extreme differences in the mass eigenstates mass ratios, is a smoking gun for confirming the determination of the absolute neutrino mass scale from cosmological observations. It also demonstrates the synergy between astrophysics and particle physics experiments.