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: We present the results of an updated fit of short-baseline neutrino oscillation data in the framework of 3+1 active-sterile neutrino mixing. We first consider v(e) and (v) over bar (e) disappearance in the light of the Gallium and reactor anomalies. We discuss the implications of the recent measurement of the reactor (v) over bar (e) spectrum in the NEOS experiment, which shifts the allowed regions of the parameter space towards smaller values of |U-e1|(2). The beta-decay constraints of the Mainz and Troitsk experiments allow us to limit the oscillation length between about 2 cm and 7 m at 3 sigma for neutrinos with an energy of 1 MeV. The corresponding oscillations can be discovered in a model-independent way in ongoing reactor and source experiments by measuring v(e) and (v) over bar (e), disappearance as a function of distance. We then consider the global fit of the data on short-baseline v(mu)((-)) -> v(e)((-)) transitions in the light of the LSND anomaly, taking into account the constraints from v(e)(( )) and v(mu)((-)) disappearance experiments, including the recent data of the MINOS and IceCube experiments. The combination of the NEOS constraints on |U-e4|(2) and the MINOS and IceCube constraints on |U-mu 4|(2) lead to an unacceptable appearance-disappearance tension which becomes tolerable only in a pragmatic fit which neglects the MiniBooNE low-energy anomaly. The minimization of the global chi(2) in the space of the four mixing parameters Delta m(41)(2), |U-e4|(2), |U-mu 4|(2) and |U-4 tau|(2) leads to three allowed regions with narrow Delta m(41)(2) widths at Delta m(41)(2) approximate to 1.7 (best-fit), 1.3 (at 2 sigma), 2.4 (at 3 sigma) eV(2). The effective amplitude of short-baseline v(mu)((-)) -> v(e)((-)) oscillations is limited by 0.00048 less than or similar to sin(2) 2 nu(e mu) less than or similar to 0.0020 at 3 sigma The restrictions of the allowed regions of the mixing parameters with respect to our previous global fits are mainly due to the NEOS constraints. We present a comparison of the allowed regions of the mixing parameters with the sensitivities of ongoing experiments, which show that it is likely that these experiments will determine in a definitive way if the reactor, Gallium and LSND anomalies are due to active-sterile neutrino oscillations or not.

Abstract: We present the cosmological bounds on the thermal axion mass in an extended cosmological scenario in which the primordial power spectrum of scalar perturbations differs from the usual power-law shape predicted by the simplest inflationary models. The power spectrum is instead modeled by means of a “piecewise cubic Hermite interpolating polynomial” (PCHIP). When using cosmic microwave background measurements combined with other cosmological data sets, the thermal axion mass constraints are degraded only slightly. The addition of the measurements of sigma(8) and Omega(m) from the 2013 Planck cluster catalog on galaxy number counts relaxes the bounds on the thermal axion mass, mildly favoring a similar to 1 eV axion mass, regardless of the model adopted for the primordial power spectrum. However, in general, such a preference disappears if the sum of the three active neutrino masses is also considered as a free parameter in our numerical analyses, due to the strong correlation between the masses of these two hot thermal relics.