# Sterile Neutrino Fits to Short-Baseline Neutrino Oscillation Measurements

Authors: J. M. Conrad; C. M. Ignarra; G. Karagiorgi; M. H. Shaevitz; J. Spitz
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

This paper reviews short-baseline oscillation experiments as interpreted within the context of one, two, and
three sterile neutrino models associated with additional neutrino mass states in the~1 eV range. Appearance and disappearance signals and limits are considered. We show that fitting short-baseline datasets to a 3 + 3 (3 + 2) model, defined by three active and three (two) sterile neutrinos, results in an overall goodness of fit of 67% (69%) and good compatibility between data sets—to be compared to a 3 + 1 model with a 55% goodness of fit. While the (3 + 3) fit yields the highest quality overall, it still finds
inconsistencies with the MiniBooNE appearance datasets; in particular, the global fit fails to account for the observed MiniBooNE low-energy excess. Given the overall improvement, we recommend using the results of
(3 + 2) and (3 + 3) fits, rather than (3 + 1) fits, for future neutrino oscillation phenomenology. These results motivate the pursuit of further short-baseline experiments, such as those reviewed in this paper.


---

## Body

## 1. Introduction

Over the past 15 years, neutrino oscillations associated with small splittings between the neutrino mass states have become well established [1–16]. Based on this, a phenomenological extension of the Standard Model (SM) has been constructed involving three neutrino mass states, over which the three known flavors of neutrinos (νe, νμ, and ντ) are distributed. This is a minimal extension of the SM requiring a lepton mixing matrix, analogous to the quark sector and introducing neutrino mass.Despite its success, the model does not address fundamental questions such as how neutrino masses should be incorporated into an SM Lagrangian or why the neutrino sector has small masses and large mixing angles compared to the quark sector. As a result, while this structure makes successful predictions, one would like to gain a deeper understanding of neutrino phenomenology. This has led to searches for other unexpected properties of neutrinos that might lead to clues towards a more complete theory governing their behavior.Recalling that the mass splitting is related to the frequency of oscillation, short-baseline (SBL) experiments search for evidence of “rapid” oscillations above the established solar (~10−5 eV2) and atmospheric (~10−3 eV2) mass splittings that are incorporated into today’s framework. A key motivation is the search for light sterile neutrinos-fermions that do not participate in SM interactions but do participate in mixing with the established SM neutrinos. Indications of oscillations between active and sterile neutrinos have been observed in the LSND [17], MiniBooNE [18], and reactor [19] experiments, though many others have contributed additional probes of the effect, which are of comparable sensitivity and/or complementary to those above.This paper examines these results within the context of models describing oscillations with sterile neutrinos. An oscillation formalism that introduces multiple sterile neutrinos is described in the next section. Following this, we review the SBL datasets used in the fits presented in this paper, which include both positive signals and stringent limits. We then detail the analysis approach, which we have developed in a series of past papers [20–22]. The global fits are presented with one, two, and three light sterile neutrinos. While groups [20, 23, 24] have explored fits with two sterile neutrinos in the past, the fits presented here represent an important step forward. In particular, we show that, for the first time, the (3 + 3) model resolves some disagreements between the datasets. Lastly, the future of SBL searches for sterile neutrinos is reviewed.


## 2. Oscillations Involving Sterile Neutrinos

### 2.1. Light Sterile States

Sterile neutrinos are additional states beyond the standard electron, muon, and tau flavors, which do not interact via the exchange ofW or Z bosons [25] and are thus “sterile” with respect to the weak interaction. These states are motivated by many Beyond Standard Model theories, where they are often introduced as gauge singlets. Traditionally, sterile neutrinos were introduced at very high mass scales within the context of grand unification and leptogenesis. For many years, sterile neutrinos with light masses were regarded as less natural. However, as recent data [17, 19, 26, 27] has indicated the potential existence of light sterile neutrinos, the theoretical view has evolved to accommodate these light mass gauge singlets [28, 29]. At this point, it is generally accepted that the mass scale for sterile neutrinos is not well predicted, and the existence of one or more sterile neutrinos accommodated by introducing extra neutrino mass states at the eV scale is possible. An excellent review of the phenomenology of sterile neutrinos, as well as the data motivating light sterile models, is provided in [23].Within the expanded oscillation phenomenology, sterile neutrinos are handled as additional noninteracting flavors, which are connected to additional mass states via an extended mixing matrix with extra mixing angles and CP violating phases. These additional mass states must be mostly sterile, with only a small admixture of the active flavors, in order to accommodate the limits on oscillations to sterile neutrinos from the atmospheric and solar neutrino data. Experimental evidence for these additional mass states would come from the disappearance of an active flavor to a sterile neutrino state or additional transitions from one active flavor to another through the sterile neutrino state.The number of light sterile neutrinos is not predicted by theory. However, a natural tendency is to introduce three sterile states. Depending on how the states are distributed in mass scale, one, two, or all three states may be involved in SBL oscillations. These are referred to as (3+N) models where the “3” refers to the three active flavors and the “N" refers to the number of sterile neutrinos.Introducing sterile neutrinos can have implications in cosmological observations, especially measurements of the radiation density in the early universe. These are compounded if the extra neutrinos have significant mass (>1 eV) and do not decay. Currently, cosmological data allow additional states and in many cases favor light sterile neutrinos [30–36]. Upcoming Planck data [37] is expected to precisely measure Neff. This parameter, however, can be considered a model-dependent one. As an example, there are a variety of classes of theories where the neutrinos do not thermalize in the early universe [23]. In these cases, the cosmological neutrino abundance would substantially decrease, rendering cosmological measurements of Neff invalid. Therefore, while the community certainly looks forward to cosmological measurements of Neff, we think that SBL experiments are a largely better approach for probing light sterile neutrinos and constraining their mixing properties. We therefore proceed with a study of the SBL data, without further reference to the cosmological results.


### 2.2. The Basic Oscillation Formalism

Before considering the phenomenology of light sterile neutrinos, it is useful to introduce the idea of oscillations within a simpler model. In this section, we first consider the two-neutrino formalism. We then extend these ideas to form the well-established three-active-flavor neutrino model. Based on these concepts, we expand the discussion to include more states in the following section.Neutrino oscillations require that (1) neutrinos have mass; (2) the difference between the masses is small; (3) the mass eigenstates are rotations of the weak interaction eigenstates. These rotations are given in a simple two-neutrino model as follows:(1)νe=cosθν1+sinθν2,νμ=-sinθν1+cosθν2,
where νi(i=1,2) is the “mass eigenstate,” να(α=e,μ) is the “flavor eigenstate,” and θ is the “mixing angle.” Under these conditions, a neutrino born in a pure flavor state through a weak decay can oscillate into another flavor as the state propagates in space, due to the fact that the different mass eigenstate components propagate with different frequencies. The mass splitting between the two states is Δm2=|m22-m12|>0. The oscillation probability for νμ→νe oscillations is then given by the following:
(2)P(νμ⟶νe)=sin22θsin2(1.27Δm2(eV2)L(km)E(GeV)),
where L is the distance from the source, and E is the neutrino energy.From (2), one can see that the probability for observing oscillations is large when Δm2~E/L. In the discussions below, we will focus on experiments with signals in the Δm2~1eV2 range. These experiments are therefore designed with E/L~1 GeV/km (or, alternatively, 1 MeV/m). Typically, neutrino source energies range from a few MeV to a few GeV. Thus, most of the experiments considered are located between a few meters and a few kilometers from the source. This is not absolutely necessary, a very high-energy experiment with a very long baseline is sensitive to oscillations in the Δm2~1eV2 range, as long as the ratio E/L~1 GeV/km is maintained. In other words, “short-baseline experiments” is something of a misnomer—what is meant is the experiments with sensitivity to Δm2~1eV2 oscillations.In the case whereE/L≪1GeV/km, such as in accelerator-based experiments with long baselines (hundreds of kilometers), one can see from (2) that the oscillations will be rapid. In the case of Δm2~1eV2, sensitivity to the mass splitting is lost because the sin2(1.27Δm2(L/E)) term will average to 1/2 due to the finite energy and position resolution of the experiment. The oscillation probability becomes P=(sin22θ)/2 in this case. Thus, the information from “long-baseline experiments” can be used to constrain the mixing angle, but not the Δm2.The exercise of generalizing to a three-neutrino model is useful, since the inclusion of more states follows from this procedure. Within a three-neutrino model, the mixing matrix is written as follows:(3)(νeνμντ)=(Ue1Ue2Ue3Uμ1Uμ2Uμ3Uτ1Uτ2Uτ3)(ν1ν2ν3).
The matrix elements are parametrized by three mixing angles, analogous to the Euler angles. As in the quark sector, the three-neutrino model can be extended to include an imaginary term that introduces a CP-violating phase. This formalism is analogous to the quark sector, where strong and weak eigenstates are rotated and the resultant mixing is described conventionally by a unitary mixing matrix.The oscillation probability for three-neutrino oscillations is typically written as the following:(4)P(να⟶νβ)=δαβ-4∑j>i‍UαiUβi*Uαj*Uβjsin2(1.27Δmij2LE),
where Δmij2=mj2-mi2, α and β are flavor-state indices (e,μ,τ), and i and j are mass-state indices (1,2,3 in the three-neutrino case, though (4) holds for n-neutrino oscillations). Although in general there will be mixing among all three flavors of neutrinos, if the mass scales are quite different (m3≫m2≫m1), then the oscillation phenomena tend to decouple and the two-neutrino mixing model is a good approximation in limited regions.Three differentΔm2 parameters appear in (4); however, only two are independent since the two small Δm2 parameters must sum to the largest. If we consider the oscillation data measured at >5σ [1–16], then two Δm2 ranges, 7×10-5eV2 (solar) and 3×10-3 eV2 (atmospheric), are already defined. These constrain the third Δm2, so that oscillation results at ~1 eV2, such as those discussed in this paper, cannot be accommodated within a three-neutrino model.


### 2.3. (3+N) Oscillation Formalism

The sterile neutrino oscillation formalism followed in this paper assumes up to three additional neutrino mass eigenstates, beyond the established three SM neutrino species. We know, from solar and atmospheric oscillation observations, that three of the mass states must be mostly active. Experimental hints point toward the existence of additional mass states that are mostly sterile, in the range ofΔm2 = 0.01–100 eV2.Introducing extra mass states results in a large number of extra parameters in the model. Approximation is required to allow for efficient exploration of the available parameters. To this end, in our model we assume that the three lowest states,ν1, ν2, and ν3, that are the mostly active states accounting for the solar and atmospheric observations, have masses so small as to be effectively degenerate with equal masses. This is commonly called the “short-baseline approximation” and it reduces the picture to two-, three-, and four-neutrino-mass oscillation models, corresponding to (3+1), (3+2), and (3+3), respectively.The active(e,μ,τ) content of the N additional mass eigenstates is assumed to be small; specifically, the Uαi elements of the extended (3+N)×(3+N) mixing matrix for i=4–6 and α=e,μ,τ, are restricted to values |Uαi|≤0.5, while the following constraints are applied by way of unitarity:
(5)∑α=e,μ,τ‍|Uαi|2≤0.3,
for each i=4–6, and
(6)∑i=4-6‍|Uαi|2≤0.3,
for each α=e,μ,τ. In our fits, since the SBL experiments considered have no ντ sensitivity, we explicitly assume that |Uτi|=0. The above restrictions therefore apply only for α=e,μ, and are consistent with solar and atmospheric neutrino experiments, which indicate that there can only be a small electron and muon flavor content in the fourth, fifth, and sixth mass eigenstates [23].In this formalism, the probabilities forνα→νβ oscillations can be deduced from the following equation:
(7)P(να⟶νβ)=δαβ-∑j<i‍(4Re{UβiUαi*Uβj*Uαj}sin2(1.27Δmij2LE)-2Im{UβiUαi*Uβj*Uαj}sin(2.53Δmij2LE)),
where Δmij2=mi2-mj2 is in eV2, L is in m, and  E is in MeV. This formalism conserves CPT, but does not necessarily conserve CP.To be explicit, for the (3+3) scenario, the mixing formalism is extended in the following way:
(8)(νeνμντνs1νs2νs3)=(Ue1Ue2Ue3Ue4Ue5Ue6Uμ1Uμ2Uμ3Uμ4Uμ5Uμ6Uτ1Uτ2Uτ3Uτ4Uτ5Uτ6Us11Us12Us13Us14Us15Us16Us21Us22Us23Us24Us25Us36Us31Us32Us33Us34Us35Us36)(ν1ν2ν3ν4ν5ν6).
The SBL approximation states that m1≈m2≈m3. With this assumption, and for the case of the (3+3) scenario, the appearance (α≠β) oscillation probability can be rewritten as the following:
(9)P(να⟶νβ)≃-4|Uα5||Uβ5||Uα4||Uβ4|cosϕ54sin2(1.27Δm542LE)-4|Uα6||Uβ6||Uα4||Uβ4|cosϕ64sin2(1.27Δm642LE)-4|Uα5||Uβ5||Uα6||Uβ6|cosϕ65sin2(1.27Δm652LE)+4(|Uα4||Uβ4|+|Uα5||Uβ5|cosϕ54+|Uα6||Uβ6|cosϕ64)×|Uα4||Uβ4|sin2(1.27Δm412LE)+4(|Uα4||Uβ4|cosϕ54+|Uα5||Uβ5|+|Uα6||Uβ6|cosϕ65)×|Uα5||Uβ5|sin2(1.27Δm512LE)+4(|Uα4||Uβ4|cosϕ64+|Uα5||Uβ5|cosϕ65+|Uα6||Uβ6|)×|Uα6||Uβ6|sin2(1.27Δm612LE)+2|Uβ5||Uα5||Uβ4||Uα4|sinϕ54sin(2.53Δm542LE)+2|Uβ6||Uα6||Uβ4||Uα4|sinϕ64sin(2.53Δm642LE)+2|Uβ6||Uα6||Uβ5||Uα5|sinϕ65sin(2.53Δm652LE)+2(|Uα5||Uβ5|sinϕ54+|Uα6||Uβ6|sinϕ64)×|Uα4||Uβ4|sin(2.53Δm412LE)+2(-|Uα4||Uβ4|sinϕ54+|Uα6||Uβ6|sinϕ65)×|Uα5||Uβ5|sin(2.53Δm512LE)+2(-|Uα4||Uβ4|sinϕ64-|Uα5||Uβ5|sinϕ65)×|Uα6||Uβ6|sin(2.53Δm612LE).
CP violation appears in (9) in the form of the three phases defined by
(10)ϕ54=arg(Ue5Uμ5*Ue4*Uμ4),ϕ64=arg(Ue6Uμ6*Ue4*Uμ4),ϕ65=arg(Ue6Uμ6*Ue5*Uμ5).
In each case, ν→ν- implies ϕ→-ϕ. In the case of disappearance (α=β), the survival probability can be rewritten as the following:
(11)P(να⟶να)≃1-4|Uα4|2|Uα5|2sin2(1.27Δm542LE)-4|Uα4|2|Uα6|2sin2(1.27Δm642LE)-4|Uα5|2|Uα6|2sin2(1.27Δm652LE)-4(1-|Uα4|2-|Uα5|2-|Uα6|2)×(|Uα4|2sin2(1.27Δm412LE)+|Uα5|2sin2×(1.27Δm512LE)+|Uα6|2sin2(1.27Δm612LE)).
This formula has no ϕij dependencies because CP violation only affects appearance.We have discussed the formulas for (3+1) and (3+2) oscillations that arise from (7) in previous papers [20–22]. To reduce to a (3+2) model, the parameters Δm612, |Ue6|, |Uμ6|, ϕ64, and ϕ65 are explicitly set to zero; consequently, we have the following appearance and disappearance formulas for a (3+2) model:
(12)P(να⟶νβ)≃-4|Uα5||Uβ5||Uα4||Uβ4|cosϕ54sin2(1.27Δm542LE)+4(|Uα4||Uβ4|+|Uα5||Uβ5|cosϕ54)×|Uα4||Uβ4|sin2(1.27Δm412LE)+4(|Uα4||Uβ4|cosϕ54+|Uα5||Uβ5|)×|Uα5||Uβ5|sin2(1.27Δm512LE)+2|Uβ5||Uα5||Uβ4||Uα4|sinϕ54sin(2.53Δm542LE)+2(|Uα5||Uβ5|sinϕ54)|Uα4||Uβ4|sin(2.53Δm412LE)+2(-|Uα4||Uβ4|sinϕ54)|Uα5||Uβ5|sin(2.53Δm512LE),P(να⟶να)≃1-4|Uα4|2|Uα5|2sin2(1.27Δm542LE)-4(1-|Uα4|2-|Uα5|2)×(|Uα4|2sin2(1.27Δm412LE)+|Uα5|2sin2(1.27Δm512LE)).For a (3+1) model, Δm612, Δm512, |Ue6|, |Uμ6|, |Ue5|, |Uμ5|, ϕ64, ϕ65, and ϕ54 should be set to zero. This further simplifies the oscillation probabilities, and one recovers the familiar two-neutrino appearance and disappearance probabilities. The appearance and disappearance formulas for a (3+1) model are then given by the following:
(13)P(να⟶νβ)≃4|Uα4|2|Uβ4|2sin2(1.27Δm412LE),P(να⟶να)≃1-4(1-|Uα4|2)|Uα4|2sin2(1.27Δm412LE).In principle, the probability for neutrino oscillation is modified in the presence of matter. “Matter effects” arise because the electron neutrino flavor experiences both Charged-Current (CC) and Neutral-Current (NC) elastic forward scattering with electrons as it propagates through matter, while theνμ and ντ experience only NC forward-scattering. The sterile component experiences no forward-scattering. In practice, SM-inspired matter effects are very small given the short baselines of the experiments, and so we do not consider them further here. Beyond-SM matter effects are beyond the scope of this paper, but are considered in [38].


## 2.1. Light Sterile States

Sterile neutrinos are additional states beyond the standard electron, muon, and tau flavors, which do not interact via the exchange ofW or Z bosons [25] and are thus “sterile” with respect to the weak interaction. These states are motivated by many Beyond Standard Model theories, where they are often introduced as gauge singlets. Traditionally, sterile neutrinos were introduced at very high mass scales within the context of grand unification and leptogenesis. For many years, sterile neutrinos with light masses were regarded as less natural. However, as recent data [17, 19, 26, 27] has indicated the potential existence of light sterile neutrinos, the theoretical view has evolved to accommodate these light mass gauge singlets [28, 29]. At this point, it is generally accepted that the mass scale for sterile neutrinos is not well predicted, and the existence of one or more sterile neutrinos accommodated by introducing extra neutrino mass states at the eV scale is possible. An excellent review of the phenomenology of sterile neutrinos, as well as the data motivating light sterile models, is provided in [23].Within the expanded oscillation phenomenology, sterile neutrinos are handled as additional noninteracting flavors, which are connected to additional mass states via an extended mixing matrix with extra mixing angles and CP violating phases. These additional mass states must be mostly sterile, with only a small admixture of the active flavors, in order to accommodate the limits on oscillations to sterile neutrinos from the atmospheric and solar neutrino data. Experimental evidence for these additional mass states would come from the disappearance of an active flavor to a sterile neutrino state or additional transitions from one active flavor to another through the sterile neutrino state.The number of light sterile neutrinos is not predicted by theory. However, a natural tendency is to introduce three sterile states. Depending on how the states are distributed in mass scale, one, two, or all three states may be involved in SBL oscillations. These are referred to as (3+N) models where the “3” refers to the three active flavors and the “N" refers to the number of sterile neutrinos.Introducing sterile neutrinos can have implications in cosmological observations, especially measurements of the radiation density in the early universe. These are compounded if the extra neutrinos have significant mass (>1 eV) and do not decay. Currently, cosmological data allow additional states and in many cases favor light sterile neutrinos [30–36]. Upcoming Planck data [37] is expected to precisely measure Neff. This parameter, however, can be considered a model-dependent one. As an example, there are a variety of classes of theories where the neutrinos do not thermalize in the early universe [23]. In these cases, the cosmological neutrino abundance would substantially decrease, rendering cosmological measurements of Neff invalid. Therefore, while the community certainly looks forward to cosmological measurements of Neff, we think that SBL experiments are a largely better approach for probing light sterile neutrinos and constraining their mixing properties. We therefore proceed with a study of the SBL data, without further reference to the cosmological results.


## 2.2. The Basic Oscillation Formalism

Before considering the phenomenology of light sterile neutrinos, it is useful to introduce the idea of oscillations within a simpler model. In this section, we first consider the two-neutrino formalism. We then extend these ideas to form the well-established three-active-flavor neutrino model. Based on these concepts, we expand the discussion to include more states in the following section.Neutrino oscillations require that (1) neutrinos have mass; (2) the difference between the masses is small; (3) the mass eigenstates are rotations of the weak interaction eigenstates. These rotations are given in a simple two-neutrino model as follows:(1)νe=cosθν1+sinθν2,νμ=-sinθν1+cosθν2,
where νi(i=1,2) is the “mass eigenstate,” να(α=e,μ) is the “flavor eigenstate,” and θ is the “mixing angle.” Under these conditions, a neutrino born in a pure flavor state through a weak decay can oscillate into another flavor as the state propagates in space, due to the fact that the different mass eigenstate components propagate with different frequencies. The mass splitting between the two states is Δm2=|m22-m12|>0. The oscillation probability for νμ→νe oscillations is then given by the following:
(2)P(νμ⟶νe)=sin22θsin2(1.27Δm2(eV2)L(km)E(GeV)),
where L is the distance from the source, and E is the neutrino energy.From (2), one can see that the probability for observing oscillations is large when Δm2~E/L. In the discussions below, we will focus on experiments with signals in the Δm2~1eV2 range. These experiments are therefore designed with E/L~1 GeV/km (or, alternatively, 1 MeV/m). Typically, neutrino source energies range from a few MeV to a few GeV. Thus, most of the experiments considered are located between a few meters and a few kilometers from the source. This is not absolutely necessary, a very high-energy experiment with a very long baseline is sensitive to oscillations in the Δm2~1eV2 range, as long as the ratio E/L~1 GeV/km is maintained. In other words, “short-baseline experiments” is something of a misnomer—what is meant is the experiments with sensitivity to Δm2~1eV2 oscillations.In the case whereE/L≪1GeV/km, such as in accelerator-based experiments with long baselines (hundreds of kilometers), one can see from (2) that the oscillations will be rapid. In the case of Δm2~1eV2, sensitivity to the mass splitting is lost because the sin2(1.27Δm2(L/E)) term will average to 1/2 due to the finite energy and position resolution of the experiment. The oscillation probability becomes P=(sin22θ)/2 in this case. Thus, the information from “long-baseline experiments” can be used to constrain the mixing angle, but not the Δm2.The exercise of generalizing to a three-neutrino model is useful, since the inclusion of more states follows from this procedure. Within a three-neutrino model, the mixing matrix is written as follows:(3)(νeνμντ)=(Ue1Ue2Ue3Uμ1Uμ2Uμ3Uτ1Uτ2Uτ3)(ν1ν2ν3).
The matrix elements are parametrized by three mixing angles, analogous to the Euler angles. As in the quark sector, the three-neutrino model can be extended to include an imaginary term that introduces a CP-violating phase. This formalism is analogous to the quark sector, where strong and weak eigenstates are rotated and the resultant mixing is described conventionally by a unitary mixing matrix.The oscillation probability for three-neutrino oscillations is typically written as the following:(4)P(να⟶νβ)=δαβ-4∑j>i‍UαiUβi*Uαj*Uβjsin2(1.27Δmij2LE),
where Δmij2=mj2-mi2, α and β are flavor-state indices (e,μ,τ), and i and j are mass-state indices (1,2,3 in the three-neutrino case, though (4) holds for n-neutrino oscillations). Although in general there will be mixing among all three flavors of neutrinos, if the mass scales are quite different (m3≫m2≫m1), then the oscillation phenomena tend to decouple and the two-neutrino mixing model is a good approximation in limited regions.Three differentΔm2 parameters appear in (4); however, only two are independent since the two small Δm2 parameters must sum to the largest. If we consider the oscillation data measured at >5σ [1–16], then two Δm2 ranges, 7×10-5eV2 (solar) and 3×10-3 eV2 (atmospheric), are already defined. These constrain the third Δm2, so that oscillation results at ~1 eV2, such as those discussed in this paper, cannot be accommodated within a three-neutrino model.


## 2.3. (3+N) Oscillation Formalism

The sterile neutrino oscillation formalism followed in this paper assumes up to three additional neutrino mass eigenstates, beyond the established three SM neutrino species. We know, from solar and atmospheric oscillation observations, that three of the mass states must be mostly active. Experimental hints point toward the existence of additional mass states that are mostly sterile, in the range ofΔm2 = 0.01–100 eV2.Introducing extra mass states results in a large number of extra parameters in the model. Approximation is required to allow for efficient exploration of the available parameters. To this end, in our model we assume that the three lowest states,ν1, ν2, and ν3, that are the mostly active states accounting for the solar and atmospheric observations, have masses so small as to be effectively degenerate with equal masses. This is commonly called the “short-baseline approximation” and it reduces the picture to two-, three-, and four-neutrino-mass oscillation models, corresponding to (3+1), (3+2), and (3+3), respectively.The active(e,μ,τ) content of the N additional mass eigenstates is assumed to be small; specifically, the Uαi elements of the extended (3+N)×(3+N) mixing matrix for i=4–6 and α=e,μ,τ, are restricted to values |Uαi|≤0.5, while the following constraints are applied by way of unitarity:
(5)∑α=e,μ,τ‍|Uαi|2≤0.3,
for each i=4–6, and
(6)∑i=4-6‍|Uαi|2≤0.3,
for each α=e,μ,τ. In our fits, since the SBL experiments considered have no ντ sensitivity, we explicitly assume that |Uτi|=0. The above restrictions therefore apply only for α=e,μ, and are consistent with solar and atmospheric neutrino experiments, which indicate that there can only be a small electron and muon flavor content in the fourth, fifth, and sixth mass eigenstates [23].In this formalism, the probabilities forνα→νβ oscillations can be deduced from the following equation:
(7)P(να⟶νβ)=δαβ-∑j<i‍(4Re{UβiUαi*Uβj*Uαj}sin2(1.27Δmij2LE)-2Im{UβiUαi*Uβj*Uαj}sin(2.53Δmij2LE)),
where Δmij2=mi2-mj2 is in eV2, L is in m, and  E is in MeV. This formalism conserves CPT, but does not necessarily conserve CP.To be explicit, for the (3+3) scenario, the mixing formalism is extended in the following way:
(8)(νeνμντνs1νs2νs3)=(Ue1Ue2Ue3Ue4Ue5Ue6Uμ1Uμ2Uμ3Uμ4Uμ5Uμ6Uτ1Uτ2Uτ3Uτ4Uτ5Uτ6Us11Us12Us13Us14Us15Us16Us21Us22Us23Us24Us25Us36Us31Us32Us33Us34Us35Us36)(ν1ν2ν3ν4ν5ν6).
The SBL approximation states that m1≈m2≈m3. With this assumption, and for the case of the (3+3) scenario, the appearance (α≠β) oscillation probability can be rewritten as the following:
(9)P(να⟶νβ)≃-4|Uα5||Uβ5||Uα4||Uβ4|cosϕ54sin2(1.27Δm542LE)-4|Uα6||Uβ6||Uα4||Uβ4|cosϕ64sin2(1.27Δm642LE)-4|Uα5||Uβ5||Uα6||Uβ6|cosϕ65sin2(1.27Δm652LE)+4(|Uα4||Uβ4|+|Uα5||Uβ5|cosϕ54+|Uα6||Uβ6|cosϕ64)×|Uα4||Uβ4|sin2(1.27Δm412LE)+4(|Uα4||Uβ4|cosϕ54+|Uα5||Uβ5|+|Uα6||Uβ6|cosϕ65)×|Uα5||Uβ5|sin2(1.27Δm512LE)+4(|Uα4||Uβ4|cosϕ64+|Uα5||Uβ5|cosϕ65+|Uα6||Uβ6|)×|Uα6||Uβ6|sin2(1.27Δm612LE)+2|Uβ5||Uα5||Uβ4||Uα4|sinϕ54sin(2.53Δm542LE)+2|Uβ6||Uα6||Uβ4||Uα4|sinϕ64sin(2.53Δm642LE)+2|Uβ6||Uα6||Uβ5||Uα5|sinϕ65sin(2.53Δm652LE)+2(|Uα5||Uβ5|sinϕ54+|Uα6||Uβ6|sinϕ64)×|Uα4||Uβ4|sin(2.53Δm412LE)+2(-|Uα4||Uβ4|sinϕ54+|Uα6||Uβ6|sinϕ65)×|Uα5||Uβ5|sin(2.53Δm512LE)+2(-|Uα4||Uβ4|sinϕ64-|Uα5||Uβ5|sinϕ65)×|Uα6||Uβ6|sin(2.53Δm612LE).
CP violation appears in (9) in the form of the three phases defined by
(10)ϕ54=arg(Ue5Uμ5*Ue4*Uμ4),ϕ64=arg(Ue6Uμ6*Ue4*Uμ4),ϕ65=arg(Ue6Uμ6*Ue5*Uμ5).
In each case, ν→ν- implies ϕ→-ϕ. In the case of disappearance (α=β), the survival probability can be rewritten as the following:
(11)P(να⟶να)≃1-4|Uα4|2|Uα5|2sin2(1.27Δm542LE)-4|Uα4|2|Uα6|2sin2(1.27Δm642LE)-4|Uα5|2|Uα6|2sin2(1.27Δm652LE)-4(1-|Uα4|2-|Uα5|2-|Uα6|2)×(|Uα4|2sin2(1.27Δm412LE)+|Uα5|2sin2×(1.27Δm512LE)+|Uα6|2sin2(1.27Δm612LE)).
This formula has no ϕij dependencies because CP violation only affects appearance.We have discussed the formulas for (3+1) and (3+2) oscillations that arise from (7) in previous papers [20–22]. To reduce to a (3+2) model, the parameters Δm612, |Ue6|, |Uμ6|, ϕ64, and ϕ65 are explicitly set to zero; consequently, we have the following appearance and disappearance formulas for a (3+2) model:
(12)P(να⟶νβ)≃-4|Uα5||Uβ5||Uα4||Uβ4|cosϕ54sin2(1.27Δm542LE)+4(|Uα4||Uβ4|+|Uα5||Uβ5|cosϕ54)×|Uα4||Uβ4|sin2(1.27Δm412LE)+4(|Uα4||Uβ4|cosϕ54+|Uα5||Uβ5|)×|Uα5||Uβ5|sin2(1.27Δm512LE)+2|Uβ5||Uα5||Uβ4||Uα4|sinϕ54sin(2.53Δm542LE)+2(|Uα5||Uβ5|sinϕ54)|Uα4||Uβ4|sin(2.53Δm412LE)+2(-|Uα4||Uβ4|sinϕ54)|Uα5||Uβ5|sin(2.53Δm512LE),P(να⟶να)≃1-4|Uα4|2|Uα5|2sin2(1.27Δm542LE)-4(1-|Uα4|2-|Uα5|2)×(|Uα4|2sin2(1.27Δm412LE)+|Uα5|2sin2(1.27Δm512LE)).For a (3+1) model, Δm612, Δm512, |Ue6|, |Uμ6|, |Ue5|, |Uμ5|, ϕ64, ϕ65, and ϕ54 should be set to zero. This further simplifies the oscillation probabilities, and one recovers the familiar two-neutrino appearance and disappearance probabilities. The appearance and disappearance formulas for a (3+1) model are then given by the following:
(13)P(να⟶νβ)≃4|Uα4|2|Uβ4|2sin2(1.27Δm412LE),P(να⟶να)≃1-4(1-|Uα4|2)|Uα4|2sin2(1.27Δm412LE).In principle, the probability for neutrino oscillation is modified in the presence of matter. “Matter effects” arise because the electron neutrino flavor experiences both Charged-Current (CC) and Neutral-Current (NC) elastic forward scattering with electrons as it propagates through matter, while theνμ and ντ experience only NC forward-scattering. The sterile component experiences no forward-scattering. In practice, SM-inspired matter effects are very small given the short baselines of the experiments, and so we do not consider them further here. Beyond-SM matter effects are beyond the scope of this paper, but are considered in [38].


## 3. Experimental Datasets

This section provides an overview of the various types of past and current neutrino sources and detectors used in SBL experiments. After introducing the experimental concepts, the specific experimental datasets used in this analysis are discussed.The data fall into two overall categories: disappearance, where the active flavor is assumed to have oscillated into a sterile neutrino and/or another flavor which is kinematically not allowed to interact or leaves no detectable signature, and appearance, where the transition is between active flavors, but with mass splittings corresponding to the mostly sterile states. Appearance and disappearance are natural divisions for testing the compatibility of datasets. If|Uα4|2 and |Uβ4|2 are shown to be small, then the effective mixing angle for appearance, 4|Uα4|2|Uβ4|2, cannot be large. This constraint that the disappearance experiments place on appearance experiments extends to (3+2) and (3+3) models also.CPT conservation, which is assumed in the analysis, demands that neutrino and antineutrino disappearance probabilities are the same after accounting for cp violating effects. To test this, we divide the data into antineutrino and neutrino sets and fit each set separately. If CP violation is already allowed in the oscillation formalism, then any incompatibility found between respective neutrino and antineutrino fits could imply effective CPT violation, as discussed in [22].Figures1, 2, and 3 provide summaries of the datasets, showing the constraints they provide in a simple two-neutrino oscillation model, which is functionally equivalent to the (3+1) scenario. Figure 1 shows the muon-to-electron flavor datasets in neutrino and antineutrino mode at 95% confidence level (CL). Figures 2 and 3 show results for νμ and ν-μ, and νe and ν-e disappearance, respectively.Figure 1
Summary ofν-μ→ν-e and νμ→νe results, shown at 95% CL. Top row: LSND, KARMEN, BNB-MB(νapp); bottom row: BNB-MB(ν-app), NuMI-MB(νapp), NOMAD. See Section 3.2 for details and references.Figure 2
Summary ofν-μ→ν-μ and νμ→νμ results, shown at 95% CL. Top row: BNB-MB(νdis), CCFR84, CDHS; bottom row: MINOS-CC, ATM. See Section 3.2 for details and references.Figure 3
Summary ofν-e→ν-e and νe→νe results, shown at 95% CL. From left: KARMEN/LSND(xsec), Bugey, and Gallium. See Section 3.2 for details and references.


### 3.1. Sources and Detectors Used in Short-Baseline Neutrino Experiments

Before considering the datasets in detail, we provide an overview of how SBL experiments are typically designed.


#### 3.1.1. Sources of Neutrinos for Short-Baseline Experiments

The neutrino sources used in SBL experiments range in energy from a few MeV to hundreds of GeV and include manmade radioactive sources, reactors, and accelerator-produced beams. While the higher energy accelerator sources are mixtures of different neutrino flavors, the <10 MeV sources rely on beta decay and are thus pure electron neutrino flavor.At the low-energy end of the spectrum, the rate of electron neutrino interactions from the beta decay of the ~1 MCi sources51Cr (half-life: 28 days) and 37Ar (half-life: 35 days) have been studied. These sources were originally produced for the low-energy (~1 MeV) calibration of solar neutrino detectors [39, 40] but have proven themselves interesting as a probe of electron neutrino disappearance.Moving up in energy by a few MeV, nuclear reactors are powerful sources of ~2−8 MeVν-e through the β+-decaying elements produced primarily in the decay chains of 235U, 239Pu, 238U, and 241Pu. While these four isotopes are the progenitors of most of the reactor flux, modern reactor simulations include all fission sources [41]. Reactor simulations convolute predictions of fission rates over time with neutrino production per fission. Recently, a reanalysis of the production cross-section per fission [23, 42, 43] has led to an increase in the predicted reactor flux. As their energy is too low for an appearance search (the neutrino energy is below the muon production kinematic threshold), reactor source antineutrinos  can only be used for ν-e disappearance searches, where the antineutrinos are detected using CC interactions with an outgoing e+.The lowest neutrino energy (up to 53 MeV) accelerator sources used in existing SBL experiments are based on pion- and muon-decay-at-rest (DAR). The neutrino flux comes from the stopped pion decay chain:π+→μ+νμ and μ+→e+ν-μνe. Pions are produced in interactions of accelerator protons with, typically, a graphite or water target. The contribution from the decay chain π-→μ-ν-μ is suppressed by designing the target such that the π- mesons are captured with high probability. The result is a source which has a well-understood neutrino flavor content and energy distribution, with a minimal (<10−3) ν-e content [44, 45]. This last point is important as  ν-μ→ν-e is the dominant channel used for oscillation searches by DAR sources.In a conventional high-energy (from ~100 MeV to hundreds of GeV) accelerator-based neutrino beam, protons impinge on a target (beryllium and carbon are typical) to produce secondary mesons. The boosted mesons enter and subsequently decay inside a long, often evacuated, pipe. Neutrinos are primarily produced byπ+ and π- decay in flight (DIF). Pion sign selection, via a large magnet placed directly in the beamlines before the decay pipe, allows for nearly pure neutrino or antineutrino running, with only a few percent “wrong sign” neutrino flux content in the case of neutrino running, and ~15% [46] in the case of antineutrino running. These beams are generally produced by protons at 8 GeV and above. At these energies, in addition to pion production, kaon production contributes to the flux of both muon and electron neutrino flavors. There is often a substantial muon DIF content as well, contributing both νe/ν-e and ν-μ/νμ to the beam. The result of the kaon and muon secondary content is that, while the neutrinos are predominantly muon flavored, the beam will always have some intrinsic electron flavor neutrino content, usually at the several percent level. Accelerator-based beams are predominantly used for νμ→νe and ν-μ→ν-e appearance searches, as well as νμ and ν-μ disappearance searches. An excellent review of methods in producing accelerator-based neutrino beams can be found in [47].In contrast to lower-energy neutrino sources (DAR, reactor, and isotope sources), high-energy accelerator-based neutrino sources are subject to significant energy-dependent neutrino flux uncertainties, often at the level of 10–15%, due to in-target meson production uncertainties. These uncertainties can affect the energy distribution, flavor content, and absolute normalization of a neutrino beam. Typically, meson production systematics are constrained with dedicated measurements by experiments such as HARP [48] and MIPP [49], which use replicated targets (geometry and material) and a wide range of proton beam energies to study meson production cross-sections and kinematics directly. Alternatively, experiments can employ a two-detector design for comparing near-to-far event rate in energy to effectively reduce these systematics. However, due to the short baselines employed for studying sterile neutrino oscillations, a two-detector search is often impractical. In situ measurements in single-detector experiments can exploit flux (multiplied by cross-section) correlations among different beam components and energies to reduce flux uncertainties, as has been done in the case of the MiniBooNE νe and ν-e appearance searches described below.


#### 3.1.2. Short-Baseline Neutrino Detectors

Because low-energy neutrino interaction cross-sections are very small, the options for SBL detectors are typically limited to designs which can be constructed on a massive scale. There are several generic neutrino detection methods in use today: unsegmented scintillator detectors, unsegmented Cerenkov detectors, segmented scintillator-and-iron calorimeters, and segmented trackers.Neutrino oscillation experiments usually require sensitivity to CC neutrino interactions, whereby one can definitively identify the flavor of the interacting neutrino by the presence of a charged lepton in the final state. However, in the case of sterile neutrino oscillation searches, NC interactions can also provide useful information, as they are directly sensitive to the sterile flavor content of the neutrino mass eigenstate,|Usi|2=1-|Uei|2-|Uμi|2-|Uτi|2.Unsegmented scintillator detectors are typically used for few-MeV-scale SBL experiments, which require efficient electron neutrino identification and reconstruction. These detectors consist of large tanks of oil-based (CnH2n) liquid scintillator surrounded by phototubes. The free protons in the oil provide a target for the inverse beta decay interaction, ν-ep→e+n. The reaction threshold for this interaction is 1.8 MeV due to the mass difference between the proton and neutron and the mass of the positron. The scintillation light from the e+, as well as light from the Compton scattering of the 0.511 MeV annihilation photons provides an initial (“prompt”) signal. This is followed by n capture on hydrogen and a 2.2 MeV flash of light, as the resulting γ Compton-scatters in the scintillator. This coincidence sequence in time (positron followed by neutron capture) provides a clean, mostly background-free interaction signature. Experiments often dope the liquid scintillator using an element with a high neutron capture cross-section for improved event identification efficiency.The CC interaction with the carbon in the oil (which produces either nitrogen or boron depending on whether the scatterer is a neutrino or antineutrino) has a significantly higher energy threshold than the free proton target-scattering process. The CC quasielastic interactionνe+C→e-+N has an energy threshold of 17.3 MeV, which arises from the carbon-nitrogen mass difference and the mass of the electron. In the case of both reactor and radioactive decay sources, the flux cuts off below this energy threshold. However, neutrinos from DAR sources are at sufficiently high energy to produce these carbon scatters.Unsegmented Cerenkov detectors make use of a target which is a large volume of clear medium (undoped oil or water is typical) surrounded by, or interspersed with, phototubes. Undoped oil has a larger refractive index, leading to a larger Cerenkov opening angle. Water is the only affordable medium once the detector size surpasses a few kilotons. In this paper, the only unsegmented Cerenkov detector that is considered is the 450-ton oil-based MiniBooNE detector. In such a detector, a track will project a ring with a sharp inner and outer edge onto the phototubes. Consider an electron produced in aνe CC quasielastic interaction. As the electron is low mass, it will multiple-scatter and easily bremsstrahlung, smearing the light projected on the tubes and producing a “fuzzy” ring. A muon produced by a CC quasielastic νμ interaction (νμn→μ-p) is heavier and will thus produce a sharper outer edge to the ring. For the same visible energy, the track will also extend farther, filling the interior of the ring and, perhaps, exiting the tank. If the muon stops within the tank and subsequently decays, the resulting electron provides an added tag for particle identification. In the case of the μ-, 18% will capture in water and, thus, have no electron tag, while only 8% will capture in the oil.Scintillator and iron calorimeters provide an affordable detection technique for ~1 GeV and higherνμ interactions. At these energies, multiple hadrons may be produced at the interaction vertex and will be observed as hadronic showers. In these devices, the iron provides the target, while the scintillator provides information on energy deposition per unit length. This information allows separation between the hadronic shower, which occurs in both NC and CC events, and the minimum-ionizing track of an outgoing muon, which occurs in CC events. Transverse information can be obtained if segmented scintillator strips are used, or if drift chambers are interspersed. The light from scintillator strips is transported to tubes by wavelength-shifting fibers. Information in the transverse plane improves separation of electromagnetic and hadronic showers. The iron can be magnetized to allow separation of neutrino and antineutrino events based on the charge of the outgoing lepton.To address the problem of running at ~1 GeV, where hadron track reconstruction is desirable, highly segmented tracking designs have been developed. The best resolution comes from stacks of wire chambers, where the material enclosing the gas provides the target. However, a more practical alternative has been stacks of thin extruded scintillator bars that are read out using wavelength-shifting fibers.


### 3.2. Data Used in the Sterile Neutrino Fits

There are many SBL datasets that can be included in this analysis. In this work, we have substantially expanded the number of datasets used beyond those in our past papers [20–22]. In the sections following, we identify and discuss new and updated datasets, as well as provide information on those used in past fits. The fit technique is described in Section 4.2.


#### 3.2.1. Experimental Results from Decay at Rest Studies

In past sterile neutrino studies [20–22], we have included the LSND and KARMEN appearance results described below. Since that work, a new study that constrains νe disappearance from the relative LSND-to-KARMEN cross-section measurements was published [50]. This new dataset is included in this analysis.LSND Appearance. LSND was a DAR experiment that ran in the 1990s, searching for ν-μ→ν-e. The beam was produced using 800 MeV protons on target from the LAMPF accelerator at Los Alamos National Laboratory, where a 1 mA beam of protons impinged on a water target. The center of the 8.75 m long, nearly cylindrical detector, was located at 29.8 m from the target, at an angle of 12° from the proton beam direction. This was an unsegmented detector with a fiducial mass of 167 tons of oil (CH2), lightly doped with b-PBD scintillator. The intrinsic ν-e content of the beam was 8×10-4 of the ν-μ content. The experiment observed a ν-e excess of 87.9±22.4±6.0 events above background, which was interpreted as oscillations with a probability of (0.264±0.067±0.045)%. Details are available in [17].This dataset is referred to as LSND in the analysis below and indicates a signal at 95% CL, as shown in Figure1. This data covers energies between 20 and 53 MeV and contributes five energy bins to the global fit. Statistical errors are taken into account by using a log-likelihood χ2 definition in the fit, while systematic errors on the background prediction are not included because these are small relative to the statistical error. Energy and baseline smearing are taken into account by averaging the oscillation probability over the energy bin width and over the neutrino flight path uncertainty.KARMEN Appearance. KARMEN was another DAR experiment searching for ν-μ→ν-e. KARMEN ran at the ISIS facility at Rutherford Laboratory, with 200 μA of protons impinging on a copper, tantalum, or uranium target. The neutrino detector was located at an angle of 100° with respect to the targeting protons to reduce background from π- DIF. The resulting intrinsic ν-e content was 6.4×10-4 of the ν-μ content.The center of the approximately cubic segmented scintillator detector was located at 17.7 m. Thus, this detector was 60% of the distance from the source compared to LSND. The liquid scintillator target volume was 56 m3 and consisted of 512 optically independent modules (17.4 cm × 17.8 cm × 353 cm) wrapped in gadolinium-doped paper. KARMEN saw no signal and set a limit on appearance. More details are available in [51].This dataset is referred to as KARMEN in the analysis below and indicates a limit at 95% CL, as shown in Figure1. This dataset contributes nine energy bins, in the range 16 to 50 MeV. As in the case of LSND, statistical errors are taken into account by using a log-likelihood χ2 definition in the fit, while systematic errors on the background prediction are not included. Energy and baseline smearing are taken into account by averaging the sin2(1.27Δm2L/E) and sin(2.53Δm2L/E) term contributions in the total signal prediction over energy bin widths. The limit which is shown here is determined using a Δχ2-based raster scan, as discussed in Section 4.2.LSND and KARMEN Cross-Section Measurements. Along with the oscillation searches, LSND and KARMEN measured νe+12C→12Ngs+e- scattering. In this two-body interaction, with a Q-value of 17.3 MeV, the neutrino energy can be reconstructed by measuring the outgoing visible energy of the electron. The 12N ground state is identified by the subsequent β decay, 12Ngs→12C+e++νe, which has a Q-value of 16.3 MeV and a lifetime of 15.9 ms.The cross-section is measured by both experiments under the assumption that theνe flux has not oscillated, leading to disappearance. The excellent agreement between the two results, as a function of energy, allows a limit to be placed on νe oscillations. The energy dependence of the cross-section, as well as the normalization, are well predicted and both constraints are used in the analysis [50].This dataset is referred to as KARMEN/LSND(xsec) in the analysis below, and indicates a limit at 95% CL, as shown in Figure 3. A total of six (for KARMEN) plus five (for LSND) bins are used in the fit, which extend approximately from 28–50 MeV in the case of KARMEN and from 38–50 MeV in the case of LSND. In calculating the oscillation probability, the signal is averaged across the lengths of the detectors. The experiments have correlated systematics arising from the flux normalization due to a shared underlying analysis for pion production in DAR experiments. This is addressed through application of pull terms as described in [50].


#### 3.2.2. The MiniBooNE Experimental Results

The MiniBooNE experiment provides multiple results from a single detector. This oil-based 450 t fiducial volume Cerenkov detector was exposed to two conventional beams, the Booster Neutrino Beam (BNB) and the off-axis NuMI beam. The primary goal of MiniBooNE was to search forνμ→νe and ν-μ→ν-e appearance, using the BNB, which provides sensitivity to Δm2~0.1–10eV2 oscillations. The NuMI beam also provides some sensitivity to νμ→νe appearance at a similar Δm2. In addition to the appearance searches, MiniBooNE also looked for νμ and ν-μ disappearance using the BNB.The MiniBooNE datasets included in our analysis have increased throughout the period that our group has been performing fits. Reference [21] used a Monte Carlo prediction for neutrinos and antineutrinos to estimate MiniBooNE’s sensitivity to sterile neutrinos. The full BNB neutrino and first published BNB antineutrino datasets from MiniBooNE form the experimental constraints in [22]. Here, we have updated the analysis to include the full BNB antineutrino datasets. A further update has been to employ a log-likelihood method for the BNB neutrino and antineutrino datasets from [18], as this was recently adopted by the MiniBooNE Collaboration [52]. We also use the updated constraints on electron neutrino flux from kaons [18]. A partial dataset from NuMI data taking was presented in [22] and has not been updated, as the result was already systematics limited. In this analysis we also introduce the MiniBooNE disappearance search [53].In our fits to MiniBooNE appearance data, when drawing allowed regions and calculating compatibilities, which make use ofΔχ2’s and not absolute χ2’s, we use MiniBooNE’s log-likelihood χ2 definition, summing over both νe and νμ bins, as described in [52]. For consistency, the absolute MiniBooNE BNB νe and ν-e appearance χ2 values quoted in our paper also correspond to the same definition, that is, fitting to both νe and νμ spectra; therefore, they differ from the ones published by MiniBooNE in [18], which are obtained by fitting only to a priori constrained νe distributions. Note that the two definitions yield consistent allowed regions and compatibility results.The Booster Neutrino Beam Appearance Search in Neutrino Running Mode. The BNB flux composition in neutrino mode consists of >90% νμ, 6%ν-μ, and 0.06% νe and ν-e combined [46]. In the MiniBooNE BNB search for νe appearance, the νe and ν-e signal was normalized to the νμ and ν-μ CC quasielastic events observed in the detector, which peaked at 700 MeV.The global fits presented here use the full statistics of the MiniBooNEνμ→νe dataset, representing 6.46×1020 protons on target. In this dataset, MiniBooNE has observed an excess of events at 200–1250 MeV, corresponding to 162.0±47.8 electron-like events [18]. The dataset is referred to as BNB-MB(νapp) in the analysis below.We include the BNB-MB(νapp) dataset in our fits in the form of the full νe CC reconstructed energy distribution, in 11 energy bins from 200 to 3000 MeV, fit simultaneously with the full νμ CC energy distribution, in eight energy bins up to 1900 MeV. We account for statistical and systematic uncertainties in each sample, as well as systematic correlations (from flux and cross-section) among the νe signal and background and νμ background distributions. The systematic correlations are provided in the form of a full 19-bin×19-bin fractional covariance matrix. By fitting the νe and νμ spectra simultaneously, we are able to exploit the high-statistics νμ CC sample as a constraint on background and signal event rates. This assumes no significant νμ disappearance; this simplification could lead to a <20% effect on appearance probability obtained in MiniBooNE only fits [18].The dataset results in a signal at 95% CL, as shown in Figure1. This has changed slightly from our past analysis [22] now that we are using updated constraints on intrinsic electron neutrinos from kaons and the log-likelihood method, but is in agreement with the equivalent analysis from the MiniBooNE Collaboration [18].The Booster Neutrino Beam Appearance Search in Antineutrino Running Mode. The BNB flux composition in antineutrino mode consists of 83% ν-μ, 0.6% νe and ν-e combined, and a significantly larger wrong-sign composition than in neutrino mode, of 16%  νμ. As in the BNB νe appearance search, the electron flavor signal was normalized to the muon flavor CC quasielastic events observed in the detector, which peaked at 500 MeV.The global fits presented here use the full statistics of the MiniBooNEν-μ dataset, representing 11.27×1020 protons on target. In this dataset, MiniBooNE has observed an excess of events at 200–1250MeV, corresponding to 78.4±28.5 electron-like events. The dataset is referred to as BNB-MB(ν-app) in the analysis below.As in neutrino mode, we fit the fullν-e CC energy distribution, in 11 energy bins from 200 to 3000 MeV, simultaneously with the full ν-μ CC energy distribution, in 8 energy bins up to 1900 MeV. The wrong-sign contamination in the beam (νμ) is assumed to not contribute to any oscillations; only ν-μ→ν-e oscillations are assumed for this dataset. We account for statistical and systematic uncertainties in each sample as well as systematic correlations among the ν-e and ν-μ distributions in the form of a full 19-bin × 19-bin fractional covariance matrix in each fit. For further information, see [18].The dataset results in a signal at 95% CL, as shown in Figure1.The NuMI Beam Appearance Search. The MiniBooNE detector is also exposed to the NuMI neutrino beam, produced from a 120 GeV proton beam impinging on a carbon target. This beam is nominally used for the MINOS long-baseline neutrino oscillation experiment. NuMI events arrive out of time with the BNB-produced events. This 200 MeV to 3 GeV neutrino energy source is dominated by kaon decays near the NuMI target, which is 110 mrad off-axis and located 745 m upstream of the MiniBooNE detector. The beam consists of 81% νμ, 13% ν-μ, 5% νe, and 1% ν-e. For more information on this data, see [54].This dataset is referred to as NuMI-MB(νapp) in the analysis below. As seen in Figure 1, the dataset provides a limit at 95% CL. In the fits presented here, this data is used to constrain electron flavor appearance in neutrino mode, with 10 bins used in the fit. Statistical and systematic errors for this dataset are added in quadrature.The Booster Neutrino Beam Disappearance Search. The MiniBooNE experiment also searched for νμ and ν-μ disappearance using the BNB. The neutrino (antineutrino) dataset corresponded to 5.6×1020 (3.4×1020) protons on target, which produced a beam covering the neutrino energy range up to 1.9 GeV. The MiniBooNE νμ disappearance result provides restrictions on sterile neutrino oscillations which are comparable to those provided by the CDHS experiment, discussed below. Therefore, we include that dataset in these fits. On the other hand, the ν-μ result was weaker due to the combination of fewer protons on target and a lower cross-section. The MINOS ν-μ CC constraint, described below, is stronger, and so we do not use the MiniBooNE ν-μ dataset.The fit to theνμ dataset uses 16 bins ranging up to 1900 MeV in reconstructed neutrino energy. A shape-only fit is performed, where the predicted spectrum given any set of oscillation parameters is renormalized so that the total number of predicted events, after oscillations, is equal to the total number of observed events. Then the normalized predicted spectrum is compared to the observed spectrum in the form of a χ2 which accounts for statistical and shape-only systematic uncertainties and bin-to-bin correlations in the form of a covariance matrix.This dataset is referred to as BNB-MB(νdis) in the analysis below. Figure 2 shows that this data sets a limit at 95% CL. It should be noted that the published MiniBooNE analysis used a Pearson χ2 method [53], and we are able to reproduce those results. However, to fold these results into our analysis, we reverted to the Δχ2 definition used consistently among all datasets included in the fits (see Section 4.2).


#### 3.2.3. Results from Multi-GeV Conventional Short-Baselineνμ Beams

The set of multi-GeV conventional SBLνμ experiments is the same as was used in previous fits. Our overview of these experiments is therefore very brief.NOMAD Appearance Search. The NOMAD experiment [55], which ran at CERN using protons from the 450 GeV SPS accelerator, employed a conventional neutrino beamline to create a wideband 2.5 to 40 GeV neutrino energy source. These neutrinos were created with a carbon-based, low-mass tracking detector located 600 m downstream of the target. This detector had fine spatial resolution and could search for muon-to-electron and muon-to-tau oscillations. No signal was observed in either channel. In this analysis, we use the νμ→νe constraint.This dataset is referred to as NOMAD in the analysis below. This dataset contributes 30 energy bins to the global fit. The statistical and systematic errors are added in quadrature. This experiment sets a limit at 95% CL, as seen in Figure1.CCFR Disappearance Search. The CCFR dataset was taken at Fermilab in 1984 [56] with a narrowband beamline, with meson energies set to 100, 140, 165, 200, and 250 GeV, yielding νμ and ν-μ beams that ranged from 40 to 230 GeV in energy. This was a two-detector disappearance search, with the near detector at 715 m and the far detector at 1116 m from the center of the 352 m long decay pipe. The calorimetric detectors were constructed of segmented iron with scintillator and spark chambers, and each had a downstream toroid to measure the muon momentum.This dataset is referred to as CCFR84 in the analysis below. The data were published as the double ratios of the observed-to-expected rates in a near-to-far ratio. For each secondary mean setting, the data are divided into three energy bins. The systematic uncertainty is assumed to be energy independent and fully correlated between the energy bins. Due to the high beam energies and short baselines, this experiment sets a limit at highΔm2 in the muon flavor disappearance search at 95% CL, as shown in Figure 2.CDHS Disappearance Search. The CDHS experiment [57] at CERN searched for νμ disappearance with a two-detector design of segmented calorimeters with iron and scintillator. The experiment used 19.2 GeV protons on a beryllium target to produce mesons that were subsequently focused into a 52 m decay channel. The detectors were located 130 m and 885 m downstream of the target.This dataset is referred to as CDHS in the analysis below. CDHS provides data and errors in 15 bins of muon energy, as seen in in Table 1 of [57]. We relate these bins to the neutrino energy distributions using the method described in [20]: the neutrino energy distribution for a given muon energy or range is determined via the NUANCE [58] neutrino event generator. The experiment has a limit at 95% CL and sets constraints that are comparable to the MiniBooNE νμ disappearance limit described above, but which extend to slightly lower Δm2. See Figure 2 for comparison.


#### 3.2.4. Reactor and Source Experiments

The reactor experiment dataset has been updated to reflect recent changes in the predicted neutrino fluxes, as discussed below. The source-based experimental datasets are both new to this paper, and have been published since our last set of fits [22].Bugey Dataset. This analysis uses energy-dependent data from the Bugey 3 reactor experiment [59]. The detector consisted of 6Li-doped liquid scintillator, with data taken at 15, 45, and 90 m from the 2.8 GW reactor source. The detectors are taken to be pointlike in the analysis.Recently, a reanalysis of reactorν-e flux predictions [23, 42, 43] has led to a reinterpretation of the Bugey data. The data has transitioned from being simply a limit on neutrino disappearance to an allowed region at 95% CL. In this analysis, we adjust the predicted Bugey flux spectra normalization according to the calculations from [23].There are many other SBL reactor datasets in existence. However, we have chosen to use only Bugey in these fits as the measurement has the lowest combined errors. Also, any global fit to multiple reactor datasets must correctly account for the correlated systematics between them, which is beyond the scope of our fits at present.This dataset is referred to as Bugey in the analysis below. As shown in Figure3, this dataset presents a signal at 95% CL. There are 60 bins in this analysis in total, each extending from 1 to 6 MeV in positron energy: the 15 m and 45 m baselines contributing 25 bins each and the 90 m baseline contributing 10 bins. The fit follows the “normalized energy spectra” fit method and χ2 definition detailed in [59]. The χ2 definition depends not only on the mass and mixing parameters we fit for, but also on five large-scale deformations of the positron spectrum due to systematic effects. Energy resolution and baseline smearing due to the finite reactor core are taken into account. To fold in the flux normalization correction mentioned above, we update the theoretical prediction for the expected ratio by an overall normalization factor of 1.06237, 1.06197, and 1.0627 for the 15 m, 45 m, and 90 m baselines, respectively [23].Gallium Calibration Dataset. Indications of νe disappearance have recently been published from calibration data taken by the SAGE [39] and GALLEX [40] experiments. These were solar neutrino experiments that used Mega-curie sources of 51Cr and 37Ar, which produce νe, to calibrate the detectors. Each of the two experiments had two calibration periods. The overall rates from these four measurements are consistent and show an overall deficit that has been reported to be consistent with electron flavor disappearance [27, 60]. We use the four ratios of calibration data to expectation, as reported in [27], Table 2: 1.00±0.10, 0.81±0.10, 0.95±0.12, and 0.79±0.10. These correspond to the two periods from GALLEX and the two periods from SAGE, respectively. Our analysis of this dataset, referred to as Gallium below, follows that of [27]; a 4-bin fit to the above measured calibration period rates is used. The predicted rates after oscillations are obtained by averaging the oscillation probabilities taking into account the detector geometry, the location of the source within the detector, and the neutrino energy distribution for each source (energy line and branching fraction). The neutrino energies are approximately 430 and 750 keV for 51Cr and 812 keV for 37Ar. The data result in a limit at 95% CL, as shown in Figure 3.


#### 3.2.5. Long-Baseline Experimental Results Contributing to the Fits

While this study concentrates mainly on results from SBL experiments, the data from experiments with baselines of hundreds of kilometers can be valuable. At such long baselines, the ability to identify theΔm2 associated with any observed oscillation has disappeared due to the rapid oscillations. However, these experiments can place strong constraints on the mixing parameters. New to this paper is the inclusion of the MINOS ν-μ CC constraint. We have included the atmospheric dataset in our previous fits [20–22].We note two long-baseline results not included in this analysis. First, we have dropped the Chooz dataset that was included in previous fits [20–22] due to the discovery that sin22θ13 is large [12–16], which significantly complicates the use of this data for SBL oscillation searches. Second, the recent muon flavor disappearance results from IceCube [61] were published too late to be included in this iteration of fits. However, the MiniBooNE and MINOS muon flavor disappearance results are more stringent than the IceCube limits, and so we do not expect this to significantly affect the results.MINOS
ν
-
μ
CC Disappearance Search. MINOS is a muon flavor disappearance experiment featuring two (near and far) iron-scintillator segmented calorimeter-style detectors in the NuMI beamline (described above) at Fermilab. The near detector is located 1 km from the target while the far detector is located 730 km away. The wideband beam peaks at about 4 GeV.MINOS ran in both neutrino and antineutrino mode. We employ the antineutrino data in our fits as it constrains the allowed region for muon antineutrino disappearance when we divide the datasets into neutrino versus antineutrino fits. The MINOS neutrino mode disappearance limit is not as restrictive as the atmospheric result, and so only the antineutrino dataset is utilized.This result is referred to as MINOS-CC in the analysis below. The data present a limit at 95% CL as discussed above and shown in Figure2. In our analysis of MINOS-CC, we fit both the antineutrino (right sign) data published by MINOS in antineutrino mode running [62] and the antineutrino (wrong sign) data published by MINOS in neutrino mode running [63]. The right sign data are considered in 12 bins from 0 to 20 GeV, and the wrong sign in 13 bins from 0 to 50 GeV. We account for possible oscillations in the near detector due to high Δm2 values by using the ratio of the oscillation probabilities at the far and near detectors for each mass and mixing model. As MINOS is sensitive to Δmatm2, we add an extra mass state to the oscillation probability using the best-fit atmospheric mass and mixing parameters from the MINOS experiment [10]. The data points and systematic errors are taken from [62, 63].Atmospheric Constraints on
ν
μ
Disappearance Used in Fits. Atmospheric neutrinos are produced when cosmic rays interact with nuclei in the atmosphere to produce showers of mesons. The neutrino path length varies from a few to 12,800 km, while neutrino energies range from sub- to few-GeV. Thus, this is a long-baseline source with sensitivity to primarily νμ disappearance and effectively no sensitivity to Δmij2. The former is a consequence of the atmospheric neutrino flux composition and the detector technology used in atmospheric experiments. Thus, atmospheric neutrino measurements and long-baseline accelerator-based νμ disappearance experiments constrain the same parameters and are treated in our fits in a similar way.As with our past fits, we include atmospheric constraints following the prescription of [64]. We refer to this dataset as ATM. This makes use of two datasets: (1) 1489 days of Super-K muon-like and electron-like events with energies in the sub- to multi-GeV range, taking into account atmospheric flux predictions from [65] and treating systematic uncertainties according to [66]; (2) νμ disappearance data from the long-baseline accelerator-based experiment K2K [6, 67, 68]. The atmospheric constraint is implemented in the form of a χ2 available as a function of the parameter dμ, which depends on the muon flavor composition of m4, m5, and m6 as follows:
(14)dμ=1-1-4A2,
where
(15)A=(1-|Uμ4|2-|Uμ5|2-|Uμ6|2)×(|Uμ4|2+|Uμ5|2+|Uμ6|2)+|Uμ4|2|Uμ5|2+|Uμ4|2|Uμ6|2+|Uμ5|2|Uμ6|2.
The atmospheric constraints set a limit at 95% CL as shown in Figure 2.


## 3.1. Sources and Detectors Used in Short-Baseline Neutrino Experiments

Before considering the datasets in detail, we provide an overview of how SBL experiments are typically designed.


### 3.1.1. Sources of Neutrinos for Short-Baseline Experiments

The neutrino sources used in SBL experiments range in energy from a few MeV to hundreds of GeV and include manmade radioactive sources, reactors, and accelerator-produced beams. While the higher energy accelerator sources are mixtures of different neutrino flavors, the <10 MeV sources rely on beta decay and are thus pure electron neutrino flavor.At the low-energy end of the spectrum, the rate of electron neutrino interactions from the beta decay of the ~1 MCi sources51Cr (half-life: 28 days) and 37Ar (half-life: 35 days) have been studied. These sources were originally produced for the low-energy (~1 MeV) calibration of solar neutrino detectors [39, 40] but have proven themselves interesting as a probe of electron neutrino disappearance.Moving up in energy by a few MeV, nuclear reactors are powerful sources of ~2−8 MeVν-e through the β+-decaying elements produced primarily in the decay chains of 235U, 239Pu, 238U, and 241Pu. While these four isotopes are the progenitors of most of the reactor flux, modern reactor simulations include all fission sources [41]. Reactor simulations convolute predictions of fission rates over time with neutrino production per fission. Recently, a reanalysis of the production cross-section per fission [23, 42, 43] has led to an increase in the predicted reactor flux. As their energy is too low for an appearance search (the neutrino energy is below the muon production kinematic threshold), reactor source antineutrinos  can only be used for ν-e disappearance searches, where the antineutrinos are detected using CC interactions with an outgoing e+.The lowest neutrino energy (up to 53 MeV) accelerator sources used in existing SBL experiments are based on pion- and muon-decay-at-rest (DAR). The neutrino flux comes from the stopped pion decay chain:π+→μ+νμ and μ+→e+ν-μνe. Pions are produced in interactions of accelerator protons with, typically, a graphite or water target. The contribution from the decay chain π-→μ-ν-μ is suppressed by designing the target such that the π- mesons are captured with high probability. The result is a source which has a well-understood neutrino flavor content and energy distribution, with a minimal (<10−3) ν-e content [44, 45]. This last point is important as  ν-μ→ν-e is the dominant channel used for oscillation searches by DAR sources.In a conventional high-energy (from ~100 MeV to hundreds of GeV) accelerator-based neutrino beam, protons impinge on a target (beryllium and carbon are typical) to produce secondary mesons. The boosted mesons enter and subsequently decay inside a long, often evacuated, pipe. Neutrinos are primarily produced byπ+ and π- decay in flight (DIF). Pion sign selection, via a large magnet placed directly in the beamlines before the decay pipe, allows for nearly pure neutrino or antineutrino running, with only a few percent “wrong sign” neutrino flux content in the case of neutrino running, and ~15% [46] in the case of antineutrino running. These beams are generally produced by protons at 8 GeV and above. At these energies, in addition to pion production, kaon production contributes to the flux of both muon and electron neutrino flavors. There is often a substantial muon DIF content as well, contributing both νe/ν-e and ν-μ/νμ to the beam. The result of the kaon and muon secondary content is that, while the neutrinos are predominantly muon flavored, the beam will always have some intrinsic electron flavor neutrino content, usually at the several percent level. Accelerator-based beams are predominantly used for νμ→νe and ν-μ→ν-e appearance searches, as well as νμ and ν-μ disappearance searches. An excellent review of methods in producing accelerator-based neutrino beams can be found in [47].In contrast to lower-energy neutrino sources (DAR, reactor, and isotope sources), high-energy accelerator-based neutrino sources are subject to significant energy-dependent neutrino flux uncertainties, often at the level of 10–15%, due to in-target meson production uncertainties. These uncertainties can affect the energy distribution, flavor content, and absolute normalization of a neutrino beam. Typically, meson production systematics are constrained with dedicated measurements by experiments such as HARP [48] and MIPP [49], which use replicated targets (geometry and material) and a wide range of proton beam energies to study meson production cross-sections and kinematics directly. Alternatively, experiments can employ a two-detector design for comparing near-to-far event rate in energy to effectively reduce these systematics. However, due to the short baselines employed for studying sterile neutrino oscillations, a two-detector search is often impractical. In situ measurements in single-detector experiments can exploit flux (multiplied by cross-section) correlations among different beam components and energies to reduce flux uncertainties, as has been done in the case of the MiniBooNE νe and ν-e appearance searches described below.


### 3.1.2. Short-Baseline Neutrino Detectors

Because low-energy neutrino interaction cross-sections are very small, the options for SBL detectors are typically limited to designs which can be constructed on a massive scale. There are several generic neutrino detection methods in use today: unsegmented scintillator detectors, unsegmented Cerenkov detectors, segmented scintillator-and-iron calorimeters, and segmented trackers.Neutrino oscillation experiments usually require sensitivity to CC neutrino interactions, whereby one can definitively identify the flavor of the interacting neutrino by the presence of a charged lepton in the final state. However, in the case of sterile neutrino oscillation searches, NC interactions can also provide useful information, as they are directly sensitive to the sterile flavor content of the neutrino mass eigenstate,|Usi|2=1-|Uei|2-|Uμi|2-|Uτi|2.Unsegmented scintillator detectors are typically used for few-MeV-scale SBL experiments, which require efficient electron neutrino identification and reconstruction. These detectors consist of large tanks of oil-based (CnH2n) liquid scintillator surrounded by phototubes. The free protons in the oil provide a target for the inverse beta decay interaction, ν-ep→e+n. The reaction threshold for this interaction is 1.8 MeV due to the mass difference between the proton and neutron and the mass of the positron. The scintillation light from the e+, as well as light from the Compton scattering of the 0.511 MeV annihilation photons provides an initial (“prompt”) signal. This is followed by n capture on hydrogen and a 2.2 MeV flash of light, as the resulting γ Compton-scatters in the scintillator. This coincidence sequence in time (positron followed by neutron capture) provides a clean, mostly background-free interaction signature. Experiments often dope the liquid scintillator using an element with a high neutron capture cross-section for improved event identification efficiency.The CC interaction with the carbon in the oil (which produces either nitrogen or boron depending on whether the scatterer is a neutrino or antineutrino) has a significantly higher energy threshold than the free proton target-scattering process. The CC quasielastic interactionνe+C→e-+N has an energy threshold of 17.3 MeV, which arises from the carbon-nitrogen mass difference and the mass of the electron. In the case of both reactor and radioactive decay sources, the flux cuts off below this energy threshold. However, neutrinos from DAR sources are at sufficiently high energy to produce these carbon scatters.Unsegmented Cerenkov detectors make use of a target which is a large volume of clear medium (undoped oil or water is typical) surrounded by, or interspersed with, phototubes. Undoped oil has a larger refractive index, leading to a larger Cerenkov opening angle. Water is the only affordable medium once the detector size surpasses a few kilotons. In this paper, the only unsegmented Cerenkov detector that is considered is the 450-ton oil-based MiniBooNE detector. In such a detector, a track will project a ring with a sharp inner and outer edge onto the phototubes. Consider an electron produced in aνe CC quasielastic interaction. As the electron is low mass, it will multiple-scatter and easily bremsstrahlung, smearing the light projected on the tubes and producing a “fuzzy” ring. A muon produced by a CC quasielastic νμ interaction (νμn→μ-p) is heavier and will thus produce a sharper outer edge to the ring. For the same visible energy, the track will also extend farther, filling the interior of the ring and, perhaps, exiting the tank. If the muon stops within the tank and subsequently decays, the resulting electron provides an added tag for particle identification. In the case of the μ-, 18% will capture in water and, thus, have no electron tag, while only 8% will capture in the oil.Scintillator and iron calorimeters provide an affordable detection technique for ~1 GeV and higherνμ interactions. At these energies, multiple hadrons may be produced at the interaction vertex and will be observed as hadronic showers. In these devices, the iron provides the target, while the scintillator provides information on energy deposition per unit length. This information allows separation between the hadronic shower, which occurs in both NC and CC events, and the minimum-ionizing track of an outgoing muon, which occurs in CC events. Transverse information can be obtained if segmented scintillator strips are used, or if drift chambers are interspersed. The light from scintillator strips is transported to tubes by wavelength-shifting fibers. Information in the transverse plane improves separation of electromagnetic and hadronic showers. The iron can be magnetized to allow separation of neutrino and antineutrino events based on the charge of the outgoing lepton.To address the problem of running at ~1 GeV, where hadron track reconstruction is desirable, highly segmented tracking designs have been developed. The best resolution comes from stacks of wire chambers, where the material enclosing the gas provides the target. However, a more practical alternative has been stacks of thin extruded scintillator bars that are read out using wavelength-shifting fibers.


## 3.1.1. Sources of Neutrinos for Short-Baseline Experiments

The neutrino sources used in SBL experiments range in energy from a few MeV to hundreds of GeV and include manmade radioactive sources, reactors, and accelerator-produced beams. While the higher energy accelerator sources are mixtures of different neutrino flavors, the <10 MeV sources rely on beta decay and are thus pure electron neutrino flavor.At the low-energy end of the spectrum, the rate of electron neutrino interactions from the beta decay of the ~1 MCi sources51Cr (half-life: 28 days) and 37Ar (half-life: 35 days) have been studied. These sources were originally produced for the low-energy (~1 MeV) calibration of solar neutrino detectors [39, 40] but have proven themselves interesting as a probe of electron neutrino disappearance.Moving up in energy by a few MeV, nuclear reactors are powerful sources of ~2−8 MeVν-e through the β+-decaying elements produced primarily in the decay chains of 235U, 239Pu, 238U, and 241Pu. While these four isotopes are the progenitors of most of the reactor flux, modern reactor simulations include all fission sources [41]. Reactor simulations convolute predictions of fission rates over time with neutrino production per fission. Recently, a reanalysis of the production cross-section per fission [23, 42, 43] has led to an increase in the predicted reactor flux. As their energy is too low for an appearance search (the neutrino energy is below the muon production kinematic threshold), reactor source antineutrinos  can only be used for ν-e disappearance searches, where the antineutrinos are detected using CC interactions with an outgoing e+.The lowest neutrino energy (up to 53 MeV) accelerator sources used in existing SBL experiments are based on pion- and muon-decay-at-rest (DAR). The neutrino flux comes from the stopped pion decay chain:π+→μ+νμ and μ+→e+ν-μνe. Pions are produced in interactions of accelerator protons with, typically, a graphite or water target. The contribution from the decay chain π-→μ-ν-μ is suppressed by designing the target such that the π- mesons are captured with high probability. The result is a source which has a well-understood neutrino flavor content and energy distribution, with a minimal (<10−3) ν-e content [44, 45]. This last point is important as  ν-μ→ν-e is the dominant channel used for oscillation searches by DAR sources.In a conventional high-energy (from ~100 MeV to hundreds of GeV) accelerator-based neutrino beam, protons impinge on a target (beryllium and carbon are typical) to produce secondary mesons. The boosted mesons enter and subsequently decay inside a long, often evacuated, pipe. Neutrinos are primarily produced byπ+ and π- decay in flight (DIF). Pion sign selection, via a large magnet placed directly in the beamlines before the decay pipe, allows for nearly pure neutrino or antineutrino running, with only a few percent “wrong sign” neutrino flux content in the case of neutrino running, and ~15% [46] in the case of antineutrino running. These beams are generally produced by protons at 8 GeV and above. At these energies, in addition to pion production, kaon production contributes to the flux of both muon and electron neutrino flavors. There is often a substantial muon DIF content as well, contributing both νe/ν-e and ν-μ/νμ to the beam. The result of the kaon and muon secondary content is that, while the neutrinos are predominantly muon flavored, the beam will always have some intrinsic electron flavor neutrino content, usually at the several percent level. Accelerator-based beams are predominantly used for νμ→νe and ν-μ→ν-e appearance searches, as well as νμ and ν-μ disappearance searches. An excellent review of methods in producing accelerator-based neutrino beams can be found in [47].In contrast to lower-energy neutrino sources (DAR, reactor, and isotope sources), high-energy accelerator-based neutrino sources are subject to significant energy-dependent neutrino flux uncertainties, often at the level of 10–15%, due to in-target meson production uncertainties. These uncertainties can affect the energy distribution, flavor content, and absolute normalization of a neutrino beam. Typically, meson production systematics are constrained with dedicated measurements by experiments such as HARP [48] and MIPP [49], which use replicated targets (geometry and material) and a wide range of proton beam energies to study meson production cross-sections and kinematics directly. Alternatively, experiments can employ a two-detector design for comparing near-to-far event rate in energy to effectively reduce these systematics. However, due to the short baselines employed for studying sterile neutrino oscillations, a two-detector search is often impractical. In situ measurements in single-detector experiments can exploit flux (multiplied by cross-section) correlations among different beam components and energies to reduce flux uncertainties, as has been done in the case of the MiniBooNE νe and ν-e appearance searches described below.


## 3.1.2. Short-Baseline Neutrino Detectors

Because low-energy neutrino interaction cross-sections are very small, the options for SBL detectors are typically limited to designs which can be constructed on a massive scale. There are several generic neutrino detection methods in use today: unsegmented scintillator detectors, unsegmented Cerenkov detectors, segmented scintillator-and-iron calorimeters, and segmented trackers.Neutrino oscillation experiments usually require sensitivity to CC neutrino interactions, whereby one can definitively identify the flavor of the interacting neutrino by the presence of a charged lepton in the final state. However, in the case of sterile neutrino oscillation searches, NC interactions can also provide useful information, as they are directly sensitive to the sterile flavor content of the neutrino mass eigenstate,|Usi|2=1-|Uei|2-|Uμi|2-|Uτi|2.Unsegmented scintillator detectors are typically used for few-MeV-scale SBL experiments, which require efficient electron neutrino identification and reconstruction. These detectors consist of large tanks of oil-based (CnH2n) liquid scintillator surrounded by phototubes. The free protons in the oil provide a target for the inverse beta decay interaction, ν-ep→e+n. The reaction threshold for this interaction is 1.8 MeV due to the mass difference between the proton and neutron and the mass of the positron. The scintillation light from the e+, as well as light from the Compton scattering of the 0.511 MeV annihilation photons provides an initial (“prompt”) signal. This is followed by n capture on hydrogen and a 2.2 MeV flash of light, as the resulting γ Compton-scatters in the scintillator. This coincidence sequence in time (positron followed by neutron capture) provides a clean, mostly background-free interaction signature. Experiments often dope the liquid scintillator using an element with a high neutron capture cross-section for improved event identification efficiency.The CC interaction with the carbon in the oil (which produces either nitrogen or boron depending on whether the scatterer is a neutrino or antineutrino) has a significantly higher energy threshold than the free proton target-scattering process. The CC quasielastic interactionνe+C→e-+N has an energy threshold of 17.3 MeV, which arises from the carbon-nitrogen mass difference and the mass of the electron. In the case of both reactor and radioactive decay sources, the flux cuts off below this energy threshold. However, neutrinos from DAR sources are at sufficiently high energy to produce these carbon scatters.Unsegmented Cerenkov detectors make use of a target which is a large volume of clear medium (undoped oil or water is typical) surrounded by, or interspersed with, phototubes. Undoped oil has a larger refractive index, leading to a larger Cerenkov opening angle. Water is the only affordable medium once the detector size surpasses a few kilotons. In this paper, the only unsegmented Cerenkov detector that is considered is the 450-ton oil-based MiniBooNE detector. In such a detector, a track will project a ring with a sharp inner and outer edge onto the phototubes. Consider an electron produced in aνe CC quasielastic interaction. As the electron is low mass, it will multiple-scatter and easily bremsstrahlung, smearing the light projected on the tubes and producing a “fuzzy” ring. A muon produced by a CC quasielastic νμ interaction (νμn→μ-p) is heavier and will thus produce a sharper outer edge to the ring. For the same visible energy, the track will also extend farther, filling the interior of the ring and, perhaps, exiting the tank. If the muon stops within the tank and subsequently decays, the resulting electron provides an added tag for particle identification. In the case of the μ-, 18% will capture in water and, thus, have no electron tag, while only 8% will capture in the oil.Scintillator and iron calorimeters provide an affordable detection technique for ~1 GeV and higherνμ interactions. At these energies, multiple hadrons may be produced at the interaction vertex and will be observed as hadronic showers. In these devices, the iron provides the target, while the scintillator provides information on energy deposition per unit length. This information allows separation between the hadronic shower, which occurs in both NC and CC events, and the minimum-ionizing track of an outgoing muon, which occurs in CC events. Transverse information can be obtained if segmented scintillator strips are used, or if drift chambers are interspersed. The light from scintillator strips is transported to tubes by wavelength-shifting fibers. Information in the transverse plane improves separation of electromagnetic and hadronic showers. The iron can be magnetized to allow separation of neutrino and antineutrino events based on the charge of the outgoing lepton.To address the problem of running at ~1 GeV, where hadron track reconstruction is desirable, highly segmented tracking designs have been developed. The best resolution comes from stacks of wire chambers, where the material enclosing the gas provides the target. However, a more practical alternative has been stacks of thin extruded scintillator bars that are read out using wavelength-shifting fibers.


## 3.2. Data Used in the Sterile Neutrino Fits

There are many SBL datasets that can be included in this analysis. In this work, we have substantially expanded the number of datasets used beyond those in our past papers [20–22]. In the sections following, we identify and discuss new and updated datasets, as well as provide information on those used in past fits. The fit technique is described in Section 4.2.


### 3.2.1. Experimental Results from Decay at Rest Studies

In past sterile neutrino studies [20–22], we have included the LSND and KARMEN appearance results described below. Since that work, a new study that constrains νe disappearance from the relative LSND-to-KARMEN cross-section measurements was published [50]. This new dataset is included in this analysis.LSND Appearance. LSND was a DAR experiment that ran in the 1990s, searching for ν-μ→ν-e. The beam was produced using 800 MeV protons on target from the LAMPF accelerator at Los Alamos National Laboratory, where a 1 mA beam of protons impinged on a water target. The center of the 8.75 m long, nearly cylindrical detector, was located at 29.8 m from the target, at an angle of 12° from the proton beam direction. This was an unsegmented detector with a fiducial mass of 167 tons of oil (CH2), lightly doped with b-PBD scintillator. The intrinsic ν-e content of the beam was 8×10-4 of the ν-μ content. The experiment observed a ν-e excess of 87.9±22.4±6.0 events above background, which was interpreted as oscillations with a probability of (0.264±0.067±0.045)%. Details are available in [17].This dataset is referred to as LSND in the analysis below and indicates a signal at 95% CL, as shown in Figure1. This data covers energies between 20 and 53 MeV and contributes five energy bins to the global fit. Statistical errors are taken into account by using a log-likelihood χ2 definition in the fit, while systematic errors on the background prediction are not included because these are small relative to the statistical error. Energy and baseline smearing are taken into account by averaging the oscillation probability over the energy bin width and over the neutrino flight path uncertainty.KARMEN Appearance. KARMEN was another DAR experiment searching for ν-μ→ν-e. KARMEN ran at the ISIS facility at Rutherford Laboratory, with 200 μA of protons impinging on a copper, tantalum, or uranium target. The neutrino detector was located at an angle of 100° with respect to the targeting protons to reduce background from π- DIF. The resulting intrinsic ν-e content was 6.4×10-4 of the ν-μ content.The center of the approximately cubic segmented scintillator detector was located at 17.7 m. Thus, this detector was 60% of the distance from the source compared to LSND. The liquid scintillator target volume was 56 m3 and consisted of 512 optically independent modules (17.4 cm × 17.8 cm × 353 cm) wrapped in gadolinium-doped paper. KARMEN saw no signal and set a limit on appearance. More details are available in [51].This dataset is referred to as KARMEN in the analysis below and indicates a limit at 95% CL, as shown in Figure1. This dataset contributes nine energy bins, in the range 16 to 50 MeV. As in the case of LSND, statistical errors are taken into account by using a log-likelihood χ2 definition in the fit, while systematic errors on the background prediction are not included. Energy and baseline smearing are taken into account by averaging the sin2(1.27Δm2L/E) and sin(2.53Δm2L/E) term contributions in the total signal prediction over energy bin widths. The limit which is shown here is determined using a Δχ2-based raster scan, as discussed in Section 4.2.LSND and KARMEN Cross-Section Measurements. Along with the oscillation searches, LSND and KARMEN measured νe+12C→12Ngs+e- scattering. In this two-body interaction, with a Q-value of 17.3 MeV, the neutrino energy can be reconstructed by measuring the outgoing visible energy of the electron. The 12N ground state is identified by the subsequent β decay, 12Ngs→12C+e++νe, which has a Q-value of 16.3 MeV and a lifetime of 15.9 ms.The cross-section is measured by both experiments under the assumption that theνe flux has not oscillated, leading to disappearance. The excellent agreement between the two results, as a function of energy, allows a limit to be placed on νe oscillations. The energy dependence of the cross-section, as well as the normalization, are well predicted and both constraints are used in the analysis [50].This dataset is referred to as KARMEN/LSND(xsec) in the analysis below, and indicates a limit at 95% CL, as shown in Figure 3. A total of six (for KARMEN) plus five (for LSND) bins are used in the fit, which extend approximately from 28–50 MeV in the case of KARMEN and from 38–50 MeV in the case of LSND. In calculating the oscillation probability, the signal is averaged across the lengths of the detectors. The experiments have correlated systematics arising from the flux normalization due to a shared underlying analysis for pion production in DAR experiments. This is addressed through application of pull terms as described in [50].


### 3.2.2. The MiniBooNE Experimental Results

The MiniBooNE experiment provides multiple results from a single detector. This oil-based 450 t fiducial volume Cerenkov detector was exposed to two conventional beams, the Booster Neutrino Beam (BNB) and the off-axis NuMI beam. The primary goal of MiniBooNE was to search forνμ→νe and ν-μ→ν-e appearance, using the BNB, which provides sensitivity to Δm2~0.1–10eV2 oscillations. The NuMI beam also provides some sensitivity to νμ→νe appearance at a similar Δm2. In addition to the appearance searches, MiniBooNE also looked for νμ and ν-μ disappearance using the BNB.The MiniBooNE datasets included in our analysis have increased throughout the period that our group has been performing fits. Reference [21] used a Monte Carlo prediction for neutrinos and antineutrinos to estimate MiniBooNE’s sensitivity to sterile neutrinos. The full BNB neutrino and first published BNB antineutrino datasets from MiniBooNE form the experimental constraints in [22]. Here, we have updated the analysis to include the full BNB antineutrino datasets. A further update has been to employ a log-likelihood method for the BNB neutrino and antineutrino datasets from [18], as this was recently adopted by the MiniBooNE Collaboration [52]. We also use the updated constraints on electron neutrino flux from kaons [18]. A partial dataset from NuMI data taking was presented in [22] and has not been updated, as the result was already systematics limited. In this analysis we also introduce the MiniBooNE disappearance search [53].In our fits to MiniBooNE appearance data, when drawing allowed regions and calculating compatibilities, which make use ofΔχ2’s and not absolute χ2’s, we use MiniBooNE’s log-likelihood χ2 definition, summing over both νe and νμ bins, as described in [52]. For consistency, the absolute MiniBooNE BNB νe and ν-e appearance χ2 values quoted in our paper also correspond to the same definition, that is, fitting to both νe and νμ spectra; therefore, they differ from the ones published by MiniBooNE in [18], which are obtained by fitting only to a priori constrained νe distributions. Note that the two definitions yield consistent allowed regions and compatibility results.The Booster Neutrino Beam Appearance Search in Neutrino Running Mode. The BNB flux composition in neutrino mode consists of >90% νμ, 6%ν-μ, and 0.06% νe and ν-e combined [46]. In the MiniBooNE BNB search for νe appearance, the νe and ν-e signal was normalized to the νμ and ν-μ CC quasielastic events observed in the detector, which peaked at 700 MeV.The global fits presented here use the full statistics of the MiniBooNEνμ→νe dataset, representing 6.46×1020 protons on target. In this dataset, MiniBooNE has observed an excess of events at 200–1250 MeV, corresponding to 162.0±47.8 electron-like events [18]. The dataset is referred to as BNB-MB(νapp) in the analysis below.We include the BNB-MB(νapp) dataset in our fits in the form of the full νe CC reconstructed energy distribution, in 11 energy bins from 200 to 3000 MeV, fit simultaneously with the full νμ CC energy distribution, in eight energy bins up to 1900 MeV. We account for statistical and systematic uncertainties in each sample, as well as systematic correlations (from flux and cross-section) among the νe signal and background and νμ background distributions. The systematic correlations are provided in the form of a full 19-bin×19-bin fractional covariance matrix. By fitting the νe and νμ spectra simultaneously, we are able to exploit the high-statistics νμ CC sample as a constraint on background and signal event rates. This assumes no significant νμ disappearance; this simplification could lead to a <20% effect on appearance probability obtained in MiniBooNE only fits [18].The dataset results in a signal at 95% CL, as shown in Figure1. This has changed slightly from our past analysis [22] now that we are using updated constraints on intrinsic electron neutrinos from kaons and the log-likelihood method, but is in agreement with the equivalent analysis from the MiniBooNE Collaboration [18].The Booster Neutrino Beam Appearance Search in Antineutrino Running Mode. The BNB flux composition in antineutrino mode consists of 83% ν-μ, 0.6% νe and ν-e combined, and a significantly larger wrong-sign composition than in neutrino mode, of 16%  νμ. As in the BNB νe appearance search, the electron flavor signal was normalized to the muon flavor CC quasielastic events observed in the detector, which peaked at 500 MeV.The global fits presented here use the full statistics of the MiniBooNEν-μ dataset, representing 11.27×1020 protons on target. In this dataset, MiniBooNE has observed an excess of events at 200–1250MeV, corresponding to 78.4±28.5 electron-like events. The dataset is referred to as BNB-MB(ν-app) in the analysis below.As in neutrino mode, we fit the fullν-e CC energy distribution, in 11 energy bins from 200 to 3000 MeV, simultaneously with the full ν-μ CC energy distribution, in 8 energy bins up to 1900 MeV. The wrong-sign contamination in the beam (νμ) is assumed to not contribute to any oscillations; only ν-μ→ν-e oscillations are assumed for this dataset. We account for statistical and systematic uncertainties in each sample as well as systematic correlations among the ν-e and ν-μ distributions in the form of a full 19-bin × 19-bin fractional covariance matrix in each fit. For further information, see [18].The dataset results in a signal at 95% CL, as shown in Figure1.The NuMI Beam Appearance Search. The MiniBooNE detector is also exposed to the NuMI neutrino beam, produced from a 120 GeV proton beam impinging on a carbon target. This beam is nominally used for the MINOS long-baseline neutrino oscillation experiment. NuMI events arrive out of time with the BNB-produced events. This 200 MeV to 3 GeV neutrino energy source is dominated by kaon decays near the NuMI target, which is 110 mrad off-axis and located 745 m upstream of the MiniBooNE detector. The beam consists of 81% νμ, 13% ν-μ, 5% νe, and 1% ν-e. For more information on this data, see [54].This dataset is referred to as NuMI-MB(νapp) in the analysis below. As seen in Figure 1, the dataset provides a limit at 95% CL. In the fits presented here, this data is used to constrain electron flavor appearance in neutrino mode, with 10 bins used in the fit. Statistical and systematic errors for this dataset are added in quadrature.The Booster Neutrino Beam Disappearance Search. The MiniBooNE experiment also searched for νμ and ν-μ disappearance using the BNB. The neutrino (antineutrino) dataset corresponded to 5.6×1020 (3.4×1020) protons on target, which produced a beam covering the neutrino energy range up to 1.9 GeV. The MiniBooNE νμ disappearance result provides restrictions on sterile neutrino oscillations which are comparable to those provided by the CDHS experiment, discussed below. Therefore, we include that dataset in these fits. On the other hand, the ν-μ result was weaker due to the combination of fewer protons on target and a lower cross-section. The MINOS ν-μ CC constraint, described below, is stronger, and so we do not use the MiniBooNE ν-μ dataset.The fit to theνμ dataset uses 16 bins ranging up to 1900 MeV in reconstructed neutrino energy. A shape-only fit is performed, where the predicted spectrum given any set of oscillation parameters is renormalized so that the total number of predicted events, after oscillations, is equal to the total number of observed events. Then the normalized predicted spectrum is compared to the observed spectrum in the form of a χ2 which accounts for statistical and shape-only systematic uncertainties and bin-to-bin correlations in the form of a covariance matrix.This dataset is referred to as BNB-MB(νdis) in the analysis below. Figure 2 shows that this data sets a limit at 95% CL. It should be noted that the published MiniBooNE analysis used a Pearson χ2 method [53], and we are able to reproduce those results. However, to fold these results into our analysis, we reverted to the Δχ2 definition used consistently among all datasets included in the fits (see Section 4.2).


### 3.2.3. Results from Multi-GeV Conventional Short-Baselineνμ Beams

The set of multi-GeV conventional SBLνμ experiments is the same as was used in previous fits. Our overview of these experiments is therefore very brief.NOMAD Appearance Search. The NOMAD experiment [55], which ran at CERN using protons from the 450 GeV SPS accelerator, employed a conventional neutrino beamline to create a wideband 2.5 to 40 GeV neutrino energy source. These neutrinos were created with a carbon-based, low-mass tracking detector located 600 m downstream of the target. This detector had fine spatial resolution and could search for muon-to-electron and muon-to-tau oscillations. No signal was observed in either channel. In this analysis, we use the νμ→νe constraint.This dataset is referred to as NOMAD in the analysis below. This dataset contributes 30 energy bins to the global fit. The statistical and systematic errors are added in quadrature. This experiment sets a limit at 95% CL, as seen in Figure1.CCFR Disappearance Search. The CCFR dataset was taken at Fermilab in 1984 [56] with a narrowband beamline, with meson energies set to 100, 140, 165, 200, and 250 GeV, yielding νμ and ν-μ beams that ranged from 40 to 230 GeV in energy. This was a two-detector disappearance search, with the near detector at 715 m and the far detector at 1116 m from the center of the 352 m long decay pipe. The calorimetric detectors were constructed of segmented iron with scintillator and spark chambers, and each had a downstream toroid to measure the muon momentum.This dataset is referred to as CCFR84 in the analysis below. The data were published as the double ratios of the observed-to-expected rates in a near-to-far ratio. For each secondary mean setting, the data are divided into three energy bins. The systematic uncertainty is assumed to be energy independent and fully correlated between the energy bins. Due to the high beam energies and short baselines, this experiment sets a limit at highΔm2 in the muon flavor disappearance search at 95% CL, as shown in Figure 2.CDHS Disappearance Search. The CDHS experiment [57] at CERN searched for νμ disappearance with a two-detector design of segmented calorimeters with iron and scintillator. The experiment used 19.2 GeV protons on a beryllium target to produce mesons that were subsequently focused into a 52 m decay channel. The detectors were located 130 m and 885 m downstream of the target.This dataset is referred to as CDHS in the analysis below. CDHS provides data and errors in 15 bins of muon energy, as seen in in Table 1 of [57]. We relate these bins to the neutrino energy distributions using the method described in [20]: the neutrino energy distribution for a given muon energy or range is determined via the NUANCE [58] neutrino event generator. The experiment has a limit at 95% CL and sets constraints that are comparable to the MiniBooNE νμ disappearance limit described above, but which extend to slightly lower Δm2. See Figure 2 for comparison.


### 3.2.4. Reactor and Source Experiments

The reactor experiment dataset has been updated to reflect recent changes in the predicted neutrino fluxes, as discussed below. The source-based experimental datasets are both new to this paper, and have been published since our last set of fits [22].Bugey Dataset. This analysis uses energy-dependent data from the Bugey 3 reactor experiment [59]. The detector consisted of 6Li-doped liquid scintillator, with data taken at 15, 45, and 90 m from the 2.8 GW reactor source. The detectors are taken to be pointlike in the analysis.Recently, a reanalysis of reactorν-e flux predictions [23, 42, 43] has led to a reinterpretation of the Bugey data. The data has transitioned from being simply a limit on neutrino disappearance to an allowed region at 95% CL. In this analysis, we adjust the predicted Bugey flux spectra normalization according to the calculations from [23].There are many other SBL reactor datasets in existence. However, we have chosen to use only Bugey in these fits as the measurement has the lowest combined errors. Also, any global fit to multiple reactor datasets must correctly account for the correlated systematics between them, which is beyond the scope of our fits at present.This dataset is referred to as Bugey in the analysis below. As shown in Figure3, this dataset presents a signal at 95% CL. There are 60 bins in this analysis in total, each extending from 1 to 6 MeV in positron energy: the 15 m and 45 m baselines contributing 25 bins each and the 90 m baseline contributing 10 bins. The fit follows the “normalized energy spectra” fit method and χ2 definition detailed in [59]. The χ2 definition depends not only on the mass and mixing parameters we fit for, but also on five large-scale deformations of the positron spectrum due to systematic effects. Energy resolution and baseline smearing due to the finite reactor core are taken into account. To fold in the flux normalization correction mentioned above, we update the theoretical prediction for the expected ratio by an overall normalization factor of 1.06237, 1.06197, and 1.0627 for the 15 m, 45 m, and 90 m baselines, respectively [23].Gallium Calibration Dataset. Indications of νe disappearance have recently been published from calibration data taken by the SAGE [39] and GALLEX [40] experiments. These were solar neutrino experiments that used Mega-curie sources of 51Cr and 37Ar, which produce νe, to calibrate the detectors. Each of the two experiments had two calibration periods. The overall rates from these four measurements are consistent and show an overall deficit that has been reported to be consistent with electron flavor disappearance [27, 60]. We use the four ratios of calibration data to expectation, as reported in [27], Table 2: 1.00±0.10, 0.81±0.10, 0.95±0.12, and 0.79±0.10. These correspond to the two periods from GALLEX and the two periods from SAGE, respectively. Our analysis of this dataset, referred to as Gallium below, follows that of [27]; a 4-bin fit to the above measured calibration period rates is used. The predicted rates after oscillations are obtained by averaging the oscillation probabilities taking into account the detector geometry, the location of the source within the detector, and the neutrino energy distribution for each source (energy line and branching fraction). The neutrino energies are approximately 430 and 750 keV for 51Cr and 812 keV for 37Ar. The data result in a limit at 95% CL, as shown in Figure 3.


### 3.2.5. Long-Baseline Experimental Results Contributing to the Fits

While this study concentrates mainly on results from SBL experiments, the data from experiments with baselines of hundreds of kilometers can be valuable. At such long baselines, the ability to identify theΔm2 associated with any observed oscillation has disappeared due to the rapid oscillations. However, these experiments can place strong constraints on the mixing parameters. New to this paper is the inclusion of the MINOS ν-μ CC constraint. We have included the atmospheric dataset in our previous fits [20–22].We note two long-baseline results not included in this analysis. First, we have dropped the Chooz dataset that was included in previous fits [20–22] due to the discovery that sin22θ13 is large [12–16], which significantly complicates the use of this data for SBL oscillation searches. Second, the recent muon flavor disappearance results from IceCube [61] were published too late to be included in this iteration of fits. However, the MiniBooNE and MINOS muon flavor disappearance results are more stringent than the IceCube limits, and so we do not expect this to significantly affect the results.MINOS
ν
-
μ
CC Disappearance Search. MINOS is a muon flavor disappearance experiment featuring two (near and far) iron-scintillator segmented calorimeter-style detectors in the NuMI beamline (described above) at Fermilab. The near detector is located 1 km from the target while the far detector is located 730 km away. The wideband beam peaks at about 4 GeV.MINOS ran in both neutrino and antineutrino mode. We employ the antineutrino data in our fits as it constrains the allowed region for muon antineutrino disappearance when we divide the datasets into neutrino versus antineutrino fits. The MINOS neutrino mode disappearance limit is not as restrictive as the atmospheric result, and so only the antineutrino dataset is utilized.This result is referred to as MINOS-CC in the analysis below. The data present a limit at 95% CL as discussed above and shown in Figure2. In our analysis of MINOS-CC, we fit both the antineutrino (right sign) data published by MINOS in antineutrino mode running [62] and the antineutrino (wrong sign) data published by MINOS in neutrino mode running [63]. The right sign data are considered in 12 bins from 0 to 20 GeV, and the wrong sign in 13 bins from 0 to 50 GeV. We account for possible oscillations in the near detector due to high Δm2 values by using the ratio of the oscillation probabilities at the far and near detectors for each mass and mixing model. As MINOS is sensitive to Δmatm2, we add an extra mass state to the oscillation probability using the best-fit atmospheric mass and mixing parameters from the MINOS experiment [10]. The data points and systematic errors are taken from [62, 63].Atmospheric Constraints on
ν
μ
Disappearance Used in Fits. Atmospheric neutrinos are produced when cosmic rays interact with nuclei in the atmosphere to produce showers of mesons. The neutrino path length varies from a few to 12,800 km, while neutrino energies range from sub- to few-GeV. Thus, this is a long-baseline source with sensitivity to primarily νμ disappearance and effectively no sensitivity to Δmij2. The former is a consequence of the atmospheric neutrino flux composition and the detector technology used in atmospheric experiments. Thus, atmospheric neutrino measurements and long-baseline accelerator-based νμ disappearance experiments constrain the same parameters and are treated in our fits in a similar way.As with our past fits, we include atmospheric constraints following the prescription of [64]. We refer to this dataset as ATM. This makes use of two datasets: (1) 1489 days of Super-K muon-like and electron-like events with energies in the sub- to multi-GeV range, taking into account atmospheric flux predictions from [65] and treating systematic uncertainties according to [66]; (2) νμ disappearance data from the long-baseline accelerator-based experiment K2K [6, 67, 68]. The atmospheric constraint is implemented in the form of a χ2 available as a function of the parameter dμ, which depends on the muon flavor composition of m4, m5, and m6 as follows:
(14)dμ=1-1-4A2,
where
(15)A=(1-|Uμ4|2-|Uμ5|2-|Uμ6|2)×(|Uμ4|2+|Uμ5|2+|Uμ6|2)+|Uμ4|2|Uμ5|2+|Uμ4|2|Uμ6|2+|Uμ5|2|Uμ6|2.
The atmospheric constraints set a limit at 95% CL as shown in Figure 2.


## 3.2.1. Experimental Results from Decay at Rest Studies

In past sterile neutrino studies [20–22], we have included the LSND and KARMEN appearance results described below. Since that work, a new study that constrains νe disappearance from the relative LSND-to-KARMEN cross-section measurements was published [50]. This new dataset is included in this analysis.LSND Appearance. LSND was a DAR experiment that ran in the 1990s, searching for ν-μ→ν-e. The beam was produced using 800 MeV protons on target from the LAMPF accelerator at Los Alamos National Laboratory, where a 1 mA beam of protons impinged on a water target. The center of the 8.75 m long, nearly cylindrical detector, was located at 29.8 m from the target, at an angle of 12° from the proton beam direction. This was an unsegmented detector with a fiducial mass of 167 tons of oil (CH2), lightly doped with b-PBD scintillator. The intrinsic ν-e content of the beam was 8×10-4 of the ν-μ content. The experiment observed a ν-e excess of 87.9±22.4±6.0 events above background, which was interpreted as oscillations with a probability of (0.264±0.067±0.045)%. Details are available in [17].This dataset is referred to as LSND in the analysis below and indicates a signal at 95% CL, as shown in Figure1. This data covers energies between 20 and 53 MeV and contributes five energy bins to the global fit. Statistical errors are taken into account by using a log-likelihood χ2 definition in the fit, while systematic errors on the background prediction are not included because these are small relative to the statistical error. Energy and baseline smearing are taken into account by averaging the oscillation probability over the energy bin width and over the neutrino flight path uncertainty.KARMEN Appearance. KARMEN was another DAR experiment searching for ν-μ→ν-e. KARMEN ran at the ISIS facility at Rutherford Laboratory, with 200 μA of protons impinging on a copper, tantalum, or uranium target. The neutrino detector was located at an angle of 100° with respect to the targeting protons to reduce background from π- DIF. The resulting intrinsic ν-e content was 6.4×10-4 of the ν-μ content.The center of the approximately cubic segmented scintillator detector was located at 17.7 m. Thus, this detector was 60% of the distance from the source compared to LSND. The liquid scintillator target volume was 56 m3 and consisted of 512 optically independent modules (17.4 cm × 17.8 cm × 353 cm) wrapped in gadolinium-doped paper. KARMEN saw no signal and set a limit on appearance. More details are available in [51].This dataset is referred to as KARMEN in the analysis below and indicates a limit at 95% CL, as shown in Figure1. This dataset contributes nine energy bins, in the range 16 to 50 MeV. As in the case of LSND, statistical errors are taken into account by using a log-likelihood χ2 definition in the fit, while systematic errors on the background prediction are not included. Energy and baseline smearing are taken into account by averaging the sin2(1.27Δm2L/E) and sin(2.53Δm2L/E) term contributions in the total signal prediction over energy bin widths. The limit which is shown here is determined using a Δχ2-based raster scan, as discussed in Section 4.2.LSND and KARMEN Cross-Section Measurements. Along with the oscillation searches, LSND and KARMEN measured νe+12C→12Ngs+e- scattering. In this two-body interaction, with a Q-value of 17.3 MeV, the neutrino energy can be reconstructed by measuring the outgoing visible energy of the electron. The 12N ground state is identified by the subsequent β decay, 12Ngs→12C+e++νe, which has a Q-value of 16.3 MeV and a lifetime of 15.9 ms.The cross-section is measured by both experiments under the assumption that theνe flux has not oscillated, leading to disappearance. The excellent agreement between the two results, as a function of energy, allows a limit to be placed on νe oscillations. The energy dependence of the cross-section, as well as the normalization, are well predicted and both constraints are used in the analysis [50].This dataset is referred to as KARMEN/LSND(xsec) in the analysis below, and indicates a limit at 95% CL, as shown in Figure 3. A total of six (for KARMEN) plus five (for LSND) bins are used in the fit, which extend approximately from 28–50 MeV in the case of KARMEN and from 38–50 MeV in the case of LSND. In calculating the oscillation probability, the signal is averaged across the lengths of the detectors. The experiments have correlated systematics arising from the flux normalization due to a shared underlying analysis for pion production in DAR experiments. This is addressed through application of pull terms as described in [50].


## 3.2.2. The MiniBooNE Experimental Results

The MiniBooNE experiment provides multiple results from a single detector. This oil-based 450 t fiducial volume Cerenkov detector was exposed to two conventional beams, the Booster Neutrino Beam (BNB) and the off-axis NuMI beam. The primary goal of MiniBooNE was to search forνμ→νe and ν-μ→ν-e appearance, using the BNB, which provides sensitivity to Δm2~0.1–10eV2 oscillations. The NuMI beam also provides some sensitivity to νμ→νe appearance at a similar Δm2. In addition to the appearance searches, MiniBooNE also looked for νμ and ν-μ disappearance using the BNB.The MiniBooNE datasets included in our analysis have increased throughout the period that our group has been performing fits. Reference [21] used a Monte Carlo prediction for neutrinos and antineutrinos to estimate MiniBooNE’s sensitivity to sterile neutrinos. The full BNB neutrino and first published BNB antineutrino datasets from MiniBooNE form the experimental constraints in [22]. Here, we have updated the analysis to include the full BNB antineutrino datasets. A further update has been to employ a log-likelihood method for the BNB neutrino and antineutrino datasets from [18], as this was recently adopted by the MiniBooNE Collaboration [52]. We also use the updated constraints on electron neutrino flux from kaons [18]. A partial dataset from NuMI data taking was presented in [22] and has not been updated, as the result was already systematics limited. In this analysis we also introduce the MiniBooNE disappearance search [53].In our fits to MiniBooNE appearance data, when drawing allowed regions and calculating compatibilities, which make use ofΔχ2’s and not absolute χ2’s, we use MiniBooNE’s log-likelihood χ2 definition, summing over both νe and νμ bins, as described in [52]. For consistency, the absolute MiniBooNE BNB νe and ν-e appearance χ2 values quoted in our paper also correspond to the same definition, that is, fitting to both νe and νμ spectra; therefore, they differ from the ones published by MiniBooNE in [18], which are obtained by fitting only to a priori constrained νe distributions. Note that the two definitions yield consistent allowed regions and compatibility results.The Booster Neutrino Beam Appearance Search in Neutrino Running Mode. The BNB flux composition in neutrino mode consists of >90% νμ, 6%ν-μ, and 0.06% νe and ν-e combined [46]. In the MiniBooNE BNB search for νe appearance, the νe and ν-e signal was normalized to the νμ and ν-μ CC quasielastic events observed in the detector, which peaked at 700 MeV.The global fits presented here use the full statistics of the MiniBooNEνμ→νe dataset, representing 6.46×1020 protons on target. In this dataset, MiniBooNE has observed an excess of events at 200–1250 MeV, corresponding to 162.0±47.8 electron-like events [18]. The dataset is referred to as BNB-MB(νapp) in the analysis below.We include the BNB-MB(νapp) dataset in our fits in the form of the full νe CC reconstructed energy distribution, in 11 energy bins from 200 to 3000 MeV, fit simultaneously with the full νμ CC energy distribution, in eight energy bins up to 1900 MeV. We account for statistical and systematic uncertainties in each sample, as well as systematic correlations (from flux and cross-section) among the νe signal and background and νμ background distributions. The systematic correlations are provided in the form of a full 19-bin×19-bin fractional covariance matrix. By fitting the νe and νμ spectra simultaneously, we are able to exploit the high-statistics νμ CC sample as a constraint on background and signal event rates. This assumes no significant νμ disappearance; this simplification could lead to a <20% effect on appearance probability obtained in MiniBooNE only fits [18].The dataset results in a signal at 95% CL, as shown in Figure1. This has changed slightly from our past analysis [22] now that we are using updated constraints on intrinsic electron neutrinos from kaons and the log-likelihood method, but is in agreement with the equivalent analysis from the MiniBooNE Collaboration [18].The Booster Neutrino Beam Appearance Search in Antineutrino Running Mode. The BNB flux composition in antineutrino mode consists of 83% ν-μ, 0.6% νe and ν-e combined, and a significantly larger wrong-sign composition than in neutrino mode, of 16%  νμ. As in the BNB νe appearance search, the electron flavor signal was normalized to the muon flavor CC quasielastic events observed in the detector, which peaked at 500 MeV.The global fits presented here use the full statistics of the MiniBooNEν-μ dataset, representing 11.27×1020 protons on target. In this dataset, MiniBooNE has observed an excess of events at 200–1250MeV, corresponding to 78.4±28.5 electron-like events. The dataset is referred to as BNB-MB(ν-app) in the analysis below.As in neutrino mode, we fit the fullν-e CC energy distribution, in 11 energy bins from 200 to 3000 MeV, simultaneously with the full ν-μ CC energy distribution, in 8 energy bins up to 1900 MeV. The wrong-sign contamination in the beam (νμ) is assumed to not contribute to any oscillations; only ν-μ→ν-e oscillations are assumed for this dataset. We account for statistical and systematic uncertainties in each sample as well as systematic correlations among the ν-e and ν-μ distributions in the form of a full 19-bin × 19-bin fractional covariance matrix in each fit. For further information, see [18].The dataset results in a signal at 95% CL, as shown in Figure1.The NuMI Beam Appearance Search. The MiniBooNE detector is also exposed to the NuMI neutrino beam, produced from a 120 GeV proton beam impinging on a carbon target. This beam is nominally used for the MINOS long-baseline neutrino oscillation experiment. NuMI events arrive out of time with the BNB-produced events. This 200 MeV to 3 GeV neutrino energy source is dominated by kaon decays near the NuMI target, which is 110 mrad off-axis and located 745 m upstream of the MiniBooNE detector. The beam consists of 81% νμ, 13% ν-μ, 5% νe, and 1% ν-e. For more information on this data, see [54].This dataset is referred to as NuMI-MB(νapp) in the analysis below. As seen in Figure 1, the dataset provides a limit at 95% CL. In the fits presented here, this data is used to constrain electron flavor appearance in neutrino mode, with 10 bins used in the fit. Statistical and systematic errors for this dataset are added in quadrature.The Booster Neutrino Beam Disappearance Search. The MiniBooNE experiment also searched for νμ and ν-μ disappearance using the BNB. The neutrino (antineutrino) dataset corresponded to 5.6×1020 (3.4×1020) protons on target, which produced a beam covering the neutrino energy range up to 1.9 GeV. The MiniBooNE νμ disappearance result provides restrictions on sterile neutrino oscillations which are comparable to those provided by the CDHS experiment, discussed below. Therefore, we include that dataset in these fits. On the other hand, the ν-μ result was weaker due to the combination of fewer protons on target and a lower cross-section. The MINOS ν-μ CC constraint, described below, is stronger, and so we do not use the MiniBooNE ν-μ dataset.The fit to theνμ dataset uses 16 bins ranging up to 1900 MeV in reconstructed neutrino energy. A shape-only fit is performed, where the predicted spectrum given any set of oscillation parameters is renormalized so that the total number of predicted events, after oscillations, is equal to the total number of observed events. Then the normalized predicted spectrum is compared to the observed spectrum in the form of a χ2 which accounts for statistical and shape-only systematic uncertainties and bin-to-bin correlations in the form of a covariance matrix.This dataset is referred to as BNB-MB(νdis) in the analysis below. Figure 2 shows that this data sets a limit at 95% CL. It should be noted that the published MiniBooNE analysis used a Pearson χ2 method [53], and we are able to reproduce those results. However, to fold these results into our analysis, we reverted to the Δχ2 definition used consistently among all datasets included in the fits (see Section 4.2).


## 3.2.3. Results from Multi-GeV Conventional Short-Baselineνμ Beams

The set of multi-GeV conventional SBLνμ experiments is the same as was used in previous fits. Our overview of these experiments is therefore very brief.NOMAD Appearance Search. The NOMAD experiment [55], which ran at CERN using protons from the 450 GeV SPS accelerator, employed a conventional neutrino beamline to create a wideband 2.5 to 40 GeV neutrino energy source. These neutrinos were created with a carbon-based, low-mass tracking detector located 600 m downstream of the target. This detector had fine spatial resolution and could search for muon-to-electron and muon-to-tau oscillations. No signal was observed in either channel. In this analysis, we use the νμ→νe constraint.This dataset is referred to as NOMAD in the analysis below. This dataset contributes 30 energy bins to the global fit. The statistical and systematic errors are added in quadrature. This experiment sets a limit at 95% CL, as seen in Figure1.CCFR Disappearance Search. The CCFR dataset was taken at Fermilab in 1984 [56] with a narrowband beamline, with meson energies set to 100, 140, 165, 200, and 250 GeV, yielding νμ and ν-μ beams that ranged from 40 to 230 GeV in energy. This was a two-detector disappearance search, with the near detector at 715 m and the far detector at 1116 m from the center of the 352 m long decay pipe. The calorimetric detectors were constructed of segmented iron with scintillator and spark chambers, and each had a downstream toroid to measure the muon momentum.This dataset is referred to as CCFR84 in the analysis below. The data were published as the double ratios of the observed-to-expected rates in a near-to-far ratio. For each secondary mean setting, the data are divided into three energy bins. The systematic uncertainty is assumed to be energy independent and fully correlated between the energy bins. Due to the high beam energies and short baselines, this experiment sets a limit at highΔm2 in the muon flavor disappearance search at 95% CL, as shown in Figure 2.CDHS Disappearance Search. The CDHS experiment [57] at CERN searched for νμ disappearance with a two-detector design of segmented calorimeters with iron and scintillator. The experiment used 19.2 GeV protons on a beryllium target to produce mesons that were subsequently focused into a 52 m decay channel. The detectors were located 130 m and 885 m downstream of the target.This dataset is referred to as CDHS in the analysis below. CDHS provides data and errors in 15 bins of muon energy, as seen in in Table 1 of [57]. We relate these bins to the neutrino energy distributions using the method described in [20]: the neutrino energy distribution for a given muon energy or range is determined via the NUANCE [58] neutrino event generator. The experiment has a limit at 95% CL and sets constraints that are comparable to the MiniBooNE νμ disappearance limit described above, but which extend to slightly lower Δm2. See Figure 2 for comparison.


## 3.2.4. Reactor and Source Experiments

The reactor experiment dataset has been updated to reflect recent changes in the predicted neutrino fluxes, as discussed below. The source-based experimental datasets are both new to this paper, and have been published since our last set of fits [22].Bugey Dataset. This analysis uses energy-dependent data from the Bugey 3 reactor experiment [59]. The detector consisted of 6Li-doped liquid scintillator, with data taken at 15, 45, and 90 m from the 2.8 GW reactor source. The detectors are taken to be pointlike in the analysis.Recently, a reanalysis of reactorν-e flux predictions [23, 42, 43] has led to a reinterpretation of the Bugey data. The data has transitioned from being simply a limit on neutrino disappearance to an allowed region at 95% CL. In this analysis, we adjust the predicted Bugey flux spectra normalization according to the calculations from [23].There are many other SBL reactor datasets in existence. However, we have chosen to use only Bugey in these fits as the measurement has the lowest combined errors. Also, any global fit to multiple reactor datasets must correctly account for the correlated systematics between them, which is beyond the scope of our fits at present.This dataset is referred to as Bugey in the analysis below. As shown in Figure3, this dataset presents a signal at 95% CL. There are 60 bins in this analysis in total, each extending from 1 to 6 MeV in positron energy: the 15 m and 45 m baselines contributing 25 bins each and the 90 m baseline contributing 10 bins. The fit follows the “normalized energy spectra” fit method and χ2 definition detailed in [59]. The χ2 definition depends not only on the mass and mixing parameters we fit for, but also on five large-scale deformations of the positron spectrum due to systematic effects. Energy resolution and baseline smearing due to the finite reactor core are taken into account. To fold in the flux normalization correction mentioned above, we update the theoretical prediction for the expected ratio by an overall normalization factor of 1.06237, 1.06197, and 1.0627 for the 15 m, 45 m, and 90 m baselines, respectively [23].Gallium Calibration Dataset. Indications of νe disappearance have recently been published from calibration data taken by the SAGE [39] and GALLEX [40] experiments. These were solar neutrino experiments that used Mega-curie sources of 51Cr and 37Ar, which produce νe, to calibrate the detectors. Each of the two experiments had two calibration periods. The overall rates from these four measurements are consistent and show an overall deficit that has been reported to be consistent with electron flavor disappearance [27, 60]. We use the four ratios of calibration data to expectation, as reported in [27], Table 2: 1.00±0.10, 0.81±0.10, 0.95±0.12, and 0.79±0.10. These correspond to the two periods from GALLEX and the two periods from SAGE, respectively. Our analysis of this dataset, referred to as Gallium below, follows that of [27]; a 4-bin fit to the above measured calibration period rates is used. The predicted rates after oscillations are obtained by averaging the oscillation probabilities taking into account the detector geometry, the location of the source within the detector, and the neutrino energy distribution for each source (energy line and branching fraction). The neutrino energies are approximately 430 and 750 keV for 51Cr and 812 keV for 37Ar. The data result in a limit at 95% CL, as shown in Figure 3.


## 3.2.5. Long-Baseline Experimental Results Contributing to the Fits

While this study concentrates mainly on results from SBL experiments, the data from experiments with baselines of hundreds of kilometers can be valuable. At such long baselines, the ability to identify theΔm2 associated with any observed oscillation has disappeared due to the rapid oscillations. However, these experiments can place strong constraints on the mixing parameters. New to this paper is the inclusion of the MINOS ν-μ CC constraint. We have included the atmospheric dataset in our previous fits [20–22].We note two long-baseline results not included in this analysis. First, we have dropped the Chooz dataset that was included in previous fits [20–22] due to the discovery that sin22θ13 is large [12–16], which significantly complicates the use of this data for SBL oscillation searches. Second, the recent muon flavor disappearance results from IceCube [61] were published too late to be included in this iteration of fits. However, the MiniBooNE and MINOS muon flavor disappearance results are more stringent than the IceCube limits, and so we do not expect this to significantly affect the results.MINOS
ν
-
μ
CC Disappearance Search. MINOS is a muon flavor disappearance experiment featuring two (near and far) iron-scintillator segmented calorimeter-style detectors in the NuMI beamline (described above) at Fermilab. The near detector is located 1 km from the target while the far detector is located 730 km away. The wideband beam peaks at about 4 GeV.MINOS ran in both neutrino and antineutrino mode. We employ the antineutrino data in our fits as it constrains the allowed region for muon antineutrino disappearance when we divide the datasets into neutrino versus antineutrino fits. The MINOS neutrino mode disappearance limit is not as restrictive as the atmospheric result, and so only the antineutrino dataset is utilized.This result is referred to as MINOS-CC in the analysis below. The data present a limit at 95% CL as discussed above and shown in Figure2. In our analysis of MINOS-CC, we fit both the antineutrino (right sign) data published by MINOS in antineutrino mode running [62] and the antineutrino (wrong sign) data published by MINOS in neutrino mode running [63]. The right sign data are considered in 12 bins from 0 to 20 GeV, and the wrong sign in 13 bins from 0 to 50 GeV. We account for possible oscillations in the near detector due to high Δm2 values by using the ratio of the oscillation probabilities at the far and near detectors for each mass and mixing model. As MINOS is sensitive to Δmatm2, we add an extra mass state to the oscillation probability using the best-fit atmospheric mass and mixing parameters from the MINOS experiment [10]. The data points and systematic errors are taken from [62, 63].Atmospheric Constraints on
ν
μ
Disappearance Used in Fits. Atmospheric neutrinos are produced when cosmic rays interact with nuclei in the atmosphere to produce showers of mesons. The neutrino path length varies from a few to 12,800 km, while neutrino energies range from sub- to few-GeV. Thus, this is a long-baseline source with sensitivity to primarily νμ disappearance and effectively no sensitivity to Δmij2. The former is a consequence of the atmospheric neutrino flux composition and the detector technology used in atmospheric experiments. Thus, atmospheric neutrino measurements and long-baseline accelerator-based νμ disappearance experiments constrain the same parameters and are treated in our fits in a similar way.As with our past fits, we include atmospheric constraints following the prescription of [64]. We refer to this dataset as ATM. This makes use of two datasets: (1) 1489 days of Super-K muon-like and electron-like events with energies in the sub- to multi-GeV range, taking into account atmospheric flux predictions from [65] and treating systematic uncertainties according to [66]; (2) νμ disappearance data from the long-baseline accelerator-based experiment K2K [6, 67, 68]. The atmospheric constraint is implemented in the form of a χ2 available as a function of the parameter dμ, which depends on the muon flavor composition of m4, m5, and m6 as follows:
(14)dμ=1-1-4A2,
where
(15)A=(1-|Uμ4|2-|Uμ5|2-|Uμ6|2)×(|Uμ4|2+|Uμ5|2+|Uμ6|2)+|Uμ4|2|Uμ5|2+|Uμ4|2|Uμ6|2+|Uμ5|2|Uμ6|2.
The atmospheric constraints set a limit at 95% CL as shown in Figure 2.


## 4. Analysis Description

The analysis method follows the formalism described in Section2.3, and fits are performed to each of the (3+1), (3+2), and (3+3) hypotheses separately.


### 4.1. Fit Parameters

The independent parameters considered in the (3+1) fit are Δm412, representing the splitting between the (degenerate) first three mass eigenstates and the fourth mass eigenstate, and |Ue4| and |Uμ4|, representing the electron and muon flavor content in the fourth mass eigenstate, which are assumed to be small. The (3+2) model introduces a fifth mass eigenstate (where Δm512≥Δm412) two additional mixing parameters” |Ue5| and |Uμ5|, and the CP-violating phase ϕ54, defined by (10). The (3+3) model includes all the previous parameters and a sixth mass eigenstate, described by Δm612, where Δm612≥Δm512≥Δm412, two additional mixing parameters, |Ue6| and |Uμ6|, and two more CP-violating phases, ϕ46 and ϕ56. The above model parameters are allowed to vary freely within the following ranges: Δm412, Δm512, and Δm612 within 0.01–100 eV2; |Uαi| within 0.01–0.5; ϕij within 0–2π, with the exception that for the |Uαi| there are additional constraints imposed on the mixing parameters in order to conserve unitarity of the full (3+N)×(3+N) mixing matrix in each scenario, as described in Section 2.3.


### 4.2. Fitting Method

The fitting method closely follows what has been done in [21]. Given an oscillation model, (3+1), (3+2), or (3+3), the corresponding independent oscillation parameters are randomly generated within their allowed range, and then varied via a Markov Chain χ2 minimization procedure [69]. Each independent parameter x is generated and varied according to
(16)x=xold+s(R-0.5)(xmin-xmax),
where xold is the value of parameter x previously tested in the χ2 minimization chain; xmin and xmax represent the boundaries on the parameter x as described in Section 4.1; R is a random number between 0 and 1, which is varied as one steps from  xold to x and  s is the stepsize, a parameter of the Markov Chain. By definition, within the Markov Chain minimization method  the point is accepted based only on the point directly preceding it. The acceptance of any new point x in the chain, where x is the new point in the oscillation parameter space, is determined by the following:
(17)P=min(1,e-(χ2-χold2)/T),
where T is the Markov Chain parameter “temperature.” The stepsize and temperature control how quickly the Markov Chain diffuses toward the minimum χ2 value. At every step in the chain, each of which corresponds to a point in the oscillation parameter space, a χ2 is calculated by summing together the individual χd2 contributed from each dataset, d, included in the fit, where d denotes a dataset as described in Section 3.2.In any given fit, we define possible signal indications at 90% and 99% CL by marginalizing over the full parameter space, and looking for closed contours formed about a global minimum,χmin2, when projected onto any two-dimensional parameter space, assuming only two degrees of freedom. We use the standard two degrees of freedom Δχ2 cuts of 4.61 for exploring allowed 90% CL regions, 5.99 for exploring allowed 95% CL regions (used only for Figures 1, 2, and 3), and 9.21 for 99% CL regions. If the null point (Uαi,Uβi=0) is allowed at >95% CL, we instead proceed with drawing one-dimensional raster scan limits, obtained with the standard Δχ2 cuts of 2.70, 3.84, and 6.63 for 90%, 95% (used only for Figures 1, 2, and 3), and 99% CL, respectively.


### 4.3. Parameter Goodness-of-Fit Test

In any given fit, in addition to a standardχ2-probability, which is quoted for the global χmin2 and number of degrees of freedom in the fit, we also report statistical compatibility comparisons using the parameter goodness-of-fit test (PG test) from [70]. This test reduces the bias imposed toward datasets with a large number of bins in the standard χ2-probability in order to calculate the compatibility between datasets simply on the basis of preferred parameters. The PG (%) can be calculated to quantify compatibility between any two or more datasets, or between combinations of datasets, according to
(18)χPG2=χmin, combined2-∑i‍χmin,d2,
where χmin,combined2 is the χ2-minimum of the combined fit of the datasets in consideration, and χmin,d2 is the χ2-minimum of each dataset fit individually. When comparing groups of datasets (i.e., appearance experiments versus disappearance experiments), each group is treated as an individual dataset. The number of degrees of freedom (ndfPG) for the PG test is given by
(19)ndfPG=∑d‍Npd-Npcombined.
Here, Npd represents the number of independent parameters involved in the fit of a particular dataset and Npcombined represents the number of independent parameters involved in the global fit.


## 4.1. Fit Parameters

The independent parameters considered in the (3+1) fit are Δm412, representing the splitting between the (degenerate) first three mass eigenstates and the fourth mass eigenstate, and |Ue4| and |Uμ4|, representing the electron and muon flavor content in the fourth mass eigenstate, which are assumed to be small. The (3+2) model introduces a fifth mass eigenstate (where Δm512≥Δm412) two additional mixing parameters” |Ue5| and |Uμ5|, and the CP-violating phase ϕ54, defined by (10). The (3+3) model includes all the previous parameters and a sixth mass eigenstate, described by Δm612, where Δm612≥Δm512≥Δm412, two additional mixing parameters, |Ue6| and |Uμ6|, and two more CP-violating phases, ϕ46 and ϕ56. The above model parameters are allowed to vary freely within the following ranges: Δm412, Δm512, and Δm612 within 0.01–100 eV2; |Uαi| within 0.01–0.5; ϕij within 0–2π, with the exception that for the |Uαi| there are additional constraints imposed on the mixing parameters in order to conserve unitarity of the full (3+N)×(3+N) mixing matrix in each scenario, as described in Section 2.3.


## 4.2. Fitting Method

The fitting method closely follows what has been done in [21]. Given an oscillation model, (3+1), (3+2), or (3+3), the corresponding independent oscillation parameters are randomly generated within their allowed range, and then varied via a Markov Chain χ2 minimization procedure [69]. Each independent parameter x is generated and varied according to
(16)x=xold+s(R-0.5)(xmin-xmax),
where xold is the value of parameter x previously tested in the χ2 minimization chain; xmin and xmax represent the boundaries on the parameter x as described in Section 4.1; R is a random number between 0 and 1, which is varied as one steps from  xold to x and  s is the stepsize, a parameter of the Markov Chain. By definition, within the Markov Chain minimization method  the point is accepted based only on the point directly preceding it. The acceptance of any new point x in the chain, where x is the new point in the oscillation parameter space, is determined by the following:
(17)P=min(1,e-(χ2-χold2)/T),
where T is the Markov Chain parameter “temperature.” The stepsize and temperature control how quickly the Markov Chain diffuses toward the minimum χ2 value. At every step in the chain, each of which corresponds to a point in the oscillation parameter space, a χ2 is calculated by summing together the individual χd2 contributed from each dataset, d, included in the fit, where d denotes a dataset as described in Section 3.2.In any given fit, we define possible signal indications at 90% and 99% CL by marginalizing over the full parameter space, and looking for closed contours formed about a global minimum,χmin2, when projected onto any two-dimensional parameter space, assuming only two degrees of freedom. We use the standard two degrees of freedom Δχ2 cuts of 4.61 for exploring allowed 90% CL regions, 5.99 for exploring allowed 95% CL regions (used only for Figures 1, 2, and 3), and 9.21 for 99% CL regions. If the null point (Uαi,Uβi=0) is allowed at >95% CL, we instead proceed with drawing one-dimensional raster scan limits, obtained with the standard Δχ2 cuts of 2.70, 3.84, and 6.63 for 90%, 95% (used only for Figures 1, 2, and 3), and 99% CL, respectively.


## 4.3. Parameter Goodness-of-Fit Test

In any given fit, in addition to a standardχ2-probability, which is quoted for the global χmin2 and number of degrees of freedom in the fit, we also report statistical compatibility comparisons using the parameter goodness-of-fit test (PG test) from [70]. This test reduces the bias imposed toward datasets with a large number of bins in the standard χ2-probability in order to calculate the compatibility between datasets simply on the basis of preferred parameters. The PG (%) can be calculated to quantify compatibility between any two or more datasets, or between combinations of datasets, according to
(18)χPG2=χmin, combined2-∑i‍χmin,d2,
where χmin,combined2 is the χ2-minimum of the combined fit of the datasets in consideration, and χmin,d2 is the χ2-minimum of each dataset fit individually. When comparing groups of datasets (i.e., appearance experiments versus disappearance experiments), each group is treated as an individual dataset. The number of degrees of freedom (ndfPG) for the PG test is given by
(19)ndfPG=∑d‍Npd-Npcombined.
Here, Npd represents the number of independent parameters involved in the fit of a particular dataset and Npcombined represents the number of independent parameters involved in the global fit.


## 5. Results

This section presents the results of the analysis for the (3+1), (3+2), and (3+3) sterile neutrino model fits. For reference, information about the datasets used in the analyses is provided in Table 1. Tables 2 and 3 summarize the results of the fits, which will be described in more detail below. Table 2 gives the fit results for the overall global fits and for various combinations of datasets. When interpreting compatibilities, one should keep in mind that, along with a high compatibility among the individual datasets in a global fit, high compatibility values among groups of datasets is also important. Finally, Table 3 provides the parameters for the best-fit points for each of the models.Table 1
Datasets used in the fits and their corresponding use in the analysis. Column 1 provides the tag for the data. Column 2 references the description in Section3.2. Column 3 lists the relevant oscillation process. Column 4 lists which datasets are included in the neutrino versus antineutrino analyses and column 5 lists which datasets are included in the appearance versus disappearance study.
Tag
Section
Process
ν versus ν-
App versus Dis
LSND
3.2.1
ν
-
μ
→
ν
-
e
ν
-
App
KARMEN
3.2.1
ν
-
μ
→
ν
-
e
ν
-
App
KARMEN/LSND(xsec)
3.2.1
ν
e
→
ν
e
ν
Dis
BNB-MB(νapp)
3.2.2
ν
μ
→
ν
e
ν
App
BNB-MB(ν-app)
3.2.2
ν
-
μ
→
ν
-
e
ν
-
App
NuMI-MB(νapp)
3.2.2
ν
μ
→
ν
e
ν
App
BNB-MB(νdis)
3.2.2
ν
μ
→
ν
μ
ν
Dis
NOMAD
3.2.3
ν
μ
→
ν
e
ν
App
CCFR84
3.2.3
ν
μ
→
ν
μ
ν
Dis
CDHS
3.2.3
ν
μ
→
ν
μ
ν
Dis
Bugey
3.2.4
ν
-
e
→
ν
-
e
ν
-
Dis
Gallium
3.2.4
ν
e
→
ν
e
ν
Dis
MINOS-CC
3.2.5
ν
-
μ
→
ν
-
μ
ν
-
Dis
ATM
3.2.5
ν
μ
→
ν
μ
ν
DisTable 2
Theχ2 values, degrees of freedom (dof), and probabilities associated with the best-fit and null hypothesis in each scenario. Also shown are the results from the parameter goodness-of-fit tests. Pbest refers to the χ2-probability at the best-fit point and Pnull refers to the χ2-probability at null.
χ
min
2 (dof)
χ
null
2 (dof)
P
best
P
null
χ
PG
2 (dof)
PG (%)
3
+
1
All
233.9 (237)
286.5 (240)
55%
2.1%
54.0 (24)
0.043%
App
87.8 (87)
147.3 (90)
46%
0.013%
14.1 (9)
12%
Dis
128.2 (147)
139.3 (150)
87%
72%
22.1(19)
28%
ν
123.5 (120)
133.4 (123)
39%
25%
26.6(14)
2.2%
ν
¯
94.8 (114)
153.1 (117)
90%
1.4%
11.8(7)
11%
App versus Dis
—
—
—
—
17.8(2)
0.013%
ν versus ν¯
—
—
—
—
15.6(3)
0.14%
3
+
2
All
221.5 (233)
286.5 (240)
69%
2.1%
63.8 (52)
13%
App
75.0 (85)
147.3 (90)
77%
0.013%
16.3 (25)
90%
Dis
122.6 (144)
139.3 (150)
90%
72%
23.6 (23)
43%
ν
116.8 (116)
133.4 (123)
77%
25%
35.0 (29)
21%
ν
¯
90.8 (110)
153.1 (117)
90%
1.4%
15.0(16)
53%
App versus Dis
—
—
—
—
23.9(4)
0.0082%
ν versus ν¯
—
—
—
—
13.9(7)
5.3%
3
+
3
All
218.2 (228)
286.5 (240)
67%
2.1%
68.9 (85)
90%
App
70.8 (81)
147.3 (90)
78%
0.013%
17.6 (45)
100%
Dis
120.3 (141)
139.3 (150)
90%
72%
24.1 (34)
90%
ν
116.7 (111)
133.4 (123)
34%
25%
39.5 (46)
74%
ν
¯
90.6 (105)
153 (117)
84%
1.4%
18.5 (27)
89%
App versus Dis
—
—
—
—
27.1(6)
0.014%
ν versus ν¯
—
—
—
—
10.9 (12)
53%Table 3
The oscillation parameter best-fit points in each scenario considered. The values ofΔm2 shown are in units of eV2.
(a)
3
+
1
Δ
m
41
2
|
U
μ
4
|
|
U
e
4
|
All
0.92
0.17
0.15
App
0.15
0.39
0.39
Dis
18
0.18
0.18
ν
7.8
0.059
0.26
ν
¯
0.92
0.23
0.13
(b)
3
+
2
Δ
m
41
2
Δ
m
51
2
|
U
μ
4
|
|
U
e
4
|
|
U
μ
5
|
|
U
e
5
|
ϕ
54
All
0.92
17
0.13
0.15
0.16
0.069
1.8π
App
0.31
1.0
0.31
0.31
0.17
0.17
1.1π
Dis
0.92
18
0.015
0.12
0.17
0.12
N/A
ν
7.6
17.6
0.05
0.27
0.18
0.052
1.8π
ν
¯
0.92
3.8
0.25
0.13
0.12
0.079
0.35π
(c)
3
+
3
Δ
m
41
2
Δ
m
51
2
Δ
m
61
2
|
U
μ
4
|
|
U
e
4
|
|
U
μ
5
|
|
U
e
5
|
|
U
μ
6
|
|
U
e
6
|
ϕ
54
ϕ
64
ϕ
65
All
0.90
17
22
0.12
0.11
0.17
0.11
0.14
0.11
1.6π
0.28π
1.4π
App
0.15
1.8
2.7
0.37
0.37
0.12
0.12
0.12
0.12
1.4π
0.32π
0.94π
Dis
0.92
7.2
18
0.013
0.12
0.019
0.16
0.15
0.069
N/A
N/A
N/A
ν
13
17
26
0.076
0.24
0.16
0.067
0.10
0.017
1.1π
1.8π
0.037π
ν
¯
7.5
9.1
18
0.024
0.28
0.098
0.11
0.18
0.029
1.8π
2.0π
0.61π


### 5.1. (3+1) Fit Results

For a (3+1) model, three parameters are determined: Δm412, |Ue4|, and |Uμ4|. A global (3+1) fit of all of the experiments (Figure 4) a yields  χ2-probability of 55% but a very low compatibility of 0.043%, indicating a low compatibility among all individual datasets. Contrasting the result from the χ2 test to the poor compatibility illustrates how the χ2 test can be misleading. As discussed above, this is due to some datasets dominating others due to the number of bins in the fit, many of which may not have strong oscillation sensitivity. It is for this reason that most groups fitting for sterile neutrinos now use the PG test as the figure of merit.Figure 4
TheΔm412 versus sin22θμe allowed space from fits to all data—neutrino and antineutrino—in a (3+1) model.In order to understand the source of the poor compatibility, the datasets are subdivided, as shown in Table2, into separate neutrino and antineutrino results. Within each of these categories, the PG compatibility values are 2.2% and 11% for neutrinos and antineutrinos, respectively. However, the two datasets favor very different oscillation parameters, as seen in Table 3 and Figure 5. This leads to a very low PG of 0.14% when the neutrino and antineutrino data are compared. The separation of the datasets into appearance and also shows a strong incompatibility, as illustrated in Figure 6, leading to an even lower PG value of 0.013%. These results imply that the (3+1) model is not sufficient to describe all datasets simultaneously.TheΔm412 versus sin22θμe allowed space from fits to neutrino (a) and antineutrino (b) data in a (3+1) model.
(a)
(b)TheΔm412 versus sin22θμe allowed space from fits to appearance (a) and disappearance (b) data in a (3+1) model.
(a)
(b)As can be seen in Table3, two different Δm412 values are preferred for neutrino versus antineutrino and for appearance versus disappearance. This leads one to suspect that the data would prefer at least two mass splittings between the mostly active and the mostly sterile states and, thus, encourages the consideration of a (3+2) interpretation [20]. Moreover, a (3+2) model allows the introduction of a CP-violating phase, which can address the differences between neutrino and antineutrino datasets [21]. Therefore, these results lead us to abandon (3+1) and move on to testing the (3+2) hypothesis. It should be noted that the shortcomings of the (3+1) model have now been established by a number of independent analyses [20, 22, 23, 71–73].


### 5.2. (3+2) Fit Results

In a (3+2) model, there are seven parameters to determine: Δm412, Δm512, |Ue4|, |Uμ4|, |Ue5|, |Uμ5|, and ϕ54. The best-fit values for these parameters from a global fit to all datasets are given in Table 3. The 90% and 99% CL contours in marginalized (Δm412,Δm512) space can be seen in Figure 7.Figure 7
TheΔm512 versus Δm412 allowed space from fits to all data in a (3+2) model.Adding a second mass eigenstate reduces the tension seen in the (3+1) fits, bringing the overall compatibility to 13% (see Table 2) and reducing the χ2 of the global fit by 12.4 units with four extra parameters introduced in the fit. For this compatibility test, the BNB-MB(νapp) dataset has the worst χ2-probability. When considered by itself, the BNB-MB(νapp) dataset gives a constrained (see Section 3.2.2) χ2 (dof) of 19.2 (4) for the global best-fit parameters, which corresponds to a χ2-probability of 0.07%. This is one of the first indications that the MiniBooNE neutrino data has some tension with the other datasets.The need to introduce a CP-violating phase was established in previous studies of global fits [22]. This term affects only fits involving appearance datasets and results in a difference in the oscillation probabilities for νμ→νe versus ν-μ→ν-e. In particular, previous studies considered CP-violating fits in an attempt to reconcile the MiniBooNE neutrino appearance results with the MiniBooNE and LSND antineutrino appearance results.Table2 gives the fit results in dataset combinations for cross-comparison. We find that the separate neutrino and antineutrino dataset fits remain in good agreement and that the compatibility between them has risen to 5.3%—a significant improvement over the (3+1) result. The best-fit values and allowed regions are shown in Table 3 and in Figure 8, respectively.TheΔm512 versus Δm412 allowed space from fits to neutrino (a) and antineutrino (b) data in a (3+2) model.
(a)
(b)While the neutrino versus antineutrino discrepancy has been somewhat reduced, Table2 points out a second important problem. The appearance and disappearance datasets still have very poor compatibility (0.0082%), even in a (3+2) model. The poor compatibility can be partially traced to a discrepancy in the preferred mass splittings for these two datasets. As reported in Table 3, the appearance datasets prefer a low (0.31 eV2) and a medium (1.0 eV2) mass splitting while the disappearance datasets prefer a medium (0.92 eV2) and a high (18 eV2) splitting. This is also illustrated in Figure 9. This suggests that three mass splittings may be required to reconcile appearance and disappearance results and motivates the consideration of a (3+3) model.TheΔm512 versus Δm412 allowed space from fits to appearance (a) and disappearance (b) data in a (3+2) model.
(a)
(b)


### 5.3.  (3+3) Fit Results

For a (3+3) model, there are 12 model parameters to be determined: Δm412, Δm512, Δm612, |Ue4|, |Uμ4|, |Ue5|, |Uμ5|, |Ue6|, |Uμ6|, ϕ54, ϕ64, and ϕ65. Adding a third mass eigenstate does not significantly change the global fit χmin2; however, the tension between the individual dataset fits is further reduced, raising the compatibility from 13%, in (3+2), to 90%. The neutrino and antineutrino compatibility rises by an order of magnitude from the (3+2) value to 53%, indicating that the (3+3) model can better accommodate the differences in these datasets.It is interesting to note that the (3+3) fit prefers an “inverted hierarchy” among the three mostly sterile states, with Δm542=m52-m42=16eV2 and Δm652=m62-m52=5.0eV2. The overall splitting relative to the three mostly active states is Δm412=m42-m12=0.90eV2.The one puzzling discrepancy for the (3 + 3) fits is the tension still exhibited by the appearance versus disappearance datasets, shown in Figures10, 11, and 12, where the compatibility remains low, at less than 0.01%. We find that important sources of this incompatibility are the BNB-MB(νapp) and BNB-MB(ν-app) datasets. The BNB-MB(νapp) dataset has fairly small statistical and systematic uncertainties and therefore has a large impact on the fits and compatibility calculations. This is shown in Figure 13 where the MiniBooNE data agrees well with the appearance-only fit but disagrees with the overall global fit. Removing both the BNB-MB(νapp) and BNB-MB(ν-app) sets raises the compatibility to 3.5%, corresponding to an improvement of over two orders of magnitude. It has been known since the first MiniBooNE publication [18] that the BNB-MB(νapp) data was fairly consistent with no oscillations above 475 MeV; however, a significant low-energy excess was present below this energy. The energy dependence of the BNB-MB(νapp) excess does not fit very well with oscillation models extracted from fits to global datasets, unless very low Δmij2 with large mixing elements |Uei| and |Uμi| are involved in the fit. This may lead to the poor compatibility when included in appearance versus disappearance comparisons. Other possible explanations for this incompatibility include downward fluctuations of the BNB-MB(νapp) data in the higher-energy region or some other process contributing to the low-energy excess such as those suggested in [28, 74].TheΔm512 versus Δm412 and Δm612 versus Δm412 correlations from fits to all data in a (3+3) model.
(a)
(b)TheΔm512 versus Δm412 and Δm612 versus Δm412 correlations from fits to appearance data in a (3+3) model.
(a)
(b)TheΔm512 versus Δm412 and Δm612 versus Δm412 correlations from fits to disappearance data in a (3+3) model.
(a)
(b)A comparison of the BNB-MB(νapp) and BNB-MB(ν-app) excess data with the global best-fit oscillation signal predictions (solid colored lines) and with the appearance only best-fit predictions (dashed colored lines) for each of the models, (3+1), (3+2), and (3+3). The error bars on the excess correspond to statistical and unconstrained background systematic errors, added in quadrature.
(a)
(b)
(c)Statistical issues could be addressed with more MiniBooNE neutrino data that may become available over the next few years. In addition, the MicroBooNE experiment, which is expected to start running in 2014, will provide more information on the low-energy excess events and will answer the question of whether the excess is associated with outgoing electrons or photons [75].


### 5.4. Summary of Results

The sterile neutrino fits to global datasets show that a (3+1) model is inadequate; multiple Δmij2 values are needed along with CP-violating effects to explain the neutrino versus antineutrino differences. A (3+2) model improves significantly on the (3+1) results but still shows some tension in the neutrino versus antineutrino compatibility and cannot explain the appearance versus disappearance differences. The (3+3) model does not seem to further improve the fit and still has poor appearance versus disappearance compatibility. The BNB-MB(νapp) (and BNB-MB(ν-app)) dataset is a prime contributor to this incompatibility and additional experimental information in this region should be available soon. Figure 13 gives a comparison of the BNB-MB(νapp) and BNB-MB(ν-app) data with the global best-fit predictions and with the appearance-only best-fit predictions for each of the three models, (3+1), (3+2), and (3+3). The global fit prediction is significantly below the data at low energy, which contributes to the poor appearance versus disappearance compatibility.In summary, out of the three sterile neutrino oscillation hypotheses considered in the analysis, we find that the (3+2) and (3+3) models provide a better description than the (3+1) model, although the MiniBooNE appearance data continue to raise issues within the fits. As has been shown before, (3+1) scenarios provide a poor fit to the data, and should not be emphasized. We therefore recommend continued investigations of (3+2) and (3+3) scenarios.


## 5.1. (3+1) Fit Results

For a (3+1) model, three parameters are determined: Δm412, |Ue4|, and |Uμ4|. A global (3+1) fit of all of the experiments (Figure 4) a yields  χ2-probability of 55% but a very low compatibility of 0.043%, indicating a low compatibility among all individual datasets. Contrasting the result from the χ2 test to the poor compatibility illustrates how the χ2 test can be misleading. As discussed above, this is due to some datasets dominating others due to the number of bins in the fit, many of which may not have strong oscillation sensitivity. It is for this reason that most groups fitting for sterile neutrinos now use the PG test as the figure of merit.Figure 4
TheΔm412 versus sin22θμe allowed space from fits to all data—neutrino and antineutrino—in a (3+1) model.In order to understand the source of the poor compatibility, the datasets are subdivided, as shown in Table2, into separate neutrino and antineutrino results. Within each of these categories, the PG compatibility values are 2.2% and 11% for neutrinos and antineutrinos, respectively. However, the two datasets favor very different oscillation parameters, as seen in Table 3 and Figure 5. This leads to a very low PG of 0.14% when the neutrino and antineutrino data are compared. The separation of the datasets into appearance and also shows a strong incompatibility, as illustrated in Figure 6, leading to an even lower PG value of 0.013%. These results imply that the (3+1) model is not sufficient to describe all datasets simultaneously.TheΔm412 versus sin22θμe allowed space from fits to neutrino (a) and antineutrino (b) data in a (3+1) model.
(a)
(b)TheΔm412 versus sin22θμe allowed space from fits to appearance (a) and disappearance (b) data in a (3+1) model.
(a)
(b)As can be seen in Table3, two different Δm412 values are preferred for neutrino versus antineutrino and for appearance versus disappearance. This leads one to suspect that the data would prefer at least two mass splittings between the mostly active and the mostly sterile states and, thus, encourages the consideration of a (3+2) interpretation [20]. Moreover, a (3+2) model allows the introduction of a CP-violating phase, which can address the differences between neutrino and antineutrino datasets [21]. Therefore, these results lead us to abandon (3+1) and move on to testing the (3+2) hypothesis. It should be noted that the shortcomings of the (3+1) model have now been established by a number of independent analyses [20, 22, 23, 71–73].


## 5.2. (3+2) Fit Results

In a (3+2) model, there are seven parameters to determine: Δm412, Δm512, |Ue4|, |Uμ4|, |Ue5|, |Uμ5|, and ϕ54. The best-fit values for these parameters from a global fit to all datasets are given in Table 3. The 90% and 99% CL contours in marginalized (Δm412,Δm512) space can be seen in Figure 7.Figure 7
TheΔm512 versus Δm412 allowed space from fits to all data in a (3+2) model.Adding a second mass eigenstate reduces the tension seen in the (3+1) fits, bringing the overall compatibility to 13% (see Table 2) and reducing the χ2 of the global fit by 12.4 units with four extra parameters introduced in the fit. For this compatibility test, the BNB-MB(νapp) dataset has the worst χ2-probability. When considered by itself, the BNB-MB(νapp) dataset gives a constrained (see Section 3.2.2) χ2 (dof) of 19.2 (4) for the global best-fit parameters, which corresponds to a χ2-probability of 0.07%. This is one of the first indications that the MiniBooNE neutrino data has some tension with the other datasets.The need to introduce a CP-violating phase was established in previous studies of global fits [22]. This term affects only fits involving appearance datasets and results in a difference in the oscillation probabilities for νμ→νe versus ν-μ→ν-e. In particular, previous studies considered CP-violating fits in an attempt to reconcile the MiniBooNE neutrino appearance results with the MiniBooNE and LSND antineutrino appearance results.Table2 gives the fit results in dataset combinations for cross-comparison. We find that the separate neutrino and antineutrino dataset fits remain in good agreement and that the compatibility between them has risen to 5.3%—a significant improvement over the (3+1) result. The best-fit values and allowed regions are shown in Table 3 and in Figure 8, respectively.TheΔm512 versus Δm412 allowed space from fits to neutrino (a) and antineutrino (b) data in a (3+2) model.
(a)
(b)While the neutrino versus antineutrino discrepancy has been somewhat reduced, Table2 points out a second important problem. The appearance and disappearance datasets still have very poor compatibility (0.0082%), even in a (3+2) model. The poor compatibility can be partially traced to a discrepancy in the preferred mass splittings for these two datasets. As reported in Table 3, the appearance datasets prefer a low (0.31 eV2) and a medium (1.0 eV2) mass splitting while the disappearance datasets prefer a medium (0.92 eV2) and a high (18 eV2) splitting. This is also illustrated in Figure 9. This suggests that three mass splittings may be required to reconcile appearance and disappearance results and motivates the consideration of a (3+3) model.TheΔm512 versus Δm412 allowed space from fits to appearance (a) and disappearance (b) data in a (3+2) model.
(a)
(b)


## 5.3.  (3+3) Fit Results

For a (3+3) model, there are 12 model parameters to be determined: Δm412, Δm512, Δm612, |Ue4|, |Uμ4|, |Ue5|, |Uμ5|, |Ue6|, |Uμ6|, ϕ54, ϕ64, and ϕ65. Adding a third mass eigenstate does not significantly change the global fit χmin2; however, the tension between the individual dataset fits is further reduced, raising the compatibility from 13%, in (3+2), to 90%. The neutrino and antineutrino compatibility rises by an order of magnitude from the (3+2) value to 53%, indicating that the (3+3) model can better accommodate the differences in these datasets.It is interesting to note that the (3+3) fit prefers an “inverted hierarchy” among the three mostly sterile states, with Δm542=m52-m42=16eV2 and Δm652=m62-m52=5.0eV2. The overall splitting relative to the three mostly active states is Δm412=m42-m12=0.90eV2.The one puzzling discrepancy for the (3 + 3) fits is the tension still exhibited by the appearance versus disappearance datasets, shown in Figures10, 11, and 12, where the compatibility remains low, at less than 0.01%. We find that important sources of this incompatibility are the BNB-MB(νapp) and BNB-MB(ν-app) datasets. The BNB-MB(νapp) dataset has fairly small statistical and systematic uncertainties and therefore has a large impact on the fits and compatibility calculations. This is shown in Figure 13 where the MiniBooNE data agrees well with the appearance-only fit but disagrees with the overall global fit. Removing both the BNB-MB(νapp) and BNB-MB(ν-app) sets raises the compatibility to 3.5%, corresponding to an improvement of over two orders of magnitude. It has been known since the first MiniBooNE publication [18] that the BNB-MB(νapp) data was fairly consistent with no oscillations above 475 MeV; however, a significant low-energy excess was present below this energy. The energy dependence of the BNB-MB(νapp) excess does not fit very well with oscillation models extracted from fits to global datasets, unless very low Δmij2 with large mixing elements |Uei| and |Uμi| are involved in the fit. This may lead to the poor compatibility when included in appearance versus disappearance comparisons. Other possible explanations for this incompatibility include downward fluctuations of the BNB-MB(νapp) data in the higher-energy region or some other process contributing to the low-energy excess such as those suggested in [28, 74].TheΔm512 versus Δm412 and Δm612 versus Δm412 correlations from fits to all data in a (3+3) model.
(a)
(b)TheΔm512 versus Δm412 and Δm612 versus Δm412 correlations from fits to appearance data in a (3+3) model.
(a)
(b)TheΔm512 versus Δm412 and Δm612 versus Δm412 correlations from fits to disappearance data in a (3+3) model.
(a)
(b)A comparison of the BNB-MB(νapp) and BNB-MB(ν-app) excess data with the global best-fit oscillation signal predictions (solid colored lines) and with the appearance only best-fit predictions (dashed colored lines) for each of the models, (3+1), (3+2), and (3+3). The error bars on the excess correspond to statistical and unconstrained background systematic errors, added in quadrature.
(a)
(b)
(c)Statistical issues could be addressed with more MiniBooNE neutrino data that may become available over the next few years. In addition, the MicroBooNE experiment, which is expected to start running in 2014, will provide more information on the low-energy excess events and will answer the question of whether the excess is associated with outgoing electrons or photons [75].


## 5.4. Summary of Results

The sterile neutrino fits to global datasets show that a (3+1) model is inadequate; multiple Δmij2 values are needed along with CP-violating effects to explain the neutrino versus antineutrino differences. A (3+2) model improves significantly on the (3+1) results but still shows some tension in the neutrino versus antineutrino compatibility and cannot explain the appearance versus disappearance differences. The (3+3) model does not seem to further improve the fit and still has poor appearance versus disappearance compatibility. The BNB-MB(νapp) (and BNB-MB(ν-app)) dataset is a prime contributor to this incompatibility and additional experimental information in this region should be available soon. Figure 13 gives a comparison of the BNB-MB(νapp) and BNB-MB(ν-app) data with the global best-fit predictions and with the appearance-only best-fit predictions for each of the three models, (3+1), (3+2), and (3+3). The global fit prediction is significantly below the data at low energy, which contributes to the poor appearance versus disappearance compatibility.In summary, out of the three sterile neutrino oscillation hypotheses considered in the analysis, we find that the (3+2) and (3+3) models provide a better description than the (3+1) model, although the MiniBooNE appearance data continue to raise issues within the fits. As has been shown before, (3+1) scenarios provide a poor fit to the data, and should not be emphasized. We therefore recommend continued investigations of (3+2) and (3+3) scenarios.


## 6.  The Future

Establishing the existence of sterile neutrinos would have a major impact on particle physics. Motivated by this, there are a number of existing and planned experiments set to probe the parameter space indicative of one or more sterile neutrinos. Such experiments are necessary in order to confirm or refute the observed anomalies in theΔm2~1eV2 region. The new experiments are being designed to have improved sensitivity, with the goal of 5σ sensitivity and the ability to observe oscillatory behavior in L and/or E within single or between multiple detectors. In these experiments, the oscillation signal needs to be clearly separated from any backgrounds.Sterile neutrino oscillation models are based on oscillations associated with mixing between active and sterile states and demand the presence of both appearance and disappearance. It is therefore imperative that the future program explore both of these oscillation types. Establishing sterile neutrinos will require that both types of measurements are compatible with sterile neutrino oscillation models. Future experiments will search for evidence of sterile neutrino(s) using a variety of neutrino creation sources: (1) pion/muon DIF (e.g., [75–80]), (2) pion or kaon DAR (e.g., [81–86]), (3) unstable isotopes (e.g., [15, 87–90]), and (4) atmospheric (see [23]) and (5) nuclear reactors (e.g., [91, 92]). All of these experiments are under development and the sensitivities are likely to change. Therefore, rather than displaying sensitivity curves for each future program, we instead focus on the conceptual ideas behind the experiments. Unless otherwise mentioned, the experiments below will provide “significant” sensitivity to a large portion of the favored sterile neutrino parameter space through searches for neutrino disappearance and/or appearance.


### 6.1. The Importance of theL/E Signature from Multiple Experiments

Ultimately, in order to determine if there are zero, one, two, or three sterile states contributing to oscillations in SBL experiments, it will be necessary to observe the expectedL/E-dependent oscillation probabilities discussed in Section 2.3. Assuming that the SBL anomalies are confirmed, a consistent L/E dependence is the only signature which is distinct for oscillations and excludes other exotic explanations such as CPT violation [22], decays [93], and Lorentz violation [94]. The ideal experiment would reconstruct the oscillation wave as a function of L/E [95]. The combined information from many experiments, however, is more suitable for covering the widest possible range in L/E as well as providing valuable flavor and neutrino versus antineutrino information.The three models, (3+1), (3+2), and (3+3), have distinct signatures as a function of L/E. To illustrate this, we consider the case of a hypothetical experiment with 10% resolution in L/E, assuming the best-fit values presented in Table 3. In the case of (3+1), as shown in Figure 14, the disappearance (appearance) probabilities shown on the left (right), have maxima and minima that evolve monotonically to P=1/2sin2(2θ), the long-baseline limit discussed in Section 2.2. This can be contrasted with Figures 15 and 16, where the structure of the oscillation wave, in the approach to the long-baseline limit, is more “chaotic" due to the interference between the various mass splitting terms.The (3+1) oscillation probabilities for the global best-fit (“all” datasets) values in Table 3 with 10% resolution in L/E.
(a)
(b)The (3+2) oscillation probabilities for the global best-fit (“all” datasets) values in Table 3 with 10% resolution in L/E.
(a)
(b)The (3+3) oscillation probabilities for the global best-fit (“all” datasets) values in Table 3 with 10% resolution in L/E.
(a)
(b)In Figures14, 15, and 16, the two curves on the disappearance plots on the left refer to muon and electron flavor, respectively. As the theory is CPT-conserving, these disappearance curves should be identical for neutrinos and antineutrinos. The appearance curves on the right also show the importance of neutrino and antineutrino running, which can lead to very different L/E dependencies for the three models, and constrain CP-violating parameters.In summary, it seems very unlikely that anysingle future experiment will be able to differentiate between the sterile neutrino models. Multiple experiments looking at different oscillation channels and covering a wide range of L/E regions are required. Thus, the consideration of many independent experiments of relatively modest size, such as those listed in Table 4, is essential.Table 4
A summary of current and future sterile neutrino oscillation experiments.
Source
App/Dis
Channel
Experiment
Reactor
Dis
ν
-
e
→
ν
-
e
Nucifer, Stereo, SCRAMM, NIST, Neutrino4, DANSS
Radioactive
Dis
ν
e
→
ν
e,
ν-e→ν-e
Baksan, LENS, Borexino, SNO+, Ricochet, CeLAND, Daya Bay
Accelerator-based isotope
Dis
ν
-
e
→
ν
-
e
IsoDAR
Pion/Kaon DAR
App and Dis
ν
-
μ
→
ν
-
e,
νe→νe
OscSNS, DAEδALUS, KDAR
ν
μ
→
ν
e
Accelerator (Pion DIF)
App and Dis
ν
μ
→
ν
e,
ν-μ→ν-e
MINOS+, MicroBooNE, LAr1kton + MicroBooNE, CERN SPS
ν
μ
→
ν
μ,
ν-μ→ν-μ
Low-energyν-Factory
App and Dis
ν
e
→
ν
μ,
ν-e→ν-μ
νSTORM
ν
μ
→
ν
μ,
ν-μ→ν-μ


### 6.2. Future Experiments

A summary of future sterile neutrino experiments is provided in Table4.


#### 6.2.1. Pion Decay in Flight

Muon neutrinos (antineutrinos) from positive (negative) pion DIF can be used to search for (anti)neutrino disappearance and electron (anti)neutrino appearance in the sterile neutrino region of interest. Given the usual neutrino energies for these experiments (𝒪(1 GeV)), the baseline for such an experiment can be considered “short” (𝒪(100–1000 m)).The BNB at Fermilab will provide pion-induced neutrinos to the MicroBooNE LArTPC-based detector starting in 2014 [75]. MicroBooNE will probe the MiniBooNE low-energy anomaly [96] with a ~90 ton active volume about 100 meters closer to the neutrino source than MiniBooNE. Some coverage of the LSND allowed region in neutrino mode is also expected, along with LArTPC development and needed precision neutrino-argon cross-section measurements [97]. A design involving two LArTPC-based detectors in a near/far configuration, with MicroBooNE as the near detector, is also being considered for deployment in the BNB at Fermilab [78]. A similar two-detector configuration in the CERN-SPS neutrino beam has recently been proposed [79]. Two identical LArTPCs, in combination with magnetized spectrometers, would measure the mostly pion DIF-induced muon neutrino composition of the beam as a function of distance (300 m, 1600 m) to probe electron neutrino appearance in the sterile neutrino parameter space.Another BNB-based idea calls for a significant upgrade to the MiniBooNE experiment in which the current MiniBooNE detector becomes the 540 m baseline far detector in a two detector configuration and a MiniBooNE-like oil-based near detector is installed at a baseline of 200 m [77]. Such a configuration could significantly reduce the now largely irreducible systematics associated with MiniBooNE-far-only, which mainly come from neutral pion background events and flux uncertainty. In conjunction with MicroBooNE, “BooNE” could provide sensitivity to LSND-like electron (anti)neutrino appearance, muon (anti)neutrino disappearance, and the MiniBooNE low-energy anomaly.A low-energy 3-4 GeV/c muon storage ring could deliver a precisely known flux of electron neutrinos for a muon neutrino appearance search in the parameter space of interest for sterile neutrinos [80]. The MINOS-like detectors, envisioned at 20–50 m (near) and ~2000 m (far), would need to be magnetized in order to differentiate muon neutrino appearance from intrinsic muon antineutrinos created from the positive muon decay. Similar to most pion DIF beams, the muon storage ring could run in both neutrino mode and antineutrino mode. Such an experiment would also provide a technological demonstration of a muon storage ring with a “simple” neutrino factory [98].


#### 6.2.2. Pion or Kaon Decay at Rest

As discussed above, neutrinos from pion DAR and subsequent daughter muon DAR, with their well-known spectrum, provide a source for an oscillation search. Notably, LSND employed muon antineutrinos from the pion daughter’s muon DAR in establishing their 3.8σ excess consistent with electron antineutrino appearance.The 1 MW Spallation Neutron Source at Oak Ridge National Laboratory, a pion and muon DAR neutrino source, in combination with an LSND-style detector could directly probe the LSND excess with a factor of 100 lower steady state background and higher beam power [83, 84]. A 1 MW source at a large liquid scintillator detector is also under consideration [95]. Such an experiment could reconstruct appearance and disappearance oscillation waves across a ~50 m length of detector.If higher energy proton beams are used, then positive kaon DAR and the resulting monoenergetic (235.5 MeV) muon neutrino can also be used to search for sterile neutrinos through an electron neutrino appearance search with a LArTPC-based device [86]. An intense >3 GeV kinetic energy proton beam is required for such an experiment so as to produce an ample number of kaons per incoming proton.


#### 6.2.3. Unstable Isotopes

The disappearance of electron antineutrinos from radioactive isotopes is a direct probe of the reactor/gallium anomaly and an indirect probe of the LSND anomaly. As such neutrinos are in the ones-of-MeV range, the baseline for these experiments is generally on the order of tens of meters or so. Oscillation waves within a single detector can be observed if the neutrinos originate from a localized source, if the oscillation length is short enough, and if the detector has precise enough vertex resolution.The IsoDAR concept [90] calls for an intense 60 MeV proton source in combination with a kiloton-scale scintillation-based detector for sensitivity to the sterile neutrino. Such a source is being developed concurrently with the DAEδALUS experiment, nominally a search for nonzero δCP [99]. Cyclotron-produced 60 MeV protons impinge on a beryllium-based target, which mainly acts as a copious source of neutrons and is surrounded by an isotopically pure shell of 7Li. 8Li, created via neutron capture on 7Li inside the shell, decays to a 6.4 MeV mean energy electron antineutrino. Placing such an antineutrino source next to an existing detector such as KamLAND [3] could quickly provide discovery-level sensitivity in the reactor anomaly allowed region. Furthermore, IsoDAR has the ability to distinguish between one and multiple sterile neutrinos.Another unstable-isotope-based idea involves the deployment of a radioactive source inside an existing kiloton-scale detector [89] such as Borexino [100], KamLAND [3], or SNO+ [101]. Electron antineutrinos from a small-extent  ~2 PBq 144Ce or 106Ru beta source can be used to probe the sterile neutrino parameter space. For currently favored parameters associated with sterile neutrino(s), such antineutrinos are expected to disappear and reappear as a function of distance and energy inside the detector, much like the IsoDAR concept described above.


#### 6.2.4. Nuclear Reactor

A nuclear reactor can be used as a source for an electron antineutrino disappearance experiment with sensitivity to sterile neutrino(s). The Nucifer detector will likely be the first reactor-based detector to test the sterile neutrino hypothesis using antineutrino energy shape rather than just rate [91]. The experiment will take data in 2012/2013. The idea is to place a 1 m3 scale Gd-doped liquid scintillator device within a few tens of meters of a small-extent 70 MW research reactor in an attempt to observe antineutrino disappearance as a function of energy. The observation of an oscillation wave consistent with high Δm2 would be unambiguous evidence for the existence of at least one sterile neutrino. Cosmic ray interactions and their products represent the largest source of background for this class of experiment.One of the challenges of a reactor-based search is the need for a relatively small reactor size given the baseline required for maximal sensitivity toΔmij2~1eV2; a large neutrino source size relative to the neutrino baseline smears L and reduces Δmij2 resolution. A sterile search at a GW-scale power reactor is possible, however. The SCRAAM experiment (see [23]) calls for a Gd-doped liquid scintillator detector at the San Onofre Nuclear Generating Station.


#### 6.2.5. Neutral Current Based Experiments

All of the future experiments discussed previously involve either disappearance or appearance of neutrinos and antineutrinos detected via the charged current. However, a NC-based disappearance experiment provides unique sensitivity to the sterile neutrino. If neutrino disappearance was observed in a NC experiment, one would know that the active flavor neutrino(s) in question had oscillated into the noninteracting sterile flavor. Particularly, such an experiment would provide a measure of|Us4|, the sterile flavor composition of the fourth neutrino mass eigenstate, and definitively prove the existence of a sterile flavor neutrino, especially when considered in combination with CC-based experiments. A full understanding of the mixing angles associated with sterile neutrino(s) will require a NC-based experiment. The Ricochet concept [88] calls for oscillometry measurements using NC coherent neutrino-nucleus scattering detected via low temperature bolometers [81, 82, 88]. Both reactor and isotope decay sources are being considered for these measurements, utilizing the as-yet-undetected coherent neutrino-nucleus scattering process. The OscSNS experiment [83, 84] will also have the capability of looking for muon neutrino disappearance via the neutral current channel. Such a measurement would directly probe the sterile neutrino content of possible extra mass eigenstates.


## 6.1. The Importance of theL/E Signature from Multiple Experiments

Ultimately, in order to determine if there are zero, one, two, or three sterile states contributing to oscillations in SBL experiments, it will be necessary to observe the expectedL/E-dependent oscillation probabilities discussed in Section 2.3. Assuming that the SBL anomalies are confirmed, a consistent L/E dependence is the only signature which is distinct for oscillations and excludes other exotic explanations such as CPT violation [22], decays [93], and Lorentz violation [94]. The ideal experiment would reconstruct the oscillation wave as a function of L/E [95]. The combined information from many experiments, however, is more suitable for covering the widest possible range in L/E as well as providing valuable flavor and neutrino versus antineutrino information.The three models, (3+1), (3+2), and (3+3), have distinct signatures as a function of L/E. To illustrate this, we consider the case of a hypothetical experiment with 10% resolution in L/E, assuming the best-fit values presented in Table 3. In the case of (3+1), as shown in Figure 14, the disappearance (appearance) probabilities shown on the left (right), have maxima and minima that evolve monotonically to P=1/2sin2(2θ), the long-baseline limit discussed in Section 2.2. This can be contrasted with Figures 15 and 16, where the structure of the oscillation wave, in the approach to the long-baseline limit, is more “chaotic" due to the interference between the various mass splitting terms.The (3+1) oscillation probabilities for the global best-fit (“all” datasets) values in Table 3 with 10% resolution in L/E.
(a)
(b)The (3+2) oscillation probabilities for the global best-fit (“all” datasets) values in Table 3 with 10% resolution in L/E.
(a)
(b)The (3+3) oscillation probabilities for the global best-fit (“all” datasets) values in Table 3 with 10% resolution in L/E.
(a)
(b)In Figures14, 15, and 16, the two curves on the disappearance plots on the left refer to muon and electron flavor, respectively. As the theory is CPT-conserving, these disappearance curves should be identical for neutrinos and antineutrinos. The appearance curves on the right also show the importance of neutrino and antineutrino running, which can lead to very different L/E dependencies for the three models, and constrain CP-violating parameters.In summary, it seems very unlikely that anysingle future experiment will be able to differentiate between the sterile neutrino models. Multiple experiments looking at different oscillation channels and covering a wide range of L/E regions are required. Thus, the consideration of many independent experiments of relatively modest size, such as those listed in Table 4, is essential.Table 4
A summary of current and future sterile neutrino oscillation experiments.
Source
App/Dis
Channel
Experiment
Reactor
Dis
ν
-
e
→
ν
-
e
Nucifer, Stereo, SCRAMM, NIST, Neutrino4, DANSS
Radioactive
Dis
ν
e
→
ν
e,
ν-e→ν-e
Baksan, LENS, Borexino, SNO+, Ricochet, CeLAND, Daya Bay
Accelerator-based isotope
Dis
ν
-
e
→
ν
-
e
IsoDAR
Pion/Kaon DAR
App and Dis
ν
-
μ
→
ν
-
e,
νe→νe
OscSNS, DAEδALUS, KDAR
ν
μ
→
ν
e
Accelerator (Pion DIF)
App and Dis
ν
μ
→
ν
e,
ν-μ→ν-e
MINOS+, MicroBooNE, LAr1kton + MicroBooNE, CERN SPS
ν
μ
→
ν
μ,
ν-μ→ν-μ
Low-energyν-Factory
App and Dis
ν
e
→
ν
μ,
ν-e→ν-μ
νSTORM
ν
μ
→
ν
μ,
ν-μ→ν-μ


## 6.2. Future Experiments

A summary of future sterile neutrino experiments is provided in Table4.


### 6.2.1. Pion Decay in Flight

Muon neutrinos (antineutrinos) from positive (negative) pion DIF can be used to search for (anti)neutrino disappearance and electron (anti)neutrino appearance in the sterile neutrino region of interest. Given the usual neutrino energies for these experiments (𝒪(1 GeV)), the baseline for such an experiment can be considered “short” (𝒪(100–1000 m)).The BNB at Fermilab will provide pion-induced neutrinos to the MicroBooNE LArTPC-based detector starting in 2014 [75]. MicroBooNE will probe the MiniBooNE low-energy anomaly [96] with a ~90 ton active volume about 100 meters closer to the neutrino source than MiniBooNE. Some coverage of the LSND allowed region in neutrino mode is also expected, along with LArTPC development and needed precision neutrino-argon cross-section measurements [97]. A design involving two LArTPC-based detectors in a near/far configuration, with MicroBooNE as the near detector, is also being considered for deployment in the BNB at Fermilab [78]. A similar two-detector configuration in the CERN-SPS neutrino beam has recently been proposed [79]. Two identical LArTPCs, in combination with magnetized spectrometers, would measure the mostly pion DIF-induced muon neutrino composition of the beam as a function of distance (300 m, 1600 m) to probe electron neutrino appearance in the sterile neutrino parameter space.Another BNB-based idea calls for a significant upgrade to the MiniBooNE experiment in which the current MiniBooNE detector becomes the 540 m baseline far detector in a two detector configuration and a MiniBooNE-like oil-based near detector is installed at a baseline of 200 m [77]. Such a configuration could significantly reduce the now largely irreducible systematics associated with MiniBooNE-far-only, which mainly come from neutral pion background events and flux uncertainty. In conjunction with MicroBooNE, “BooNE” could provide sensitivity to LSND-like electron (anti)neutrino appearance, muon (anti)neutrino disappearance, and the MiniBooNE low-energy anomaly.A low-energy 3-4 GeV/c muon storage ring could deliver a precisely known flux of electron neutrinos for a muon neutrino appearance search in the parameter space of interest for sterile neutrinos [80]. The MINOS-like detectors, envisioned at 20–50 m (near) and ~2000 m (far), would need to be magnetized in order to differentiate muon neutrino appearance from intrinsic muon antineutrinos created from the positive muon decay. Similar to most pion DIF beams, the muon storage ring could run in both neutrino mode and antineutrino mode. Such an experiment would also provide a technological demonstration of a muon storage ring with a “simple” neutrino factory [98].


### 6.2.2. Pion or Kaon Decay at Rest

As discussed above, neutrinos from pion DAR and subsequent daughter muon DAR, with their well-known spectrum, provide a source for an oscillation search. Notably, LSND employed muon antineutrinos from the pion daughter’s muon DAR in establishing their 3.8σ excess consistent with electron antineutrino appearance.The 1 MW Spallation Neutron Source at Oak Ridge National Laboratory, a pion and muon DAR neutrino source, in combination with an LSND-style detector could directly probe the LSND excess with a factor of 100 lower steady state background and higher beam power [83, 84]. A 1 MW source at a large liquid scintillator detector is also under consideration [95]. Such an experiment could reconstruct appearance and disappearance oscillation waves across a ~50 m length of detector.If higher energy proton beams are used, then positive kaon DAR and the resulting monoenergetic (235.5 MeV) muon neutrino can also be used to search for sterile neutrinos through an electron neutrino appearance search with a LArTPC-based device [86]. An intense >3 GeV kinetic energy proton beam is required for such an experiment so as to produce an ample number of kaons per incoming proton.


### 6.2.3. Unstable Isotopes

The disappearance of electron antineutrinos from radioactive isotopes is a direct probe of the reactor/gallium anomaly and an indirect probe of the LSND anomaly. As such neutrinos are in the ones-of-MeV range, the baseline for these experiments is generally on the order of tens of meters or so. Oscillation waves within a single detector can be observed if the neutrinos originate from a localized source, if the oscillation length is short enough, and if the detector has precise enough vertex resolution.The IsoDAR concept [90] calls for an intense 60 MeV proton source in combination with a kiloton-scale scintillation-based detector for sensitivity to the sterile neutrino. Such a source is being developed concurrently with the DAEδALUS experiment, nominally a search for nonzero δCP [99]. Cyclotron-produced 60 MeV protons impinge on a beryllium-based target, which mainly acts as a copious source of neutrons and is surrounded by an isotopically pure shell of 7Li. 8Li, created via neutron capture on 7Li inside the shell, decays to a 6.4 MeV mean energy electron antineutrino. Placing such an antineutrino source next to an existing detector such as KamLAND [3] could quickly provide discovery-level sensitivity in the reactor anomaly allowed region. Furthermore, IsoDAR has the ability to distinguish between one and multiple sterile neutrinos.Another unstable-isotope-based idea involves the deployment of a radioactive source inside an existing kiloton-scale detector [89] such as Borexino [100], KamLAND [3], or SNO+ [101]. Electron antineutrinos from a small-extent  ~2 PBq 144Ce or 106Ru beta source can be used to probe the sterile neutrino parameter space. For currently favored parameters associated with sterile neutrino(s), such antineutrinos are expected to disappear and reappear as a function of distance and energy inside the detector, much like the IsoDAR concept described above.


### 6.2.4. Nuclear Reactor

A nuclear reactor can be used as a source for an electron antineutrino disappearance experiment with sensitivity to sterile neutrino(s). The Nucifer detector will likely be the first reactor-based detector to test the sterile neutrino hypothesis using antineutrino energy shape rather than just rate [91]. The experiment will take data in 2012/2013. The idea is to place a 1 m3 scale Gd-doped liquid scintillator device within a few tens of meters of a small-extent 70 MW research reactor in an attempt to observe antineutrino disappearance as a function of energy. The observation of an oscillation wave consistent with high Δm2 would be unambiguous evidence for the existence of at least one sterile neutrino. Cosmic ray interactions and their products represent the largest source of background for this class of experiment.One of the challenges of a reactor-based search is the need for a relatively small reactor size given the baseline required for maximal sensitivity toΔmij2~1eV2; a large neutrino source size relative to the neutrino baseline smears L and reduces Δmij2 resolution. A sterile search at a GW-scale power reactor is possible, however. The SCRAAM experiment (see [23]) calls for a Gd-doped liquid scintillator detector at the San Onofre Nuclear Generating Station.


### 6.2.5. Neutral Current Based Experiments

All of the future experiments discussed previously involve either disappearance or appearance of neutrinos and antineutrinos detected via the charged current. However, a NC-based disappearance experiment provides unique sensitivity to the sterile neutrino. If neutrino disappearance was observed in a NC experiment, one would know that the active flavor neutrino(s) in question had oscillated into the noninteracting sterile flavor. Particularly, such an experiment would provide a measure of|Us4|, the sterile flavor composition of the fourth neutrino mass eigenstate, and definitively prove the existence of a sterile flavor neutrino, especially when considered in combination with CC-based experiments. A full understanding of the mixing angles associated with sterile neutrino(s) will require a NC-based experiment. The Ricochet concept [88] calls for oscillometry measurements using NC coherent neutrino-nucleus scattering detected via low temperature bolometers [81, 82, 88]. Both reactor and isotope decay sources are being considered for these measurements, utilizing the as-yet-undetected coherent neutrino-nucleus scattering process. The OscSNS experiment [83, 84] will also have the capability of looking for muon neutrino disappearance via the neutral current channel. Such a measurement would directly probe the sterile neutrino content of possible extra mass eigenstates.


## 6.2.1. Pion Decay in Flight

Muon neutrinos (antineutrinos) from positive (negative) pion DIF can be used to search for (anti)neutrino disappearance and electron (anti)neutrino appearance in the sterile neutrino region of interest. Given the usual neutrino energies for these experiments (𝒪(1 GeV)), the baseline for such an experiment can be considered “short” (𝒪(100–1000 m)).The BNB at Fermilab will provide pion-induced neutrinos to the MicroBooNE LArTPC-based detector starting in 2014 [75]. MicroBooNE will probe the MiniBooNE low-energy anomaly [96] with a ~90 ton active volume about 100 meters closer to the neutrino source than MiniBooNE. Some coverage of the LSND allowed region in neutrino mode is also expected, along with LArTPC development and needed precision neutrino-argon cross-section measurements [97]. A design involving two LArTPC-based detectors in a near/far configuration, with MicroBooNE as the near detector, is also being considered for deployment in the BNB at Fermilab [78]. A similar two-detector configuration in the CERN-SPS neutrino beam has recently been proposed [79]. Two identical LArTPCs, in combination with magnetized spectrometers, would measure the mostly pion DIF-induced muon neutrino composition of the beam as a function of distance (300 m, 1600 m) to probe electron neutrino appearance in the sterile neutrino parameter space.Another BNB-based idea calls for a significant upgrade to the MiniBooNE experiment in which the current MiniBooNE detector becomes the 540 m baseline far detector in a two detector configuration and a MiniBooNE-like oil-based near detector is installed at a baseline of 200 m [77]. Such a configuration could significantly reduce the now largely irreducible systematics associated with MiniBooNE-far-only, which mainly come from neutral pion background events and flux uncertainty. In conjunction with MicroBooNE, “BooNE” could provide sensitivity to LSND-like electron (anti)neutrino appearance, muon (anti)neutrino disappearance, and the MiniBooNE low-energy anomaly.A low-energy 3-4 GeV/c muon storage ring could deliver a precisely known flux of electron neutrinos for a muon neutrino appearance search in the parameter space of interest for sterile neutrinos [80]. The MINOS-like detectors, envisioned at 20–50 m (near) and ~2000 m (far), would need to be magnetized in order to differentiate muon neutrino appearance from intrinsic muon antineutrinos created from the positive muon decay. Similar to most pion DIF beams, the muon storage ring could run in both neutrino mode and antineutrino mode. Such an experiment would also provide a technological demonstration of a muon storage ring with a “simple” neutrino factory [98].


## 6.2.2. Pion or Kaon Decay at Rest

As discussed above, neutrinos from pion DAR and subsequent daughter muon DAR, with their well-known spectrum, provide a source for an oscillation search. Notably, LSND employed muon antineutrinos from the pion daughter’s muon DAR in establishing their 3.8σ excess consistent with electron antineutrino appearance.The 1 MW Spallation Neutron Source at Oak Ridge National Laboratory, a pion and muon DAR neutrino source, in combination with an LSND-style detector could directly probe the LSND excess with a factor of 100 lower steady state background and higher beam power [83, 84]. A 1 MW source at a large liquid scintillator detector is also under consideration [95]. Such an experiment could reconstruct appearance and disappearance oscillation waves across a ~50 m length of detector.If higher energy proton beams are used, then positive kaon DAR and the resulting monoenergetic (235.5 MeV) muon neutrino can also be used to search for sterile neutrinos through an electron neutrino appearance search with a LArTPC-based device [86]. An intense >3 GeV kinetic energy proton beam is required for such an experiment so as to produce an ample number of kaons per incoming proton.


## 6.2.3. Unstable Isotopes

The disappearance of electron antineutrinos from radioactive isotopes is a direct probe of the reactor/gallium anomaly and an indirect probe of the LSND anomaly. As such neutrinos are in the ones-of-MeV range, the baseline for these experiments is generally on the order of tens of meters or so. Oscillation waves within a single detector can be observed if the neutrinos originate from a localized source, if the oscillation length is short enough, and if the detector has precise enough vertex resolution.The IsoDAR concept [90] calls for an intense 60 MeV proton source in combination with a kiloton-scale scintillation-based detector for sensitivity to the sterile neutrino. Such a source is being developed concurrently with the DAEδALUS experiment, nominally a search for nonzero δCP [99]. Cyclotron-produced 60 MeV protons impinge on a beryllium-based target, which mainly acts as a copious source of neutrons and is surrounded by an isotopically pure shell of 7Li. 8Li, created via neutron capture on 7Li inside the shell, decays to a 6.4 MeV mean energy electron antineutrino. Placing such an antineutrino source next to an existing detector such as KamLAND [3] could quickly provide discovery-level sensitivity in the reactor anomaly allowed region. Furthermore, IsoDAR has the ability to distinguish between one and multiple sterile neutrinos.Another unstable-isotope-based idea involves the deployment of a radioactive source inside an existing kiloton-scale detector [89] such as Borexino [100], KamLAND [3], or SNO+ [101]. Electron antineutrinos from a small-extent  ~2 PBq 144Ce or 106Ru beta source can be used to probe the sterile neutrino parameter space. For currently favored parameters associated with sterile neutrino(s), such antineutrinos are expected to disappear and reappear as a function of distance and energy inside the detector, much like the IsoDAR concept described above.


## 6.2.4. Nuclear Reactor

A nuclear reactor can be used as a source for an electron antineutrino disappearance experiment with sensitivity to sterile neutrino(s). The Nucifer detector will likely be the first reactor-based detector to test the sterile neutrino hypothesis using antineutrino energy shape rather than just rate [91]. The experiment will take data in 2012/2013. The idea is to place a 1 m3 scale Gd-doped liquid scintillator device within a few tens of meters of a small-extent 70 MW research reactor in an attempt to observe antineutrino disappearance as a function of energy. The observation of an oscillation wave consistent with high Δm2 would be unambiguous evidence for the existence of at least one sterile neutrino. Cosmic ray interactions and their products represent the largest source of background for this class of experiment.One of the challenges of a reactor-based search is the need for a relatively small reactor size given the baseline required for maximal sensitivity toΔmij2~1eV2; a large neutrino source size relative to the neutrino baseline smears L and reduces Δmij2 resolution. A sterile search at a GW-scale power reactor is possible, however. The SCRAAM experiment (see [23]) calls for a Gd-doped liquid scintillator detector at the San Onofre Nuclear Generating Station.


## 6.2.5. Neutral Current Based Experiments

All of the future experiments discussed previously involve either disappearance or appearance of neutrinos and antineutrinos detected via the charged current. However, a NC-based disappearance experiment provides unique sensitivity to the sterile neutrino. If neutrino disappearance was observed in a NC experiment, one would know that the active flavor neutrino(s) in question had oscillated into the noninteracting sterile flavor. Particularly, such an experiment would provide a measure of|Us4|, the sterile flavor composition of the fourth neutrino mass eigenstate, and definitively prove the existence of a sterile flavor neutrino, especially when considered in combination with CC-based experiments. A full understanding of the mixing angles associated with sterile neutrino(s) will require a NC-based experiment. The Ricochet concept [88] calls for oscillometry measurements using NC coherent neutrino-nucleus scattering detected via low temperature bolometers [81, 82, 88]. Both reactor and isotope decay sources are being considered for these measurements, utilizing the as-yet-undetected coherent neutrino-nucleus scattering process. The OscSNS experiment [83, 84] will also have the capability of looking for muon neutrino disappearance via the neutral current channel. Such a measurement would directly probe the sterile neutrino content of possible extra mass eigenstates.


## 7. Conclusions

This paper has presented results of SBL experiments discussed within the context of oscillations involving sterile neutrinos. Fits to (3+1), (3+2), and (3+3) models have been presented. We have examined whether the (3+3) model addresses tensions observed with (3+1) and (3+2) fits.Several issues arise when comparing datasets in (3+1) and (3+2) models. In a (3+1) model, the compatibility of the neutrino versus antineutrino datasets is poor (0.14%), and the compatibility among all datasets is only 0.043%. In a (3+2) model, there is a striking disagreement between appearance and disappearance datasets, with a compatibility of 0.0082%.A 3 + 3 (3 + 2) model fit has aχ2-probability for the best-fit of 67% (69%), compared to 2.1% for the no oscillation scenario. Though these values are on the order of the χ2-probabilities found for the 3 + 1 model, the 3 + 3 (3 + 2) fits resolve the incompatibility issues seen in the 3 + 1 model, with the exception of the MiniBooNE appearance datasets. Therefore, we argue that the 3 + 2 and 3 + 3 fits should be the main focus of sterile neutrino phenomenological studies in the future.While the indications of sterile neutrino oscillations have historically been associated with only appearance-based SBL experiments, the recently realized suppression in observedν-e in disappearance reactor experiments provides further motivation for these models. As we have shown, one can consistently fit most results under the (3+3) hypothesis with improved compatibility. However, the need for additional information from both appearance and disappearance experiments provides strong motivation for pursuing the future experiments discussed in this paper.



---

*Source: 163897-2013-03-06.xml*