Advances in High Energy Physics

Volume 2012 (2012), Article ID 586413, 41 pages

doi:10.1155/2012/586413

Review Article

Neutrinos from Cosmic Accelerators including Magnetic Field and Flavor Effects

Walter Winter

Institut für Theoretische Physik und Astrophysik, Universität Würzburg, 97074 Würzburg, Germany

Received 2 February 2012; Accepted 5 April 2012

Academic Editor: Guey-Lin Lin

Copyright © 2012 Walter Winter. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

We review the particle physics ingredients affecting the normalization, shape, and flavor composition of astrophysical neutrinos fluxes, such as different production modes, magnetic field effects on the secondaries (muons, pions, and kaons), and flavor mixing, where we focus on interactions. We also discuss the interplay with neutrino propagation and detection, including the possibility to detect flavor and its application in particle physics, and the use of the Glashow resonance to discriminate from interactions in the source. We illustrate the implications on fluxes and flavor composition with two different models: (1) the target photon spectrum is dominated by synchrotron emission of coaccelerated electrons and (2) the target photon spectrum follows the observed photon spectrum of gamma-ray bursts. In the latter case, the multimessenger extrapolation from the gamma-ray fluence to the expected neutrino flux is highlighted.

1. Introduction

In addition to gamma-ray and cosmic ray instruments, neutrino telescopes, such as IceCube [1] or ANTARES [2], provide interesting data on the sources of the highest-energetic particles found in the universe, so-called “cosmic accelerators”; see [3–6] for reviews. In particular, neutrinos are a prominent way to search for the origin of the cosmic rays or to discriminate between leptonic and hadronic models describing the observed spectral energy distribution of photons. There are numerous possible sources, see [7] for an overview and [8] for the general theory. Very interesting extragalactic candidates for neutrino and cosmic ray production may be gamma-ray bursts (GRBs) [9] and active galactic nuclei (AGNs) [10–12]. The most stringent bounds for these sources, which are expected to be roughly uniformly distributed over the sky, so far come from IceCube, which has recently released data on time-integrated [13] and time-dependent [14] point source searches, GRB neutrino searches [15], and diffuse flux searches [16]. So far, no astrophysical neutrinos have been detected, which has been for a long time consistent with generic Waxman-Bahcall [17] and Mannheim-Protheroe-Rachen [18] bounds. However, data from IC40 and IC59, referring to the 40 and 59 string configuration of IceCube, respectively, start to significantly exceed these bounds, see [15, 16, 19], which is in tension with the corresponding neutrino production models, such as [9, 20, 21] for GRBs. For example, neutrino data may soon challenge the paradigm that GRB fireballs are the sources of the ultra-high-energy cosmic rays (UHECRs) [22]. For constraints to AGN models, see, for example, [23]. As a consequence, the age of truth has come for neutrino astrophysics, which is also the age of precision. Especially since data are available now, it is necessary to critically review the underlying assumptions from both the astrophysics and particle physics perspectives and to develop the models from rough analytical estimates into more accurate numerical predictions.

In this paper, we focus on the minimal set of particle physics ingredients for the neutrino production, which must be present in virtually all sources, using several specific examples. We focus on photohadronic () interactions for the meson production, with the exception of Section 4.3. We do not only discuss the predicted neutrino flux, but also the flavor and neutrino-antineutrino composition at source and detector. The discussed effects include the following.

(i) Additional pion production modes, such as -channel (direct) and multipion production; see, for example, [24–30].(ii) Neutron and kaon production; see, for example, [29–36] (in general, there is also an additional contribution from charmed meson production, see [32, 33] for a detailed comparison).(iii) The cooling and decay of secondaries (pions, muons, and kaons); see, for example, [29, 30, 33, 34, 37–39].(iv) Flavor mixing and possible new physics effects; see [40] for a review.(v) The helicity-dependence of the muon decays; see, for example, [28, 38] (this effect has been discussed earlier in the context of atmospheric neutrinos, see, for instance, [41–43]).(vi) Spectral effects, such as the energy dependence of the mean free path of the protons, and their impact on the prediction; see, for example, [44, 45].(vii) The impact of the maximal proton energy on the neutrino spectrum; see, for example, [36].(viii) Deviations from the frequently used neutrino flux assumption; see, for example, [46].

While many of these effects have been studied elsewhere in the literature, we mainly show examples generated with the NeuCosmA (Neutrinos from Cosmic Accelerators) software in this paper to present them in a self-consistent way.

The structure of this paper is as follows: in Section 2, we give a simplified picture for the connection among neutrinos, cosmic rays, and gamma-rays. Then in Section 3, we review the minimal set of ingredients for neutrino production from the particle physics perspective. In Section 4, we discuss neutrino propagation and detection, including the possibility to detect flavor and the use of the Glashow resonance, where we illustrate how new physics can be tested in the neutrino propagation in Section 4.4. We furthermore present two specific applications: a generic AGN-like model in Section 5 and a model for GRBs in Section 6, where the main difference is the model for the target photons. Then we finally summarize in Section 7.

2. Neutrinos and the Multimessenger Connection

Here we outline a simplified picture of the neutrino or cosmic ray source, as often used in the literature, whereas we add extra ingredients in the next section. In this approach, charged mesons originate from or interactions, where we focus on (photohadronic) interactions in this work; see, for example, [47, 48] for interactions, which may be dominant for particular source classes, such as supernova remnants. In the simplest possible picture, charged pions are produced by the -resonance

While this process is not sufficient for state-of-the-art models for neutrino production, it is very useful to illustrate a few qualitative points common to many cosmic ray and neutrino production models. The protons on the l.h.s. of (2.1) are typically assumed to be injected into the interaction volume with an spectrum (here primed parameters refer to the shock rest frame (SRF), whereas unprimed parameters to the observer’s frame) coming from Fermi shock acceleration, where . They interact with the photons on the l.h.s. of (2.1) with energy . While the assumptions for the injected protons are similar for most models (except from the minimal and maximal energies), the target photons are typically described in a model- and source-dependent way, such as by

(1) synchrotron emission from coaccelerated electrons or positrons,(2) thermal emission, such as from an accretion disk,(3) a more complicated combination of radiation processes,(4) an estimate inferred from the gamma-ray observation,

just to name a few examples. While any realistic simulation of a particular source will imply option (3), the other options typically rely on fewer parameters and may be good approximations in many cases. In particular, a reliable prediction for the photon density in the source may be obtained from the gamma-ray observation, option (4), if the photons can escape. In fact, we will use option (1) in Section 5 and option (4) in Section 6.

After an interaction between proton and photon, the particles on the r.h.s. of (2.1) are produced with the given branching ratios. The neutrinos then originate from decays via the decay chain

where in this standard picture are produced in the ratio 1 : 2 : 0 if the polarities (neutrinos and antineutrinos) are added. In addition, high-energy gamma-rays are produced by

These are typically emitted from the source at lower energies due to electromagnetic cascades, in addition to gamma-rays escaping from the interaction volume (the ones contributing on the l.h.s. of (2.1)).

From (2.1), we can also illustrate the production of cosmic ray protons, ignoring for the moment that the composition of cosmic rays may be heavier at high energies [49]. First of all, some of the protons injected into the interaction volume on the l.h.s. of (2.1) may escape, leading to cosmic ray production. However, even if the protons are magnetically confined, the neutrons on the r.h.s. of (2.1), which are electrically neutral, can easily escape if the source is optically thin to neutron escape. After decay (typically outside the source)

they lead to cosmic ray flux and an additional neutrino flux which is an unavoidable consequence of the interactions in (2.1). The cosmic ray protons with energies above  eV interact with the cosmic microwave background (CMB) photons by (2.1), leading to the so-called Greisen-Zatsepin-Kuzmin (GZK) cutoff [50, 51]. However, according to (2.1), charged pions are produced in these interactions as well, which means that an additional neutrino flux should come with that, which is often called “cosmogenic neutrino flux.”

In summary, the photohadronic interaction in (2.1) offers a self-consistent picture for a cosmic ray source, with a possible connection among cosmic ray, neutrino, and gamma-ray escape. In specific models, however, that does not mean that a large neutrino flux is guaranteed for every cosmic ray source. For instance, the interaction rate for the process in (2.1), which depends on the photon density, may be low.

3. Simulation of Neutrino Sources

Here we give a more detailed generic picture of the simulation of neutrino sources from the particle physics perspective with the minimal set of ingredients, without using a specific model. A flowchart summarizing the contents of this section, which can be followed during the reading, is given in Figure 6.

3.1. Photohadronic Interactions

In order to describe the processes within an interaction volume (one zone in the simplest case), two kinds of spectra are needed: (in units of GeV⁻¹ cm⁻³ s⁻¹) describes the number of particles injected or ejected per volume element and energy interval, and (in units of GeV⁻¹ cm⁻³) describes the particle density per energy interval. The secondary meson injection rate for a pion or kaon species produced in photohadronic interactions is given by (following [28])

Here is the fraction of energy going into the secondary, (here can be related to the center-of-mass energy by , where is the angle between proton and photon momentum; corresponds to heads-on collisions), and is the “response function.” If many interaction types are considered, the response function can be quite complicated. However, if it is known from particle physics, (3.1) can be used to compute the secondary injection for arbitrary proton and photon spectra. The important point here is that the secondary production depends on the product normalization of the proton density and the target photon density within the interaction volume. Thus a higher proton density can be compensated by a lower photon density, and vice versa (strictly speaking, this degeneracy only holds in the absence of any other radiation process, e.g., inverse Compton scattering depends on the photon density individually, which means that the input spectral shapes (especially in (3.1)) will be modified if that process contributes significantly). Another implication of (3.1) is that the secondary production depends on the densities within the source , not the injection rates . Of course, one cannot look into the source, but can only observe cosmic messengers escaping from the source. As we will demonstrate later, the observed/ejected photon or cosmic ray spectrum is only directly representative for the corresponding density spectrum within the source if “trivial” escape is the leading process, that is, with and the size of the interaction region. For this section, (3.1) is used as a starting point for the computation of the neutrino fluxes, where we do not discuss the origin of spectral shape and normalization of and . In practice, typically an injection spectrum is assumed for the protons, as mentioned above, where the maximal energy is limited by synchrotron and adiabatic losses. The photon density may be a consequence of a complicated interplay of radiation processes. In either case, the derivation of these densities depends on the model, and we will show several examples in Sections 5 and 6.

Compared to a numerical approach, (2.1) has limitations to describe the meson production. First of all, additional pion production modes contribute, such as higher resonances, direct (-channel), and multipion production, which will also lead to production (cf. Figure 6). These are not as easy to describe as the Δ-resonance because of different shapes of the cross sections and more complicated kinematics. The Monte Carlo package SOPHIA [25] can deal with these interactions. In order to increase the efficiency, often parameterizations of SOPHIA are used, such as [28, 52]. In the following, we use [28] (model Sim-B), because the secondary muons and pions are needed explicitly. We show the impact of the resonant production (including higher resonances) on production in Figure 1 for a typical GRB, an AGN, and a (10 eV) black body (BB) target photon field. As one can read off from this figure, the resonances always give a reasonable first estimate for the actual pion production, but quantitatively they only dominate at the breaks. In addition, multipion processes can change the spectral shape significantly, such as for the GRB example, which is a consequence of the cross section dependence on the center-of-mass energy. As a further limitation, note that (2.1) does not describe kaon production and subsequent decay into neutrinos, where the leading modes are given by

(a)

(b)

(c)

Figure 1: Contributions of different production modes to the total spectrum for a typical GRB, an AGN, and a (10 eV) black body (BB) target photon density, figure taken from [28].

The branching ratio for the leading channel in (3.3) is about 64%. The second-most-important decay mode is (20.7%). The other decay modes account for 16%, no more than about 5% each. Because interesting effects can only be expected in the energy range with the most energetic neutrinos, we only use the direct decays from the leading mode.

In the literature, the Δ-resonance approximation in (2.1) is, even in analytical approaches, typically not taken literally. For example, a simple case is the approximation by Waxman and Bahcall [9] (“WB Δ-approx.”), for which one can write the response function as

which implies that charged pions are produced in 50% of all cases, and these take 20% of the proton energy. In addition, the width of the Δ-resonance is taken into account (note that this description is slightly more accurate than [9], which uses an additional integral approximation). In fact, this function peaks at  GeV, which is higher than the threshold for photohadronic interactions—and it is even a little bit higher in the numerical calculation. The reason is that for the threshold often head-on collisions are assumed (), whereas these only contribute a small part to the total number of interactions. Using (3.4) in (3.1) and rewriting the integral over in one over , it is easy to show that for power law spectra and

This means that the pion spectral index depends on both the proton and photon spectra, where it is inversely proportional to the photon spectral index. As we will see below, the neutrino spectrum follows the pion spectrum, which means that the assumption of an spectrum for the neutrinos, as it is often used in data analysis and many models in the literature, is only a valid assumption for —which is roughly observed for GRBs below the break. On the other hand, if the target photons come from synchrotron emission, such a hard photon spectrum is not possible, and may be more plausible if the electrons are injected with a spectral index similar to the protons. As a consequence, the neutrino spectrum becomes harder. In addition, multipion processes in the photohadronic interactions will act in the same direction and make the neutrino spectrum even harder, compare Figure 1; see also [29] for a detailed comparison between the approximation in (3.4) and the numerics. Note that for interactions with “cold” (non-relativistic) protons, the assumption may be plausible [47]. We discuss the implications of the assumption for the detector response in Section 4.

3.2. Decays of Secondaries

The weak decays of pions and muons are described in detail in [38]. In general, in case of ultra-relativistic parents of type , the distribution of the daughter particle of type takes a scaling form in order to obtain for the energy spectra

summed over all parent species. The functions for pion-, kaon-, and helicity-dependent muon decays can be read off from [38] (Section 4). Consider the simplified case of a -function for . For instance, for decays of neutrons in (2.4), one may approximate

with . If the neutrons do not interact, (see below), and we find for the neutrino injection  (3.6)

In this case, the neutrino spectrum follows the neutron spectrum. If the neutrinos originate from pion decays, the neutrino injection follows the pion injection spectrum by similar arguments.

Since the decays of muons are helicity dependent, that is, is different for left- and right-handed muons, it is necessary to keep track of these two species separately (cf. Figure 6). Although the effects of the helicity-dependent muon decays on the fluxes are probably small, the flavor composition is slightly affected, depending on the parameters of the source. We illustrate this effect in Figure 2 for the GRB, AGN, and black-body examples in Figure 1. In this figure, the horizontal line corresponds to the standard assumption, that is, neutrinos being produced in the flavor composition of 1 : 2 : 0. Including the scaling of the secondary decays, the dotted (blue) curves are obtained if the helicity of the muons is averaged over. From the comparison with the light (red) solid curves, it is clear that this assumption is implemented in SOPHIA. On the other hand, as pointed out in [38], keeping track of the muon helicity slightly changes the flavor composition, see dark (black) solid curves, which are significantly different from the standard assumption and the helicity-averaged version. However, it is also clear from Figure 2 that the deviation from the standard prediction depends on the input spectra, as it is smaller for the AGN than the GRB example. The dashed (green) curves show the contributions of the neutron and kaon decays, which affect the flavor composition at very low and high energies, respectively.

(a)

(b)

(c)

Figure 2: Flux ratio at the source between electron and muon neutrinos as a function of the neutrino energy, for the same objects as in Figure 1. No secondary cooling is included in this figure, figure taken from [28].

3.3. Cooling of Secondaries

In order to describe the cooling of the secondary pions, muons, and kaons, we use the steady state approach; that is, we do not allow for an explicit time-dependence since the statistics of neutrino observations is typically expected to be low. The steady state equation for the particle spectrum, assuming continuous energy losses, is given by

with the characteristic escape time, with the rate characterizing energy losses. This differential equation balances the particle injection on the l.h.s. with energy losses and escape on the r.h.s. of the equation. Note again that the steady density is needed for the photohadronic interactions in (3.1), not the actual injection spectrum. In addition, note that if there are no energy losses (), one has immediately from (3.9), which we have already used above. If decay is the dominant escape mechanism, one finds .

While the primary proton and photon spectra in (3.1) could be affected by a number of radiation processes, for the neutrino fluxes and flavor compositions, at least the processes of the secondaries (pions, muons, kaons) are important. We illustrate the synchrotron cooling and decay rates for pions, muons, and kaons as a function of energy in Figure 3. As one can easily see in the figure, for any species, decay dominates at low energies, while synchrotron cooling dominates at high energies. Other cooling or escape processes are often subdominant, as illustrated by the gray-shaded region for an adiabatic cooling components in GRBs. The two curves meet at a critical energy for each species, which is different depending on the particle physics parameters. As a consequence of (3.9), the corresponding steady spectra are loss-steepened by two powers above

where synchrotron and decay rates are equal. These critical energies depend on the particle physics properties of the parent, that is, the mass and the rest frame lifetime , and the magnetic field as the only astrophysical parameter. It is therefore a very robust prediction and might allow for the only direct measurement of . Rescaling the magnetic field shifts the critical energies by a constant amount on the horizontal axis, but does not change the spacing in the logarithmic picture. In Figure 3, we also show the estimated range for adiabatic cooling as additional cooling component (shaded region), which may have some impact, especially on the muons, in extreme cases. In these cases, the height of the spectral peaks will be somewhat reduced, but the qualitative picture does not change.

Figure 3: Synchrotron cooling and decay rates for pions, muons, and kaons as a function of energy for a magnetic field of  kG. The gray-shaded region shows the estimated range for the escape (or adiabatic cooling) rate for GRBs for a variability timescale ranging from 0.01 s to 1 s, figure taken from [30].

We show in Figure 4 the consequences for the spectral shape at one GRB example. Here the injection spectrum of electron neutrinos Figure 4(a) and muon neutrinos Figure 4(b) is shown, including the individual contributions from the parents. First of all, one can read off from this figure that in the case of no losses (dashed curves) the spectral shapes of all contributions are very similar, and the neutrino fluxes add in a trivial manner. A change of the primary spectra and in (3.1) may change the shape of the dashed curves, but almost in the same way for all curves. If the synchrotron losses are switched on (solid curves), the spectral split predicted by (3.10) (see also Figure 3) among the neutrino spectra coming from different parent species can be clearly seen in Figure 4(b). One can also see a small pile-up effect coming from the muon decays; that is, a small region where the cooled muons coming from higher energies pile up and dominate, and lead to a higher flux than in the “no losses” case. In Figure 4(a), only two spectra are shown, since only muon or neutron decays may produce the electron flavor. Because of the very small in (3.8), the neutron decays only show up at low energies. Since the decay of pions, muons, or kaons, which dominate in different energy ranges, leads to different neutrino flavor compositions, the flavor composition at the source will be changed as a function of energy. In the literature, often the following source classes in terms of are distinguished (see [53] for a summary, including other options leading to a neutron beam or muon beam-like source).

(a)

(b)

Figure 4: Injection spectra for electron neutrinos (a) and muon neutrinos (b) at the source for a GRB example, where the individual contributions from the parents are shown (neutrinos and antineutrinos added), figure taken from [30].

Pion Beam Source

Pion decays and muon decays equally contribute; flavor composition 1 : 2 : 0.

Muon Damped Source

Strong muon cooling, which means that pion decays dominate; flavor composition 0 : 1 : 0.

Muon Beam Source

Pile-up muons dominate; flavor composition 1 : 1 : 0.

Neutron Beam Source

Neutron decays dominate; flavor comp. 1 : 0 : 0.

Undefined Source

Several of these processes compete; flavor composition ().

Note that in none of these cases are produced in significant quantities, which means that it is sufficient to use the ratio between the and production at the source, sometimes called “flavor ratio.” From the preceding discussion, it is clear that the above classifications can only hold in specific energy ranges. However, all of these sources can be recovered in a numerical simulation. In Figure 5, four examples for energy-dependent flavor ratios at the source are shown for different parameter sets, where in this model is generated by synchrotron losses of coaccelerated electrons. In this figure, also the relevant flavor ratio ranges for the different sources introduced above are shown. Figure 5(b) shows the classical pion beam source, which is typically found for low magnetic fields. Nevertheless, the contribution of neutron decays at low energies can be clearly seen. Figure 5(d) shows a pion beam evolving in a muon damped source at high energies, as, for example, in [37]. Figure 5(a) depicts a muon beam to muon damped source. In this case, the cooled muons pile up at lower energies, where the muon decays dominate. And Figure 5(c) shows an undefined source, where several processes compete.

(a)

(b)

(c)

(d)

Figure 5: Energy-dependent flavor ratio at the source for several selected parameter sets (test points “TP”; cf. Figure 10) for a model where is generated by synchrotron losses of coaccelerated electrons. The (dark) gray-shaded areas mark the regions where the fluxes peak, figure taken from [34].

Figure 6: Flowchart describing the model (in SRF). The functions denote (injection) spectra per time frame [(GeV cm³ s)⁻¹] and steady spectra [(GeV cm³)⁻¹] derived from the balance between injection and losses or escape. Dashed arrows stand for solving the steady state differential equation (3.9), figure taken from [30].

Of course, not only the secondaries are affected by synchrotron (or adiabatic) losses, but also the primary protons. Depending on the model, one can use these losses to determine the maximal proton energy, or one can put in the maximal proton energy by hand. From Figure 4, it is interesting to discuss the impact of the maximal proton energy on the neutrino fluxes. In the “no losses” case (dashed curves), the maximal neutrino energy is directly determined by the maximal proton energy . There is, however, one exception: the neutrinos from neutron decays are limited by , cf., (3.8). If the synchrotron losses are switched on (solid curves), the neutrino spectrum from neutron decays still follows the proton spectrum, since the neutrons are electrically neutral, whereas the maximal neutrino energies for the other production modes are determined by (3.10). As a consequence, the neutron decay spectrum strongly depends on the assumptions for , whereas the other spectral shapes are entirely unaffected by as long as for kaons (see (3.10), factor six from kaon production and decay kinematics). It is therefore not surprising that for strong enough magnetic fields one can find parameter sets for which the neutron decays dominate (cf. [35, 36]), but one should keep in mind that this depends on the assumptions for the maximal proton energy (and the inclusion of multipion processes, etc., which may mask this effect).

3.4. Transformation into Observer’s Frame

The transformation of the injection spectrum of the neutrinos from the source to the observable flux of (in units of GeV⁻¹ cm⁻² s⁻¹) at the Earth is given by

where a simple Lorentz boost is used (instead of a viewing angle-dependent Doppler factor). Here is the transition probability , discussed in Section 4, and is a (model-dependent) normalization factor. For example, if an isotropically emitting spherical zone is boosted with towards the observer, then since the emission is boosted into a cone with opening angle . For a relativistically expanding fireball, it is simpler to perform the transformation in a different way, see Section 6. Furthermore, is the luminosity distance, and is the comoving distance. From (3.11), one can read off that the redshift dependence of the neutrino luminosity scales as independent of the model (for fixed ), which is expected. Note that in (3.11) the neutrino and antineutrino fluxes are often added if the detector cannot distinguish these.

3.5. Summary of Ingredients, and Limitations of the Approach

We summarize the generic neutrino production in Figure 6. As one can see in this flowchart, the starting points are the proton and photon densities within the source. Once these (and ) are fixed, the rest is just a particle physics consequence. Therefore, for the computation of specific neutrino fluxes, the main effort is actually to determine , , and .

Of course, there are some processes not taken into account in this picture, which may add to the ingredients discussed above for specific source classes. For instance, secondary neutrons, produced in the photohadronic interactions, may interact again if the source is optically thick to neutrons, the secondary pions, muons, and kaons may be reaccelerated [54], synchrotron photons of the secondaries may add to , and so forth. In addition, the neutrino spectrum may be more complicated in multizone models, since these naturally allow for more freedom.

From a particle physics perspective, additional kaon and charmed meson production modes may be added, and also the secondaries may interact again (see [32, 33]). However, compared to analytical computations, the numerical approach described in this section already takes into account the secondary cooling in a self-consistent way, additional neutrino production modes can be easily included, and the full energy dependencies can be accounted for. For example, it has been demonstrated in [44] that all necessary ingredients to reproduce the analytical GRB fireball (neutrino) calculations in [9, 15, 20, 21] are contained. It should represent the minimal set of of ingredients for neutrino production which are present in every source in the spirit of constructing the simplest possible model first. Of course, if is small, the secondary cooling effects will be small as well, which is automatically included.

4. Neutrino Propagation and Detection

In this section, we discuss several aspects of neutrino propagation and detection from the theoretical perspective.

4.1. Neutrino Propagation and Observables

It is well known that neutrinos may change flavor from the production to the detection point. While this phenomenon is in general described by neutrino oscillations, astrophysical neutrinos are typically assumed to suffer from decoherence over very long distances. This means that effectively (in most practical cases) only flavor mixing enters the astrophysical neutrino propagation (see [55] for a more detailed discussion). In that case, in (3.11) becomes

for three active neutrinos, where are the usual PMNS mixing matrix elements in the standard parameterization; see, for example, [56] for recent values of the mixing angles. This implies that neutrino oscillations, that is, the -dependence, are averaged out. An initial flavor composition of 1 : 2 : 0 will therefore evolve (approximately) into 1 : 1 : 1 at the detector, see, for example, [57, 58]. In Section 4.4, we will see that (4.1) can significantly change in the presence of new physics effects, which opens new possibilities to test such effects. However, (4.1) also implies that there could be some sensitivity to standard flavor mixing, which may be complementary to Earth-based experiments, see discussions in [59–71]. In the light of the current bounds for astrophysical neutrino fluxes from IceCube, however, such applications might be unlikely.

The main observable in neutrino telescopes is muon tracks from charged current interactions of muon neutrinos, producing Cherenkov light, which can be detected in so-called digital optical modules (DOMs). Because of the long muon range which is increasing with energy, the muon track does not have to be fully contained in the detector volume, which leads to excellent statistics increasing with energy. In addition, muon tracks have a very good directional resolution (order one degree). Additional event topologies include electromagnetic (mostly from electron neutrinos) and hadronic (from tau neutrinos) cascades, as well as neutral current cascades for all flavors. For even higher energies, the tau track may be separated, leading to so-called double bang or lollipop events; see [72, 73] for an overview. In practice, the main “flavor” analysis so far performed by the IceCube collaboration has been a cascade analysis [74]. To see that, consider that electromagnetic (from ) and hadronic (from ) cascades cannot be distinguished. A useful observable is therefore the ratio of muon tracks to cascades [59]

Note that neutral current events will also produce cascades, and will also produce muon tracks in 17% of all cases, which, in practice, have to be included as backgrounds. In [74], the contribution of the different flavors to the cascade rate for a extragalactic test flux with equal contributions of all flavors at the Earth was given as electron neutrinos 40%, tau neutrinos 45%, and muon neutrinos 15% (after all cuts). This implies that charged current showers dominate and that electron and tau neutrinos are detected with comparable efficiencies, that is, that (4.2) is a good first approximation to discuss flavor at a neutrino telescope. The benefit of this flavor ratio is that the normalization of the source drops out. In addition, it represents the experimental flavor measurement with the simplest possible assumptions. For a pion beam source, one finds at the detector, for a muon damped source , and for a neutron beam source , with some dependence on the mixing angles. In principle, this and other observables, such as the ratio between electromagnetic and hadronic showers, can therefore be used to determine the flavor composition of the sources, see, for example, [75–79] (unless there are new physics effects present, see for example, [80], which may yield similar ratios at the detector for different flavor compositions at the source).

We show the ratio as a function of energy for a GRB neutrino flux in Figure 7. At the “plateau” of the flux, where the largest number of events is expected, one can clearly see the flavor transition between a pion beam and a muon damped source at higher energies. In addition, the neutron decays have an impact at very low energies. In this figure, the uncertainty on coming from the current mixing angle uncertainties is shown, as well as the expectation for the next generation of reactor and long-baseline experiments (2015), and for a high-precision neutrino oscillation facility (2025). Obviously, the current uncertainties on the mixing angles are still too large to allow for a clear identification of the flavor composition at the source, whereas already the knowledge from the next generation will allow for a flavor ratio discrimination—at least in principle.

Figure 7: Flavor ratio as a function of at the detector for a GRB neutrino flux with the Waxman-Bahcall (WB) “plateau” (in ) between and  GeV (cf. [17]). The shaded regions show the impact of the mixing angle uncertainties now (current) [81], in about 2015 (next generation, dominated by Daya Bay and T2K [82]), and in about 2025 (for a neutrino factory [83, 84]). Here is assumed for the sake of simplicity, figure taken from [30].

4.2. Detector Response and Impact of Spectral Shape

For time-integrated point source searches in IceCube [13], the simplest possible approach to describe the event rate of muon tracks is to use the exposure , where is the neutrino effective area and is the observation (exposure) time. Here is a function of the flavor or interaction type (which we do not show explicitly), the incident neutrino energy , and the declination of the source . The neutrino effective area already includes Earth attenuation effects (above PeV energies) and event selection cuts to reduce the backgrounds, which depend on the type of source considered, the declination, and the assumptions for the input neutrino flux, such as the spectral shape. Normally, the cuts are optimized for an flux. The total event rate of a neutrino telescope can be obtained by folding the input neutrino flux with the exposure as

Here is, for point sources, given in units of GeV⁻¹ cm⁻² s⁻¹ for neutrinos and antineutrinos added. If backgrounds are negligible, the 90% (Feldman-Cousins) sensitivity limit for an arbitrarily normalized input flux used in (4.3) can be estimated as [85]. This implies that a predicted flux at the level of the sensitivity limit, irrespective of the spectral shape, would lead to the same number (2.44) of events. The 90% confidence level differential limit in terms of can be defined as , see, for example, [86].

The interplay between spectral shape and detector response has, for instance, been studied in [46]. In order to illustrate that, we show in Figure 8 the limits for selected muon neutrino spectra and two different source declinations, corresponding to downgoing and upgoing events, for the time-integrated point source search in IC40. The limits are computed with (4.3), where the normalization is obtained from and the spectral shape is predicted by the model. In this figure, the thick horizontal lines show the limits for an flux (in the dominant energy range), and the thick curves the differential limits. One can easily see that the differential limits are useful, since any “well behaved” neutrino flux will stay below these limits. In addition, the differential limit shows the energy range where the instrument is most sensitive to a flux, whereas, for the horizontal lines, representing an flux limit, only the contribution close to the differential limit minimum contributes to (4.3). In Figure 8, the minimal and maximal energies for the horizontal lines are indeed arbitrarily chosen, since the limit hardly depends on these (for the chosen example, this holds as long as the energy range around the differential limit minimum within about 1.5 orders of magnitude in energy is included. Sometimes the range where 90% of the events come from is shown). From the different fluxes in Figure 8 it is clear that the interplay between spectral shape and detector response is important. For example, spectrum no. 4 will be better constrained in Figure 8(a) (downgoing events) than in Figure 8(b) (upgoing events) because the differential limit peaks at higher energies—in spite of the lower absolute performance at the differential limit minimum. The reason is the coincidence between spectral peak and differential limit minimum. This picture changes completely in the presence of strong magnetic field effects on the secondaries, see spectrum no. 2, where even the imprint of these effects becomes important. The comparison to the test flux (horizontal line) clearly demonstrates that the assumption is insufficient for all practical purposes, since it does not take into account the energy range of the flux. In addition, as we discussed around (3.5) earlier, the assumption only holds for a very special case for the target photon density. Therefore, detector response and source model are intimately connected, and it is not quite clear if some sources may be missed just because the detector response does not match the source prediction. We will come back to this issue in Section 5.

(a)

(b)

Figure 8: Limits for selected muon neutrino spectra (including flavor mixing) from [34] for different declination bands for the time-integrated point source search in IC40 (90% CL) [13]. The thick lines show the limits for an flux (in the dominant energy range), and the thick curves the differential limits, figure taken from [46]. The different spectra correspond to the test points in Figure 10 (same numbers used).

Figure 8 is also useful to illustrate the impact of the source declination, which is, for IceCube, a measure for the direction of the muon track in terms of the nadir angle because of the location at the South Pole. Note, however, that for other neutrino telescopes, such as ANTARES, this relationship is nontrivial. Obviously, the differential limit is different in the left and right panels of this figure, corresponding to downgoing and upgoing events, respectively. First of all, note that the downgoing events have to fight the atmospheric muon background, which leads to a worse performance at low energies because of appropriate cuts. However, the upgoing events suffer from the Earth attenuation for energies ≫PeV, which leads to the better sensitivity of the downgoing events for high energies, and the shift of the differential limit minimum to higher energies. From the discussion above it is clear that both event types (upgoing and downgoing) are complementary not only because they test a different part of the sky, but also because they test different energy ranges.

Finally, let us briefly comment on the statistics expected for neutrinos. First of all, it is clear from the discussion above that per definition any current IC40 limit, obtained over about one year of data taking, is compatible with about 2.4 events at the current 90% confidence limit. Assuming that the effective area increases by about a factor of three to four from IC40 to IC86 [87], one can extrapolate that the current limits are compatible with events over ten years of full IceCube operation if the current bounds are saturated. Any further nonobservation of events will reduce this maximal expectation. Therefore, one can already say right now that any conclusion about the astrophysical neutrino sources, the sources of the cosmic rays, or leptonic versus hadronic models for -ray observations will most likely be based on the information from many sources of one class. A typical example is the stacking of GRBs using their gamma-ray counterparts, such as in [15]. The aggregation of fluxes, no matter if diffused or stacked, will however imply new systematics and model-dependent ingredients, see discussion in [30].

4.3. Glashow Resonance to Discriminate from ?

A useful observable may be the Glashow resonance at around 6.3 PeV [34, 57, 88–92] to distinguish between neutrinos and antineutrinos in the detector, since this process is only sensitive to . For photohadronic () interactions, however, mostly and therefore are produced at the source, see (2.1) and (2.2), which means that no excess of events should be seen at this resonance energy—at least in the absence of flavor mixing. On the other hand, for interactions in the source, and are produced in about equal ratios, which increases the production rate. Therefore, the Glashow resonance is frequently proposed as a discriminator between and sources. Note that this argument is especially interesting for neutrino fluxes, whereas any other, (significantly different) spectral index may be a clear indicator for interactions.

In [34], the electron neutrino-antineutrino ratio at source (Section 3.3) and detector (Section 4.3) has been computed explicitly. At the source, the following observations have been made.

(i) Additional production modes (such as direct and multipion production) in photohadronic interactions produce in addition to . The pure source therefore does not exist, and up to 20% contamination from at the source has to be accepted even in this case. As a consequence, the Glashow resonance must be seen to some degree, even for the source.(ii) Since the photohadronic interaction in (2.1) also exists for neutrons producing in this case, any optical thickness to neutron escape will lead to a - equilibration. This means that the optically thick source cannot be distinguished from the source by the Glashow resonance.(iii) Neutron decays, which are inherently present in any photohadronic source, lead to flux, faking a contribution and therefore a source in particular energy intervals (determined by the maximal proton energy).

As a consequence, only the “ optically thin” (to neutron escape) source might be uniquely identified if less than 20% of contamination are found. On the other hand, one cannot uniquely identify a source. Note that modern approaches take the composition between and interactions as a variable, see [90–92]. However, these approaches cannot typically describe the contamination from neutron decays, because this depends very much on the model.

At the detector, the use of the Glashow resonance becomes even more complicated for the following reasons.

(i) Flavor mixing reintroduces a component from produced in decays even for pure production, compare (2.2).(ii) As a consequence, the flavor composition at the source is important and has to be determined at the same time. For example, a muon damped source may be easily mixed up with a pion beam source.(iii) The Glashow resonance occurs at a specific energy (6.3 PeV), which means that for this process the energy dependence is not important and only one particular energy matters. However, the transformation of the energy from source to detector depends on redshift and a possible Lorentz boost, compare (3.11), which have to be known to draw any conclusions.

In summary, the discovery of a neutrino signal with significantly suppressed Glashow resonance events will be interesting and may allow for possible conclusions about the source if flavor composition, , and are known. On the other hand, the detection of a pronounced Glashow resonance events is probably nonconclusive for the physics of the source.

4.4. Testing New Physics in the Neutrino Propagation

In the Standard Model, the transition probability in (3.11) is described by the usual flavor mixing in (4.1), which is independent of energy. New physics effects may lead to deviations from this picture, where the effects discussed in the literature include sterile neutrinos, neutrino decay, quantum decoherence, Lorentz invariance violation, among others, see [53] for a review. In some of these cases, may even be a function of energy . We choose neutrino decay as an example in this section, see [38, 80, 90, 93–100], where in [98–100] energy-dependent effects have been considered.

Following [100], for decay into invisible particles (see, e.g., [90, 101] for the systematical discussion/classification of possible scenarios. In general, neutrinos may decay into particles invisible for the detector, or other visible active flavors. Furthermore, the decay may be complete, i.e., all particles have decayed, or incomplete, where the energy dependence of the decay can be seen. For the case of complete decays, there are only eight different effective decay scenarios (each mass eigenstate can be stable or unstable), no matter if the decays are visible or invisible or if intermediate (unstable) invisible states are involved [90]. Likewise, there are only eight scenarios for invisible incomplete decays, whereas the treatment of incomplete visible decays is more complicated, see, e.g., [101, 102]), the transition probability can be described by a modified version of (4.1):

as the damping coefficient [103]. Here with the rest frame lifetime for mass eigenstate . Typically the neutrino lifetime is quoted as since is unknown, see [53] for a review. Note that this transition probability is energy dependent, compared to (4.1). In this case, the flavor ratio in (4.2) can be rewritten as

if is the ratio between electron neutrinos and muon neutrinos ejected at the source (assuming that hardly any tau neutrinos are produced). Equation (4.5) carries now two energy dependencies: the energy dependence of the flavor composition at the source , and the energy dependence of the new physics effect , which have to be intrinsically disentangled. On the other hand, the energy dependence of the new physics effect may provide a unique signature [98, 99], and the energy dependence of the flavor composition at the source may help to disentangle different scenarios [100].

In general, there are decay scenarios for invisible incomplete decays, since every (active) mass eigenstate may be either stable or unstable. In Figure 9(a), is shown as a function of the initial flavor composition for these different (complete) decay scenarios, where filled disks correspond to stable mass eigenstates and unfilled disks to unstable mass eigenstates (the scenario with only unstable states is of course not shown, since no neutrinos can be detected in this case). In this panel, the different types of sources are also marked. One can clearly see that especially for a pion beam source, three of the scenarios cannot be distinguished, whereas these can be disentangled in principle for any other type of source. In addition, for three scenarios (the scenarios where only one mass eigenstate is stable), is independent of the initial flavor composition. Note that the scenario with the largest ( and unstable) faces the strongest constraints because of the observation of neutrinos from supernova 1987A.

(a)

(b)

Figure 9: (a) as a function of the initial flavor composition for different decay scenarios (complete decays), where filled disks correspond to stable mass eigenstates and unfilled disks to unstable mass eigenstates. (b) as a function of energy for a pion beam to muon damped source (cf. Figure 5, TP 13) for a specific decay parameter (same for all unstable states), as given in the panel. The dashed curves stand for complete decays, figures adopted from [100].

In Figure 9(b), is shown as a function of energy for a particular choice of the decay parameter times distance and a specific source with a flavor transition from pion beam to muon damped source at about  GeV. First of all note that because of the exponential damping in (4.4), the decays are practically complete for energies  GeV, and the neutrinos are stable for energies  GeV for the chosen decay parameters. This is why the scenarios deviate from the complete decay curves (dashed) at very high energies and converge into the standard scenario there (all neutrinos stable). The pion beam (low energies) cannot distinguish three of the scenarios at low energies, as expected. However, above  GeV, where the flavor transition occurs into a muon damped source, the corresponding curves split up (and the neutrino decays are still complete), before they converge into the standard case. This example illustrates that the energy-dependent flavor transition of a specific source may be a useful for new physics tests, provided that enough statistics can be collected.

5. Application to Generic (AGN-Like) Sources

As it is illustrated in Figure 6 (see also (3.1)), the proton and photon densities within the source control the secondary meson, and hence the neutrino production. In this section, we follow the ansatz in [34]: we assume that the target photons are produced in a self-consistent way, by the synchrotron emission of coaccelerated electrons (positrons); see, for example, [104] for a corresponding specific (BL Lac) AGN blazar model. The main purpose of this model is the prediction of spectral shape and flavor composition of the neutrino source, while it cannot predict the normalization of the neutrino flux in the form presented here. In addition, no neutrinos have been observed yet, which means that an ansatz tailor-made for neutrinos may be useful to test the interplay between detector response and source model; after all, there may be sources for which the optical counterpart is absorbed, so-called “hidden sources,” see, for example, [105–108]. The ingredients to this model are comparable to the conventional GRB neutrino models, discussed in Section 6, where however the origin of the target photons is different.

5.1. Additional Model Ingredients

The primaries, in this case protons and electrons, are assumed to be injected with an injection spectrum, where the (universal) injection index is one of the important model parameters. They are assumed to lose energy dominantly by synchrotron radiation, controlled by , and adiabatic cooling, controlled by , which also determine their maximal energies. This means that (3.9) is applied to the primaries, which implies that the acceleration and radiation zones are different in this model. As a consequence, the electron spectrum becomes loss-steepened by one power.

The synchrotron photons (the target photons in (3.1)) are computed in the Melrose-approximation [109] averaged over the pitch angle. The power radiated per photon energy by one particle with energy , mass and charge in a magnetic field is given by

We have to convolute this with the spectrum of radiating electrons as

The number of produced photons per time can be computed with

Assuming that the photons escape (and hardly interact again), the steady photon spectrum, which is needed for the computation of photohadronic interactions, can be estimated by multiplying with escape time of the photons. The synchrotron spectral index obtained from this approach is , which means that the dependence on the primary injection index is small. The rest of the computation follows Section 3.

The main parameters in this model are , , and . For the sake of simplicity, we assume in the following that the source is only moderately Lorentz boosted with respect to the observer’s frame, that is, and . In this case, a convenient description of the parameter space of interest is the Hillas plot [110]. In order to confine a particle in a magnetic field at the source, the Larmor radius has to be smaller than the extension of the acceleration region . This can be translated into the Hillas condition for the maximal particle energy

Here is the charge (number of unit charges) of the accelerated particle, is the magnetic field in Gauss, and can be interpreted as an efficiency factor or linked to the characteristic velocity of the scattering centers. Potential cosmic ray sources are then often shown in a plot as a function of and , as it is illustrated in Figure 10(a), by the numbered disks (see legend for possible source correspondences). Assuming that a source produces the highest energetic cosmic rays with  eV, one can interpret (5.4) as a necessary condition excluding the region below the line in Figure 10 (for protons with ). If one allows for relativistic boosts between source and observer, this condition is relaxed, as one can read off from the figure (in this case, and in the plot have to interpreted as and ). However, this method does not take into account energy loss mechanisms, which lead to a qualitatively different picture, see, for example, [111, 112], and which are implied in our model. In the following, we will study the complete parameter space covered by Figure 10(a) without any prejudice. Since the location of the sources in Figure 10 cannot be taken for granted, we will refer to the individual sources as “test points” (TP), and leave the actual interpretation to the reader.

(a)

(b)

Figure 10: (a) Possible acceleration sites in Hillas plot as a function of and (version adopted from M. Boratav). (b) Classification of sources for injection index in this plot (see main text). Some points from left plot are shown for orientation, as well as two new points (12 and 13), figure taken from [34].

Concerning the limitations of the model, it certainly does not apply exactly to all types of sources. For example, in supernova remnants, (proton-proton) or (proton-nucleus) interactions may dominate the neutrino production, which would require additional parameters to describe the target protons or nucleons. In addition, at ultra-high energies, heavier nuclei may be accelerated. The spirit of this model is different: it is developed as the simplest (minimal) possibility including nontrivial magnetic field and flavor effects. Another ingredient is the target photon density, which is assumed to come from synchrotron emission of coaccelerated electrons here. In more realistic models, typically a combination of different radiation processes is at work. However, in many examples with strong magnetic fields, the specific shape of the photon spectrum is less important for the neutrino spectral shape than the cooling and decay of the secondaries, which depend on particle physics only. To check this, we have tested the hypothesis that acceleration and radiation zones of the electrons are identical (i.e., the electron spectrum is not loss-steepened by synchrotron losses, which is actually a simpler version of this model). Thus, while it is unlikely that the model applies exactly to a particular source, it may be used as a good starting hypothesis.

5.2. Flavor Composition at the Source

We discussed the flavor composition at the source already in Section 3.3, where we also showed several examples for the energy dependence in Figure 5. Let us now approach this from a systematical point of view as a function of and (for ). In this case, a qualitative classification of the sources can be found in Figure 10(b), where it is implied that a certain flavor ratio can be clearly identified over one order of magnitude in energy close enough to the spectral peak; see [34]. From this figure, we find that for  G, all charged species lose energy so rapidly that the neutron decays dominate the neutrino flux. For  G, several processes compete, leading to an undefined source. For  kG, the sources behave as classical pion beams, which typically applies to sources on galactic scales. In the intermediate range, , the source classification somewhat depends on the spectral shape, since this affects possible muon pile-up effects, the energy range close to the spectral peak, and the competition of several effects. Depending on , muon beam to muon damped sources (as a function of energy), muon damped sources, and pion beam to muon damped sources are found. There is some dependence on , which affects the spectral shapes. For instance, for a similar pion spectrum as for GRBs is obtained, which leads to a pion beam to muon damped sources for typical GRB parameters in consistency with [37, 38]. In summary, the pion beam source assumption is safe if  kG, whereas for stronger magnetic fields the magnetic field effects on the secondaries have to be taken into account. However, for given parameters, this flavor composition can be predicted. For instance, may be estimated from the time variability of the source, from energy equipartition arguments, and from the observed photon spectrum.

5.3. Interplay between Spectral Shape and Detector Response

We discussed already earlier in Section 4.2 the interplay between spectral shape and detector response, see [46]. From Figure 8, it is clear that spectral shapes with a peak at the position of the differential limit minimum can be better limited than others. This is especially clear if two event types or detectors are compared for the same spectrum. However, how should one quantify that for different spectra seen by the same detector? Consider, for example, the fluxes no. 2 and no. 11 in Figure 8(a), both leading to the same event rate by definition. Which of the two neutrino sources is the detector more sensitive to in terms of the physics of the source? In order to address this question, it is useful to assign a single number to each spectrum which measures how much energy in neutrinos can be tested for a specific spectrum and event type. We choose the energy flux density

as this quantity, which we show in units of erg cm⁻² s⁻¹ for point sources in order to distinguish it from in units of GeV cm⁻² s⁻¹ ( GeV).

This quantity measures the total energy flux in neutrinos, and it is useful as a performance indicator measuring the efficiency of neutrino production in the source. In order to see that, consider the alternative derivation of from (3.11) (neglecting a possible Lorentz boost of the source)

where is the “neutrino luminosity” and is the volume of the interaction region. Since the neutrinos originate mostly from pion decays and take a certain fraction of the pion energy (about 1/4 per produced neutrino for each charged pion), the neutrino luminosity is directly proportional to the (internal) luminosity of protons (or the proton energy dissipated within a certain time frame ) and the fraction of the proton energy going into pion production, commonly denoted by (if the energy losses of the secondaries can be neglected). Since a possibly emitted photon flux can be often linked to by energy partition arguments, one has , and is a measure for of the source (if no photon counterpart is observed), or even itself (if a photon counterpart is observed).

Regions for a specific sensitivity to are shown as a function of and in Figure 11 for (Figure 11(a)) and (Figure 11(b)), for several source declinations in IceCube and for Earth-skimming tau neutrinos in Auger (in this case, for a diffuse flux). There are several conclusions from this figure: first of all, IceCube responds very well to the usual suspects, such as AGNs (Figure 11(a), points 3 to 6). Even for , which produce a muon track in only 17% of all cases, most of this parameter space can be covered (Figure 11(b)). Thus, it is clear that most sources will be also detected if the partition between and was heavily disturbed, such as by a new physics effect in the neutrino propagation. The sensitive region, however, somewhat depends on the source declination. For very large values of , the neutrino energies are lower and instruments such as the DeepCore array may respond better to the flux (Figure 11). On the other hand, the region where the UHECRs are expected in this model (lower right corner) is better covered by Auger (Figure 11(b)). This is not a big surprise: the neutrino spectrum follows the proton spectrum in the absence of strong magnetic field effects on the secondaries, which means that the spectra (cf. Figure 8) peak at high energies. A very interesting region may be the gap between the IceCube and Auger regions: perhaps a future instrument such as KM3NeT should optimize their geometry to be complementary in terms of energy range coverage.

(a)

(b)

Figure 11: Regions where the sensitivity in exceeds  erg cm⁻² s⁻¹ [sr⁻¹] for (a) and (b) for several selected data samples (90% CL, ). The dashed regions “UHECR” indicate where  eV cosmic ray protons are expected to be produced in the model. The different sources (test points) correspond to the ones in Figure 10, figure taken from [46].

6. GRB Neutrinos and the Multimessenger Connection

Recall again that, as it is illustrated in Figure 6 (cf. (3.1)), the proton and photon densities control neutrino production. In the previous section, we assumed that is produced by synchrotron emission of coaccelerated electrons. Here we emphasize the multimessenger connection, assuming that can be derived from the gamma-ray observation [30, 44]; see also [15, 20, 21, 113] for analytical approaches. The motivation is also different from the previous section: whereas we were interested in a systematic parameter space study of spectral shape and flavor ratio before, the main emphasis here is the prediction of spectral shape, flavor composition, and absolute neutrino flux normalization for a specific set of GRBs observed in gamma rays.

6.1. Additional Model Ingredients

Here we describe the key ingredients of the conventional (numerical) fireball GRB model for neutrino emission, following [30], where we focus on the normalization. For a detailed comparison to the analytical calculations, see [29, 44]. Because the model does not describe the neutrino production in a time-resolved way, it makes sense to relate the neutrino production to the (bolometric) gamma-ray fluence (in units of erg cm⁻²) during the burst. The isotropic bolometric equivalent energy (in erg) in the source (engine) frame can be then obtained as

It can easily be boosted into the SRF by . Assuming energy equipartition between photons and electrons, the photons carry a fraction (fraction of energy in electrons assuming equipartition between kinetic energy in electrons and energy in photons) of the total energy as

In order to compute the photon and proton densities in the SRF, it turns out to be useful to define an “isotropic volume”

where is the collision radius, is the shell thickness of the colliding shells, and is the (observed) variability timescale. Thus, can be estimated from , , and and can be regarded as the volume of the interaction region assuming isotropic emission by the source. Because of the intermittent nature of GRBs, the total fluence is assumed to be coming from such interaction regions, where is the duration of the burst (time during which 90% of the total energy is observed). Now one can determine the normalization of the photon spectrum in (3.1) from

assuming that the spectral shape is determined by the observed spectrum. Similarly, one can compute the normalization of the proton spectrum in (3.1) by

where is the ratio between energy in electrons and protons (: baryonic loading). Note that, in the end, we will obtain the neutrino flux per time frame per interaction region from (3.11) with . Assuming that the magnetic field carries a fraction of , one has in addition

Typical values used in the literature are (see, e.g., [21]). An explicit calculation for yields

In summary, once the photon (from observation) and proton (typically ) spectral shapes are determined, the proton and photon densities in the source and can be calculated with the above formulas from the observables (gamma-ray fluence, , , ). Equation (6.3) implies that for fixed , the larger , the larger the interaction region, and the smaller the photon density in (6.4), which directly enters the fraction of fireball proton energy lost into pion production [20], or, consequently, . Therefore, the main contribution to the neutrino flux is often believed to come from bursts with small Lorentz factors, see discussions in [20, 114, 115]. We follow this conventional fireball approach in the following. Note that this numerical approach contains all the ingredients of [20, 21] explicitly—such as the cooling of the secondaries. The neutrino emission can be then easily computed as shown in Figure 6. In addition, note that there are alternatives to the model. For instance, if the bursts are alike in the comoving frame, as suggested in [116], one has [30]. One can read off from (6.3) that the correlation is expected in that case, since will be similar for the bursts. See [30] for a discussion of the neutrino flux for different model hypotheses, and how the neutrino flux can in principle be used to discriminate among these.

6.2. Systematics in the Interpretation of Aggregated Fluxes

Since the number of neutrinos expected from a single GRB is small, dedicated aggregation methods are needed. For instance, one may search for the diffuse flux from GRBs, which, however, has to fight the background from atmospheric neutrinos. Another possibility is to use the gamma-ray observation to infer on time window and direction of the neutrino signal, which effectively leads to significantly reduced backgrounds. In addition, a procedure such as the one in Section 6.1 may be used to predict the absolute neutrino flux, its shape (see, e.g., Figure 4), or its flavor composition (see, e.g., Figure 7) on a burst-by-burst basis. Summing over many observed bursts, such an analysis is also called stacking analysis, see, for example, [15] for a recent example. The diffuse limit can be extrapolated from such a stacked flux, which is also called the quasi-diffuse limit; see [7] for details. Here we discuss some of the implications when a stacked neutrino flux is translated into a quasi-diffuse flux.

In order to illustrate this problem, following [30], let us consider the redshift distribution of the GRBs as an example. We show in Figure 12(a), a population of 10 000 bursts, representative for the number of GRBs in the visible universe over about 10 years (cf. histogram). These bursts are assumed to follow the star formation rate from Hopkins and Beacom [117] with the correction from Kistler et al. [118]. For the sake of simplicity, assume that all bursts have the same isotropic luminosity. From the discussion after (3.11), we have , which means that closer GRBs, which are however rarer, will lead to a larger neutrino flux. Thus, it is the product which determines the main contribution to the neutrino flux, shown as solid curve in Figure 12. While the peak contribution in terms of the GRB distribution is at z ~ 2–4, this contribution function peaks at . This observation has several implications: first of all, if the redshift is not measured, the neutrino flux may be overestimated if is assumed; compare (6.1). Second, the number of bursts contributing in the region is rather small, which means that large statistical fluctuations are expected in quasi-diffuse flux estimates based on the stacking of a few bursts only.

(a)

(b)

Figure 12: (a) Distribution of 10 000 bursts as a function of redshift (histograms) and relative contribution of the individual GRBs (solid curves). The dashed curve shows the exact distribution function. Here it is assumed that the GRBs follow the star formation rate from Hopkins and Beacom [117] with the correction from Kistler et al. [118]. (b) Probability that the quasi-diffuse flux extrapolated from a low statistics sample with bursts is larger than a certain fraction of the diffuse limit (see legend for different values of ). This function corresponds to (one minus) the cumulative distribution function of the probability density, figure taken from [30].

We quantify this systematical errors in Figure 12(b). In this panel, we show the probability that the quasi-diffuse flux extrapolated from a low statistics sample with bursts is larger than a certain fraction of the diffuse limit (see legend for different values of ). This function corresponds to (one minus) the cumulative distribution function of the probability density, and the step function corresponds to the diffuse limit. One can read off from this figure that for , corresponding to the analysis in [15], the quasi-diffuse extrapolation will be within 50% of the diffuse limit in the probability range corresponding to 90% of all cases (between 0.05 and 0.95). This means that a 50% error on the quasi-diffuse flux can be estimated from the redshift distribution only, while additional parameter variations increase this error [30].

6.3. Neutrino Flux Predictions from Gamma-Ray Observations

Since IceCube has not observed any GRB neutrino flux yet, there has been increasing tension between the model predictions [9, 20, 21] and the observation [15, 19] (recently, the superluminal propagation of neutrinos has also been proposed as a reason why no neutrinos have been seen, see, e.g., [119]). On the other hand, a direct comparison with the photomeson production in [9] was performed in [29] (see also Figure 1), demonstrating that the neutrino flux is actually underestimated in the analytical approaches. Therefore, the differences between the numerical approach in Section 6.1 and the analytical models in [20, 21] (based on [9]) have been identified in [44] by a recomputation of the analytical models and the analytical computation of a simplified version of the numerical code. As far as the astrophysics ingredients are concerned, these approaches can be shown to be equivalent, based on the same logic; see Section 6.1. The main differences are that magnetic field and flavor-dependent effects are explicitly included in numerical approach, additional pion, neutron, and kaon production modes are computed, and the full energy dependencies are taken into account.

We illustrate in Figure 13 the comparison between analytical and numerical approaches, where three different recent Fermi-measured GRBs have been chosen as examples: GRB 080916C, GRB 090902B, and GRB 091024. GRB 080916C has been selected, because it is one of the brightest bursts ever seen, although at a large redshift, and one of the best studied Fermi-LAT bursts. The gamma-ray spectrum of GRB 090902B can be fit by a Band function and a cutoff power law (CPL), which means that it can be used to illustrate the difference. GRB 091024 can be regarded as a typical example representative for many Fermi-GBM bursts [120], except from the long duration. Note that the first two bursts have an exceptionally large , whereas for the third burst. All three bursts have in common that the required parameters for the neutrino flux computation can be taken from the literature, in particular, the properties of the gamma-ray spectrum (including fluence), , , , and , see figure caption for the references.

Figure 13: Expected (time-integrated) neutrino flux for three different GRBs (parameters: see [120, 123] for GRB 080916C, [120, 124] for GRB 090902B, and [120, 125] for GRB 091024). First row: revision of IceCube analytical method conventional fireball calculation “CFB” (correction of shape, normalization from pion production efficiency , and normalization from neutrino versus proton spectral shape ), leading to the revised fireball calculation “RFB.” Second row: comparison analytical (CFB, RFB) methods with simplified numerical method “WB Δ-approximation” (3.4) and full interactions. Third row: impact of adiabatic cooling on protons and secondaries. Fourth and fifth row: corresponding inverse timescales (rates) for protons and secondaries, respectively, courtesy of Hümmer [126].

In Figure 13, we show in the first row the computation of the predicted neutrino flux with the IceCube analytical method [21], called CFB (conventional fireball calculation) here. As described in [44], the corrections of shape, normalization from pion production efficiency , and normalization from neutrino versus proton spectral shape (see also [45]) lead to a revised analytical calculation RFB. From the figure, one can easily read off that these revisions strongly depend on the burst parameters, especially the photon spectral shape. Comparing the predicted CFB fluxes with the bursts used for the IC40 analysis (Figure 1 in [15]), one can easily see that the expected fluxes of the first two bursts are about a factor 50 below that of the most luminous bursts in that analysis, and the third example about a factor of 5 below. This is expected from the scaling of the pion production efficiency in that approach.

In the second row of Figure 13, we show the comparison of the analytical (CFB, RFB) methods with a simplified numerical method “WB Δ-approximation”, compare (3.4), and the full interactions. In most cases, the simplified numerical approach matches the method RFB rather well, which proves the validity of the derived corrections. However, for GRB 090902B (middle panel), the spectrum below the first break is different from the analytical estimate because the scalings of the weak decays limit the steepness of the spectrum there. The final numerical calculation including all production modes is then significantly enhanced again, especially due to multipion production [29]. Using a cutoff power law for GRB 090902B for the gamma-ray spectral fit (middle panel), the normalization of the prediction slightly reduces in that example because the photon density above the photon break is suppressed. Comparing the original CFB method with the final numerical computation “full ”, it is interesting that this can significantly deviate in both normalization and shape, but this deviation depends on the burst parameters. Very interestingly, a similar neutrino flux normalization for all three bursts is obtained, which means that it is probably not warranted to say that the neutrino flux from high- bursts such as GRB 080916C is expected to be small. However, note that such extreme bursts only make up for a small fraction of the observed bursts, and the conclusions from neutrinos will be determined by the statistical properties of the burst sample in the stacking analysis. Finally, note that the numerical calculations in Figure 13 do not depend on any approximations, whereas different analytical methods lead to different predictions, similar to CFB. Therefore, the numerical computations should be regarded as the benchmark which defines the corrections, not vice versa. Within the simplest fireball neutrino model, there are only small model dependencies within the numerical approach. For instance, the integral limits in (6.5) (the minimal and maximal proton energies) have to be specified (note, however, that for an injection spectrum, the energy partition only logarithmically depends on the minimal and maximal proton energies), and a bolometric correction may have to be applied to (6.4)—which typically has small effects.

One of the effects included in the final result of the second row of Figure 13 is the adiabatic cooling of primaries (protons) and secondaries (muons, pions, kaons). We illustrate the effect of this cooling component in the lower three rows of the figure, where we show the impact on the final result (third row) and the respective inverse timescales (rates) for protons (fourth row) and secondaries (fifth row). For the protons (fourth row), synchrotron losses are assumed to determine the maximal energies in the absence of adiabatic cooling, whereas the larger of the synchrotron or adiabatic cooling loss rates determines the maximal proton energy otherwise. From the fourth row, one can also easily read off that the energy losses due to interactions are typically subdominant. The comparison with the third row illustrates that, depending on the burst parameters, the kaon hump (the rightmost one) can be suppressed by the adiabatic cooling of the protons, whereas the normalization of the spectra is hardly affected. From the fifth row, one can also read off that adiabatic cooling may have a small effect on the muons, for which it sometimes dominates in a small energy range. This leads to a small suppression of the first hump, coming from the muon decays (see left and middle panels of third row).

We show in Figure 14 the predicted quasi-diffuse neutrino flux from the above numerical method applied to the IC40 bursts for the same bursts and parameters (solid black curve) used in that analysis, which is about one order of magnitude lower than the IC40 limit and a factor of two below the current limit. In this figure a number of systematical errors are shown as well, such as the statistical error discussed in the previous subsection, and the estimated astrophysical uncertainty (by varying the unknown parameters, such as proton injection index , variability timescale by one order of magnitude around the IceCube standard values,  s for long bursts, from 200 to 500, and the ratio from 0.1 to 10). In addition, note that has only been measured for a few bursts used in the IC40 analysis, whereas has been assumed for the long bursts with unknown . As we illustrated in the previous subsection, this is potentially problematic, which means that a solid lower limit for the prediction can be only obtained for bursts with measured (dashed black curve in Figure 14). Note that our prediction varies not as strongly as one may expect from [20]. The synchrotron losses of the secondaries damp this variation [115]: for smaller , the energy densities in the source increase because of the energy equipartition, and consequently in (6.6). This increases the energy losses of the secondaries, which means that less energy goes into the neutrinos.

Figure 14: Numerical prediction (NFC) of the quasi-diffuse flux expected from the bursts used in the IC40-analysis, including the estimates for several model- or method-specific systematical uncertainties (see main text). In addition, the IC40 limit is shown, and two expectations are shown for comparison (IC59+40 from [19] and IC86 extrapolated for from IC40; see, e.g., [87]), figure taken from [44].

From the extrapolated IceCube limit for IC86 in Figure 14, it is obvious that IceCube will finally test the parameter space of the GRB fireball models, whereas the current limits already start to enter the meaningful parameter space. On the other hand, it is not clear what values of and the dominant bursts actually have, and what would contribute most to the neutrino flux. For example, a theoretical study of the dominant to the diffuse neutrino flux has been performed in [30] for different model hypotheses, with a result strongly depending on the model and underlying distribution. In the conventional fireball approach, which is presented above, bursts with lower than (as used by IceCube) tend to dominate the expected neutrino flux, such as GRB 091024. In this case, one would expect the prediction rather in the upper half of the shaded region of Figure 14, and a near-future detection in IceCube may be rather likely. On the other hand, recent studies for seem to prefer larger average  s than the  s used by IceCube for long bursts [121], which points towards the lower half of the shaded region. Note that in practice, bursts with very large and (or collision radii) will hardly contribute to the quasi-diffuse limit, which means that the expectation depends on the (weighted) population of bursts with small and compared to the total population. Independent of these parameters, the neutrino flux prediction is proportional to , which means that the final IceCube results will limit the baryonic loading as key parameter.

Finally, note that similar considerations as for the GRBs may apply to AGN models, such as [10, 12, 122], and it is yet to be seen what the impact on the cosmic ray connection is, see [22].

7. Summary and Conclusions

We have discussed the impact of particle physics effects on the neutrino fluxes from cosmic accelerators, based on a generic numerical approach of the photohadronic charged meson production in sources optically thin to neutron escape; compare Figure 6. The starting point has been the proton and photon densities (spectra) within the source, which lead to the secondary meson production. We have included additional -channel, higher resonance, and multipion production modes, neutrinos from kaon and neutron decays, the synchrotron cooling of the secondaries (pions, muons, and kaons), the helicity dependence of the muon decays, and the full spectral dependencies. While this approach may not be able to describe every source accurately, it can be regarded as the minimal approach to neutrino production by interactions using at least the well-known particle physics ingredients, to be extended by further components if applied to particular source classes. We have applied this framework to GRBs in Section 6.1, where the target photons are inferred from the gamma-ray observation. We have demonstrated that it contains all the ingredients in frequently-used analytical models, such as [9, 15, 20, 21], which means that it can be regarded as a numerical refined version of the conventional fireball neutrino model. In addition, we have discussed a generic AGN-like application where the target photons are produced by the synchrotron emission of coaccelerated electrons in Section 5. From the presentation in this paper, it must be clear that this application is at a similar level as the GRB approach in terms of the ingredients. Applications to specific neutrino flux normalization predictions are, however, yet to be performed.

The main consequences of the particle physics effects discussed in this paper can be summarized as follows.

(i) Magnetic field effects and flavor mixing change the flavor composition. The assumption of flavor composition at the source of 1 : 2 : 0 (pion beam) can only be justified for  kG, while it will normally change as a function of energy for larger values of .(ii) The neutrino spectral shape for interactions depends on proton and photon injection spectra and the magnetic field and flavor effects. Roughly speaking, the assumption only holds for the special case for the photon spectrum. This means that this assumption is too simple for many sources, since the detector response depends on the interplay between spectral shape and differential limit. For example, the detector may miss a neutrino signal because the spectrum does not peak at the right energy, although the fraction of energy going into pion production in the source is the same as in another detectable case.(iii) Additional pion production processes increase the neutrino production significantly and lead to an impact on the spectral shape. Approximations using the Δ-resonance approximation are useful for analytical estimates, but they do not take into account the cross section dependence at high center-of-mass energies. The impact on the neutrino flux is a factor of a few.(iv) Spectral effects, such as the energy dependence of the proton interaction length or the energy dependence of the photon spectrum, leads a significant reduction of the GRB neutrino flux prediction of about one order of magnitude in comparison to analytical estimates. As a consequence, IceCube has just entered the relevant part of the parameter space to test the simplest GRB fireball model. Similar effects on AGN flux models are yet to be tested.(v) Neutrino fluxes from kaon and neutron decays are generically expected in any source. However, the neutrinos from neutron decays show typically up below the peak, the neutrinos from kaon decays above the peak. The relative importance of these effects depends on the strength of (which may separate the spectra from the different parents) and the maximal proton energy (which controls the neutron decay contribution).(vi) Using the Glashow resonance for the discrimination of and interactions in the source is challenging. Especially, it is very hard to infer any physics information from the discovery of Glashow resonant events, whereas the absence of Glashow resonant events (in presence of a neutrino signal) can be interpreted as a source optically thin to neutron escape.(vii) Already from the current IceCube limits one can extrapolate that conclusions for astrophysical sources will most likely be based on source classes, not individual sources. The aggregation of fluxes from many sources, however, introduces new systematical errors. For instance, the extrapolation from the neutrino fluxes from 117 GRBs to a quasi-diffuse flux implies an error of at least 50% (at the 90% confidence level), from the redshift distribution only.(viii) In the presence of a high statistics neutrino flux (close to the current bound), useful information on possible new physics effects in the neutrino propagation may be inferred from the ratio between muon tracks and cascades. In this case, the energy dependence of the flavor ratio at the source and the energy dependence of the new physics effects may provide the necessary information to identify the effect.

In conclusion, in the presence of data, neutrino astrophysics is now at the point where rough analytical estimates for neutrino fluxes are not sufficient anymore, since in that case no reliable conclusions can be drawn for the astrophysical models and thus the origin of the cosmic rays. Especially the particle physics effects discussed in this paper can be relatively easily taken into account, and they are well known. We have demonstrated with numerous examples that these effects on flux normalization, spectral shape, flavor composition, and neutrino-antineutrino composition cannot be neglected. For example, a correction of the GRB neutrino flux prediction of one order of magnitude has been identified. Therefore, a reliable treatment of the particle physics of the neutrino production should be the state-of-the-art of any neutrino data interpretation in the multimessenger context.

Acknowledgments

The author would like to thank Philipp Baerwald, Mauricio Bustamante, Svenja Hümmer, and Guey-Lin Lin for useful comments, and Svenja Hümmer for providing some figure from her Ph.D. thesis. This publication was funded by the German Research Foundation (DFG) and the University of Würzburg in the funding program Open Access Publishing. It has been supported by Deutsche Forschungsgemeinschaft, Grants WI 2639/3-1 and WI 2639/4-1.

References

J. Ahrens, J. N. Bahcall, X. Bai et al., “Sensitivity of the IceCube detector to astrophysical sources of high energy muon neutrinos,” Astroparticle Physics, vol. 20, no. 5, pp. 507–532, 2004.

ANTARES, E. Aslanides, et al., “A deep sea telescope for high energy neutrinos,” http://www.arxiv.org/abs/astro-ph/9907432/

J. G. Learned and K. Mannheim, “High-energy neutrino astrophysics,” Annual Review of Nuclear and Particle Science, vol. 50, no. 1, pp. 679–749, 2000.

F. Halzen and D. Hooper, “High-energy neutrino astronomy: the cosmic ray connection,” Reports on Progress in Physics, vol. 65, no. 7, pp. 1025–1078, 2002.

T. Chiarusi and M. Spurio, “High-energy astrophysics with neutrino telescopes,” European Physical Journal C, vol. 65, no. 3, pp. 649–701, 2009.

U.F. Katz and C. Spiering, “High-energy neutrino astrophysics: status and perspectives,” Progress in Particle and Nuclear Physics, vol. 67, no. 3, pp. 651–704, 2012.

J.K. Becker, “High-energy neutrinos in the context of multimessenger astrophysics,” Physics Reports, vol. 458, no. 4-5, pp. 173–246, 2008.

J. P. Rachen and P. Mészáros, “Photohadronic neutrinos from transients in astrophysical sources,” Physical Review D, vol. 58, no. 12, article 123005, 1998.

E. Waxman and J. Bahcall, “High energy neutrinos from cosmological gamma-ray burst fireballs,” Physical Review Letters, vol. 78, no. 12, pp. 2292–2295, 1997.

F. W. Stecker, C. Done, M. H. Salamon, and P. Sommers, “High-energy neutrinos from active galactic nuclei,” Physical Review Letters, vol. 66, no. 21, pp. 2697–2700, 1991.

K. Mannheim, “The proton blazar,” Astronomy and Astrophysics, vol. 269, no. 1-2, pp. 67–76, 1993.

F. W. Stecker, “Note on high-energy neutrinos from active galactic nuclei cores,” Physical Review D, vol. 72, no. 10, Article ID 107301, 2005.

R. Abbasi, Y. Abdou, T. Abu-Zayyad et al., “Time-integrated searches for point-like sources of neutrinos with the 40-string IceCube detector,” Astrophysical Journal Letters, vol. 732, no. 1, article 18, 2011.

R. Abbasi, Y. Abdou, T. Abu-Zayyad et al., “Time-dependent searches for point sources of neutrinos with the 40-string and 22-string configurations of icecube,” Astrophysical Journal, vol. 744, no. 1, article 1, 2012.

R. Abbasi, Y. Abdou, T. Abu-Zayyad et al., “Limits on neutrino emission from gamma-ray bursts with the 40 string IceCube detector,” Physical Review Letters, vol. 106, no. 14, Article ID 141101, 2011.

R. Abbasi, Y. Abdou, T. Abu-Zayyad et al., “Search for a diffuse flux of astrophysical muon neutrinos with the IceCube 40-string detector,” Physical Review D, vol. 84, no. 8, Article ID 082001, 2011.

E. Waxman and J. Bahcall, “High energy neutrinos from astrophysical sources: an upper bound,” Physical Review D, vol. 59, no. 2, pp. 1–8, 1999.

K. Mannheim, R.J. Protheroe, and J.P. Rachen, “Cosmic ray bound for models of extragalactic neutrino production,” Physical Review D, vol. 63, no. 2, Article ID 023003, 2001.

IceCube, R. Abbasi, et al., “The IceCube neutrino observatory I: point source searches,” in Proceedings of the 32nd International Cosmic Ray Conference (ICRC'11), Beijing, China, 2011.

D. Guetta, D. Hooper, J. Alvarez-Muñiz, F. Halzen, and E. Reuveni, “Neutrinos from individual gamma-ray bursts in the BATSE catalog,” Astroparticle Physics, vol. 20, no. 4, pp. 429–455, 2004.

R. Abbasi, Y. Abdou, T. Abu-Zayyad et al., “Search for muon neutrinos from gamma-ray bursts with the IceCube neutrino telescope,” Astrophysical Journal Letters, vol. 710, no. 1, pp. 346–359, 2010.

M. Ahlers, M. C. Gonzalez-Garcia, and F. Halzen, “GRBs on probation: testing the UHE CR paradigm with IceCube,” Astroparticle Physics, vol. 35, no. 2, pp. 87–94, 2011.

C. A. Argüelles, M. Bustamante, and A. M. Gago, “IceCube expectations for two high-energy neutrino production models at active galactic nuclei,” Journal of Cosmology and Astroparticle Physics, vol. 2010, no. 12, Article ID 005, 2010.

A. Mücke, J. P. Rachen, R. Engel, R. J. Protheroe, and T. Stanev, “On photohadronic processes in astrophysical environments,” Publications of the Astronomical Society of Australia, vol. 16, article 160, 1999.

A. Mücke, R. Engel, J. P. Rachen, R. J. Protheroe, and T. Stanev, “Monte Carlo simulations of photohadronic processes in astrophysics,” Computer Physics Communications, vol. 124, no. 2-3, pp. 290–314, 2000.

A. Mücke and R. J. Protheroe, “A proton synchrotron blazar model for flaring in Markarian ~501,” Astroparticle Physics, vol. 15, no. 1, pp. 121–136, 2001.

K. Murase and S. Nagataki, “High energy neutrino emission and neutrino background from gamma-ray bursts in the internal shock model,” Physical Review D, vol. 73, no. 6, Article ID 063002, 2006.

S. Hümmer, M. Rüger, F. Spanier, and W. Winter, “The Astrophysical Journal,” 721, article 630, p. 2010.

P. Baerwald, S. Hümmer, and W. Winter, “Magnetic field and flavor effects on the Gamma-ray burst neutrino flux,” Physical Review D, vol. 83, no. 6, Article ID 067303, 2011.

P. Baerwald, S. Hümmer, and W. Winter, “Systematics in the interpretation of aggregated neutrino flux limits and flavor ratios from gamma-ray bursts,” Astroparticle Physics, vol. 35, no. 8, pp. 508–529, 2012.

K. Asano and S. Nagataki, “Very high energy neutrinos originating from kaons in Gamma-ray bursts,” Astrophysical Journal Letters, vol. 640, no. 1, pp. L9–L12, 2006.

M. Kachelrieß and R. Tomàs, “High-energy neutrino yields from astrophysical sources: weakly magnetized sources,” Physical Review D, vol. 74, no. 6, Article ID 063009, 2006.

M. Kachelrieß, S. Ostapchenko, and R. Tomàs, “High energy neutrino yields from astrophysical sources. II. Magnetized sources,” Physical Review D, vol. 77, no. 2, Article ID 023007, 2008.

S. Hümmer, M. Maltoni, W. Winter, and C. Yaguna, “Energy dependent neutrino flavor ratios from cosmic accelerators on the Hillas plot,” Astroparticle Physics, vol. 34, no. 4, pp. 205–224, 2010.

R. Moharana and N. Gupta, “Tracing cosmic accelerators with decaying neutrons,” Physical Review D, vol. 82, no. 2, Article ID 023003, 2010.

R. Moharana and N. Gupta, “Neutrinos from decaying muons, pions, kaons and neutrons in gamma ray bursts,” in press, http://www.arxiv.org/abs/1107.4483/

T. Kashti and E. Waxman, “Astrophysical neutrinos: flavor ratios depend on energy,” Physical Review Letters, vol. 95, no. 18, Article ID 181101, 2005.

P. Lipari, M. Lusignoli, and D. Meloni, “Flavor composition and energy spectrum of astrophysical neutrinos,” Physical Review D, vol. 75, no. 12, Article ID 123005, 2007.

M. M. Reynoso and G. E. Romero, “Magnetic field effects on neutrino production in microquasars,” Astronomy and Astrophysics, vol. 493, no. 1, pp. 1–11, 2009.

S. Pakvasa, “Neutrino flavor goniometry by high energy astrophysical beams,” Modern Physics Letters A, vol. 23, no. 17–20, pp. 1313–1324, 2008.

S. Barr, T. K. Gaisser, P. Lipari, and S. Tilav, “Ratio of ve vμ in atmospheric neutrinos,” Physics Letters B, vol. 214, no. 1, pp. 147–150, 1988.

G. Barr, T. K. Gaisser, and T. Stanev, “Flux of atmospheric neutrinos,” Physical Review D, vol. 39, no. 11, pp. 3532–3534, 1989.

P. Lipari, “Lepton spectra in the earth's atmosphere,” Astroparticle Physics, vol. 1, no. 2, pp. 195–227, 1993.

S. Hummer, P. Baerwald, and W. Winter, “Neutrino emission from gamma-ray burst fireballs, revised,” Physical Review Letters, vol. 108, no. 23, 5 pages, 2012.

Z. Li, “Note on the normalization of predicted gamma-ray burst neutrino flux,” Physical Review D, vol. 85, no. 2, Article ID 027301, 2012.

W. Winter, “Interpretation of neutrino flux limits from neutrino telescopes on the Hillas plot,” Physical Review D, vol. 85, no. 2, Article ID 023013, 2012.

S. R. Kelner, F. A. Aharonian, and V. V. Bugayov, “Energy spectra of gamma rays, electrons, and neutrinos produced at proton-proton interactions in the very high energy regime,” Physical Review D, vol. 74, no. 3, Article ID 034018, 2006.

F. Vissani and F. Aharonian, “Galactic sources of high-energy neutrinos: highlights,” In press, http://www.arxiv.org/abs/1112.3911/

Pierre Auger Observatory Collaboration, J. Abraham, et al., “Measurement of the depth of maximum of extensive air showers above 10¹⁸ eV,” Physical Review Letters, vol. 104, Article ID 091101, 7 pages, 2010.

K. Greisen, “End to the cosmic-ray spectrum?” Physical Review Letters, vol. 16, no. 17, pp. 748–750, 1966.

G. Zatsepin and V. Kuzmin, “Upper limit of the spectrum of cosmic rays,” Journal of Experimental and Theoretical Physics Letters, vol. 4, no. 3, p. 78, 1966.

S. R. Kelner and F. A. Aharonian, “Energy spectra of gamma rays, electrons, and neutrinos produced at interactions of relativistic protons with low energy radiation,” Physical Review D, vol. 78, no. 3, Article ID 034013, 2008.

S. Pakvasa, “Neutrino flavor detection at neutrino telescopes and its uses,” In press, http://www.arxiv.org/abs/1004.5413/

H. B. Koers and R. A. Wijers, “Enhanced high-energy neutrino emission from choked gamma-ray bursts due to meson and muon acceleration,” In press, http://www.arxiv.org/abs/0711.4791/

Y. Farzan and A.Yu. Smirnov, “Coherence and oscillations of cosmic neutrinos,” Nuclear Physics B, vol. 805, no. 1-2, pp. 356–376, 2008.

T. Schwetz, M. Tortola, and J.W.F. Valle, “Where we are on theta;13: addendum to 'Global neutrino data and recent reactor fluxes: Status of three-flavor oscillation parameters',” New Journal of Physics, vol. 13, Article ID 109401, 2011.

J.G. Learned and S. Pakvasa, “Detecting ν T oscillations at PeV energies,” Astroparticle Physics, vol. 3, no. 3, pp. 267–274, 1995.

W. Rodejohann, “Neutrino mixing and neutrino telescopes,” Journal of Cosmology and Astroparticle Physics, no. 1, Article ID 029, 2007.

P.D. Serpico and M. Kachelrieß, “Measuring the 13 neutrino mixing angle and the CP phase with neutrino telescopes,” Physical Review Letters, vol. 94, no. 21, Article ID 211102, 2005.

P.D. Serpico, “Probing the 2-3 leptonic mixing at high-energy neutrino telescopes,” Physical Review D, vol. 73, no. 4, Article ID 047301, 2006.

W. Winter, “How astrophysical neutrino sources could be used for early measurements of neutrino mass hierarchy and leptonic CP phase,” Physical Review D, vol. 74, no. 3, Article ID 033015, 2006.

Z.-Z. Xing, “Neutrino telescopes as a probe of broken μ-τ symmetry,” Physical Review D, vol. 74, no. 1, Article ID 013009, 2006.

D. Majumdar and A. Ghosal, “Probing deviations from tribimaximal mixing through ultrahigh energy neutrino signals,” Physical Review D, vol. 75, no. 11, Article ID 113004, 2007.

K. Blum, Y. Nir, and E. Waxman, “Probing CP violation in neutrino oscillations with neutrino telescopes,” In press, http://www.arxiv.org/abs/0706.2070/

R.L. Awasthi and S. Choubey, “Confusing sterile neutrinos with deviation from tribimaximal mixing at neutrino telescopes,” Physical Review D—Particles, Fields, Gravitation and Cosmology, vol. 76, no. 11, Article ID 113002, 2007.

G.-R. Hwang and K. Siyeon, “Conflict between the identification of a cosmic neutrino source and the sensitivity to mixing angles in a neutrino telescope,” Physical Review D, vol. 78, no. 9, Article ID 093008, 2008.

S. Pakvasa, W. Rodejohann, and T. J. Weiler, “Flavor ratios of astrophysical neutrinos: implications for precision measurements,” Journal of High Energy Physics, vol. 2008, no. 2, Article ID 005, 2008.

A. Donini and O. Yasuda, “Signatures of sterile neutrino mixing in high-energy cosmic neutrino flux,” In press, http://www.arxiv.org/abs/0806.3029/

C. Quigg, “Cosmic neutrinos,” In press, http://www.arxiv.org/abs/0802.0013/

S. Choubey, V. Niro, and W. Rodejohann, “Probing θ ₂₃ in neutrino telescopes,” Physical Review D, vol. 77, no. 11, Article ID 113006, 2008.

Z.-Z. Xing and S. Zhou, “Implications of leptonic unitarity violation at neutrino telescopes,” Physics Letters. Section B, vol. 666, no. 2, pp. 166–172, 2008.

J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa, and T. J. Weiler, “Measuring flavor ratios of high-energyastrophysical neutrinos,” Physical Review D, vol. 68, no. 9, Article ID 093005, 2003.

J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa, and T.J. Weiler, “Erratum: Measuring flavor ratios of high-energy astrophysical neutrinos,” Physical Review D, vol. 72, no. 1, Article ID 019901, pp. 1–3, 2005.

R. Abbasi, Y. Abdou, T. Abu-Zayyad et al., “First search for atmospheric and extraterrestrial neutrino-induced cascades with the IceCube detector,” Physical Review D, vol. 84, no. 7, Article ID 072001, 2011.

Z.-Z. Xing and S. Zhou, “Determination of the initial flavor composition of ultrahigh-energy neutrino fluxes with neutrino telescopes,” Physical Review D, vol. 74, no. 1, Article ID 013010, 8 pages, 2006.

K.-C. Lai, G.-L. Lin, and T. C. Liu, “Determination of the neutrino flavor ratio at the astrophysical source,” Physical Review D, vol. 80, no. 10, Article ID 103005, 2009.

S. Choubey and W. Rodejohann, “Flavor composition of ultrahigh energy neutrinos at source and at neutrino telescopes,” Physical Review D, vol. 80, no. 11, Article ID 113006, 2009.

A. Esmaili and Y. Farzan, “An analysis of cosmic neutrinos: flavor composition at source and neutrino mixing parameters,” Nuclear Physics B, vol. 821, no. 1-2, pp. 197–214, 2009.

K.-C. Lai, G.-L. Lin, and T. C. Liu, “Flavor transition mechanisms of propagating astrophysical neutrinos: a model independent parametrization,” Physical Review D, vol. 82, no. 10, Article ID 103003, 2010.

M. Bustamante, A. M. Gagoa, and C. Peña-Garayc, “Energy-independent new physics in the flavour ratios of high-energy astrophysical neutrinos,” Journal of High Energy Physics, vol. 2010, no. 4, Article ID 066, 2010.

T. Schwetz, M. Tórtola, and J. W. F. Valle, “Three-flavour neutrino oscillation update,” New Journal of Physics, vol. 10, Article ID 113011, 2008.

P. Huber, M. Lindner, T. Schwetz, and W. Winter, “First hint for CP violation in neutrino oscillations from upcoming superbeam and reactor experiments,” Journal of High Energy Physics, vol. 2009, no. 11, Article ID 044, 2009.

P. Huber and W. Winter, “Neutrino factories and the ‘magic’ baseline,” Physical Review D, vol. 68, no. 3, Article ID 037301, 2003.

J. Tang and W. Winter, “Physics with near detectors at a neutrino factory,” Physical Review D, vol. 80, no. 5, Article ID 053001, 2009.

G. J. Feldman and R. D. Cousins, “Unified approach to the classical statistical analysis of small signals,” Physical Review D, vol. 57, no. 7, pp. 3873–3889, 1998.

J. Abraham, P. Abreu, M. Aglietta et al., “Limit on the diffuse flux of ultrahigh energy tau neutrinos with the surface detector of the Pierre Auger Observatory,” Physical Review D, vol. 79, no. 10, Article ID 102001, 2009.

A. Karle and F. T. I. Collaboration, “IceCube,” in Proceedings of the 31st International Cosmic Ray Conference (ICRC'10), Lodz, Poland, 2010.

L. Anchordoqui, H. Goldberg, F. Halzen, and T. J. Weiler, “Neutrinos as a diagnostic of high energy astrophysical processes,” Physics Letters. Section B, vol. 621, no. 1-2, pp. 18–21, 2005.

P. Bhattacharjee and N. Gupta, “Probing neutrino mixing angles with ultrahigh energy neutrino telescopes,” in press, http://www.arxiv.org/abs/hep-ph/0501191/

M. Maltoni and W. Winter, “Testing neutrino flavor mixing plus decay with neutrino telescopes,” Journal of High Energy Physics, vol. 2008, no. 7, Article ID 064, 2008.

Z.-Z. Xing and S. Zhou, “Glashow resonance as a discriminator of ultrahigh-energy cosmic neutrinos originating from pγ and pp collisions,” Physical Review D, vol. 84, no. 3, Article ID 033006, 2011.

A. Bhattacharya, R. Gandhi, W. Rodejohann, and A. Watanabe, “The Glashow resonance at IceCube: signatures, event rates and pp vs. pγ interactions,” Journal of Cosmology and Astroparticle Physics, vol. 2011, no. 10, Article ID 017, 2011.

Y. Farzan and A.Yu. Smirnov, “Leptonic unitarity triangle and CP violation,” Physical Review D, vol. 65, no. 11, Article ID 113001, 2002.

J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa, and T. J. Weiler, “Decay of high-energy astrophysical neutrinos,” Physical Review Letters, vol. 90, no. 18, Article ID 181301, 4 pages, 2003.

J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa, and T. J. Weiler, “Sensitivity to θ ₁₃ and δ in the decaying astrophysical neutrino scenario,” Physical Review D, vol. 69, no. 1, Article ID 017303, 2004.

D. Meloni and T. Ohlsson, “Neutrino flux ratios at neutrino telescopes: the role of uncertainties of neutrino mixing parameters and applications to neutrino decay,” Physical Review D, vol. 75, no. 12, Article ID 125017, 2007.

D. Majumdar, “Unparticle decay of neutrinos and it's effect on ultra high energy neutrinos,” In press, http://www.arxiv.org/abs/0708.3485/

A. Bhattacharya, S. Choubey, R. Gandhi, and A. Watanabe, “Diffuse ultra-high energy neutrino fluxes and physics beyond the standard model,” Physics Letters, Section B, vol. 690, no. 1, pp. 42–47, 2010.

A. Bhattacharya, S. Choubey, R. Gandhi, and A. Watanabe, “Ultra-high energy neutrino fluxes as a probe for non-standard physics,” Journal of Cosmology and Astroparticle Physics, vol. 2010, no. 9, Article ID 009, 2010.

P. Mehta and W. Winter, “Interplay of energy dependent astrophysical neutrino flavor ratios and new physics effects,” Journal of Cosmology and Astroparticle Physics, vol. 2011, no. 3, Article ID 041, 2011.

M. Lindner, T. Ohlsson, and W. Winter, “A combined treatment of neutrino decay and neutrino oscillations,” Nuclear Physics B, vol. 607, no. 1-2, pp. 326–354, 2001.

M. Lindner, T. Ohlsson, and W. Winter, “Decays of supernova neutrinos,” Nuclear Physics B, vol. 622, no. 1-2, pp. 429–456, 2002.

M. Blennow, T. Ohlsson, and W. Winter, “Damping signatures in future neutrino oscillation experiments,” Journal of High Energy Physics, no. 6, pp. 1095–1125, 2005.

A. Mücke, R. J. Protheroe, R. Engel, J. P. Rachen, and T. Stanev, “BL Lac objects in the synchrotron proton blazar model,” Astroparticle Physics, vol. 18, no. 6, pp. 593–613, 2003.

S. Razzaque, P. Mészáros, and E. Waxman, “TeV neutrinos from core collapse supernovae and hypernovae,” Physical Review Letters, vol. 93, no. 18, Article ID 181101, 4 pages, 2004.

S. Ando and J. F. Beacom, “Revealing the supernova-gamma-ray burst connection with TeV neutrinos,” Physical Review Letters, vol. 95, no. 6, Article ID 061103, 2005.

S. Razzaque, P. Mészáros, and E. Waxman, “High energy neutrinos from a slow jet model of core collapse supernovae,” Modern Physics Letters A, vol. 20, no. 31, pp. 2351–2367, 2005.

S. Razzaque and A.Yu. Smirnov, “Flavor conversion of cosmic neutrinos from hidden jets,” Journal of High Energy Physics, vol. 2010, no. 3, Article ID 031, 2010.

D.B. Melrose, Plasma Astrophysics: Nonthermal Processes in Diffuse Magnetized Plasmas, vol. 1 of The Emission, Absorption and Transfer of Waves in Plasmas, Gordon and Breach, 1980.

A. M. Hillas, “The origin of ultra-high-energy cosmic rays,” Annual Review of Astronomy and Astrophysics, vol. 22, pp. 425–444, 1984.

M.V. Medvedev, “Constraint on electromagnetic acceleration of highest energy cosmic rays,” Physical Review E, vol. 67, no. 4 2, Article ID 045401, pp. 454011–454014, 2003.

R.J. Protheroe, “Effect of energy losses and interactions during diffusive shock acceleration: Applications to SNR, AGN and UHE cosmic rays,” Astroparticle Physics, vol. 21, no. 4, pp. 415–431, 2004.

J.K. Becker, M. Stamatikos, F. Halzen, and W. Rhode, “Coincident GRB neutrino flux predictions: implications for experimental UHE neutrino physics,” Astroparticle Physics, vol. 25, no. 2, pp. 118–128, 2006.

D. Guetta, M. Spada, and E. Waxman, “Efficiency and spectrum of internal gamma-ray burst shocks,” Astrophysical Journal Letters, vol. 557, no. 1, pp. 399–407, 2001.

D. Guetta, M. Spada, and E. Waxman, “On the neutrino flux from gamma-ray bursts,” Astrophysical Journal Letters, vol. 559, no. 1 PART 1, pp. 101–109, 2001.

G. Ghirlanda, L. Nava, G. Ghisellini et al., “Gamma-ray bursts in the comoving frame,” Monthly Notices of the Royal Astronomical Society, vol. 420, no. 1, pp. 483–494, 2012.

A. M. Hopkins and J. F. Beacom, “On the normalization of the cosmic star formation history,” Astrophysical Journal Letters, vol. 651, no. 1 I, pp. 142–154, 2006.

M. D. Kistler, H. Yüksel, J. F. Beacom, A. M. Hopkins, and J. S. B. Wyithe, “The star formation rate in the reionization era as indicated by gamma-ray bursts,” Astrophysical Journal Letters, vol. 705, no. 2 PART 2, pp. L104–L108, 2009.

D. Autiero, P. Migliozzi, and A. Russo, “The neutrino velocity anomaly as an explanation of the missing observation of neutrinos in coincidence with GRB,” Journal of Cosmology and Astroparticle Physics, vol. 2011, no. 11, Article ID 026, 2011.

L. Nava, G. Ghirlanda, G. Ghisellini, and A. Celotti, “Spectral properties of 438 GRBs detected by Fermi /GBM,” Astronomy and Astrophysics, vol. 530, Article ID A21, 2011.

G. A. MacLachlan, A. Shenoy, E. Sonbas, et al., “Minimum variability time scales of long and short GRBs,” In press, http://www.arxiv.org/abs/1201.4431/

K. Mannheim, “High-energy neutrinos from extragalactic jets,” Astroparticle Physics, vol. 3, no. 3, pp. 295–302, 1995.

J. Greiner, C. Clemens, T. Krühler et al., “The redshift and afterglow of the extremely energetic gamma-ray burst GRB 080916C,” Astronomy and Astrophysics, vol. 498, no. 1, pp. 89–94, 2009.

A. A. Abdo, M. Ackermann, M. Ajello et al., “Fermi observations of GRB 090902b: a distinct spectral component in the prompt and delayed emission,” Astrophysical Journal Letters, vol. 706, no. 1 PART 2, pp. L138–L144, 2009.

D. Gruber, T. Krühler, S. Foley et al., “Fermi/GBM observations of the ultra-long GRB 091024: a burst with an optical flash,” Astronomy and Astrophysics, vol. 528, Article ID A15, 2011.

S. Hümmer, Ph.D. thesis, Würzburg University, 2012.