Abstract: We present a consistent framework to set limits on properties of light sterile neutrinos coupled to all three active neutrinos using a combination of the latest cosmological data and terrestrial measurements from oscillations, beta-decay, and neutrinoless double-beta-decay (0 nu beta beta) experiments. We directly constrain the full 3 + 1 active-sterile mixing matrix elements vertical bar U-alpha 4 vertical bar(2) , with alpha is an element of (e,mu,tau), and the mass-squared splitting Delta m(41)(2) (math) m(4)(2) – m(1)(2). We find that results for a 3 + 1 case differ from previously studied 1 + 1 scenarios where the sterile is coupled to only one of the neutrinos, which is largely explained by parameter space volume effects. Limits on the mass splitting and the mixing matrix elements are currently dominated by the cosmological datasets. The exact results are slightly prior dependent, but we reliably find all matrix elements to be constrained below vertical bar U-alpha 4 vertical bar(2) less than or similar to 10(-3) . Short-baseline neutrino oscillation hints in favor of eV-scale sterile neutrinos arc in serious tension with these bounds, irrespective of prior assumptions. We also translate the bounds from the cosmological analysis into constraints on the parameters probed by laboratory searches, such as m(beta) or m(beta)(beta), the effective mass parameters probed by beta-decay and 0 nu beta beta searches, respectively. When allowing for mixing with a light sterile neutrino, cosmology leads to upper bounds of m(beta) < 0.09 eV and m(beta)(beta )< 0.07 eV at 95% CL, more stringent than the limits from current laboratory experiments.

Abstract: We present an updated global fit of neutrino oscillation data in the simplest three-neutrino framework. In the present study we include up-to-date analyses from a number of experiments. Concerning the atmospheric and solar sectors, besides the data considered previously, we give updated analyses of IceCube DeepCore and Sudbury Neutrino Observatory data, respectively. We have also included the latest electron antineutrino data collected by the Daya Bay and RENO reactor experiments, and the long-baseline T2K and NO nu A measurements, as reported in the Neutrino 2020 conference. All in all, these new analyses result in more accurate measurements of theta (13), theta (12), Delta m212 and Delta m312. The best fit value for the atmospheric angle theta (23) lies in the second octant, but first octant solutions remain allowed at similar to 2.4 sigma. Regarding CP violation measurements, the preferred value of delta we obtain is 1.08 pi (1.58 pi) for normal (inverted) neutrino mass ordering. The global analysis still prefers normal neutrino mass ordering with 2.5 sigma statistical significance. This preference is milder than the one found in previous global analyses. These new results should be regarded as robust due to the agreement found between our Bayesian and frequentist approaches. Taking into account only oscillation data, there is a weak/moderate preference for the normal neutrino mass ordering of 2.00 sigma. While adding neutrinoless double beta decay from the latest Gerda, CUORE and KamLAND-Zen results barely modifies this picture, cosmological measurements raise the preference to 2.68 sigma within a conservative approach. A more aggressive data set combination of cosmological observations leads to a similar preference for normal with respect to inverted mass ordering, namely 2.70 sigma. This very same cosmological data set provides 2 sigma upper limits on the total neutrino mass corresponding to Sigma m(nu)< 0.12 (0.15) eV in the normal (inverted) neutrino mass ordering scenario. The bounds on the neutrino mixing parameters and masses presented in this up-to-date global fit analysis include all currently available neutrino physics inputs.