New Clues About Light Sterile Neutrinos: Preference for Models with Damping Effects in Global Fits

Abstract: This article reports global fits of short-baseline neutrino data to oscillation models involving light sterile neutrinos. In the commonly-used 3+1 plane wave model, there is a well-known 4.9$\sigma$ tension between data sets sensitive to appearance versus disappearance of neutrinos. We find that models that damp the oscillation prediction for the reactor data sets, especially at low energy, substantially improve the fits and reduce the tension. We consider two such scenarios. The first scenario introduces the quantum mechanical wavepacket effect that accounts for the source size in reactor experiments into the 3+1 model. We find that inclusion of the wavepacket effect greatly improves the overall fit compared to a 3$\nu$ model by $\Delta \chi^2/$DOF$=61.1/4$ ($7.1\sigma$ improvement) with best-fit $\Delta m^2=1.4$ eV$^2$ and wavepacket length of 67fm. The internal tension is reduced to 3.4$\sigma$. If reactor-data only is fit, then the wavepacket preferred length is 91 fm ($>20$ fm at 99\% CL). The second model introduces oscillations involving sterile flavor and allows the decay of the heaviest, mostly sterile mass state, $\nu_4$. This model introduces a damping term similar to the wavepacket effect, but across all experiments. Compared to a three-neutrino fit, this has a $\Delta \chi^2/$DOF$=60.6/4$ ($7\sigma$ improvement) with preferred $\Delta m^2=1.4$ eV$^2$ and decay $\Gamma = 0.35$ eV$^2$. The internal tension is reduced to 3.7$\sigma$. For many years, the reactor event rates have been observed to have structure that deviates from prediction. Community discussion has focused on an excess compared to prediction observed at 5 MeV; however, other deviations are apparent. This structure has $L$ dependence that is well-fit by the damped models. Before assuming this points to new physics, we urge closer examination of systematic effects that could lead to this $L$ dependence.