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: We present a new global fit of neutrino oscillation parameters within the simplest three-neutrino picture, including new data which appeared since our previous analysis[1]. In this update we include new long-baseline neutrino data involving the antineutrino channel in T2K, as well as new data in the neutrino channel, data from NO nu A, as well as new reactor data, such as the Daya Bay 1230 days electron antineutrino disappearance spectrum data and the 1500 live days prompt spectrum from RENO, as well as new Double Chooz data. We also include atmospheric neutrino data from the IceCube DeepCore and ANTARES neutrino telescopes and from Super-Kamiokande. Finally, we also update our solar oscillation analysis by including the 2055-day day/night spectrum from the fourth phase of the Super-Kamiokande experiment. With the new data we find a preference for the atmospheric angle in the upper octant for both neutrino mass orderings, with maximal mixing allowed at Delta chi(2)= 1.6 (3.2) for normal (inverted) ordering. We also obtain a strong preference for values of the CP phase delta in the range [pi, 2 pi], excluding values close to pi/2at more than 4 sigma. More remarkably, our global analysis shows a hint in favorof the normal mass ordering over the inverted one at more than 3 sigma. We discuss in detail the status of the mass ordering, CP violation and octant sensitivities, analyzing the interplay among the different neutrino data samples.

Abstract: We propose a realization of the so-called bimodal/schizophrenic model proposed recently. We assume 54, the permutation group of four objects as flavor symmetry giving tri-bimaximal lepton mixing at leading order. In these models the second massive neutrino state is assumed quasi-Dirac and the remaining neutrinos are Majorana states. In the case of inverse mass hierarchy, the lower bound on the neutrinoless double beta decay parameter m(ee) is about two times that of the usual lower bound, within the range of sensitivity of the next generation of experiments.

Abstract: We revisit the joint constraints in the mixed hot dark matter scenario in which both thermally produced QCD axions and relic neutrinos are present. Upon recomputing the cosmological axion abundance via recent advances in the literature, we improve the state-of-the-art analyses and provide updated bounds on axion and neutrino masses. By avoiding approximate methods, such as the instantaneous decoupling approximation, and limitations due to the limited validity of the perturbative approach in QCD that forced to artificially divide the constraints from the axion-pion and the axion-gluon production channels, we find robust and self-consistent limits. We investigate the two most popular axion frameworks: KSVZ and DFSZ. From Big Bang Nucleosynthesis (BBN) light element abundances data we find for the KSVZ axion Delta N-eff < 0.31 and an axion mass bound m(a) < 0.53 eV (i.e., a bound on the axion decay constant f(a) > 1.07 x 10(7) GeV) both at 95% CL. These BBN bounds are improved to Delta N-eff < 0.14 and m(a) < 0.16 eV (f(a) > 3.56 x 10(7) GeV) if a prior on the baryon energy density from Cosmic Microwave Background (CMB) data is assumed. When instead considering cosmological observations from the CMB temperature, polarization and lensing from the Planck satellite combined with large scale structure data we find Delta N-eff < 0.23, m(a) < 0.28 eV (f(a) > 2.02 x 10(7) GeV) and Sigma m(nu) < 0.16 eV at 95% CL. This corresponds approximately to a factor of 5 improvement in the axion mass bound with respect to the existing limits. Very similar results are obtained for the DFSZ axion. We also forecast upcoming observations from future CMB and galaxy surveys, showing that they could reach percent level errors for m(a) similar to 1 eV.