Pushing the Energy and Cosmic Frontiers with High-Energy Astrophysical Neutrinos
Mauricio Bustamante (Bohr Inst. & Copenhagen U.)

TL;DR
High-energy astrophysical neutrinos detected by IceCube serve as unique probes for fundamental particle physics at energy scales and distances unreachable by other methods, enabling tests of physics beyond the Standard Model.
Contribution
This paper reviews how IceCube neutrinos enable new tests of fundamental physics and demonstrates the transition from theoretical proposals to experimental reality.
Findings
Neutrinos travel up to Gpc distances, providing long-baseline probes.
Detection of TeV to PeV neutrinos opens new energy frontiers.
Stringent tests of physics beyond the Standard Model are possible.
Abstract
The astrophysical neutrinos recently discovered by the IceCube neutrino telescope have the highest detected neutrino energies --- from TeV to PeV --- and travel the longest distances --- up to a few Gpc, the size of the observable Universe. These features make them naturally attractive probes of fundamental particle-physics properties, possibly tiny in size, at energy scales unreachable by any other means. The decades before the IceCube discovery saw many proposals of particle-physics studies in this direction. Today, those proposals have become a reality, in spite of prevalent astrophysical unknowns. We showcase examples of studying fundamental neutrino physics at these scales, including some of the most stringent tests of physics beyond the Standard Model.
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
Pushing the Energy and Cosmic Frontiers with High-Energy Astrophysical Neutrinos111Based on an invited talk given at the 6th Symposium on Prospects in the Physics of Discrete Symmetries (DISCRETE 2018), Vienna, 26–30 November 2018
Mauricio Bustamante
Niels Bohr International Academy & DARK, Niels Bohr Institute, University of Copenhagen, Blegdamsvej 17, Copenhagen, Denmark [email protected]
Abstract
The astrophysical neutrinos recently discovered by the IceCube neutrino telescope have the highest detected neutrino energies — from TeV to PeV — and travel the longest distances — up to a few Gpc, the size of the observable Universe. These features make them naturally attractive probes of fundamental particle-physics properties, possibly tiny in size, at energy scales unreachable by any other means. The decades before the IceCube discovery saw many proposals of particle-physics studies in this direction. Today, those proposals have become a reality, in spite of prevalent astrophysical unknowns. We showcase examples of studying fundamental neutrino physics at these scales, including some of the most stringent tests of physics beyond the Standard Model.
1 Introduction
High-energy astroparticle physics has two complementary facets. One facet is oriented to astrophysics [1]. Its goal is to find and study the most energetic astrophysical sources and discover the long-sought origin of ultra-high-energy cosmic rays (UHECRs), cosmic rays with energies in excess of GeV. The other facet is oriented to particle physics [2]. Its goal is to probe the properties and interactions of fundamental particles at the highest energies, well beyond the reach of present-day particle accelerators. Here we focus on the particle physics facet; in particular, on the unique strengths of high-energy astrophysical neutrinos to explore it.
2 The fundamental-physics reach of high-energy astrophysical neutrinos
Probes of fundamental physics that use high-energy particles of astrophysical origin crucially complement collider tests performed at lower energies. In many proposed beyond-Standard-Model (BSM) theories, the size of the new effects grows with the energy of the interaction.
The Large Hadron Collider probes physics at a center-of-momentum energy of TeV. In comparison, UHECR interactions probe it at hundreds of TeV. However, UHECR interactions on the cosmic microwave background make the Universe opaque to UHECRs, so the most energetic of them rarely reach Earth. In addition, there are numerous unknowns related to UHECRs that complicate extracting fundamental physics from them, though it is still possible to do it [3]. High-energy astrophysical gamma rays are detected with energies of 100 TeV, and are believed to be produced with even higher energies. However, because the Universe is also opaque to PeV photons, they are degraded to GeV–TeV upon reaching Earth.
In contrast, high-energy astrophysical neutrinos — so far detected in the TeV–PeV range — reach Earth without being attenuated. There are at least three reasons why they are fitting probes of fundamental physics:
They have the highest neutrino energies: The have energies that are orders of magnitude higher than accelerator neutrinos, which reach a few hundred GeV at most. Thus, they can probe physics in the neutrino sector at , comparable to Tevatron energies. In the future, the detection of EeV neutrinos would allow to probe physics at , comparable to energies envisioned for the Future Circular Collider. 2. 2.
They have the longest baselines: Even though the sources of high-energy astrophysical neutrinos remain undiscovered, the isotropy in their arrival directions hints at an origin in extragalactic sources located at Gpc-scale distances from Earth (1 Gpc km), essentially the size of the observable Universe. During the long trip, even tiny new-physics effects could compound to observable levels by the time the neutrinos reach Earth. 3. 3.
They are weakly interacting: While the production mechanism and emission spectrum of neutrinos is unknown — though there are solid theoretical expectations [4, 5] — once neutrinos are emitted, they are not expected to interact during their propagation and until detected, due to their tiny cross sections. However, new-physics effects could change this. The absence of standard interaction effects during propagation removes one layer of uncertainty when looking for new physics.
There is one additional reason, arguably of practical importance. A significant part of the science program of neutrino telescopes, like IceCube, is related to searching for the highest-energy, non-thermal astrophysical sources [1]. The data that are recorded for that purpose can be directly repurposed for searching for fundamental physics [2].
Figure 1 shows the distance and energy scales of neutrinos from different sources. High-energy neutrinos reach the highest values on both axes. The figure also shows a representative sample of BSM models that can be probed at high energies, complementing probes at lower energies.
Numerous new-physics effects grow as , where is a model-specific coupling strengths, is the neutrino energy, is the propagation distance, and is an integer that determines the energy-dependence of the effect. For instance, for neutrino decays, , for CPT-odd violation of Lorentz-invariance, , and for CPT-even violation of Lorentz invariance, . Using astrophysical PeV neutrinos coming from sources located Gpc away, we can probe tiny couplings, i.e., , an improvement of several orders of magnitude over limits obtained using atmospheric neutrinos [6, 7].
3 High-energy astrophysical neutrinos: basics and current status
3.1 Producing neutrinos
Astrophysical neutrinos with TeV–PeV energies were discovered by IceCube only a few years ago [8, 9]. To date, their origin remains unknown. We believe that they are made in cosmic accelerators that take protons and ions up to energies of EeV and above. This may happen via collisionless shocks, wherein charged particles repeatedly cross shock fronts that move at relativistic speeds. This produces an energy spectrum of protons (and ions) with an expected shape close to , where is the proton energy.
The high-energy protons interact with either photon fields [10, 11] () or surrounding matter [12] () in the sources [13, 14, 15, 16, 17]. To be concrete, below we consider interactions, in a simplified form. When the proton energy and photon energy satisfy the kinematic condition GeV2, their interaction produces a short-lived resonance. This happens, for instance, when protons of 0.1–100 PeV interact with MeV photons, which several classes of astrophysical sources have been seen to emit. The resonance decays into high-energy pions: , two-thirds of the time, and , one-third of the time.
The final-state neutrons escape the sources and decay into protons that are detected as UHECRs. The final-state pions decay promptly: neutral pions make TeV–PeV gamma rays via , while charged pions make TeV–PeV neutrinos via , followed by . Gamma rays carry of the parent proton energy, and neutrinos carry ; thus, cosmic rays of tens of PeV make neutrinos and gamma rays of a few PeV. The resulting neutrino energy spectrum is a power law whose slope, spectral breaks, and high-energy cut-off is inherited from the parent protons and photons [15]. Neutrino production in interactions proceeds in a similar way, though in this case the neutrino spectrum is expected to be a power law with fewer features. The fact that high-energy cosmic rays, gamma rays, and neutrinos are made in the same processes is at the core of multi-messenger astronomy.
Unlike cosmic rays and gamma rays, neutrinos travel to Earth without losing energy to interactions with cosmic photon backgrounds. Further, while cosmic rays are deflected by extragalactic and Galactic magnetic fields, neutrinos point back to their sources. The bulk of the high-energy astrophysical neutrinos detected is compatible with coming from an isotropic distribution of extragalactic sources, of which several candidate source classes exist [4, 5]. These candidates have been observed to emit gamma rays or are believed to harbor the necessary conditions to accelerate UHECRs. The contribution from neutrinos produced in sources inside the Milky Way, concentrated towards the Galactic Center, is of order 10%, at most [18, 19].
3.2 Detecting neutrinos
The tiny cross sections that make neutrinos fitting probes of fundamental physics also make them challenging to detect. IceCube monitors 1 km3 of clear Antarctic ice using photomultipliers (PMTs) buried deep in it to search for the dim Cherenkov light that is created as a result of the interaction of a high-energy neutrino with the ice.
When a TeV–PeV neutrino deep-inelastic-scatters off a nucleon in the ice, the interaction can be charged-current (CC), i.e., , or neutral-current (NC), i.e., , where is either a proton or neutron, are final-state hadrons, and . The final-state charged particles initiate showers that emit Cherenkov light that is recorded by the PMTs. Anti-neutrinos undergo the same interactions, charge-conjugated.
The CC interaction of a or , and the NC interaction of a neutrino of any flavor, creates a “cascade” — a particle shower whose light profile expands approximately radially out from the interaction vertex for several hundred meters. Because the showers made by and look alike, it is difficult to tell these two flavors apart. On the other hand, the CC interaction of a creates a “track” — the final-state high-energy muon of the interaction propagates for up to several kilometers, leaving an elongated, clearly distinguishable track of Cherenkov light.
When looking for an astrophysical neutrino flux, the main background is from atmospheric neutrinos and muons made when cosmic rays interact in the upper atmosphere of Earth. The energy spectrum of atmospheric neutrinos follows that of the cosmic rays: it is roughly . Beyond 100 TeV, the harder spectrum of astrophysical neutrinos becomes dominant.
After about 7.5 years of running, IceCube has detected 103 showers and tracks that start inside the instrumented volume, 60 of which have energies above 60 TeV [20]. These contained events constitute evidence of astrophysical neutrinos with a statistical significance of . A separate data set of about 1000 astrophysical through-going tracks, that are born outside the detector and cross part of it, provides independent evidence with a similar statistical significance [21]. The energy spectrum of astrophysical neutrinos is either [22] or [23], depending on what data set is used.
4 Neutrino observables to probe fundamental physics
New-physics effects may affect one or more of the following high-energy neutrino observables:
- •
Energy: The neutrino energy spectrum is expected to be a relatively featureless power law. However, BSM effects can introduce spectral features that are not expected from astrophysical production mechanisms, e.g., dips and peaks due to new resonant interactions between high-energy neutrinos and low-energy relic neutrinos [24, 25, 26], or to interactions between neutrinos and dark matter [27].
- •
Arrival direction: The bulk of the diffuse neutrino flux is expected to come from isotropically distributed extragalactic sources, and so it is expected to be isotropic itself. Thus, anisotropies in the distribution of arrival directions towards large accumulations of matter — the Galactic Center or dwarf satellite galaxies of the Milky Way — could be indicative of BSM neutrino-matter or neutrino-dark matter interactions.
- •
Flavor composition: The decay of pions produces neutrinos with flavor ratios . En route to Earth, the standard expectation is that neutrino oscillations transform the flavor composition into nearly [28]. Currently, measurements of flavor composition at IceCube are compatible with this expectation [29, 23, 30]. Yet, the measurements have large uncertainties, due to how challenging it is to measure flavor [29, 23]. Numerous BSM models acting during neutrino propagation may introduce large deviations in the flavor composition [31, 32, 33, 34, 35, 36, 37, 26, 38, 39, 40]. Representative possibilities are Lorentz- and CPT-invariance violation, neutrino decay, and new long-range neutrino interactions with matter, which we expand upon below.
- •
Timing: Modifications of the energy-momentum relation introduced by Lorentz-invariance violation cause neutrinos and photons of different energies, produced simultaneously in transient astrophysical phenomena, to have different speeds and arrive on Earth at different times [41, 42, 43, 44]. Similar delays could be induced by BSM scattering of high-energy neutrinos off the low-energy relic neutrino background [45].
There are significant astrophysical and cosmological uncertainties associated to high-energy astrophysical neutrinos. Regardless, it is possible, already today, to extract fundamental physics from the above neutrino observables. We explore this below by way of a few examples.
5 Current status, open questions
Figure 3 shows a possible classification of new-physics models that are testable using high-energy astrophysical neutrinos [46]. Models are classified according to what neutrino observable from Section 4 they affect, and on whether they affect neutrinos during their production, propagation, or detection at Earth. The list of models in Fig. 3 is not exhaustive, but representative.
Below, we briefly expand on a few selected fundamental-physics measurements; see Ref. [47] for a more complete presentation.
- •
Neutrino-nucleon cross section: Until recently, the neutrino-nucleon cross section had been measured only up to neutrino energies of 350 GeV [48]. The detection of high-energy neutrinos in IceCube has allowed us to measure the cross section at energies in the TeV–PeV range [49, 50, 51]. Figure 3 shows that the measurements are compatible with precise SM predictions [52]. Yet, because the measurement uncertainties are still large, there is still room for small contributions from new physics. At higher neutrino energies, still undetected, we could probe larger deviations from the SM predictions due to new physics — e.g., large extra dimensions [53] or electroweak sphalerons [54] — and the internal structure of nucleons — e.g., color glass condensates [55].
- •
Testing fundamental symmetries: The SM is an effective theory. At energies above its energy scale, its symmetries may break down and new ones may appear. Breaking CPT or Lorentz invariance [56], two of the linchpins of the SM, would affect multiple neutrino observables [57, 58, 59, 60, 34, 61, 62, 63, 64, 65, 66, 67], increasing in strength with neutrino energy. So far, tests of the violation of Lorentz invariance have been performed with intermediate-energy [6, 7] and high-energy [68] atmospheric neutrinos. Using astrophysical neutrinos of higher energies will increase the sensitivity further.
- •
Neutrino decay: Neutrino decay channels in the SM have associated lifetimes that are at least thirty orders of magnitude longer than the age of the Universe; for all practical purposes, in the SM neutrinos are stable [69, 70, 71]. However, BSM models may introduce new, faster decay modes; e.g., via the emission of a Majoron [72, 73]. Regardless of the specific BSM decay channel, the decay of the heavier neutrinos into the lighter ones, acting across cosmological distances, could leave sizable imprints in the high-energy neutrino energy spectrum and flavor composition [74, 75, 26, 76, 77]. The sensitivity of IceCube to these imprints yields neutrino lifetime constraints that are several orders of magnitude better than the ones obtained using neutrinos coming from closer sources [76].
- •
Dark matter: The decay or self-annihilation of dark matter into neutrinos [78, 79, 80, 81], or the interaction of high-energy neutrinos with dark matter in the Galactic Center [82], could induce line-like features and low-energy pile-ups in the neutrino energy spectrum and an excess or deficit of events in the direction of large Galactic [83, 84] and extragalactic [85] concentrations of dark matter. The absence of these features in the data has placed significant constraints on the properties of dark matter, from light to superheavy [86, 27].
- •
Hidden neutrino interactions: BSM mediators may introduce new, “secret” interactions of high-energy astrophysical neutrinos with the relic neutrino background [24, 25, 26, 87], with dark backgrounds [88, 89, 90], or with cosmological distributions of ordinary matter [91] that may significantly affect the neutrino energy spectrum and flavor composition.
- •
Sterile neutrinos: The mixing of high-energy active and sterile neutrinos during propagation over cosmological distances may induce deviations from the expectation of equi-flavor composition at Earth. However, the stringent limits on active-sterile mixing angles from oscillation experiments constrain the deviations to be small, except in the unlikely scenario where the sources emit a significant amount of sterile neutrinos [92].
6 Final remarks, and the future
In the TeV–PeV range, IceCube, its planned extension IceCube-Gen2 [93], and the planned KM3NeT [94] should provide a continuous stream of data in coming years. The larger event statistics will improve measurements of neutrino cross sections and flavor composition. This will be complemented by improvements in neutrino energy, direction, and flavor reconstruction.
In the EeV range, neutrinos are predicted [95, 96], but remain undetected [97, 98, 99]. There is a guaranteed flux of “cosmogenic” EeV-scale neutrinos expected from the photohadronic interaction of UHECRs with the cosmic microwave background. However, due to prevailing UHECR unknowns — their mass composition and injection spectrum, and the redshift evolution of the number density of their sources — the predictions of the cosmogenic neutrino flux shape and normalization are uncertain [100, 101, 102, 103, 104]. Yet, if the EeV neutrino flux is high enough, a suite of envisioned next-generation ultra-high-energy neutrino telescopes [3] could probe new physics at even higher energy scales and with tinier couplings, of EeV1-n.
Astroparticle physics lies at the root of the history of particle physics: the positron, muon, and pion were discovered in cosmic rays. During the XX century, progress in particle physics was mainly driven by a progression of increasingly energetic and luminous particle colliders. It is perhaps fitting that, at a time when colliders are at an impasse in finding how to extend the Standard Model, we seek again guidance in astroparticle physics.
\ack
MB is supported by the Villum Fonden project no. 13164.
References
- [1]
Ackermann M et al. 2019 (Preprint 1903.04334)
- [2]
Ackermann M et al. 2019 (Preprint 1903.04333)
- [3]
Alves Batista R et al. 2019 (Preprint 1903.06714)
- [4]
Anchordoqui L A et al. 2014 JHEAp 1-2 1 (Preprint 1312.6587)
- [5]
Ahlers M and Halzen F 2018 Prog. Part. Nucl. Phys. 102 73–88 (Preprint 1805.11112)
- [6]
Abbasi R et al. (IceCube) 2010 Phys. Rev. D 82 112003 (Preprint 1010.4096)
- [7]
Abe K et al. (Super-Kamiokande) 2015 Phys. Rev. D 91 052003 (Preprint 1410.4267)
- [8]
Aartsen M G et al. (IceCube) 2013 Phys. Rev. Lett. 111 021103 (Preprint 1304.5356)
- [9]
Aartsen M G et al. (IceCube) 2013 Science 342 1242856 (Preprint 1311.5238)
- [10]
Mücke A, Engel R, Rachen J P, Protheroe R J and Stanev T 2000 Comput. Phys. Commun. 124 290 (Preprint astro-ph/9903478)
- [11]
Hummer S, Ruger M, Spanier F and Winter W 2010 Astrophys. J. 721 630 (Preprint 1002.1310)
- [12]
Kelner S R, Aharonian F A and Bugayov V V 2006 Phys. Rev. D 74 034018 [Erratum: Phys. Rev. D 79, 039901 (2009)] (Preprint astro-ph/0606058)
- [13]
Berezinsky V and Zatsepin G 1969 Phys. Lett. B 28 423
- [14]
Berezinsky V S and Smirnov A Yu 1975 Astrophys. Space Sci. 32 461
- [15]
Stecker F W 1979 Astrophys. J. 228 919
- [16]
Hill C T and Schramm D N 1983 Phys. Lett. B 131 247 [, 495 (1983)]
- [17]
Yoshida S and Teshima M 1993 Prog. Theor. Phys. 89 833
- [18]
Ahlers M, Bai Y, Barger V and Lu R 2016 Phys. Rev. D 93 013009 (Preprint 1505.03156)
- [19]
Aartsen M G et al. (IceCube) 2017 Astrophys. J. 849 67 (Preprint 1707.03416)
- [20]
Taboada I 2018 Presented at the Neutrino 2018 Conference, Heidelberg, Germany, June 04–09, 2018 http://doi.org/10.5281/zenodo.1286918
- [21]
Aartsen M G et al. (IceCube) 2017 (Preprint 1710.01191)
- [22]
Aartsen M G et al. (IceCube) 2016 Astrophys. J. 833 3 (Preprint 1607.08006)
- [23]
Aartsen M G et al. (IceCube) 2015 Astrophys. J. 809 98 (Preprint 1507.03991)
- [24]
Ioka K and Murase K 2014 PTEP 2014 061E01 (Preprint 1404.2279)
- [25]
Ng K C Y and Beacom J F 2014 Phys. Rev. D 90 065035 [Erratum: Phys. Rev. D 90, 089904 (2014)] (Preprint 1404.2288)
- [26]
Shoemaker I M and Murase K 2016 Phys. Rev. D 93 085004 (Preprint 1512.07228)
- [27]
Bhattacharya A, Esmaili A, Palomares-Ruiz S and Sarcevic I 2019 (Preprint 1903.12623)
- [28]
Learned J G and Pakvasa S 1995 Astropart. Phys. 3 267 (Preprint hep-ph/9405296)
- [29]
Aartsen M G et al. (IceCube) 2015 Phys. Rev. Lett. 114 171102 (Preprint 1502.03376)
- [30]
Bustamante M and Ahlers M 2019 (Preprint 1901.10087)
- [31]
Beacom J F, Bell N F, Hooper D, Pakvasa S and Weiler T J 2003 Phys. Rev. D 68 093005 [Erratum: Phys. Rev. D 72, 019901 (2005)] (Preprint hep-ph/0307025)
- [32]
Pakvasa S, Rodejohann W and Weiler T J 2008 JHEP 02 005 (Preprint 0711.4517)
- [33]
Bustamante M, Gago A M and Jones Perez J 2011 JHEP 05 133 (Preprint 1012.2728)
- [34]
Bustamante M, Gago A M and Pena-Garay C 2010 JHEP 04 066 (Preprint 1001.4878)
- [35]
Mehta P and Winter W 2011 JCAP 1103 041 (Preprint 1101.2673)
- [36]
Bustamante M, Beacom J F and Winter W 2015 Phys. Rev. Lett. 115 161302 (Preprint 1506.02645)
- [37]
Argüelles C A, Katori T and Salvado J 2015 Phys. Rev. Lett. 115 161303 (Preprint 1506.02043)
- [38]
Gonzalez-Garcia M C, Maltoni M, Martinez-Soler I and Song N 2016 Astropart. Phys. 84 15 (Preprint 1605.08055)
- [39]
Rasmussen R W, Lechner L, Ackermann M, Kowalski M and Winter W 2017 Phys. Rev. D 96 083018 (Preprint 1707.07684)
- [40]
Ahlers M, Bustamante M and Mu S 2018 Phys. Rev. D 98 123023 (Preprint 1810.00893)
- [41]
Longo M J 1987 Phys. Rev. D 36 3276
- [42]
Wang Z Y, Liu R Y and Wang X Y 2016 Phys. Rev. Lett. 116 151101 (Preprint 1602.06805)
- [43]
Diaz J S 2016 Symmetry 8 105 (Preprint 1609.09474)
- [44]
Boran S, Desai S and Kahya E O 2018 (Preprint 1807.05201)
- [45]
Murase K and Shoemaker I M 2019 (Preprint 1903.08607)
- [46]
Argüelles C, Bustamante M, Conrad J, Kheirandish A, Palomares-Ruiz S, Salvadó J and Vincent A In preparation
- [47]
Ahlers M, Helbing K and Pérez de los Heros C 2018 Eur. Phys. J. C 78 924 (Preprint 1806.05696)
- [48]
Tzanov M et al. (NuTeV) 2006 Phys. Rev. D 74 012008 (Preprint hep-ex/0509010)
- [49]
Aartsen M G et al. (IceCube) 2017 Nature 551 596 (Preprint 1711.08119)
- [50]
Bustamante M and Connolly A 2019 Phys. Rev. Lett. 122 041101 (Preprint 1711.11043)
- [51]
Aartsen M G et al. (IceCube) 2019 Phys. Rev. D 99 032004 (Preprint 1808.07629)
- [52]
Cooper-Sarkar A, Mertsch P and Sarkar S 2011 JHEP 08 042 (Preprint 1106.3723)
- [53]
Alvarez-Muniz J, Halzen F, Han T and Hooper D 2002 Phys. Rev. Lett. 88 021301 (Preprint hep-ph/0107057)
- [54]
Ellis J, Sakurai K and Spannowsky M 2016 JHEP 05 085 (Preprint 1603.06573)
- [55]
Henley E M and Jalilian-Marian J 2006 Phys. Rev. D 73 094004 (Preprint hep-ph/0512220)
- [56]
Colladay D and Kostelecky V A 1998 Phys. Rev. D 58 116002 (Preprint hep-ph/9809521)
- [57]
Kostelecky V A and Mewes M 2004 Phys. Rev. D 69 016005 (Preprint hep-ph/0309025)
- [58]
Hooper D, Morgan D and Winstanley E 2005 Phys. Rev. D 72 065009 (Preprint hep-ph/0506091)
- [59]
Kostelecky V A and Russell N 2011 Rev. Mod. Phys. 83 11 (Preprint 0801.0287)
- [60]
Bazo J L, Bustamante M, Gago A M and Miranda O G 2009 Int. J. Mod. Phys. A 24 5819 (Preprint 0907.1979)
- [61]
Kostelecky A and Mewes M 2012 Phys. Rev. D 85 096005 (Preprint 1112.6395)
- [62]
Gorham P W et al. 2012 Phys. Rev. D 86 103006 (Preprint 1207.6425)
- [63]
Borriello E, Chakraborty S, Mirizzi A and Serpico P D 2013 Phys. Rev. D 87 116009 (Preprint 1303.5843)
- [64]
Stecker F W and Scully S T 2014 Phys. Rev. D 90 043012 (Preprint 1404.7025)
- [65]
Tomar G, Mohanty S and Pakvasa S 2015 JHEP 11 022 (Preprint 1507.03193)
- [66]
Amelino-Camelia G, Guetta D and Piran T 2015 Astrophys. J. 806 269
- [67]
Liao J and Marfatia D 2018 Phys. Rev. D 97 041302 (Preprint 1711.09266)
- [68]
Aartsen M G et al. (IceCube) 2018 Nature Phys. 14 961 (Preprint 1709.03434)
- [69]
Pal P B and Wolfenstein L 1982 Phys. Rev. D 25 766
- [70]
Hosotani Y 1981 Nucl. Phys. B 191 411 [Erratum: Nucl. Phys. B 197, 546 (1982)]
- [71]
Nieves J F, Pal P B and Unger D G 1983 Phys. Rev. D 28 908
- [72]
Chikashige Y, Mohapatra R N and Peccei R D 1980 Phys. Rev. Lett. 45 1926 [, 921 (1980)]
- [73]
Gelmini G B, Nussinov S and Roncadelli M 1982 Nucl. Phys. B 209 157
- [74]
Beacom J F, Bell N F, Hooper D, Pakvasa S and Weiler T J 2003 Phys. Rev. Lett. 90 181301 (Preprint hep-ph/0211305)
- [75]
Baerwald P, Bustamante M and Winter W 2012 JCAP 1210 020 (Preprint 1208.4600)
- [76]
Bustamante M, Beacom J F and Murase K 2017 Phys. Rev. D 95 063013 (Preprint 1610.02096)
- [77]
Denton P B and Tamborra I 2018 Phys. Rev. Lett. 121 121802 (Preprint 1805.05950)
- [78]
Beacom J F, Bell N F and Mack G D 2007 Phys. Rev. Lett. 99 231301 (Preprint astro-ph/0608090)
- [79]
Yuksel H, Horiuchi S, Beacom J F and Ando S 2007 Phys. Rev. D 76 123506 (Preprint 0707.0196)
- [80]
Feng J L 2010 Ann. Rev. Astron. Astrophys. 48 495 (Preprint 1003.0904)
- [81]
Murase K and Beacom J F 2012 JCAP 1210 043 (Preprint 1206.2595)
- [82]
Argüelles C A, Kheirandish A and Vincent A C 2017 Phys. Rev. Lett. 119 201801 (Preprint 1703.00451)
- [83]
Adrian-Martinez S et al. (ANTARES) 2015 JCAP 1510 068 (Preprint 1505.04866)
- [84]
Aartsen M G et al. (IceCube) 2017 Eur. Phys. J. C 77 627 (Preprint 1705.08103)
- [85]
Aartsen M G et al. (IceCube) 2013 Phys. Rev. D 88 122001 (Preprint 1307.3473)
- [86]
Chianese M, Miele G and Morisi S 2017 Phys. Lett. B 773 591 (Preprint 1707.05241)
- [87]
Kelly K J and Machado P A N 2018 JCAP 1810 048 (Preprint 1808.02889)
- [88]
Anchordoqui L, Barger V, Goldberg H and Marfatia D 2008 Phys. Lett. B 660 529 (Preprint 0711.4055)
- [89]
Klop N and Ando S 2018 Phys. Rev. D 97 063006 (Preprint 1712.05413)
- [90]
Capozzi F, Shoemaker I M and Vecchi L 2018 JCAP 1807 004 (Preprint 1804.05117)
- [91]
Bustamante M and Agarwalla S K 2019 Phys. Rev. Lett. 122 061103 (Preprint 1808.02042)
- [92]
Brdar V, Kopp J and Wang X P 2017 JCAP 1701 026 (Preprint 1611.04598)
- [93]
Aartsen M G et al. (IceCube) 2014 (Preprint 1412.5106)
- [94]
Adrian-Martinez S et al. (KM3Net) 2016 J. Phys. G 43 084001 (Preprint 1601.07459)
- [95]
Greisen K 1966 Phys. Rev. Lett. 16 748
- [96]
Zatsepin G T and Kuzmin V A 1966 JETP Lett. 4 78 [Pisma Zh. Eksp. Teor. Fiz. 4, 114 (1966)]
- [97]
Zas E (Pierre Auger) 2018 PoS ICRC2017 972
- [98]
Aartsen M G et al. (IceCube) 2018 Phys. Rev. D 98 062003 (Preprint 1807.01820)
- [99]
Gorham P W et al. (ANITA) 2019 (Preprint 1902.04005)
- [100]
Kotera K, Allard D and Olinto A V 2010 JCAP 1010 013 (Preprint 1009.1382)
- [101]
Ahlers M and Halzen F 2012 Phys. Rev. D 86 083010 (Preprint 1208.4181)
- [102]
Romero-Wolf A and Ave M 2018 JCAP 1807 025 (Preprint 1712.07290)
- [103]
Alves Batista R, de Almeida R M, Lago B and Kotera K 2019 JCAP 1901 002 (Preprint 1806.10879)
- [104]
Heinze J, Fedynitch A, Boncioli D and Winter W 2019 (Preprint 1901.03338)
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Ackermann M et al. 2019 ( Preprint 1903.04334 )
- 2[2] Ackermann M et al. 2019 ( Preprint 1903.04333 )
- 3[3] Alves Batista R et al. 2019 ( Preprint 1903.06714 )
- 4[4] Anchordoqui L A et al. 2014 JHE Ap 1-2 1 ( Preprint 1312.6587 )
- 5[5] Ahlers M and Halzen F 2018 Prog. Part. Nucl. Phys. 102 73–88 ( Preprint 1805.11112 )
- 6[6] Abbasi R et al. (Ice Cube) 2010 Phys. Rev. D 82 112003 ( Preprint 1010.4096 )
- 7[7] Abe K et al. (Super-Kamiokande) 2015 Phys. Rev. D 91 052003 ( Preprint 1410.4267 )
- 8[8] Aartsen M G et al. (Ice Cube) 2013 Phys. Rev. Lett. 111 021103 ( Preprint 1304.5356 )
