Abstract: Sterile neutrinos with a mass in the eV range have been invoked as a possible explanation of a variety of short baseline (SBL) neutrino oscillation anomalies. However, if one considers neutrino oscillations between active and sterile neutrinos, such neutrinos would have been fully thermalised in the early universe, and would be therefore in strong conflict with cosmological bounds. In this study we first update cosmological bounds on the mass and energy density of eV-scale sterile neutrinos. We then perform an updated study of a previously proposed model in which the sterile neutrino couples to a new light pseudoscalar degree of freedom. Consistently with previous analyses, we find that the model provides a good fit to all cosmological data and allows the high value of H-0 measured in the local universe to be consistent with measurements of the cosmic microwave background. However, new high l polarisation data constrain the sterile neutrino mass to be less than approximately 1 eV in this scenario. Finally, we combine the cosmological bounds on the pseudoscalar model with a Bayesian inference analysis of SBL data and conclude that only a sterile mass in narrow ranges around 1 eV remains consistent with both cosmology and SBL data.

Abstract: The atmospheric neutrino passing fraction, or self-veto, is defined as the probability for an atmospheric neutrino not to be accompanied by a detectable muon from the same cosmic-ray air shower. Building upon previous work, we propose a redefinition of the passing fractions by unifying the treatment for muon and electron neutrinos. Several approximations have also been removed. This enables performing detailed estimations of the uncertainties in the passing fractions from several inputs: muon losses, cosmic-ray spectrum, hadronic-interaction models and atmosphere-density profiles. We also study the passing fractions under variations of the detector configuration: depth, surrounding medium and muon veto trigger probability. The calculation exhibits excellent agreement with passing fractions obtained from Monte Carlo simulations. Finally, we provide a general software framework to implement this veto technique for all large-scale neutrino observatories.

Abstract: We revisit the topic of invisible neutrino decay in the precision cosmological context, via a first-principles approach to understanding the cosmic microwave background and large-scale structure phenomenology of such a non-standard physics scenario. Assuming an effective Lagrangian in which a heavier standard-model neutrino nu(H) couples to a lighter one nu(l) and a massless scalar particle phi via a Yukawa interaction, we derive from first principles the complete set of Boltzmann equations, at both the spatially homogeneous and the firstorder inhomogeneous levels, for the phase space densities of nu(H), nu(l), and phi in the presence of the relevant decay and inverse decay processes. With this set of equations in hand, we perform a critical survey of recent works on cosmological invisible neutrino decay in both limits of decay while nu(H) is ultra-relativistic and non-relativistic. Our two main findings are: (i) in the non-relativistic limit, the effective equations of motion used to describe perturbations in the neutrino-scalar system in the existing literature formally violate momentum conservation and gauge invariance, and (ii) in the ultra-relativistic limit, exponential damping of the anisotropic stress does not occur at the commonly-used rate Gamma(T) = (1/tau(0))( m(nu H)/E-nu H)(3), but at a rate similar to (1/ tau(0))(m(nu H)/E-nu H)(5). Both results are model-independent. The impact of the former finding on the cosmology of invisible neutrino decay is likely small. The latter, however, implies a significant revision of the cosmological limit on the neutrino lifetime tau(0) from tau(old)(0) greater than or similar to 1.2 x 10(9) s (m(nu H)/50 meV)(3) to tau(0) greater than or similar to (4 x 10(5) -> 4 x 10(6)) s (m(nu H)/50 meV)(5).