Acta Polytechnica doi:10.14311/AP.2013.53.0770 Acta Polytechnica 53(Supplement):770–775, 2013 © Czech Technical University in Prague, 2013 available online at http://ojs.cvut.cz/ojs/index.php/ap ICECUBE OBSERVATORY: NEUTRINOS AND THE ORIGIN OF COSMIC RAYS Paolo Desiatia,b,∗, for the IceCube Collaborationc a Wisconsin IceCube Particle Astrophysics Center (WIPAC), University of Wisconsin, Madison, WI 53706, U.S.A. b Department of Astronomy, University of Wisconsin, Madison, WI 53706, U.S.A. c http: // icecube. wisc. edu ∗ corresponding author: desiati@icecube.wisc.edu Abstract. The completed IceCube Observatory, the first km3 neutrino telescope, is already providing the most stringent limits on the flux of high energy cosmic neutrinos from point-like and diffuse galactic and extra-galactic sources. The non-detection of extra-terrestrial neutrinos has important consequences on the origin of the cosmic rays. Here the current status of astrophysical neutrino searches, and of the observation of a persistent cosmic ray anisotropy above 100 TeV, are reviewed. Keywords: neutrinos, cosmic rays, anisotropy. 1. Introduction One hundred years after their discovery, the origin of the cosmic rays is still a mystery. The current leading model is that cosmic rays are accelerated in diffusive shocks. In this case Supernova Remnants (SNRs) in our Galaxy could be the major source of cosmic rays up to about 1015 ÷1017 eV. The SNR energy output in the Galaxy can provide the energy budget necessary to maintain the presently observed population of galactic cosmic-rays. In particular, in order to achieve such high energies it is expected that acceleration occurs during the relatively short period in the SNR evolution between the end of free expansion and the beginning of the so-called Sedov phase. This period is about 103 years from the explosion when the shock velocity is high enough to allow for efficient acceleration. At energies in excess of about 1017 eV, Active Galactic Nuclei (AGN) and Gamma Ray Bursts (GRB) could play an important role in the origin of the extra- galactic cosmic rays. Since cosmic rays are deflected by magnetic fields, it is not possible to associate them to their sources. However, if hadronic particles are accelerated, a frac- tion of them would interact within their sources or in surrounding molecular clouds to produce mesons. The mesons eventually decay into high energy γ-rays and neutrinos with an energy spectrum ∼ E−2 of the accel- erated cosmic rays. The remaining hadronic particles propagate until their detection on Earth. Detection of γ-rays and neutrinos from individual galactic or extra- galactic source candidates of cosmic rays, or from extended molecular clouds, is therefore a method to indirectly probe the origin of cosmic rays. During the last decade, detection of γ-rays from galactic sources has been successfully achieved by satel- lite experiments such as AGILE and Fermi up to 10 and 100 GeV, respectively. Imaging Cherenkov Tele- scope Arrays such as MAGIC, VERITAS and H.E.S.S., and water Cherenkov detectors such as Milagro have made measurements up to O(10 TeV). High energy di- rect emission from old SNRs appears to be inconsistent with hadronic acceleration1. It is interesting, however, that delayed secondary γ-ray emissions can be pro- duced by the most energetic particles that escaped the acceleration region when they propagate through molecular clouds that surround the star forming re- gions [1]. With this mechanism, indirect evidence of hadronic acceleration is present even when SNR are several 104 years old. In fact, the detection of an extended emission of TeV γ-rays from the Galactic Center by H.E.S.S., which is attributed to cosmic rays accelerated by SNR G0.9+0.1 interacting with the surrounding clouds, might provide the first evidence of hadronic acceleration [2]. The most compelling evi- dence currently comes from low energy γ-ray emission from the regions surrounding the intermediate-age SNR W44. AGILE observations in the energy range of 50 MeV ÷ 10 GeV [3] and Fermi observations up to 100 GeV [4] show that while leptonic models fail to describe simultaneously γ and radio emissions with- out requiring too large circumstellar densities, the hadronic models are consistent with experimental con- straints from radio, optical, X and γ-rays observations. Although the γ-ray energy spectrum is consistent with a proton spectral index of 3 and a low energy cut-off of approximately 10 GeV 2, the hadronic origin of the observed emission is considered likely. The observed steep spectrum and low energy cut-off may be caused by suppression of efficient particle acceleration in the dense environment of this source [5]. Ion-neutral col- 1Most probably SNR older than several thousand years no longer efficiently accelerate cosmic rays. 2This is the reason why such a source was not observed at TeV energy. 770 http://dx.doi.org/10.14311/AP.2013.53.0770 http://ojs.cvut.cz/ojs/index.php/ap http://icecube.wisc.edu vol. 53 supplement/2013 IceCube Observatory: Neutrinos and the Origin of Cosmic Rays lisions in the weakly ionized dense gas surrounding the remnant lead to a softer spectrum as well as to damping of the plasma Alfvén waves that form the shock. The resulting poor particle confinement leads to a low energy cutoff [6]. Other than the specific properties of single objects, evidence of an instance of hadronic acceleration is a very important step towards the discovery of the origin of cosmic rays. However, this would not mean that all galactic cosmic rays are necessarily accelerated in SNR. If cosmic ray acceleration occurs predominantly on a larger scale, such as in superbubbles [7] or in the Galaxy cluster medium where particles could be accelerated to ultra-high energies [8], the search for the origin of cosmic rays should concentrate on extended sources or diffuse fluxes. While the TeV γ-ray horizon is limited within our Galaxy, because of absorption in the infrared and mi- crowave cosmic background, the GeV γ-emissions can be observed within about 100 Mpc making it possible to search for extragalactic sources of cosmic rays. On the other hand, detection of neutrinos from individual sources are an efficient and unambiguous probe for the high energy hadronic acceleration mechanism, and therefore for the sources of cosmic rays. However, the very same property that makes neutrinos an excellent cosmic messenger also makes them difficult to detect. Thus large instrumented volume of target matter is required to capture sufficient event statistics. The IceCube Neutrino Observatory (see Fig. 1), completed in December 2010, is currently the only km3 scale neutrino telescope collecting data. The observatory consists of an array of 5,160 optical sen- sors arranged along 86 cables (or strings) between 1,450 and 2,450 meters below the geographic South Pole, where the antarctic ice is particularly transpar- ent. IceCube includes a surface shower array, IceTop, and a dense instrumented core with a lower energy threshold, DeepCore. The surface array, IceTop, is 81 stations each consisting of two tanks of frozen clean water with each tank containing two optical sensors. IceTop, using events in coincidence with the deep IceCube array, provides the measurement of the spec- trum and mass composition of cosmic rays at the knee and up to about 1018 eV. The DeepCore sub-array, consisting of 6 densely instrumented strings located at the bottom-center of IceCube, lowers the obser- vatory neutrino energy threshold to about 10 GeV. DeepCore uses the surrounding IceCube instrumented volume as a veto for the background of cosmic ray induced through-going muon bundles, thus enhancing the detection of down-going neutrinos within the Deep Core volume. Veto rejection power in excess of 108 has been achieved [9]. The basic detection compo- nent of IceCube is the Digital Optical Module (DOM) which consists of a 10-inch Hamamatsu photomulti- plier tube (PMT) and its own data acquisition (DAQ) circuitry enclosed in a pressure-resistant glass sphere. The DOMs detect, digitize and timestamp the signals Figure 1. A schematic view of the IceCube Obser- vatory with the surface array IceTop and the densely instrumented DeepCore. from Cherenkov radiation photons. Their main DAQ board is connected to the central DAQ in the IceCube Laboratory at the surface, where the global trigger is determined [11]. The construction of IceCube started in 2004 and physics quality data taking commenced in 2006. With this early data the observatory is provid- ing the most stringent limits on the flux of high energy neutrinos from extra-terrestrial origin, and therefore strong constraints on the models of individual sources of cosmic rays and unidentified diffuse sources. At the same time, IceCube has accumulated a large number of cosmic ray induced neutrinos produced in the at- mosphere, making it possible to probe the combined effect of hadronic interaction models, cosmic ray spec- trum and composition on the neutrino spectrum up to a few hundred TeV [10]. In the search for high energy neutrinos, the large exposure of IceCube makes it possible to collect an unprecedented number of events in the form of bundles of high energy muons generated in the cosmic ray induced extensive air showers. Although these events represent an overwhelming background in the neutrino searches, they make it possible, for the first time, to determine the degree of anisotropy of cosmic rays from a few TeV to several PeV of particle energy. The persistence of a cosmic ray anisotropy at high energy raises the question of the responsible mechanism. The notion that cosmic ray anisotropy might be connected to the distribution of nearby and recent supernovae is intriguing, and might thus provide a new probe into the origin of the cosmic rays. On the other hand the complex energy-dependent topology suggests that non-diffusive processes in the local interstellar medium most probably play an important role. 2. Physics Results If the signals from detected Cherenkov photons satisfy specific trigger conditions, an event is defined and recorded by the surface data acquisition system. On- 771 Paolo Desiati, for the IceCube Collaboration Acta Polytechnica strings year mean µrate final νµ rate 22 2007 500 Hz 18/day 40 2008 1100 Hz 40/day 59 2009 1700 Hz 130/day 79 2010 2000 Hz 170/day 86 2011 2100 Hz 200/day Table 1. Mean rate of muon bundles and atmospheric neutrinos after final event selection for different string configurations of the IceCube Observatory (numbers in italic are predictions). line data filtering at the South Pole reduces the event volume to about 10 % of the trigger rate, based on a series of reconstruction and filter algorithms aimed to select events based on directionality, topology and energy [15]. The filter makes it possible to transfer data via satellite from the experimental site for prompt physics analyses. 2.1. Atmospheric neutrinos Of the events that trigger IceCube, the vast majority are muon bundles produced by the impact of primary cosmic rays in the atmosphere. Only a small fraction of the detected events (∼ 10−5) are muons produced by the charged current interaction of atmospheric muon neutrinos. The easiest way to reject the down- going muon bundle background is to exclusively select well reconstructed up-going events, since these can only be produced by neutrinos crossing the Earth and interacting in the matter surrounding the detector. Depending on the detector configuration and on the specific reconstruction algorithms and event selection utilized, the atmospheric neutrino sample is charac- terized by a directional resolution of better than 1° above 1 TeV. The corresponding resolution in the es- timation of the muon energy is about 0.2 ÷ 0.3 (in log10 of the energy) for crossing track-like events, and about 0.1 or better for contained cascade-like events. Typically, 30 % ÷ 40 % of the up-going events survive the selection with a background contamination of less than about 1 % (see Tab. 1). The atmospheric neutrino sample collected by Ice- Cube over the years is the largest ever recorded and currently reaches energies near 400 TeV (see Fig. 2). For the first time the precision of this measurement is providing a powerful tool to constrain the effects of high energy hadronic interaction models that represent our present knowledge of the cosmic ray induced ex- tensive air showers and the spectrum and composition of primary cosmic rays [10]. 2.2. Search for astrophysical ν’s Atmospheric neutrinos represent an irreducible back- ground for the search of high energy astrophysical neutrinos. If hadronic acceleration is the underlying process of high energy cosmic ray production and γ-ray observations in galactic and extra-galactic sources, the Figure 2. Collection of theoretical calculations and experimental measurements of the atmospheric neu- trino spectrum. Shown is the predicted conventional νµ + ν̄µ (blue line) and νe + ν̄e (red line) flux from [16], and the predicted prompt flux of neutrinos (magenta band) from [17]. The unfolded energy spectrum [18] (black filled circles) and forward folded spectrum [19] (gray band) from the 40-string IceCube configuration, unfolded spectrum [20] (black open circles) and for- ward folded spectrum [21] (ecru band) from AMANDA are presented. The results from Super-K [22] (aqua band) and that from Fréjus [23] (black filled squares for νµ + ν̄µ and black open squares for νe + ν̄e) are also presented. charged mesons could produce enough neutrinos to be observed in a detector the size of IceCube. Figure 3 shows the sensitivity (90 % CL) of IceCube for the full-sky search of steady point sources of E−2 muon neutrinos as a function of declination, along with that of other experiments. The extension of the point source search to the southern hemisphere is made possible by a high energy event selection that rejects the background down-going events by five or- ders of magnitude, and restricts neutrino energies to above 100 TeV. Still dominated by high energy large muon bundles, this makes the southern hemisphere poor in atmospheric neutrinos yielding a low neu- trino detection sensitivity. Nevertheless, this provides IceCube with a full-sky view that complements cover- age of the neutrino telescopes in the Mediterranean. The figure shows the sensitivities from IceCube and other observatories (interpreted as the median upper limit we expect to observe from individual sources across the sky) along with upper limits from selected sources. The sensitivity is reaching the level of current predictions for flux from astrophysical sources (i.e. be- low 1012 × E−2 TeV cm−2 s−1) although the discovery potential, defined to be 5σ for 50 % of the trials, is typically a factor of three higher than the sensitivity. Therefore constraints on the parameters of hadronic acceleration models are starting to develop. Searches for neutrinos from transient [30] and peri- 772 vol. 53 supplement/2013 IceCube Observatory: Neutrinos and the Origin of Cosmic Rays Figure 3. Sensitivity (90 % CL) for a full-sky search of steady point sources of muon neutrinos with an E−2 energy spectrum as a function of declination angle for IceCube and other experiments. Note that for IceCube, events with δ < 0° are down-going, coming from the southern hemisphere, and events with δ > 0° are up- going and come from the northern hemisphere. Figure 4. Upper limits (90 % CL) for neutrino searches in coincidence with Gamma Ray Bursts with 40 strings of IceCube and the combined 40- and 59- string detector configurations [32]. Also shown is the Waxman & Bahcall predicted average flux [33]. odic [31] sources have also been performed. In particu- lar, a time window scan for transient sources (with no external triggers) shows that the discovery potential drops by a factor of 2 if searching for 1 day duration flares. A particular search for transient sources is that for neutrinos from GRB. For the first time, the Ice- Cube Observatory has provided a definitive test of the GRB models with the most stringent constraints. Fig- ure 4 shows the upper limits obtained with the data collected by the 40-string configuration o IceCube and by the combined data of the 40- and 59-string config- urations [32]. For each detector configuration, a list of GRBs detected during the corresponding physics runs was compiled and the predicted neutrino flux was calculated based on the γ-ray spectrum shown in [34]. The corresponding stacked neutrino flux was used to search for events collected within the time window in which 5 % to 95 % of the fluence is recorded (i.e. T90). The upper limit is about 3 times below the predicted flux of the Waxman & Bahcall model, challenging the hypothesis that GRB are the sources of Ultra High Energy Cosmic Rays (UHECR). This re- sult has profound consequences for the predicted flux of neutrinos produced by the interaction of UHECR with the cosmic microwave background, the so-called cosmogenic neutrinos, as well as for the GeV–TeV γ-ray background flux (see for instance [35, 36]). It is important to note that it was recently shown that the fireball model with refined assumptions yields a 10 times smaller predicted flux (see [37, 38]). There is the possibility that the bulk of cosmic rays does not originate from individual sources, but from large-scale acceleration processes in superbubbles or even Galaxy clusters. In addition, unresolved sources of cosmic rays over cosmological times are expected to have produced detectable fluxes of diffuse neutrinos. Since shock acceleration is expected to provide an ∼ E−2 energy spectrum, harder than the ∼ E−3.7 of the atmospheric neutrinos, the diffuse flux is expected to dominate at high energy where the sensitivity is strongly dependent on the experimental quality of the selected events. Figure 5 shows a collection of sensitivities and up- per limits (90 % CL) for an E−2 flux of νµ + ν̄µ, from AMANDA, Antares and various IceCube configura- tions compared to the experimental and theoretical flux of the atmospheric neutrinos and various models of astrophysical neutrinos. The most recent results lie below the Waxmann & Bahcall neutrino bound [42], again indicating IceCube’s potential for discovering the origin of cosmic rays. In the Ultra High Energy range (UHE), above ∼ 106 GeV, IceCube is reaching a competitive sen- sitivity as well. At this level one begins to reach current models of cosmogenic neutrino production (see Fig. 6) that are simultaneously constrained by the current observations of UHECRs and the GeV γ-rays by Fermi-LAT [49]. Taking into account that UHECR mass composition is a key ingredient for the absolute flux and spectral shape of cosmogenic neutri- nos [35], its large uncertainty still weighs profoundly on current models. This means that although the IceCube sensitivity to UHE neutrinos is currently the best ever achieved below 1010 GeV it might be still far from the actual flux. From this point of view, the current developments toward a radio array in Antarc- tica, such as Askaryan Radio Array (ARA) [51] is a natural extension toward the highest energies. It is worth noting that the preliminary sensitivity for an arbitrary spectrum, shown in Fig. 6, has a minimum just above 1 PeV, where no significant cos- mogenic neutrino flux is expected. In the experimental analysis performed on data collected during 2010–12, where events with a large number of detected photons were selected, two events were found on a background of conventional atmospheric neutrinos of 0.3. The events deposited an energy in the detector of about 1 PeV, and further study is underway to determine their nature. One possible hypothesis is that these 773 Paolo Desiati, for the IceCube Collaboration Acta Polytechnica [GeV]) ν log10(E 3 4 5 6 7 8 ] -1 s r -1 s -2 [ G e V c m ν /d E Φ d 2 ν E -9 10 -8 10 -710 -6 10 -5 10 2000-2003 90%CL limitµνAMANDA 07-09 90%CL limitµνANTARES 90%CL limitµνIC40 IC59 diffuse sensitivity IC59 diffuse 90%CL limit IC40 atmospheric unfolding (HKKM07)µνconventional atmospheric µνconventional (HKKM07) + prompt (Enberg et al.) Waxman-Bahcall upper bound (2011) Mannheim 1995 BBR I 2005 steep spectra sources Stecker AGN (Seyfert) 2005 High Peaked BL Lac (max) Mucke 2003 Prompt GRB Razzaque et al. 2008 Preliminary Figure 5. Experimental upper limits (90 % CL) for the diffuse muon neutrino flux (including the pre- liminary result from the 59-string configuration of IceCube) along with atmospheric neutrino observa- tions and theoretical models of atmospheric and extra- terrestrial neutrino fluxes. From top to bottom in the legend [16, 18, 39–47]. events represent an upper fluctuation of the prompt neutrino production in the atmosphere from the decay of heavy charm mesons. 2.3. Cosmic ray anisotropy The large number of muon bundle events collected by IceCube (about 1010–1011 each year, depending on the detector configuration) makes it possible to study the arrival direction distribution of the cosmic rays at a level of about 10−5. The bundles of highly collimated atmospheric muons share the same direction as the parent cosmic ray particle. Since this study does not require highly well reconstructed muon directions, all collected and reconstructed events with a median angular resolution of about 3° are used. Using full simulation of cosmic ray induced extensive air shower we find that the median particle energy of the IceCube data sample is about 20 TeV. With these data IceCube provides the first high statistics determination of the anisotropy of galactic cosmic rays in the southern hemisphere in the multi-TeV energy range. The large scale anisotropy observed by IceCube [52] appears to complement the observations in the north- ern hemisphere, providing for the first time an all-sky view of TeV cosmic ray arrival directions. The sky map obtained by subtracting an averaged map (over a scale of 30°–60°) from the data [53], shows signifi- cant small angular scale structures in the cosmic ray anisotropy, similarly to observations in the northern hemisphere [54, 55]. Figure 6. Preliminary sensitivity (90 % CL) for the detection of UHE neutrinos, compared to other ex- perimental results and to predictions [48–50]. The sensitivity curves are evaluated at each decade of en- ergy. Another interesting result obtained by IceCube is the persistence of the anisotropy at an energy in ex- cess of 100 TeV. At such energies a different structure is observed that can be interpreted in terms of a differ- ent phase [56] as already reported by the EAS-TOP shower array in the northern hemisphere for the first time [58]. The observation at high energy was re- cently confirmed by the preliminary result from the IceTop shower array [57]. The change of the anisotropy pattern at about 100 TeV may suggest that the helio- sphere could have an effect in flipping the apparent direction of the anisotropy. In fact, at about 100 TeV the cosmic rays’ gyro-radius in the 3 µG local inter- stellar magnetic field is of the order of magnitude of the elongated heliosphere. Below this energy scale the scattering processes on the heliospheric perturba- tions at the boundary with the interstellar magnetic field might be the dominant processes affecting the global cosmic ray arrival distribution and the small angular structure as well (see [59] where a review of other proposed models is also given). The Milagro observation of a likely harder than average cosmic ray spectrum from the localized excess region toward the direction of the heliotail, the so-called region B in [54] and also observed by ARGO-YBJ shower ar- ray [55], have triggered astrophysics interpretations (see [60–62]). However, this may suggest that some type of re-acceleration mechanism associated with cos- mic ray propagation in the turbulent heliospheric tail might occur [63, 64]. On the other hand, the TeV cosmic ray anisotropy is a tracer of the local interstel- lar magnetic field, and it might indicate cosmic ray 774 vol. 53 supplement/2013 IceCube Observatory: Neutrinos and the Origin of Cosmic Rays streaming along the magnetic field lines due to the Loop I shell expanding from the Scorpion-Centaurus Association [65]. If the local propagation effects on the cosmic ray anisotropy below 100 TeV are dominant, at higher energy it is reasonable to believe that the persistent anisotropy might be a natural consequence of the stochastic nature of cosmic ray galactic sources, in particular nearby and recent SNRs [66–68]. References [1] Gabici S., Aharonian F.A.: 2007, ApJ 665, L131 [2] Aharonian F.A., et al.: 2006, Nature 439 (7077), 695 [3] Giuliani A., et al.: 2011, arXiv:1111.4868 [4] Uchiyama Y., et al.: 2012, arXiv:1203.3234 [5] Uchiyama Y., et al.: 2010, ApJ 723, L122 [6] Malkov M.A., et al.: 2011, arXiv:1004.4714 [7] Butt, Y. :2009, Nature 460, 701 [8] Lazarian A. and Brunetti G.: 2011, arXiv:1108.2268 [9] Hyon Ha, C., et al.: 2012, J. Phys.: Conf. Ser. 375, 052034 [10] Fedynitch A., et al.: 2012, submitted to PRD, arXiv:1206.6710 [11] Abbasi R., et al.: 2009, NIM A 601, 294 [12] Abbasi R., et al.: 2010 NIM A 618, 139 [13] Andrés E., et al.: 2006 J. Geophys. Res 111, D13203 [14] Abbasi R., et al.: in preparation [15] Ahrens, J., et al.: 2004, NIM A 524, 169 [16] Honda M., et al.: 2007, Phys. Rev. D75, 043006 [17] Enberg R., et al.:2008, Phys. Rev. D78, 043005 [18] Abbasi R., et al.: 2011, Phys. Rev. D83, 012001 [19] Abbasi R., et al: 2011, Phys. Rev. D84, 082001 [20] Abbasi R., et al.: 2010, Astropart. Phys. 34, 48 [21] Abbasi R., et al.: 2009, Phys. Rev. D79, 102005 [22] Gonzalez-Garcia M.C., Maltoni M, and Rojo: 2006, J. High Energy Phys. 0610, 075 [23] Daum K., et al.: 1995, Z. Phys. C66, 417 [24] Thrane E., et al.: 2009, ApJ 704, 503 [25] Adrián-Martínez S., et al.: 2011, arXiv:1108.0292 [26] Abbasi R., et al.: 2011, Astrop. Phys. 34, 420 [27] Abbasi R., et al.: 2011, ApJ 732, 18 [28] Abbasi R., et al.: 2011, in Proc. of ICRC 2011, Beijing, China, arXiv:1111.2741 [29] Kooijman P., et al.: 2011, in Proc. of 32nd ICRC, Vol. 4, p.351, Beijing, China [30] Abbasi R., et al.: 2012, ApJ 744, 1 [31] Abbasi R., et al.: 2012, ApJ 748, 118 [32] Abbasi R., et al.: 2012, Nature 484, 351 [33] Waxman, E. & Bahcall, J.N.: 1997, Phys. Rev. Lett. 78, 2292 [34] Guetta D., et al.: 2004, Astropart. Phys. 20, 429 [35] Anchordoqui L.A., et al., 2007, Phys. Rev. D76, 123008 [36] Ahlers M., et al.: 2011, Astrop. Phys. 35-2, 87 [37] Baerwald, et al.: 2011, Phys. Rev. D83, 067303 [38] Hümmer, et al.: 2011, arXiv:1112.1076 [39] Achterberg A., et al.: 2007, Phys. Rev. D76, 042008 [40] Biagi S., et al.: 2010, NEUTRINO 2010 Athens, Greece [41] Abbasi R., et al.: 2011, Phys. Rev. D84, 082001 [42] Waxman E. & Bahcall J.N.: 1998, Phys. Rev. D59, 023002 [43] Waxman E.: 2011, to be published in Astronomy at the Frontiers of Science [44] Mannheim K.: 1995, Astropart. Phys. 3, 295 [45] Muecke A., et al.: 2003, Astropart. Phys. 18, 593 [46] Becker J.K., et al.: 2005, Astropart. Phys. 23, 355 [47] Stecker F.W.: 2005, Phys. Rev. D72, 107301 [48] Kotera K., et al.: 2010, JCAP 10, 013 [49] Ahlers M., et al.: 2010, Astrop. Phys. 34, 106 [50] Yoshida S., et al.: 1997, ApJ 479, 547 [51] Allison P., et al.: 2011, arXiv:1105.2854 [52] Abbasi R., et al.: 2010, ApJ 718, L194 [53] Abbasi R., et al.: 2011, ApJ, 740, 16 [54] Abdo A.A., et al.: 2008, Phys. Rev. Lett. 101, 221101 [55] Di Sciascio G., et al.: 2012, arXiv:1202.3379 [56] Abbasi R., et al.: 2012, ApJ 746, 33 [57] Santander M., et al.: 2012, arXiv:1205.3969 [58] Aglietta M., et al.: 2009, ApJ 692, L130 [59] Desiati P. and Lazarian A.: 2012, submitted to ApJ [60] Salvati M. & Sacco B.: 2008, Astron. & Astrophys 485, 527 [61] Drury L.O’C. & Aharonian F.A.: 2008, Astropart. Phys. 29, 420 [62] Malkov M.A., et al.: 2010, Astrophys. J. 721, 750 [63] Lazarian A. and Desiati P.: 2010, ApJ 722, 188 [64] Desiati P. and Lazarian A.: 2012, Nonlinear Proc. in Geophys. 19, 351 [65] Frisch P.C., et al.: 2012, submitted to ApJ, arXiv:1206.1273 [66] Erlykin A.D. & Wolfendale A.W.: 2006, Astrop. Phys. 25, 183 [67] Blasi P. & Amato E.: 2012, JCAP 1, 11 [68] Biermann P., et al.: 2012 Discussion Peter Grieder — Do the anisotropies which you have established with IceCube muons correlate with known magnetic structures in our part of the Galaxy? Paolo Desiati — There are a few articles from Frisch et al. on the properties of the local interstellar magnetic field in relation to the evolved Loop I sub-shell expanding from the Scorpion-Centaurus Association. There seems to be some link between the direction of the magnetic field and the Loop I shell that is correlated with the TeV cosmic ray anisotropy. This is a compelling argument that I mentioned in the text of this paper as well. 775 Acta Polytechnica 53(Supplement):770–775, 2013 1 Introduction 2 Physics Results 2.1 Atmospheric neutrinos 2.2 Search for astrophysical neutrino's 2.3 Cosmic ray anisotropy References