# Measurement of Neutrino Oscillations with Neutrino Telescopes

Authors: Jürgen Brunner
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

IceCube and ANTARES are the world-largest neutrino telescopes. They are successfully taking data, producing a wealth of scientific results. Whereas their main goal is the detection of cosmic neutrinos with energies in the TeV-PeV range, both have demonstrated their capability to measure neutrino oscillations by studying atmospheric neutrinos with energies of 10–50 GeV. After recalling the methods of these measurements and the first published results of these searches, the potential of existing, and planned low-energy extensions of IceCube and KM3Net are discussed. These new detectors will be able to improve the knowledge of the
atmospheric neutrino oscillation parameters, and in particular they might help to understand the neutrino mass hierarchy. Such studies,
which use atmospheric neutrinos, could be complemented by measurements in a long-baseline neutrino beam, which is discussed as a long-term future option.


---

## Body

## 1. Introduction

The main goal of neutrino telescopes such as IceCube [1] and ANTARES [2] is the observation of high-energy neutrinos from nonterrestrial sources. These telescopes are optimised for the detection of Cherenkov light induced by charged particles issued from neutrino interactions at TeV energies. To cope with the expected feeble neutrino fluxes at these energies, target masses of up to one Gton are equipped with low-density arrays of photomultipliers. Such a low-density of detector elements implies a rather high energy threshold for neutrino detection of about 50 GeV.Nevertheless both IceCube (exploiting its low-energy infill DeepCore [3]) and ANTARES developed dedicated analysis methods to access neutrino energies as low as 20 GeV, where the effect of neutrino oscillations starts to be measurable and both collaborations published recently first results on the measurement of parameters for atmospheric neutrino oscillations.Two principal event signatures can be distinguished. Track-like events are characterized by the presence of a long muon track, which can be identified by the light pattern from Cherenkov light emitted along its path. Such events are the result of muon neutrino charged current (CC) interactions or tau neutrino CC interactions followed by a muonic tau decay in conjunction with an analysis-dependent condition on the minimal muon energy. Cascade-like events are distinguished by the absence of such a long muon track. The low density of detector elements does not allow to separate hadronic from electromagnetic cascades. Therefore all remaining reaction channels such as neutral current (NC) reactions, electron and tau neutrino CC interactions (the latter followed by nonmuonic tau decays) are commonly called cascade-like events or just “cascades”.First, published results of ANTARES and IceCube are presented, followed by the discussion of very recent and still preliminary studies and their long-term potential. The second part of the paper is dedicated to the planned low-energy extensions of neutrino telescopes in the Mediterranean Sea as well as in the Antarctic ice. Their potential of studying atmospheric neutrinos is discussed. The final section evokes the possibilities of sending neutrino beams towards these future neutrino telescopes.


## 2. Atmospheric Neutrino Oscillations atEν>15 GeV

At energies ofEν>15 GeV, Earth matter effects can be ignored because the resonance energies for typical Earth densities are significantly lower (see Section 7 and Figure 6). Oscillation effects are prominently seen in the suppression of muon neutrinos. The survival probability of atmospheric νμ in the framework of three-flavour mixing is given as
(1)P(νμ⟶νμ)=1-4|Uμ1|2|Uμ2|2sin2(1.27Δm212LEν)-4|Uμ1|2|Uμ3|2sin2(1.27Δm312LEν)-4|Uμ2|2|Uμ3|2sin2(1.27Δm322LEν),
where L is the travel path (in km) of the neutrino through the earth and Eν, its energy (in GeV). Uαi is the 3×3 PMNS matrix which describes the mixing between flavour eigenstates νe,νμ,ντ and mass eigenstates ν1, ν2, and ν3 and Δmij2=mi2-mj2 (in eV2) is the difference of the squares of the masses of the corresponding neutrino mass eigenstates. The numbering of the mass eigenstates ν1, ν2, and ν3 is chosen such that
(2)|Ue1|2>|Ue2|2>|Ue3|2.
For atmospheric neutrinos, L≤D, where D is the Earth diameter, is always satisfied. If the condition Eν>15 GeV is taken into account, the term sin2(1.27Δm212L/Eν) does not exceed 5·10-3 when using Δm212 from [4]. This term can safely be ignored as well as differences between Δm312 and Δm322 and (1) simplifies to
(3)P(νμ⟶νμ)=1-4(1-|Uμ3|2)×|Uμ3|2sin2(1.27Δm322LEν).
Results could in principle be extracted in terms of |Uμ3|2 and |Δm322| which are the two oscillation parameters in (3). To maintain compatibility with earlier results, a mixing angle sin2θ23=|Uμ3|2 is defined, ignoring the 2.4% deviation from 1 of cos2θ13=0.976 [5, 6]. This leads to the usual two-flavour description
(4)P(νμ⟶νμ)=1-sin22θ23sin2(16200Δm322cosΘEν),
where L has been replaced by D·cosΘ with Θ the zenith angle. The transition probability P depends now on only two oscillation parameters, |Δm322| and sin22θ23, which determine the behaviour for the atmospheric neutrino oscillations. With |Δm322|=2.32·10-3 eV2 and sin22θ23=1 from [4] one expects the first oscillation maximum; that is, P(νμ→νμ)=0 for vertical upgoing neutrinos (cosΘ=1) at Eν=24 GeV. Muons induced by 24 GeV neutrinos can travel up to 120 m in sea water or Antarctic ice. Both ANTARES and IceCube have performed studies of neutrino oscillations using the suppression of muon neutrino events in the vicinity of this first oscillation maximum. Results of these studies are presented in the following.


## 3. The ANTARES Detector

A detailed description of the ANTARES detector can be found in [2]. The detector consists of 12 lines and equipped with photosensors and a junction box which distributes the power and clock synchronization signals to the lines and collects the data. The junction box is connected to the shore by a 42 km electro-optical cable. The length of the detection lines is 450 m, of which the lowest 100 m are not instrumented. Their horizontal separation is about 65 m and they are arranged to form regular octagon on the sea floor. They are connected to the junction box with the help of a submarine using wet-mateable connectors. Each line comprises 25 storeys each separated by a vertical distance of 14.5 m. The lines are kept taut by a buoy at the top of the line and an anchor on the seabed. The movement of the line elements due to the sea currents is continuously measured by an acoustic calibration system with an accuracy of 10 cm [7]. Each storey contains three 45° downward-looking 10′′ photomultiplier tubes (PMT) inside pressure resistant glass spheres—the optical modules [8]. Some of the storeys contain supplementary calibration equipment such as acoustic hydrophones or optical beacons [9].The signals of each photomultiplier are read out by two ASICs. The charges and arrival times of the PMT signals are digitised and stored for transfer to the shore station [10]. The time stamps are synchronised by a clock signal which is sent at regular intervals from the shore to all electronic cards. The overall time calibration is better than 0.5 ns [11]. Therefore the time resolution of the signal pulses is limited by the transit time spread of the photomultipliers (σ~1.3 ns) [12] and by chromatic dispersion for distant light sources. All data are sent to the shore station. With the observed optical background rate of 70 kHz per PMT at the single photon level, this produces a data flow of several Gbit/s to the shore. In the shore station a PC farm performs data filtering to reduce the data rate by at least a factor of 100 [13]. Several trigger algorithms are applied depending on the requested physics channel and on the observed optical noise.


## 4. ANTARES Neutrino Oscillation Measurement

This analysis is based on data taken with the ANTARES detector between March 2007 and December 2010. Until December 2007, ANTARES operated in a 5-line configuration, followed by several months of operation with 10 installed detector lines. The detector construction was completed in May 2008. All physics runs taken under normal conditions have been used. The events selected by two tight trigger conditions are used [2]. The analysed sample consists of 293 million triggers, dominated by atmospheric muons, corresponding to a detector live time of 863 days.Downgoing atmospheric muons were simulated with the program MUPAGE [14, 15] which provides parameterised muon bundles at the detector. Alternatively, also CORSIKA [16] was used for cross checks and systematic studies. Upgoing neutrinos were simulated according to the parameterisations from [17–19] in the energy range from 10 GeV to 10 PeV. The Cherenkov light, produced inside or in the vicinity of the detector instrumented volume, was propagated taking into account light absorption and scattering in sea water [20]. The angular acceptance, quantum efficiency, and other characteristics of the PMTs were taken from [8] and the overall geometry corresponded to the layout of the ANTARES detector [2]. The optical noise was simulated from counting rates observed in the data. At the same time, the definition of active and inactive channels has been applied from data runs as well. The generated statistics correspond to an equivalent observation time of 100 years for atmospheric neutrinos and ten months for atmospheric muons.The algorithm used to determine the zenith angleΘR of the muon track is described in [21]. A first step of this procedure is a strict hit selection, which aims at the selection of those Cherenkov photon hits from the muon track, which did not undergo any substantial scattering on their path. If the selected hits occur only on one detector line, a single-line fit is performed. No azimuth angle is determined in this case due to the rotational symmetry of the problem. This does not affect the present measurement, as the oscillation probability does not depend on the azimuth angle (see (4)). If selected hits occur instead on several detector lines, a multiline fit is performed which provides both the zenith and azimuth angles of the track. The inclusion of single-line events is a special feature of the reconstruction method and allows to significantly lower the energy threshold of the final atmospheric neutrino sample. Whereas for multiline events, the threshold energy of the final neutrino sample is about 50 GeV due to the 65 m horizontal spacing between lines, single-line events are reconstructed down to 20 GeV for nearly vertical tracks, accessing events with Eν/L values close to the first oscillation maximum.The neutrino energy is estimated from the observed muon range in the detector (ER). The presence of a hadronic shower at the neutrino vertex is ignored, as well as the fact that the muon might leave or enter the detector, making only a fraction of its actual range available for measurement.Downgoing atmospheric muons might contaminate the event sample of upgoing atmospheric neutrinos if misreconstructed. A cut on the quality of the reconstructed tracks is needed to reduce this contamination and derive reliably neutrino oscillation parameters from the data set. The cut is chosen in such a way that the contamination of misreconstructed atmospheric muons remains below 5%. The selected event sample consists of 2126 events.The oscillation parameters are extracted from aχ2 minimisation of an event distribution in ER/cosΘR. Systematic uncertainties lead to correlated distortion of this distribution. They are implemented as pull factors ϵ and η which act as independent free normalisations for the single-line, and multiline sample, respectively.Figure1 shows the event distribution in ER/cosΘR and the fraction of measured and simulated events with respect to the nonoscillation Monte Carlo hypothesis. A clear event deficit for ER/cosΘR<60 GeV can be seen as expected assuming atmospheric neutrino oscillations.(a) Distribution ofER/cosΘR for selected events. Black crosses are data with statistical uncertainties, whereas the blue histogram shows simulations of atmospheric neutrinos without neutrino oscillations plus the residual background from atmospheric muons. The red histogram shows the result of the fit. (b) The fraction of events with respect to the nonoscillation hypothesis as function of ER/cosΘR. Same color code as for the left figure (Figure from [22]).
(a)
(b)Restricting theχ2 minimisation to parameters in the physically allowed region yields the red curve of Figure 1 with |Δm322|=3.1·10-3eV2 and sin22θ23=1.00. The corresponding pull factors are ϵ=-0.138 and η=-0.142. The fit yields χ2/NDF=17.1/21. For the nonoscillation hypothesis, that is, sin22θ23=0, χ2/NDF=31.1/23 is obtained. The pull factors in this case are ϵ=-0.302 and η=-0.196. The event deficit in the single-line channel is seen here as ϵ becoming lower than η. Requiring in addition ϵ=η the χ2 increases further to χ2/NDF=40.0/24, which has a probability of only 2.1%. This measurement is converted into a 90% C.L. contour of the oscillation parameters and is shown in Figure 5. More details on this analysis can be found in [22].


## 5. The IceCube Detector

IceCube is a cubic-kilometer neutrino detector installed in the ice at the geographic South Pole [1]. It consists of 5,160 digital optical modules (DOMs). Each of them contains a photomultiplier tube, supporting hardware and electronics with in situ pulse digitization [23, 24] inside a pressure glass sphere. These optical modules are arranged on 86 strings frozen into the antarctic ice at depths from 1450 m to 2450 m each holding 60 DOMs. The primary (high-energy) detector has a spacing of 17 m between sensors and an average horizontal distance of 125 m between neighboring strings. The low-energy infill array DeepCore consists of eight dedicated strings with a typical spacing of 70 m deployed near the center of the IceCube array. On the dedicated Deep- Core strings, the sensors are concentrated in the clearest deep ice, with a denser 7 m vertical spacing.Optical properties of South Pole ice are depth dependent. Typical optical absorption lengths are 100–140 m and typical optical effective scattering lengths are 25–35 m. The latter limits the time resolution of the signal pulses from Cherenkov light. Unlike sea water, Antarctic ice has a low concentration in radioactive isotopes. In particular40K, abundant in salty sea water, is absent in ice. Furthermore, no light-producing biological activity occurs. Therefore, the observed optical background rate at the single photon level is only 500 Hz per PMT in IceCube.Beginning with the installation of the first string in 2005, IceCube has been operated in approximately year-long data taking seasons. During the Antarctic summer seasons in December and January new strings had been deployed at a pace up to 20 strings per year. IceCube construction was finished in December 2010. From May 2010 to April 2011 79 strings had been active (IC79), after this period data taking proceeds with the complete detector (IC86). The results which will be presented in the following have been obtained from analysing these two data-taking periods.


## 6. IceCube Neutrino Oscillation Measurements

Downgoing atmospheric muons were simulated with CORSIKA [16]. The flux of atmospheric neutrinos has been taken from [19]. Both the flux normalisation and the spectral index are allowed to vary and they are included as nuisance parameters in the analyses described below.


### 6.1. IC79 Zenith Angle Analysis

A first analysis [25] exploited two samples of upgoing muon neutrino events which had been collected during a 318.9-day period with the IC79 configuration, excluding periods of calibration runs, partial detector configurations, and detector downtime. The first sample was obtained from high-energy events using data from the entire IceCube detector. The second sample, selected from events starting in the DeepCore volume, was using the surrounding IceCube array as an active veto to reject atmospheric muon background and high-energy (>100 GeV) neutrinos [26].The directions of the neutrino-induced muon tracks in the high-energy sample were determined with the standard maximum likelihood muon track reconstruction of IceCube [27]. For low-energy events, the same method was applied as an initial step. However, the standard hypothesis of a throughgoing track is not appropriate at low energies. In a subsequent step, the finite length of the low-energy muons which might start or stop inside the equipped detector volume was taken into account. Quality cuts on variables related to the track fit eliminated misreconstructed downgoing atmospheric muons. The final event sample consisted of 719 low-energy events, while the high-energy sample contained 39,638 events after final cuts. The simulated neutrino energy distribution of the two samples is shown in Figure 2. The high-energy sample is not affected by standard oscillations and it serves as a reference sample which absorbs systematic effects related to the overall normalisation. Instead neutrino oscillations are expected to deplete the low-energy sample.Figure 2
Expected distribution of the neutrino energy of atmospheric neutrinos in the low-energy (DeepCore) and in the high-energy (IceCube) samples according to simulations (Figure from [25]).The oscillation parameters are extracted from the zenith angle distribution, shown in Figure3 of the events in both samples. Within the two samples, no further use of energy is made. A χ2-minimization including various systematic effects is applied to the 20 bins in cosΘ in both samples.Data and Monte Carlo expectation at world average oscillation parameters [4] and at the nonoscillation scenario for the low-energy sample (a) and for the high-energy sample (b). The shaded bands illustrate separately those systematic effects which affect the normalisation (norm) and the shape (shape) of the distributions (figure from [25]).
(a)
(b)The fit yields|Δm322|=2.3-0.6+0.5·10-3eV2 and sin22θ23>0.93 at 68% C.L. The nonoscillation hypothesis can be rejected at 5.6 σ. All pull factors related to systematic effects are found within one σ from their nominal values. The resulting 90% C.L. contour is shown in Figure 5.Both analyses, which have been discussed so far, are final and their results have been published. Two more studies of IceCube data are under way. They have reported preliminary results recently.


### 6.2. IC79 Two-Dimensional Analysis

As this analysis [28] is well based on data in the IC79 configuration, here 312 days are used. The event selection requires that neutrino event candidates should have their interaction vertices inside the instrumented volume of the Deep-Core detector, which reduces background from high-energy neutrino events and downgoing atmospheric muons. When using the surrounding IceCube strings as a veto, their contributions can be partly discarded. Apart from the veto cut, only weak selection criteria are applied, which leads to a final event sample of 8117 events with an estimated purity in νμ events of 70%. The oscillation parameters are fitted on a two-dimensional distribution in zenith angle and energy. The zenith angle is reconstructed using standard IceCube methods [27]. The energy is derived from the measured track length [26]. For the likelihood, the standard Poisson formulation is used. For each bin (i,j), the probability is calculated to observe dij events in the measured data, given sij events in the simulated data. In addition to fitting the oscillation parameters, nuisance parameters qk are also left free in the fit which absorb certain systematic uncertainties. Gaussian priors are used for each nuisance parameter k with a mean value 〈qk〉 and a width σk which are added to the likelihood. The full likelihood expression has the form
(5)-LLH=∑i,j‍(sij-dijlnsij)+12∑k‍(qk-〈qk〉σk)2.
Other systematics, which are not directly implemented in the fit, are evaluated by separate simulations. As a preliminary result of this analysis, the nonoscillation hypothesis is rejected by a likelihood ratio corresponding to 5.1 standard deviations.  The corresponding 90% C.L. contour is also shown in Figure 5.


### 6.3. IC86 Analysis

The second new analysis [29] uses a complementary approach. A strict selection of events with clusters of essentially unscattered Cherenkov photon hits is applied, following the procedure which had previously been used for the ANTARES oscillation search [21]. The described hit selection and fitting procedure has been adapted to IceCube/DeepCore. A lower energy threshold and a better zenith angle distribution are observed compared to the performances of the above mentioned analyses. Also the dependence of the result on uncertainties in the description of the ice properties is reduced as scattered hits are largely discarded in this analysis. The neutrino energy for the νμ CC sample is estimated from the sum of the fitted cascade energy at the neutrino vertex and the muon energy obtained from the measured muon range. Figure 4 shows the performance of the reconstructions used for the final event selection. A zenith angle resolution of 7 degrees, comparable with the kinematic angle between the neutrino, and the muon in the sample, is obtained. The energy is resolved with an error of 0.25 in log10(E). The analysis was applied to the first year of data taken by the complete IceCube detector (IC86), from May 2011 until April 2012. In 343 days of lifetime, 1487 neutrino events were found. For the likelihood, the standard Poisson formulation with added nuisance terms is also used here (see (5)). The best fit values for single parameters obtained from the likelihood profile are sin2(2θ23)=1(>0.93 at 68% C.L.) and |Δm322|=(2.4±0.4)·10-3 eV2. The nuisance parameters at the best fit point are all well within their one-sigma ranges. The resulting zenith angle and energy distributions for the best fit values are shown in Figure 4. The corresponding 90% C.L. contour is also overlaid in Figure 5 and represents currently the best such measurement provided by a neutrino telescope.(a) Performance of reconstructions for theνμ CC component of the IC86 analysis in zenith angle (top) and energy (bottom). (b) Zenith angle and energy distributions of data and best fit simulation. Statistical errors are added to the data; systematic errors are attributed to the simulation (figure from [29]).
(a)
(b)Figure 5
90% C.L. contours from recent measurements of the neutrino oscillation parameters by ANTARES [22] and IceCube [25, 28, 29] compared to earlier measurements by MINOS [30] and Super-Kamiokande [31] (Figure kindly provided by S. Euler, Aachen University).Figure 6
Neutrino oscillation probabilitiesP(νμ→νμ) as a function of neutrino energy for different zenith angles cosΘ=0.2, 0.4, 0.6, and 0.8 which correspond to baselines of 2548; 5097; 7645; 10194 km. Oscillation parameters from a global fit are used. The group of red lines is for NH, blue for IH, and ϕCP varied in steps of 30° between 0°  and 330°.Whereas the ANTARES oscillation measurement was based on four years of data, the current IceCube analyses did all use just one year of data. All of them are currently statistically limited. It can be expected that the measurement accuracy improves with roughlyT (T being the detector lifetime) for a couple of years. But to reach a sensitivity, which would be competitive to the current MINOS [30] and SuperKamiokande [31] results, increased statistics alone are not sufficient.  An improved energy estimator combined with a clean high-statistics low-energy event sample is needed. It has to be seen whether this can be accomplished within the DeepCore/IceCube setup or whether new, denser detectors are required.


## 6.1. IC79 Zenith Angle Analysis

A first analysis [25] exploited two samples of upgoing muon neutrino events which had been collected during a 318.9-day period with the IC79 configuration, excluding periods of calibration runs, partial detector configurations, and detector downtime. The first sample was obtained from high-energy events using data from the entire IceCube detector. The second sample, selected from events starting in the DeepCore volume, was using the surrounding IceCube array as an active veto to reject atmospheric muon background and high-energy (>100 GeV) neutrinos [26].The directions of the neutrino-induced muon tracks in the high-energy sample were determined with the standard maximum likelihood muon track reconstruction of IceCube [27]. For low-energy events, the same method was applied as an initial step. However, the standard hypothesis of a throughgoing track is not appropriate at low energies. In a subsequent step, the finite length of the low-energy muons which might start or stop inside the equipped detector volume was taken into account. Quality cuts on variables related to the track fit eliminated misreconstructed downgoing atmospheric muons. The final event sample consisted of 719 low-energy events, while the high-energy sample contained 39,638 events after final cuts. The simulated neutrino energy distribution of the two samples is shown in Figure 2. The high-energy sample is not affected by standard oscillations and it serves as a reference sample which absorbs systematic effects related to the overall normalisation. Instead neutrino oscillations are expected to deplete the low-energy sample.Figure 2
Expected distribution of the neutrino energy of atmospheric neutrinos in the low-energy (DeepCore) and in the high-energy (IceCube) samples according to simulations (Figure from [25]).The oscillation parameters are extracted from the zenith angle distribution, shown in Figure3 of the events in both samples. Within the two samples, no further use of energy is made. A χ2-minimization including various systematic effects is applied to the 20 bins in cosΘ in both samples.Data and Monte Carlo expectation at world average oscillation parameters [4] and at the nonoscillation scenario for the low-energy sample (a) and for the high-energy sample (b). The shaded bands illustrate separately those systematic effects which affect the normalisation (norm) and the shape (shape) of the distributions (figure from [25]).
(a)
(b)The fit yields|Δm322|=2.3-0.6+0.5·10-3eV2 and sin22θ23>0.93 at 68% C.L. The nonoscillation hypothesis can be rejected at 5.6 σ. All pull factors related to systematic effects are found within one σ from their nominal values. The resulting 90% C.L. contour is shown in Figure 5.Both analyses, which have been discussed so far, are final and their results have been published. Two more studies of IceCube data are under way. They have reported preliminary results recently.


## 6.2. IC79 Two-Dimensional Analysis

As this analysis [28] is well based on data in the IC79 configuration, here 312 days are used. The event selection requires that neutrino event candidates should have their interaction vertices inside the instrumented volume of the Deep-Core detector, which reduces background from high-energy neutrino events and downgoing atmospheric muons. When using the surrounding IceCube strings as a veto, their contributions can be partly discarded. Apart from the veto cut, only weak selection criteria are applied, which leads to a final event sample of 8117 events with an estimated purity in νμ events of 70%. The oscillation parameters are fitted on a two-dimensional distribution in zenith angle and energy. The zenith angle is reconstructed using standard IceCube methods [27]. The energy is derived from the measured track length [26]. For the likelihood, the standard Poisson formulation is used. For each bin (i,j), the probability is calculated to observe dij events in the measured data, given sij events in the simulated data. In addition to fitting the oscillation parameters, nuisance parameters qk are also left free in the fit which absorb certain systematic uncertainties. Gaussian priors are used for each nuisance parameter k with a mean value 〈qk〉 and a width σk which are added to the likelihood. The full likelihood expression has the form
(5)-LLH=∑i,j‍(sij-dijlnsij)+12∑k‍(qk-〈qk〉σk)2.
Other systematics, which are not directly implemented in the fit, are evaluated by separate simulations. As a preliminary result of this analysis, the nonoscillation hypothesis is rejected by a likelihood ratio corresponding to 5.1 standard deviations.  The corresponding 90% C.L. contour is also shown in Figure 5.


## 6.3. IC86 Analysis

The second new analysis [29] uses a complementary approach. A strict selection of events with clusters of essentially unscattered Cherenkov photon hits is applied, following the procedure which had previously been used for the ANTARES oscillation search [21]. The described hit selection and fitting procedure has been adapted to IceCube/DeepCore. A lower energy threshold and a better zenith angle distribution are observed compared to the performances of the above mentioned analyses. Also the dependence of the result on uncertainties in the description of the ice properties is reduced as scattered hits are largely discarded in this analysis. The neutrino energy for the νμ CC sample is estimated from the sum of the fitted cascade energy at the neutrino vertex and the muon energy obtained from the measured muon range. Figure 4 shows the performance of the reconstructions used for the final event selection. A zenith angle resolution of 7 degrees, comparable with the kinematic angle between the neutrino, and the muon in the sample, is obtained. The energy is resolved with an error of 0.25 in log10(E). The analysis was applied to the first year of data taken by the complete IceCube detector (IC86), from May 2011 until April 2012. In 343 days of lifetime, 1487 neutrino events were found. For the likelihood, the standard Poisson formulation with added nuisance terms is also used here (see (5)). The best fit values for single parameters obtained from the likelihood profile are sin2(2θ23)=1(>0.93 at 68% C.L.) and |Δm322|=(2.4±0.4)·10-3 eV2. The nuisance parameters at the best fit point are all well within their one-sigma ranges. The resulting zenith angle and energy distributions for the best fit values are shown in Figure 4. The corresponding 90% C.L. contour is also overlaid in Figure 5 and represents currently the best such measurement provided by a neutrino telescope.(a) Performance of reconstructions for theνμ CC component of the IC86 analysis in zenith angle (top) and energy (bottom). (b) Zenith angle and energy distributions of data and best fit simulation. Statistical errors are added to the data; systematic errors are attributed to the simulation (figure from [29]).
(a)
(b)Figure 5
90% C.L. contours from recent measurements of the neutrino oscillation parameters by ANTARES [22] and IceCube [25, 28, 29] compared to earlier measurements by MINOS [30] and Super-Kamiokande [31] (Figure kindly provided by S. Euler, Aachen University).Figure 6
Neutrino oscillation probabilitiesP(νμ→νμ) as a function of neutrino energy for different zenith angles cosΘ=0.2, 0.4, 0.6, and 0.8 which correspond to baselines of 2548; 5097; 7645; 10194 km. Oscillation parameters from a global fit are used. The group of red lines is for NH, blue for IH, and ϕCP varied in steps of 30° between 0°  and 330°.Whereas the ANTARES oscillation measurement was based on four years of data, the current IceCube analyses did all use just one year of data. All of them are currently statistically limited. It can be expected that the measurement accuracy improves with roughlyT (T being the detector lifetime) for a couple of years. But to reach a sensitivity, which would be competitive to the current MINOS [30] and SuperKamiokande [31] results, increased statistics alone are not sufficient.  An improved energy estimator combined with a clean high-statistics low-energy event sample is needed. It has to be seen whether this can be accomplished within the DeepCore/IceCube setup or whether new, denser detectors are required.


## 7. Projects for Low-Energy Extensions of Neutrino Telescopes

All three mixing angles and both mass square differences of the neutrino mass eigenstates are known after the measurement of the mixing angleθ13 [5, 6]. Global fits of all experimental input [32] provide a coherent picture of the oscillation parameters. In the following, the best fit values from this analysis are used if not stated otherwise. One yet unknown feature of the oscillation scheme is the ordering of the mass eigenstates of the neutrinos, the neutrino mass hierarchy (MH). Whereas m2>m1 is determined through matter effects in the sun, which influence the solar neutrino flux measurements, the options m3>m2 (normal hierarchy: NH) and m1>m3 (inverted hierarchy: IH) are both allowed. They could be distinguished by measuring matter effects in the Earth, which is considered to be in reach for the next generation of experiments. In the approximation of a constant electron density ne and for single Δm322 dominance as derived for (3), its main effect is the modification of the mixing strength related to θ13, which has to be replaced by the effective matter mixing angle θ13m [33]:
(6)sin2θ13m=sin(2θ13)(22GFneEν/Δm322-cos2θ13)2+(sin2θ13)2
with GF the Fermi constant. The sign of ne is positive for neutrino interactions and negative for antineutrinos. A resonance, that is, sin2θ13m=1, occurs for neutrinos and Δm322>0 (NH) and for antineutrinos and Δm322<0 (IH). For typical densities in the Earth mantle, the resonance is observed at Eν=6–8 GeV.Neutrino oscillation probabilities are now calculated in a full three-flavour scheme including matter effects in the Earth. The Globes package [34, 35] is used for this purpose. Figure 6 shows the survival probabilities for upgoing atmospheric νμ for different zenith angles as function of neutrino energy. The Earth density profile is taken from the PREM model [36]. The two colors refer to the two MH choices. A striking difference of the oscillation probabilities for the two options can be seen for zenith angles cosΘ>0.4 and energies close to the resonance energy of the matter effect in the Earth mantle. Neutrino interactions in this energy range can occur in various channels: quasielastic, via nuclear resonances or as deep-inelastic reactions. Moreover nuclear effects should be taken into account. Most of the studies mentioned below are based on the event generator GENIE [37] which comprises all mentioned effects.It has been proposed to determine the MH by measuring atmospheric neutrinos [38] which has motivated proposals for low-energy extensions of existing and future neutrino telescope projects. To measure reliably neutrinos with energies of few GeVs, higher densities of the photosensor arrays are needed. Two such detectors have been proposed, which are presented in the following.


### 7.1. ORCA

KM3NeT [39] will be the next generation neutrino telescope in the Mediterranean Sea with an effective volume of several cubic kilometers. Funds for a first phase of the project, corresponding to about 20% of the total envisaged budget, are meanwhile available. A dedicated “Oscillation Research with Cosmics in the Abyss” (ORCA) is currently under evaluation as a low energy part of the KM3Net detector [40]. ORCA might consist of at least 50 detector lines, each equipped with 20 optical modules. Each optical module hosts 31 3-inch phototubes, as designed in the KM3NeT TDR [39]. Optical modules with many small phototubes have an enhanced photon counting capability and provide better directional information than larger tubes. This can be exploited in reconstruction algorithms and for background rejection. The readout will be largely simplified compared to ANTARES. For each detected pulse, only a timestamp and the “time over threshold” will be recorded. Trigger algorithms will be applied in an onshore computing farm, as successfully done for ANTARES. The distance of the detector lines on the sea bed will be approximately 20 m, a value determined by deployment safety. The vertical distance between modules on a detector line will be chosen between 5 m and 10 m, the precise value being subject to optimization. This leads to an equipped water mass of about 2 Mtons. A newly developed and simplified deployment scheme [39] will allow deploying several detector lines during one sea operation. It can be envisaged to complete the detector installation within 5 years.Studies have been carried out, to evaluate the achievable resolution in zenith angle and energy. The reconstruction methods have been derived from ANTARES and have been adapted to the denser geometry of the ORCA detector. Preliminary results have been published recently [40] and they are shown in Figure 7. For events which interact inside the instrumented volume, a zenith angle resolution close to the kinematical limit between neutrino and muon direction can be reached. 68% of events with a muon of a true energy between 3 and 8 GeV can be reconstructed with a precision of 3 GeV from range estimation of the fitted track. An effective mass of 1.8 Mtons is observed for νμ CC events with Eν>5 GeV, having their interaction vertex inside the instrumented volume.(a) Median of the difference between the true neutrino zenith angle and the reconstructed zenith angle calculated for all the upgoing reconstructed events with the interaction vertex inside the instrumented volume. For reference, the median of the difference between the true neutrino zenith angle and the true muon zenith angle is shown. (b) Reconstructed muon energy as a function of the muon true energy. The drawn black line gives the median value and the red lines indicate the 16% and 84% quantiles—the equivalent of a one-sigma range for a Gaussian distribution (figure from [40]).
(a)
(b)


### 7.2. PINGU

PINGU stands for “Precision IceCube Next Generation Upgrade”, and it is planned as a dense core extension of the existing IceCube/DeepCore ensemble [41]. Several geometries are under consideration. The two most promising configurations include additional 20–40 detector strings, each holding 60 to 120 light detectors, deployed within the DeepCore volume (see Figure 8 for the example of the 20-lines configuration). For the 20-lines configuration the dependence of the effective mass for a sample of νμ CC events and a condition about a minimal number of optical modules hit is shown in Figure 8. Different reconstruction algorithms are currently tested for PINGU, some of them derived from IceCube tools (see above) and some others are under development specifically for PINGU. Currently a median in zenith angle resolution of 8–15 degrees is reported for neutrino energies between 2 and 20 GeV whereas the median of the absolute neutrino energy resolution can be conveniently parametrised by 0.7GeV+0.2·Eν for an energy reconstruction which exploits the total brightness of fully contained events.(a) Overhead view of the PINGU geometry layout with 20 additional detector strings. (b) Effective volume for PINGU in the baseline configuration (20 strings) in comparison with DeepCore. A fiducial volume defined by a cylinder with radius of 75 m and a height of 332 m is assumed. A threshold of 20 hits per event was assumed (figure from [41]).
(a)
(b)


### 7.3. Sensitivity for Neutrino Mass Hierarchy Determination with Atmospheric Neutrinos

Despite the seemingly easy distinction of the two MH hypotheses for a given oscillation channel (as illustrated in Figure6), the distinction of the two hierarchy hypotheses with atmospheric neutrinos is challenging due to the cancellation of contributions from neutrinos and antineutrinos. Further cancellation effects occur between νe and νμ assuming an unavoidably imperfect flavour tagging method. The question of flavour tagging, that is, distinguishing cascade from track signatures and thereby isolating the νμ CC, signal has still to be improved—in ORCA as well as in PINGU. The reachable sensitivity for MH is further attenuated by the finite energy and angular resolutions of the considered detectors which are discussed above. Last but not least MH effects can be partly masked by readjusting oscillation parameters (in particular Δm322 and sin22θ23) within their present uncertainties.Both ORCA and PINGU have yet to provide a sensitivity calculation for MH, including all these effects plus detector and flux related sources of systematic uncertainties. The most detailed sensitivity study to date is based on input from PINGU performance figures and can be found in [42]. Its main result is shown in Figure 9. A 1.3 σ to 2.9 σ significance might be reached for three years running with the 20-lines PINGU detector depending on the true MH choice and the true values for other oscillation parameters. The result is based on a rather optimistic assumption of a cascade misidentification fraction of 5% and it assumes that the neutrino energy can be determined with a perfectly Gaussian resolution of 25%.Figure 9
Number ofσ for the hierarchy discovery of a 20-lines PINGU detector as a function of time for two extreme cases of the true parameter values (figure from [42]).The cancellation of matter effects between neutrinos and anti-neutrinos can be partly avoided by determining for eachνμ CC event the hadronic and muon part separately. This would allow to measure the inelasticity variable Bjorken-y which in turn could be used to statistically separate neutrinos from anti-neutrinos. The best separation would be reached for events of high inelasticity (large y), but these are hardly distinguishable from cascades due to the low-energy muon tracks. The potential of the y determination for the MH sensitivity is discussed in [43]. When assuming a realistic experimental smearing of the measured quantities, an enhancement of 25% is quoted in [43] by measuring hadrons and muons separately.The study of the neutrino mass hierarchy is not the only physics goal of PINGU and ORCA. The knowledge of the oscillation parametersΔm322 and sin2θ23 itself can be improved. The octant of the mixing angle can be determined, if its value is sufficiently distinct from maximal mixing and |Δm322| can be measured with a precision better than 1% based on one year of data with PINGU alone [41] under the assumption that systematic effects such as the uncertainty on the ice model do not deteriorate this measurement too strongly. Further, indirect dark matter searches from annihilation processes of gravitationally confined particles can be extended toward lower masses compared to current results from ANTARES and IceCube.


## 7.1. ORCA

KM3NeT [39] will be the next generation neutrino telescope in the Mediterranean Sea with an effective volume of several cubic kilometers. Funds for a first phase of the project, corresponding to about 20% of the total envisaged budget, are meanwhile available. A dedicated “Oscillation Research with Cosmics in the Abyss” (ORCA) is currently under evaluation as a low energy part of the KM3Net detector [40]. ORCA might consist of at least 50 detector lines, each equipped with 20 optical modules. Each optical module hosts 31 3-inch phototubes, as designed in the KM3NeT TDR [39]. Optical modules with many small phototubes have an enhanced photon counting capability and provide better directional information than larger tubes. This can be exploited in reconstruction algorithms and for background rejection. The readout will be largely simplified compared to ANTARES. For each detected pulse, only a timestamp and the “time over threshold” will be recorded. Trigger algorithms will be applied in an onshore computing farm, as successfully done for ANTARES. The distance of the detector lines on the sea bed will be approximately 20 m, a value determined by deployment safety. The vertical distance between modules on a detector line will be chosen between 5 m and 10 m, the precise value being subject to optimization. This leads to an equipped water mass of about 2 Mtons. A newly developed and simplified deployment scheme [39] will allow deploying several detector lines during one sea operation. It can be envisaged to complete the detector installation within 5 years.Studies have been carried out, to evaluate the achievable resolution in zenith angle and energy. The reconstruction methods have been derived from ANTARES and have been adapted to the denser geometry of the ORCA detector. Preliminary results have been published recently [40] and they are shown in Figure 7. For events which interact inside the instrumented volume, a zenith angle resolution close to the kinematical limit between neutrino and muon direction can be reached. 68% of events with a muon of a true energy between 3 and 8 GeV can be reconstructed with a precision of 3 GeV from range estimation of the fitted track. An effective mass of 1.8 Mtons is observed for νμ CC events with Eν>5 GeV, having their interaction vertex inside the instrumented volume.(a) Median of the difference between the true neutrino zenith angle and the reconstructed zenith angle calculated for all the upgoing reconstructed events with the interaction vertex inside the instrumented volume. For reference, the median of the difference between the true neutrino zenith angle and the true muon zenith angle is shown. (b) Reconstructed muon energy as a function of the muon true energy. The drawn black line gives the median value and the red lines indicate the 16% and 84% quantiles—the equivalent of a one-sigma range for a Gaussian distribution (figure from [40]).
(a)
(b)


## 7.2. PINGU

PINGU stands for “Precision IceCube Next Generation Upgrade”, and it is planned as a dense core extension of the existing IceCube/DeepCore ensemble [41]. Several geometries are under consideration. The two most promising configurations include additional 20–40 detector strings, each holding 60 to 120 light detectors, deployed within the DeepCore volume (see Figure 8 for the example of the 20-lines configuration). For the 20-lines configuration the dependence of the effective mass for a sample of νμ CC events and a condition about a minimal number of optical modules hit is shown in Figure 8. Different reconstruction algorithms are currently tested for PINGU, some of them derived from IceCube tools (see above) and some others are under development specifically for PINGU. Currently a median in zenith angle resolution of 8–15 degrees is reported for neutrino energies between 2 and 20 GeV whereas the median of the absolute neutrino energy resolution can be conveniently parametrised by 0.7GeV+0.2·Eν for an energy reconstruction which exploits the total brightness of fully contained events.(a) Overhead view of the PINGU geometry layout with 20 additional detector strings. (b) Effective volume for PINGU in the baseline configuration (20 strings) in comparison with DeepCore. A fiducial volume defined by a cylinder with radius of 75 m and a height of 332 m is assumed. A threshold of 20 hits per event was assumed (figure from [41]).
(a)
(b)


## 7.3. Sensitivity for Neutrino Mass Hierarchy Determination with Atmospheric Neutrinos

Despite the seemingly easy distinction of the two MH hypotheses for a given oscillation channel (as illustrated in Figure6), the distinction of the two hierarchy hypotheses with atmospheric neutrinos is challenging due to the cancellation of contributions from neutrinos and antineutrinos. Further cancellation effects occur between νe and νμ assuming an unavoidably imperfect flavour tagging method. The question of flavour tagging, that is, distinguishing cascade from track signatures and thereby isolating the νμ CC, signal has still to be improved—in ORCA as well as in PINGU. The reachable sensitivity for MH is further attenuated by the finite energy and angular resolutions of the considered detectors which are discussed above. Last but not least MH effects can be partly masked by readjusting oscillation parameters (in particular Δm322 and sin22θ23) within their present uncertainties.Both ORCA and PINGU have yet to provide a sensitivity calculation for MH, including all these effects plus detector and flux related sources of systematic uncertainties. The most detailed sensitivity study to date is based on input from PINGU performance figures and can be found in [42]. Its main result is shown in Figure 9. A 1.3 σ to 2.9 σ significance might be reached for three years running with the 20-lines PINGU detector depending on the true MH choice and the true values for other oscillation parameters. The result is based on a rather optimistic assumption of a cascade misidentification fraction of 5% and it assumes that the neutrino energy can be determined with a perfectly Gaussian resolution of 25%.Figure 9
Number ofσ for the hierarchy discovery of a 20-lines PINGU detector as a function of time for two extreme cases of the true parameter values (figure from [42]).The cancellation of matter effects between neutrinos and anti-neutrinos can be partly avoided by determining for eachνμ CC event the hadronic and muon part separately. This would allow to measure the inelasticity variable Bjorken-y which in turn could be used to statistically separate neutrinos from anti-neutrinos. The best separation would be reached for events of high inelasticity (large y), but these are hardly distinguishable from cascades due to the low-energy muon tracks. The potential of the y determination for the MH sensitivity is discussed in [43]. When assuming a realistic experimental smearing of the measured quantities, an enhancement of 25% is quoted in [43] by measuring hadrons and muons separately.The study of the neutrino mass hierarchy is not the only physics goal of PINGU and ORCA. The knowledge of the oscillation parametersΔm322 and sin2θ23 itself can be improved. The octant of the mixing angle can be determined, if its value is sufficiently distinct from maximal mixing and |Δm322| can be measured with a precision better than 1% based on one year of data with PINGU alone [41] under the assumption that systematic effects such as the uncertainty on the ice model do not deteriorate this measurement too strongly. Further, indirect dark matter searches from annihilation processes of gravitationally confined particles can be extended toward lower masses compared to current results from ANTARES and IceCube.


## 8. Neutrino Beams and Neutrino Telescopes

The above-mentioned complications of determining MH with atmospheric neutrinos are avoided by counting beam related events of a specific flavour. Neutrino beams can be rather pureνμ beams, that is, with low contamination from other flavours or polarities. Assuming a moderate capability of the target neutrino telescope to distinguish track signatures (νμ CC events or ντ CC events with subsequent muonic decay) from cascade signatures (everything else), two basic methods have been proposed to determine MH.


### 8.1. Muon Counting

Track signatures are isolated and beam related events are counted. This has been recently proposed [44] as a simple method to measure MH via the disappearance channel P(μ→μ). Assuming a conventional narrow band beam with a mean neutrino energy of 6–8 GeV, one can examine Figure 6 to determine the optimal baseline for such a measurement. The cleanest distinction of the two MH hypothesis is seen for cosΘ=0.6 or equivalently a baseline of 7645 km. The ignorance about the CP-phase does not affect the distinction as one can immediately verify from Figure 6. Taking into account locations of existing particle accelerators and existing or potential future neutrino telescopes, a beamline from Fermilab to a KM3Net site close to Sicily is proposed [44] which yields a baseline of 7800 km. For 1020 protons on target, about 1000 νμ CC events are expected in a Megaton detector such as ORCA for NH. This number is reduced by 30% for IH, providing a 9 σ separation of the two hypotheses on purely statistical ground.


### 8.2. Electron Counting

An alternative approach exploits the appearance channelP(νμ→νe) by counting cascade-like events [45]. Figure 10 shows the oscillation probability P(νμ→νe) for two very long baselines. The left plot corresponds to the distance Fermilab-Sicily (discussed above) the right plot is typical for a beamline from any accelerator in the Northern hemisphere to the South Pole, as discussed in [42]. In both cases, a clear separation of the hierarchy hypotheses is seen and it could be exploited in a counting experiment of cascade-like events. To find the optimal baseline for such a measurement, the dependence of the event rate N from the baseline L has to be taken into account: N∝L-2. For a typical neutrino beam, which provides neutrinos in the energy range 2–8 GeV, the event rate dependence from baseline is shown in Figure 11. An optimal separation between NH and IH is found for cosΘ=0.2 or equivalently a baseline of about 2600 km. A beam with such a baseline could be built from the Institute of High Energy Physics (IHEP) in Protvino near Moscow (located at 54∘52′N,37∘11′E [46]) towards one of the potential sites for the ORCA detector. The IHEP hosts the U70 proton accelerator [47] which provides protons with energies up to 70 GeV. It is operational since 1967 and it had been the world largest proton accelerator at its time of commissioning. For the calculation of event numbers below it is assumed that after an intensity upgrade, Npot=1.5·1021 can be provided within few years, similar to what is planned for the NOVA experiment [48]. If the ORCA detector would be built close to the ANTARES site (42∘48′N,6∘10′E), a neutrino beam from Protvino to this location would result in a baseline of 2588 km. An alternative site in the Ionian Sea off the Sicilian coast results incidentally in an identical baseline (within 1%) for a beam from Protvino. A moderate downward inclination of 11.7° is needed.Neutrino oscillation probabilitiesP(νμ→νe) parameters from a global fit as a function of neutrino energy for baselines of 7675 km (a) and 11468 km (b) which corresponds to cosΘ=0.6 and 0.9, respectively. The group of red lines is for NH, blue for IH and ϕCP varied in steps of 30° between 0°  and 330°.
(a)
(b)Event rates of oscillatedνe-CC events as function of cosΘ, which is proportional to the baseline. (a) Comparison of event rates for both hierarchy hypotheses (same color code as in earlier figures). (b) Event rate differences from left plot.
(a)
(b)The corresponding oscillation probabilitiesP(νμ→νe,νμ,ντ) with the best fit parameters from [32] are shown in Figure 12.  Most of the neutrino path will be in the outer Earth mantle with a maximal depth of 134 km. A constant density of 3.3 g/cm3 is used for the calculation of the oscillation probabilities as given in the PREM model [36] for the outer mantle.Neutrino oscillation probabilitiesP(νμ→νe), P(νμ→νμ) and P(νμ→ντ) for a baseline of 2600 km and oscillation parameters from a global fit as a function of neutrino energy. The group of red lines is for NH, blue for IH and ϕCP varied in steps of 30∘ between 0∘ and 330∘. The special values ϕCP=0∘,180∘ (CP conservation) are indicated in magenta, ϕCP=90∘,270∘ (maximal CP violation) in orange on the left plot.
(a)
(b)ForP(νμ→νe), a significant difference between both hierarchies is observed in the range 3GeV<Eν<8GeV. The variation of the CP-phase ϕCP (different lines of the same color) leads instead only to moderate changes of the oscillation probabilities. Peak values of 11% for P(νμ→νe) are observed in NH. This is smaller than that for larger baselines (see Figure 10), but fully compensated by the higher flux due to the smaller L (see Figure 11). P(νμ→νμ) depends only weakly on MH and not on ϕCP. The first vacuum oscillation maximum at 5 GeV is the dominating feature in the shown energy range. Counting νμ events can serve as a flux normalisation or it could be used to improve the measurement of the atmospheric oscillation parameters |Δm322| and θ23.To calculate event numbers, simplified cross-section formulas are used. The total cross sections forνμ and ν-μ CC interactions are taken from [4] in the parton scaling approximation:
(7)σνμCC(Eν)=0.68·(Eν/GeV)10-38cm2σν-μCC(Eν)=0.34·(Eν/GeV)10-38cm2.
Deviations from this linear behaviour due to quasielastic or resonant interactions are ignored. Their contribution would not alter the result of this study in a significant way. The relevant neutrino energies are significantly larger than mμ and me; therefore, from flavour universality, σνeCC=σνμCC and σν-eCC=σν-μCC. However for ντ CC interactions, the mass of the τ-lepton cannot be neglected. We use the calculation from [49]. For a threshold energy E0=5 GeV and E0<Eν<30 GeV, we find the following simple parametrisation:
(8)σντCC(Eν)=0.29log(EνE0)σνμCC(Eν).
The cross sections for ν(ν-) NC interactions for all flavours are approximated as 1/3σνeCC(Eν) and 1/3σν-eCC(Eν), respectively.As a first step, the impact of the MH on the total event count is evaluated, regardless of its topology in the detector. Whereas the sum of the oscillation probabilities∑α‍P(μ→α) equals unity due to the unitarity of the flavour mixing matrix, the cross section weighted sum
(9)Pμσ(Eν)=σνμNC(Eν)+∑αP(μ⟶α)σναCC(Eν)σνμNC(E0)+σνCC(E0)
(arbitrarily normalised at E0=1 GeV) may help identifying the optimal energy range to separate the two MH hypotheses. Pμσ is shown on the left plot of Figure 13 and it can be interpreted as the event rate per neutrino energy seen in a detector with an energy independent detection efficiency for a pure νμ flux flat in energy. No distinction is made neither between NC and CC interactions nor between different flavours. Nevertheless, a clear separation of the two MH hypotheses is observed for energies above 3 GeV. It can be attributed to the kinematical suppression of ντ CC interactions. The size of the effect is quantified in the right plot of Figure 13 which shows the ratio between IH and NH integrals from the left plot (for different values of ϕCP) taken between an assumed threshold energy of 2.5 GeV and a variable maximal energy Emax. An optimal value of Emax=6 GeV can be read from the figure, which yields a suppression of the IH event rate of 11–14% compared to the expected rate for NH. Extending the energy range to higher values reduces the relative size of the separation of the MH hypotheses. It can be concluded that a measurement of the mass hierarchy will be possible even without any flavour tagging capabilities using a neutrino beam in a limited energy range of 2–6 GeV and a large detector which can reliably count beam related neutrino interactions. Nevertheless, flavour tagging methods are discussed below and they are used to improve the significance of the measurement.(a) Summed oscillation probabilitiesP(μ→α) weighted by cross-sections (red: NH, blue: IH). (b) Ratio of integrals NIH/NNH from left plot between Emin=2.5 GeV and Emax.
(a)
(b)In the past, a neutrino beam was provided to several experiments, among them the bubble chamber SKAT [47]. Secondary hadrons were produced in an aluminum target and focused by parabolic lenses. Neutrinos were produced in a 140 m long decay line. The bubble chamber was situated 270 m behind the target and at a distance of lSKAT=245m downstream the beginning of the decay tube. This value will be used in the following to scale the beam intensity to the remote location. The neutrino fluxes dΦν/dEν, as they were delivered to the SKAT experiment for focusing of positively charged hadrons and for all four flavours present in the beam, are shown in Figure 14. A parametrisation of the beam intensity as function of the neutrino energy at the SKAT detector has been obtained from [47, 50] and it will be used in the following. As seen from Figure 14, a very clean νμ beam was provided with energies dominantly in the range 2–8 GeV. Contaminations from other flavours were on the subpercent level.Figure 14
Neutrino fluxdΦν/dEν per proton as seen by the SKAT bubble chamber. The different lines correspond to the different neutrino flavours.The rate of detected CC events of flavourα can now be calculated as follows:
(10)dNαdEν=Npot(lSKATlLBL)2Meff(Eν)mp×[σναCC(dΦνμdEνPμα+dΦνedEνPeα)+σν-αCC(dΦν-μdEνPμα¯+dΦν-edEνPeα¯)]
with mp the proton mass and the abbreviation Pβα=P(νβ→να) and Pβα¯=P(ν-β→ν-α) for the oscillation probabilities. The effective mass Meff(Eν) is taken from Section 7.1 and it is assumed to be flavour independent. This overestimates the relative ντ CC detection efficiency due to the escaping neutrino(s) of the tau decay. As ντ interactions are a background in the present analysis, this is a conservative approximation. For Npot=1.5·1021, one obtains dNe/dEν=1621±255 for NH but only dNe/dEν=497±100 for IH. The errors correspond to the uncertainty of the CP-phase. The statistical separation of both samples is better than 20 σ. dNμ/dEν depends only weakly on MH (10317 versus 10015) and it can serve as reference sample.The last step consists in adding background channels and considering misidentification for the two topologies. Contributions from atmospheric neutrinos and misreconstructed downgoing atmospheric muons can be ignored. A pulsed beam with a typical duty cycle lower than10-6 allows to safely discard these events. The beam itself is a source of background events. ντ CC events with a muonic τ-decay produce “track” events, whereas NC events and ντ CC events with a nonmuonic τ-decay have a genuine cascade signature. The rate of detected NC events is given by
(11)dNNCdEν=Npot(lSKATlLBL)2Meff(Eν/2)mp×[σνNC(dΦνμdEν+dΦνedEν)+σν-NC(dΦν-μdEν+dΦν-edEν)].
No dependence on oscillation parameters enters here. As, on average, 50% of the neutrino energy is transferred to the outgoing neutrino, the effective mass Meff is evaluated at Eν/2.In addition, the probabilityϵ(Eν) to misidentify a cascade as track is introduced, while the misidentification of a track as cascade will be called η(Eν). Both functions are parametrised by
(12)ϵ(Eν)=η(Eν)=1(Eν/GeV);Eν>2GeV.
For Eν=2 GeV, ϵ=η=0.5, which means that the two topologies cannot be distinguished and the attribution of an event to one of them is random. A 5 GeV neutrino produces a muon with an average range of 15 m in sea water which exceeds already by far the typical longitudinal size of a hadronic or electromagnetic shower. Correspondingly, the misidentification probability is assumed to drop to 20%. For Eν=10 GeV, ϵ and η further decrease to 10%. With these two quantities, the total background for the two event samples can be written as
(13)dNbgtrackdEν=ϵdNsigcascdEν+[ϵ(1-BRτμ)+(1-η)BRτμ]×dNτdEν+ϵdNNCdEνdNbgcascdEν=ηdNsigtrackdEν+[(1-ϵ)(1-BRτμ)+ηBRτμ]×dNτdEν+(1-ϵ)dNNCdEν.BRτμ stands for the muonic branching ratio of the tau decay (17.4%). The total number of observed events in each channel is given by adding the (reduced) number of signal events and the background contribution from (13) as follows:
(14)dNtottrackdEν=(1-η)dNsigtrackdEν+dNbgtrackdEνdNtotcascdEν=(1-ϵ)dNsigcascdEν+dNbgcascdEν.
The resulting rates are shown in Figure 15 and quoted in Table 1. The event rates between the two MH hypotheses differ by 9–18% with a statistical uncertainty of 1.2%. The statistical significance of the MH hypothesis test is still better than 7 σ and it remains at the level of 3 σ even when adding an additional systematic uncertainty of 3-4% (depending on the true value of the CP-phase) for the determination of the total cascade event rate.  Figure 15 illustrates the total event numbers as function of neutrino energy, detailing the different contributions. NC events are added to the Eν/2 bins according to their lower light yield in the detector. Their main contribution is found for energies below 6 GeV. Background from CC events is instead mainly seen above 5 GeV. The separation of the MH hypotheses is most pronounced in the range 4–8 GeV.Table 1
Event numbers forNpot=1.5·1021 in the track and cascade channel for both mass hierarchy schemes and varying ϕCP values.
Channel
Tracks NH
Tracks IH
Cascades NH
Cascades IH
No oscil
26315
—
Signal
8990
8735
1134–1547
350–519
Misreco
232–329
47–79
1326
1280
ν
τ
324–332
351–355
978–998
1057–1068
NC
1092
1092
3640
3640
BG Total
1655–1745
1494–1522
5944–5964
5977–5988
Total
10645–10736
10229–10257
7099–7491
6338–6496Observed events forNpot=1.5·1021 as function of the neutrino energy. (a) Track-like events. (b) Cascade-Like events. In green the background from NC events is given. The blue histogram shows the total background including misreconstructed CC events. Magenta histograms are for the total event rate (signal plus background) for NH (solid) and IH (dashed) and various values of ϕCP.
(a)
(b)


## 8.1. Muon Counting

Track signatures are isolated and beam related events are counted. This has been recently proposed [44] as a simple method to measure MH via the disappearance channel P(μ→μ). Assuming a conventional narrow band beam with a mean neutrino energy of 6–8 GeV, one can examine Figure 6 to determine the optimal baseline for such a measurement. The cleanest distinction of the two MH hypothesis is seen for cosΘ=0.6 or equivalently a baseline of 7645 km. The ignorance about the CP-phase does not affect the distinction as one can immediately verify from Figure 6. Taking into account locations of existing particle accelerators and existing or potential future neutrino telescopes, a beamline from Fermilab to a KM3Net site close to Sicily is proposed [44] which yields a baseline of 7800 km. For 1020 protons on target, about 1000 νμ CC events are expected in a Megaton detector such as ORCA for NH. This number is reduced by 30% for IH, providing a 9 σ separation of the two hypotheses on purely statistical ground.


## 8.2. Electron Counting

An alternative approach exploits the appearance channelP(νμ→νe) by counting cascade-like events [45]. Figure 10 shows the oscillation probability P(νμ→νe) for two very long baselines. The left plot corresponds to the distance Fermilab-Sicily (discussed above) the right plot is typical for a beamline from any accelerator in the Northern hemisphere to the South Pole, as discussed in [42]. In both cases, a clear separation of the hierarchy hypotheses is seen and it could be exploited in a counting experiment of cascade-like events. To find the optimal baseline for such a measurement, the dependence of the event rate N from the baseline L has to be taken into account: N∝L-2. For a typical neutrino beam, which provides neutrinos in the energy range 2–8 GeV, the event rate dependence from baseline is shown in Figure 11. An optimal separation between NH and IH is found for cosΘ=0.2 or equivalently a baseline of about 2600 km. A beam with such a baseline could be built from the Institute of High Energy Physics (IHEP) in Protvino near Moscow (located at 54∘52′N,37∘11′E [46]) towards one of the potential sites for the ORCA detector. The IHEP hosts the U70 proton accelerator [47] which provides protons with energies up to 70 GeV. It is operational since 1967 and it had been the world largest proton accelerator at its time of commissioning. For the calculation of event numbers below it is assumed that after an intensity upgrade, Npot=1.5·1021 can be provided within few years, similar to what is planned for the NOVA experiment [48]. If the ORCA detector would be built close to the ANTARES site (42∘48′N,6∘10′E), a neutrino beam from Protvino to this location would result in a baseline of 2588 km. An alternative site in the Ionian Sea off the Sicilian coast results incidentally in an identical baseline (within 1%) for a beam from Protvino. A moderate downward inclination of 11.7° is needed.Neutrino oscillation probabilitiesP(νμ→νe) parameters from a global fit as a function of neutrino energy for baselines of 7675 km (a) and 11468 km (b) which corresponds to cosΘ=0.6 and 0.9, respectively. The group of red lines is for NH, blue for IH and ϕCP varied in steps of 30° between 0°  and 330°.
(a)
(b)Event rates of oscillatedνe-CC events as function of cosΘ, which is proportional to the baseline. (a) Comparison of event rates for both hierarchy hypotheses (same color code as in earlier figures). (b) Event rate differences from left plot.
(a)
(b)The corresponding oscillation probabilitiesP(νμ→νe,νμ,ντ) with the best fit parameters from [32] are shown in Figure 12.  Most of the neutrino path will be in the outer Earth mantle with a maximal depth of 134 km. A constant density of 3.3 g/cm3 is used for the calculation of the oscillation probabilities as given in the PREM model [36] for the outer mantle.Neutrino oscillation probabilitiesP(νμ→νe), P(νμ→νμ) and P(νμ→ντ) for a baseline of 2600 km and oscillation parameters from a global fit as a function of neutrino energy. The group of red lines is for NH, blue for IH and ϕCP varied in steps of 30∘ between 0∘ and 330∘. The special values ϕCP=0∘,180∘ (CP conservation) are indicated in magenta, ϕCP=90∘,270∘ (maximal CP violation) in orange on the left plot.
(a)
(b)ForP(νμ→νe), a significant difference between both hierarchies is observed in the range 3GeV<Eν<8GeV. The variation of the CP-phase ϕCP (different lines of the same color) leads instead only to moderate changes of the oscillation probabilities. Peak values of 11% for P(νμ→νe) are observed in NH. This is smaller than that for larger baselines (see Figure 10), but fully compensated by the higher flux due to the smaller L (see Figure 11). P(νμ→νμ) depends only weakly on MH and not on ϕCP. The first vacuum oscillation maximum at 5 GeV is the dominating feature in the shown energy range. Counting νμ events can serve as a flux normalisation or it could be used to improve the measurement of the atmospheric oscillation parameters |Δm322| and θ23.To calculate event numbers, simplified cross-section formulas are used. The total cross sections forνμ and ν-μ CC interactions are taken from [4] in the parton scaling approximation:
(7)σνμCC(Eν)=0.68·(Eν/GeV)10-38cm2σν-μCC(Eν)=0.34·(Eν/GeV)10-38cm2.
Deviations from this linear behaviour due to quasielastic or resonant interactions are ignored. Their contribution would not alter the result of this study in a significant way. The relevant neutrino energies are significantly larger than mμ and me; therefore, from flavour universality, σνeCC=σνμCC and σν-eCC=σν-μCC. However for ντ CC interactions, the mass of the τ-lepton cannot be neglected. We use the calculation from [49]. For a threshold energy E0=5 GeV and E0<Eν<30 GeV, we find the following simple parametrisation:
(8)σντCC(Eν)=0.29log(EνE0)σνμCC(Eν).
The cross sections for ν(ν-) NC interactions for all flavours are approximated as 1/3σνeCC(Eν) and 1/3σν-eCC(Eν), respectively.As a first step, the impact of the MH on the total event count is evaluated, regardless of its topology in the detector. Whereas the sum of the oscillation probabilities∑α‍P(μ→α) equals unity due to the unitarity of the flavour mixing matrix, the cross section weighted sum
(9)Pμσ(Eν)=σνμNC(Eν)+∑αP(μ⟶α)σναCC(Eν)σνμNC(E0)+σνCC(E0)
(arbitrarily normalised at E0=1 GeV) may help identifying the optimal energy range to separate the two MH hypotheses. Pμσ is shown on the left plot of Figure 13 and it can be interpreted as the event rate per neutrino energy seen in a detector with an energy independent detection efficiency for a pure νμ flux flat in energy. No distinction is made neither between NC and CC interactions nor between different flavours. Nevertheless, a clear separation of the two MH hypotheses is observed for energies above 3 GeV. It can be attributed to the kinematical suppression of ντ CC interactions. The size of the effect is quantified in the right plot of Figure 13 which shows the ratio between IH and NH integrals from the left plot (for different values of ϕCP) taken between an assumed threshold energy of 2.5 GeV and a variable maximal energy Emax. An optimal value of Emax=6 GeV can be read from the figure, which yields a suppression of the IH event rate of 11–14% compared to the expected rate for NH. Extending the energy range to higher values reduces the relative size of the separation of the MH hypotheses. It can be concluded that a measurement of the mass hierarchy will be possible even without any flavour tagging capabilities using a neutrino beam in a limited energy range of 2–6 GeV and a large detector which can reliably count beam related neutrino interactions. Nevertheless, flavour tagging methods are discussed below and they are used to improve the significance of the measurement.(a) Summed oscillation probabilitiesP(μ→α) weighted by cross-sections (red: NH, blue: IH). (b) Ratio of integrals NIH/NNH from left plot between Emin=2.5 GeV and Emax.
(a)
(b)In the past, a neutrino beam was provided to several experiments, among them the bubble chamber SKAT [47]. Secondary hadrons were produced in an aluminum target and focused by parabolic lenses. Neutrinos were produced in a 140 m long decay line. The bubble chamber was situated 270 m behind the target and at a distance of lSKAT=245m downstream the beginning of the decay tube. This value will be used in the following to scale the beam intensity to the remote location. The neutrino fluxes dΦν/dEν, as they were delivered to the SKAT experiment for focusing of positively charged hadrons and for all four flavours present in the beam, are shown in Figure 14. A parametrisation of the beam intensity as function of the neutrino energy at the SKAT detector has been obtained from [47, 50] and it will be used in the following. As seen from Figure 14, a very clean νμ beam was provided with energies dominantly in the range 2–8 GeV. Contaminations from other flavours were on the subpercent level.Figure 14
Neutrino fluxdΦν/dEν per proton as seen by the SKAT bubble chamber. The different lines correspond to the different neutrino flavours.The rate of detected CC events of flavourα can now be calculated as follows:
(10)dNαdEν=Npot(lSKATlLBL)2Meff(Eν)mp×[σναCC(dΦνμdEνPμα+dΦνedEνPeα)+σν-αCC(dΦν-μdEνPμα¯+dΦν-edEνPeα¯)]
with mp the proton mass and the abbreviation Pβα=P(νβ→να) and Pβα¯=P(ν-β→ν-α) for the oscillation probabilities. The effective mass Meff(Eν) is taken from Section 7.1 and it is assumed to be flavour independent. This overestimates the relative ντ CC detection efficiency due to the escaping neutrino(s) of the tau decay. As ντ interactions are a background in the present analysis, this is a conservative approximation. For Npot=1.5·1021, one obtains dNe/dEν=1621±255 for NH but only dNe/dEν=497±100 for IH. The errors correspond to the uncertainty of the CP-phase. The statistical separation of both samples is better than 20 σ. dNμ/dEν depends only weakly on MH (10317 versus 10015) and it can serve as reference sample.The last step consists in adding background channels and considering misidentification for the two topologies. Contributions from atmospheric neutrinos and misreconstructed downgoing atmospheric muons can be ignored. A pulsed beam with a typical duty cycle lower than10-6 allows to safely discard these events. The beam itself is a source of background events. ντ CC events with a muonic τ-decay produce “track” events, whereas NC events and ντ CC events with a nonmuonic τ-decay have a genuine cascade signature. The rate of detected NC events is given by
(11)dNNCdEν=Npot(lSKATlLBL)2Meff(Eν/2)mp×[σνNC(dΦνμdEν+dΦνedEν)+σν-NC(dΦν-μdEν+dΦν-edEν)].
No dependence on oscillation parameters enters here. As, on average, 50% of the neutrino energy is transferred to the outgoing neutrino, the effective mass Meff is evaluated at Eν/2.In addition, the probabilityϵ(Eν) to misidentify a cascade as track is introduced, while the misidentification of a track as cascade will be called η(Eν). Both functions are parametrised by
(12)ϵ(Eν)=η(Eν)=1(Eν/GeV);Eν>2GeV.
For Eν=2 GeV, ϵ=η=0.5, which means that the two topologies cannot be distinguished and the attribution of an event to one of them is random. A 5 GeV neutrino produces a muon with an average range of 15 m in sea water which exceeds already by far the typical longitudinal size of a hadronic or electromagnetic shower. Correspondingly, the misidentification probability is assumed to drop to 20%. For Eν=10 GeV, ϵ and η further decrease to 10%. With these two quantities, the total background for the two event samples can be written as
(13)dNbgtrackdEν=ϵdNsigcascdEν+[ϵ(1-BRτμ)+(1-η)BRτμ]×dNτdEν+ϵdNNCdEνdNbgcascdEν=ηdNsigtrackdEν+[(1-ϵ)(1-BRτμ)+ηBRτμ]×dNτdEν+(1-ϵ)dNNCdEν.BRτμ stands for the muonic branching ratio of the tau decay (17.4%). The total number of observed events in each channel is given by adding the (reduced) number of signal events and the background contribution from (13) as follows:
(14)dNtottrackdEν=(1-η)dNsigtrackdEν+dNbgtrackdEνdNtotcascdEν=(1-ϵ)dNsigcascdEν+dNbgcascdEν.
The resulting rates are shown in Figure 15 and quoted in Table 1. The event rates between the two MH hypotheses differ by 9–18% with a statistical uncertainty of 1.2%. The statistical significance of the MH hypothesis test is still better than 7 σ and it remains at the level of 3 σ even when adding an additional systematic uncertainty of 3-4% (depending on the true value of the CP-phase) for the determination of the total cascade event rate.  Figure 15 illustrates the total event numbers as function of neutrino energy, detailing the different contributions. NC events are added to the Eν/2 bins according to their lower light yield in the detector. Their main contribution is found for energies below 6 GeV. Background from CC events is instead mainly seen above 5 GeV. The separation of the MH hypotheses is most pronounced in the range 4–8 GeV.Table 1
Event numbers forNpot=1.5·1021 in the track and cascade channel for both mass hierarchy schemes and varying ϕCP values.
Channel
Tracks NH
Tracks IH
Cascades NH
Cascades IH
No oscil
26315
—
Signal
8990
8735
1134–1547
350–519
Misreco
232–329
47–79
1326
1280
ν
τ
324–332
351–355
978–998
1057–1068
NC
1092
1092
3640
3640
BG Total
1655–1745
1494–1522
5944–5964
5977–5988
Total
10645–10736
10229–10257
7099–7491
6338–6496Observed events forNpot=1.5·1021 as function of the neutrino energy. (a) Track-like events. (b) Cascade-Like events. In green the background from NC events is given. The blue histogram shows the total background including misreconstructed CC events. Magenta histograms are for the total event rate (signal plus background) for NH (solid) and IH (dashed) and various values of ϕCP.
(a)
(b)


## 9. Conclusion

Neutrino telescopes have been designed to identify cosmic sources of high-energy neutrinos. Both ANTARES and IceCube have demonstrated their capability to contribute as well to the measurement of neutrino oscillation parameters. With several years of IceCube/DeepCore sensitivities comparable to MINOS might be reachable. Low-energy extensions such as ORCA and PINGU will allow to lower the effective detection threshold to values of few GeV. This will open the possibility to perform competitive measurements of neutrino oscillation parameters and to pin down yet unknown parameters such as the neutrino mass hierarchy or the octant of the mixing angleθ23. The significance of these future measurements with atmospheric neutrinos could be complemented by using accelerator neutrinos from yet to be built long-baseline beams pointing towards one of the future neutrino telescopes.



---

*Source: 782538-2013-11-27.xml*