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: Cosmological observables show a dependence with the neutrino mass, which is partially degenerate with parameters of extended models of gravity. We study and explore this degeneracy in Horndeski generalized scalar-tensor theories of gravity. Using forecasted cosmic microwave background and galaxy power spectrum datasets, we find that a single parameter in the linear regime of the effective theory dominates the correlation with the total neutrino mass. For any given mass, a particular value of this parameter approximately cancels the power suppression due to the neutrino mass at a given redshift. The extent of the cancellation of this degeneracy depends on the cosmological large-scale structure data used at different redshifts. We constrain the parameters and functions of the effective gravity theory and determine the influence of gravity on the determination of the neutrino mass from present and future surveys.

Abstract: We have studied the consequences of breaking the CPT symmetry in the neutrino sector, using the expected high-energy neutrino flux from distant cosmological sources such as active galaxies. For this purpose we have assumed three different hypotheses for the neutrino production model, characterised by the flavour fluxes at production phi(0)(e) : phi(0)(mu) : phi(0)(tau) = 1 : 2 : 0, 0 : 1 : 0, and 1 : 0 : 0, and studied the theoretical and experimental expectations for the muon-neutrino flux at Earth, phi(mu), and for the flavour ratios at Earth, R = phi(mu)/phi(e) and S = phi(tau)/phi(mu). CPT violation (CPTV) has been implemented by adding an energy-independent term to the standard neutrino oscillation Hamiltonian. This introduces three new mixing angles, two new eigenvalues and three new phases, all of which have currently unknown values. We have varied the new mixing angles and eigenvalues within certain bounds, together with the parameters associated to pure standard oscillations. Our results indicate that, for the models 1 : 2 : 0 and 0 : 1 : 0, it might be possible to find large deviations of phi(mu), R, and S between the cases without and with CPTV, provided the CPTV eigenvalues lie within 10(-29) – 10(-27) GeV, or above. Moreover, if CPTV exists, there are certain values of R and S that can be accounted for by up to three production models. If no CPTV were observed, we could set limits on the CPTV eigenvalues of the same order. Detection prospects calculated using IceCube suggest that for the models 1 : 2 : 0 and 0 : 1 : 0, the modifications due to CPTV are larger and more clearly separable from the standard-oscillations predictions. We conclude that IceCube is potentially able to detect CPTV but that, depending on the values of the CPTV parameters, there could be a mis-determination of the neutrino production model.