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 describe the main features of a new and updated version of the program PArthENoPE, which computes the abundances of light elements produced during Big Bang Nucleosynthesis. As the previous first release in 2008, the new one, PArthENoPE2.0, is publicly available and distributed from the code site, http://parthenope.na.infn.it . Apart from minor changes, which will be also detailed, the main improvements are as follows. The powerful, but not freely accessible, NAG routines have been substituted by ODEPACK libraries, without any significant loss in precision. Moreover, we have developed a Graphical User Interface (GUI) which allows a friendly use of the code and a simpler implementation of running for grids of input parameters. New Version program summary Program Title: PArthENoPE2.0 Program Files doi : http://dx.doi.org/10.17632/wvgr7d8yt9.1 Licensing provisions: GPLv3 Programming language: Fortran 77 and Python Supplementary material: User Manual available on the web page http://parthenope.na.infn.it Journal reference of previous version: Comput. Phys. Commun. 178 (2008) 956 971 Does the new version supersede the previous version?: Yes Reasons for the new version: Make the code more versatile and user friendly Summary of revisions: (1) Publicly available libraries (2) GUI for configuration Nature of problem: Computation of yields of light elements synthesized in the primordial universe Solution method: Livermore Solver for Ordinary Differential Equations (LSODE) for stiff and nonstiff systems

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.