Probing Particle Physics with IceCube

The IceCube observatory located at the South Pole is a cubic-kilometre optical Cherenkov telescope primarily designed for the detection of high-energy astrophysical neutrinos. IceCube became fully operational in 2010, after a seven-year construction phase, and reached a milestone in 2013 by the first observation of cosmic neutrinos in the TeV-PeV energy range. This observation does not only mark an important breakthrough in neutrino astronomy, but it also provides a new probe of particle physics related to neutrino production, mixing, and interaction. In this review we give an overview of the various possibilities how IceCube can address fundamental questions related to the phenomena of neutrino oscillations and interactions, the origin of dark matter, and the existence of exotic relic particles, like monopoles. We will summarize recent results and highlight future avenues.


Introduction
Not long after the discovery of the neutrino by Cowan and Reines in 1956 [1], the idea emerged that it represented the ideal astronomical messenger [2]. Neutrinos are only weakly interacting with matter and can cross cosmic distances without being absorbed or scattered. However, this weak interaction is also a challenge for the observation of these particles. Early estimates of the expected flux of high-energy neutrinos associated with the observed flux of extra-galactic cosmic rays indicated that neutrino observatories require gigaton masses as a necessary condition to observe a few neutrino interactions per year [3]. These requirements can only be met by special a e-mail: markus.ahlers@nbi.ku.dk b e-mail: helbing@uni-wuppertal.de c e-mail: cph@physics.uu.se experimental setups that utilise natural resources. Not only that -the detector material has to be suitable so that these few interactions can be made visible and separated from large atmospheric backgrounds.
Despite these obstacles, there exist a variety of experimental concepts to detect high-energy neutrinos. One particularly effective method is based on detecting the radiation of optical Cherenkov light produced by relativistic charged particles. This requires the use of optically transparent detector media like water or ice, where the Cherenkov emission can be read out by optical sensors deployed in the medium. This information then allows to reconstruct the various Cherenkov light patterns produced in neutrino events and infer the neutrino flavour, arrival direction, and energy. The most valuable type of events for neutrino astronomy are charged current interactions of muon-neutrinos with matter near the detector. These events produce muons that can range into the detector and allow the determination of the initial muon-neutrino direction within a precision of better than one degree.
Presently the largest optical Cherenkov telescope is the IceCube Observatory, which uses the deep Glacial Ice at the geographic South Pole as its detector medium. The principal challenge of any neutrino telescope is the large background of atmospheric muons and neutrinos produced in cosmic ray interactions in the atmosphere. High-energy muons produced in the atmosphere have a limited range in ice and bedrock. Nevertheless IceCube, at a depth of 1.5 kilometres, observes about 100 billion atmospheric muon events per year. This large background can be drastically reduced by only looking for up-going events, i.e., events that originate below the horizon. This cut leaves only muons produced by atmospheric neutrinos at a rate of about 100,000 per year. While these large backgrounds are an obstacle for neutrino astronomy they provide a valuable probe for cosmic ray physics in general and for neutrino oscillation and interaction studies in particular.
In this review we want to highlight IceCube's potential as a facility to probe fundamental physics. There 2 exist a variety of methods to test properties of the Standard Model (SM) and its possible extensions. The flux of atmospheric and astrophysical neutrinos observed in IceCube allows to probe fundamental properties in the neutrino sector related to the standard neutrino oscillations (neutrino mass differences, mass ordering, and flavour mixing) and neutrino-matter interactions. It also provides a probe for exotic oscillation effects, e.g., related to the presence of sterile neutrinos or non-standard neutrino interactions with matter. The ultra-long baselines associated with the propagation of cosmic neutrinos observed beyond 10 TeV allow for various tests of feeble neutrino oscillation effects that can leave imprints on the oscillation-averaged flavour composition.
One of the fundamental questions in cosmology is the origin of dark matter that today constitutes one quarter of the total energy density of the Universe. Candidate particles for this form of matter include weakly interacting massive particles (WIMPs) that could have been thermally produced in the early Universe. IceCube can probe the existence of these particles by the observation of a flux of neutrinos produced in the annihilation or decay of WIMPs gravitationally clustered in nearby galaxies, the halo of the Milky Way, the Sun, or the Earth. In the case of compact objects, like Sun and Earth, neutrinos are the only SM particles that can escape the dense environments to probe the existence of WIMPs.
Neutrino telescopes can also probe exotic particles leaving direct or indirect Cherenkov signals during their passage through the detector. One important example are relic magnetic monopoles, topological defects that could have formed during a phase transition in the early Universe. Light exotic particles associated with extensions of the Standard Model can also be produced by the interactions of high-energy neutrinos or cosmic rays. Collisions of neutrinos and cosmic rays with nucleons in the vicinity of the Cherenkov detector can reach center-ofmass energies of the order √ s 1 TeV (neutrino energy E ν 10 15 eV) or even √ s 100 TeV (cosmic ray energy E CR 10 20 eV), respectively, only marginally probed by collider experiments. The outline of this review is as follows. We will start in sections 2 and 3 with a description of the IceCube detector, atmospheric backgrounds, standard event reconstructions, and event selections. In section 4 we summarise the phenomenology of three-flavour neutrino oscillation and IceCube's contribution to test the atmospheric neutrino mixing. We will cover standard model neutrino interactions in section 5 and highlight recent measurements of the inelastic neutrino-nucleon cross sections with Ice-Cube. We then move on to discuss IceCube's potential to probe non-standard neutrino oscillation with atmospheric and astrophysical neutrino fluxes in section 6. In section 7 we highlight IceCube results on searches for dark matter and section 8 is devoted to magnetic monopoles while section 9 covers other massive exotic particles and Big Bang relics.
Any review has its limitations, both in scope and timing. We have given priority to present a comprehensive view of the activity of IceCube in areas related to the topic of this review, rather than concentrating on a few recent results. We have also chosen at times to include older results for completeness, or when it was justified as an illustration of the capabilities of the detector on a given topic. The writing of any review develops along its own plot and updated results on some analyses have been made public while this paper was in preparation, and could not be included here. This only reflects on the lively activity of the field.

The IceCube Neutrino Observatory
The IceCube Neutrino Observatory [4] consists of an inice array (simply "IceCube" hereafter) and a surface air shower array, IceTop [5]. IceCube utilises one cubic kilometre of the deep ultra-clear glacial ice at the South Pole as its detector medium (see left panel of Fig. 1). This volume is instrumented with 5,160 Digital Optical Modules (DOMs) that register the Cherenkov light emitted by relativistic charged particles passing through the detector. The DOMs are distributed on 86 read-out and support cables ("strings") and are deployed between 1.5 km and 2.5 km below the surface. Most strings follow a triangular grid with a width of 125 m, evenly spaced over the volume (see green markers in right panel of Fig. 1).
Eight strings are placed in the centre of the array and are instrumented with a denser DOM spacing and typical inter-string separation of 55 m (red markers in right panel of Fig. 1). They are equipped with photomultiplier tubes with higher quantum efficiency. These strings, along with the first layer of the surrounding standard strings, form the DeepCore low-energy sub-array [6]. Its footprint is depicted by a blue dashed line in Fig. 1. While the original IceCube array has a neutrino energy threshold of about 100 GeV, the addition of the denser infill lowers the energy threshold to about 10 GeV. The DOMs are operated to trigger on single photo-electrons and to digitise in-situ the arrival time of charge ("waveforms") detected in the photomultiplier. The dark noise rate of the DOMs is about 500 Hz for standard modules and 800 Hz for the highquantum-efficiency DOMs in the DeepCore sub-array.
Some results highlighted in this review were derived from data collected with the AMANDA array [7], the predecessor of IceCube built between 1995 and 2001 at the same site, and in operation until May 2009. AMANDA was not only a proof of concept and a hardware test-bed for the IceCube technology, but a full fledged detector which obtained prime results in the field.

Neutrino Event Signatures
As we already highlighted in the introduction, the main event type utilised in high-energy neutrino astronomy are charged current (CC) interactions of muon neutrinos with nucleons (N), ν µ + N → µ − + X. These interactions produce high-energy muons that lose energy by ionisation, bremsstrahlung, pair production and photo-nuclear interactions in the ice [8]. The combined Cherenkov light from the primary muon and secondary relativistic charged particles leaves a track-like pattern as the muon passes through the detector. An example is shown in the left panel of Fig. 2. In this figure, the arrival time of Cherenkov light in individual DOMs is indicated by colour (earlier in red and later in blue) and the size of each DOM is proportional to the total Cherenkov light it detected. 1 . Since the average scattering angle between the incoming neutrino and the outgoing muon decreases with energy, Ψ ν→µ ∼ 0.7°(E ν /TeV) −0.7 [9], an angular resolution below 1°can be achieved for neutrinos with energies above a few TeV, only limited by the detector's intrinsic angular resolution. This changes at low energies, where muon tracks are short and their angular resolution deteriorates rapidly. For neutrino energies of a few tens of GeVs the angular resolution reaches a median of ∼ 40°.
All deep-inelastic interactions of neutrinos, both neutral current (NC), ν α + N → ν α + X and charged current, ν α + N → − α + X, create hadronic cascades X that are visible by the Cherenkov emission of secondary charged particles. However, these secondaries can not produce elongated tracks in the detector due to their rapid scattering or decay in the medium. Because of the large separation of the strings in IceCube and the scattering of light in the ice, the Cherenkov light distribution from particle cascades in the detector is rather spherical, see right panel of Fig. 2. For cascades or tracks fully contained in the detector, the energy resolution is significantly better since the full energy is deposited in the detector and it is proportional to the detected light. The ability to distinguish these two light patterns in any energy range is crucial, since cascades or tracks can contribute to background or signal depending on the analysis performed.
The electrons produced in charged current interactions of electron neutrinos, ν e + N → e − + X, will contribute to an electromagnetic cascade that overlaps with the hadronic cascade X at the vertex. At energies of E ν 6.3 PeV, electron anti-neutrinos can interact resonantly with electrons in the ice via a W-resonance ("Glashow" resonance) [10]. The W-boson decays either into hadronic states with a branching ratio (BR) of 67%, or into leptonic states (BR 11% for each flavour). This type of event can be visible by the appearance of isolated muon tracks starting in the detector or by spectral features in the event distribution [11].
Also the case of charged current interactions of tau neutrinos, ν τ + N → τ + X, is special. Again, the hadronic cascade X is visible in Cherenkov light. The tau has a lifetime (at rest) of 0.29 ps and decays to leptons as τ − → µ − + ν µ + ν τ (BR 18%) and τ − → e − + ν e + ν τ (BR 18%) or to hadrons (mainly pions and kaons, BR 64%) as τ − → ν τ + mesons. With tau energies below 100 TeV these charged current events will also contribute to track and cascade events. However, the delayed decay of taus at higher energies can become visible in Ice-Cube, in particular above around a PeV when the decay length becomes of the order of 50 m. This allows for a variety of characteristic event signatures, depending on the tau energy and decay channel [12,13]. The left plot shows a muon track from a ν µ interaction crossing the detector. Each coloured dot represents a hit DOM. The size of the dot is proportional to the amount of light detected and the colour code is related to the relative timing of light detection: read denotes earlier hits, blue, later hits. The right plot shows a ν e or ν τ charged-current (or any flavour neutral-current) interaction inside the detector.

Event Selection and Reconstruction
In this review we present results from analyses which use different techniques tailored to the characteristics of the signals searched for. It is therefore impossible to give a description of a generic analysis strategy which would cover all aspects of every approach. There are, however, certain levels of data treatment and analysis techniques that are common for all analyses in IceCube, and which we cover in this section.

Event Selection
Several triggers are active in IceCube in order to preselect potentially interesting physics events [14]. They are based on finding causally connected spatial hit distributions in the array, typically requiring a few neighbour or next-to-neighbour DOMs to fire within a predefined time window. Most of the triggers aim at finding relativistic particles crossing the detector and use time windows of the order of a few microseconds. In order to extend the reach of the detector to exotic particles, e.g., monopoles catalysing nucleon-decay, which can induce events lasting up to milliseconds, a dedicated trigger sensitive to non-relativistic particles with velocities down to β −4 has also been implemented.
When a trigger condition is fulfilled the full detector is read out. IceCube triggers at a rate of 2.5 kHz, collecting about 1 TB/day of raw data. To reduce this amount of data to a more manageable level, a series of software filters are applied to the triggered events: fast reconstructions [15] are performed on the data and a first event selection carried out, reducing the data stream to about 100 GB/day. These reconstructions are based on the position and time of the hits in the detector, but do not include information about the optical properties of the ice, in order to speed up the computation. The filtered data is transmitted via satellite to several IceCube institutions in the North for further processing.
Offline processing aims at selecting events according to type (tracks or cascades), energy, or specific arrival directions using sophisticated likelihood-based reconstructions [16,17]. These reconstructions maximise the likelihood function built from the probability of obtaining the actual temporal and spatial information in each DOM ("hit") given a set of track parameters (vertex, time, energy, and direction). For low-energy events, where the event signature is contained within the volume of the detector, a joint fit of muon track and an hadronic cascade at the interaction vertex is performed. For those events the total energy can be reconstructed with rather good accuracy, depending on further details of the analysis. Typically, more than one reconstruction is performed for each event. This allows, for example, to estimate the probability of each event to be either a track or a cascade. Each analysis will then use complex classification methods based on machine-learning techniques to further separate a possible signal from the background. Variables that describe the quality of the reconstructions, the time development and the spatial distribution of hit DOMs in the detector are usually used in the event selection.

Effective Area and Volume
After the analysis-dependent event reconstruction and selection, the observed event distribution in energy and arrival direction can be compared to the sum of background and signal events. For a given neutrino flux, φ ν , the total number of signal events, µ s , expected at the detector can be expressed as where T is the exposure time and A eff ν α the detector effective area for neutrino flavour α. The effective area encodes the trigger and analysis efficiencies and depends on the observation angle and neutrino energy.
In practice, the figure of merit of a neutrino telescope is the effective volume, V eff , the equivalent volume of a detector with 100% detection efficiency of neutrino events. 5 This quantity is related to the signal events as where φ is the neutrino flux after taking into account Earth absorption and regeneration effects, n is the local target density, and σ the neutrino cross section for the relevant neutrino signal. The effective volume allows to express the event number by the local density of events, i.e., the quantity within [·]. This definition has the practical advantage that the effective volume can be simulated from a uniform distribution of neutrino events: if n gen (E ν , Ω) is the number of Monte-Carlo events generated over a large geometrical generation volume V gen by neutrinos with energy E ν injected into the direction Ω, then the effective volume is given by where n s is the number of remaining signal events after all the selection cuts of a given analysis.

Background Rejection
There are two backgrounds in any analysis with a neutrino telescope: atmospheric muons and atmospheric neutrinos, both produced in cosmic-ray interactions in the atmosphere. The atmospheric muon background measured by IceCube [18] is much more copious than the atmospheric neutrino flux, by a factor up to 10 6 depending on declination. Note that cosmic ray interactions can produce several coincident forward muons ("muon bundle") which are part of the atmospheric muon background. Muon bundles can be easily identified as background in some cases, but they can also mimic bright single tracks (like monopoles for example) and are more difficult to separate from the signal in that case. Even if many of the IceCube analyses measure the atmospheric muon background from the data, the CORSIKA package [19] is generally used to generate samples of atmospheric muons that are used to cross-validate certain steps of the analyses. The large background of atmospheric muons can be efficiently reduced by using the Earth as a filter, i.e., by selecting up-going track events, at the expense of reducing the sky coverage of the detector to the Northern Hemisphere (see Fig. 3). Still, due to light scattering in the ice and the emission angle of the Cherenkov cone, a fraction of the down-going atmospheric muon tracks can be misreconstructed as up-going through the detector. This typically leads to a mismatch between the predicted atmospheric neutrino rate and the data rate at the final level of many analyses. There are analyses where a certain atmospheric muon contamination can be tolerated and it does not affect the final result. These are searches that look for a difference in the shape of the energy and/or angular  spectra of the signal with respect to the background, and are less sensitive to the absolute normalisation of the latter. For others, like monopole searches, misreconstructed atmospheric muons can reduce the sensitivity of the detector. We will describe in more detail how each analysis deals with this background when we touch upon specific analyses in the rest of this review. The atmospheric neutrino flux constitutes an irreducible background for any search in IceCube, and sets the baseline to define a discovery in many analyses. It is therefore crucial to understand it both quantitatively and qualitatively. The flux of atmospheric neutrinos is dominated by the production and decay of mesons produced by cosmic ray interactions with air molecules [20]. The behaviour of the neutrino spectra can be understood from the competition of meson (m) production and decay in the atmosphere: At high energy, where the meson decay rate is much smaller than the production rate, the meson flux is calorimetric and simply follows the cosmic ray spectrum, Φ m ∝ E −Γ . Below a critical energy m , where the decay rate becomes comparable to the production rate, the spectrum becomes harder by one power of energy, Φ m ∝ E 1−Γ . The corresponding neutrino spectra from the decay of mesons are softer by one power of energy, Φ ν ∝ Φ m /E due to the energy dependence of the meson decay rate.
The neutrino flux arising from pion and kaon decay is reasonably well understood, with an uncertainty in the range 10%-20% [20]. Figure 4 shows the atmospheric 6 10 10 2 10 3 10 4 10 5 10 6 [21,22] and ν µ + ν µ [23] spectra. The green data show the inferred bin-wise spectrum of the four-year highenergy starting event (HESE) analysis [24,25]. The green line and green-shaded area indicate the best-fit and 1σ uncertainty range of a power-law fit to the HESE data. Note that the HESE analysis vetoes atmospheric neutrinos, and the true background level is much lower as indicated in the plot. The green line and green-shaded area indicate the best-fit and 1σ uncertainty range of a power-law fit of the up-going muon neutrino analysis [26]. neutrino fluxes measured by IceCube. The atmospheric muon neutrino spectrum (ν µ + ν µ ) was obtained from one year of IceCube data (April 2008 to May 2009) using up-going muon tracks [23]. The atmospheric electron neutrino spectra (ν µ + ν µ ) were analysed by looking for contained cascades observed with the low-energy infill array DeepCore between June 2010 and May 2011 in the energy range from 80 GeV to 6 TeV [21]. This agrees well with a more recent analysis using contained events observed in the full IceCube detector between May 2011 and May 2012 with an extended energy range from 100 GeV to 100 TeV [22]. All measurements agree well with model prediction of "conventional" atmospheric neutrinos produced in pion and kaon decay. IceCube uses the public Monte Carlo software GENIE [27] and the internal software NUGEN (based on [28]) to generate samples of atmospheric neutrinos for its analyses, following the flux described in [29].
Kaons with an energy above 1 TeV are also significantly attenuated before decaying and the "prompt" component, arising mainly from very short-lived charmed mesons (D ± , D 0 , D s and Λ c ) is expected to dominate the spectrum. The prompt atmospheric neutrino flux, however, is much less understood, because of the uncertainty on the cosmic ray composition and relatively poor knowledge of QCD processes at small Bjorken-x [30,31,32,33,34]. In IceCube analyses the normalisation of the prompt atmospheric neutrino spectrum is usually treated as a nuisance parameter, while the energy distributions follows the model prediction of Ref. [30].
For high enough neutrino energies (O(10) TeV), the possibility exists of rejecting atmospheric neutrinos by selecting starting events, where an outer layer of DOMs acts as a virtual veto region for the neutrino interaction vertex. This technique relies on the fact that atmospheric neutrinos are accompanied by muons produced in the same air shower, that would trigger the veto [35,36]. The price to pay is a reduced effective volume of the detector for down-going events and a different sensitivity for upgoing and down-going events. This approach has been extremely successful, extending the sensitivity of IceCube to the Southern Hemisphere including the Galactic centre. There is not a generic veto region defined for all Ice-Cube analyses, but each analysis finds its optimal definition depending on its physics goal. Events that present more than a predefined number of hits within some time window in the strings included in the definition of the veto volume are rejected. A reduction of the atmospheric muon background by more than 99%, depending on analysis, can be achieved in this way (see for example [36,26]).
This approach has been also the driver behind one of the most exciting recent results in multi-messenger astronomy: the first observation of high-energy astrophysical neutrinos by IceCube. The first evidence of this flux could be identified from a high-energy starting event (HESE) analysis, with only two years of collected data in 2013 [37,24,25]. The event sample is dominated by cascade events, with only a rather poor angular resolution of about 10°. The result is consistent with an excess of events above the atmospheric neutrino background observed in up-going muon tracks from the Northern Hemisphere [38,26]. Figure 4 summarises the neutrino spectra inferred from these analyses. Based on different methods for reconstruction and energy measurement, their results agree, pointing at extra-galactic sources whose flux has equilibrated in the three flavours after propagation over cosmic distances [39] with ν e : ν µ : ν τ ∼ 1 : 1 : 1. While both types of analyses have now reached a significance of more than 5σ for an astrophysical neutrino flux, the origin of this neutrino emission remains a mystery (see, e.g., Ref. [40]).

Standard Neutrino Oscillations
Over the past decades, experimental evidence for neutrino flavour oscillations has been accumulating in solar (ν e ), atmospheric (ν e,µ & ν e,µ ), reactor (ν e ), and accelerator (ν µ & ν µ ) neutrino data (for a review see [41]). These oscillation patterns can be convincingly interpreted as a nontrivial mixing of neutrino flavour and mass states with a small solar and large atmospheric mass splitting. Neutrinos ν α with flavour α = e, µ, τ refer to those neutrinos that 7 couple to leptons α in weak interactions. Flavour oscillations are based on the effect that these flavour states are a non-trivial superposition of neutrino mass eigenstates ν j (j = 1, 2, 3) expressed as where the U αj 's are elements of the unitary neutrino massto-flavour mixing matrix, the so-called Pontecorvo-Maki-Nagakawa-Sakata (PMNS) matrix [42,43,44]. In general, the mixing matrix U has nine degrees of freedom, which can be reduced to six by absorbing three global phases into the flavour states ν α . The neutrino mixing matrix U is then conveniently parametrised [41] by three Euler rotations θ 12 , θ 23 , and θ 13 , and three CP-violating phases δ, α 1 and α 2 , Here, we have made use of the abbreviations sin θ ij = s ij and cos θ ij = c ij . The phases α 1/2 are called Majorana phases, since they have physical consequences only if the neutrinos are Majorana spinors, i.e., their own antiparticles. Note, that the phase δ (Dirac phase) appears only in combination with non-vanishing mixing sin θ 13 . Neutrino oscillations can be derived from plane-wave solutions of the Hamiltonian, that coincide with mass eigenstates in vacuum, exp(−i(ET − pL)). To leading order in m/E, the neutrino momentum is p E − m 2 /(2E) and a wave packet will travel a distance L T. Therefore, the leading order phase of the neutrino at distance L from its origin is exp(−im 2 L/(2E)). From this expression we see that the effect of neutrino oscillations depend on the difference of neutrino masses, ∆m 2 ij ≡ m 2 i − m 2 j . After traveling a distance L an initial state ν α becomes a superposition of all flavours, with probability of transition to flavour β given by P ν α →ν β = | ν β |ν α | 2 . This can be expressed in terms of the PMNS matrix elements as [41] where the oscillation phase ∆ ij can be parametrised as Note, that the third term in Eq. (6) comprises CP-violating effects, i.e., this term can change sign for the process P ν α →ν β , corresponding to the exchange U ↔ U * in Eq. (6). For the standard parametrisation (5) the single CP-violating contribution can be identified as the Dirac phase δ; oscillation experiments are not sensitive to Majorana phases.
The first compelling evidence for the phenomenon of atmospheric neutrino oscillations was observed with Super-Kamiokande (SK) [45]. The simplest and most direct interpretation of the atmospheric data [46] is oscillations of muon neutrinos, most likely converting into tau neutrinos. The survival probability of ν µ can be approximated as an effective two-level system with P ν µ →ν µ = 1 − sin 2 2θ atm sin 2 ∆ atm (8) The angular distribution of contained events in SK shows that for E ν ∼ 1 GeV, the deficit comes mainly from L atm ∼ 10 2 − 10 4 km. The corresponding oscillation phase must be nearly maximal, ∆ atm ∼ 1, which requires a mass splitting ∆m 2 atm ∼ 10 −4 − 10 −2 eV 2 . Moreover, assuming that all up-going ν µ 's which would yield multi-GeV events oscillate into a different flavour while none of the down-going ones do, the observed up-down asymmetry leads to a mixing angle very close to maximal, sin 2 2θ atm > 0.92 at 90%C.L. These results were later confirmed by the KEK-to-Kamioka (K2K) [47] and the Main Injector Neutrino Oscillation Search (MINOS) [48] experiments, which observed the disappearance of accelerator ν µ 's at a distance of 250 km and 735 km, respectively, as a distortion of the measured energy spectrum.
Furthermore, solar neutrino data collected by SK [49], the Sudbury Neutrino Observatory (SNO) [50] and Borexino [51] show that solar ν e 's produced in nuclear processes convert to ν µ or ν τ . For the interpretation of solar neutrino data it is crucial to account for matter effects that can have a drastic effect on the neutrino flavour evolution. The coherent scattering of electron neutrinos off background electrons with a density N e introduces a unique 2 potential term V mat = √ 2G F N e , where G F is the Fermi constant [52]. In the effective two-level system for the survival of electron neutrinos, the effective matter oscillation parameters (∆m 2 eff & θ eff ) relate to the vacuum values (∆m 2 & θ ) as tan 2θ eff tan 2θ where the resonance density is given by The effective oscillation parameters in the case of electron anti-neutrinos are the same as (9) and (10) Table 1 Results of a global analysis [58] of mass splittings, mixing angles, and Dirac phase for normal and inverted mass ordering. We best-fit parameters are shown with 1σ uncertainty.
The previous mixing and oscillation parameters are derived under the assumption of a constant electron density N e . If the electron density along the neutrino trajectory is only changing slowly compared to the effective oscillation frequency, the effective mass eigenstates will change adiabatically. Note that the oscillation frequency and oscillation depth in matter exhibits a resonant behaviour [53,54,52]. This Mikheyev-Smirnov-Wolfenstein (MSW) resonance can have an effect on continuous neutrino spectra, but also on monochromatic neutrinos passing through matter with slowly changing electron densities, like the radial density gradient of the Sun. Once these matter effect is taken into account, the observed intensity of solar electron neutrinos at different energies compared to theoretical predictions can be used to extract the solar mixing parameters. In addition to solar neutrino experiments, the KamLAND Collaboration [55] has measured the flux of ν e from distant reactors and find that ν e 's disappear over distances of about 180 km. This observation allows a precise determination of the solar mass splitting ∆m 2 consistent with solar data.
The results obtained by short-baseline reactor neutrino experiments show that the remaining mixing angle θ 13 is small. This allows to identify the mixing angle θ 12 as the solar mixing angle θ and θ 23 as the atmospheric mixing angle θ atm . Correspondingly, the mass splitting can be identified as ∆m 2 ∆m 2 21 and ∆m 2 atm |∆m 2 32 | |∆m 2 31 |. However, observations by the reactor neutrino experiments Daya-Bay [56] and RENO [57] show that the small reactor neutrino mixing angle θ 13 is larger than zero. As pointed out earlier, this is important for the observation of CP-violating effects parametrised by the Dirac phase δ in the PMNS matrix (5).
The global fit to neutrino oscillation data is presently incapable to determine the ordering of neutrino mass states. The fit to the data can be carried out under the assumption of normal (m 1 < m 2 < m 3 ) or inverted (m 3 < m 1 < m 2 ) mass ordering. A recent combined analysis [58] of solar, atmospheric, reactor, and accelerator neutrino data gives the values for the mass splittings, mixing angles, and CP-violating Dirac phase for normal or inverted mass ordering shown in Table 1. Note that, presently, the Dirac phase is inconsistent with δ = 0 at the 3σ level, independent of mass ordering. Ref. [64] Neutrino oscillation measurements are only sensitive to the relative neutrino mass differences. The absolute neutrino mass scale can be measured by studying the electron spectrum of tritium ( 3 H) β-decay. Present upper limits (95% C.L.) on the (effective) electron anti-neutrino mass are at the level of m ν e < 2 eV [59,60]. The KATRIN experiment [61] is expected to reach a sensitivity of m ν e < 0.2 eV. Neutrino masses are also constrained by their effect on the expansion history of the Universe and the formation of large-scale structure. Assuming standard cosmology dominated at late times by dark matter and dark energy, the upper limit (95% C.L.) on the combined neutrino masses is ∑ i m i < 0.23 eV [62].
The mechanism that provides neutrinos with their small masses is unknown. The existence of right-handed neutrino fields, ν R , would allow to introduce a Dirac mass term of the form m D ν L ν R + h.c., after electroweak symmetry breaking. Such states would be "neutral" with respect to the standard model gauge interactions, and therefore sterile [44]. However, the smallness of the neutrino masses would require unnaturally small Yukawa couplings. This can be remedied in seesaw models (see, e.g., Ref. [63]). Being electrically neutral, neutrinos can be Majorana spinors, i.e., spinors that are identical to their charge-conjugate state, ψ c ≡ Cψ T , where C is the charge-conjugation matrix. In this case, we can introduce Majorana mass terms of the form m L ν L ν c L /2 + h.c. and the analogous term for ν R . In seesaw models the individual size of the mass terms are such that m L 0 and m D m R . After diagonalization of the neutrino mass matrix, the masses of active neutrinos are then proportional to m i m 2 D /m R . This would explain the smallness of the effective neutrino masses via a heavy sector of particles beyond the Standard Model.

Atmospheric Neutrino Oscillations with IceCube
The atmospheric neutrino "beam" that reaches IceCube allows to perform high-statistics studies of neutrino oscillations at higher energies, and therefore is subject to  [65]. Systematic uncertainties are split into the normalisation contribution (dashed areas) and the shape contribution (filled areas) for each assumption shown. Right Panel: Significance contours at 68% and 90% C.L. for the best-fit values of the IceCube analysis (red curves), compared with results of the ANTARES [66], MINOS [67] and Super-Kamiokande [68] experiments. Figures reprinted with permission from Ref. [69] (Copyright 2013 APS) different systematic uncertainties, than those typically available in reactor-or accelerator-based experiments. Atmospheric neutrinos arrive at the detector from all directions, i.e., from travelling more than 12,700 km (vertically up-going) to about 10 km (vertically down-going), see Fig. 3. The path length from the production point in the atmosphere to the detector is therefore related to the measured zenith angle θ zen . Combined with a measurement of the neutrino energy, this opens the possibility of measuring ν µ disappearance due to oscillations, exploiting the dependence of the disappearance probability with energy and arrival angle.
Although the three neutrino flavours play a role in the oscillation process, a two-flavour approximation as in Eq. (8) is usually accurate to the percent level with ∆ atm ∆ 23 and θ atm θ 23 . The survival probability of muon neutrinos as a function of path length through the Earth and neutrino energy is shown in Fig. 5. It can be seen that, for the largest distance travelled by atmospheric neutrinos (the diameter of the Earth), Eq. (8) shows a maximum ν µ disappearance at about 25 GeV. This is precisely within the energy range of contained events in DeepCore. Simulations show that the neutrino energy response of DeepCore spans from about 6 GeV to about 60 GeV, peaking at 30 GeV. Neutrinos with higher energies will produce muon tracks that are no longer contained in the DeepCore volume.
Given this relatively narrow energy response of Deep-Core compared with the wide range of path lengths, it is possible to perform a search for ν µ disappearance through a measurement of the rate of contained events as a function of arrival direction, even without a precise energy determination. This is the approach taken in Ref. [69]. Events starting in DeepCore were selected by using the rest of the IceCube strings as a veto. A "high-energy" sample of events not contained in DeepCore was used as a reference, since ν µ disappearance due to oscillations at higher energies (O(100) GeV) is not expected. The atmospheric muon background is reduced to a negligible level by removing tracks that enter the DeepCore fiducial volume from outside, and by only considering up-going events, i.e., events that have crossed the Earth (cos θ zen ≤ 0), although a contamination of about 10%-15% of ν e events misidentified as tracks remained, as well as ν τ from ν µ oscillations. These two effects were included as background.
After all analysis cuts, a high-purity sample of 719 events contained in DeepCore were detected in a year. The left panel of Fig. 6 shows the angular distribution of the remaining events compared with the expected event rate without oscillations (red-shaded area) and with oscillations using current world-average values for sin 2 θ 23 and |∆m 2 32 | [65] (grey-shaded area). A statistically significant deficit of events with respect to the nonoscillation scenario can be seen near the vertical direction (−0.6 < cos θ zen < −1.0), while no discrepancy was observed in the reference high-energy sample (see Fig. 2 in [69]). The discrepancy between the data and the nonoscillation case can be used to fit the oscillation parameters, without assuming any a priori value for them. The right panel of Fig. 6 shows the result of that fit, with 68% (1σ) and 90% contours around the best-fit values found: sin 2 (2θ 23 ) = 1 and |∆m 2 32 | = 2.3 +0.6 −0.5 × 10 −3 eV 2 . The next step in complexity in an oscillation analysis with IceCube is to add the measurement of the neutrino energy, so the quantities L and E ν in equation (7) can be calculated separately. This is the approach followed in Ref. [75], where the energy of the neutrinos is obtained by using contained events in DeepCore and the assumption that the resulting muon is minimum ionising. Once the vertex of the neutrino interaction and the muon decay point have been identified, the energy of the muon   plane compared with results of Super-Kamiokande [70], T2K [71], MINOS [72] and NOvA [73]. A normal mass ordering is assumed. Figures from Ref. [74] can be calculated assuming constant energy loss, and it is proportional to the track length. The energy of the hadronic particle cascade at the vertex is obtained by maximising a likelihood function that takes into account the light distribution in adjacent DOMs. The neutrino energy is then the sum of the muon and cascade energies, E ν = E cascade + E µ . The most recent oscillation analysis from IceCube [74] improves on the mentioned techniques in several fronts. It is an all-sky analysis and also incorporates some degree of particle identification by reconstructing the events under two hypotheses: a ν µ charged-current interaction which includes a muon track, and a particle-shower only hypothesis at the interaction vertex. This latter hypothesis includes ν e and ν τ chargedcurrent interactions, although these two flavours can not be separately identified. The analysis achieves an energy resolution of about 25% (30%) at ∼20 GeV for muon-like (cascade-like) events and a median angular resolution of 10°(16°). Full sensitivity to lower neutrino energies, for example to reach the next oscillation minimum at ∼6 GeV, can only be achieved with a denser array, like the proposed PINGU low-energy extension [76].
In order to determine the oscillation parameters, the data is binned into a two-dimensional histogram where each bin contains the measured number of events in the corresponding range of reconstructed energy and arrival direction. The expected number of events per bin depend on the mixing angle, θ 23 , and the mass splitting, ∆m 2 32 , as shown in Fig. 5. This allows to determine the mixing angle θ 23 and the mass splitting ∆m 2 32 as the maximum of the binned likelihood. The fit also includes the likelihood of the track and cascade hypotheses. Systematic uncertainties and the effect of the Earth density profile are included as nuisance parameters. In this analysis, a full three-flavour oscillation scheme is used and the rest of the oscillation parameters are kept fixed to ∆m 2 21 = 7.53 × 10 −5 eV 2 , sin 2 θ 12 = 3.04 × 10 −1 , sin 2 θ 13 = 2.17 × 10 −2 and δ CP = 0. The effect of ν µ disappearance due to oscillations is clearly visible in the left panel of Fig. 7, which shows the number of events as a function of the reconstructed L/E ν , compared with the expected event distribution, shown as a dotted magenta histogram, if oscillations were not present. The results of the best fit to the data are shown in the right panel of Fig. 7. The best-fit values are ∆m 2 32 = 2.31 +0.11 −0.13 × 10 −3 eV 2 and sin 2 2θ 23 = 0.51 +0.07 −0.09 , assuming a normal mass ordering. The results of the two analyses mentioned above are compatible within statistics but, more importantly, they agree and are compatible in precision with those from dedicated oscillation experiments.

Flavour of Astrophysical Neutrinos
The neutrino oscillation phase in equation (7) depends on the ratio L/E ν of distance travelled, L, and neutrino energy, E ν . For astrophysical neutrinos we have to consider ultra-long oscillation baselines L corresponding to many oscillation periods between source and observer. The initial mixed state of neutrino flavours has to be averaged over ∆L, corresponding to the size of individual neutrino emission zones or the distribution of sources for diffuse emission. In addition, the observation of neutrinos can only decipher energies within an experimental energy resolution ∆E ν . The oscillation phase in (7) has therefore an absolute uncertainty that is typically much larger than π for astrophysical neutrinos. As a consequence, only the oscillation-averaged flavour ratios can be observed.
The flavour-averaged survival and transition probability of neutrino oscillations in vacuum, can be derived from Eq. (6) by replacing sin 2 ∆ ij → 1/2 and sin 2∆ ij → 0.  [77]. The best-fit flavour ratio is indicated by a white "×", with 68% and 96% confidence levels indicated by white lines. The expected oscillation-averaged composition is indicated for three different initial compositions, corresponding to standard pion decay (1 : 2 : 0), muon-damped pion decay (0 : 1 : 0), and neutron decay (1 : 0 : 0). The white "+" indicate the best-fit from a previous analysis [39]. (From Ref. [77]) The resulting expression can be expressed as To a good approximation, neutrinos are produced in astrophysical environments as a mixed state involving ν e , ν e , ν µ , and ν µ . Due to the similarity of neutrino and antineutrino signals in Cherenkov telescopes we consider in the following only flavour ratios of the sum of neutrino and anti-neutrino fluxes φ ν+ν with flavour ratios N e : N µ : N τ . Note, that the mixing angles shown in Table 1 are very close to the values for "tri-bi-maximal" mixing [78] corresponding to sin 2 θ 12 ∼ 1/3, sin 2 θ 23 ∼ 1/2 and sin 2 θ 13 ∼ 0. If we use this approximation then the oscillation-averaged spectrum will be close to a flavour ratio where x e = N e /N tot is the electron neutrino fraction on production. For instance, pion decays π + → µ + + ν µ followed by muon decay µ + → e + + ν e + ν µ produces an initial electron fraction of x e = 1/3. The resulting flavour ratio is then close to 1 : 1 : 1. It is also feasible that the muon from pion decay loses energy as a result of synchrotron radiation in strong magnetic fields ("muondamped" scenario) resulting in x e 0 and a flavour ratio of 4 : 7 : 7. Radioactive decay, on the other hand, will produce an initial electron neutrino fraction x e 1 and a flavour ratio 5 : 2 : 2. Figure 8 shows a visualisation of the observable neutrino flavour. Each location in the triangle corresponds to a unique flavour composition indicated by the three axis. The coloured markers correspond to the oscillationaveraged flavour ratios from the three scenarios (x e = 1/3, x e = 0, and x e = 1) discussed earlier, where the best-fit oscillation parameters have been used (instead of "tri-bi-maximal" mixing). The blue-shaded regions show the relative flavour log-likelihood ratio of a global analysis of IceCube data [77]. The best-fit is indicated as a white cross. IceCube's observations are consistent with the assumption of standard neutrino oscillations and the production of neutrino in pion decay (full or "muondamped"). Neutrino production by radioactive decay is disfavoured at the 2σ level.

Standard Model Interactions
The measurement of neutrino fluxes requires a precise knowledge of the neutrino interaction probability or, equivalently, the cross section with matter. At neutrino energies of less than a few GeV the cross section is dominated by elastic scattering, e.g., ν x + p → ν x + p, and quasi-elastic scattering, e.g., ν e + p → e + + n. In the energy range of 1-10 GeV, the neutrino-nucleon cross section is dominated by processes involving resonances, e.g. ν e + p → e − + ∆ ++ . At even higher energies neutrino scattering with matter proceeds predominantly via deep inelastic scattering (DIS) off nucleons, e.g., ν µ + p → µ − + X, where X indicates a secondary particle shower. The neutrino cross sections have been measured up to neutrino energies of a few hundreds of GeV. However, the neutrino energies involved in scattering of atmospheric and astrophysical neutrinos off nucleons far exceed this energy scale and we have to rely on theoretical predictions.
We will discuss in the following the expected cross section of high-energy neutrino-matter interactions. In weak interactions with matter the left-handed neutrino couples via Z 0 and W ± exchange with the constituents of a proton or neutron. Due to the scale-dependence of the strong coupling constant, the calculation of this process involves both perturbative and non-perturbative aspects due to hard and soft processes, respectively.

Deep Inelastic Scattering
The gauge coupling of quantum chromodynamics (QCD) increases as the renormalisation scale µ decreases, a behaviour which leads to the confinement of quarks and gluons at distances smaller that the characteristic size Λ −1 QCD (200MeV) −1 1 fm. In nature (except in high temperature environments (T Λ QCD ) as in the early universe) the only manifestations of coloured representations are composite gauge singlets such as mesons and baryons. These bound states consist of valence quarks, which determine the overall spin, isospin, and flavour of the hadron, and a sea of gluons and quark-anti-quark The kinematics of deep inelastic scattering.
pairs, which results from QCD radiation and pair-creation. These constituents of baryons and mesons are also called "partons". Due to the strength of the QCD coupling at small scales the neutrino-nucleon interactions cannot be described in a purely perturbative way. However, since the QCD interaction decreases as the renormalisation scale increases (asymptotic freedom) the constituents of a nucleon may be treated as loosely bound objects within sufficiently small distance and time scales (Λ −1 QCD ). Hence, in a hard scattering process of a neutrino involving a large momentum transfer to a nucleon the interactions between quarks and gluons may factorise from the sub-process (see Fig. 9). Due to the renormalisation scale dependence of the couplings this factorisation will also depend on the absolute momentum transfer Q 2 ≡ −q 2 . Figure 9 shows a sketch of a general lepton-nucleon scattering process. A nucleon N with mass M scatters off the lepton by a t-channel exchange of a boson. The final state consist of a lepton and a hadronic state H with centre of mass energy (P + q) 2 = W 2 . This scattering process probes the partons, the constituents of the nucleon with a characteristic size M −1 at length scales of the order of Q −1 . Typically, this probe will be "deep" and "inelastic", corresponding to Q M and W M, respectively. The sub-process between lepton and parton takes place on time scales which are short compared to those of QCD interactions and can be factorised from the soft QCD interactions. The intermediate coloured states, corresponding to the scattered parton and the remaining constituents of the nucleus, will then softly interact and hadronise into the final state H.
The kinematics of a lepton-nucleon scattering is conveniently described by the Lorentz scalars x = Q 2 /(2q · P), also called Bjorken-x, and inelasticity y = (q · P)/(k · P) (see Fig. 9 for definitions). In the kinematic region of deep inelastic scattering (DIS) where Q M and W M we also have Q 2 2q · p and thus x (q · p)/(q · P). The scalars x and y have simple interpretations in particular reference frames. In a reference frame where the nucleon is strongly boosted along the neutrino 3-momentum k the relative transverse momenta of the partons is negligible. The parton momentum p in the boosted frame is approximately aligned with P and the scalar x expresses the momentum fraction carried by the parton. In the rest frame of the nucleus the quantity y is the fractional energy loss of the lepton, y = (E − E )/E, where E and E are the lepton's energy before and after scattering, respectively.
From the previous discussion we obtain the following recipe for the calculation of the total (anti-)neutrinonucleon cross section σ(ν(ν)N). The differential leptonparton cross section may be calculated using a perturbative expansion in the weak coupling. The relative contribution of this partonic sub-process with Bjorken-x and momentum transfer Q 2 in the nucleon N is described by structure functions, which depend on the particular parton distribution functions (PDFs) of quarks ( f q (x, Q 2 )) and gluons ( f g (x, Q 2 )) . These functions must be measured in fixed target and accelerator experiments, that only access a limited kinematic region in x and Q 2 . Figure 10 shows the regions in the kinematical x-Q 2 -plane which have been covered in electron-proton (HERA), antiproton-proton (Tevatron), and proton-proton (LHC) collisions as well as in fixed target experiments with neutrino, electron, and muon beams (see, e.g., Ref. [79] and references therein).

Charged and Neutral Current Interactions
The parton level charged current interactions of neutrinos with nucleons are shown as the top two diagrams (a) and (b) of Fig. 11. The leading-order contribution is given by where G F 1.17 × 10 −5 GeV −2 is the Fermi coupling constant. The effective parton distribution functions are

Fig. 11
The parton level W (a/b) and Z (c/d) boson exchange between neutrinos and light quarks.
For antineutrino scattering we simply have to replace all f q by f q . These structure functions f q are determined experimentally via deep inelastic lepton-nucleon scattering or hard scattering processes involving nucleons. The corresponding relation of neutron structure function are given by the exchange u ↔ d and u ↔ d due to approximate isospin symmetry. In neutrino scattering with matter one usually makes the approximation of an equal mix between protons and neutrons. Hence, for an iso-scalar target, i.e., averaging over isospin, Analogously, the parton level neutral current (NC) interactions of the neutrino with nucleons are shown in the bottom two diagrams (c) and (d) of Fig. 11. The leading-order double differential neutral current cross section can be expressed as Here, the structure functions are given by The weak couplings after electro-weak symmetry breaking depend on the combination I 3 − q sin 2 θ W , where I 3 is the weak isospin, q the electric charge, and θ W the Weinberg angle. More explicitly, the couplings for left-handed (I 3 = ±1/2) and right-handed (I 3 = 0) quarks are given by As in the case of charged current interactions, the relation of neutron structure function f q are given by the exchange u ↔ d and u ↔ d and for an iso-scalar target one takes

High-Energy Neutrino-Matter Cross Sections
The expressions for the total charged and neutral current neutrino cross sections are derived from Eqs. (14) and (15) after integrating over Bjorken-x and momentum transfer Q 2 (or equivalently inelasticity y). The evolution of PDFs with respect to factorisation scale µ can be calculated by a perturbative QCD expansion and results in the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations [80,81,82,83]. The solution of the (leading-order) DGLAP equations correspond to a re-summation of powers (α s ln(Q 2 /µ 2 )) n which appear by QCD radiation in the initial state partons. However, these radiative processes will also generate powers (α s ln(1/x)) n and the applicability of the DGLAP formalism is limited to moderate values of Bjorken-x (small ln(1/x)) and large Q 2 (small α s ). If these logarithmic contributions from a small x become large, a formalism by Balitsky, Fakin, Kuraev, and Lipatov (BFKL) may be used to re-sum the α s ln(1/x) terms [84,85]. This approach applies for moderate values of Q 2 , since contributions of α s ln(Q 2 /µ 2 ) have to be kept under control.
There are unified forms [86] and other improvements of the linear DGLAP and BFKL evolution for the problematic region of small Bjorken-x and large Q 2 . The extrapolated solutions of the linear DGLAP and BFKL equations predict an unlimited rise of the gluon density at very small x. It is expected that, eventually, non-linear effects like gluon recombination g + g → g dominate the evolution and screen or even saturate the gluon density [87,88,89].
Note, that neutrino-nucleon scattering in charged (14) and neutral (15) current interactions via t-channel exchange of W and Z bosons, respectively, probe the parton content of the nucleus effectively up to momentum transfers of Q 2 M 2 Z/W (see Fig. 11). The present range of Bjorken-x probed by experiments only extends down to x 10 −4 at this Q-range, and it is limited to 10 −6 for arbitrary Q values. On the other hand, the Bjorken-x probed by neutrino interactions is, roughly, (20)  using HERAPDF1.5 at NLO are shown in Fig. 13. The general trend of the uncertainties can be understood by noting that as one moves to higher neutrino energy one also moves to lower x where the PDF uncertainties are increasing. The PDF uncertainties are smallest at 10 −2 < ∼ x < ∼ 10 −1 , corresponding to s ∼ 10 5 GeV 2 . Moving to smaller neutrino energies brings us into the high x region where PDF uncertainties increase again. This effect is greater for the HERAPDF1.5 because the HERA data have less statistics at high x than the fixed target data which are included in CT10; however these data have further uncertainties that are not fully accounted for in CT10, e.g. heavy target corrections, deuterium corrections and assumptions regarding higher twist effects. When the high x region becomes important the neutrino and antineutrino cross-sections are different because the valence contribution to xF 3 is now significant. This is seen in Fig. 13, as is the onset of the linear dependence of the cross-sections for s < M 2 W . Note that our predictions are made for Q 2 > 1 GeV 2 since perturbative QCD cannot sensibly be used at lower values. Moreover for s below ∼ 100 GeV 2 , there can be contributions to the cross-section of O(10%) from even lower values of Q 2 which are not accounted for here; hence we do not show results for E ν below 50 GeV where there are other contributions to the neutrino cross-section and the use of a code such as GENIE [59] is appropriate. For higher energies, we intend to upgrade ANIS [57] to use the HERAPDF1.5 (differential) cross-sections. Meanwhile we have provided the total DIS cross-sections for CC and NC scattering of neutrinos and antineutrinos on isoscalar targets in Tables I and II and recommend these as a benchmark for use by experimentalists. These cross-sections as well as those for isoscalar targets are available from a webpage [60]; differential cross sections are available upon request. Any measured deviation from these values would signal the need for new physics beyond the DGLAP formalism.

VII. ACKNOWLEDGEMENTS
We are grateful to James Ferrando for providing us with an up-to-date version of DISPRED and for discussions. We also thank Mike Whalley and Voica Radescu for the speedy implementation of HERAPDF1.5 in LHAPDF. PM and SS thank their colleagues in the Auger and IceCube collaborations for stimulating exchanges and acknowledge partial  13. The general trend of the uncertainties can be understood by noting one also moves to lower x where the PDF uncertainties are increasing.
x < ∼ 10 −1 , corresponding to s ∼ 10 5 GeV 2 . Moving to smaller neutrino ere PDF uncertainties increase again. This effect is greater for the less statistics at high x than the fixed target data which are included certainties that are not fully accounted for in CT10, e.g. heavy target ptions regarding higher twist effects. When the high x region becomes s-sections are different because the valence contribution to xF 3 is now nset of the linear dependence of the cross-sections for s < M 2 W . Note eV 2 since perturbative QCD cannot sensibly be used at lower values. be contributions to the cross-section of O(10%) from even lower values e we do not show results for E ν below 50 GeV where there are other d the use of a code such as GENIE [59] is appropriate. For higher se the HERAPDF1.5 (differential) cross-sections. Meanwhile we have and NC scattering of neutrinos and antineutrinos on isoscalar targets benchmark for use by experimentalists. These cross-sections as well a webpage [60]; differential cross sections are available upon request. uld signal the need for new physics beyond the DGLAP formalism. This shows that high-energy neutrino-nucleon interactions beyond 100 PeV strongly rely on extrapolations of the structure functions. Figure 12 shows the results of Ref. [90,91] which used an update of the PDF fit formalism of the published ZEUS-S global PDF analysis [92]. The total cross sections in the energy range 10 7 ≤ (E ν /GeV) ≤ 10 12 can be approximated to within ∼ 10% by the relations [91], Neutrino electron interactions can often be neglected with respect to neutrino nucleon interactions due to the electron's small mass. There is, however, one exception with ν e + e − interactions, because of the intermediateboson resonance formed in the neighbourhood of E res ν = M 2 W /2m e 6.3 PeV, generally referred to as the Glashow resonance [10]. The total cross section for the resonant where B in = Br(W − → ν e + e − ) and B out = Br(W − → X) are the corresponding branching ratios of W decay and Γ W 2.1 GeV the W decay width. The branching ratios into ν α + α are 10.6% and into hadronic states 67.4% [41].

Neutrino Cross Section Measurement with IceCube
Similar to the study of neutrino oscillations, that can be inferred from the low-energy atmospheric neutrino flux that reaches IceCube from different directions, highenergy atmospheric neutrinos can be used to measure the neutrino-nucleon cross section at energies beyond what is currently reached at accelerators. The technique is based on measuring the amount of atmospheric muonneutrinos as a function of zenith angle θ zen , and compare it with the expected number from the known atmospheric flux assuming the Standard Model neutrino cross sections. Neglecting regeneration effects, the number of events scales as where X(θ zen ) is the integrated column depth along the line of sight (n(θ zen )) from the location of IceCube (r IC ), The neutrino-matter cross section σ νN increases with neutrino energy, and above 100 TeV the Earth becomes opaque to vertically up-going neutrinos, i.e., neutrinos that traverse the whole Earth, see Fig. 13. Therefore, any deviation from the expected absorption pattern of atmospheric neutrinos can be linked to deviations from the assumed cross section, given all other inputs are known with sufficient precision.
IceCube has performed such an analysis [93] by a maximum likelihood fit of the neutrino-matter cross section. The data, binned into neutrino energy, E ν , and neutrino arrival direction, cos θ zen , was compared to the expected event distribution from atmospheric and astrophysical 15 neutrinos. Deviations from the Standard Model cross section σ SM were fitted by the ratio R = σ νN /σ SM . The analysis assumes priors on the atmospheric and astrophysical neutrino flux based on the baseline models in Refs. [94,30,24]. In practice, the likelihood maximisation uses the product of the flux and the cross section, keeping the observed number of events as a fixed quantity. Thus, trials with higher cross sections must assume lower fluxes (or vice-versa) in order to preserve the total number of events. The procedure is thus sensitive to neutrino absorption in the Earth alone, and not to the total number of observed events. Since the astrophysical flux is still not known to a high precision, the uncertainties in the normalisation and spectral index were included as nuisance parameters in the analysis. Other systematics considered are the Earth density and core radius as obtained from the Preliminary Earth Model [95], the effects of temperature variations in the atmosphere, which impact the neutrino flux during the year, and detector systematics.
The analysis results in a value of . This is compatible with the Standard Model prediction (R = 1) within uncertainties but, most importantly, it is the first measurement of the neutrino-nucleon cross section at an energy range (few TeV to about 1 PeV) unexplored so far with accelerator experiments [41]. This is illustrated in Fig. 14 which shows current accelerator measurements (within the yellow shaded area) and the results of the IceCube analysis as the light brown shaded area. The authors of Ref. [96] performed a similar analysis based on six years of high-energy starting event data. Their results are also consistent with perturbative QCD predictions of the neutrino-matter cross section.

Probe of Cosmic Ray Interactions with IceCube
On a slightly different topic, but still related to the products of cosmic ray interactions in the atmosphere, the high rate of atmospheric muons detected by IceCube can be used to perform studies of hadronic interactions at high energies and high momentum transfers. Muons are created from the decays of pions, kaons and other heavy hadrons. For primary energies above about 1 TeV, muons with a high transverse momentum, p t 2 GeV, can be produced alongside the many particles created in the forward direction, the "core" of the shower. This will show up in IceCube as two tracks separated by a few hundred meters: one track for the main muon bundle following the core direction, and another track for the high-p t muon. The muon lateral distribution in cosmic-ray interactions depends on the composition of the primary flux and details of the hadronic interactions [97,98]. If the former is sufficiently well known, the measurement of high-p t muons can be used to probe hadronic processes involving nuclei and to calibrate existing Monte-Carlo codes at energies not accessible with particle accelerators.
The lateral separation, d t , of high p t muons from the core of the shower is given by d t = p t H/E µ cos θ zen , where H is the interaction height of the primary with a zenith angle θ zen . The initial muon energy E µ is close to that at ground level due to minimal energy losses in the atmosphere. That is, turning the argument around, the identification of single, laterally separated muons at a given d t accompanying a muon bundle in IceCube is a measurement of the transverse momentum of the muon's parent particle, and a handle into the physics of the primary interaction. Given the depth of IceCube, only muons with an energy above ∼400 GeV at the surface can reach the depths of the detector. This, along with the inter-string separation of 125 m, sets the level for the minimum p t accessible in IceCube. However, since the exact interaction height of the primary is unknown and varies with energy, a universal p t threshold can not be given. For example, a 1 TeV muon produced at 50 km height and detected at 125 m from the shower core has a transverse momentum p t of 2.5 GeV.
Our current understanding of lateral muon production in hadronic interactions shows an exponential behaviour at low p t , exp(−p t /T), typically below 2 GeV, due to soft, non-perturbative interactions, and a powerlaw behaviour at high p t values, (1 + p t /p 0 ) −n , reflecting the onset of hard processes described by perturbative QCD. The approach traces back to the QCD inspired "modified Hagedorn function" [99,100]. The parameters T, p 0 and n can be obtained from fits to proton-proton or heavy ion collision data [100,101]. This is also the behaviour seen by IceCube. Fig. 15 shows the muon lateral distribution at high momenta obtained from a selection of events reconstructed with a two-track hypothesis in the 59-string detector [102], along with a fit to a compound exponential plus power-law function. Due to the size of the 59-string detector and the short live time of the analysis (1 year of data), the statistics for large separations is low and fluctuations in the data appear for track separations beyond 300 m. Still, the presence of an expected hard component at large lateral distances (high p t ) that can be described by perturbative quantum chromodynamics (a power-law behaviour) is clearly visible.
Significant discrepancies also exist between interaction models on the expected p t distribution of muons [103]. The high-p t muon yield per collision depends, both, on the primary composition (protons typically producing higher rates of high-p t muons than iron interactions) and the hadronic interaction model. Future IceCube analyses with the larger, completed, 86-string detector using Ice-Cube/IceTop coincident events can extend the range of the cos θ zen distribution as well as provide an updated comparison of the muon p t distribution with predictions from existing hadronic models with higher statistics than that shown in Fig. 15. There is definitely complementary information from neutrino telescopes to be added to the efforts in understanding hadronic interactions using air shower arrays and heavy-ion and hadronic accelerator experiments [104,105].

Non-Standard Neutrino Oscillations and Interactions
In the previous two sections we have summarised the phenomenology of weak neutrino interactions and standard oscillations based on the mixing between three active neutrino flavour states and the eigenstates of the Hamiltonian (including matter effects). However, the Standard Model of particle physics does not account for neutrino masses and is therefore incomplete. The neces-sary extensions of the Standard Model that allow for the introduction of neutrino mass terms can also introduce non-standard oscillation effects that are suppressed in a low-energy effective theory. This is one motivation to study non-standard neutrino oscillations. In the following, we will discuss various extensions to the Standard Model that can introduce new neutrino oscillation effects and neutrino interactions. The large energies and very long baselines associated with atmospheric and cosmic neutrinos, respectively, provide a sensitive probe for these effects.
For the following discussion it is convenient to introduce neutrino oscillations via the evolution of the density operator ρ of a mixed state. For a given Hamiltonian H, the time evolution of the neutrino density operator is governed by the Liouville equation, The solution to the Liouville equation allows to describe neutrino oscillation effects in a basis-independent way. For instance, a neutrino that is produced at times t = 0 in the flavour state ν α can be represented by ρ(0) = Π α , where Π α = |ν α ν α | is the projection operator onto the flavour state |ν α . The transition probability between two flavour states ν α and ν β can then be recovered as the expectation value of the projector Π β , which is given by trace For standard neutrino oscillations the Hamiltonian is composed of two terms, H SM = H 0 + V mat , describing the free evolution in vacuum and coherent matter effects, respectively. The free Hamiltonian in vacuum can be written as where the sum runs over the available flavours and we have introduced the projection operators Π i ≡ |ν i ν i | and Π i ≡ |ν i ν i | onto neutrino and anti-neutrino mass eigenstates. Matter effects introduced earlier can be cast into the form where Π e ≡ |ν e ν e | and Π e ≡ |ν e ν e | are the projections onto electron neutrino and anti-neutrinos states, respectively, and N e is the number density of electrons in the background.
In the following, we will discuss non-standard contributions to the Liouville equation that we assume to take on the forṁ The second term on the right hand side of Eq. (30) accounts for additional terms in the Hamiltonian, that affect the phenomenology of neutrino oscillations. These 17 terms are therefore parametrised by a sequence of Hermitian operators, H n = H † n , that preserve the unitarity of the density operator. The third term corresponds to non-unitary effects, related to dissipation or decoherence, that can not only affect the neutrino flavour composition, but also their spectra. In particular, we will focus in the following on cases where the sequence of terms D n can be expressed in terms of a set of Lindblad operators L j as These Lindblad operators act on the Hilbert space of the open quantum system, H, and satisfy ∑ j L † j L j ∈ B(H), where B(H) indicates the space of bounded operators acting on H [106]. We will discuss in the following various instances of these Lindblad operators for neutrino decoherence, neutrino decay, or hidden neutrino interactions with cosmic backgrounds.

Effective Hamiltonians
Non-trivial mixing terms of the effective Hamiltonian can be generated in various ways, including non-standard interactions with matter and Standard Model extensions that violate the equivalence principle, Lorentz invariance, or CPT symmetry (see Ref. [107] for references). We will first study oscillations that are induced by additional terms to the Hamiltonian that can be parametrised in the form [107] Here we introduced the projectors Π a = |ν a ν a | and Π a = |ν a ν a | onto a new set of states |ν a that are related to the usual flavor states by a new unitary mixing matrix The expansion parameters δ a have mass dimensions 1 − n with integer n.
It is convenient to group these contributions into CPTeven terms obeying the relation δ a = δ a and CPT-odd terms with δ a = −δ a . For a CPT-symmetric process we have P(ν α → ν β ) = P(ν β → ν α ). In particular, the survival probability between neutrinos and anti-neutrinos is the same. CPT-odd terms break this symmetry. Indeed, we have already encountered such a CPT-odd term as the CP-violating matter effect contributing by the effective potential of Eq. (29). The corresponding expansion in terms the effective Hamiltonian of Eq. (32) is given by U = 1, δ e = −δ e = √ 2G F N e , δ µ,τ = 0, and n = 0. On the other hand, the free Hamiltonian in Eq. (28) is a CPT-even term with U = U PMNS , δ i = δ i = m 2 i /2, and n = −1. The contribution of effective Hamiltonians can be tested by the oscillation of atmospheric muon neutrinos [107]. For simplicity, we will assume that this can be treated effectively as a two-level system, analogous to the case of standard atmospheric neutrino oscillations. The most general form of the unitary matrix U is then The ν µ survival probability can then be expressed as [107] where we have introduced the scaling parameter and with ∆δ = δ 3 − δ 2 . The effective rotation angle θ eff is given by sin 2 2θ eff = 1 R 2 sin 2 2θ atm + R 2 sin 2 2ξ + 2R sin 2θ atm sin 2ξ cos η . (38) Note that the survival probability in Eq. (35) reduces to the familiar expression of atmospheric oscillation in the limit R → 0.  Table 1). The coloured contours show the range of oscillation-averaged flavor ratios based on the unitarity of the neutrino mixing matrix. Figure 16 shows the results of an analysis of atmospheric muon neutrino data in the range 100 GeV to 10 TeV taken by AMANDA-II in the years 2000 to 2006 [109]. The data was binned into two-dimensional histograms in terms of the number of hit optical modules (as a measure of energy) and the zenith angle (cos θ zen ). The predicted effect of non-standard oscillation parameters can be compared to the data via a profile likelihood method. No evidence of non-standard neutrino oscillations was found and the statistically allowed region of the [∆δ, sin 2 2ξ]plane is shown as 90%, 95%, 99% C.L. These results are derived under the assumption that η = π/2 in the unitary mixing matrix. The red dashed line shows the 90% limit of a combined analysis by Super-Kamiokande and K2K [108] which is compatible with the AMANDA-II bound. The projected IceCube 90% sensitivity after ten years of data taking is given as a yellow dotted line and may improve the limit on ∆δ by one order of magnitude [107].
Effective Hamiltonians with positive energy dependence n > 0 can dominate the standard Hamiltonian H 0 at sufficiently high energy. For astrophysical neutrinos this can lead to the situation that the oscillation-averaged flavour transition probability is completely dominated by the new unitary transition matrix leading to This effect can therefore be tested by looking for statistically significant deviations from the predictions of standard oscillations in Eq. (12). However, the transition probabilities described by Eq.  a non-trivial set of unitarity conditions: The resulting oscillation-averaged flavour composition from these unitarity bounds are shown as the contours in Fig. 17. We assume three initial flavour composition from pion decay (1:2:0), muon-damped pion decay (0:1:0), and neutron decay (1:0:0).

Violation of Lorentz Invariance
One of the foundations of the Standard Model of particle physics is the principle of Lorentz symmetry: The fundamental laws in nature are thought to be independent of the observer's inertial frame. However, some extensions of the Standard Model, like string theory or quantum gravity, allow for the spontaneous breaking of Lorentz symmetry, that can lead to Lorentz-invariance violating (LIV) effects in the low-energy effective theory. There also exist a deep connection between the appearance of LIV effects with the violation of CPT-invariance 3 in local quantum field theories [111]. Such effects were incorporated in the Standard Model Extension (SME), an effective-field Lorentz-violating extension of the Standard Model, which includes CPT-even and CPT-odd terms [112]. The SME provides a benchmark for experiments to gauge possible Lorentz violating processes in nature, by expressing experimental results in terms of the parameters of the model. The size of LIV effect is expected to be suppressed by Planck scale M P 10 19 GeV (or Planck length λ P 10 −33 cm), consistent with the strong experimental limits on the effect [113,114].
Oscillations of atmospheric neutrinos with energies above 100 GeV provide a sensitive probe of LIV effects. For instance, LIV in the neutrino sector can lead to small differences in the maximal attainable velocity of neutrino states [116]. Since the "velocity eigenstates" are different from the flavour eigenstates, a new oscillation pattern  [115] can arise. The effect can be described in the framework of effective Hamiltonians described in the previous section by CPT-even states with n = 1. For the approximate two-level system with survival probability described by Eq. (35) we can identify effective Hamiltonian parameters as the velocity difference ∆δ = ∆c, together with a new mixing angle ξ and a phase η. Higher order contributions n > 1 have been considered for non-renormalisable LIV effects caused by quantum mechanical fluctuations of the space-time metric and topology [117]. Both the ∆δ ∝ E (n = 1) and the ∆δ ∝ E 3 (n = 3) cases have been examined in the context of violations of the equivalence principle (VEP) [118,119,120]. The 90% C.L. upper limits on the corresponding coefficients ∆δ/E n ν [GeV 1−n ] are shown in Tab. 2.
Violations of Lorentz invariance can also manifest themselves by the dependence of neutrino oscillations on a preferred orientation of the neutrino arrival direction. This effect can also be described by effective Hamiltonians similar to Eq. (32) where the expansion parameters not only depend on energy but also on direction. The contribution of these terms were studied in the analysis [121] using data collected by the 40-string configuration of Ice-Cube between April 2008 and May 2009. Due to IceCube's unique position at the geographic South Pole, the field of view of the observatory is constant in time. The atmospheric neutrino data extracted from this period was binned in terms of sidereal time. If neutrino oscillations depend on the relative neutrino momenta in the cosmic rest frame, then the muon neutrino data in this reference frame is expected to show oscillation patterns. The nonobservation of this effect allowed to constrain these LIV parameters.
A more recent search for isotropic LIV effects with the complete IceCube detector was the subject of the analysis presented in [115], based on data collected during the period from May 2010 to May 2012. The LIV effects were parametrised by two parameters ρ d and θ d , which are related to the strength of the LIV effect and the a combination of other LIV parameters of the SME [112], respectively. These parameters accompany the expansion of the Hamiltonian in powers of energy, and the subscript d refers to the power of the corresponding operator in the Hamiltonian. Most experiments are sensitive to effects of dimension d = 3 and d = 4, but the energy reach of IceCube makes it possible to extend the search to LIV effects induced by operators of dimension up to d = 8. The analysis used binned atmospheric neutrino data into a horizontal (cos θ zen > −0.6) and a vertical (cos θ zen < −0.6) sample which allowed to study the LIV effect by the energy-dependent ratio of the samples using a likelihood analysis. No LIV effects were identified, the best fit values of all the ρ d parameters being compatible with zero. This allowed to place limits in the [ρ d -cos θ d ] parameter space. The results are shown in Fig. 18 for The violation of Lorentz invariance associated with Planck-scale physics can also affect neutrino spectra over long baselines. The LIV effects can result in modified dispersion relations, e.g., E 2 − p 2 = m 2 − E 2 , that introduce non-trivial maximal particle velocities [116]. Whereas at low energies the Lorentz invariance is recovered, E 2 − p 2 m 2 , at high energies we can observe sub-or superluminal maximal particle velocities. These can allow otherwise forbidden neutrino decays, in particular, vacuum pair production, ν α → ν α + e + + e − , and vacuum

−2∆LLH
90% CL 90% CI 90% CL Super-K 90% CI Salvado et al Fig. 19 Confidence limits on the NSI parameter µτ using the event selection from [75,125] shown as solid vertical red lines. Similarly, dashed vertical red lines show the 90% credibility interval using a flat prior on µτ and where we have profiled over the nuisance parameters. The light blue vertical lines show the Super-Kamiokande 90% confidence limit [126]. The light green lines show the 90% credibility region from [127]. Finally, the horizontal dash-dot line indicates the value of −2∆LLH that corresponds to a 90% confidence interval according to Wilks' theorem. Figure from [128].

Non-Standard Matter Interactions
We have already discussed coherent scattering of neutrinos and anti-neutrinos in dense matter, that can be accounted for by an effective matter potential in the standard Hamiltonian. Since neutrino oscillations are only sensitive to non-universal matter effects, only the unique charged-current interactions of electron neutrinos and anti-neutrinos with electrons are expected to contribute. However, non-standard interactions (NSI) can change this picture.
A convenient theoretical framework to study NSI contributions is given by the Hamiltonian, The dimensionless coefficients αβ measure the NSI contribution relative to the strength of the standard matter potential. Hermiticity of the Hamiltonian requires that αβ = * βα . Non-standard interactions with nucleons are not necessarily universal for electrons and valence quarks in matter [129]. Therefore, the parameters αβ can be considered as effective parameters of the form where are the couplings to up-and down-type valence quarks in matter. In charge-neutral matter we must have N p = N e . The average neutron abundance according to Earth density models is N n 1.051 [95]. Therefore, we have In the full three-flavour neutrino system, the hermiticity condition leaves nine degrees of freedom for the αβ coefficients, that reduce to eight after absorption of a global phase. The case of atmospheric muon neutrino oscillations can again be approximated by oscillations in the ν µ -ν τsystem. The oscillations induced by NSI interactions has three degrees of freedom that can be parametrized by ττ − µµ = and µτ = exp(iα). We can make use of the fact that the NSI Hamiltonian can be recast as an effective Hamiltonian as in Eq. (32). The NSI Hamiltonian in the two level system can be diagonalised by η = α and tan 2ξ = 2 / . After removing a global phase, the effective splitting is ∆δ = √ 2G F N e √ 1 + 4 2 / 2 with σ = −1.
The effect of NSI in atmospheric neutrino oscillations with IceCube data were studied in Refs. [128] and [130]. The IceCube analysis of Ref. [128] is based on three years of data collected by the low-energy extension DeepCore, that was also used in the standard neutrino oscillation analysis discussed in Section 4. The data was binned into a two-dimensional histogram in terms of reconstructed neutrino energy and zenith angle, cos θ zen , and analysed via a profile likelihood method. The analysis focused on NSI interactions with d-quarks assuming ττ . The constraints for µτ are shown in Fig. 19.
A two-dimensional analysis on the parameters µτ and was carried out by the authors of Ref. [130]. The analysis was performed on the atmospheric neutrino measurement done with one year of IceCube-40 data [23] and with two years of public DeepCore data, which was initially used for the first oscillation analysis by IceCube [69]. A combined confidence level region of allowed values for µτ and is obtained from this analysis, see Fig. 20. Although the analysis is only based on the arrival direction distribution of muon-neutrinos and not their energy, it is competitive with previous results from Super-Kamiokande, and even produces more restrictive limits.

Neutrino Decoherence
The Hamiltonian evolution in Eq. (26) is a characteristic of physical systems isolated from their surroundings. The time evolution of such a quantum system is given by the continuous group of unitarity transformations, U t = e −iHt , where t is the time. The hermiticity of the Hamiltonian guarantees the reversibility of the processes, For open quantum systems, the introduction of dissipative effects lead to modifications of Eq. (26) that account for the irreversible nature of the evolution. The transformations responsible for the time evolution of these systems are defined by the operators of the Lindblad quantum dynamical semi-groups [106]. Since this does not admit an inverse, such a family of transformations has the property of acting only forward in time. The monotonic increase of the von Neumann entropy, S(ρ) = −Tr (ρ ln ρ), implies the hermiticity of the Lindblad operators, L j = L † j [131]. In addition, the conservation of the average value of the energy can be enforced by taking [H, L j ] = 0 [132].
In the three-flavour basis it is convenient to expand the density operator ρ, the Hamiltonian H and the Lindblad operators L i in terms of Hermitian 3 × 3 matrices F µ (µ = 0, . . . , 8), given by, where λ i (i = 1, . . . , 8) are the Gell-Mann matrices [133]. This set of basis matrices satisfy the orthogonality condition Tr(F † µ F ν ) = δ µν /2. The Liouville equation (30) can now be expressed in terms of the expansion coefficients ρ i and h i of the density operator and Hamiltonian, respectively, as the matrix equatioṅ The coefficients f ijk are the SU(3) structure constants defined by [F i , F j ] = i ∑ k f ijk F k [133]. The decoherence effects are cast into a matrix where l (n) i are the expansion coefficient of the Lindblad operator L n .
In atmospheric neutrino data we are interested in the survival probability of muon neutrinos (and antineutrinos). Typically, it is assumed that the 8 × 8 matrix D ij takes on a diagonal form, i.e., D ij δ ij D i , with positive diagonal elements. The three-level system including decoherence and oscillations can be further simplified assuming universal decoherence parameters D i = D. The muon neutrino survival probability can then be expressed as [133], Note that this reduces to the standard survival probability in the case D → 0 and leads to a flavour ratio 1 : 1 : 1 in the case of DL → ∞, independent of initial condition. Table 2 shows the limits on the universal decoherence parameter D/E n ν [GeV n−1 ] for various energy dependencies (n = 1, 2, 3) derived from AMANDA-II data of the years 2000 to 2006 [109].
On the other hand, keeping decoherence parameters D i independent, we can derive the oscillation-averaged flavour survival/appearance probability in the form [133] P avg ν α →ν β = Only the parameters D 3 and D 8 contribute in this case. Again, for D 3 L → ∞ and D 8 L → ∞ we recover a flavour ratio 1 : 1 : 1, as expected for full decoherence.

Neutrino Decay
The flavour ratio of neutrinos from astrophysical sources can be significantly altered if neutrinos can decay en route to Earth [134,135,136,137,138]. Two-body decays of the form ν i → ν j + X with a massless light scalar X (e.g., a Majoron [139,140,141]) are only weakly limited by solar neutrino data [142] with τ i 10 −4 s(m i /eV). In the following, we will assume that the decay of neutrino mass eigenstates can be written as a dissipation term in Eq. (30) of the form The first term describes the disappearance of a neutrino mass eigenstate |ν i with total decay width Γ i (E i ). The second term describes the appearance of neutrino mass 22 eigenstates from the transition |ν i and |ν j with differential production rate γ i→j (E i , E j ). Note, that the corresponding set of Lindblad operators representing the expression in Eq. (50) can be written as where the index i runs over active neutrino states and the index j over active and sterile states. In contrast to quantum decoherence discussed in the previous section, this set of Lindblad operators are non-Hermitian and do not commute with the Hamiltonian. As a further approximation, we will consider neutrino propagation over cosmological distances, i.e., all oscillation terms from neutrino mass differences are averaged. As a representative example, we consider transitions that leave the neutrino energy constant, i.e., γ i→j Br i→j Γ i δ(E i − E j ), and normal ordering of neutrino masses with transitions 3 → 2, 3 → 1, and 2 → 1. In the mass eigenbasis the evolution of the density matrix has then the forṁ In the oscillation-averaged solution only the diagonal elements ρ ii have non-zero contributions. We can then express the neutrino survival/appearance probability as where the elements of transition matrix, M ji , can be written in the form [143] The sum in the previous equation runs over all production chains c 1 → . . . → c n c with c 1 = i and c n c = j and In this particular scenario, the final neutrino flavour ratio for LΓ i 1 will thus correspond to the composition of the lowest mass eigenstate.
Note, that the decay rates Γ i are expected to decrease with energy due to relativistic boosting of the mass eigenstate's lifetime. Therefore, the neutrino flavour composition can experience strong energy dependencies. On the other hand, active neutrino decay into sterile neutrinos can introduce spectral cutoffs due to the energy dependence of the neutrino lifetime. This process is limited by the observation of IceCube's TeV-PeV neutrino flux and could be responsible for a tentative cutoff [136]. Neutrino decay has also been considered as a possibility to alleviate a mild tension in the best-fit power-law spectra between cascade-and track-dominated IceCube data [138].
Astrophysical neutrinos propagating over cosmic distances are also susceptible to feeble interactions with cosmic backgrounds. In particular, feeble interactions with the cosmic neutrino background (CνB) that can be enhanced by resonant interactions, e.g., ν α + ν α → Z → ν β + ν β have been discussed as a source for absorption features [144,145,146,147]. The evolution of the neutrino density matrix is identical to that for neutrino decay with dissipation term as in Eq. (50). However, in this case the interaction rates have a non-trivial dependence on redshift via the density evolution of the CνB.

Sterile Neutrinos
It is conceivable that the three neutrino flavour states, ν e , ν µ , and ν τ , are augmented by additional states. The number of light "active" neutrinos, i.e., neutrinos that take part in the Standard Model weak interactions is limited to the known flavour states from the observed decay width of the Z-boson. However, it is feasible that there are light "sterile" neutrino states with no Standard Model interactions. These sterile states would imply the existence of additional neutrino mass eigenstates that could impact the three-flavour neutrino oscillation phenomenology, due to an extended PMNS mixing matrix U, The term "sterile" referring to neutrinos was introduced by Pontecorvo already in 1968 when discussing the, at that time hypothetical, possibility of vacuum neutrino oscillations [44]. Many extensions of the Standard Model that relate to the appearance of neutrino mass terms predict the existence of sterile neutrinos. As discussed earlier, the righthanded neutrino field, that can provide a Dirac mass term m D ν L ν R + h.c., does not interact via weak interactions and is therefore sterile. However, in the minimal type I seesaw models (see, e.g., Ref. [63]) these sterile neutrinos have a large Majorana mass term, Mν R Cν T R /2 + h.c., with m D M, that give rise to a large effective neutrino mass after diagonalising the mass matrix. These massive sterile states are practically unobservable in low-energy oscillation experiments. On the other hand, for values of M m D ("pseudo-Dirac" case), the active and sterile state mass states are degenerate after diagonalisation, leading to maximal mixing between the left (active) and right (sterile) states.
The minimal sterile neutrino model is a "3+1" model where, in addition to the three standard weakly-interacting neutrino flavours, one additional heavier sterile neutrino state is added. Such a simple extension of the neutrino sector has been advocated to explain certain tensions between experimental results from accelerator [155,156,157], reactor [158] and radio-chemical [159] experiments, and the predictions from the standard three active  [148,149,150,151], and the 99% C.L. allowed region from global fits to appearance experiments including MiniBooNE and LSND, assuming |U e4 | 2 =0.0023 [152] and |U e4 | 2 =0.0027 [153], respectively. flavour scenario. In the most general case, the introduction of one sterile neutrino adds six new parameters to the neutrino oscillation phenomenology [152]: three mixing angles θ 14 , θ 24 , θ 34 , two CP-violating phases, δ 14 and δ 34 , and one mass difference, ∆m 2 41 , where the indexes '1-3' stand for the known neutrino mass states and '4' for the sterile state.
Although sterile neutrinos can not be detected directly, their existence can leave an imprint on the oscillation pattern of active neutrinos. The sterile neutrino modifies the oscillation pattern of the standard neutrinos since these can now undergo vacuum oscillations into the new state, with a probability that is proportional to the new mixing angles. The period of these oscillations can be small, smaller than the directional resolution of IceCube, and the net effect is then to distort the overall ν µ flux normalisation with respect to the three-flavour case. An additional effect arises from the different interactions of flavours with matter when traversing the Earth [160,161]. The new possibility to oscillate to a state that does not interact results in energy and angular dependent oscillation amplitudes that depend on the mixing angles, but also on the new ∆m 2 41 . More precisely, the comparison between the oscillation pattern in the (energy, zenith) phase space predicted by the sterile "3+1" model and the pattern seen in experimental data is what IceCube exploits to set limits on the sterile mixing parameters. As for the standard oscillation case described in section 4, searches for sterile neutrinos in IceCube are based on the detection of the disappearance of muon neutrinos and are more sensitive to θ 24 . Therefore the choice of a simplified minimal mixing scenario where only θ 24 is not zero is justified, and is indeed the approach followed in the analyses described below. Nonzero values of the weakly constrained θ 34 within current limits would not significantly affect the results presented here [162].
There have been several searches in the past at sub-TeV energies with neutrinos from the Sun, reactors, and accelerator setups (see, e.g., Refs. [152,163,164,165] and references therein). The large flux of atmospheric neutrinos that reach IceCube and the wide range of energy response of the detector allow to search for signatures of anomalous ν µ + ν µ disappearance caused by oscillations to an sterile neutrino, ν s , at TeV energies, an energy not probed before. Furthermore, the DeepCore detector can extend the search down to about 6 GeV, improving previous limits from smaller detectors.
The analysis techniques are very similar to the search for standard oscillations and are based on comparing the equivalent oscillogram from Fig. 5 for oscillations with an sterile component, to the measured data. Although the low-energy tracks provide a quite short lever arm to reconstruct their direction with precision, and the angular resolution of the analysis varies between 6°and 12°, depending on energy. This is enough, though, to perform the analysis (see Fig. 5 in Ref. [125]) since the differences in the oscillation pattern in the presence of a sterile neutrino can still be distinguished from the non-sterile case with such angular resolution. At higher energies the analyses necessarily follow a slightly different approach since the tracks originate outside the detector and the energy deposited in the hadronic shower at the interaction vertex is not accessible. Such analyses use measured muon energies instead of neutrino energies, like the one presented in Ref. [154]. The muon energy is reconstructed from the light emission profile from stochastic energy losses of the muon along its trajectory [17], achieving an energy resolution of σ log 10 (E µ /GeV) ∼ 0.5. Since at the energies of this analysis the muon track can be well reconstructed, the angular resolution reaches values between 0.2°and 0.8°, depending on incoming angle.
As with searches for standard oscillations, the aim is to reach a pure sample of atmospheric neutrinos to be able to compare the measured number of events as a function of angle and energy with the model prediction. None of the analyses detected a deviation from the expected standard oscillation scenario and, therefore, limits on the oscillation parameters which depend on the additional neutrino can be set. The low-energy analysis uses a 8×8 binned grid in the (energy, cos θ zen ) plane while, due to its higher energy resolution, the high-energy analysis uses a 10×20 binned grid. A log-likelihood approach is used to find the best fit to the data given the model parameters. Confidence levels are calculated from the difference between the profile log-likelihood and the loglikelihood at the best fit point. The results are illustrated in Figs. 21 and 22. Figure 21 shows the result of the highenergy analysis expressed in the (∆m 2 14 , sin 2 2θ 24 ) plane. The figure shows the 90% confidence level contour (red line) compared with 90% exclusions from previous disappearance searches. The exclusion is compatible with the sensitivity (green and yellow areas) calculated under the assumption of the no-sterile neutrino hypothesis. The result therefore disfavours much of the parameter space of the "3+1" model. On the other hand, Fig. 22 shows the results of an analysis using low-energy events (6 GeV to 60 GeV) contained in DeepCore. The results are shown in the (|U τ4 | 2 , |U µ4 | 2 ) plane, where |U τ4 | is defined as sin 2 θ 34 cos 2 θ 34 and |U µ4 | is just sin 2 θ 24 . The best-fit value of these elements of the mixing matrix are |U µ4 | 2 = 0.0 and |U τ4 | 2 = 0.08, while the 90% confidence level limits on their values are still compatible with 0, |U µ4 | 2 ≤ 0.11 and |U τ4 | 2 ≤ 0.15. Combinations of the parameters shown in the axes above the contours are disfavoured. The figure includes a recent exclusion contour from Super-Kamiokande [148] as a comparison.
As in the previous cases, the presence of sterile neutrinos can also have an effect on neutrino spectra from astrophysical sources. In the case of a pseudo-Dirac scenario, where the mass splitting between active and sterile neutrinos is very small, there exist the possibility that the corresponding neutrino oscillation effects are still visible, i.e., L∆m 2 /E 1 [166,137].

Indirect Dark Matter Detection
There is a large corpus of evidence that supports the existence of a non-baryonic, non-luminous component of matter in the cosmos. A way to understand the rotation curves of galaxies, the peculiar velocities of galaxies in clusters, and the formation of first galaxies growing out of small density perturbations imprinted in the cosmic microwave background, is to introduce a "dark matter" component in the energy budget of the Universe [167,168]. Attractive candidates for dark matter consists of stable relic particles whose present density is determined by the thermal history of the early universe [169,170,171,172]. The present abundance of dark matter can be naturally explained by physics beyond the Standard Model providing stable, weakly-interacting massive particle (WIMP) in the few GeV-TeV mass range. For thermally produced WIMPs, the upper mass limit arises from theoretical arguments in order to preserve unitarity [173], although higher masses can be accommodated in models where the dark matter candidates are not produced thermally [174].
There is a vast ongoing experimental effort to try to identify the nature of dark matter through different strategies: production at colliders [175] or through the detection of nuclear recoils in a selected target in "direct detection" experiments [176]. A complementary, "indirect", approach is based on searching for the products of the annihilation of dark matter particles gravitationally trapped in the halo of galaxies or accumulated in heavy celestial objects like the Sun or Earth [177,178,179,180,181,182,183,184]. In this latter case, neutrinos are the only possible messengers, since other particles produced in the annihilations will be absorbed. These search techniques are competitive since they can set limits on the same physical quantities (the dark matter-nucleon cross section for example). But they are also complementary since they are subject to different backgrounds (the gamma-ray sky is very different from the proton or neutrino sky), different astrophysical inputs (dark matter density and velocity distribution) and different systematics (nucleon and nuclear form factors of different targets).
The strength of the expected neutrino flux emitted from a celestial object depends, among other factors, on the capture rate of WIMPs, which is proportional to the WIMP-nucleon cross section, and the annihilation rate, which is proportional to the velocity-averaged WIMP-WIMP annihilation cross section, σ A v . The evolution of the WIMP number density n DM in compact celestial objects follows the balance equatioṅ where Q C is the capture rate per unit volume and the numerical factors account for the fact that annihilations remove two WIMPs per interaction but there are only 1/2 of the WIMPS available to form pairs. If the capture rate remains constant over a long time, the WIMP density reaches an equilibrium solution The WIMP capture rate from interaction with baryonic matter can have spin-dependent, σ SD , and spinindependent, σ SI , contributions [185]. Since the Sun is 25 primarily a proton target (75% of H and 24% of He in mass) [186] the capture of WIMPs from the halo can be considered to be driven mainly via the spin-dependent scattering. Heavier elements constitute less than 2% of the mass of the Sun, but can still play a role when considering spin-independent capture, since the spin-independent cross section is proportional to the square of the atomic mass number. Note that these heavy elements can also take part in the spin-dependent capture process if WIMPs present momentum-dependent interactions with normal matter [187]. The situation for the Earth is rather different, since the most abundant isotopes of the main components of the Earth inner core, mantle and crust, 56 Fe, 28 Si and 16 O [188], are spin 0 nuclei. Searches for dark matter accumulated in the Earth are then more sensitive to the σ SI component of the WIMP-nucleon cross section. Another difference with respect to solar searches is that equilibrium between the capture and annihilation rates can not be assumed, and to be able to extract a limit on σ SI , an assumption on the value of the WIMP self-annihilation cross section must be made.
Dark matter searches from our own Galaxy, nearby galaxies or galaxy clusters present some distinct features with respect to searches from the Sun or Earth which are advantageous. Firstly, capture is not an issue since the presence of dark matter over-densities has been an essential part in the process of galaxy formation. What can be measured then is the velocity-averaged WIMP self-annihilation cross section, σ A v . Secondly, the products of the annihilations are not necessarily absorbed at the production site, and other indirect signatures (photons, anti-protons, etc.) can also be searched for in γ-ray and cosmic-ray observatories. These multi-wavelength and/or multi-messenger searches can increase the sensitivity of dark matter searches. Neutrinos remain, however, an attractive signature since they do not suffer from uncertainties in their propagation (as charged particles do) and no background or foreground from astrophysical objects is present (as in the case of γ-rays). Note, that some of the sources are extended (the Galactic halo for example) and point-source analysis techniques have to be modified. On the other hand we expect that the flux of secondaries from these distant objects is much lower than that predicted from WIMP annihilations in the Sun and, furthermore, there are new systematics effecting the calculations. For example, the assumed shape of the dark matter halo profile effects significantly the interpretation of the results since the annihilation rate depends on the square of the dark matter number density. We will discuss these issues in section 7.4.

Neutrinos from WIMP Annihilation and Decay
The annihilation of dark matter into Standard Model particles can be probed by γ-ray, cosmic ray, and neutrino emission. We focus in the following on neutrino emission with a spectral production rate per unit volume given by where ρ DM = m X n DM is the WIMP mass density at a given position r. The energy distribution dN ν α /dE ν of neutrinos is normalised to the total number of neutrinos of flavour ν α expected from the annihilation. The factor 1/2 compensates for the symmetry of WIMP combinations. The differential flux observed from the solid angle Ω can then be expressed as the line-of-sight integral from Earth's location 4 r ⊕ , where r(l, Ω) ≡ r ⊕ + l n(Ω) and n(Ω) is a unit vector in the direction Ω. The factor P α→β (E ν ) accounts for flavour oscillations 5 . The factor in parenthesis [·] encodes the cosmological and astrophysical dependence on the emission from the dark matter density distribution. This is sometimes called J-factor [189]. The decay of dark matter can be treated analogously. Here, the spectral production rate is given by where Γ X denotes the dark matter lifetime. As before, the distribution dN ν α /dE ν of neutrinos is normalized to the total number of neutrinos of flavour ν α and it depends implicitly on the branching fraction of the decay. The differential flux observed from the solid angle Ω can then be expressed as the line-of-sight integral Again, the factor in parenthesis [·] encodes the cosmological and astrophysical dependence and is sometimes called D-factor. IceCube analyses are performed in the most modelindependent way possible, but the allowed parameter space of the underlying theoretical models force some assumptions to be made. Since the exact branching ratios of WIMP annihilation into different channels is modeldependent, two annihilation channels which give extreme neutrino emission spectra dN/dE are chosen, and the analysis is optimised for each of them separately. Annihilation into bb is taken as a representative case for a soft neutrino energy spectrum. This is in part due to the bquark interaction with the medium (for dense objects like the Sun) and the hadronisation of the final state quarks leading to neutrinos from meson decays. In contrast, the annihilation into W + W − or τ + τ − follows a hard spectrum. Assuming a 100% branching ratio to each of these channels brackets the expected neutrino spectrum from any realistic model with branching to more channels. The WIMP mass is left free, and independent searches are performed for a few benchmark masses. The total number of signal events, µ s , expected at the detector from a given model is then given by, where T live is the exposure time, A eff ν α (Ω, E ν ) the detector effective area for neutrino flavour α, that depends on the detector's response with respect to observation angle and neutrino energy.
In the absence of a signal, the 90% confidence level limit on the number of signal events, µ 90 s , can then be directly translated into a limit on either the velocityaveraged annihilation cross section σ A v or the dark matter lifetime Γ X . The interplay between the total number of observed events in a given data sample, n obs , and the estimated number of background events, n bg , is the basis to perform a simple event-counting statistical analysis to constrain µ s , i.e., to constrain a given model. This was the approach followed in early IceCube publications, e.g., Refs. [190,191].
In order to improve the power of the statistical test, a distribution-shape analysis can be used instead. The angular distribution of the expected signal events is expected to peak towards the source, and is very different from the flatter atmospheric neutrino background. If Ω is the solid angle of a reconstructed track position, we can weight the event by signal and background probability distributions S(Ω) and B(Ω), respectively. The likelihood that the data sample contains µ s signal events out of n obs observed events is then defined as From this likelihood, one can determine confidence intervals for µ s by using the likelihood-ratio test statistic as proposed in [192]. This method produces stronger limits than just event counting, and it is less prone to unsimulated background contamination in the final data sample. There are systematic uncertainties in the translation of the number of detected events into capture cross section values due to uncertainties in the element composition of the Sun [193], the effect of planets on the capture of WIMPS from the halo [194] and the uncertainty on the values of the nuclear form factors needed in the rather complex capture calculations [185,195,196]. These effects Fig. 23 IceCube limits on the spin-dependent WIMP-proton scattering cross section, σ SD χ−p as a function of WIMP mass, compared to results from other neutrino detectors and direct detection experiments [211,203,212,213]. The IceCube limits have been scaled up to the upper edge of the 1σ systematic uncertainty band. The coloured points correspond to models from a scan of the pMSSM. The model points are shown colour-coded according to the "hardness" of the resultant neutrino spectrum. Red points correspond to models that annihilate predominantly into harder channels (such as τ + τ − ) and can hence be excluded by the IceCube red line, while blue points correspond to models that favour annihilations into softer channels (such as bb) and are probed by the blue lines. Similar coding applies for intermediate colours. Figure from [214].
can be of relevance when comparing results from different search techniques [197].
The experimental effort to detect neutrinos as a signature of dark matter annihilations in celestial bodies came of age in the mid 90's with underground detectors like MACRO, Baksan, Kamiokande and Super-Kamiokande. These detectors provided the first limits on the flux of neutrinos from dark matter annihilations in the Earth or the Sun [198,199,200,201]. Baksan and Super-Kamiokande continue to be competitive in the field today [202,203,204]. Baikal [205] and AMANDA were the first largescale neutrino detectors with an open geometry to perform dark matter searches in the late 90's, soon followed by ANTARES [206]. Early results of these experiments can be found in Refs. [205,207,208,209,190,210].

Dark Matter Signals from the Sun
The interplay between the dark matter capture and annihilation determines the number of WIMPs accumulated in the Sun. Losses through evaporation due to WIMPnucleus scattering have been shown to be negligible for WIMP masses above a few GeV [178,215,216] and can therefore be neglected in IceCube analyses. Given the age of the Sun (4.5 Gyr), the estimated dark matter density (ρ local ∼ 0.4 GeV/cm 3 ) and a weak-scale interaction between dark matter and baryons, many models predict that dark matter capture and annihilation in the Sun have reached equilibrium with density following Eq. 58. In this case, the J-factor in Eq. (60) becomes proportional to 1/ σ A v and total neutrino emission only depends  24 IceCube limits on the spin-dependent LKP-proton scattering cross section, σ SD χ−p as a function of LKP mass, compared with limits from direct detection experiments [221,222,223]. The theoretically allowed phase space is indicated by the green shaded [224]. The region below m LKP = 300 GeV is excluded by collider experiments [225], and the upper bound on m LKP arises from arguments to avoid over-closure in the early universe [226]. The lighter blue region is allowed when considering a wide range of dark matter relic density 0.05 < Ω CDM h 2 < 0.20, and the darker blue region corresponds to the 1σ WMAP relic density, 0.1037 < Ω CDM h 2 < 0.1161 [227]. on the capture Q C related to the WIMP-nucleon scattering cross section. Under the assumption that the capture rate is fully dominated either by the spin-dependent or spin-independent scattering, conservative limits can be extracted on either the spin-dependent, σ SD , or spinindependent, σ SI , WIMP-proton scattering cross section. The number of events in the detector observed in a given live time can therefore be translated into a limit on a physical quantity that can be used to compare with other experiments or to test predictions of a specific particle physics model. Different dark matter scenarios can be probed through the predicted neutrino flux in Eq. (63). Indeed, IceCube has set limits to the muon flux from annihilations of the lightest Kaluza-Klein mode arising in models of universal extra dimensions, as well as to its scattering cross section with protons [217] (see Fig. 24). Other, non-standard, scenarios like strongly interacting dark matter or self-interacting dark matter can also be tested since the IceCube event selections are quite generic and model-independent [218,219,220].
Traditionally, solar WIMP searches with IceCube have used the muon channel since it gives better pointing and, in the end, dark matter searches from the Sun are really point-source searches. The first analyses used the Earth as a filter of atmospheric muons and "looked" at the Sun only in the austral winter, when the Sun is below the horizon at the South Pole [217,228,229]. With the completion of IceCube-79 and DeepCore, it was possible to define effective veto regions to efficiently reject incoming atmospheric muons from above [230]. Since then the IceCube solar WIMP searches cover also the austral sum-mer, doubling the exposure of the detector per calendar year. DeepCore has also allowed to extend the search for neutrinos from WIMPs with masses as low as 20 GeV/c 2 , whereas past IceCube searches have only been sensitive above 50 GeV/c 2 . Additionally, all-flavour analyses are being developed [231], since the addition of ν e and ν τ events triples the expected signal. Improved low-energy reconstruction techniques allow to reconstruct electron and tau neutrino interactions with sufficiently good angular resolution to be useful in solar and Earth dark matter searches.
Signal and background differ, though, not only in their angular distribution, but also in their energy spectra. This information can be encapsulated in a likelihood function that includes a number-counting (normalisation) term and an angular and energy (spectral) terms, both for the signal and background p.d.f.'s. This is the approach taken in the latest IceCube solar WIMP analyses [232,214]. Allflavour analyses also benefit from this technique since the better energy resolution of cascade events benefits from the use of energy information in the likelihood. The general form of such an extended likelihood is, where the second term is analogous to Eq. (64) and the prefactor is the Poisson number likelihood for observing n obs events given a total of µ tot predicted signal and background events. The signal and background probability distribution functions, S and B, introduced in Eq. (64) now include the event's reconstructed solid angle Ω i and energy E i . For instance, the signal probability distribution is given by In general, the reconstructed neutrino energy E and solid angle Ω are different from the true energy E and true arrival direction Ω . The relation between these quantities on a statistical basis must be obtained from simulations, and it is specific to the annihilation channel under  Fig. 25 Improvement on the dark-matter proton cross section due to the use of an event-level likelihood analysis (full lines) compared to a traditional cut-and-count analysis (dashed lines) and an analysis based on the difference in shape of the space-angle distribution for signal and background (dotted lines tagged 'PRL' and originally from [230]). 100% annihilation to bb, W + W − and τ + τ − is assumed in each case shown. Figure from [232]. study. The true energy E can be related to the number of hit DOMs or can be estimated from more elaborate energy reconstructions. The quantity Q(E i , Ω i |E , Ω ) is the probability density (in effective units of inverse steradian and proxy energy) for reconstructing E i and Ω i for the ith event when the true values are E and Ω , respectively.
Note that systematic uncertainties on the signal and/or background prediction or on the angular or energy resolutions can be easily incorporated in a likelihood approach as nuisance parameters by marginalising over them. The only knowledge needed is the functional form of the nuisance parameters. P is a function of energy and angle, which in a simplified approach can be decomposed in an angle-dependent part (the PSF of the detector) and an energy-dependent part (the energy dispersion of the detector). More generally, the angular response of the detector can depend on energy, and then this decomposition is not valid. For point-source searches, due to the restricted angular region in the sky considered, the PSF and energy dispersion can be taken to be uncorrelated. Figure 23 shows the limits on the spin-dependent WIMP-proton cross section using three years of IceCube data, including both angular and energy terms in the likelihood [214]. The plot shows the latest IceCube results as full lines for three benchmark annihilation channels (in different colours), compared with the results of Super-Kamiokande [203] and ANTARES [211], as well as results from the PICO direct detection experiment [212,213]. The dots correspond to a parameter scan of the phenomenological Minimal Supersymmetric extension of the Standard Model (pMSSM) where the color code represents the leading annihilation channel: channels producing a soft neutrino spectra are marked as blue, and are probed by the soft experimental limits, while harder neutrino spectra are marked in red and are probed by the hard ex-perimental limits. Since large-volume neutrino telescopes are high-energy neutrino detectors, the limits for annihilation channels leading to a soft neutrino spectrum can be more than two orders of magnitude less restrictive than those resulting in harder spectra. Even in the latter case the limits can decrease rapidly for low WIMP masses. Direct search experiments do not reach σ SD values much below 10 −38 cm 2 at their best point, worsening rapidly away from it, while, together, the limits from IceCube and Super-Kamiokande cover the WIMP mass range between from a few GeV to 100 TeV and reach cross section values of about 10 −40 cm 2 . In the case of the σ SD cross section, direct-search experiments have the advantage of dedicated spinless targets, and the limits from neutrino telescopes lie about three orders of magnitude above the best limit from LUX at a WIMP mass of about 50 GeV [236].
The improvement due to using a full event-based likelihood in comparison to just an angular shape analysis is illustrated in Fig. 25, taken from [232]. The figure shows the limit on the spin-dependent WIMP-proton cross section as a function of WIMP mass obtained with two analyses performed on the same data set taken with IceCube-79. Shown are the limits using an event-count likelihood (dashed lines), the limits obtained using an analysis based on the difference in shape of the space-angle distribution for signal and background (tagged 'PRL' and originally from [230]) and the limits obtained using a full likelihood like in Eq. (66). Including the event-level energy information has the most impact at high WIMP mass, due to the relatively good energy resolution of IceCube at high neutrino energies. Note that the full likelihood analysis in [232] used a rather simple energy proxy based on the number of hit DOMs. Better energy reconstruction algorithms being developed within IceCube, particularly at low energies, will further improve the performance of this method [76].
There is an additional step in complexity when using neutrino telescope results to probe specific WIMP models, that avoids the need to simulate specific annihilation benchmark channels [237,238,235]. For a given model, this approach takes into account the full neutrino spectrum from WIMP annihilations including all allowed annihilation channels, i.e., the full d 2 P/dE dφ (E, φ, ξ) is calculated. The expected number of signal events is then obtained through Eq. (63) and, through the use of the likelihood function in Eq. (65), the model under consideration can be accepted or rejected at a given confidence level. To explore a wide parameter space, which is typical for supersymmetric extension of the Standard Model, each allowed combination of the free parameters of the model needs to be tested, and the procedure becomes computationally demanding. Even ad-hoc models like the constrained MSSM (cMSSM) [239] with just 7 parameters pose a computational challenge. But it is a powerful way of assigning a statistically meaningful weight to different areas of the model parameter space. The au-  thors in [237] and [235] have performed sensitivity studies of IceCube and IceCube-DeepCore under different model assumptions. Figure 26 shows the results of the model scan in [235]. The figure shows the spin-dependent neutralino-proton cross section versus neutralino mass. Each point in the plot represents an allowed combination of the cMSSM parameters. The colour code indicates which models can be probed by IceCube and disfavoured at the 1σ (green), 2σ (light blue) and 5σ (dark blue) level. Figure 27 illustrates the complementarity between accelerator and neutrino telescope searches for supersymmetry. The plot shows the gluino-squark mass parameter space with colour-coded exclusion levels by IceCube (same colour coding as in Fig. 26). The brown line shows the 95% CL exclusion region from searches for coloured sparticles in jets+missing transverse energy with, at the time, 4.71 fb −1 of data at centre-of-mass energies √ s 7 TeV from ATLAS [240]. As can be seen, there is a wealth of models not excluded by ATLAS that are under the reach of IceCube (blue dots). Note that more recent results from ATLAS and CMS (see, e.g., the summary on the experimental status of Supersymmetry in [41]) do not change the picture of complementarity between the parameter space reach of IceCube and the LHC.

Dark Matter Signals from the Earth
The rationale for searching for dark matter accumulated in the Earth follows a similar line than that of the Sun: WIMPs gravitationally accumulated in the centre of the Earth can annihilate, giving rise to a flux of neutrinos. A signal from dark matter annihilations in the centre of the Earth will produce a unique signature in IceCube as vertically up-going muons, where each string can act as a more or less independent detector. Searches from the vertical direction pose, however, some challenges in Ice-Cube due to its geometry. While in any other point source search an off-source region at the same declination of the source can be defined to measure the background, this is not possible for the vertical direction, and one needs to rely on accurate simulations of the background components (atmospheric neutrinos and muons). Through the use of advanced classification methods to separate signal and background, the amount of misreconstructed atmospheric muons can be reduced to a negligible level, and a likelihood shape analysis can be performed using Eq. (64). The results of such an analysis using 327 days of live time with the IceCube-79 configuration [241] are shown in Fig. 28. The shape of the IceCube limits reflects the resonant capture of WIMPs of certain masses that nearly match the mass of the main isotopes of the Earth (the peaks show the Iron, Silicon and Oxygen resonances at 53 GeV/c 2 , 26 GeV/c 2 and 15 GeV/c 2 , respectively). A self-annihilation cross section, σ A v , of 3×10 −26 cm 3 s −1 has been assumed, a typical value for a particle species to be a thermal relic. The lack of a signal can also be used to set limits on the dependence of σ SI on the annihilation cross section, as illustrated in Fig. 29. Values above the full lines are disfavoured at the 90% confidence level since the combination of capture and annihilation would have produced a detectable muon flux in IceCube. LUX limit m χ =50GeV IC86−I Earth limit : m χ =1TeV, χχ →W + W − IC86−I Earth limit : m χ =50GeV, χχ →τ + τ −

Fig. 29 Upper limits at 90% confidence level on σ SI
X+N as a function of the velocity-averaged annihilation cross section σ A v for 50 GeV WIMPs annihilating into τ + τ − and for 1 TeV WIMPs annihilating into W + W − . Limits from LUX [236] are shown as dashed lines for comparison. The red vertical line indicates the thermal annihilation cross section for a particle species to constitute the dark matter. Also indicated are Ice-Cube limits on the annihilation cross section for the respective masses and annihilation channels [244]. Figure from [241]

Dark Matter Signals from Galaxies and Galaxy Clusters
The Milky Way centre and halo, as well as nearby dwarf galaxies and galaxy clusters provide natural large-scale regions of increased dark matter density. Since dark matter played a significant role in the formation of such structures from primordial density fluctuations, the issue of capture is not relevant, and what neutrino telescopes can prove when considering such objects is the WIMP selfannihilation cross section, σ A v . In order to predict the rate of annihilation of dark matter particles in galactic halos, the precise size and shape of the halo needs to be known. There is still some controversy on how dark halos evolve and which shape they have. There are different numerical simulations, observational fits, and parametrisations of the dark matter density around visible galaxies, including the Navarro-Frenk-White (NFW) profile [245], the Kravtsov profile [246], the Moore profile [247], and the Burkert [248]. The common feature of these profiles is a denser spherically symmetric region of dark matter in the centre of the galaxy, with decreasing density as the radial distance to the centre increases. Where they diverge is in the predicted shape of the central region. Profiles obtained from N-body simulations of galaxy formation and evolution tend to predict a steep power-law type behaviour of the dark matter component in the central region, while profiles based on observational data (stellar velocity fields) tend to favour a constant dark matter density near the core. This is the core-cusp problem [249], and it is an unresolved issue which affects the signal prediction from dark matter annihilations in neutrino telescopes. Note that the shape of the dark halo can depend on local characteristics of any given galaxy, like the size of the galaxy [250] or on its evolution history [251,252] The shape of the dark matter halo is important because the expected annihilation signal depends on the line-of-sight integral from the observation point (the Earth) to the source, and involves an integration over the square of the dark matter density. This is included in the J-factor of Eq. (60), which is galaxy-dependent, and absorbs all the assumptions on the shape of the specific halo being considered. In the case of our Galaxy, the expected signal from the Galactic Centre assuming one halo parametrisation or another can differ by as much as four orders of magnitude depending on the halo model used (see, e.g., Fig. 2 in Ref. [260]).
As before, a few benchmark annihilation channels can be chosen to bracket the model expectations, and then dN ν /dE represents the neutrino flux assuming 100% annihilation to each of the benchmark scenarios. The analyses use the same likelihood approach as described by Eq. (64). The signal hypothesis (excess of events at small angular distances ψ to the Galactic Centre) can then be tested against the background-only hypothesis (an event distribution isotropic in the sky). There is, however, an additional effect to take into account when dealing with

Fig. 30
Comparison of upper limits on σ A v versus WIMP mass, for the annihilation channel χχ → τ + τ − . IceCube results obtained with different detector configurations [253,254,244,255] are compared to ANTARES [256] and the latest upper limits from γ-ray searches from H.E.S.S. [257] and from a combination of Fermi-LAT and MAGIC results [258]. Figure from [259]. extended sources, like the Galactic halo. Since the signal is allowed to come from anywhere in the halo, the background distribution, which is usually taken from data scrambled over azimuth angle, B = D, is necessarily contaminated by a potential signal. Thus, the background distribution B depends indirectly on the number of signal events µ S and needs to be corrected. The effective background distribution can be written as where S and D are the probability distribution functions of the azimuth-scrambled arrival directions of signal simulation and data events respectively. There is another, complementary, analysis approach for extended sources that naturally incorporates the diffuse character of the signal over a large region in the sky. It is based on a multipole expansion of the sky map of event arrival directions. Dark matter annihilations in the halo would produce a diffuse flux of neutrinos with a characteristic large scale structure following the shape of the halo, while the atmospheric neutrino background presents small anisotropies at smaller scales. The sky map of reconstructed event arrival directions can be constructed as where (θ i , φ i ) are the reconstructed coordinates (declination and right ascension respectively) of event i, n obs is the total number of events in the data sample and δ D is the Dirac-delta-distribution. Such distribution can be mapped onto an expansion in spherical harmonics on the sphere, where the expansion coefficients a m are given by the sum over events with reconstructed arrival directions Ω i , Note that this expansion depends on the coordinate system. In particular, in the equatorial coordinate system where the event distribution has strong dependence on declination angle (equivalent to zenith angle) from the detector acceptance, the expansion is dominated by the m = 0 coefficients. It is possible to design a test statistic of the remaining m = 0 components to separate signal from background (see Eqs. (8) and (9) in [253]). Results from the searches performed by IceCube with different techniques on the Galactic Centre, halo and dwarf spheroid galaxies [253,254,244,255] are shown in Fig. 30, compared with other experiments and theory interpretations. All sources considered showed results compatible with the background-only hypothesis yielding limits on the velocity-averaged annihilation cross section at the level of few 10 −23 cm 3 s −1 . Recent results from γ-ray telescopes on dwarf spheroids currently provide the best limits on σ A v , due to their accurate pointing and lack of foreground or background from these kind of sources.
The high-energy diffuse astrophysical neutrino flux discovered by IceCube opens a new possibility of probing the galactic dark matter distribution through neutrinodark matter interactions [262,263,261]. Indeed dark matter couplings to standard model particles are commonly assumed to exist, and this is the basis of direct detection experiments. Such a coupling can be extended to neutrinos if one assumes the existence of a mediator φ, which can be either bosonic of fermionic in nature, which couples to dark matter with a coupling g. For simplicity one can assume that the strength of the ν − φ coupling is also given by g. Under these assumptions, dark matterneutrino interactions could distort the isotropy of the astrophysical neutrino flux, resulting in an attenuation of the flux towards the Galactic Centre, where the density of dark matter is higher. Therefore, an analysis of the isotropy of the high-energy astrophysical neutrino data of IceCube can be translated into a limit on the strength of the neutrino-dark matter coupling, g. Such an analysis has been performed by the authors in [261], where the mass of the dark matter candidate, the mass of the mediator and the strength of the coupling are left as free parameters in a likelihood calculation which aims at evaluating the suppression of astrophysical neutrino events from the direction of the Galactic Centre. The results are shown in Fig. 31. The left panel shows contours of the  Fig. 31 Contours of the maximum allowed value of the coupling of dark matter to neutrinos, g max , as a function of the dark matter mass, m χ , and the mediator mass m φ . Left Panel: fermionic dark matter coupled through a vector mediator. Right Panel: scalar dark matter coupled through a fermionic mediator. The magenta line in both plots delimits the region where cosmological observations from large scale structure become more restrictive than the IceCube limits. The plots have been obtained by marginalising over the atmospheric and astrophysical fluxes, allowing the astrophysical spectral index to vary between 2 and 3. The diamond and the stars refer to specific models studied in [261]. Reprinted with permission from [261]. Copyright 2017 by the American Physical Society. maximum allowed value of the coupling of fermionic dark matter coupled to neutrinos through a vector mediator, while the right panel shows the case of a scalar dark matter coupled through a fermionic mediator. Interestingly IceCube is sensitive to a region of parameter space complementary to results derived from cosmological arguments alone [264,265,266], as indicated by the magenta line. This line delimits the region where limits from analyses using large-scale structure data become more restrictive than the IceCube limits shown in the plot.

TeV-PeV Dark Matter Decay
The origin of the TeV-PeV diffuse flux observed with Ice-Cube is so far unknown. Whereas most models assume an astrophysical origin of the emission, it is also feasible that the emission is produced via the decay of dark matter, as first proposed in Ref. [267]. Various studies have argued that heavy dark matter decay can be responsible for various tentative spectral features in the inferred neutrino spectrum of the HESE analysis [268,269,270,271,272,273,274] and also its low energy extension [275,276,277]. Some authors have also discussed the necessary condition for PeV dark matter production in the early Universe, e.g., via a secluded dark matter sector [278], resonantly-enhanced freeze-out [279,280], or freeze-in [281,282,283,284].
If dark matter decay is responsible for the high-energy neutrino emission, the arrival directions of TeV-PeV neutrino events observed with IceCube should correlate with the line-of-sight integral of the dark matter distribution ("D-factor"). The contribution of neutrinos from dark matter decay in the Galactic halo can be similar to the isotropic extra-galactic contribution. Half of the neutrino events from dark matter decay in the halo are predicted to fall within 60°around the Galactic Centre. This introduces a weak large scale anisotropy that can be tested against the observed event distribution [289,290,291]. The neutrino emission from galaxies and galaxy clusters could also be identified as (extended) point-source emission in future IceCube searches [292].
It can be expected that the secondary emission from decaying dark matter scenarios will also include other standard model particles that can be constrained by multimessenger observations. In particle the production of PeV γ-rays, that have a pair production length of O(10) kpc in the CMB, would be a smoking gun for a Galactic contributions [292]. However, also the the secondary GeV-TeV emission from electro-magnetic cascades initiated by PeV γ-rays, electrons, and positrons can constrain the heavy dark matter decay scenario [292,290,293,294,295].
From an experimental point of view, the search for a neutrino signal from heavy dark matter decay follows closely the strategies used for the search for annihilation signatures. A given analysis can be used for both constraining dark matter lifetime and annihilation, as can be inferred from the similarities between Eqs. (60) and (62). A feature of searches for dark matter decay in the 100 TeV mass region is that the recently discovered astrophysical neutrino flux becomes a background to the search, and its current uncertainties on normalisation and the presence, or not, of a cut-off effect the interpretation of any search for an additional component due to dark matter decay. IceCube has performed a search for signatures of heavy dark matter decay by assuming that the background is a combination of atmospheric and astrophysical neutrinos, where the normalisation and energy dependence of the astrophysical flux is allowed to flow freely [288]. This allows for deviations from a strict power law which could be interpreted as an additional contribution from dark matter decays. The signal consists of both a galactic component and a diffuse component  Right Panel: IceCube lower limits on dark matter lifetime versus dark matter mass assuming a 100% branching ratio decay to bb (full line), µ + µ − (dashed line), compared with limits from HAWC, both for Dwarf Spheroidal Galaxies [285] and the Galactic Centre/Halo [286], and Fermi/LAT [287]. Figures from [288].
from the contribution of distant galaxies. Typical neutrino spectra from the decay of a 2 PeV dark matter candidate, assuming 100% decay to each channel, are shown in the left panel of Fig. 32, featuring the characteristic peak at the dark matter mass from those decays where the neutrino takes most of the initial available energy. Recent results from IceCube on the dark matter lifetime as a function of the dark matter mass are shown in the right panel of Fig. 32, compared with results from HAWC and Fermi/LAT. Values below the lines are disfavoured at 90% confidence level, with IceCube limits showing the strongest constraint for masses above about 10 TeV.

Magnetic Monopoles
Maxwell's equations of classical electrodynamics appear to be asymmetric due to the absence of magnetic charges. However, this is merely by choice. We can always redefine electric and magnetic fields by a suitable duality transformation, E + iB = e iφ (E + iB), such that Maxwell's equations in terms of the new fields are completely symmetric. However, this duality transformation requires that for every particle the ratio between magnetic and electric charges are the same. If this is not the case, then it is necessary to include a source term for magnetic charges q m (i.e., magnetic monopoles), that create a magnetic field of the form B = q m e r /4πr 2 . P. Dirac was the first to speculate about the existence of magnetic monopoles, guided 7 by "new mathematical features" appearing in the quantum-mechanical description of electrodynamics [297]. He showed that the presence of a magnetic monopole with a minimum charge q m would explain why the electric charge is always quantised, i.e., in integer multiples of an elementary charge.
His argument was based on the observation that the magnetic potential of the static magnetic monopole field, B = rot A, has to be singular along a semi-infinite line. The wave function of a particle with charge q e encircling this line will pick up a phase that is equal to ∆φ = ±q e q m . To be unobservable, this phase should be a multiple of 2π independent of the type of monopole or electrically charged particle. This leads to the condition 8 q m q e = 2πn (n ∈ Z) .
At the time of Dirac's seminal work, the smallest electric charge unit was considered to be that of the electron, q e = e, and hence the smallest possible (Dirac) magnetic charge would be g D = 2π/e = e/2α 68.5e .
This elementary monopole charge derives from basic considerations still valid in modern physics. In the Standard Model the smallest charge is provided by the down-type quarks, increasing the minimal magnetic charge to 3g D . Note that the large magnetic charge of monopoles predicted by the Dirac quantisation condition (72) leads to distinct electromagnetic features as monopoles pass through matter. We will return to this point shortly.
In solids, structures have been found which resemble poles. These are sometimes mistakenly called magnetic monopoles, although the poles can only occur in pairs and do not exist as free particles. For distinction, recently the term magnetricity has been coined for the field that exhibits theses poles. Fundamental magnetic monopoles have not been observed so far.
Although, as Dirac showed, magnetic monopoles can be consistently described in quantum theory, they do not appear automatically in that framework. As was first found independently by 't Hooft [299] and Polyakov [300], this is different in Grand Unified Theories (GUT) which embed the Standard Model interactions into a larger gauge group. These theories are motivated by the observation that the scale-dependent Standard Model gauge couplings seem to unify at very high energies. Generally, a 't Hooft-Polyakov monopole can arise via spontaneously breaking of the GUT group via the Higgs mechanism. The stability of the monopole is due to the Higgs field configuration which cannot smoothly be transformed to a spatially uniform vacuum configuration. An early candidate for a GUT theory studied in this context has been the Georgi-Glashow model [301], where matter is unified in SU(5) representations, that spontaneously break to the Standard Model gauge representations below the unification scale Λ GUT 10 15 GeV. This is related to the mass of the monopole as where α GUT = g 2 /4π is the gauge coupling constant at the unification scale. If the original unified group undergoes secondary symmetry breaking at lower energies also monopoles with lower masses may be generated. Depending on details of the GUT model, the monopole mass can range from 10 7 GeV to 10 17 GeV [298,302].  Fig. 33). The presence of virtual particles within these "shells" influences the monopole's interaction with matter. For instance, within the monopole core, the GUT gauge symmetry is restored, and can mediate baryon-number violating processes. At large distances, only electromagnetic interactions are visible by the magnetic monopole field.
Although the greatest interest has been in the supermassive monopoles that are motivated by GUTs, with the advent of the LHC and the MoEDAL experiment [303] the possibility of lighter monopoles lately received renewed attention. In the electroweak theory, with SU(2) × U(1) broken to U(1), there are no topologically nontrivial configurations of the Higgs field, and hence no topologically stable monopole solutions exist. However, there exist specific modifications with a non-minimal coupling of the Higgs field that would allow electroweak monopoles with TeV-scale masses, e.g., Refs. [304,305].

Cosmological Bounds
While inaccessible to collider experiments, very heavy GUT monopoles could have been produced in the early universe if the temperature exceeds T cr ∼ Λ GUT . In this case, the expected monopole density would be roughly one per correlated volume, which corresponds to the horizon size in a second-order phase transition [306,307]. This gives a naive monopole energy density (relative to the critical density) of This prediction is in clear conflict with the observed spatial flatness of the Universe (Ω tot 1) and known as the "monopole problem". An elegant solution to this problem is an inflationary universe, i.e., a universe that underwent an exponential expansion of the scale factor, diluting any initial monopole abundance to an (almost) unobservable level. This inflationary mechanism is a very powerful idea since it simultaneously explains why our Universe has been extremely flat at early times (flatness problem), e.g., Ω − 1 10 −16 at the epoch of big bang nucleosynthesis, and why the Universe appears to be so homogenous over causally disconnected distances (Horizon problem), e.g., temperature fluctuations in the CMB of only 10 −5 .
The requirement that the contribution of monopoles in Eq. (74) does not exceed today's dark matter abundance, Ω GUT Ω m , results in the overclosure bound on the integrated isotropic flux

Parker Bound
Magnetic monopoles are accelerated in magnetic fields -analogously to charged particle acceleration in electric fields. Therefore, relic monopoles that are initially nonrelativistic are expected to gain energy while they travel along galactic and intergalactic magnetic fields. The requirement that monopoles have to be rare not to shortcircuit these magnetic fields gives the so-called Parker bound [308]. The galactic magnetic field with a strength of a few µG can be generated by a dynamo action on a time scale that is comparable to the Milky Way's rotation period, τ 10 8 yr. A monopole with magnetic charge q m will gain an energy of ∆E kin = ∆ Bq m after it travels a distance ∆ along magnetic field lines. The power density of the galactic dynamo ∼ B 2 /τ should be larger than the energy drained by the magnetic monopoles.
During its acceleration the monopole encounters different magnetic field orientations coherent over a length scale λ c which are much smaller than the typical size r of the magnetic halo. If the monopoles stay non-relativistic, i.e., M m q m Bλ c , the energy gain is always large compared to the kinetic energy and the particle will be accelerated. The processes can be estimated by a random walk with N r/λ c encounters with coherent field regions. For weak magnetic fields, i.e., M m q m Bλ c this process loses its efficiency since the monopole does not follow the potential drop along the field lines. A careful analysis for our own galaxy gives the following bound for masses M m 10 17 GeV [309] Very massive monopoles M m 10 17 GeV will not be significantly deflected by the galactic magnetic field, since the acceleration across the galaxy does not change much of the initial virial velocity v ∼ 10 −3 c. The energy drain of the field by these monopoles depends thus on the initial monopole trajectory with respect to the field lines. To first order, the effect from the motion of these heavy monopoles and their anti-monopoles across the galaxy will cancel. However, the effect will be visible at second order which introduces a mass dependence ∆E kin ∝ M m . The Parker bound beyond M m 10 17 GeV is hence weaker than (77) by a factor (M m /10 17 GeV). Applying Parker's arguments to the seed magnetic fields of galaxies or galaxy clusters leads to tighter bounds [310,311]. As these bounds are less secure, and for consistency with other literature on the subject, we compare the experimental results with the original Parker bound (see Fig. 37).   35 Light yield for the different radiation mechanisms of magnetic monopoles [313]. For comparison the direct Cherenkov light emitted by a minimum ionising muon is shown.

Nucleon Decay Catalysis
The central core of a GUT monopole (see Fig. 33) contains the fields of the superheavy gauge bosons that mediate baryon number violation, so one might expect that baryon number conservation could be violated in baryonmonopole scattering and the possibility that a proton or a neutron in contact with a GUT monopole can decay. This feature was pointed out by Callan [314] and Rubakov [315]. The cross sections for the catalysis processes such as p + M → e + + π 0 + M (Fig. 34) are essentially geometric: The correction F(β) = β β 0 γ takes into account an additional angular momentum of the monopole-nucleus system. Both parameters γ and β 0 depend on the nucleus [316]. Current estimates for the catalysis cross sections are of the order of 10 −27 cm 2 to 10 −21 cm 2 [317].

Monopole Searches with IceCube
The experimental search for magnetic monopoles has a long history. Searches were pursued at accelerators, in cosmic rays, and for bound monopoles in matter. Detection methods include induction in SQUIDs, the observation of excessive energy loss of the monopole compared to Standard Model particles and particles describing nonhelical paths in a uniform magnetic field, or other unusual trajectories like non-relativistic velocities combined with a high stopping power and long ranges. In indirect searches at accelerators virtual monopole processes are assumed to influence the production rates of final states. In a direct search, evidence of the passage of a monopole through material is sought. Here we primarily address direct searches for primordial monopoles with IceCube.  [318]. The total time of the event is 5000 ns. Less light has been detected in a horizontal plane [4] roughly at the half height of the detector due to dust in the ice. Middle Panel: monopole at 0.3c moving from the top of the detector to the bottom emitting luminescence light. The simulated track of the particle is indicated in red. Only a few noise hits contribute throughout the detector due to the short time frame of the event [318]. Right Panel: Monopole with β = 10 −3 catalysing nucleon decay with λ cat = 1 cm with superimposed background noise [319]. The black line represents the simulated monopole track.
Different light production mechanisms induced by monopoles dominate depending on their velocity (see Fig. 35): (i) Direct Cherenkov light is produced at highly relativistic velocities above 0.76 c as with any other charged Standard Model particle. Due to the high relative Dirac charge, as shown in Eq. (72), several thousand times more light is radiated with a monopole than from a minimum ionising singly electrically charged particle like a muon [320].
(ii) Indirect Cherenkov light from secondary knock-off δ-electrons is relevant at mildly relativistic velocities ( 0.5 c to 0.76 c). The high-energy δ-electrons in turn can have velocities above the Cherenkov threshold themselves. The energy transfer of the monopole to the δelectrons can be inferred from the differential cross section calculated by Kasama, Yang and Goldhaber (KYG) [321,322].
(iii) Luminescence light from excitation of the ice dominates at low relativistic velocities ( 0.1 c to 0.5 c). The observables of luminescence, such as the wavelength spectrum and decay times, are dependent on the properties of the medium, in particular, temperature and purity. The signature is relatively dim in comparison to muon signatures. Pending further laboratory measurements in ice and water [323], the efficiency of luminescence photon production per deposited energy is in the range of dN γ /dE = 0.2 γ/MeV and 2.4 γ/MeV [324,325]. Even for the lower plausible light yield, luminescence is a viable signature due to the high excitation of the medium induced by a monopole [313].
(iv) At velocities well below 0.1 c luminescence is expected to fall off (see Fig. 35). The catalysis of nucleon decays is a plausible scenario for GUT monopoles (see Sec. 8.3) and may be observed if its mean free path is small compared to the detector size. The Cherenkov light from secondaries emitted in nucleon decays along the mo-nopole trajectory can lead to a characteristic slow moving event pattern across the detector [319].
For each of these speed ranges, searches for magnetic monopoles at the IceCube experiment are either in progress (luminescence) or have already set the world's best upper limits on the flux of magnetic monopoles over a wide range of velocities. Examples of magnetic monopole passing through the detector at different velocities are shown in Fig. 36.
These unique signatures from monopoles lead to the following general analysis strategies in the different velocity regions: (i) Relativistic monopoles are selected based on their brightness, arrival direction, and velocity [320]. The high energy astrophysical neutrino flux is an important background to this signal. At similar brightness, monopoles show less stochastic energy loss than Standard Model particles leading to a smoother light yield distribution along the track.
(ii) Due to its lower rest mass a Standard Model particle with a velocity below the speed of light in vacuum, c, would not be able to traverse the whole detector. However, the discrimination power of the reconstruction of the velocity is insufficient for the suppression of the vast air shower backgrounds against the identification of mildly relativistic monopoles. Instead variables describing the topology, smoothness, and time distribution of the events are processed in a Boosted Decision Tree (BDT) machine learning [320]. (iii) Searching for mildly relativistic monopoles using luminescence light can be performed using analysis techniques that combine the non-relativistic reconstructed particle velocity and the continuous but dim light production of a through-going track in the detector [313,318].  log 10 γ Fig. 37 Upper limits on the flux of magnetic monopoles as a function of the velocity of the monopole at the detector [313]. Shown are limits from IceCube [319,326,320], Baikal [324], ANTARES [327,328], RICE [329], ANITA [330] and Auger [331].
of electromagnetic and hadronic cascades. Because of the high light yield this channel is typically used as a benchmark in analyses. Due to their low speed, the duration of the event is in the order of 10 milliseconds. As obvious from the right event display in Fig. 36 at such timescales random noise pulses are a significant contribution. Various effects contribute to subtle temporal correlations on long time scales of this noise [332,4] complicating an adequate description in simulation. Instead, a background model is established from reshuffling experimental data. For signal identification, time-isolated local coincidences in neighbouring DOMs are searched for along a monopole trajectory hypothesis consistent with a straight particle track of constant speed. A Kalman filter is used to separated noise from monopole signals and a combination of observables are fed to a BDT to further improve the signal purity [319]. Figure 37 shows a compilation of current flux upper limits of relic monopoles from various experiments. Only in the past decade astrophysical experiments have been able to improve upon the original Parker bound which is shown for comparison. These recent experiments have employed large scale detectors for cosmic rays to achieve the highest sensitivities in the whole β range in which GUT magnetic monopoles are expected. Typically it is assumed that the flux at the respective detector site is isotropic implying sufficient kinetic energy in order to cross the Earth or the overburden above the detector due to their large rest masses. This assumption is justified for monopoles of masses in excess of 10 10 GeV. The monopole flux limits commonly assume a single Dirac magnetic charge q m = g D (see Eq. (72)) with no additional electric charge. In most detectors and velocity ranges, the detection efficiency for larger magnetic charges or for electrically charged monopoles (dyons) is expected to increase.
Operational until 2000, MACRO searched for magnetic monopoles using three types of subdetectors -liquid scintillation counters, limited streamer tubes and nuclear track detectors. No monopole was found, with an upper flux limit at the 90% confidence level of 1.4 × 10 −16 cm −2 s −1 sr −1 for monopoles with velocity between 4 × 10 −5 c to c and magnetic charge with n ≥ 1 [333,334]. Under the assumption that monopoles are gravitationally accumulated in the centre of the Sun, Super-Kamiokande [335] could impose stringent limits for nonrelativistic velocities in an indirect search for neutrinos from the direction of the Sun. Baikal [324] has investigated the direct Cherenkov light from relativistic monopoles. The analysis by ANTARES includes also the mildly relativistic regime employing similar techniques like IceCube [328]. The reduced scattering in water at the ANTARES site compared to ice leads to a better velocity reconstruction which helps with the Standard Model background suppression. This partly compensates the higher noise level in the detector.
Intermediate and low mass monopoles may acquire highly relativistic velocities in intergalactic magnetic fields reaching Lorentz factors of γ 10 10 for the example of a PeV-mass [302]. Ultra-relativistic particles with magnetic charge (or large electric charge) dramatically loose energy in their passage through matter, initiating a large number of bright showers along the track. At high Lorentz boost factors the photo-nuclear effect is the dominant energy loss mechanism generating hadronic showers. While these showers are continuously produced, they may overlap with each other. In the atmosphere of the Earth this leads to a built-up such that the energy deposit increases with slant depth. The Auger experiment has used this feature to distinguish monopoles from the background of electrically singly charged ultrahigh-energy cosmic rays like protons [331]. The RICE [329] and ANITA [330] exper-iments have searched for such multiple sub-shower signatures in the Antarctic ice sheet with the Radio-Cherenkov technique, the discriminant against conventional cosmic rays here being primarily the rapid succession of several radio pulses received from each sub-shower.
The extrapolation of IceCube's limit towards highly relativistic velocities by a constant line in Fig. 37 is a very conservative approach. It not only neglects the increase in signal detection efficiency with more energy deposited, but it also ignores the onset of the photo-nuclear effect and pair production. These effects would produce showers with light emission orders of magnitude brighter than from the Cherenkov effect considered here, hence visible also from far outside the instrumented detector volume.
While these flux limits reflect the cosmic density, at the electroweak scale monopoles may be created in accelerator collisions, which is studied at the MoEDAL experiment. Also cosmic ray collisions or high energy neutrino interactions in the Earth may produce monopoles. This might be an additional detection opportunity, also for IceCube.

Other Exotic Signals
The detection principle of Cherenkov telescopes is very general in the sense that it applies to any flux of particles that can penetrate the detector shielding and produce light signals inside the detector volume. We have already covered the possibility to observe relic magnetic monopoles in the previous section. In this section we will discuss the detection potential for other exotic candidates like Q-balls and strangelets. We will also address the possibility that long-lived charged particles produced in cosmic ray or neutrino interactions may be discovered via their Cherenkov emission. All these particles have in common that their passage through the IceCube detector produces observable features that can be extracted from backgrounds.

Q-Balls
There exist non-topological solitonic solutions of a field theory, so-called Q-balls [336]. Whereas the stability of topological solitons, e.g., monopoles is guaranteed by the conservation of a topological charge (winding number) associated with a degeneracy of the vacuum state, a Q-ball can be stable by the conservation of a charge associated with a global symmetry of the theory. This can happen if its energy configuration is lower than the corresponding multi-particle Fock state. For a single complex scalar field φ carrying the charge Q, this implies that there exists a non-trivial field value |φ 0 | > 0 where the scalar potential U(φ † φ) obeys the condition U(φ † 0 φ 0 ) < m 2 φ φ † 0 φ 0 , where m φ is the mass of the scalar field (for a review see Ref. [337]).
The appearance of these flat scalar potentials is generic in supersymmetric (SUSY) field theories, predicting scalar partners for fermions and gauge bosons, that may carry global charges [338]. Supersymmetry has to be broken at low energies to provide mass terms for the SUSY partners (for a review see, e.g., Ref. [339] and references therein). Since a naive spontaneous SUSY breaking by renormalisable terms in the visible sector predicts SUSY masses that are in general too light (below TeV), it is typically assumed that the breaking occurs in a hidden sector at some unobservable high-energy scale. The mediation of SUSY breaking terms to the visible sector, e.g., by gauge interactions or supergravity generates then soft (renormalisable) SUSY breaking masses for SUSY partners of the Standard Model matter and gauge bosons. Generically, SUSY extensions of the Standard Model predict flat directions for some combination of scalar fields. For instance, gaugemediated SUSY breaking is expected to introduce a scalar potential with V ∼ m 4 F (log(φ † φ/M 2 m )) 2 (for |φ| M m ) with mass scale m and messenger mass M m m. In this case, the mass of a Q-ball with Q 1 can be estimated as [340] m Q 4π Naturally, the effective mass scale m is expected to lie within the range 100 GeV < m < 100 TeV. If baryogenesis proceeded via the Affleck-Dine mechanism [341] 9 , stable Q-balls with 10 12 Q B 10 30 could have been formed copiously as a dark matter contribution by the fragmentation of the Affleck-Dine condensate [342]. It is an appealing property of this scenario that the baryon and dark-matter abundance, Ω B and Ω DM , respectively, are related. For Q B 10 26 one obtains Ω DM 10Ω B close to the observed values.
Analogously to monopole densities, Q-ball densities are also limited by the dark matter abundance. The maximal contribution of Q-balls to the observed dark matter Ω Q h 2 Ω m h 2 results in a bound on the integrated isotropic flux since 4πF Q v M Ω m ρ cr /m Q and Here the mass term (79) from gauge-mediated SUSY breaking is used. If we assume that dark matter Q-balls cluster in our galaxy with ρ halo 10 −24 g cm −3 we obtain the limit corresponding to a few events per year and square kilometer in a 2π sky coverage. The global charge Q associated with the Q-ball can be the same as baryon number (B) or lepton number (L) or some combination of them if these symmetries are connected to a global U(1) symmetry. Since the global symmetry is spontaneously broken in the interior of the Q-ball with φ = 0, the soliton could catalyse nucleon decay traversing the detector volume. This is analogous to the case of monopole-catalysed nucleon decay. Even if the charge is not related to B or L, it is possible that the vacuum state associated with the Q-ball interior catalyses nucleon decay. This can happen if the scalar potential is very flat such that Interactions that lead to nucleon decay induced by new physics at an ultra-violet scale Λ, for instance, in grand unified theories are typically suppressed by powers of φ /Λ and can become large in the Q-ball environment [343].
The cross section of the catalyzed nucleon decay via Q-balls with baryon number Q B passing through matter can be approximated by its size. In gauge-mediated SUSY breaking, it can be estimated as [340] σ cat 10 −20 Q B 10 26 The experimental signature of this process would thus be similarly spectacular as in the case of nucleon-decay catalysed by monopoles. Hence, the upper limits discussed for slow monopoles can be reinterpreted in terms of electrically neutral Q-balls also for neutrino telescopes. An example of such a recalculation is available for the Baikal experiment in [344,345] leading to a flux upper limit of F = 3.9 × 10 −16 cm −2 sr −1 s −1 at β = 10 −3 for an assumed cross section of σ cat = 1.9 × 10 −22 cm 2 . Since the detection efficiency for Q-balls is comparable to that for slow monopoles catalysing nucleon decay with the same cross section it is to be expected that a reinterpretation of the above mentioned IceCube slow monopole limits would also present an improvement of about two orders of magnitude with respect to other limits.

Strange Quark Matter
Using energy and symmetry arguments, it has been speculated that strange quark matter (SQM), a hypothetical form of matter with roughly equal numbers of up (u), down (d), and strange (s) quarks, could be the true ground state of QCD [346,347,348]. For a plasma of quarks in thermodynamical equilibrium it might be energetically preferable to condense into a phase containing strange quarks instead of ordinary matter with neutrons (udd) and protons (uud). An approximate thermodynamical calculation with massless quarks and neglecting strong interactions shows that the average kinetic energy per quark in ordinary bulk matter could be reduced in bulk SQM by a factor of about 0.89 [347]. Therefore, it is feasible that the extra "penalty" paid by the presence of more massive strange quarks is over-compensated by the reduction in energy density. However, the meta-stable state of protons, neutrons, and composite nuclei would be very long-lived, since conversion to the SQM ground state would proceed via weak interactions.
Lumbs of SQM, so-called strangelets, can have a large atomic mass number A and charge Z. Classical strangelets have a quark charge Z ∼ 0.1A for low mass numbers (A 700). For total quark charges exceeding Z ∼ α −1 ∼ 137 strong field QED corrections lead to screening and Z ∼ 8A 1/3 (A 700) [348]. It has also been speculated that colour and flavour symmetries at high baryon densities might be broken simultaneously by the condensation of quark Cooper pairs [349]. In this scenario, the "colour-flavour-locked" strangelets have charges of Z ∼ 0.3A 2/3 [350].
Stable strangelets can absorb ordinary matter in exothermic reactions involving u → s + e +ν e or u + d → s + u [348]. If the strangelet carries a positive charge the Coulomb barrier will usually prevent a strangeletnucleus system from collapse. However, neutron-rich environments like neutron stars are not protected by this mechanism. In fact, if strange quark matter is stable then all compact stars like white dwarfs or neutron stars are likely to consist of it. Even the capture of a single strangelet would be sufficient to convert a neutron-rich environment very rapidly.
Slowly moving strangelets -so called nuclearities [351] -lose their energy in matter due to atomic collisions. The excessive energy released will heat the medium and create thermal shocks. The hot expanding plasma will emit Planck radiation from its surface over a wide range of frequencies. IceCube is sensitive to optical photons energies of about 2 ÷ 4 eV (300 nm to 600 nm). At nuclearite velocities expected for cold dark matter the fraction of total energy loss emitted in optical photons is of order 10 −5 c. Since this leads to signatures similar to slow magnetic monopoles or Q-Balls the detection efficiency is comparable again. Nuclearites with masses in excess of 10 14 GeV and typical velocities of order of 10 −3 c reach underground detectors. Correspondingly, the sensitivity to the flux of nuclearites is roughly of the same order as that to slow monopole fluxes catalysing nucleon decays. Results for the MACRO experiments [352] and ongoing studies for the ANTARES detector [353] addressing fluxes of order 10 −16 − 10 −15 cm −2 sr −1 s −1 underline the high potential for a corresponding reinterpretation of IceCube analyses.

Long-Lived Charged Massive Particles
Many extensions of the Standard Model predict the existence of long-lived charged massive particles (CHAMPs). These particles occur naturally in scenarios where the decay of charged particles is limited by (approximate) discrete symmetries and involves final states that have only very weak couplings. Analogously to muons, these CHAMPs have a reduced electro-magnetic energy-loss in matter due to the suppression of bremsstrahlung by the rest mass. Still, they may be detected by their Cherenkov emission and, due to the long range, even with an enlarged effective detection area.
In the following we will consider SUSY breaking scenarios where the right-handed stau is the next-to-lightest SUSY particle (NLSP). However, most of the arguments apply equally well to other scenarios of CHAMPs with a decay length larger than other experimental scales (see, e.g., Ref. [354]). If R-parity is conserved the stau NLSP can only decay into final states containing the LSP. Depending on the mass and coupling of this particle the lifetime of the stau NLSP can be very long and, in some cases, it can be considered as practically stable on experimental time-scales.
In the case of a neutralino LSP, the stau NLSP can be very long-lived if its mass is nearly degenerated with that of the neutralino [355]. However, there are also superweakly interacting candidates for the LSP in extensions of the MSSM, which provide the long life-time of the NLSP more naturally. Possible scenarios include SUSY extensions of gravity with a gravitino G LSP, SUSY versions of the Peccei-Quinn axion and the corresponding axino LSP [356], and the MSSM with right-handed chiral neutrinos and a right-handed sneutrino LSP.
Staus produced in SUSY interactions of EHE neutrinos could traverse Cherenkov telescopes at the level of a few per year, assuming that extragalactic neutrino fluxes are close to the existing bounds [357,358,359,360,361,362] or prompt atmospheric neutrino fluxes close to upper theoretical estimates [361]. The leading-order SUSY contribution consists of chargino χ or neutralino χ 0 exchange between neutrinos and quarks, analogous to the partonlevel SM contributions shown in Fig. 11. The reactions produce sleptons and squarks, which promptly decay into lighter stau NLSPs.
Cosmic ray interactions in the atmosphere could also produce a long-lived stau signal in neutrino telescopes [363] (see also Ref. [361]). At energies above 10 4 GeV the decay length of weakly decaying nucleons, charged pions and kaons is much larger than their interaction length in the air. Therefore, as they propagate in the atmosphere the probability that a long-lived hadron (h) collides with a nucleon to produce SUSY particles is just σ SUSY h /σ SM h , where σ SUSY h is the cross section to produce the SUSY particles X (gluinos or squarks) and σ SM h the Standard Model cross section of the hadron with nuclei in the atmosphere (which is also the total cross section to a good approximation). However, since σ SM h is above 100 mb, it is apparent that this probability will be very small and that it would be much larger for a neutrino propagating in matter. The energy of charged particle tracks observed in neutrino detectors is determined by measuring their specific energy loss, dE/dx. For muons with energies above an energy of about 500 GeV, the energy loss rises linearly with energy. However, since the energy loss depends on the particle Lorentz boost, a high-energy stau is practically indistinguishable from a muon with reduced energy, E µ /E τ ∼ m µ /m τ . A smoking-gun signal for pair-produced stau NLSPs are parallel tracks in the Cherenkov detector [357].
The detection efficiency of stau pairs, i.e., coincident parallel tracks, depends on the energies and directions of the staus, as in the case of single events, but also on their separation. A large fraction of staus reaching the detector will be accompanied by their stau partner from the same interaction. However, not all of the stau pairs might be seen as separable tracks in a Cherenkov telescope if they emerge from interactions too close to or also too far from the detector. The opening angle θ τ τ between staus can be estimated by the initial opening angle between the SUSY particles in deep inelastic scattering. The separation of stau tracks in the detector is then given as x 2 ∆ tan(θ τ τ /2), where ∆ is the distance of the interaction to the detector center.
Double tracks with low track separation are difficult to distinguish from single muons copiously produced either directly in air showers or from atmospheric neutrinos. A required minimum track separation in IceCube of 150 m was found to be necessary to suppress these muons, due to the geometry of the detector [364]. A possible background to the stau pair signal produced in neutrino interactions consists of parallel muon pair events from random coincidences produced by upgoing atmospheric neutrinos. However, this is bounded from above by the number of muons arriving within a coincidence time-window requiring them to be almost parallel and is several orders of magnitude below the stau pair event rate with N µ+µ O(10 −12 )N µ (Ref. [365]). Muon pairs from charged current muon-neutrino interactions involving final state hadrons that promptly decay into a second muon are expected to be more likely. This has been estimated in Ref. [360] for the production and decay of charmed hadrons.
The rate of stau pairs from neutrino production is largely uncertain and depends on the SUSY mass spectrum. If the SUSY mass spectrum close to observational bounds (see, e.g., Ref. [41]), the rate might reach a few events per decade in cubic-kilometer Cherenkov telescopes [357,359,360,361].
Experimentally, the situation over the years has become ever more challenging. The SUSY models typically studied to be in reach of colliders, today have lower limits of their mass scales of the order of TeV from the nonobservation at the LHC [41]. This leads to significantly reduced expected opening angles (respectively, track separations) of the double tracks and low production cross sections in these models. However the possibility exists that a model beyond the current Standard Model is realised in nature with properties that escape observation at current collider detectors. Hence, the flux of double tracks has been probed generically with IceCube [364]. It turns out, a flux of hundreds of double tracks with a separation of above 150 m may pass undetected (Fig. 38). This is largely due to the challenge of requiring two very faint tracks to be reconstructed with high accuracy while the detector is optimised for higher light yields and single tracks or cascades. Due to generally higher noise rates, a similar analysis in sea water neutrino telescopes may be even more challenging. Being the first such exploration, the potential for improvements is manifold including the implementation of a dedicated double track trigger, IceCube detector upgrades that allow identification of smaller track separations and the analysis of all years of available data.

Fractional Electric Charges
The Standard Model intrinsically does not constrain the elementary charge but observationally it appears as a physical constant, i.e., all observed colour-singlet particles have integer multiples of elementary charge. As outlined in section 8 magnetic monopoles would provide a mechanism to explain electrical charge quantisation. In GUT theories the resulting quantisation is driven by the minimum possible magnetic charge rather than directly. Free fractionally charged states hence are often predicted in multiples of 1/2 (e.g., in the Pati-Salam model [366]) or 1/3, but other and smaller fractions are possible. Beyond this, fractional charges may exist in composite objects with large (10 12 GeV) confinement scales, probably also contributing to dark matter [367].
Experimentally, fractionally charged particles in cosmic rays can be observed though the anomalously low energy loss and lower light emission. Discarding one event from the signal region, the lowest upper limits on the flux of such particles have been posed by a study of the MACRO experiment [368] assuming a simple chargesquared scaling of the energy loss. Due to its size, it is expected that IceCube can reach considerably improved sensitivity to the signature of fractionally charged particles in future analyses.

Summary
IceCube opened a new window to study the non-thermal universe in 2013 through the discovery of a high-energy neutrino flux of astrophysical origin, and a first identification of a high-energy neutrino point source may be possible through joint multi-messenger observations. While these novel results can be taken as the beginning of neutrino astronomy, the potential exists to use IceCube to probe physics topics beyond astrophysics in particle physics, not least due to its sheer size.
The observation of secondary particles produced in interactions of neutrinos or cosmic rays with matter provides a probe of Standard Model interactions at energies only marginally covered or inaccessible by particle accelerator experiments. The continuous flux of atmospheric neutrinos allows to study standard neutrino flavour oscillations at a precision that is compatible with those of dedicated oscillation experiments. Moreover, the very long oscillation baselines (thousands of kilometres for atmospheric neutrinos or giga-parsec in astrophysical neutrinos) and the very high neutrino energies (up to PeV) can probe feeble deviations from the standard three-flavour oscillation scenario, that are otherwise undetectable. Ice-Cube has indeed provided strict limits on the allowed parameter space for an additional light "sterile" neutrino state with no Standard Model interactions, reducing considerably the range of allowed values of the new oscillation variables sin 2 2θ 24 and ∆m 2 41 . Similar analyses can set limits on the degree of Lorentz invariance violation, an effect that can be factorised in terms of operators proportional to powers of the neutrino energy. The energy reach of IceCube allows to probe higher dimension operators than previous experiments.
Neutrinos are also valuable indirect messengers of dark matter annihilation and decay in the Earth, Sun, Milky Way, local galaxies, or galaxy clusters. In general, neutrino emission does not suffer from large astrophysical fore-and backgrounds like electromagnetic emission and does not suffer from deflections in magnetic fields like cosmic rays. Neither indirect dark matter detection with neutrinos shares the same systematic uncertainties from astrophysical or particle physics inputs with other search techniques. In this way, indirect limits on dark matter properties from neutrino observations are complementary to indirect searches with other messengers or direct searches with accelerator or scattering experiments. Furthermore, neutrinos can be visible from very distant dark matter sources like galaxy clusters and also probe the interior of compact sources like our Sun. Besides probing dark mater capture or self annihilation cross sections, IceCube has current world leading limits on the dark matter lifetime, extending the range of the dark matter masses probed up to two orders of magnitude with respect to results from Cherenkov telescopes.
Neutrino telescopes are also sensitive to a variety of exotic signatures produced by rare particles, like Big Bang relics, passing through the detector and emitting direct or indirect Cherenkov light, as well as luminescence. Probably the most interesting signal consists of relic monopoles that could be deciphered from atmospheric and astrophysical backgrounds as extremely bright tracks and/or anomalously slow particles. Other heavy exotic particles that could be visible in this are relic Q-balls and strangelets. Cherenkov emission of long-lived supersymmetric particles or fractionally charged particles can also be considered. These are just a few examples of the many possibilities how neutrino observatories can be uses as multi-purpose particle detectors.
In this review we have summarised the many possibilities how the IceCube Observatory can probe fundamental questions of particle physics. Proposed future extensions of IceCube will enhance the sensitivity of these searches [369]. A low-energy in-fill, such as PINGU [370], would provide highly competitive measurements of the atmospheric neutrino oscillation parameters, the neutrino mass ordering, or the rate of tau neutrino appearance. It would also be more sensitive to indirect signals of lowmass dark matter. On the other hand, high-energy extensions would allow to study the astrophysical flux of neutrinos with better precision and over a wider energy range. This would reduce systematic uncertainties regarding neutrino spectra and flavour composition and help to establish astrophysical neutrinos as a probe of neutrino interactions and oscillations over ultra-long baselines.