Neutrino physics is nowadays receiving more and more attention as a
possible source of information for the long-standing problem of new
physics beyond the Standard Model. The recent measurement of the
third mixing angle θ13 in the standard mixing oscillation scenario encourages us to pursue the still missing results on leptonic CP violation
and absolute neutrino masses. However, several puzzling measurements
exist which deserve an exhaustive evaluation. We will illustrate the
present status of the muon disappearance measurements at small L/E
and the current CERN project to revitalize the neutrino field in Europe
with emphasis on the search for sterile neutrinos. We will then illustrate the achievements that a double muon spectrometer can make with
regard to discovery of new neutrino states, using a newly developed
analysis.
1. Introduction
The unfolding of the physics of the neutrino is a long and exciting story spanning the last 80 years. Over this time the interchange of theoretical hypotheses and experimental facts has been one of the most fruitful demonstrations of the progress of knowledge in physics. The work of the last decade and a half finally brought a coherent picture within the Standard Model (SM) (or some small extensions of it), namely, the mixing of three neutrino flavour states with three ν1, ν2, and ν3 mass eigenstates. The last unknown mixing angle, θ13, was recently measured [1–4], but still many questions remain unanswered to completely settle the scenario: the absolute masses, the Majorana/Dirac nature, and the existence and magnitude of leptonic CP violation. Answers to these questions will beautifully complete the (standard) three-neutrino model, but they will hardly provide an insight into new physics Beyond the Standard Model (BSM). Many relevant questions will stay open: the reason for neutrinos, the relation between the leptonic and hadronic sectors of the SM, the origin of Dark Matter, and overall where and how to look for BSM physics. Neutrinos may be an excellent source of BSM physics and their story is supporting that at length.
There are actually several experimental hints for deviations from the “coherent” picture described above. Many unexpected results, not statistically significant on a single basis, appeared also in the last decade and a half, bringing attention to the hypothesis of sterile neutrinos [5]. A recent White Paper [6] contains a comprehensive review of these issues. In this paper we will focus on one of the most intriguing and long-standing unresolved result: the unexpected oscillation of neutrinos at relatively small values of the ratio L/E (distance in km, energy in GeV), corresponding to a scale of 𝒪(1) eV2, incompatible with the much smaller values related to the atmospheric |Δm322|≃2.4×10-3 eV2 and to the solar Δm212≃8×10-5 eV2 scales.
The first unexpected measurement came from an excess of ν¯e originating from an initial ν¯μ beam from Decay At Rest (LNSD [7]). The LSND experiment saw a 3.8 σ effect. The subsequent experiment with νμ(ν¯μ) beam from accelerator, MiniBooNE [8–11], despite confirming an independent 3.8 σ effect after sustained experimental work, was unable to draw conclusive results on the origin of the LSND effect having observed an excess at higher L/E in an energy region where background is high.
In recent years many phenomenological studies were performed by analyzing the LSND effect together with similar unexpected results coming from the measurement of lower than expected rates of ν¯e and νe interactions (disappearance), either from (a) near-by nuclear reactors (ν¯e) [12–14] or (b) from Mega-Curie K-capture calibration sources in the solar νe gallium experiments [15–21]. These νe(ν¯e) disappearance measurements, all at the statistical level of 3-4 σ, could also be interpreted [22] as oscillations between neutrinos at large Δm2≃1 eV2. Several attempts were then tried to reach a coherent picture in terms of mixing between active and sterile neutrinos, in 3+1 and 3+2 [23, 24] or even 3+1+1 [25, 26] or 3+3 [27] models, as extensions of the standard three-neutrino model. We refer to [28–30] as the most recent and industrious works where a very crucial issue is raised: “a consistent interpretation of the global data in terms of neutrino oscillations is challenged by the non-observation of a positive signal in νμ disappearance experiments” [29]. In fact, in any of the above models, essential information comes from the disappearance channels (νμ or νe), which is one of the cleanest channels to measure the oscillation parameters (see Section 3 for details).
The presence of additional sterile states introduces quite naturally appearance and disappearance phenomena involving the flavour states in all channels. In particular, a νμ disappearance effect has to be present and possibly measured. It turns out that only old experiments and measurements are available for Charged Current (CC) νμ interactions at small L/E [31]. The CDHS experiment reported in 1984 the nonobservation of νμ oscillations in the Δm2 range 0.3 eV2 to 90 eV2. Their analyzed region of oscillation did not span, however, low values of the mixing parameter down to around 0.1 in sin2(2θ). More recent results are available on νμ disappearance from MiniBooNE [32], a joint MiniBooNE/SciBooNE analysis [33, 34] and the Long-Baseline MINOS experiment [35]. These results slightly extend the νμ disappearance exclusion region, however, still leaving out the small-mixing region. Similar additional constraints on νμ disappearance could possibly come from the analysis of atmospheric neutrinos in IceCube [36, 37].
Despite this set of measurements being rather unsatisfactory when compared with the corresponding LSND allowed region that lies at somewhat lower values of the mixing angle, they are still sufficient to introduce tensions in all the phenomenological models developed so far (see, e.g., [6, 28–30, 38] for comprehensive and recent reviews). Therefore it is mandatory to set up a new experiment able to improve the small-mixing angle region exclusion by at least one order of magnitude with respect to the current results. In such a way one could also rule out the idea that the mixing angle extracted from LSND is larger than the true value due to a data overfluctuation. Once again, the νμ disappearance channel should be the optimal one to perform a full disentangling of the mechanism given the strong tension between the νe appearance and νμ disappearance around Δm2≃1 eV2. In fact, whereas the LSND effect may be confirmed by a more accurate νe oscillation measurement, only the presence of a νμ oscillation pattern could shed more light on the nature and the interpretation of the effect.
Further, in the paper we will briefly discuss the newly proposed CERN experiment, following a detailed analysis of possible ways to measure the νμ disappearance channel. In particular the evaluation of the νμ disappearance rates at two different sites leads us naturally to take their ratio, elucidating the possibility to observe depletions and/or excesses. Throughout the paper we will focus on the need for a new νμ CC measurement corresponding to an increase of one order of magnitude in the sensitivity to the mixing parameter, exploiting the determination of the muon charge for a proper evaluation of the different models and the separation of νμ and ν¯μ in the same neutrino beam. Finally, we will pay particular attention to the consistency and robustness (in statistical terms) of the results.
2. The CERN Experimental Proposal
The need for a definitive clarification on the possible existence of a neutrino mass scale around 1 eV has brought up several proposals and experimental suggestions exploiting the sterile neutrino option by using different interaction channels and refurbished experiments. In the light of the considerations discussed in the previous section there are essentially two sets of experiments which must be redone: (a) the measurement of ν¯e neutrino fluxes at reactors (primarily the ILL one) together with refined and detailed computations of the flux simulations (see, e.g., [39, 40]); (b) the appearance/disappearance oscillation measurements at low L/E with a standard muon neutrino beam with its intrinsic electron neutrino component.
There is actually another interesting option which comes from the Neutrino Factory studies and the very recently submitted EOI from νSTORM [41]. We refer to [42] for a comprehensive review of the corresponding possible νe and νμ disappearance effects. It is interesting to note that our figures of merit about νμ disappearance are rather similar to or even slightly more competitive than those illustrated in [42] (e.g., compare the exclusion regions in Figure 6 of [42] and those in Figure 7 of this paper despite the use of different C.L.), not forgetting the rather long time needed to set up the νSTORM project.
Coming to experiments with standard beams, investigations are underway at CERN where two physics proposals [43, 44] were submitted in October 2011 and later merged into a single technical proposal (ICARUS-NESSiE, [45]). CERN has subsequently set up working groups for the proton beam extraction from the SPS, the secondary beam line, and the needed infrastructure/buildings for the detectors. The work was reported in a recent LOI [46].
The experiment is based on two identical Liquid Argon (LAr)—Time Projection Chambers (TPC) [44] complemented by magnetized spectrometers [43] detecting electron and muon neutrino events at far and near positions, 1600 m and 460 m away from the proton target, respectively. The project will exploit the ICARUS T600 detector, the largest LAr-TPC ever built of about 600 ton mass, now presently in the LNGS underground laboratory where it was exposed to the CNGS beam. It is supposed to be moved at the CERN far position. An additional 1/4 of the T600 detector (T150) would be constructed from scratch as a clone of the original one, except for the dimensions, and located in the near site. Two spectrometers would be placed downstream of the two LAr-TPC detectors to greatly enhance the physics reach. The spectrometers will exploit a bipolar magnet with instrumented iron slabs, and a newly designed air-core magnet, to perform charge identification and muon momentum measurements in an extended energy range (from 0.5 GeV or less to 10 GeV), over a transverse area larger than 50 m2.
While the LAr-TPCs will mainly perform a direct measurement of electron neutrinos [47], the spectrometers will allow an extended exploitation of the muon neutrino component, with neutrino/antineutrino discrimination on an event-by-event basis.
2.1. The Neutrino Beam
The proposed new neutrino beam will be constructed in the SPS North Area [46]. The setup is based on a 100 GeV proton beam with a fast extraction scheme providing about 3.5·1013 protons/pulse in two pulses of 10.5 μs durations separated by 50 ms for a sample of about 4.5·1019 protons on target (p.o.t.) per year (pulses of 10.5 μs duration are normally put in coincidence with the fast response of the spectrometers’ detectors and efficiently used to reject the cosmic ray background. See, e.g., [48] where a time resolution of less than 2 ns is reported for the detectors used in the OPERA experiment). A target station will be located next to the so-called TCC2 target zone, 11 m underground, followed by a cylindrical He-filled decay pipe with a length of about 110 m and a diameter of 3 m. The beam dump of 15 m in length will be composed of iron blocks with a graphite inner core. Downstream of the beam dump a set of muon chambers stations will act as beam monitors. The beam will point upward, with a slope of about 5 mrad, resulting in a depth of 3 m for the detectors in the far site.
The current design of the focusing optics includes a pair of pulsed magnetic horns operated at relatively low currents. A graphite target of about 1 m in length is deeply inserted into the first horn allowing a large acceptance for the focusing of low momentum pions emitted at large angles. This design allows production of a spectrum peaking at about 2 GeV thus matching the most interesting domain of Δm2 with the detector locations at 460 and 1600 m from the target.
The charged current event rates for νμ and ν-μ at the near and far detectors are shown in Figure 1 for the positive and negative focusing configuration.
Expected neutrino CC interactions in the no-oscillation hypothesis for positive polarity (a) and negative polarity (b) for the new proposed CERN neutrino beam (elaborated from [46]) and for an integrated luminosity of 1 year. In the first (second) row the rates are shown at the near (far) position. In the third row the rates expected at the near site are displayed as function of log10(1/E), used in the paper (see Section 3). Note that the distributions of the antineutrino rates in case of positive polarity are multiplied by a factor 5 to allow a better visual inspection.
A relevant contamination of νμ in the negative polarity configuration is visible especially at high energy. This component arises as a result of the decays of high energy poorly defocused mesons produced at small angles. The charge discrimination of the magnetic system described below will allow an efficient discrimination of these two components with a charge confusion below or of the order of 1% from sub-GeV (0.3–0.5 GeV) up to momenta around 8–10 GeV [43].
2.2. Spectrometer Requirements
The main purpose of the spectrometers placed downstream of the LAr-TPC is to provide charge identification and momentum reconstruction for the muons produced in neutrino interactions occurring in the LAr volume or in the magnetized iron of the spectrometers. In order to perform this measurement with high precision and in a wide energy range, from sub-GeV to multi-GeV, a massive iron-core dipole magnet (ICM) is coupled to an air-core magnet (ACM) in front of it [43]. Low momentum muons will be measured by the ACM, while the ICM will be employed at higher momenta.
As considered in the previous sections the definition of two sites, near and far, constitutes a fundamental issue for each physics program which aims to perform any sterile neutrino search. The two layouts have to be as similar as possible in order to allow an almost complete cancellation of the systematic uncertainties when comparing the measurements made at the near and far sites. Hence the near spectrometer will be an exact clone of the far one, with identical thickness along the beam axis but a scaled transverse size. A sketch of the possible far site NESSiE detector is shown in Figure 2.
Sketch of the far site spectrometer that could be built by extensively reusing materials available from the OPERA spectrometers [51]. The basic elements of the new ACM concept are also depicted (courtesy of the NESSiE Collaboration). Neutrinos are traveling from right to left.
The key feature of the ACM is the large geometric and momentum acceptance. The need of a low momentum threshold for muon reconstruction can be met using a magnet in air. The only dead material along the muon path is given by the conductors needed to generate the magnetic field and the position detectors instrumenting the magnet itself. For the conductors the use of aluminum instead of copper is preferable due to the lower Z and a density lower by a factor 3. The magnetic field needed in the low momentum range covered by the ACM is in the range 0.1–0.15 T. A spatial resolution in the range of 0.1–1 mm can be reached using drift tubes [49] as high precision trackers in combination with scintillator strip detectors [50]. These could provide the external trigger needed by the drift tubes and a coarse spatial measurement in the nonbending direction, that is, the direction parallel to the drift tubes. Silicon photomultiplier (SiPM) devices may eventually be used to read out some scintillator planes embedded in the magnetic field.
The general layout of the two OPERA [51] iron spectrometers fulfills the requirements of the ICM detector, and they could eventually be used for the CERN project putting the two magnets one after the other in order to obtain a total of 4 m longitudinal thickness of iron.
The OPERA spectrometers are built assembling vertical iron plates (slabs) in a planar structure of 875 cm (width) × 800 cm (height). Each passive plane is made out of seven adjacent iron slabs. Resistive Plate Chambers (RPC) [52] are sandwiched between iron planes; 21 RPC detectors, arranged in seven rows and three columns, are used in each active plane. Each RPC detector has a rectangular shape and covers an area of about 3.2 m2. The RPCs provide tracking measurements with about 1 cm resolution using the digital read-out of strips with a 2.6 cm pitch in the bending direction and 3.5 cm in the nonbending direction. The magnet is made of two arms with 22 RPC layers alternating with iron layers. Top and bottom iron yokes are connecting the arms. Copper coils surrounding the yokes are used to generate a magnetic field of about 1.5 T in the iron circuit. A total of 924 RPC chambers are needed to instrument each spectrometer. They could be recovered from the OPERA spectrometers and reused. These RPC chambers are standard 2 mm gap chambers with bakelite electrodes with resistivity in the range from 1011 to 5×1012 Ω·cm at T=20°C. In OPERA they are operated in the streamer regime with a gas mixture made of Ar/C2H2F4/I-C4H10/SF6 in the volume ratios of 75.4/20/4/0.6 at five refills/day. The high-amplitude streamer signals produced by charged tracks crossing the gas volume allow to house the front-end discriminators in racks placed on top of the spectrometer.
3. The νμ Disappearance Analysis
The use of at least two sites where neutrino interactions can occur is mandatory to observe the oscillation pattern, as established by almost all the most recent neutrino experiments. Indeed the disappearance probability due to an additional sterile neutrino is given by the usual two-flavour formula:
(1)P(να⟶να)=1-sin2(2θ)·sin2(1.267·Δm2·LE),(Linkm,EinGeV),
where the disappearance of flavour α is due to the oscillation of neutrino mass states at the Δm2 scale and at an effective mixing angle θ that can be simply parametrized as a function of the elements of a 3+1 extended mixing matrix. As L is fixed by the experiment location, the oscillation is naturally driven by the neutrino energy, with an amplitude determined by the mixing parameter.
The disappearance of muon neutrinos due to the presence of an additional sterile state depends only on terms of the extended PMNS [53, 54] mixing matrix (Uαi with α=e, μ, τ and i=1,…,4) involving the νμ flavor state and the additional mass eigenstate 4. In a 3+1 model at Short Baseline (SBL) we have
(2)P(νμ⟶νμ)SBL3+1=1-4|Uμ4|2(1-|Uμ4|2)sin2Δm412L4E.
In contrast, appearance channels (i.e. νμ→νe) are driven by terms that mix up the couplings between the initial and final flavour states and the sterile state yielding a more complex picture:
(3)P(νμ⟶νe)SBL3+1=4|Uμ4|2|Ue4|2sin2Δm412L4E.
This also holds in extended 3+n models.
It is interesting to notice that the appearance channel is suppressed by two more powers in |Uα4|. Furthermore, since νe or νμ appearance requires |Ue4|>0 and |Uμ4|>0, it should be naturally accompanied by a corresponding νe and νμ disappearance. In this sense the disappearance searches are essential for providing severe constraints on the models of the theory (a more extensive discussion on this issue can be found, e.g., in Section 2 of [42]).
It must also be noted that the number of νe neutrinos depends on the νe→νs disappearance and νμ→νe appearance and, naturally, from the intrinsic νe contamination in the beam. On the other hand, the amount of νμ neutrinos depends only on the νμ→νs disappearance and νe→νμ appearance, but the latter is much smaller due to the fact that the νe contamination in νμ beams is usually at the percent level. Therefore in the νμ disappearance channel the oscillation probabilities in both near and far detectors can be measured without any interplay of different flavours, that is, by the same probability amplitude.
By taking into account the lengths L of the locations given in the CERN proposal, LNear=460 m and LFar=1600 m, we may plot the disappearance probability as a function of the more convenient variable log10(1/E), using E in MeV to avoid singularities. As an example, in Figure 3 the disappearance probability is shown for the near and far sites, by assuming an amplitude that corresponds to the averaged reactor anomaly ν¯e disappearance, R=0.927±0.023 [12–14] and Δm2=1 eV2 (it is worthwhile to note that the disappearance amplitude affects the sensitivity to the mixing angle that in turn depends mainly on the statistical extent of the data. Therefore the same evidence can be obtained for smaller amplitudes, that is, smaller effective mixing angles, by increasing the data sample. Finally the ultimate limit of an experiment is set by the intrinsic systematic errors corresponding to the kinematics of the neutrino interaction and to the muon reconstruction sensitivity. Figure 8 below illustrates the issue). It appears that at Δm2≈1 eV2 oscillations are already visible at the near site. Such a behaviour may turn out to help the actual measurement by plotting the ratio between the far and near sites data as a function of log10(1/E) (Figure 4). The Far/Near ratio of the oscillation probabilities well illustrates the pattern that could be observed at the CERN Short-Baseline experiment: a decrease above 1 GeV and an increase below 1 GeV of the neutrino energy!
Disappearance probabilities in the two-flavour limit at near (a) and far (b) sites, at 460 and 1600 m, respectively, by using the amplitude provided by the reactor anomaly, 0.146, and the mass scale Δm2=1 eV2. The x-axis corresponds to log10(1/Eν), with E in MeV.
The ratio of the disappearance probabilities at near and far sites, 460 and 1600 m, respectively, in the two-flavour limit, by using the amplitude provided by the reactor anomaly, 0.146, and the mass scale Δm2=1 eV2. The x-axis corresponds to log10(1/Eν), with E in MeV.
The disappearance patterns have to be evaluated after folding cross-sections and efficiencies with the neutrino fluxes and detector geometry. The Far/Near pattern does not depend at first order either on cross-sections and efficiencies or on the fluxes as far as they are similar at the two sites and within the use of two identical detector systems. However, due to the different ν and ν¯ cross-sections, errors are enhanced when neutrinos and antineutrinos are not identified and separated on an event-by-event basis. Therefore the use of magnetic detectors is mandatory to keep systematics under control. Then the largest source of error consists of the smearing introduced by the reconstruction procedure. Eventually this issue will be addressed by considering a conservative systematic effect in the muon momentum reconstruction.
Figure 5(a) illustrates a realistic observation of the nonoscillation probability when convoluted with the estimated CERN νμ beam flux and the neutrino cross-sections, by using the detector systems described above and one year of data taking. The flatness is lost as expected. However a specific case of oscillation is drawn, confirming that the workable behaviour is still present. The correlated estimator, the double ratio parameter, defined as (Far/Near)oscillated/(Far/Near)unoscillated, is also shown in the bottom of the figure. It can be used for the statistical analysis. In the following we will demonstrate how the oscillation behaviour of the Far/Near parameter will allow an increase of an order of magnitude in the sensitivity to the mixing angle, taking into account the systematic errors, too.
The oscillated and unoscillated νμ event distributions, parametrized in log10(1/Eν), for a luminosity of 4.5×1019 p.o.t. of the CERN νμ beam flux over the two-magnet system described in the text, are plotted with their statistical errors. (a) The far/near collected ratios are drawn. The no-oscillation shape is due to the convolution of the neutrino cross-section and the fluxes at the near and far site, while the oscillated shape is due to the sterile inclusion (with Δm2=1 eV2 and mixing amplitude equal to 0.146). (b) The double ratio (Far/Near)oscillated/(Far/Near)unoscillated is shown for the same sterile assumption.
Another parameter can be used in case measurement of Neutral Current (NC) events becomes available, either with the proposed LAr detectors or even with a less refined and massive detector. The NC/CC ratio becomes then a valuable estimator [43] that may enter into the evaluation of the significance of the sterile neutrino observation. A study that includes the information from Neutral Current events is out of the scope of the present paper. However for completeness we would like to briefly outline here the issue as it was originally described in [43].
The observation of a depletion in NC events would be a direct and manifest signature of the existence of sterile neutrinos. In fact the NC event rates are unaffected by standard neutrino mixing being flavour-blind such that their disappearance could only be explained by νμ(νe)→νs transitions. Even if NC events, either from νe or νμ, are efficiently detected by a Liquid Argon detector or a less refined detector, the transition rate measured with NC events has to agree with the νμ CC disappearance rate once the νμ→ντ and νμ→νe contributions have been subtracted. Therefore the NC disappearance is measured at best by the double ratio: (NC/CC)Far/(NC/CC)Near. The double ratio is the most robust experimental quantity to detect NC disappearance, once CCNear and CCFar are precisely measured thanks to the spectrometers, at the near and far locations. For that it is mandatory to disentangle νμ and ν¯μ contributions.
4. Simulation and Results
The magnetic detector system that has to be developed for the νμ disappearance measurement should take into account all the considerations depicted in previous sections. The system developed by the NESSiE Collaboration [43] is actually well suited since it couples a very powerful high-Z magnet for the momentum measurement via range to a low-Z magnet, to extend the useful muon momentum interval as low as possible, to allow charge discrimination on an event-by-event basis, and to allow NC event measurement whether coupled to an adequate (but not necessarily highly performant and with large mass) detector to identify NC events.
In the following we concentrate on the results achievable by the Far/Near estimators described in the previous section, by exploiting a full neutrino beam simulation, an up-to-date neutrino interaction model (GENIE 2.6 [55]), with detailed interactions in low-Z and high-Z magnets and a generic geometry. The Iron Core Magnets have been simulated by a massive cube of iron with a mass equivalent to the OPERA magnets and by using a 90% inner fiducial volume (770 tons and 330 tons for the far and near locations, resp.). Single hits have been extracted at the position of the RPC detectors and a realistic resolution has been included (0.75 cm for the single hit). The ACM has been simulated as an empty 1 m deep region with 0.1 T magnetic field, limited by two vertical aluminum slabs. The overall simulation corresponds to the GLoBES software [56] (version 3.0.11), to which it was partially compared and cross checked.
The muon reconstruction requires a minimal penetration length that corresponds to a 0.5 GeV cut. A conservative resolution of 10% was used for the momentum evaluation, while the charge misidentification was fixed at the 1% level. Systematic effects due to miscalibration of the detectors have been conservatively taken at 0.5% and 1% levels, which include systematics due to the relative flux at the near and far positions.
In Figure 6 the double ratio Far/Near is shown as obtained by using the reconstructed muon momentum, for different values of Δm2 scales. The relevant behaviour of the estimator is confirmed for Δm2 values above ≈2 eV2 also in this very realistic and conservative approach. In fact, by using the muon reconstructed momentum in CC events instead of the true νμ energy, the oscillations properties will be affected by the kinematics and the interaction processes (it is worthwhile to note that there are mostly quasi-elastic interactions in the accounted-for energy range).
The double ratio of the disappearance probabilities at the near and far sites, 460 and 1600 m, respectively, by using the amplitude provided by the reactor anomaly, 0.146, and mass scales Δm2 ranging from 0.1 to 5 eV2 as a function of log10(1/pμ), where “p” is the reconstructed muon momentum in MeV/c. We apply cuts on the minimum muon path length (25 cm of range in iron) and the minimum momentum (500MeV/c). A 90% fiducial acceptance of the magnet volume and an uncorrelated 1% systematic error (from OPERA experience [57]) are also considered. Data collection corresponds to 1 year, that is, 4.5·1019 p.o.t. The error bars correspond to the statistical errors. The blue band for Δm2=2 eV2 is obtained by including the systematic error summed quadratically to the statistical one, the same result being obtained for the other mass scales.
The estimated limits at 95% C.L. for νμ disappearance at a Short-Baseline beam at CERN for several luminosity running periods and different beam polarities, with a two-site massive spectrometer (770 tons and 330 tons, resp.) with 90% inner fiducial volume. (a) Refers to the positive polarity beam. The continuous (dashed) lines correspond to the sensitivity limits obtained with the log10(1/p)(p) variable. 3 years correspond to 13.5·1019 p.o.t., 5 years to 22.5·1019 p.o.t., and 10 years to 45.0·1019 p.o.t. The exclusion limit from combined MiniBooNE and SciBooNE νμ disappearance result at 90% C.L. from [33] is shown for comparison by the black curve in the right. (b) Refers to the negative polarity beam. Sensitivity limits are evaluated with the log10(1/p) variable. Clearly the negative polarity run allows the contemporary analysis of the ν¯μ and νμ disappearance exclusion regions thanks to the disentangling of the muon charge on an event-by-event basis. The black curve in the right shows for comparison the central value of the sensitivity at 90% C.L. from combined MiniBooNE and SciBooNE ν¯μ disappearance result [34].
Systematic effects for the estimated limits at 95% C.L. for νμ disappearance at a Short-Baseline beam at CERN for 5 (positive polarity, (a)) and 10 (negative polarity, (b)) years of data taking, with a two-site massive spectrometer (770 tons and 330 tons, resp.) with 90% inner fiducial volume. The systematic uncertainty, corresponding to each of the exclusion contours, ranging from 0% to 1%, is indicated in the legend as “sys”. The exclusion regions are evaluated by using the log10(1/p) (continuous lines) and p (dashed lines) variables in case of the positive polarity beam data taking.
Finally, in Figure 7 the estimated limits at 95% C.L. on νμ disappearance that can be achieved via the Far/Near estimator are shown for different data periods (3, 5, and 10 years, corresponding to 13.5·1019, 22.5·1019, and 45.0·1019 p.o.t., resp.). The different results for νμ and ν¯μ beams were evaluated using the two variables, p and log10(1/p). In negative polarity runs the muon charge identification allows an independent, simultaneous, and similar-sensitivity measurements of the ν¯μ and νμ disappearance rates, due to the large νμ contamination in the ν¯μ beam.
The parameterization used for the neutrino energy, log10(1/E), which we conservatively prefer to address as log10(1/pμ-rec), provides slightly better limits in almost all the (Δm2, sin2(2θ)) excluded regions. This is partially due to the Gaussian shape sensitivity of 1/p when cuts are applied to the corresponding variable. Moreover, the elucidation of the behaviour of the Far/Near ratio in terms of depletions and excesses in different regions of the spectra allows a better comparison between unoscillated and oscillated hypotheses.
As a result more than an order of magnitude improvement can be obtained in the sensitivity to the mixing parameter space between standard neutrino and sterile ones with respect to today’s limits over the whole Δm2 range investigated. Specifically the mixing angle sensitivity (at 95% C.L.) reaches 2.2×10-3 at Δm2=1 eV2, while at full mixing a Δm2=0.03 eV2 sensitivity at 95% C.L. is obtained already with only 3 years of data collection with the neutrino beam.
Slightly worse results are obtained by collecting data with a negative polarity beam. The limiting mixing angle sensitivity at 95% C.L. is around 10-2 at Δm2=1 eV2. At full mixing a Δm2=0.1 eV2 sensitivity at 95% C.L. is obtained with 3 years of data collection with the antineutrino beam, the limiting factor being the intensity of the ν¯μ beam. However it is worthwhile to note that by operating in negative polarity similar sensitivities can be obtained at the same time also for the neutrino component due to the large contamination of ν in the anti-ν beams.
The relevance of statistics in presence of systematic effects is depicted in Figure 8. It is evident that the overall error, dominated by the sensitivity on the measurement of the muon momentum, puts an intrinsic limitation to the data statistics which can be collected. As expected the useful luminosity that can be collected is asymptotically limited by systematics. The positive effect of the variable log10(1/pμ-rec) is more evident when larger errors are in place.
5. Conclusions
Neutrino physics is receiving more and more attention as a venue for the long standing search for new physics beyond the Standard Model. The current anomalies which do not fit into the established standard scenario with 3 neutrinos deserve refined studies and experiments. The CERN proposal for a new Short-Baseline experimental project is a very valuable one. We illustrated the current critical tensions in the muon neutrino disappearance field and the achievements that can be obtained within the CERN project. Specifically, by considering 𝒪(1) kton massive spectrometers, an improvement by an order of magnitude can be obtained in the sensitivity to the mixing parameter space between standard neutrinos and sterile ones with respect to today’s limits. Conversely, a possible νμ disappearance signal will be essential to measure the relevant physical parameters and to fully disentangle the different sterile models.
An effective analysis can be performed with a two-site experiment by using muon spectrometers with a low-Z part that allows clean charge identification on an event-by-event basis and with a massive part allowing clean momentum measurement through range. Such a kind of spectrometer is under study by the NESSiE Collaboration and might be available with a limited investment. The performances of these kinds of detectors, enlightened by the type of analysis developed in this paper, are suitable to put a definitive result on the sterile neutrino issue at the eV mass scale.
Acknowledgments
The authors are heartily dependent of the contributions of the NESSiE, ICARUS, and CERN-CENF groups in developing the CERN project and the experimental proposals. The authors wish also to warmly thank the encouragements of Marzio Nessi and Carlo Rubbia in supporting such studies. The authors are finally in debt to Maury Goodman for suggested criticisms and a careful proof-reading of the paper.
AnF. P.BaiJ. Z.BalantekinA. B.Observation of electron-antineutrino disappearance at Daya Bay2012108177171803AhnJ.ChebotaryovS.ChoiJ. H.Observation of reactor electron antineutrino disappearance in the RENO experiment2012108196191802AbeY.AberleC.AkiriT.Indication of reactor v̅e disappearance in the double chooz experiment20121087131801AbeK.AbgrallN.AjimaY.Indication of electron neutrino appearance from an accelerator-produced off-axis muon neutrino beam201110748041801PontecorvoB.Neutrino experiments and the problem of conservation of leptonic charge1967531717translation in Soviet Physics—Journal of Experimental and Theoretical Physics, vol. 26, p. 984, 1968AbazajianK. N.AceroM. A.AgarwallaS. K.Light sterile neutrinos: a white paperhttp://arxiv.org/abs/1204.5379AguilarA.AuerbachL. B.BurmanR. L.Evidence for neutrino oscillations from the observation of v̅e appearance in a v̅μ beam2001641122112007Aguilar-ArevaloA. A.BazarkoA. O.BriceS. J.Search for electron neutrino appearance at the Δm2 ~ 1 eV2 scale20079823710.1103/PhysRevLett.98.231801231801Aguilar-ArevaloA. A.AndersonC. E.BriceS. J.BrownB. C.BugelL.ConradJ. M.DharmapalanR.DjurcicZ.FlemingB. T.FordR.GarciaF. G.GarveyG. T.MirabalJ.GrangeJ.GreenJ. A.ImlayR.JohnsonR. A.KaragiorgiG.KatoriT.KobilarcikT.LindenS. K.LouisW. C.MahnK. B. M.MarshW.MaugerC.MetcalfW.MillsG. B.MooreC. D.MousseauJ.NelsonR. H.NguyenV.NienaberP.NowakJ. A.OsmanovB.PavlovicZ.PerevalovD.PollyC. C.RayH.RoeB. P.RussellA. D.SchiratoR.ShaevitzM. H.SorelM.SpitzJ.StancuI.StefanskiR. J.TayloeR.TzanovM.Van De WaterR. G.WasckoM. O.WhiteD. H.WilkingM. J.ZellerG. P.ZimmermanE. D.Event excess in the MiniBooNE search for ν̅μ→ν̅e oscillations2010105182-s2.0-7804923382810.1103/PhysRevLett.105.181801181801PollyC. C.Improved statistics in the search for ν̅μ→ν̅e oscillations in MiniBooNE2013235-2361207213Aguilar-ArevaloA.BrownB. C.BugelL.Improved search for v̅μ→v̅e oscillations in the MiniBooNE experiment20131106161801MentionG.FechnerM.LasserreTh.MuellerA. Th.LhuillierD.CribierM.LetourneauA.Reactor antineutrino anomaly2011832007300610.1103/PhysRevD.83.073006LhuillierD.Reactor flux calculations2013235-236111610.1016/j.nuclphysbps.2013.03.004LasserreT.Testing the reactor and gallium anomalies with intense (anti)neutrino emitters2013235-23621421910.1016/j.nuclphysbps.2013.04.013AbdurashitovJ. N.GavrinV. N.GorbachevV. V.Measurement of the solar neutrino capture rate with gallium metal. III. Results for the 2002–2007 data-taking period200980116015807AbdurashitovJ. N.GavrinV. N.GinnS. V.GorbachevV. V.IbragimovaT. V.KalikhovA. V.KhairnasovN. G.KnodelT. V.KornoukhovV. N.MirmovI. N.ShikhinA. A.VeretenkinE. P.VermulV. M.YantsV. E.ZatsepinG. T.BowlesT. J.NicoJ. S.TeasdaleW. A.WarkD. L.CherryM. L.KaraulovV. N.LevitinV. L.MaevV. I.NazarenkoP. I.Shkol'NikV. S.SkorikovN. V.ClevelandB. T.DailyT.DavisR.Jr.LandeK.LeeC. K.WildenhainP. W.KhomyakovY. S.ZvonarevA. V.ElliottS. R.WilkersonJ. F.The Russian-American Gallium Experiment (SAGE) Cr neutrino source measurement19967723470847112-s2.0-4243772839AbdurashitovJ. N.GavrinV. N.GirinS. V.Measurement of the response of a gallium metal solar neutrino experiment to neutrinos from a 51Cr source1999592246226310.1103/PhysRevC.59.2246AbdurashitovJ. N.GavrinV. N.GirinS. V.Measurement of the response of a Ga solar neutrino experiment to neutrinos from a 37Ar source20067312045805KaetherF.HampelW.HeusserG.KikoJ.KirstenT.Reanalysis of the Gallex solar neutrino flux and source experiments2010685147542-s2.0-7584912462210.1016/j.physletb.2010.01.030AnselmannP.FockenbrockR.HampelW.First results from the 51Cr neutrino source experiment with the GALLEX detector199534244045010.1016/0370-2693(94)01586-2HampelW.HeusserG.KikoJ.Final results of the 51Cr neutrino source experiments in GALLEX19984201-211412610.1016/S0370-2693(97)01562-1GiuntiC.LavederM.Statistical significance of the gallium anomaly20118362-s2.0-7996109805910.1103/PhysRevC.83.065504065504KoppJ.MaltoniM.SchwetzT.Are there sterile neutrinos at the eV scale?201110792-s2.0-8005195701710.1103/PhysRevLett.107.091801091801GiuntiC.LavederM.3+1 and 3+2 sterile neutrino fits20118471010.1103/PhysRevD.84.073008073008KuflikE.McDermottS. D.ZurekK. M.Neutrino phenomenology in a 3+1+1 framework2012863033015NelsonA. E.Effects of CP violation from neutral heavy fermions on neutrino oscillations, and the LSND/MiniBooNE anomalies20118452-s2.0-8005346993710.1103/PhysRevD.84.053001053001ConradJ. M.IgnarraC. M.KaragiorgiG.ShaevitzM. H.SpitzJ.Sterile neutrino fits to short-baseline neutrino oscillation measurements201320132610.1155/2013/163897163897KoppJ.MachadoP. A. N.MaltoniM.SchwetzT.Sterile neutrino oscillations: the global picture20132013, article 05010.1007/JHEP05(2013)050SchwetzT.Status of sterile neutrino oscillations2013235-236229235GiuntiC.LavederM.LiY. F.LiuQ. Y.LongH. W.Update of short-baseline electron neutrino and antineutrino disappearance2012861114113014DydakF.FeldmanG. J.GuyotC.A search for νμ oscillations in the Δm2 range 0.3–90 eV2198413428110.1016/0370-2693(84)90688-9Aguilar-ArevaloA. A.AndersonC. E.BazarkoA. O.Search for muon neutrino and antineutrino disappearance in MiniBooNE20091036510.1103/PhysRevLett.103.061802061802MahnK. B. M.NakajimaY.Aguilar-ArevaloA. A.Dual baseline search for muon neutrino disappearance at 0.5 eV2< Δm2< 40 eV2201285310032007ChengG.HuelsnitzW.Aguilar-ArevaloA. A.Dual baseline search for muon antineutrino disappearance at 0.1 eV2< Δm2< 100 eV2201286514052009AdamsonP.AndreopoulosC.AutyD. J.Search for sterile neutrino mixing in the MINOS long-baseline experiment2010815052004RazzaqueS.SmirnovA. Y.Searches for sterile neutrinos with IceCube DeepCore201285093010EsmailiA.HalzenF.PeresO.Constraining sterile neutrinos with AMANDA and IceCube atmospheric neutrino data2012201211, article 041PalazzoA.Phenomenology of light sterile neutrinos: a brief review20132871110.1142/S0217732313300048HuberP.On the determination of anti-neutrino spectra from nuclear reactors201184216024617HuberP.Reactor neutrinos fluxesProceedings of the 15th International Workshop on Neutrino Telescopes (NeuTel '13)March 2013Venice, ItalyAdeyD.AgarwallaS. K.AnkenbrandtC. M.Neutrinos from stored muons νSTORM: expression of interesthttp://arxiv.org/abs/1305.1419WinterW.Optimization of a very low energy neutrino factory for the disappearance into sterile neutrinos201285119113005BernardiniP.BertolinA.BozzaC.Prospect for charge current neutrino interactions measurements at the CERN-PShttp://arxiv.org/abs/1111.2242RubbiaC.AntonelloM.BaglianiD.A comprehensive search for “anomalies” from neutrino and anti-neutrino oscillations at large mass differences (Δm2 ~ 1 eV2) with two LAr-TPC imaging detectors at different distances from the CERN-PSSPSC-P-345, 2011AntonelloM.BaglianiD.BaibussinovB.Search for “anomalies” from neutrino and anti-neutrino oscillations at Δm2 ~ 1 eV2 with muon spectrometers and large LAr-TPC imaging detectorshttp://arxiv.org/abs/1203.3432CERN Neutrinos Study Group, ICARUS- NESSiE Collaborations, LAGUNA-LBNO ConsortiumLetter of Intent for the new CERN Neutrino Facility (CENF)https://edms.cern.ch/nav/P:CERN-0000077383:V0/P:CERN-0000096728:V0/TAB3RubbiaC.GuglielmiA.PietropaoloF.SalaP.Sterile neutrinos: the necessity for a 5 sigma definitive clarificationhttp://arxiv.org/abs/1304.2047BalsamoE.BergnoliA.BertolinA.The OPERA RPCs front end electronics; a novel application of LVDS line receiver as low cost discriminator201271110.1088/1748-0221/7/11/P11007P11007ZimmermannR.EbertJ.HagnerC.The precision tracker of the OPERA detector20055551-2435450Erratum in Nuclear Instruments and Methods in Physics Research A, vol. 557, p. 690, 200610.1016/j.nima.2005.09.003AdamT.BaussanE.BorerK.CampagneJ.-E.Chon-SenN.de La TailleC.DickN.DracosM.GaudiotG.GoeltzenlichterT.GornushkinY.GraptonJ.-N.GuyonnetJ.-L.HessM.IgersheimR.Janicsko CsathyJ.JolletC.JugetF.KocherH.KrasnoperovA.KrumsteinZ.Martin-ChassardG.MoserU.NozdrinA.OlchevskiA.PorokhovoiS.RauxL.SadovskiA.SchulerJ.SchützH.-U.SchwabC.SmolnikovA.Van BeekG.VilainP.WälchliT.WilquetG.WurtzJ.The OPERA experiment target tracker200757735235392-s2.0-3425084957810.1016/j.nima.2007.04.147AcquafreddaR.AdamT.AgafonovaN.The OPERA experiment in the CERN to Gran Sasso neutrino beam20094410.1088/1748-0221/4/04/P04018P04018BergnoliA.BorsatoE.BrugneraR.BuccheriE.CandelaA.CarraraE.CiesielskiR.CorradiG.D'InceccoM.Dal CorsoF.Degli EspostiL.Di GiovanniA.Di MarcoN.Di TroiaC.DusiniS.FaninC.FeliciG.GambararaA.GarfagniniA.GorinA.GriantiF.GustavinoC.LindozziM.LonghinA.MastropietroF.MengucciA.MonacelliP.PaoloniA.PaoluzziG.PapalinoG.SpinettiM.StancoL.StipcevicM.SugonyaevV.TatananniE.TerminielloL.TerranovaF.VenturaM.VotanoL.ZaunerB.ZaytchenkoA.Tests of OPERA RPC detectors2005526296329702-s2.0-3314448505110.1109/TNS.2005.862902PontecorvoB.Neutrino experiments and the question of leptonic-charge conservation196826984MakiZ.NakagawaM.SakataS.Remarks on the unified model of elementary particles196228587088010.1143/PTP.28.870AndreopoulosC.BellA.BhattacharyaD.CavannaF.DobsonJ.DytmanS.GallagherH.GuzowskiP.HatcherR.KehayiasP.MeregagliaA.NaplesD.PearceG.RubbiaA.WhalleyM.YangT.The GENIE neutrino Monte Carlo generator20106141871042-s2.0-7644910826310.1016/j.nima.2009.12.009HuberP.KoppJ.LindnerM.RolinecM.WinterW.New features in the simulation of neutrino oscillation experiments with GLoBES 3.0. (General Long Baseline Experiment Simulator)200717754324382-s2.0-3444751396210.1016/j.cpc.2007.05.004AgafonovaN.AleksandrovA.AltinokO.Study of neutrino interactions with the electronic detectors of the OPERA experiment20111305305110.1088/1367-2630/13/5/053051