Ultra-high Energy Cosmic Rays and Neutrinos from Gamma-Ray Bursts, Hypernovae and Galactic Shocks1footnote 11footnote 1Based on a talk given at the Origin of Cosmic Rays: Beyond the Standard Model conference in San Vito di Cadore, Dolomites, Italy, 16-22 March 2014. This is not a comprehensive review of the topics in the title; it is weighted towards work in which I have been more personally involved.

Ultra-high Energy Cosmic Rays and Neutrinos from Gamma-Ray Bursts, Hypernovae and Galactic Shocks111Based on a talk given at the Origin of Cosmic Rays: Beyond the Standard Model conference in San Vito di Cadore, Dolomites, Italy, 16-22 March 2014. This is not a comprehensive review of the topics in the title; it is weighted towards work in which I have been more personally involved.

P. Mészáros  Center for Particle and Gravitational Astrophysics, Dept. of Astronomy & Astrophysics, Dept. of Physics, 525 Davey Laboratory, Pennsylvania State University, University Park, PA 16802, U.S.A.
Abstract

I review gamma-ray burst models (GRBs) and observations, and discuss the possible production of ultra-high energy cosmic rays and neutrinos in both the standard internal shock models and the newer generation of photospheric and hadronic GRB models, in the light of current constraints imposed by IceCube, Auger and TA observations. I then discuss models that have been proposed to explain the recent astrophysical PeV neutrino observations, including star-forming and star-burst galaxies, hypernovae and galaxy accretion and merger shocks.

keywords:
Cosmic rays, Neutrinos,

1 Introduction

The origin of the cosmic rays above the knee () and up to the range of ultra-high energy cosmic rays (UHECRs, ) remains a mystery. Attempts at correlating the arrival directions of UHECRs with known AGNs have so far yielded no convincing results Auger+10corr (); Auger+13anis2 (); Abuzayad+13TAcorr (). Partly for this reason, other high energy sources, which are distributed among, or connected with, more common galaxies, have been the subject of much interest. These include gamma-ray bursts (GRBs), hypernovae (HNe) and galactic shocks, the latter being due either to accretion onto galaxies (or clusters) or galaxy mergers.

An important clue for the presumed sources of UHECR would be the detection of ultra-high energy neutrinos (UHENUs) resulting from either photohadronic () or hadronuclear () interactions of the UHECR within the host source environment and/or during propagation towards the observer. The value of this is of course that neutrinos travel essentially unabsorbed along straight lines (or geodesics) to the observer, thus pointing back at the source. Such interactions leading to neutrinos, arising via charged pions, also result in a comparable number of neutral pions leading to high energy gamma-rays, which are however more prone to subsequent degradation via cascades against low energy ambient or intergalactic photons.

The prospect of tagging UHECRs via their secondary neutrinos has recently become extremely interesting because of the announcement by IceCube IC3+13pevnu2 () of the discovery of an isotropic neutrino background (INB) at PeV and sub-PeV energies, which so far cannot be associated with any known sources, but whose spectrum is clearly well above the atmospheric neutrino background, and is almost certainly astrophysical in origin.

2 Gamma-Ray Bursts

There are at least two types of GRBs Kouveliotou+93 (), the long GRBs (LGRBs), whose -ray light curve lasts , and the short GRBs (SGRBs), whose light curve lasts . The spectra of both peak in the MeV range, with power law extensions below and above the peak of (photon number) slopes and , the peak energy of the SGRBs being generally harder (few MeV) than those of the LGRBs () Fishman+95cgro (). This broken power law spectral shape, known as a Band spectrum, is accompanied in some cases by a lower energy (tens of keV) and less prominent black-body hump, and/or by a second, higher energy power law component, in the sub-GeV to GeV range, whose photon number slope is appreciably harder then the super-MeV slope, e.g. Gehrels+12sci (); Meszaros+12raa () (Fig. 1).

The MeV light curves exhibit short timescale variability down to ms, extensively charted along with the MeV spectra by the CGRO BATSE, the Swift BAT and more recently by the Fermi GBM instruments, while the GeV light curves and spectra have in the last several years been charted by the Fermi LAT instrument, e.g. Gehrels+12sci ().

An extremely interesting property shown by most of the LAT-detected bursts is that the light curves at GeV energies start with a time lag of several seconds for LGRBs, and fractions of a second for SGRB, relative to the start of the lightcurves at MeV energies, as seen in Fig. 2.

The GeV emission amounts to about 10% and 30-50% of the total energy budget of LGRBs and SGRBs respectively, and is detected in roughly 10% of the LGRBs, and in a somewhat larger percentage of SGRBs, although the GeV detection is ubiquitous in the brighter bursts and the non-detections may be due to being below the LAT sensitivity threshold Pelassa+11-latgrbrev ().

The huge energies involved in GRB led to the view that it involves a fireball of electrons, photons and baryons which expands relativistically Cavallo+78 (); Paczynski86 (); Goodman86grb (); Shemi+90 (); Narayan+92merg (); Meszaros+92tidal (), produced by a cataclysmic stellar event. The observational and theoretical work over the past twenty years has resulted in a generally accepted view of LGRBs as originating from the core collapse of massive () stars Woosley93col (); Paczynski98grbhn (); Woosley+12coll (), whose central remnant quickly evolves to a few solar mass black hole (BH), which for a fast enough rotating core results in a brief accretion episode powering a jet which breaks through the collapsing stellar envelope. This view is observationally well supported, the LGRBs arising in star-forming regions, sometimes showing also the ejected stellar envelope as a broad-line Ic supernova, a “hypernova”, whose kinetic energy is , an order of magnitude higher than that of an ordinary SN Ic or garden variety supernova.

For SGRBs, the leading paradigm is that they arise from the merger of a compact double neutron star (NS-NS) or neutron star-black hole (NS-BH) binary Paczynski86 (); Narayan+92merg (); Meszaros+92tidal (), resulting also in an eventual central BH and a briefer accretion-fed episode resulting in a jet. Observationally this is supported by the lack of an observable supernova, and by the fact that they are observed both in star-forming and in old population galaxies, often off-set from the optical image, as expected if in the merger the remnant has been kicked off and had time to move appreciably. While the SGRB origin is less firmly established than that of LGRBs, compact mergers are nonetheless widely considered the most likely explanation, which are also of great importance as a guaranteed source of gravitational waves (GWs), being the object of scrutiny by LIGO, VIRGO and other GW detectors.

The MeV radiation providing the detector trigger as well as the slightly delayed GeV radiation are jointly called the prompt emission of the GRB. In a fraction of bursts, a prompt optical flash is also detected by ground-based robotic telescopes Akerlof+99 () or by rapidly slewing multi-wavelength GRB missions such as Swift Gehrels+09araa (). The most widely accepted view of the GRB emission is that it is produced by shocks in the relativistic outflows, the simplest example of which are the external shocks where the outflow is decelerated in the external interstellar medium or in the stellar wind of its progenitor Rees+92fball (); Meszaros+93impact (). In such shocks magnetic field can be amplified, and electrons can be Fermi accelerated into a relativistic power law energy distribution, leading to broken power law spectra peaking initially in the MeV range. Both a forward and reverse shock are expected to be present, the latter producing synchrotron radiation in the optical range, while inverse Compton (IC) radiation in the shocks also produces a GeV component Meszaros+93multi (). The fast time variability of the MeV light curves is however better explained through what is called the standard internal shock model Rees+94unsteady (), which are expected to occur in the optically thin region outside the scattering photosphere of the outflow, but inside the radius of the external shocks. The radii of the photosphere, the internal shocks and the external shocks are, respectively,

 rph ≃(LσT/4πmpc3η3)∼4×1012Lγ,52η−32.5~{}cm (1) ris ≃Γ2ctv∼3×1013η22.5tv,−2~{}cm (2) res ≃(3E0/4πnextmpc2η2)1/3 (3) ∼2×1017(E53/n0)1/2η2/32~{}cm, (4)

where are the burst total energy, luminosity, initial dimensionless entropy, coasting bulk Lorentz factor, external density and intrinsic variability timescale, e.g. Zhang+04grbrev (); Meszaros06grbrev (). If the prompt emission is due to an internal shock, the external shock can naturally result, via inverse Compton, in a delayed GeV component Meszaros+94gev ().

The above simple picture of internal and external shocks served well in the CGRO, HETE and Beppo-SAX satellite eras extending into the first half of the 2000 decade, including the discovery of X-ray and optical afterglows as well as the prompt optical emission, which were predicted by the models. It also accommodated fairly well the observed fact that the jet is collimated and when the Lorentz factor drops below the inverse of the opening angle the light curves steepen in a predictable manner.

It was however realized that simple internal shocks radiating via electron synchrotron had low radiative efficiency, and many bursts showed low energy slopes incompatible with synchrotron Preece+98death (); Ghisellini+99grbspec (). Attempts at resolving this included different radiation mechanisms, e.g. Medvedev00jitter (), which addresses the spectrum, or invoking a larger role for the scattering photosphere Meszaros+00phot (); Rees+05photdis (); Peer+06phot (), which addresses both the spectrum and efficiency issues. It is worth stressing that the need for such ”non standard” internal shocks or photospheres is important (a fact not widely recognized) when considering IceCube neutrino fluxes expected from GRBs.

The Swift satellite launched in 2004 had gamma-ray, X-ray and optical detectors, which revealed new features of the GRB afterglows, including an initial steep decay followed by a flatter decay portion of the X-ray afterglow, interspersed by X-ray spikes, finally blending into the previously known standard power law decay behavior. These features could be represented through the high latitude emission Kumar+00naked (), a continued or multi-Lorentz factor outflow, and continued internal shocks, e.g. Zhang+06ag (); Nousek+06ag ().

The Fermi satellite, launched in late 2008 and sensitive between , extended the MeV studies and opened wide the detailed study of bursts in the GeV band, which can last for and whose spectra extend in some cases up to in the source frame. The observed GeV-MeV photon delays from bursts at redshifts led to an interesting constraint on quantum gravity theories, excluding the first order term in of the usual effective field theory series expansion formulations Abdo+09-090510 (). This limit is only reinforced by the presence of additional astrophysical mechanisms for such delays.

In general, the GeV emission of all but the first few time bins is well represented by a forward shock synchrotron radiation Ghisellini+10grbrad (); Kumar+10fsb (). This holds also for the brightest GeV bursts ever discovered, GRB130427A. However, the first few time bins of the GeV emission He+11-090510 () may need to be ascribed to the prompt emission, which is also responsible for the MeV emission - for which, as mentioned, a self-consistent analysis must consider models going beyond the standard simple internal shock, c.f. below.

3 GRB UHE Neutrinos and Cosmic Rays

The pioneering works of Waxman95cr (); Waxman+97grbnu () have served as the basis for most of the thinking on UHE cosmic ray acceleration and VHE neutrino production in GRBs. These first-generation models, as one may call them, were based on a simplified “standard” internal shock (IS) model, where the bulk Lorentz factor and the variability timescale entering equ. (4) are either assumed of inferred from -ray observations, the photon spectrum is assumed to be a standard Band function and interactions occur via the -resonance. More detailed calculations of a diffuse neutrino flux, still based on this simple IS model but using specific electromagnetically (henceforth: EM) observed bursts serving to calibrate the neutrino to photon (or relativistic proton to electron, ) luminosity ratio were made by Guetta+04grbnu ().

The first IceCube data on GRBs using 40 strings and then 56 strings as the array completion progressed were presented in Ahlers+11-grbprob (); Abbasi+11-ic40nugrb (); Abbasi+12grbnu-nat (). The results using 215 EM-detected GRBs with fluxes normalized to the -ray fluxes indicated that the diffuse neutrino upper limits were a factor below this IS model predictions (Fig. 3), unless the proton to electron ratio was much less than . Both this and a model independent analysis using a broken power law photon spectrum with variable break energy and -resonance interaction indicated an inconsistency between the upper limits and a significant contribution of GRB to the UHE cosmic ray flux observed by Auger and HiRes. This was a very important first cut in constraining models with IceCube.

Subsequent investigations pointed out that the IS model fluxes used for this comparison were overestimated Li12grbnu (); Hummer+12nu-ic3 (). More careful consideration of the interaction in this model beyond the -resonance, including multi-pion and Kaon channels with the entire target photon spectrum yielded substantially lower predicted fluxes in the TeV-PeV energy range considered Hummer+12nu-ic3 (), indicating that years of observation with the full 86 string array may be needed to rule out the simple IS model (Fig. 4).

The internal shock radius depends on both the bulk Lorentz factor and the time variability of the outflow , see eq.(4). Both factors also influence the comoving magnetic field in the shock, the photon spectral peak and the photon luminosity, thus affecting the neutrino spectral flux, see Fig. (5).

Another simplification affecting the results is that the internal shocks in the above were assumed to have a constant radius, whereas they advance and expand with the flow. Calculating numerically such time-dependent IS models which accelerate CRs, including the full range of interactions and the observed -ray luminosity function and variability distributions, the current IceCube 40+59 strings upper limits are in fact compatible with GRBs contributing a significant fraction of the UHECR flux Asano+14grbcr () (Fig. 6), but the IceCube PeV neutrino flux.

More importantly, the use of the standard internal shock model, which is favored by observers for its simplicity and ease of computation, needs to be reconsidered. This model has been known for the past decade to have problems explaining the low energy -ray spectral slopes and the radiative efficiency (§2, Meszaros06grbrev ()), and alternatives free of the -ray inconsistencies have been investigated, e.g. photospheric models and modified internal shock models. The neutrino emission of baryonic photosphere models Murase08grbphotnu (); Wang+09grbphotnu () and modified IS models Murase+12reac () differs qualitatively from that of the standard IS models.

In the case of photospheric models, it is worth stressing that the spectrum is likely to deviate from a blackbody; a non-thermal of broken power law can be produced by dissipative effects, such as sub-photospheric scattering Peer+06phot (), inelastic nuclear collisions Beloborodov10pn (), photospheric shocks or magnetic dissipation Meszaros+11gevmag (), etc. The spectrum and luminosity normalization also depend on whether the dynamics of the expansion is dominated by baryonic inertia, in which case the bulk Lorentz factor initially accelerates as until it reaches the saturation value at the coasting radius ; the photospheric radius occurs generally beyond the saturation radius and is given by the first line of eq.(4).

Alternatively the dynamics might be dominated by magnetic stresses. In this case the photospheric radius depends on the value of the magnetization index , where and for extreme magnetic domination, e.g. Meszaros+11gevmag (); Gao+12photnu (). For such magnetic cases, the photospheric radius is generally in the accelerating phase , and is given by

 rph=r0η1/μT(ηT/η)1/(1+2μ) (5)

where is the launch radius () and . Fits to determine the degree of magnetic domination have been done using Fermi GBM and LAT data Veres+13fit (); Burgess+14therm (), indicating that a degree of magnetic domination does exist, which differs between bursts. A related point is that if magnetic stresses are significant in a GRB jet, this reduces the comoving photon density in the jet, allowing heavy nuclei in the jet to survive photo-dissociation Horiuchi+12nucjet (), a point of interest in view of the Auger Auger+11anis-comp () data pointing towards a heavy composition of UHECR at high energies.

The diffuse neutrino flux from both baryonic and extreme magnetic photospheric models has been computed (see Fig. 5, Gao+12photnu ()), where the extreme magnetic photospheric model is shown as red dashed lines and the baryonic photospheric model as blue dotted lines. They appear compliant with the IceCube 40+59 string upper limits, which however were calculated for a canonical Band spectral shapes, and a more spectral-specific comparison is necessary. This has been done for baryonic photospheres IC3-ICRC13nuphot (), see Fig. 7. As the observations accumulate, these constraints are getting tighter, at least for the simple IS and the simple baryonic photosphere models.

Concerning the IceCube non-detection of the EM brightest burst ever observed, GRB 130427A, a detailed calculation Gao+13-130427nu () shows that for this particular burst, except for the extreme magnetic photosphere model, the standard IS model, the baryonic photosphere model and a model-independent analysis are compatible with its non-detection.

4 Hypernovae, SFGs/SBGs, Galactic/Cluster Shocks, AGNs as UHECR/UHENU Sources

Hypernovae (henceforth HNe) are Type Ic core collapse supernovae with unusually broad lines, denoting a much higher ejecta velocity component than in usual SNeIc. This indicates a component of the ejecta reaching up to semi-relativistic velocities, with , and a corresponding inferred ejecta kinetic energy , one order of magnitude higher than that of normal SNeIc and normal SNe in general Woosley+06araa (); Nomoto+10hn-nar (). Their rate may be 1%-5% of the normal SNIc rate, i.e. as much as 500 times as frequent as GRBs Guetta+07grbhnrate (). While core collapse (collapsar) type long GRBs appear to be accompanied by HNe, the majority of HNe appear not to have a detected GRB, e.g. Soderberg+10hn (). The semi-relativistic velocity component may be due to an accretion powered jet forming in the core collapse, as in GRBs, which only for longer accretion episodes is able to break through the collapsing envelopes, while for shorter accretion episodes it is unable to break out. In both cases the jet accelerates the envelope along the jet axis more forcefully (jet-driven supernova), causing an anisotropic expansion, e.g. Campana+06-060218 (), whereas in the majority of core collapses a slow core rotation or short accretion times lead to no jet or only weak jets and a “normal” quasi-spherical SNeIc Lazzati+12hn ().

The dominant fraction of GRB-less HNe, if indeed due to a non-emerging (choked) jet, would be effectively a failed GRB, which could be detected via a neutrino signal produced in the choked, non-exiting jet, or a neutrino precursor in those collapses where the jet did emerge to produce a successful GRB Meszaros+01choked (); Horiuchi+08choked (); Murase+13choked (). Searches with IceCube have so far not found them, e.g. Taboada11choked (); Daughhetee12choked ().

An interesting aspect of HNe is that the higher bulk Lorentz factor of the ejecta leads to estimates of the maximum UHECR energy accelerated which, unlike for normal SNe, is now in the GZK range,

 εmax≃ZeBRβ=4×1018Z~{}eV (6)

especially if heavy nuclei (e.g. Fe, ) are accelerated Wang+07crhn (); Budnik+08hn (); Wang+08crnuc (). The photon field is dilute enough so that the heavy nuclei avoid photo-dissociation Wang+08crnuc (); Horiuchi+12nucjet (). The HNe kinetic energy and occurrence rate is sufficient then to explain the observed UHECR diffuse flux at GZK energies, without appearing to violate the IceCube upper limits. The HNe, as other core collapse SNe and long GRBs, occur in early type galaxies, with a larger rate in star-forming galaxies (SFGs) and even larger rate in star-burst galaxies (SBGs).

Magnetars, another type of high energy source expected from some core collapse supernovae in SFGs and SBGs, are a sub-class of fast-rotating neutron stars with an ultra-strong magnetic field, which have been considered as possible sources of UHECR and UHENU Arons03crmag (); Zhang+03magnu (); Ghisellini08corr (); Murase09magnu (); Kotera11crmag ().

SFGs make up of all galaxies, while SBGs make up . AGNs make up of all galaxies, most AGNs being radio-quiet, i.e. without an obvious jet, while radio-loud AGNs (with a prominent jet) represent of all galaxies. Radio-loud AGNs have long been considered possible UHECR and UHENU sources, e.g. Rachen+93agncr (); Berezinsky+02diffprop (). However, the lack of an angular correlation between Auger or TA UHECR events and AGNs Auger+10corr (); Abuzayad+13TAcorr () may be suggesting that more common galaxies, e.g. SFGs or SBGs, may be hosting the UHECR sources, which could be HNe, GRBs or magnetars, all of which appear capable of accelerating UHECRs at a rate sufficient to give the observed diffuse UHECR flux.

Another possibility is that UHECR are accelerated in shocks near the core of radio-quiet AGNs, where they would produce UHENUs Stecker+91nuagn (); Alvarez+04agn (); Peer+09agn (), or alternatively UHECR could be accelerated in stand-off shocks caused by the infall of intergalactic gas onto clusters of galaxies Keshet+03igshock (); Inoue+05igshock (); Murase+08clusternu (). Galactic merger shocks (GMSs) also appear capable of accelerating UHECR, with a similar energy input rate into the IGM Kashiyama+14pevmerg (); see below.

5 The PeV Neutrino Background

In 2013 the IceCube collaboration announced the discovery of the first PeV and sub-PeV neutrinos which, to a high confidence level, are of astrophysical origin IC3+13pevnu1 (); IC3+13pevnu2 (). The majority of these are cascades, whose angular resolution is , ascribed to , while a minority are Cherenkov tracks with an angular resolution due to . Their spectrum stands out above that of the diffuse atmospheric spectrum by at least , with a best fit spectrum . There is no statistically significant evidence for a concentration either towards the galactic center or the galactic plane, being compatible with an isotropic distribution. No credible correlation has been so far established with any well-defined extragalactic objects, such as AGNs, but the working assumption of an extragalactic origin is widely accepted.

A flux of PeV neutrinos from starburst galaxies at a level close to that observed level was predicted by Loeb+06nustarburst (). The actual accelerators could be hypernovae; the maximum energy of protons, from eq.(6), is sufficient for the production of PeV neutrinos Fox+13pev (), and statistically, (1) of the observed events could be due to a hypernova (or at most a few) located in the bulge of the Milky Way. However, the bulk of the observed events must come from an isotropic distribution, and hypernovae in ultra-luminous infrared galaxies (ULIRGs) or SFGs/SBGs could be responsible He+13pevhn (); Liu+13pevnuhn ().

More generally Murase+13pev (), one can ask whether hadronuclear () interactions may be responsible for this isotropic neutrino background (INB) at PeV energies, without violating the constraints imposed by the isotropic gamma-ray background (IGB) Abdo+10igbfermi () measured by Fermi. As shown by Murase+13pev (), this requires the accelerated protons to reach at least and to have an energy distribution with an index no steeper than (Fig. 8). An important point is that most events are cascades, involving electron flavor neutrinos, and the cross section is resonant at CM energies comparable to the meson mass (Glashow resonance), at around 6.3 PeV in the lab frame. Since events are not seen at this energy, but they would be expected if the proton (and neutrino) slopes where , one concludes that the proton distribution steepens or cuts off at energies PeV. Such a cutoff may be expected in scenarios where the acceleration occurs in galaxy cluster shocks or in SFG/SBGs, where this energy may correspond to that where the escape diffusion time out of the acceleration region becomes less than the injection time or the time Murase+13pev (). Broadly similar conclusions are reached by He+13pevhn (); Liu+13pevnuhn (); Chang+14pevnugam (); Anchordoqui+14pevnu ().

Suggestively, Anchordoqui+14pevnu () find a weak correlation between five known SFGs (M82, NGC253, NGC4945, SMC and IRAS18293) and the very wide, error boxes of some cascade events, but not correlation so far with any track events; they estimate that 10 years may be needed with IceCube to find track correlations with SFGs at confidence level.

Another type of large scale shocks in galaxies are the galaxy merger shocks (GMSs), which occur every time two galaxies merge. Every galaxy merged at least once during the last Hubble time, and probably more then once; in fact mergers are the way galaxies grow over cosmological time. Such galaxy mergers were considered in the PeV neutrino background context by Kashiyama+14pevmerg (). They estimate that individual major mergers involving galaxies with have an average kinetic energy of , occurring at a rate , with a relative shock velocity . For a CR acceleration fraction the UHECR energy injection rate into the Universe is (which is also the observationally inferred rate UHECR energy injection rate), with a maximum CR energy of . The interactions in the shocks and in the host galaxies lead to PeV neutrinos and GeV -rays (Fig. 9).

Individual GMSs from major mergers at would yield in IceCube on average muon events/year, or an isotropic neutrino background (INB) of of the IceCube observed PeV-sub-PeV flux. Minor mergers, whose rate is more uncertain, might contribute up to 70-100% of the INB (Fig. 9). The -ray flux from individual GMS expected is , possibly detectable by the future CTA, while the corresponding isotropic gamma background (IGB) is , about 20-30% of the observed Fermi IGB, or a somewhat larger percentage due to minor mergers Kashiyama+14pevmerg ().

6 Conclusion, Prospects

In conclusion, the sources of UHECR and the observed extragalactic UHENU are still unknown. For UHECR, an exotic physics explanation is almost certainly ruled out, mainly because any such mechanisms would produce a high energy photon component in UHECR which can be observationally ruled out, e.g. Auger+11photnu (). Anisotropy studies from Auger, which initially suggested a correlation with AGNs Auger07agncorr () have more recently, together with Telescope Array observations, yielded no significant correlation with any specific types of galaxies Auger+10corr (); Auger+13anis2 (); Abuzayad+13TAcorr (). This might favor some of the more common types of galaxies, such as possibly radio-quiet AGNs, or alternatively stellar type events such as GRBs, hypernovae or magnetars, as discussed in §§2,3,4.

The indications for a heavy UHECR composition at higher energies Auger+11anis-comp (); Auger+13icrc (); Kampert+13uhecr () would appear to disfavor AGN jets, where the composition is closer to solar, and favor evolved stellar sources, such as GRBs, hypernovae and magnetars, where a heavy composition is more natural, if the nuclei can avoid photo-dissociation (§4). These sources would also reside in more common galaxies, avoiding the anisotropy constraints.

The PeV and sub-PeV neutrinos discovered by IceCube IC3+13pevnu1 (); IC3+13pevnu2 () are an exciting development in the quest for finding the neutrino smoking gun pointing at UHECR sources, even if not at the highest energies. Standard IS GRBs appear to be ruled out as the sources for this observed diffuse neutrino flux, given the upper limits for GRBs from IceCube Abbasi+11-ic40nugrb (); Abbasi+12grbnu-nat (). Note however that these limits were obtained for simplified internal shock models, and more careful comparison needs to be made to more realistic models (see §2). Nonetheless, the normal high luminosity, electromagnetically detected GRBs, even if able to contribute to the GZK end of the UHECR distribution Asano+14grbcr (), appear inefficient as PeV neutrino sources. It is possible that low luminosity GRBs (in the electromagnetic channel) could yield appreciable PeV neutrinos Liu+13pevnugrb (), and also choked GRBs Murase+13choked () would be electromagnetically non-detected but might provide significant PeV neutrino fluxes. The fluxes, however, remain uncertain.

More attractive candidates for the PeV neutrinos are the star-forming and starburst galaxies, hosting an increased rate of hypernovae, or accretion shocks onto galaxies or clusters, or else galaxy mergers, all of which are capable of accelerating CRs up to and produce PeV neutrinos via interactions is discussed in §4.

It is also remarkable that the PeV neutrino flux is essentially at the Waxman-Bahcall (WB) bound level Waxman+97grbnu (); Bahcall+01bound () for UHECR near the GZK range, which is also comparable to the GeV range CR flux Katz+13uhecr (), Fig. 10. This suggests the intriguing prospect that the same sources may be responsible for the entire GeV-100 EeV energy range, a possibility whose testing would require much further work.

We can look forward to much further progress with continued observations from IceCube, Auger, TA and their upgrades, as well as HAWC, CTA and ground-based Cherenkov arrays and other instruments. UHECR composition and UHECR/UHENU clustering will provide important clues, as well as GeV and TeV photon observations to provide much needed additional constraints, especially if UHENU source localization is achieved..

This work was partially supported by NASA NNX 13AH50G. I am grateful to the organizers of the Origin of Cosmic Rays conference for their kind hospitality and for stimulating discussions, also held with K. Kashiyama, P. Baerwald, S, Gao and N. Senno.

References

• (1) The Pierre AUGER Collaboration, P. Abreu, M. Aglietta, E. J. Ahn, D. Allard, I. Allekotte, J. Allen, J. Alvarez Castillo, J. Alvarez-Muñiz, M. Ambrosio, et al., Update on the correlation of the highest energy cosmic rays with nearby extragalactic matter, Astroparticle Physics 34 (2010) 314–326.
• (2) Pierre Auger Collaboration, P. Abreu, M. Aglietta, M. Ahlers, E. J. Ahn, I. F. M. Albuquerque, D. Allard, I. Allekotte, J. Allen, P. Allison, et al., Constraints on the Origin of Cosmic Rays above 10 eV from Large-scale Anisotropy Searches in Data of the Pierre Auger Observatory, Astrophys.J.Lett. 762 (2013) L13.
• (3) T. Abu-Zayyad, the Telescope Array collaboration, Search for correlations of the arrival directions of ultra-high energy cosmic ray with extragalactic objects as observed by the telescope array experiment, ArXiv e-printsarXiv:1306.5808.
• (4) IceCube Collaboration, Evidence for High-Energy Extraterrestrial Neutrinos at the IceCube Detector, Science 342.
• (5) C. Kouveliotou, C. A. Meegan, G. J. Fishman, N. P. Bhat, M. S. Briggs, T. M. Koshut, W. S. Paciesas, G. N. Pendleton, Identification of two classes of gamma-ray bursts, Astrophys.J.Lett. 413 (1993) L101–L104.
• (6) G. J. Fishman, C. A. Meegan, Gamma-Ray Bursts, Annu.Rev.Astron.Astrophys. 33 (1995) 415–458.
• (7) N. Gehrels, P. Mészáros, Gamma-Ray Bursts, Science 337 (2012) 932–.
• (8) P. Mészáros, N. Gehrels, Gamma-ray bursts and their links with supernovae and cosmology, Research in Astronomy and Astrophysics 12 (2012) 1139–1161.
• (9) M. Ackermann, the Fermi collaboration, Detection of a Spectral Break in the Extra Hard Component of GRB 090926A, Astrophys.J. 729 (2011) 114–+.
• (10) A. Abdo, the Fermi Collaboration, Fermi Observations of High-Energy Gamma-Ray Emission from GRB 080916C, Science 323 (2009) 1688.
• (11) V. Pelassa, Observations of Gamma-Ray Bursts at high energies by the Fermi Large Area Telescope–two years review, in: J. E. McEnery, J. L. Racusin, & N. Gehrels (Ed.), American Institute of Physics Conference Series, Vol. 1358 of American Institute of Physics Conference Series, 2011, pp. 41–46.
• (12) G. Cavallo, M. J. Rees, A qualitative study of cosmic fireballs and gamma-ray bursts, M.N.R.A.S. 183 (1978) 359–365.
• (13) B. Paczýnski, Gamma-ray bursters at cosmological distances, Astrophys.J.Lett. 308 (1986) L43–L46.
• (14) J. Goodman, Are gamma-ray bursts optically thick?, Astrophys.J.Lett. 308 (1986) L47–L50.
• (15) A. Shemi, T. Piran, The appearance of cosmic fireballs, Astrophys.J.Lett. 365 (1990) L55–L58.
• (16) R. Narayan, B. Paczynski, T. Piran, Gamma-ray bursts as the death throes of massive binary stars, Astrophys.J.Lett. 395 (1992) L83–L86.
• (17) P. Meszaros, M. J. Rees, Tidal heating and mass loss in neutron star binaries - Implications for gamma-ray burst models, Astrophys.J. 397 (1992) 570–575.
• (18) S. E. Woosley, Gamma-ray bursts from stellar mass accretion disks around black holes, Astrophys.J. 405 (1993) 273–277.
• (19) B. Paczyński, Gamma-ray bursts as hypernovae, in: C. A. Meegan, R. D. Preece, T. M. Koshut (Eds.), Gamma-Ray Bursts, 4th Hunstville Symposium, Vol. 428 of American Institute of Physics Conference Series, 1998, pp. 783–787.
• (20) S. E. Woosley, A. Heger, Long Gamma-Ray Transients from Collapsars, Astrophys.J. 752 (2012) 32.
• (21) C. Akerlof, R. Balsano, S. Barthelmy, J. Bloch, P. Butterworth, D. Casperson, T. Cline, S. Fletcher, F. Frontera, G. Gisler, J. Heise, J. Hills, R. Kehoe, B. Lee, S. Marshall, T. McKay, R. Miller, L. Piro, W. Priedhorsky, J. Szymanski, J. Wren, Observation of contemporaneous optical radiation from a -ray burst, Nature 398 (1999) 400–402.
• (22) N. Gehrels, E. Ramirez-Ruiz, D. B. Fox, Gamma-Ray Bursts in the Swift Era, Annu.Rev.Astron.Astrophys. 47 (2009) 567–617.
• (23) M. J. Rees, P. Meszaros, Relativistic fireballs - Energy conversion and time-scales, M.N.R.A.S 258 (1992) 41P–43P.
• (24) P. Mészáros, M. J. Rees, Relativistic fireballs and their impact on external matter - Models for cosmological gamma-ray bursts, Astrophys.J. 405 (1993) 278–284.
• (25) P. Meszaros, M. J. Rees, Gamma-Ray Bursts: Multiwaveband Spectral Predictions for Blast Wave Models, Astrophys.J.Lett. 418 (1993) L59+.
• (26) M. J. Rees, P. Meszaros, Unsteady outflow models for cosmological gamma-ray bursts, Astrophys.J.Lett. 430 (1994) L93–L96.
• (27) B. Zhang, P. Mészáros, Gamma-Ray Bursts: progress, problems and prospects, International Journal of Modern Physics A 19 (2004) 2385–2472.
• (28) P. Meszaros, Gamma-ray bursts, Rept. Prog. Phys. 69 (2006) 2259–2322.
• (29) P. Mészáros, M. J. Rees, Delayed GEV Emission from Cosmological Gamma-Ray Bursts - Impact of a Relativistic Wind on External Matter, M.N.R.A.S. 269 (1994) L41+.
• (30) R. D. Preece, M. S. Briggs, R. S. Mallozzi, G. N. Pendleton, W. S. Paciesas, D. L. Band, The Synchrotron Shock Model Confronts a “Line of Death” in the BATSE Gamma-Ray Burst Data, Astrophys.J.Lett. 506 (1998) L23–L26.
• (31) G. Ghisellini, A. Celotti, Quasi-thermal comptonization and GRBs, Astron.Astrophys.Supp. 138 (1999) 527–528.
• (32) M. V. Medvedev, Theory of “Jitter” Radiation from Small-Scale Random Magnetic Fields and Prompt Emission from Gamma-Ray Burst Shocks, Astrophys.J. 540 (2000) 704–714.
• (33) P. Mészáros, M. J. Rees, Steep Slopes and Preferred Breaks in Gamma-Ray Burst Spectra: The Role of Photospheres and Comptonization, Astrophys.J. 530 (2000) 292–298.
• (34) M. J. Rees, P. Mészáros, Dissipative Photosphere Models of Gamma-Ray Bursts and X-Ray Flashes, Astrophys.J. 628 (2005) 847–852.
• (35) A. Pe’er, P. Mészáros, M. J. Rees, The Observable Effects of a Photospheric Component on GRB and XRF Prompt Emission Spectrum, Astrophys.J. 642 (2006) 995–1003.
• (36) P. Kumar, A. Panaitescu, Afterglow Emission from Naked Gamma-Ray Bursts, Astrophys.J.Lett. 541 (2000) L51–L54.
• (37) B. Zhang, Y. Z. Fan, J. Dyks, S. Kobayashi, P. Mészáros, D. N. Burrows, J. A. Nousek, N. Gehrels, Physical Processes Shaping Gamma-Ray Burst X-Ray Afterglow Light Curves: Theoretical Implications from the Swift X-Ray Telescope Observations, Astrophys.J. 642 (2006) 354–370.
• (38) J. A. Nousek, C. Kouveliotou, D. Grupe, K. L. Page, J. Granot, E. Ramirez-Ruiz, S. K. Patel, D. N. Burrows, V. Mangano, S. Barthelmy, A. P. Beardmore, S. Campana, M. Capalbi, G. Chincarini, G. Cusumano, A. D. Falcone, N. Gehrels, P. Giommi, M. R. Goad, O. Godet, C. P. Hurkett, J. A. Kennea, A. Moretti, P. T. O’Brien, J. P. Osborne, P. Romano, G. Tagliaferri, A. A. Wells, Evidence for a Canonical Gamma-Ray Burst Afterglow Light Curve in the Swift XRT Data, Astrophys.J. 642 (2006) 389–400.
• (39) A. A. Abdo, the Fermi Collaboration, A limit on the variation of the speed of light arising from quantum gravity effects, Nature 462 (2009) 331–334.
• (40) G. Ghisellini, G. Ghirlanda, L. Nava, A. Celotti, GeV emission from gamma-ray bursts: a radiative fireball?, M.N.R.A.S. 403 (2010) 926–937.
• (41) P. Kumar, R. Barniol Duran, External forward shock origin of high-energy emission for three gamma-ray bursts detected by Fermi, M.N.R.A.S. 409 (2010) 226–236.
• (42) H.-N. He, X.-F. Wu, K. Toma, X.-Y. Wang, P. Mészáros, On the High-energy Emission of the Short GRB 090510, Astrophys.J. 733 (2011) 22.
• (43) K. Toma, X.-F. Wu, P. Mészáros, An Up-Scattered Cocoon Emission Model of Gamma-Ray Burst High-Energy Lags, Astrophys.J. 707 (2009) 1404–1416.
• (44) K. Toma, X.-F. Wu, P. Mészáros, Photosphere-internal shock model of gamma-ray bursts: case studies of Fermi/LAT bursts, M.N.R.A.S. 415 (2011) 1663–1680.
• (45) I. Vurm, R. Hascoet, A. M. Beloborodov, Pair-dominated GeV-optical flash in GRB 130427A, ArXiv e-printsarXiv:1402.2595.
• (46) P. Mészáros, M. J. Rees, GeV Emission from Collisional Magnetized Gamma-Ray Bursts, Astrophys.J.Lett. 733 (2011) L40+.
• (47) A. M. Beloborodov, Collisional mechanism for gamma-ray burst emission, M.N.R.A.S. 407 (2010) 1033–1047.
• (48) S. Razzaque, C. D. Dermer, J. D. Finke, Synchrotron radiation from ultra-high energy protons and the Fermi observations of GRB 080916C, ArXiv e-printsarXiv:0908.0513.
• (49) K. Murase, K. Asano, T. Terasawa, P. Mészáros, The Role of Stochastic Acceleration in the Prompt Emission of Gamma-Ray Bursts: Application to Hadronic Injection, Astrophys.J. 746 (2012) 164.
• (50) E. Waxman, Cosmological Gamma-Ray Bursts and the Highest Energy Cosmic Rays, Physical Review Letters 75 (1995) 386–389.
• (51) E. Waxman, J. Bahcall, High energy neutrinos from cosmological gamma-ray fireballs, Phys. Rev. Lett.  78 (1997) 2292–2295.
• (52) D. Guetta, D. Hooper, J. Alvarez-Muñiz, F. Halzen, E. Reuveni, Neutrinos from individual gamma-ray bursts in the BATSE catalog, Astroparticle Physics 20 (2004) 429–455.
• (53) R. Abbasi, Y. Abdou, T. Abu-Zayyad, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, D. Altmann, K. Andeen, J. Auffenberg, et al., An absence of neutrinos associated with cosmic-ray acceleration in -ray bursts, Nature 484 (2012) 351–354.
• (54) M. Ahlers, M. C. Gonzalez-Garcia, F. Halzen, GRBs on probation: Testing the UHE CR paradigm with IceCube, Astroparticle Physics 35 (2011) 87–94.
• (55) R. Abbasi, Y. Abdou, T. Abu-Zayyad, J. Adams, J. A. Aguilar, M. Ahlers, K. Andeen, J. Auffenberg, X. Bai, M. Baker, et al., Limits on Neutrino Emission from Gamma-Ray Bursts with the 40 String IceCube Detector, Physical Review Letters 106 (14) (2011) 141101–+.
• (56) S. Hümmer, P. Baerwald, W. Winter, Neutrino Emission from Gamma-Ray Burst Fireballs, Revised, Physical Review Letters 108 (23) (2012) 231101.
• (57) Z. Li, Note on the normalization of predicted gamma-ray burst neutrino flux, Phys.Rev.D 85 (2) (2012) 027301.
• (58) S. Gao, K. Asano, P. Meszaros, High Energy Neutrinos from Dissipative Photospheric Models of Gamma Ray Bursts, Jour. Cosmology and Astro-Particle Phys. 11 (2012) 58.
• (59) K. Asano, P. Meszaros, Neutrino and Cosmic-Ray Release from Gamma-Ray Bursts: Time-Dependent Simulations, ArXiv e-printsarXiv:1402.6057.
• (60) K. Murase, Prompt high-energy neutrinos from gamma-ray bursts in photospheric and synchrotron self-Compton scenarios, Phys.Rev.D 78 (10) (2008) 101302.
• (61) X.-Y. Wang, Z.-G. Dai, Prompt TeV Neutrinos from the Dissipative Photospheres of Gamma-ray Bursts, Astrophys.J.Lett. 691 (2009) L67–L71.
• (62) P. Veres, B.-B. Zhang, P. Mészáros, Magnetically and Baryonically Dominated Photospheric Gamma-Ray Burst Model Fits to Fermi-LAT Observations, Astrophys.J. 764 (2013) 94.
• (63) J. M. Burgess, R. D. Preece, F. Ryde, P. Veres, P. Mészáros, V. Connaughton, M. Briggs, A. Pe’er, S. Iyyani, A. Goldstein, M. Axelsson, M. G. Baring, P. N. Bhat, D. Byrne, G. Fitzpatrick, S. Foley, D. Kocevski, N. Omodei, W. S. Paciesas, V. Pelassa, C. Kouveliotou, S. Xiong, H.-F. Yu, B. Zhang, S. Zhu, An Observed Correlation between Thermal and Non-thermal Emission in Gamma-Ray Bursts, Astrophys.J.Lett. 784 (2014) L43.
• (64) S. Horiuchi, K. Murase, K. Ioka, P. Mészáros, The Survival of Nuclei in Jets Associated with Core-collapse Supernovae and Gamma-Ray Bursts, Astrophys.J. 753 (2012) 69.
• (65) Pierre Auger Collaboration, Anisotropy and chemical composition of ultra-high energy cosmic rays using arrival directions measured by the Pierre Auger Observatory, Jour. Cosmology and Astro-Particle Phys. 6 (2011) 22–+.
• (66) IceCube Collaboration, M. G. Aartsen, R. Abbasi, Y. Abdou, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, D. Altmann, J. Auffenberg, et al., The IceCube Neutrino Observatory Part I: Point Source Searches, ArXiv e-prints (2013) p48, arXiv:1309.6979.
• (67) S. Gao, K. Kashiyama, P. Mészáros, On the Neutrino Non-detection of GRB 130427A, Astrophys.J.Lett. 772 (2013) L4.
• (68) S. E. Woosley, J. S. Bloom, The Supernova Gamma-Ray Burst Connection, Annu.Rev.Astron.Astrophys. 44 (2006) 507–556.
• (69) K. Nomoto, M. Tanaka, N. Tominaga, K. Maeda, Hypernovae, gamma-ray bursts, and first stars, New Ast.Rev. 54 (2010) 191–200.
• (70) D. Guetta, M. Della Valle, On the Rates of Gamma-Ray Bursts and Type Ib/c Supernovae, Astrophys.J.Lett. 657 (2007) L73–L76.
• (71) A. M. Soderberg, S. Chakraborti, G. Pignata, R. A. Chevalier, P. Chandra, A. Ray, M. H. Wieringa, A. Copete, V. Chaplin, V. Connaughton, S. D. Barthelmy, M. F. Bietenholz, N. Chugai, M. D. Stritzinger, M. Hamuy, C. Fransson, O. Fox, E. M. Levesque, J. E. Grindlay, P. Challis, R. J. Foley, R. P. Kirshner, P. A. Milne, M. A. P. Torres, A relativistic type Ibc supernova without a detected -ray burst, Nature 463 (2010) 513–515.
• (72) S. Campana, V. Mangano, A. J. Blustin, P. Brown, D. N. Burrows, G. Chincarini, J. R. Cummings, G. Cusumano, M. Della Valle, D. Malesani, P. Mészáros, J. A. Nousek, M. Page, T. Sakamoto, E. Waxman, B. Zhang, Z. G. Dai, N. Gehrels, S. Immler, F. E. Marshall, K. O. Mason, A. Moretti, P. T. O’Brien, J. P. Osborne, K. L. Page, P. Romano, P. W. A. Roming, G. Tagliaferri, L. R. Cominsky, P. Giommi, O. Godet, J. A. Kennea, H. Krimm, L. Angelini, S. D. Barthelmy, P. T. Boyd, D. M. Palmer, A. A. Wells, N. E. White, The association of GRB 060218 with a supernova and the evolution of the shock wave, Nature 442 (2006) 1008–1010.
• (73) D. Lazzati, B. J. Morsony, C. H. Blackwell, M. C. Begelman, Unifying the Zoo of Jet-driven Stellar Explosions, Astrophys.J. 750 (2012) 68.
• (74) P. Mészáros, E. Waxman, TeV Neutrinos from Successful and Choked Gamma-Ray Bursts, Physical Review Letters 87 (17) (2001) 171102–+.
• (75) S. Horiuchi, S. Ando, High-energy neutrinos from reverse shocks in choked and successful relativistic jets, Phys.Rev.D 77 (6) (2008) 063007.
• (76) K. Murase, K. Ioka, TeV-PeV Neutrinos from Low-Power Gamma-Ray Burst Jets inside Stars, Physical Review Letters 111 (12) (2013) 121102.
• (77) I. Taboada, Multi-messenger observations of gamma-ray bursts, in: J. E. McEnery, J. L. Racusin, N. Gehrels (Eds.), American Institute of Physics Conference Series, Vol. 1358 of American Institute of Physics Conference Series, 2011, pp. 365–370.
• (78) J. Daughhetee, Search for Choked GRBs using Icecube’s DeepCore, in: APS April Meeting Abstracts, 2012, p. C7008.
• (79) X.-Y. Wang, S. Razzaque, P. Meszaros, Z.-G. Dai, High-energy Cosmic Rays and Neutrinos from Semi-relativistic Hypernovae, Phys.Rev.D. in press (astro-ph/0705.0041) 705.
• (80) R. Budnik, B. Katz, A. MacFadyen, E. Waxman, Cosmic Rays from Transrelativistic Supernovae, Astrophys.J. 673 (2008) 928–933.
• (81) X.-Y. Wang, S. Razzaque, P. Mészáros, On the Origin and Survival of Ultra-High-Energy Cosmic-Ray Nuclei in Gamma-Ray Bursts and Hypernovae, Astrophys.J. 677 (2008) 432–440.
• (82) J. Arons, Magnetars in the Metagalaxy: An Origin for Ultra-High-Energy Cosmic Rays in the Nearby Universe, Astrophys.J. 589 (2003) 871–892.
• (83) B. Zhang, Z. G. Dai, P. Mészáros, E. Waxman, A. K. Harding, High-Energy Neutrinos from Magnetars, Astrophys.J. 595 (2003) 346–351.
• (84) G. Ghisellini, G. Ghirlanda, F. Tavecchio, F. Fraternali, G. Pareschi, Ultra-high energy cosmic rays, spiral galaxies and magnetars, M.N.R.A.S. 390 (2008) L88–L92.
• (85) K. Murase, P. Mészáros, B. Zhang, Probing the birth of fast rotating magnetars through high-energy neutrinos, Phys.Rev.D 79 (10) (2009) 103001–+.
• (86) K. Kotera, Ultrahigh energy cosmic ray acceleration in newly born magnetars and their associated gravitational wave signatures, Phys.Rev.D 84 (2) (2011) 023002–+.
• (87) J. P. Rachen, T. Stanev, P. L. Biermann, Extragalactic ultrahigh-energy cosmic rays. 2. comparison with experimental data, Astron. Astrophys. 273 (1993) 377.
• (88) V. Berezinsky, A. Z. Gazizov, S. I. Grigorieva, Signatures of AGN model for UHECR, ArXiv Astrophysics e-printsarXiv:arXiv:astro-ph/0210095.
• (89) F. W. Stecker, C. Done, M. H. Salamon, P. Sommers, High-energy neutrinos from active galactic nuclei, Physical Review Letters 66 (1991) 2697–2700.
• (90) J. Alvarez-Muñiz, P. Mészáros, High energy neutrinos from radio-quiet active galactic nuclei, Phys.Rev.D 70 (12) (2004) 123001–+.
• (91) A. Pe’er, K. Murase, P. Mészáros, Radio-quiet active galactic nuclei as possible sources of ultrahigh-energy cosmic rays, Phys.Rev.D 80 (12) (2009) 123018–+.
• (92) U. Keshet, E. Waxman, A. Loeb, V. Springel, L. Hernquist, Gamma Rays from Intergalactic Shocks, Astrophys.J. 585 (2003) 128–150.
• (93) S. Inoue, F. A. Aharonian, N. Sugiyama, Hard X-Ray and Gamma-Ray Emission Induced by Ultra-High-Energy Protons in Cluster Accretion Shocks, Astrophys.J.Lett. 628 (2005) L9–L12.
• (94) K. Murase, S. Inoue, S. Nagataki, Cosmic Rays above the Second Knee from Clusters of Galaxies and Associated High-Energy Neutrino Emission, Astrophys.J.Lett. 689 (2008) L105–L108.
• (95) K. Kashiyama, P. Mészáros, Galaxy Mergers as a Source of Cosmic Rays, Neutrinos, and Gamma Rays, Astrophys.J.Lett. 790 (2014) L14.
• (96) M. G. Aartsen, R. Abbasi, Y. Abdou, M. Ackermann, J. Adams, J. A. Aguilar, M. Ahlers, D. Altmann, J. Auffenberg, X. Bai, et al., First Observation of PeV-Energy Neutrinos with IceCube, Physical Review Letters 111 (2) (2013) 021103.
• (97) A. Loeb, E. Waxman, The cumulative background of high energy neutrinos from starburst galaxies, Jour. Cosmology and Astro-Particle Phys. 5 (2006) 3.
• (98) D. B. Fox, K. Kashiyama, P. Mészarós, Sub-PeV Neutrinos from TeV Unidentified Sources in the Galaxy, Astrophys.J. 774 (2013) 74.
• (99) H.-N. He, T. Wang, Y.-Z. Fan, S.-M. Liu, D.-M. Wei, Diffuse PeV neutrino emission from ultraluminous infrared galaxies, Phys.Rev.D 87 (6) (2013) 063011.
• (100) R.-Y. Liu, X.-Y. Wang, S. Inoue, R. Crocker, F. Aharonian, Diffuse PeV neutrinos from EeV cosmic ray sources: Semirelativistic hypernova remnants in star-forming galaxies, Phys.Rev.D 89 (8) (2014) 083004.
• (101) K. Murase, M. Ahlers, B. C. Lacki, Testing the hadronuclear origin of PeV neutrinos observed with IceCube, Phys.Rev.D 88 (12) (2013) 121301.
• (102) A. A. Abdo, Fermi LAT Collaboration, Spectrum of the Isotropic Diffuse Gamma-Ray Emission Derived from First-Year Fermi Large Area Telescope Data, Physical Review Letters 104 (10) (2010) 101101.
• (103) X.-C. Chang, X.-Y. Wang, The diffuse gamma-ray flux associated with sub-PeV/PeV neutrinos from starburst galaxies, ArXiv e-printsarXiv:1406.1099.
• (104) L. A. Anchordoqui, T. C. Paul, L. H. M. da Silva, D. F. Torres, B. J. Vlcek, What IceCube data tell us about neutrino emission from star-forming galaxies (so far), ArXiv e-printsarXiv:1405.7648.
• (105) The Pierre AUGER Collaboration, V. Scherini, Search for primary photons and neutrinos in the ultra-high energy cosmic rays with the Pierre Auger Observatory, Nuclear Physics B Proceedings Supplements 212 (2011) 115–120.
• (106) The Pierre Auger Collaboration, Correlation of the Highest-Energy Cosmic Rays with Nearby Extragalactic Objects, Science 318 (2007) 938–.
• (107) The Pierre Auger Collaboration, A. Aab, P. Abreu, M. Aglietta, M. Ahlers, E.-J. Ahn, I. Albuquerque, I. Allekotte, J. Allen, P. Allison, et al., The Pierre Auger Observatory: Contributions to the 33rd International Cosmic Ray Conference (ICRC 2013), ArXiv e-printsarXiv:1307.5059.
• (108) K.-H. Kampert, Ultrahigh-Energy Cosmic Rays: Results and Prospects, Brazilian Journal of Physics 43 (2013) 375–382.
• (109) R.-Y. Liu, X.-Y. Wang, Diffuse PeV Neutrinos from Gamma-Ray Bursts, Astrophys.J. 766 (2013) 73.
• (110) B. Katz, E. Waxman, T. Thompson, A. Loeb, The energy production rate density of cosmic rays in the local universe is $sim10^44-45rm erg~Mpc^-3~yr^-1$ at all particle energies, ArXiv e-printsarXiv:1311.0287.
• (111) J. Bahcall, E. Waxman, High energy astrophysical neutrinos: The upper bound is robust, Phys.Rev.D 64 (2) (2001) 023002–+.
You are adding the first comment!
How to quickly get a good reply:
• Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
• Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
• Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters