The Results of MINOS and the Future with MINOS+

The Results of MINOS and the Future with MINOS+

A.Timmons,111For the MINOS collaboration. University of Manchester,
Department of Physics and Astronomy, Oxford Road, Manchester, M13 9PL, United Kingdom

The MINOS experiment took data from 2005 up until 2012. This was superseded by MINOS+, the continuation of the two-detector, on-axis, long-baseline experiment based at Fermilab, and at the Soudan Underground Laboratory in northern Minnesota. By searching for the deficit of muon neutrinos at the Far Detector, MINOS/MINOS+ is sensitive to the atmospheric neutrino oscillation parameters and . By using the full MINOS data set looking at both disappearance and appearance in both neutrino and anti-neutrino configurations at the NuMI beam along with atmospheric neutrino data recorded at the FD, MINOS has made the most precise measurement of . Using a full three-flavour framework and searching for appearance MINOS/MINOS+ gains sensitivity to , the mass hierarchy and the octant of . Exotic phenomenon is also explored with the MINOS detectors looking for non-standard interactions and sterile neutrinos. The current MINOS+ era goals are to build on the previous MINOS results improving the precision on the three-flavour oscillation parameter measurements and strengthening the constraints placed on the sterile neutrino parameter space.

1 Introduction

Over the last couple of decades physicists across the world have obtained model independent evidence for neutrino oscillations. It was in 1998 that Super-Kamiokade ref:SuperK () observed muon neutrinos changing flavour as they transversed the atmosphere. Later in 2001 the Sudbury Neutrino Observatory experiment ref:SNO () observed neutrinos oscillating between flavours which originated from the sun. Evidence for reactor anti-neutrinos was seen in 2002 with the KamLAND ref:KamLAND () experiment. There are now multiple generations of experiments designed to confirm and probe the nature of neutrino oscillations using a neutrino source from an accelerator ref:OPERA (); ref:T2K (), solar neutrinos ref:SAGE (); ref:GALLEX (), and nuclear reactor anti-neutrinos ref:Chooz (); ref:Reno (); ref:DoubleChooz (); ref:DayaBay (). The implication of neutrino flavour change is indicative that the neutrino must have a non zero mass and violate lepton number conservation; a clear observation of new physics beyond the standard model.

It was during this time of discovery that the MINOS ref:MINOSPro () experiment was proposed. MINOS was designed with a long-baseline and two detectors 1.04 km and 735 km from the neutrino production target respectively; with the goal to measure the atmospheric neutrino oscillation parameters. A previously constructed long-baseline, two detector experiment K2K ref:K2K () in Japan had the same goal. MINOS was unique in that its magnetised detectors allowed one to distinguish and interactions on an event by event basis. The MINOS detectors encountered a higher flux of neutrinos than K2K and over the next few years MINOS contributed to the the era of precision measurements of the fundamental parameters governing this quantum mechanical effect of neutrino oscillations. The MINOS experiment helped show how effectively a two-detector experiment can minimise the large systematic uncertainties associated with a neutrino interaction experiment. In 2012 the MINOS experiment ended and it is the continuation of the detectors taking data in the upgraded accelerator for the NOref:Nova () experiment that the MINOS+ ref:MINOSPlusProposal () experiment was born. MINOS+ began taking data in September 2013, with a higher flux of neutrinos at high energies MINOS+ becomes sensitive to beyond the standard model neutrinos physics such as sterile neutrinos and large extra dimensions.

The theory of neutrino oscillation describes the change in neutrino flavour composition seen in data. These oscillations arise due to a mixture between mass and flavour eigenstates; three active flavours of neutrino and three mass eigenstates are required to fully describe the neutrino oscillations observed in data. The energy dependance of these oscillations are governed by the difference of the square of the mass eigenstates, and , while the degree of mixing (the amplitude of the oscillations) is governed by three mixing angles , and a violating phase . These parameters make up a rotation matrix known as the PMNS rotation matrix ref:PMNS ().

The difference between the two mass splittings is almost two orders of magnitude, has been shown to be small by measurements of this parameter by reactors experiments ref:DayaBay (); ref:Renoth13 (); ref:DoubleChoozth13 () and so one can approximately decouple the two frequencies into two distinct regimes: the “solar” oscillation regime is driven by and and mostly determines the flavour composition of particles propagating from within the sun towards Earth. The “atmospheric” oscillation regime is driven primarily by and and their values govern neutrino oscillation observed in neutrinos decaying from secondary hadrons due to cosmic rays interacting within the Earth’s atmosphere. The MINOS experiment was designed to probe the atmospheric sector by using neutrinos from a man made source. MINOS also has sensitivity to the parameter through observing and appearance, allowing for a full analysis combining both disappearance and appearance data.

1.1 Oscillation Physics at MINOS

The neutrino oscillation regime MINOS is most sensitive to is driven by the larger of the two mass-splitting differences ; consequently a two-flavour approximation can be used to describe the data using a single mass splitting and effective mixing angle . Using this approximation one can express the muon neutrino survival probability as


where is the neutrino propagation distance and is the neutrino energy. Previous analyses by MINOS rely on this two-favour approximation. However, the neutrino community is entering an era of precision measurement; the error on the atmospheric mass splitting is down to a few percent level and with the discovery of a non zero in 2012 by Daya bay ref:DayaBay (), has now become one of the most precisely measured angles, and so the need to move to a fuller treatment of neutrino oscillations is ever present.

Within a three-flavour framework, the oscillations are driven by the mass splittings and , where . For exact calculations of the oscillation probabilities one must consider all parameters due to interferences. One can modify the two-flavour oscillation probability parameters in Eq.(1) as follows:


Only by moving to a three-flavour framework can the degeneracies between the octant of and determination of the mass hierarchy (sign of ) be broken. These equations only account for neutrino oscillation within a vacuum, however, when neutrinos traverse through matter the eigenstates become modified due to the MSW effect ref:MSW1 (); ref:MSW2 (). To account for this one can replace with a modified mixing angle , given by ref:thetamod (), such that


where the magnitude of determines the size of the matter effect and can be expressed as , where is the Fermi weak coupling constant and is the density of electrons in the medium. The sign of A is positive (negative) for neutrinos (anti-neutrinos). The magnitude of in equation 4 influences the amount of mixing between . This MSW mechanism has an effect on a disappearance analysis through mixing, therefore, in order to perform a precision measurement using disappearance one must also take into account appearance.

MINOS is also sensitive to interactions, the probability for appearance to second order can be approximated to ref:MatterMod ():


where , and . The second term in equation 5 will have a plus (minus) sign for neutrinos (anti-neutrinos). Being sensitive to and appearance allows one to probe the mass hierarchy (as the sign of will change) as well as the CP violating phase.

2 The MINOS experiment

2.1 The NuMI Beam

The Neutrinos at the Main Injector (NuMI) neutrino beam ref:NUMI () was built at Fermilab, to provide neutrinos for the MINOS experiment. The NuMI beam typically has a beam power of 350 kW with a design specification of up to 400 kW. A high-intensity beam is required to achieve a meaningful event rate at the MINOS Far Detector (FD) placed several hundred kilometres away. Such a distance significantly reduces the neutrino flux as it falls with the square of the distance from the decay point.

Figure 1: The NuMI beam.

To produce such a powerful neutrino beam the first stage is to create protons from H- ions. To achieve this the ions are accelerated by a Radio Frequency (RF) quadrupole up to an energy of 750 keV. From there a linear accelerator then accelerates the ions to an energy of 400 MeV which have subsequently passed through a thin carbon foil stripping the electrons off the ions to leave a beam of protons. The protons are fed into a rapid cycling synchrotron (Booster) and accelerated in batches up to energies of 8 GeV. Subsequently they are fed into the Main Injector where they are accelerated to 120 GeV. The Main Injector has a circumference seven times larger than the Booster and so up to six batches can be inserted into the Main Injector at once. It is the interaction between these high energy protons upon a fixed graphite target which results in plethora of charged hadrons (predominantly pions, with a significant kaon component at higher energies). These charged hadrons pass through two parabolic, magnetic horns which focus either positive or negative hadrons depending on the direction of the electric current being pulsed through the horns. The focused hadrons travel along a 675 m decay pipe. Its the decay of these hadrons within the pipe that form the predominately muon flavour neutrino beam. By focusing the positive hadrons a beam of predominately is created (-dominated beam mode), by focusing the negativley charged hadrons the component can be increased (-enhanced beam mode) . Figure 1 shows a diagram of the charged hadrons being focused by the two horns into the decay pipe.

Figure 2: The reconstructed neutrino energy spectra at the MINOS Near Detector. The top picture shows the energy spectrum for positively focused hadrons producing a predominately beam. Below is when negatively charged hadron are focused which increases the amount of events seen at the detectors. Note how in anti-neutrino mode the event rate is significantly less due to different cross sections between neutrinos and anti-neutrinos.

Figure 2 shows the composition of the NuMI beam for charged current neutrino interactions observed in the MINOS Near Detector. The significant difference in composition and event rate between these beam modes arises mainly from the fact that the interaction cross section is a factor of approximately two lower than the interaction cross section.

Figure 3: The true energy distribution of near detector events from the three different NuMI configurations. The target position (distance upstream of a nominal position) and horn current of each configuration are shown in the caption. “0.1 m, 185 kA" is the “low-energy" configuration in which most of the MINOS data has been taken. “1.0 m, 200 kA" is the “medium-energy" configuration which is the beam configuration in the NOA era and therefore is the beam setup in which MINOS+ is currently taking data in. The final configuration is the “pseudo-high-energy" configuration.

The neutrino energy spectrum provided by the NuMI beam is tunable, through changing the relative positions between the target and the focusing horns. Three of the possible configurations are shown in figure 3. The main goal of MINOS was to measure the atmospheric oscillation parameters, this would require a large flux at the oscillation dip which would be located around the 2 GeV region. In the MINOS era the configuration of the NuMI was set to the low energy setting. For the MINOS+ era the NOA experiment is mrad off axis and requires the NuMI beam to be set at the medium energy configuration. MINOS/MINOS+ is an on axis experiment (the detectors line up with the beam axis) and therefore MINOS+ observes a high flux of neutrinos at higher energies compared to the low energy configuration.

Figure 4: Showing the beam configuration for the MINOS and MINOS+ experiments. Most of the data for MINOS was taken in the low energy -beam (in green) and (orange). Special runs where the beam was configured to higher energies or the magnetic horn turn off is indicated by the red. The Magenta coloured runs are the beginning of the MINOS+ data taking in the -beam during the NOA era.

2.2 The MINOS Detectors

The MINOS experiment has two steel-scintillator calorimeters ref:MINOSNIM () designed with the same materials and to operate in an identical manner, known as functionally-identical. The calorimeters measure the energy deposition and event topologies of neutrino interaction events. The detectors are shown in figure 6. Both detectors are made of alternating layers of 1.00 cm thick plastic-scintillator and 2.54 cm thick steel planes. As neutrinos travel through the detector they interact with the iron-nuclei, and the charged final-state particles travel through the scintillator depositing energy which is read out as light.

Figure 5: A strip of MINOS scintillator with the WLS fiber installed along the center. The image in the lower right corner depicts a charged particle (red line) depositing energy by producing light (blue line) as it passes through the strip. The light travels along the WLS fiber and is read out by a PMT.

The light travels along WaveLength Shifting (WLS) fiber and is read out by a series of a PhotoMultiplier Tubes (PMTs). Figure 5 shows a strip of the scintillator used in the MINOS detectors, along the middle a groove is made so that a wavelength shifting fiber can be installed. It is from the light patterns that a neutrino candidate event can be reconstructed, so that information about the topology of the event can be extracted. The steel planes are magnetised by a coil aligned to the longitudinal axes of each detector giving a magnetic field of approximately 3 T. The trajectories of the charged particles are therefore curved (the direction depends on the polarity of the current in the coil) and thus and CC interactions can be distinguished.

Figure 6: The MINOS detectors. Left: the Near Detector at Fermilab; right: the Far Detector at the Soudan Underground Laboratory.

The Near Detector (ND) is situated 1.04 km downstream from the neutrino target at Fermilab. With a mass of 0.98 kton the ND measures the reconstructed neutrino energy spectrum before oscillations have occurred. The ND has two distinct sections so that it can take advantage of the high neutrino flux at this location to define a relatively small target fiducial volume for selection of events for the near/far comparison. The section closest to the target is used to define the interaction vertex and measure the energy of the neutrino-induced hadronic shower; every plate is instrumented with plastic scintillator to act as a calorimeter. The second section is used as a muon spectrometer to measure the momenta of energetic muons where one in every five plates is instrumented with scintillator. The scintillator planes are made up of 4 cm-wide strips. The strips on adjacent planes are oriented perpendicular to each other to allow three-dimensional reconstruction of events. The planes are oriented to the vertical defining a co-ordinate system referred to as the and directions.

The Far Detector is 735 km downstream from the neutrino production target, 705 m underground in a mineshaft in northern Minnesota. The FD is significantly larger than the ND to compensate for the decrease in the neutrino flux. With a mass of 5.4 kton the FD measures the reconstructed neutrino energy spectrum and will observe a different neutrino flavour composition of the beam due to neutrino oscillation. The geometry is similar to that of the ND in that it is spilt into two “super modules" of 239 and 247 planes, however the FD does not have a vertex and spectrometer section due to its distance away from the beam, it observes significantly less neutrino interactions. A veto shield composed of layers of scintillator covers the top and sides of the FD to better identify incoming cosmic ray muons that may enter the fiducial volume of the detector helping to obtain a high pure sample of downward-going atmospheric events.

The two detector method is a very powerful experimental setup, since it allows the cancellation of large uncertainties that beset any neutrino oscillation experiment. The uncertainty in neutrino flux and cross sections are only known to tens of percent. Therefore, by looking at the disappearance and appearance in the FD relative to the ND the uncertainties can be reduced significantly as only the relative uncertainty between the two detectors will affect the final measured result.

3 Neutrino Interactions in the MINOS Detectors

There are three neutrino interactions that are of interest to MINOS as shown in figure 7.

Figure 7: Neutrino interaction topologies observed in the MINOS detectors. Left: A charged current interaction. Middle: A neutral current interaction. Right: A charged current interaction. Each coloured pixel represents a scintillator strip with energy deposited from a charged particle. The colour scale displays the amount of light: purple and blue are low light levels, through to orange and red for the highest light levels.

The main channel is the charged current (CC) () interaction

The cascade of hadrons, , produces a diffuse shower of energy deposits near the interaction vertex. MINOS was constructed with steel planes so that it can contain a signficant proportion of the final state muons. A muon produces a long track that curves due to the magnetic field. It is the direction of curvature that allows MINOS to identify the incoming neutrino as a or a .

All active neutrino flavours undergo the neutral current (NC) interaction through the process

Only the hadronic shower is observed, producing a diffuse pattern of energy deposits within the detector. It is not possible to determine the flavour of neutrino.

Finally, electron neutrinos undergo CC interactions through the process

The electron gives rise to an electromagnetic shower, which produces a much denser, more compact shower of energy deposits. This interaction is difficult to detect due to the steel plate thickness, meaning only a few events leave any energy deposits in the plastic-scintillator.

Above a few GeV the dominant process is deep inelastic scattering (DIS). Here the neutrino has sufficient energy that it can resolve the individual quark constituents of the nucleon which manifests in the creation of a hadronic shower. However, at MINOS the oscillation dip observed in muon neutrino disappearance at the far detector occurs just below GeV; these neutrino interactions provide a large source of signal events for a neutrino oscillation analysis. At this energy neutrinos can elastically scatter off an entire nucleon liberating a nucleon (or multiple nucleons) from the target. In the case of charged current scattering this process is referred to as “quasi-elastic scattering" (CCQE). A detailed review on the current state of neutrino cross sections can be found here ref:NuCrossSectionReview ().

The energy of a neutrino event is calculated by summing the shower energy deposits and the muon track energy. The energy resolution of contained muon tracks is 4.6% ref:MuonTrk1 (). If a muon track exits either one of MINOS detectors then the curvature of the track is used to calculate the energy. The curvature of a muon track is directly proportional to the ratio of its electric charge to its momentum; at the peak of the neutrino beam at around 3 GeV the resolution is 11% ref:MuonTrk1 ().

All three interaction processes can result in a shower of energy deposited in the detectors. The MINOS detectors are too coarse to reliably use shower topology information and so the energy is reconstructed calorimetricaly. The final calorimetric hadronic and electromagnetic shower resolutions are well modelled by simulation and the resolution is parametrised as for hadronic showers and for electromagnetic showers, where E is the particle energy in GeV ref:MikeThesis (); ref:TriciaThesis ().

A calorimeter response is different for hadronic and electromagnetic showers. A fraction of the energy deposited by the showering particle is invisible, i.e. it does not contribute to the calorimeter signal, this can cause some undesirable properties in a calorimeter causing non-linearities. For CC () interactions a more sophisticated method is implemented to measure shower energies ref:BackhouseThesis (). A -nearest-neighbour algorithm ref:KNNref () uses a multivariate analysis of a broader range of event-level information to provide an estimate of shower energy. The variables are by the algorithm: the event length, the average energy deposited per scintillator plane along the track, the transverse energy deposition, and the fluctuation of the energy deposition along the track. The -nearest-neighbour algorithm uses the 400 nearest neighbours and comparing them to a simulated library of events the total shower energy of a event can be estimated. Figure 8 shows the distribution of reconstructed over true shower energy for different ranges of true shower hadronic shower energy using calorimetric energy and -nearest-neighbour algorithm estimator. It can be seen that the -nearest-neighbour algorithm gives a better estimate of the shower energy at lower energies. This improves the hadronic energy resolution of from 55% to  43% for showers between 1.0 - 1.5 GeV.

Figure 8: Distribution of reconstructed over true shower energy for different ranges of true shower energy (indicated in the corner of each panel). The calorimetric estimator is shown in black, and the -nearest-neighbour estimator in red. In each panel, the dashed line indicates the position of . Figure taken from reference ref:BackhouseThesis ().

3.1 Selection of charged-current and interactions

MINOS was designed to measure the neutrino oscillation parameters in the atmospheric region ( and ). To achieve this one needs a sample of () CC events. Figure 7 shows that this can be achieved by selecting neutrino candidate events with a muon track.

There are three main backgrounds that have an effect on signal purity. At low energies, NC interactions can result in a charged hadron producing a track, thus being mistaken for a -CC event with a low energy muon track. Atmospheric events represent a potential source of muon neutrinos distinct from those in the muon beam and hence are not useful when only considering a beam disappearance oscillation analysis. Wrong-sign events can become a background when the the muon charge deduced from curvature is measured incorrectly if one wishes to separate and samples to test if they oscillate with different probabilities.

To reduce the NC background one needs to separate the NC and CC candidate events. MINOS uses a -nearest-neighbour algorithm; a simulated high statistics data set is created with two known classes of events, one with a muon track and one without. Four variables are used to create a discriminating variable to be applied to all track-like events. The variables are: the number of MINOS detector planes associated with a muon track (muon tracks tend to extend much further than NC showers), the average energy deposited per scintillator plane along the track, the transverse energy deposition profile, and the variation of the energy deposited along the muon track. Figure 9 shows the distribution of the -nearest-neighbour algorithm as a single variable ref:RustemThesis ().

Figure 9: The -nearest-neighbour discrimination variable used to separate CC interactions from track-like hadronic backgrounds. Events with a parameter value greater than 0.3 are selected as CC interactions for analysis.

The event selection is identical for and CC events due to their similar interaction topologies. The NuMI beam can be configured to produce an anti-neutrino enhanced beam. From the curvature of the muon track reconstructed by a Kalman Filter ref:JohnMarshallThesis () algorithm the lepton number of the neutrino can be deduced. After track identification all the remaining hits which are in proximity to one another are grouped into showers.

3.2 Selection of charged-current interactions

The selection for events relies on looking for events with a dense shower arising from electromagnetic interactions from an electron as it passes through the MINOS detectors. The dominant background comes from NC events generating a dense hadronic shower. Such behaviour can be caused by a neutral pion decaying into a pair of photons. The majority of appearance expected at the FD occurs in neutrino events with energy in the range 1-8 GeV and so only this range is considered.

Figure 10: Showing an example of a simulated event compared to a good and bad match from the LEM library.

The granularity of the MINOS detectors makes resolving any topology from electromagnetic showers almost impossible, therefore all candidate events with a shower have their energy deposition patterns compared to a large library of order simulated events containing signal (40%) and background events (60%). This technique is called Library Event Matching (LEM) ref:PedroThesis (); ref:RuthThesis (); ref:AdamSThesis (). The 50 simulated events that match the event with a similar pattern of energy deposited in each scintillator strip excited by the shower are chosen. Data and library events are not spatially translated to align them for the best matching. The matching procedure provides a quantitative means of determining the likelihood that two different charge topologies were created by the same primary deposition. For an arbitrary energy deposit, the mean expected charge on a photomultiplier tube will be some value . Consider strip i in the plane of the detector and events A and B, where the detector’s response in event A was photoelectrons and the response in event B was photoelectrons. The likelihood , of a data event corresponding to the same physical shower topology as a simulated library event can therefore be calculated as


An example of a good and bad match from LEM can be seen in figure 10.

Figure 11: Distributions of . The plots in the left column correspond to the neutrino dominated beam mode. The plots in the right column correspond to the anti-neutrino enhanced beam mode. The top row shows the distributions for ND selected events with a band about the simulation representing the systematic uncertainty. The bottom row shows the distributions for the predicted FD background and signal multiplied by 10

Three variables are constructed from the 50 best simulated events (signal or background), these are: the fraction of the events that are true CC events, the average inelasticity (this is the amount of energy that goes into the hadronic shower) of the true CC events, and the average fraction of charge that overlaps between the data event and each CC library event. These three variables along with the reconstructed energy of the data event, are fed into a neural network which calculates a classification of how signal-like the data event is. A single variable, , is formed to quantify this, as shown in figure 11. it’s an output of an artificial neural network with several variables coming from the event comparisons. Events with are selected for analysis, this number was optimised to maximise the sensitivity to appearance ref:MINOSNuE2012 (). Candidate -CC and -CC events are required to fall within a fiducial volume and to be coincident in time () and direction with the NuMI beam. Events are required to have shower-like topologies by rejecting events with tracks that are longer than 25 planes or extend more than 15 planes from a shower edge.

With the absence of a -CC and -CC signal in the ND, the signal-selection efficiency can not be extrapolated from the ND events in the same way as the background estimate. By using real data, well-identified -CC events are selected. By removing the energy deposited by the muon track ref:AnnaThesis () these events can be used to calculated the signal selection efficiency by inserting the energy deposition from an electron with identical momentum to that of the removed muon. This allows one to effectively convert a well identified sample of -CC and -CC data events into a sample of -CC and -CC data events. Using this method, the -CC identification efficiency is found to be in the neutrino-dominated beam, and in the antineutrino-enhanced beam.

3.3 Selection of neutral-current interactions

The signal for a NC event is a diffuse hadronic shower, however -CC interactions can also have large hadronic showers, if the inelasticity of the event is high ( most of the energy is given to the shower). The muon track may then be difficult to detect if it does not significantly extend beyond the hadronic shower. To achieve a high purity of NC events a number of selection cuts are used ref:GemmaThesis (). An event is classified as a NC event if it has no tracks or if a track does not extend more than six planes past the end of the shower. The NC identification efficiency is 89%, with 61% purity, this is determined by taking an average over the energy spectrum for simulated events. However, this selection will identify 97% of CC interactions as NC events, which needs to be taken into account when searching for NC disappearance at the FD.

3.4 Selection of atmospheric neutrinos

Figure 12: Distribution of reconstructed zenith angle for muons with good timings and topology. In the range the observed rate of muons is dominated by the cosmic-ray background and falls steeply as the mean rock overburden increases. To minimise the background from cosmic-ray muons, events for analysis are required to satisfy .

Atmospheric neutrinos are selected as -CC events in the MINOS FD outside of the 10  window period when the NuMI beam is producing neutrinos ref:MINOSAtmosPRD (). The atmospheric neutrino signal is separated from the cosmic-ray background using two characteristic signatures of atmospheric neutrino interactions: either a reconstructed vertex inside the fiducial volume; or a reconstructed upward-going or horizontal muon trajectory.

The FD timing resolution on a hit by hit basis is 2.5 , which is enough to calculate the direction of a muon track inside the FD. For upward and horizontal angles, where the rock overburden exceeds 14,000 m water-equivalent, the absorption of cosmic-ray muons by the earth is sufficiently high that the observed flux of muons is dominated by atmospheric muon neutrino interactions ref:atmosOverburden (). At the Soudan mine (the location of the FD), upward-going tracks with a zenith angle are defined as upward-going and horizontal ref:atmosSoudan (). Therefore, upward-going and horizontal provide a signature for atmospheric neutrinos. To further reduce the background the analysis requires as shown in figure 12.

For tracks where the end point lies inside the fiducial volume, the muon momentum is reconstructed from the measured track length; for exiting tracks, the momentum is obtained from the fitted track curvature. In both cases, the fitted curvature is used to determine the muon charge sign.

4 Muon neutrino and antineutrino disappearance

MINOS can measure the atmospheric neutrino oscillation parameters and by fitting the energy dependance of -CC and -CC disappearance. The FD sees the neutrino beam as a point source whereas the ND subtends a relatively large angle to the beam. Thus once a neutrino parent decays the ND would see a large spread in energies compared to the FD which preferentially selects higher energy neutrinos at a smaller angle from the direction of the parent hadron. MINOS pioneered the technique of resolving this kinematic issue by constructing a beam matrix that allows one to convert an energy spectrum observed in the ND to the corresponding spectrum observed at the FD ref:JustinThesis (). By using the ND data as a constraint many systematics which affect both detectors cancel to first order. Only systematics affecting both detectors differently become a significant source of uncertainty; primarily reconstruction efficiencies and miscalibrations of the neutrino energy measurement in the detectors  ref:StephenThesis (). Figure 13 shows the reconstructed energy spectrum at the FD for and neutrino events compared to two predictions, if there were no oscillations and a best fit to the FD data.

Figure 13: The FD reconstructed neutrino energy spectrum for (left) and (right). The light grey light shows the predicted FD energy spectrum if neutrinos did not oscillate. The red line indicates the best fit to the data. These energy spectra were fit along with with the atmospheric and samples discussed in the later subsections. For details on the fit and the result of the neutrino oscillation parameters obtained see section 6.

There is an uncertainty on the relative normalisation of the selected Near and Far Detector event samples which is dominated by differences in the reconstruction and selection efficiencies between the two detectors, as well as relative uncertainties on fiducial mass and live time; this uncertainty is found to be 1.6%. There are two different uncertainties on the measurement of hadronic shower energy ref:BackhouseThesis (), these are the relative mis-modelling of the energy scale between the two detectors as well as the absolute mis-modelling. It was found that the relative uncertainty for the ND is 1.9% and the FD is 1.1%. The absolute mis-modelling comes from the uncertainties on the modelling of hadronic showers which is fully correlated bin to bin in reconstructed energy and has an energy dependence of the form  ref:shower ().

A previous two-flavour analysis ref:MINOS2flav2013 () of and disappearance using the combined accelerator and atmospheric data from MINOS yielded and . A symmetry in CPT requires that neutrinos and anti-neutrinos oscillate in an identical way, thus their oscillation parameters should be identical. With the ability to distinguish the lepton number of a neutrino MINOS can measure the oscillation parameters for anti-neutrinos and neutrinos separately. Using both atmospheric and beam anti-neutrinos MINOS measures the oscillation parameters to be and  ref:MINOS2flav2013 () which is in good agreement with the parameters measured from oscillations as shown in figure 14.

Figure 14: The allowed region for antineutrino oscillation parameters (black line), compared to the region measured with neutrinos alone (red line) and the region measured using both neutrinos and antineutrinos under the assumption they have the same parameters (blue line).

4.1 Atmospheric and disappearance

Figure 15: The distribution of reconstructed atmospheric and CC interactions over various energy ranges observed at the FD, compared to the expectation with and without neutrino oscillation. Two different fits are performed for the normal and inverted hierarchy of . The energy of the neutrino event is the sum of muon momentum (from range and/or curvature) and total shower energy.

MINOS is the first experiment to probe the resonance predicted to occur in multi-GeV, upward-going atmospheric neutrinos which travel through the earth’s mantle for both neutrino and anti-neutrino events on an event by event basis. The atmospheric events are separated into samples of contained-vertex and non-fiducial muons for neutrinos and anti-neutrinos, figure 15 shows the atmospheric samples containing events with a contained-vertex. By measuring the -CC and -CC interactions separately this allows MINOS to gain sensitivity to the mass hierarchy and octant. The difference between the inverted and normal mass hierarchy in figure 15 is very marginal.

The neutrino events are binned as a function of and , where E is the reconstructed energy of the event in GeV and is the zenith angle of the muon track, this binning gives enhanced sensitivity to the MSW resonance. A sample of contained-vertex showers are also selected from the data, composed mainly of NC, -CC and -CC interactions. They are used to constrain the overall flux normalization. For atmospheric neutrinos, the earth is modelled by four layers of constant electron density using the PREM model ref:PREM (). Comparisons to a more detailed 52 layer model yielded very similar results and so the extra computational time was avoided by using the simple four layer model.

5 and appearance

By searching for and appearance at the FD MINOS can perform a measurement of . The POT exposure for this data set is protons-on-target (POT) using a -dominated beam and POT using a -enhanced beam. Neutrino events with are selected for analysis in the -domainted beam mode and in the -enhanced beam. Neutrino events with are considered background-like and therefore insensitive to and appearance. The background consists of three components: NC interactions, CC- and interactions, and the intrinsic component in the beam. The relative contribution between the ND and FD is different for all of these components, since they are affected differently by oscillation, and the kinematics of the production in the beam are different. Each background must be individually measured. By changing the configuration of the NuMI beam (low, medium or pseudo-high) one can measure these backgrounds. The relative contributions of the background components changes in a well understood way ref:JoaoThesis (). Neutrino events with are used to provide validation to the analysis procedure; MINOS uses ND neutrino events with to predict FD event yields. Comparing these yields to that found in the FD agree to within () of the statistical uncertainty for the data sample in the () beam modes.

Figure 16: The reconstructed neutrino energy distributions for three ranges. The events with energy greater than 5 GeV are combined into a single bin for the fits. The vertical bars through the data points denote statistical uncertainties. The signal prediction assumes , (normal hierarchy), and (maximal mixing). The left corresponds to data taken in -dominated beam mode samples and the right is from -enhanced beam mode samples.

For an appearance analysis one has to consider a three-flavour neutrino oscillation probability that includes matter effects. The fit is done simultaneously for data from both -dominated beam mode and -enhanced beam mode samples. A total of 127.7 background events are expected at the FD in the neutrino-dominated beam, and 17.5 events in the antineutrino-enhanced beam. In the data, 152 and 20 events are observed, respectively and their energy distributions can be seen in figure 16.

MINOS cannot distinguish and events directly, however, the relative number of neutrino and antineutrino interactions in the neutrino-dominated and anti-neutrino-enhanced beams is well known. The parameter measured is while the values and remain fixed. The fit is ran over both hierarchies and all possible values of . For more details on this analysis see reference ref:MINOSNuE2012 () . The fit is performed using the 15 bins formed by three bins of and five bins of energy as seen in figure 16, note final three bins in energy from 5-8 GeV were merged for the fit.

Figure 17: The allowed regions for .

MINOS finds that the data allow for a value of for the normal hierarchy with and for the inverted hierarchy with , in both cases , this can be seen in figure 17. This the first ever appearance search in a long-baseline beam.

6 A Combined Three Flavour Analysis

To gain maximum sensitivity to the mass hierarchy and , one needs to perfrom a full three flavour analysis combining and disappearance with and appearance. For this one needs to account for both disappearance and appearance of neutrinos oscillations at both detectors.

Figure 18: The left panels show the 68% and 90% confidence limits on and for the normal hierarchy (top) and inverted hierarchy (bottom). At each point in this parameter space, the likelihood function is maximised with respect to , and the systematic parameters constraints in the fit. The likelihood surface is calculated relative to the overall best fit, which is indicated by the star. The right panels show the 1D likelihood profiles as a function of and for each hierarchy. The horizontal dotted lines indicate the 68% and 90% confidence limits.

For a full three-flavour fit MINOS uses an accelerator neutrino data set comprising exposures of protons-on-target (POT) using a -dominated beam and POT using a -enhanced beam. Both sets were acquired in the low energy NuMI beam configuration. MINOS also collected kt-years of atmosphere neutrino data.

The oscillation parameters are determined by applying a maximum likelihood fit to the data. The parameters , , and are varied in the fit with an external constraint on the mixing angle calculated from a weighted average of the latest published results from the reactor experiments Daya Bay ref:DayaBay (), RENO ref:Reno () and Double Chooz ref:DoubleChooz (). The constraints are included by adding a Gaussian prior penalty term onto the likelihood during the fit. The solar parameters are kept at the fixed values of and  ref:Fogli (). To test the impact of the solar parameters on the fit they were varied within their uncertainties and the effect was negligible on the final results.

Figure 19: The 1D likelihood profile as a function of for each combination of mass hierarchy and octant. For each value of the likelihood was maximised with respect to , and .

Figure 18 shows the 2D confidence limits on , as well as the 1D profiled confidence limits on and separately. One can see the sensitivity to the octant of through the inclusion of the atmospheric neutrino sample. The 68% (90%) confidence limits (C.L.) on these parameters are calculated by taking the range of negative log-likelihood values with < 1.00 (2.71) relative to the overall best fit. This yields at 68% C.L and at 90% C.L in the normal hierarchy; at 68% C.L and at 90% C.L in the inverted hierarchy. The data disfavour maximal mixing by < 1.54. These results give the most precise measurement made on to date.

Figure 19 shows the 1D C.L of the profiled likelihood surface for value of for each of the four possible combinations (hierarchy and octant of ). The data disfavour 36% (11%) of the parameter space defined by , the octant, and the mass hierarchy at 68% (90%) C.L.

7 NSI Interactions

Non-standard interactions (NSI) ref:NSI1 (); ref:NSI2 (); ref:NSI3 () could occur between muon or tau neutrinos and matter, and could alter the flavour composition of a neutrino beam as it propagates through the earth’s crust. Searches for NSI have already been performed with atmospheric neutrinos ref:NSIsuperK (). Using a two-flavour approach one can write the probability for muon neutrino survival as:

Figure 20: 2D contours for both 68% and 90%, showing the allowed regions of parameter space for from the fit to the data through muon disappearance.

where is the neutrino path length and is defined as the NSI matter oscillation length defined as:


where .The signs in equation 8 arise from the matter potential, , which is positive for neutrinos and negative for antineutrinos. The parameter is real-valued and carries its own sign. A positive value of implies that the neutrino disappearance probability is greater than the antineutrino disappearance probability, and vice versa.

The results presented here are based on an exposure of protons on target (POT) in neutrino mode, combined with a POT exposure in antineutrino mode. Due to the opposite sign of the matter potential in equation 8 for neutrinos and antineutrinos, NSI, if present, will alter the survival probability of neutrinos and antineutrinos in opposite directions. This analysis has identical event selection as for the disappearance analysis, however the fit takes into account perturbations from the standard three-flavour formalism brought about by NSI. The best fit parameters from this procedure are found to be , and with the allowed region (90% C.L.). The systematic uncertainties incorporated into the penalty terms when maximising the likelihood have a negligible effect on the fit. Figure 20 shows the allowed regions of the fit parameters, these are 2D C.L contours, produced by profiling the likelihood with respect to the other parameters.

This is the first direct search for nonstandard interactions with high-purity samples of both neutrinos and antineutrinos conducting a simultaneous fit to neutrino and antineutrino energy spectra of conventional oscillations with an additional NSI matter effect. This result is consistent with the null hypothesis of no NSI.

8 Sterile Neutrinos

Figure 21: Showing the reconstructed energy spectrum for CC interactions in the FD (top) and the ND (bottom). The grey hatched histogram is the NC background and the red line is the three flavour prediction with systematic uncertainties included.

There have been several anomalous results within the neutrino community that have questioned our understanding of neutrinos. The Liquid Scintillator Neutrino Detector (LSND) and MiniBooNE short-baseline experiments saw an excess in the data of electron antineutrinos that can not be explained using the current three-favour model ref:LSND (); ref:MiniBooNE (); a reinterpretation of reactor fluxes has also lead to a discrepancy in neutrino-oscillation reactor experiments ref:ReactorNeutrinoAnomaly (). One explanation is the addition of one or more neutrino types which would oscillate with the three active neutrino flavours; a comprehensive overview of this explanation to account for the above discrepancies can be found in reference ref:SterileWhitePaper ().

MINOS has sensitivity to sterile oscillation signatures by looking for perturbations from the three-flavour oscillation formalism in CC events and a deficit in NC events. This MINOS analysis ref:MINOSsterile () considers a 3+1 sterile neutrino model. By adding an additional neutrino the PMNS is extended to a matrix which introduces three additional mixing angles , and . This extra mass state also introduces an extra mass splitting, where such that . An additional deficit of muon neutrinos at the FD aside from the expected loss due to three-flavour oscillations would be an indication of interference from these additional parameters.

Figure 22: Showing the reconstructed energy spectrum for neutral current events at the FD (top) and the ND (bottom). The grey hatched histogram is the CC background and the red line is the three-flavour prediction with systematic uncertainties included. The blue dashed line represents an additional three-flavour prediction with a different with value

Figures 21 and 22 show the NC and CC reconstructed neutrino energy spectrum for the FD and ND and are in good agreement with the expectation of a null sterile neutrino hypothesis. This agreement can be quantified by using the test statistic R for the number of NC events observed at the FD:


where is the number of events observed, is the predicted number of CC background interactions, and is the predicted number of NC interactions in the detector. A value of would indicate no mixing in the data with sterile neutrinos. The test statistic R is simply based on the integrated number of events, over the full energy range MINOS obtains (stat) (syst) which is in good agreement with the null hypothesis ref:MINOSsterile ().

Figure 23: Showing the C.L.s for the three mixing angles , , . In the fit the sterile mass splitting was kept fxed at . The black and blue contours are for a different fixed value of of ad respectively.

MINOS is insensitive to the mixing angle which is primarily involved in appearance mixing; by looking at muon disappearance MINOS has set C.L.s limits on the other sterile mixing angles yielding , and  ref:MINOSsterile (). Figure 23 shows the contours for the mixing angles , , at a particular value of .

MINOS sets a limit on the sterile-active neutrino coupling by constructing a quantity , the fraction of that have oscillated into , expressed as:


For neutrino events around the oscilltion maximum with an energy of 1.4 GeV (the energy for the highest probability of muon neutrino disappearance) a large number of test values are (randomly sampled from Gaussian distributions from the sterile ) selected for the mixing angles , , and . The value of that is larger than 90% of the test cases is used as the limit, which yields at a 90% C.L ref:MINOSsterile ().

9 The Future with MINOS+

MINOS+ ref:MINOSPlusProposal () is the continuation of the MINOS detectors taking data during the NuMI beam run in the medium energy configuration. Being on-axis the neutrino flux at the detectors significantly increases at higher energies as shown in figure 24. The beam peak in the medium configuration shifts from 3 GeV to 7 GeV allowing MINOS+ to observe around 4,000 -CC interactions in the FD each year. MINOS+ has been taking data since September 2013, and with the additional statistics will provide a useful contribution to the high-precision test of the three-flavour oscillation formalism and will improve on the world-leading measurements of muon neutrino disappearance made by MINOS.

Figure 24: The energy spectrum observed for the MINOS, MINOS+ and NOA experiments.
Figure 25: The sensitivity of MINOS+ to the existence of sterile neutrinos, when combined with data from the Bugey ref:Bugey () reactor neutrino experiment. This figure assumes two years of MINOS+ running with a neutrino-dominated beam.

With more statistics at high energies MINOS+ is in the unique position to probe and significantly extend the reach of its searches for sterile neutrino signatures in the regions of parameter space favoured by LSND and MiniBooNE. Figure 25 shows a combination between the Bugey ref:Bugey () reactor experiment combined with the sensitivity of data taken with MINOS+ assuming two years of MINOS+ running with a neutrino-dominated beam. A combination with an experiment sensitive to (such as Bugey) is required with the MINOS+ data (sensitive primarily to ) in order to set a limit in the LSND style parameter space. The combined 90% MINOS+-Bugey C.L. excludes a significant amount of the parameter space where sterile neutrinos in a 3+1 model could explain the anomalies seen in past experiments.

With the increased statistics and higher flux in neutrinos, MINOS+ can also probe non-standard interaction neutrino physics. Figure 26 shows the increase in sensitivity MINOS+ is to the previous MINOS result with three different possible exposure combinations.

Figure 26: Showing the C.L. contours the MINOS+ experiment can achieve with three different projection running periods. The black contour is the current MINOS result, and the other ones are the MINOS+ sensitivity using simulation.

10 Conclusion

The MINOS/MINOS+ experiment has been contributing to the neutrino oscillation community for over a decade. In this time the most precise measurement of has been made. Such a precise measurement is an example of how powerful and necessary a two-detector setup will be for future neutrino oscillation experiments to overcome the large systematics from flux and cross section uncertainties. With the ability to differentiate between neutrinos and anti-neutrinos MINOS has measured oscillation parameters for both and show that they lie in good agreement. Since the discovery of to be non zero MINOS has been the first experiment to set constraints on the CP violating phase as well as incorporating both disappearance and appearance using a full three-flavour framework. MINOS+ will continue taking data and will improve on the MINOS results as well probe at higher energies to investigate the tension in the sterile neutrino debate.


The work of the MINOS and MINOS+ collaborations is supported by the US DoE, the UK STFC, the US NSF, the State and University of Minnesota, the University of Athens in Greece, and Brazil’s FAPESP and CNPq. We are grateful to the Minnesota Department of Natural Resources, the crew of the Soudan Underground Laboratory, and the personnel of Fermilab, for their vital contributions.


  • (1) Y. Fukuda et al. (Super-Kamiokande), Phys. Rev. Lett. 81, 1562 (1998).
  • (2) Q.R. Ahmad et al. (SNO), Phys. Rev. Lett. 87, 071301 (2001); Phys. Rev. Lett. 89, 011301 (2002).
  • (3) K. Eguchi et al. (KamLAND), Phys. Rev. Lett. 90, 021802 (2002).
  • (4) R. Acquafredda et al. (OPERA), New. J. Phys. 8, 303 (2006).
  • (5) K. Abe et al. (T2K), Nucl. Instrum. Meth. A659, 106-135 (2011).
  • (6) J.N Abdurashitov et al. (SAGE), Phys. Rev. C60, (1999).
  • (7) W. Hampel et al. (GALLEX), Phys. Lett. B447 (1999).
  • (8) M. Apollonio et al. (CHOOZ), Phys. Lett. B420, 397-404 (1998).
  • (9) J. K. Ahn et al. (RENO), Phys. Rev. Lett. 108, 191802 (2012).
  • (10) Y. Abe et al. (Double Chooz), Phys. Rev. D86, 052008 (2012).
  • (11) F.P. An et al. (Daya Bay), Chin. Phys. C37, 011001 (2013).
  • (12) P875: A Long-baseline Neutrino Oscillation Experiment at Fermilab, February 1995
  • (13) M. H. Ahn et al. (K2K), Phys. Rev. D 74, 072003 (2006).
  • (14) D. S. Ayres et al. (NO A), FERMILAB-PROPOSAL-0929, (2004).
  • (15) G. Tzanakos et al. (MINOS+), FERMILAB-PROPOSAL-1016 (2011).
  • (16) B. Pontecorvo, JETP 34, 172 (1958); V.N. Gribov and B. Pontecorvo, Phys. Lett. B28, 493 (1969); Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28, 870 (1962).
  • (17) S-B Kim et al. (RENO), Phys. Rev. Lett. 108, 191802 (2012).
  • (18) Y. Abe et al.(Double Chooz), Phys. Rev. D86, 052008 (2012).
  • (19) L. Wolfenstein, Phys. Rev. D17, 2369 (1978).
  • (20) S.P. Mikheyev, and A.Yu. Smirnov, Nuovo Cimento C, 9(1):12-26 (1986).
  • (21) C. Giunti, C.W. Kim, and M. Monteno, Nucl. Phys. B 521, 3 (1998).
  • (22) A. Cervera et al., Nucl. Phys. B 579, 17 (2000).
  • (23) K. Anderson et al., FERMILAB-DESIGN-1998-01 (1998).
  • (24) D. Michael et al. (MINOS), Nucl. Instrum. Meth. A596, 190 (2008).
  • (25) J.A. Formaggio and G.P. Zeller, Rev. Mod. Phys. 84, 1307 (2012).
  • (26) P. Adamson et al., The MINOS calibration detector. Nucl.Instrum.Meth., A556:119Ð133, 2006.
  • (27) M.A. Kordosky, PhD thesis, University of Texas at Austin (2004), FERMILAB-THESIS-2004-34.
  • (28) P.L. Vahle, PhD thesis, University of Texas at Austin (2004), FERMILAB-THESIS-2004-35.
  • (29) C. Backhouse, DPhil thesis, University of Oxford (2011), FERMILAB-THESIS-2011-17.
  • (30) S. Dytman, H. Gallagher, and M. Kordosky, arXiv:0806.2119 (2008).
  • (31) T.M. Cover and P.E. Hart, IEE Trans. Inform. Theory 13, 21 (1967).
  • (32) R. Ospanov, PhD thesis, University of Texas at Austin (2008).
  • (33) J.S. Marshall, PhD thesis, University of Cambridge (2008), FERMILAB-THESIS-2008-20.
  • (34) J.P. Ochoa, PhD thesis, Caltech (2009), FERMILAB-THESIS-2009-44.
  • (35) R. Toner, PhD thesis, University of Cambridge (2011), FERMILAB-THESIS-2011-53.
  • (36) A. Schreckenberger , PhD thesis, University of Minnesota (2013), FERMILAB-THESIS-2013-04.
  • (37) P. Adamson et al., (MINOS), Phys. Rev. Lett. 110, 171801 (2013).
  • (38) A. Holin, PhD thesis, University College London (2010), FERMILAB-THESIS-2010-41.
  • (39) G. Tinti, DPhil thesis, University of Oxford (2010), FERMILAB-THESIS-2010-44.
  • (40) P. Adamson et al. (MINOS), Phys. Rev. D86, 052007 (2012).
  • (41) M. Crouch, Proc. 20th Int. Cosmic Ray Conf. 6, 165 (1987).
  • (42) D. DeMuth et al. (Soudan 2), J. Astropart. Phys. 20, 533 (2004).
  • (43) J.J. Evans, DPhil thesis, University of Oxford (2008), FERMILAB-THESIS-2009-14.
  • (44) S.J. Coleman, PhD thesis, College of William & Mary (2011), FERMILAB-THESIS-2011-14.
  • (45) P. Adamson et al. (MINOS), Phys. Rev. Lett. 110, 251801 (2013).
  • (46) P. Adamson et al. (MINOS), Phys. Rev. Lett. 108, 191801 (2012).
  • (47) A.M Dziewonski, and D.L. Anderson, Phys. Earth Planet. Interiors. 25, 297 (1981).
  • (48) J.A.B. Coelho, PhD thesis, Universidade Estadual de Campinas (2012).
  • (49) G.L. Fogli et al., Phys. Rev. D86, 013012 (2012).
  • (50) J.W.F. Valle, Phys. Lett. B199, 432 (1987).
  • (51) M. C. GonzalezGarcia et al., Phys. Rev. Lett. 82, 3202 (1999).
  • (52) A. Friedland, C. Lunardini, and M. Maltoni, Phys. Rev. D70, 111301 (2004).
  • (53) G. Mitsuka et al. (Super-Kamiokande), Phys. Rev. D84, 113008 (2011).
  • (54) A. Aguilar et al. (LSND), Phys. Rev. D64, 112007 (2001).
  • (55) A.A. Aguilar-Arevalo et al. (MiniBooNE), Phys. Rev. Lett. 103, 111801 (2009); Phys. Rev. Lett. 110, 161801 (2013).
  • (56) G. Mention et al., Phys. Rev. D83, 073006 (2011).
  • (57) K.N. Abazajian et al., arXiv:1204.5379.
  • (58) P. Adamson et, al. (MINOS), Phys. Rev. Lett. 107, 011802 (2011).
  • (59) Nucl. Phys. B434, 503 (1995). The Bugey
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description