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: We present a full Bayesian analysis of the combination of current neutrino oscillation, neutrinoless double beta decay and Cosmic Microwave Background observations. Our major goal is to carefully investigate the possibility to single out one neutrino mass ordering, namely Normal Ordering or Inverted Ordering, with current data. Two possible parametrizations (three neutrino masses versus the lightest neutrino mass plus the two oscillation mass splittings) and priors (linear versus logarithmic) are exhaustively examined. We find that the preference for NO is only driven by neutrino oscillation data. Moreover, the values of the Bayes factor indicate that the evidence for NO is strong only when the scan is performed over the three neutrino masses with logarithmic priors; for every other combination of parameterization and prior, the preference for NO is only weak. As a by-product of our Bayesian analyses, we are able to (a) compare the Bayesian bounds on the neutrino mixing parameters to those obtained by means of frequentist approaches, finding a very good agreement; (b) determine that the lightest neutrino mass plus the two mass splittings parametrization, motivated by the physical observables, is strongly preferred over the three neutrino mass eigenstates scan and (c) find that logarithmic priors guarantee a weakly-to-moderately more efficient sampling of the parameter space. These results establish the optimal strategy to successfully explore the neutrino parameter space, based on the use of the oscillation mass splittings and a logarithmic prior on the lightest neutrino mass, when combining neutrino oscillation data with cosmology and neutrinoless double beta decay. We also show that the limits on the total neutrino mass Sigma m(nu) can change dramatically when moving from one prior to the other. These results have profound implications for future studies on the neutrino mass ordering, as they crucially state the need for self-consistent analyses which explore the best parametrization and priors, without combining results that involve different assumptions.

Abstract: In this paper, we explore a scenario where the dark matter is a mixture of interacting and non interacting species. Assuming dark matter-photon interactions for the interacting species, we find that the suppression of the matter power spectrum in this scenario can mimic that expected in the case of massive neutrinos. Our numerical studies include present limits from Planck Cosmic Microwave Background data, which render the strength of the dark matter photon interaction unconstrained when the fraction of interacting dark matter is small. Despite the large entangling between mixed dark matter and neutrino masses, we show that future measurements from the Dark Energy Instrument (DESI) could help in establishing the dark matter and the neutrino properties simultaneously, provided that the interaction rate is very close to its current limits and the fraction of interacting dark matter is at least of O (10%). However, for that region of parameter space where a small fraction of interacting DM coincides with a comparatively large interaction rate, our analysis highlights a considerable degeneracy between the mixed dark matter parameters and the neutrino mass scale.

Abstract: The extraction of the neutrino mass ordering is one of the major challenges in particle physics and cosmology, not only for its implications for a fundamental theory of mass generation in nature, but also for its decisive role in the scale of future neutrinoless double beta decay experimental searches. It has been recently claimed that current oscillation, beta decay and cosmological limits on the different observables describing the neutrino mass parameter space provide robust decisive Bayesian evidence in favor of the normal ordering of the neutrino mass spectrum [1]. We further investigate these strong claims using a rich and wide phenomenology, with different sampling techniques of the neutrino parameter space. Contrary to the findings of Jimenez et al. [1], no decisive evidence for the normal mass ordering is found. Neutrino mass ordering analyses must rely on priors and parameterizations that are ordering-agnostic: robust results should be regarded as those in which the preference for the normal neutrino mass ordering is driven exclusively by the data, while we find a difference of up to a factor of 33 in the Bayes factors among the different priors and parameterizations exploited here. An ordering-agnostic prior would be represented by the case of parameterizations sampling over the two mass splittings and a mass scale, or those sampling over the individual neutrino masses via normal prior distributions only. In this regard, we show that the current significance in favor of the normal mass ordering should be taken as 2.7 sigma (i.e. moderate evidence), mostly driven by neutrino oscillation data. Let us stress that, while current data favor NO only mildly, we do not exclude the possibility that this may change in the future. Eventually, upcoming oscillation and cosmological data may (or may not) lead to a more significant exclusion of IO.

Abstract: Current neutrino experiments are measuring the neutrino mixing parameters with an unprecedented accuracy. The upcoming generation of neutrino experiments will be sensitive to subdominant neutrino oscillation effects that can in principle give information on the yet-unknown neutrino parameters: the Dirac CP-violating phase in the PMNS mixing matrix, the neutrino mass ordering and the octant of.23. Determining the exact values of neutrino mass and mixing parameters is crucial to test various neutrino models and flavor symmetries that are designed to predict these neutrino parameters. In the first part of this review, we summarize the current status of the neutrino oscillation parameter determination. We consider the most recent data from all solar neutrino experiments and the atmospheric neutrino data from Super-Kamiokande, IceCube, and ANTARES. We also implement the data from the reactor neutrino experiments KamLAND, Daya Bay, RENO, and Double Chooz as well as the long baseline neutrino data from MINOS, T2K, and NO.A. If in addition to the standard interactions, neutrinos have subdominant yet-unknown Non-Standard Interactions (NSI) with matter fields, extracting the values of these parameters will suffer from new degeneracies and ambiguities. We review such effects and formulate the conditions on the NSI parameters under which the precision measurement of neutrino oscillation parameters can be distorted. Like standard weak interactions, the non-standard interaction can be categorized into two groups: Charged Current (CC) NSI and Neutral Current (NC) NSI. Our focus will bemainly on neutral current NSI because it is possible to build a class of models that give rise to sizeable NC NSI with discernible effects on neutrino oscillation. These models are based on new U(1) gauge symmetry with a gauge boson of mass. 10 MeV. The UV complete model should be of course electroweak invariant which in general implies that along with neutrinos, charged fermions also acquire new interactions on which there are strong bounds. We enumerate the bounds that already exist on the electroweak symmetric models and demonstrate that it is possible to build viable models avoiding all these bounds. In the end, we review methods to test these models and suggest approaches to break the degeneracies in deriving neutrino mass parameters caused by NSI.