Advances in High Energy Physics

Volume 2013 (2013), Article ID 439532, 14 pages

http://dx.doi.org/10.1155/2013/439532

Review Article

Evidence and Search for Sterile Neutrinos at Accelerators

W. C. Louis

Los Alamos National Laboratory, Los Alamos, NM 87545, USA

Received 11 June 2013; Accepted 13 August 2013

Academic Editor: Leslie Camilleri

Copyright © 2013 W. C. Louis. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The LSND short-baseline neutrino experiment has published evidence for antineutrino oscillations at a mass scale of ~1 eV². The MiniBooNE experiment, designed to test this evidence for oscillations at an order of magnitude higher neutrino energy and distance, observes excesses of events in both neutrino mode and antineutrino mode. While the MiniBooNE neutrino excess has a neutrino energy spectrum that is softer than expected from LSND, the MiniBooNE antineutrino excess is consistent with neutrino oscillations and with the LSND oscillation signal. When combined with oscillation measurements at the solar and atmospheric mass scales, assuming that the LSND and MiniBooNE signals are due to neutrino oscillations, these experiments imply the existence of more than three neutrino mass states and, therefore, one or more sterile neutrinos. Such sterile neutrinos, if proven to exist, would have a big impact on particle physics, nuclear physics, and astrophysics and would contribute to the dark matter of the universe. Future experiments under construction or proposed at Fermilab, ORNL, CERN, and in Japan will provide a definitive test of short-baseline neutrino oscillations and will have the capability of proving the existence of sterile neutrinos.

1. Introduction

Neutrino oscillations have been clearly observed at the solar mass scale of   and the atmospheric mass scale of [1]. However, evidence for neutrino oscillations at the   mass scale has also been given by the LSND [2–6] and MiniBooNE [7–9] experiments. If the signals reported by LSND and MiniBooNE are indeed due to neutrino oscillations, then there is a problem with the three-neutrino paradigm, because it is not possible to explain oscillations at the three different mass scales with only three types of neutrinos. In order to solve this problem, additional “sterile” neutrinos have been proposed [10–14]. Other explanations include, for example, Lorentz violation [15, 16] and sterile neutrino decay [17, 18]. These additional neutrinos would need to be sterile to the weak interaction due to measurements of the width of the boson, which determine that there are three and only three active neutrinos [19]. The evidence for neutrino oscillations from LSND and MiniBooNE will be discussed in the following sections, followed by a discussion of global fits to the world data and future experiments, which will have the capability of proving whether short-baseline neutrino oscillations at the ~1  mass scale and light, sterile neutrinos exist in nature.

2. The LSND Experiment

The LSND experiment [20] was designed to search for oscillations with high sensitivity and to measure cross-sections. A photograph of the inside of the detector tank is shown in Figure 1. The main characteristics of the LSND experiment are given in Table 1. LSND had the advantage of a very high proton intensity, a large detector mass, and good particle identification. LSND made use of a high-intensity, 798 MeV proton beam that interacted in an absorber to produce a large number of pions.

Table 1: The main characteristics of the LSND and MiniBooNE experiments. Also the estimated number of neutrino background events and signal events for 100% transmutation are shown. The requirement is assumed for the LSND events.

Property
LSND
MiniBooNE
MiniBooNE

Proton energy
798 MeV
8000 MeV
8000 MeV

Proton intensity
1000 A
4 A
4 A

Proton beam power
798 kW
32 kW
32 kW

Protons on target
28,896 C
104 C
181 C

Duty factor

Total mass
167 t
806 t
806 t

Neutrino distance
30 m
541 m
541 m

100% events
12,987
89,615
38,462

background events
16.9
797.7
398.2

Figure 1: A photograph of the inside of the LSND detector tank.

2.1. Neutrino Beam

LSND made use of the LAMPF accelerator, which was an intense source of low energy neutrinos produced with a proton current of 1 mA at 798 MeV kinetic energy. For the 1993–1995 running period, the production target consisted of a 30 cm long water target (20 cm in 1993) followed by a water-cooled Cu beam dump, while for the 1996–1998 running period, the production target was reconfigured with the water target replaced by a close-packed, high- target. The resulting decay-at-rest (DAR) neutrino fluxes are well understood because almost all detectable neutrinos arise from or decay; and that stop are readily captured in the Fe of the shielding and Cu of the beam stop [21, 22]. The production of kaons or heavier mesons is negligible at these proton energies. The flux is calculated to be only ~8 × 10⁻⁴ as large as the flux in the  MeV energy range, so that the observation of a event rate significantly above the calculated background would be evidence for oscillations. Figure 2 shows the neutrino energy spectra from and DAR.

(a)

(b)

Figure 2: The neutrino energy spectra from and DAR.

2.2. Detector

The LSND detector [20] consisted of an approximately cylindrical tank 8.3 m long by 5.7 m in diameter. A schematic drawing of the detector is shown in Figure 3. The center of the detector was 30 m away from the neutrino source. On the inside surface of the tank, 1220 8-inch Hamamatsu PMTs covered 25% of the area with photocathode. The tank was filled with 167 t of liquid scintillator consisting of mineral oil and 0.031 g/L of b-PBD. This low scintillator concentration allows the detection of both Cherenkov light and scintillation light and yields an attenuation length of more than 20 m for wavelengths greater than 400 nm [23]. A typical 45 MeV electron created in the detector produced a total of ~1500 photoelectrons, of which ~280 photoelectrons were in the Cherenkov cone. PMT time and pulse-height signals were used to reconstruct particle tracks with an average RMS position resolution of ~14cm, an angular resolution of ~12°, and an energy resolution of ~7% at the Michel endpoint of 52.8 MeV. The Cherenkov cone for relativistic particles and the time distribution of the light, which is broader for nonrelativistic particles [20], gave excellent separation between electrons and particles below Cherenkov threshold. Identification of neutrons was accomplished through the detection of the  MeV from neutron capture on a free proton. The veto shield enclosed the detector on all sides except the bottom. Additional counters were placed below the veto shield after the 1993 run to reduce cosmic-ray background entering through the bottom support structure. The main veto shield [24] consisted of a 15 cm layer of liquid scintillator in an external tank and 15 cm of lead shot in an internal tank. This combination of active and passive shielding tagged cosmic-ray muons that stopped in the lead shot. A veto inefficiency of <10⁻⁵ was achieved for incident charged particles.

Figure 3: A schematic drawing of the LSND detector.

2.3. Oscillation Results

The primary oscillation search in LSND is for oscillations, where the arise from DAR in the beam stop and the are identified through the reaction . This reaction allows a two-fold signature of a positron with a 52.8 MeV endpoint and a correlated 2.2 MeV from neutron capture on a free proton. There are only two significant neutrino backgrounds with a positron/electron and a correlated neutron. The first background is from DAR in the beam stop followed by scattering in the detector. This background is highly suppressed due to the requirements that a be produced, the decays in flight, and the decays at rest prior to capture. The second background is from decay in flight (DIF) in the beam stop followed by scattering in the detector. Additional contributions are from and scattering. This background will mimic the oscillation reaction if the is sufficiently low in energy that it is below the threshold of 18 hit PMTs, corresponding to MeV. Table 2 shows the estimated number of events in the MeV energy range satisfying the electron selection criteria for 100% transmutation and for the two beam-related backgrounds with neutrons.

Table 2: The LSND estimated number of events in the  MeV energy range due to 100% transmutation and the two beam-related backgrounds with neutrons, DAR in the beam stop followed by scattering in the detector and DIF in the beam stop followed by scattering. The events must satisfy the electron selection criteria, but no correlated requirement is imposed.

Neutrino source
Reaction
Number of events

DAR
100%  

DAR

DIF

Table 3 shows the LSND statistics for events that satisfy the selection criteria for the primary oscillation search. An excess of events is observed over that expected from beam-off and neutrino background that is consistent with neutrino oscillations [2–6]. A fit to the distribution, as shown in Figure 4, gives ( DOF), which leads to a beam on-off excess of events with a correlated neutron. is the likelihood that the is a correlated from neutron capture divided by the likelihood that the is an accidental , while is the fraction of electron candidate events with a correlated . Subtracting the neutrino background from DAR followed by scattering ( events) and DIF followed by scattering ( events) (this background also includes contributions from and ) leads to a total excess of events. This excess corresponds to an oscillation probability of ()%, where the first error is statistical and the second error is the systematic error arising from uncertainties in the backgrounds, neutrino flux (7%), efficiency (7%), and efficiency (7%).

Table 3: Numbers of LSND beam-on events that satisfy the selection criteria for the primary oscillation search with , , and . Also the beam-off background, the estimated neutrino background, and the excess of events that is consistent with neutrino oscillations are shown.

Selection
Beam-on
Beam-off
background
Event excess

events
background

205

86

27

Figure 4: The LSND distribution for events that satisfy the selection criteria for the primary oscillation search.

A clean sample of oscillation candidate events can be obtained by requiring , where, as shown in Table 3, the beam on-off excess is events, while the estimated neutrino background is only events. Figure 5 displays the energy distribution of events with . The shaded regions show the combination of neutrino background plus neutrino oscillations at low . The data agree well with the oscillation hypothesis. Finally, Figure 6 shows the distribution for events with and  MeV, where is the distance travelled by the neutrino in meters and is the neutrino energy in MeV determined from the measured positron energy and angle with respect to the neutrino beam. The data agree well with the expectation from neutrino background plus neutrino oscillations at low ( D.O.F.) or high ( D.O.F.).

Figure 5: The LSND electron energy distribution for events with . The shaded region shows the expected distribution from a combination of neutrino background plus neutrino oscillations at low .

Figure 6: The LSND distribution for events with and MeV, where is the distance travelled by the neutrino in meters and is the neutrino energy in MeV. The data agree well with the expectation from neutrino background and neutrino oscillations at low .

The likelihood () fitter is applied to beam-on events in the final oscillation sample and calculates a likelihood in the plane in order to extract the favored oscillation parameters. The product in the plane is formed over the individual beam-on events that pass the oscillation cuts. This three-dimensional contour is sliced to arrive finally at the LSND allowed oscillation region. The beam-related backgrounds are determined from Monte Carlo (MC) event samples for each individual background contribution. The MC contains the trigger simulation and generally very well reproduces the tank response to all particles of interest. Agreement between the data and MC is excellent.

The oscillation parameter fit for the entire data sample, MeV, is shown in Figure 7. The fit includes both and oscillations (although the latter contribution is very small), as well as all known neutrino backgrounds. The inner and outer regions correspond to 90% and 99% CL allowed regions, while the curves are 90% CL limits from the Bugey reactor experiment [25] and the KARMEN experiment at ISIS [26]. Note that the Bugey limit has been relaxed recently due to the reactor neutrino anomaly [27]. The most favored allowed region is the band from to , although a region around is also possible.

Figure 7: The oscillation parameter fit for the entire LSND data sample,  MeV. The inner and outer regions correspond to 90% and 99% CL allowed regions, while the curves are 90% CL limits from the Bugey reactor experiment and the KARMEN experiment at ISIS.

3. The MiniBooNE Experiment

The MiniBooNE experiment was designed to test the LSND neutrino oscillation signal with a neutrino energy and distance that were an order of magnitude higher than LSND but with an ratio that was approximately the same as for LSND. Due to the higher neutrino energies, the backgrounds and systematic errors are completely different between the two experiments.

3.1. Neutrino Beam

A schematic drawing of the MiniBooNE experiment at FNAL is shown in Figure 8. The experiment is fed by 8 GeV kinetic energy protons from the Booster that interact in a 71 cm long Be target located at the upstream end of a magnetic focusing horn. The horn pulses with a current of 174 kA and, depending on the polarity, either focuses and and defocuses and to form a pure neutrino beam or focuses and and defocuses and to form a somewhat pure antineutrino beam. The produced pions and kaons decay in a 50 m long pipe, and a fraction of the neutrinos and antineutrinos [28] interact in the MiniBooNE detector, which is located 541 m downstream of the Be target. For the MiniBooNE results presented here, a total of  POT were collected in neutrino mode and  POT were collected in antineutrino mode.

Figure 8: A schematic drawing of the MiniBooNE experiment.

3.2. Detector

The MiniBooNE detector [29] consists of a 12.2 m diameter spherical tank filled with approximately 806 tons of mineral oil (CH₂). A schematic drawing of the MiniBooNE detector is shown in Figure 9, while a photograph of the inside of the tank is shown in Figure 10. There are a total of 1280 8-inch detector phototubes (covering 10% of the surface area) and 240 veto phototubes. The fiducial volume has a 5 m radius and corresponds to approximately 450 tons. A comparison between the LSND and MiniBooNE experiments is given in Table 1.

Figure 9: A schematic drawing of the MiniBooNE detector.

Figure 10: A photograph of the inside of the MiniBooNE detector.

3.3. Oscillation Results

Figure 11 shows the estimated neutrino fluxes for neutrino mode and antineutrino mode, respectively. The fluxes are fairly similar (the intrinsic electron-neutrino background is approximately 0.6% for both modes of running), although the wrong-sign contribution to the flux in antineutrino mode (~18%) is much larger than in neutrino mode (~6%). The average plus energies are 0.96 GeV in neutrino mode and 0.77 GeV in antineutrino mode, while the average plus energies are 0.79 GeV in neutrino mode and 0.66 GeV in antineutrino mode. The estimated backgrounds in the two modes are very similar, especially at low energy.

(a)

(b)

Figure 11: The estimated MiniBooNE neutrino fluxes for neutrino mode (a) and antineutrino mode (b).

Figure 12 shows the reconstructed neutrino and antineutrino energy distributions, , for candidate and data events (points with error bars) compared to the MC simulation (histogram) [7–9], while Figure 13 shows the event excesses in neutrino and antineutrino mode as a function of . assumes that the events are due to charged-current quasi-elastic scattering. The magnitude and energy shape of the excess in antineutrino mode are very similar to what is expected from neutrino oscillations based on the LSND signal, as shown by the pink and green LSND reference histograms. In neutrino mode, the magnitude of the excess is consistent with LSND; however, the energy shape is softer than the pink and green LSND reference histograms.

(a)

(b)

Figure 12: The MiniBooNE neutrino mode (a) and antineutrino mode (b) neutrino energy distributions for charged-current quasi-elastic data (points with statistical errors) and background (histogram with systematic errors).

(a)

(b)

Figure 13: The MiniBooNE neutrino mode (a) and antineutrino mode (b) event excesses as a function of neutrino energy. Also the expectations from the best two-neutrino oscillation fits with neutrino energy from 200 to 3000 MeV and from two reference values in the LSND allowed region are shown. All known systematic errors are included in the systematic error estimate.

Table 4 shows the number of data, background, and excess events for neutrino mode, antineutrino mode, and combined in the neutrino oscillation energy range from 200 to 1250 MeV. The uncertainties include both statistical and constrained systematic errors, where the charged-current data sample is used to normalize the backgrounds. All known systematic errors are included in the systematic error estimate. In neutrino mode, an excess of events, corresponding to , is observed, while in antineutrino mode, the observed excess is events, corresponding to . Combining the data in neutrino and antineutrino modes, the total excess is events, corresponding to .

Table 4: The number of data, fitted (constrained) backgrounds, and excess events in the and analyses for neutrino mode, antineutrino mode, and combined in the neutrino oscillation energy range from 200 to 1250 MeV. The uncertainties include both statistical and constrained systematic errors. All known systematic errors are included in the systematic error estimate.

Mode
Data
Background
Excess

Neutrino mode
952

Antineutrino mode
478

Combined
1430

Figure 14 shows the MiniBooNE oscillation allowed regions [7–9] from fits to the antineutrino data (a) and neutrino data (b). In antineutrino mode, the probability of the best fit is good (66%, corresponding to a ), and the MiniBooNE allowed region overlaps well with the LSND allowed region for , while the KARMEN experiment [26] rules out most of the LSND allowed region for . Therefore, the MiniBooNE antineutrino data are consistent with simple two-neutrino oscillations and with the LSND oscillation signal. In neutrino mode, the probability of the best fit is not as good (6.1%, corresponding to a ), and the MiniBooNE allowed region overlaps a smaller fraction of the LSND allowed region at low , which is already ruled out by the ICARUS experiment [30]. Therefore, in neutrino mode, the data are marginally compatible with a simple two-neutrino oscillation formalism. Indeed, the initial MiniBooNE oscillation publication [34] showed that there is hardly any excess of events above 475 MeV in neutrino mode, ruling out simple two-neutrino oscillations as an explanation of the LSND signal. However, expanded models with several sterile neutrinos can reduce the incompatibility by allowing for CP violating effects between neutrino and antineutrino oscillations [10, 11].

(a)

(b)

Figure 14: MiniBooNE allowed regions in antineutrino mode (a) and neutrino mode (b) for events with neutrino energy greater than 200 MeV within a two-neutrino oscillation model. Also the ICARUS [30] and KARMEN [26] appearance limits for neutrinos and antineutrinos are shown, respectively. The shaded areas show the 90% and 99% CL LSND allowed regions. The black stars show the MiniBooNE best fit points, while the circles show two reference values in the LSND allowed region.

Figure 15 shows the MiniBooNE oscillation allowed regions from a combined fit to the antineutrino plus neutrino data samples [35]. As in neutrino mode, the probability of the best fit is not very good (6.7%). However, as discussed above, expanded models with several sterile neutrinos can improve the fit probability by allowing for CP violating effects between neutrino and antineutrino oscillations [10, 11].

Figure 15: MiniBooNE allowed regions in combined neutrino and antineutrino mode for events with neutrino energy from 200 to 3000 MeV within a two-neutrino and oscillation model. Also shown is the limit from the KARMEN experiment [26]. The shaded areas show the 90% and 99% CL LSND allowed regions. The black star shows the best fit point.

4. Global Neutrino Oscillation Analyses

Several global oscillation analyses have been performed on the world neutrino data that include both long-baseline and short-baseline oscillation experiments, including both and appearance data and and disappearance data [12]. In general, the analyses assume a model with 3 light, mostly active neutrinos, heavier, mostly sterile neutrinos, and neutrino mass eigenstates. Precision measurements of the width have determined that the number of active neutrinos is approximately 3 ( [19]), so that additional neutrinos beyond the three active neutrinos would need to be sterile to the weak interaction. Note that neutrino oscillations depend on the superposition of neutrino mass eigenstates and do not depend on whether neutrinos are active or sterile, so that the neutrino mass eigenstates would result in independent values. Note also that, for , there would exist CP violation parameters for both short-baseline and long-baseline experiments. Finally, note that, for models, neutrinos in general would be a superposition of both active and sterile components.

Table 5 shows the best fit values from recent global fits to the world neutrino data from [10, 11]. (Additional global fits can be found in [12].) Reference [10] includes data from the LSND, MiniBooNE, KARMEN, NOMAD, CCFR, CDHS, MINOS, Bugey, SAGE, and GALLEX experiments, while [11] also includes data from the E776 and ICARUS experiments. For both global fits, the best fit value of is ~0.9  and the goodness of fit probabilities (GOF) are reasonable at 55% and 19%, respectively. However, the parameter goodness of fit probabilities (PGOF) are very low at and , respectively. The PGOF [36] is based on the difference between the overall best fit chi square and the sum of the best fit chi squares for the appearance-only and disappearance-only experiments. Although there is tension between the appearance and disappearance experiments, we believe that fake data studies would need to be performed in order to calculate the PGOF true probabilities [37].

Table 5: The best fit values from global fits to the world neutrino data from [10, 11]. Also the goodness of fit probability (GOF) and the parameter goodness of fit probability (PGOF) [36] are shown.

Ref.

GOF
PGOF

[10]
0.92
0.15
0.17
55%
0.043%

[11]
0.93
0.15
0.17
19%
0.012%

Table 6 shows the best fit values from recent global fits to the world neutrino data from [10, 11]. (Additional global fits can be found in [12].) The GOF probabilities are improved from the fits, although the PGOF probabilities are worse. As mentioned above, we believe that fake data studies would need to be performed in order to calculate the PGOF true probabilities [37]. Figures 16 and 17 show the best fit versus values for the two global fits. Figure 17 also shows the best fit, where the two mostly sterile neutrinos have masses less than and greater than the masses of the mostly active neutrinos.

Table 6: The best fit values from global fits to the world neutrino data from [10, 11]. Also the goodness of fit probability (GOF) and the parameter goodness of fit probability (PGOF) are shown.

Ref.

GOF
PGOF

[10]
0.92
17
0.15
0.17
0.069
0.16

69%
0.0082%

[11]
0.47
0.87
0.13
0.15
0.14
0.13

23%
0.0034%

Figure 16: The best fit versus values from [10].

Figure 17: The and best fit versus values from [11].

5. Future Short-Baseline Accelerator Neutrino Experiments

Future accelerator neutrino experiments will have the capability of testing the current evidence for sterile neutrinos. At Fermilab, MINOS+ and MicroBooNE provide a well-coordinated program of experiments that should help determine in the near future whether or not sterile neutrinos can explain the existing short-baseline neutrino anomalies. MINOS+ involves running the MINOS detectors with the NuMI beam in medium energy mode for the NoVA experiment, which provides a much higher flux of neutrinos above a few GeV. High disappearance into sterile neutrinos can be searched for by comparing the event rates in the near and far detectors at higher energies, where the effect of neutrino oscillations at the atmospheric scale is small [38]. MINOS+ should begin taking data later in 2013. The MicroBooNE experiment is under construction in the Booster Neutrino Beam (BNB) just upstream of MiniBooNE at a distance of 0.47 km from the neutrino source. MicroBooNE [39] consists of a 170-ton (86-ton active volume) liquid argon TPC, that will be able to determine whether the excess of events observed by MiniBooNE is due to electron events, as expected for oscillations, or to photon events, which would indicate some other process. Figure 18 shows the MicroBooNE sensitivity to oscillations. MicroBooNE should begin taking data in 2014.

Figure 18: The MicroBooNE sensitivity to oscillations. Also the allowed regions from LSND and from the MiniBooNE combined fit are shown.

NuSTORM [31] is a far future proposal that involves building a muon storage ring with central momentum of 3.8 GeV/c, and it can be built at either Fermilab or CERN. NuSTORM produces fluxes of neutrinos with well-understood energy spectra: and for a stored beam and and for a stored beam. Two magnetized iron neutrino detectors will be built downstream of the storage ring at distances of ~30 m and ~1500 m in order to search for appearance and and disappearance. The allowed LSND and MiniBooNE oscillation regions can be covered at the level, as shown in Figure 19.

Figure 19: NuSTORM can cover the allowed LSND and MiniBooNE oscillation regions at the level (from [31]).

ICARUS-NESSIE [32] is another far future proposal for an experiment at CERN to test the LSND and MiniBooNE oscillation signals. A new neutrino beamline would be built from the SPS with both near and far liquid argon TPC detectors (ICARUS) and muon spectrometers (NESSIE). The near detector location would be at a distance of 300 m from the neutrino source and would contain a 150-ton (T150) liquid argon TPC followed by a muon spectrometer with a dipole magnetic field and iron slabs. The far detector location would be at a distance of 1600 m from the neutrino source and would contain a 600-ton LAr TPC (T600) followed by a muon spectrometer. ICARUS-NESSIE will be able to search for both disappearance and appearance. Figure 20 shows the expected sensitivities for oscillations (a) and oscillations (b) after one year and two years of data taking, respectively. ICARUS-NESSIE will fully cover the LSND and MiniBooNE oscillation allowed regions.

(a)

(b)

Figure 20: The ICARUS-NESSIE expected sensitivities for oscillations (a) and oscillations (b) after one year and two years of data taking, respectively (from [32]).

The OscSNS and IsoDAR proposals offer the unique opportunity to observe actual oscillations in the detector (i.e., an oscillation of the number of events in the detector as a function of distance from the neutrino source, correcting for the falloff of the neutrino flux), which would be unmistakable proof for short baseline neutrino oscillations. The OscSNS proposal [40] involves building a cylindrical detector approximately 60 m from the SNS spallation target [41]. The 1.4 MW beam power of the SNS is a prodigious source of neutrinos from the decay of and at rest. These decays produce a well specified flux of neutrinos via ,  s and ,  s. The low duty factor of the SNS (~695 ns beam pulses at 60 Hz, ) is more than 1000 times less than that found at LAMPF. This smaller duty factor provides a reduction in backgrounds due to cosmic rays, and allows the induced events from decay to be separated from the and induced events from decay. The OscSNS detector will be centered at a location 60 meters from the SNS target, in the backward direction. The cylindrical detector design is based upon the LSND and MiniBooNE detectors and will consist of an 800-ton tank of mineral oil (with a small concentration of b-PBD scintillator dissolved in the oil) that is covered by approximately 3500 8-inch phototubes, yielding a photocathode coverage of 25%. Figure 21 shows the expected sensitivity for appearance after two and six calendar years of run time. The LSND allowed region is fully covered by more than .

(a)

(b)

Figure 21: The OscSNS sensitivity curves for the simulated sensitivity to oscillations after two (a) and six (b) calendar years of operation. Note that OscSNS has more than sensitivity to the LSND allowed region after 2 years of data collection.

The IsoDAR proposal [33] involves the construction of a high-intensity proton cyclotron near an underground detector to search for disappearance. The proposed cyclotron would accelerate ions up to 60 MeV/amu at an intensity of 5 mA. The 60 MeV proton beam then interacts on a ⁹Be target that is surrounded by a ⁷Li sleeve. Neutrons produced in the Be target interact in the sleeve to produce ⁸Li, which decay to produce with an average energy of 6.4 MeV. The can then interact in a kiloton underground detector, such as KamLAND [42] in Japan. By fitting for oscillations in the detector, IsoDAR would have excellent sensitivity to disappearance, as shown in Figure 22.

Figure 22: The sensitivity of the IsoDAR experiment to disappearance in a five-year physics run (from [33]).

6. Conclusions

The LSND and MiniBooNE experiments both observe excesses of candidate events at the level. Combining these two experiments results in a total excess that is over in significance. Whereas no standard model explanation for these excesses has yet been found, oscillation models with 3 light, mostly active neutrinos and heavy, mostly sterile neutrinos provide a fairly good fit to the world neutrino oscillation data, although there is some tension at present between appearance and disappearance oscillation experiments. Future experiments at Fermilab, CERN, ORNL, and in Japan will test this neutrino oscillation evidence and have the potential to prove the existence of light, sterile neutrinos. Such sterile neutrinos would have a big impact on particle physics, nuclear physics, and astrophysics.

Acknowledgments

The author thanks the Department of Energy for support. This work was made possible by the dedicated efforts of the LSND and MiniBooNE collaborations.

References

G. L. Fogli, E. Lisi, A. Marrone, D. Montanino, A. Palazzo, and A. M. Rotunno, “Addendum to ‘Solar neutrino oscillation parameters after first KamLAND results’,” Physical Review D, vol. 69, no. 1, Article ID 017301, 2004.

C. Athanassopoulos, L. B. Auerbach, D. A. Bauer, et al., “Candidate events in a search for oscillations,” Physical Review Letters, vol. 75, no. 14, pp. 2650–2653, 1995.

C. Athanassopoulos, L. B. Auerbach, R. L. Burman, et al., “Evidence for oscillations from the LSND experiment at the los alamos meson physics facility,” Physical Review Letters, vol. 77, pp. 3082–3085, 1996.

C. Athanassopoulos, L. B. Auerbach, R. L. Burman, et al., “Results on neutrino oscillations from the LSND experiment,” Physical Review Letters, vol. 81, pp. 1774–1777, 1998.

C. Athanassopoulos, L. B. Auerbach, R. L. Burman, et al., “Results on oscillations from pion decay in flight neutrinos,” Physical Review C, vol. 58, pp. 2489–2511, 1998.

A. Aguilar, L. B. Auerbach, R. L. Burman, et al., “Evidence for neutrino oscillations from the observation of appearance in a beam,” Physical Review D, vol. 64, no. 11, Article ID 112007, 2001.

A. A. Aguilar-Arevalo, C. E. Anderson, A. O. Bazarko, et al., “Unexplained excess of electronlike events from a 1-GeV neutrino beam,” Physical Review Letter, vol. 102, no. 10, Article ID 101802, 2009.

A. A. Aguilar-Arevalo, C. E. Anderson, A. O. Bazarko, et al., “Event excess in the MiniBooNE search for oscillations,” Physical Review Letter, vol. 105, no. 10, Article ID 181801, 2010.

A. A. Aguilar-Arevalo, C. E. Anderson, A. O. Bazarko, et al., “Improved search for oscillations in the MiniBooNE experiment,” Physical Review Letter, vol. 110, no. 10, Article ID 161801, 2013.

J. M. Conrad, C. M. Ignarra, G. Karagiorgi, M. H. Shaevitz, and J. Spitz, “Sterile Neutrino Fits to Short-Baseline Neutrino Oscillation Measurements,” Advances in High Energy Physics, vol. 2013, Article ID 163897, 26 pages, 2013.

K. Joachim, A. N. Pedro Machado, M. Michele, and T. Schwetz, “Sterile neutrino oscillations: the global picture,” Journal of High Energy Physics, vol. 2013, article 50, 2013.

K. N. Abazajian, M. A. Acero, S. K. Agarwalla, et al., “Light sterile neutrinos: a white paper,” http://arxiv.org/abs/1204.5379

P. Langacker, The Standard Model and Beyond, Taylor and Francis, Boca Raton, Fla, USA, 2010.

A. de Gouvea and W.-C. Huang, “Constraining the (low-energy) type-I Seesaw,” Physical Review D, vol. 85, Article ID 053006, 12 pages, 2012.

V. A. Kostelecky and M. Mewes, “Lorentz and CPT violation in neutrinos,” Physical Review D, vol. 69, no. 1, Article ID 016005, 2004.

T. Katori, V. A. Kostelecky, and R. Tayloe, “Global three-parameter model for neutrino oscillations using Lorentz violation,” Physical Review D, vol. 74, Article ID 105009, 2006.

S. N. Gninenko, “MiniBooNE anomaly and heavy neutrino decay,” Physical Review Letter, vol. 103, no. 24, Article ID 241802, 2009.

S. N. Gninenko and D. S. Gorbunov, “MiniBooNE anomaly, the decay and heavy sterile neutrino,” Physical Review D, vol. 81, no. 7, Article ID 075013, 2010.

ALEPH Collaboration, DELPHI Collaboration, L3 Collaboration, OPAL Collaboration, and SLD Collaboration, “Precision electroweak measurements on the Z resonance,” Physical Report, vol. 427, no. 5-6, p. 257, 2006.

C. Athanassopoulos, L. B. Auerbach, D. Bauer, et al., “The liquid scintillator neutrino detector and LAMPF neutrino source,” Nuclear Instruments and Methods in Physics Research A, vol. 388, no. 1-2, pp. 149–172, 1997.

R. L. Burman, M. E. Potter, and E. S. Smith, “Monte Carlo simulation of neutrino production by medium-energy protons in a beam stop,” Nuclear Instruments and Methods in Physics Research A, vol. 291, no. 3, pp. 621–633, 1990.

R. L. Burman, A. C. Dodd, and P. Plischke, “Neutrino flux calculations for the ISIS spallation neutron facility,” Nuclear Instruments and Methods in Physics Research A, vol. 368, no. 2, pp. 416–424, 1996.

R. A. Reeder, B. D. Dieterle, C. Gregory, F. Schaefer, and K. Schum, “Dilute scintillators for large-volume tracking detectors,” Nuclear Instruments and Methods in Physics Research A, vol. 334, pp. 353–366, 1993.

J. J. Napolitano, S. J. Freedman, G. T. Garvey, et al., “Construction and performance of a large area liquid scintillator cosmic ray anticoincidence detector,” Nuclear Instruments and Methods in Physics Research Section A, vol. 274, no. 1-2, pp. 152–164, 1989.

B. Achkara, R. Aleksanb, M. Aveniera, et al., “Search for neutrino oscillations at 15, 40 and 95 meters from a nuclear power reactor at Bugey,” Nuclear Physics B, vol. 434, pp. 503–534, 1995.

B. Armbruster, I. M. Blair, B. A. Bodmann, et al., “Upper limits for neutrino oscillations from muon decay at rest,” Physical Review D, vol. 65, Article ID 112001, 2002.

G. Mention, M. Fechner, T. Lasserre et al., “Reactor antineutrino anomaly,” Physical Review D, vol. 83, no. 7, Article ID 073006, 2011.

A. A. Aguilar-Arevalo, C. E. Anderson, A. O. Bazarko, et al., “Neutrino flux prediction at MiniBooNE,” Physical Review D, vol. 79, no. 7, Article ID 072002, 2009.

A.A. Aguilar-Arevaloe, C.E. Andersonp, L.M. Bartoszek, et al., “The MiniBooNE detector,” Nuclear Instruments and Methods in Physics Research A, vol. 599, no. 1, pp. 28–46, 2009.

M. Antonello, B. Baibussinov, P. Benetti, et al., “Experimental search for the LSND anomaly with the ICARUS detector in the CNGS neutrino beam,” The European Physical Journal C, vol. 73, Article ID 2345, 2013.

D. Adey, S. K. Agarwalla, C. M. Ankenbrandt, et al., “Neutrinos from stored Muons nuSTORM: expression of interest,” http://arxiv.org/abs/1305

M. Antonello, D. Bagliani B, and B. Baibussinov, “Baibussinov Search for, “anomalies” from neutrino and anti-neutrino oscillations at with muon spectrometers and large LAr-TPC imaging detectors,” http://arxiv.org/abs/1203.3432

A. Bungau, A. Adelmann, J. R. Alonso, et al., “Proposal for an electron antineutrino disappearance search using high-rate ⁸Li production and decay,” Physical Review Letter, vol. 109, no. 14, Article ID 141802, 2012.

A. A. Aguilar-Arevalo, A. O. Bazarko, S. J. Brice, et al., “Search for electron neutrino appearance at the scale,” Physical Review Letter, vol. 98, no. 23, Article ID 231801, 2007.

A. A. Aguilar-Arevalo, B. C. Brown, L. Bugel, et al., “A combined to and to oscillation analysis of the MiniBooNE excesses,” Physical Review Letters, vol. 110, no. 16, Article ID 161801, 2013.

M. Maltoni, T. Schwetz, M. Tortola, and J. Valle, “Ruling out four-neutrino oscillation interpretations of the LSND anomaly?” Nuclear Physics B, vol. 643, no. 1–3, pp. 321–328, 2002.

M. J. Conrad, C. W. Louis, and H. M. Shaevitz, “The LSND and MiniBooNE oscillation searches at high ,” Annual Review of Nuclear and Particle Science, vol. 63, p. 45, 2013.

P. Adamson, D. J. Auty, D. S. Ayres, et al., “Active to sterile neutrino mixing limits from neutral-current interactions in MINOS,” Physical Review Letter, vol. 107, no. 1, Article ID 011802, 2011.

G. Karagiorgi, “MicroBooNE: searching for new physics in the neutrino sector with a 100-ton-scale liquid argon TPC,” Journal of Physics, vol. 375, Article ID 042067, 2012.

G. T. Garvey, A. Green, C. Green, et al., “Measuring active-sterile neutrino oscillations with a stopped pion neutrino source,” Physical Review D, vol. 72, no. 9, Article ID 092001, 2005.

“The spallation neutron source (SNS) is an accelerator-based source built in Oak Ridge, Tennessee, by the U.S. DOE,” http://sns.gov/

S. Abe, T. Ebihara, S. Enomoto, et al., “Precision measurement of neutrino oscillation parameters with KamLAND,” Physical Review Letter, vol. 100, no. 22, Article ID 221803, 5 pages, 2008.