Combining sterile neutrino fits to short-baseline data with IceCube data

Recent global fits to short-baseline neutrino oscillation data have been performed finding preference for a sterile neutrino solution (3+1) over null. In the most recent iteration, it was pointed out that an unstable sterile neutrino (3+1+decay) may be a better description of the data. This is due to the fact that this model significantly reduces the tension between appearance and disappearance datasets. In this work, we add a 1-year IceCube dataset to the global fit obtaining new results for the standard 3+1 and 3+1+decay sterile neutrino scenarios. We find that the 3+1+decay model provides a better fit than the 3+1, even in the presence of IceCube, with reduced appearance to disappearance tension. The 3+1+decay model is a 5.4σ improvement over the null hypothesis and a 2.8σ improvement over the standard 3+1 model.