We review the neutrino flux from gamma-ray bursts, which is estimated from gamma-ray observations and used for the interpretation of recent IceCube data, from a particle physics perspective. We numerically calculate the neutrino flux for the same astrophysical assumptions as the analytical fireball neutrino model, including the dominant pion and kaon production modes, flavor mixing, and magnetic field effects on the secondary muons, pions, and kaons. We demonstrate that taking into account the full energy dependencies of all spectra, the normalization of the expected neutrino flux reduces by about one order of magnitude and the spectrum shifts to higher energies, where we can pin down the exact origin of the discrepancies by the re-computation of the analytical models. We also reproduce the IceCube-40 analysis for exactly the same bursts and same assumptions and illustrate the impact of uncertainties. We conclude that the baryonic loading of the fireballs, which is an important control parameter for the emission of cosmic rays, can be constrained significantly with the full-scale experiment after about ten years.

Gamma-ray burst analyses at neutrino telescopes are typically based on diffuse or stacked (i.e., aggregated) neutrino fluxes, because the number of events expected from a single burst is small. The interpretation of aggregated flux limits implies new systematics not present for a single burst, such as by the integration over parameter distributions (diffuse fluxes), or by the low statistics in small burst samples (stacked fluxes). We simulate parameter distributions with a Monte Carlo method computing the spectra burst by burst, as compared to a conventional Monte Carlo integration. With this approach, we can predict the behavior of the flux in the diffuse limit as well as in low statistics stacking samples, such as used in recent IceCube data analyses. We also include the flavor composition at the detector (ratio between muon tracks and cascades) into our considerations. We demonstrate that the spectral features, such as a characteristic multi-peak structure coming from photohadronic interactions, flavor mixing, and magnetic field effects, are typically present even in diffuse neutrino fluxes if only the redshift distribution of the sources is considered, with z \simeq 1 dominating the neutrino flux. On the other hand, we show that variations of the Lorentz boost can only be interpreted in a model-dependent way, and can be used as a model discriminator. For example, we illustrate that the observation of spectral features in aggregated fluxes will disfavor the commonly used assumption that bursts with small Lorentz factors dominate the neutrino flux, whereas it will be consistent with the hypothesis that the bursts have similar properties in the comoving frame.

We discuss neutrino fluxes and energy dependent flavor ratios of cosmic accelerators as a function of the size of the acceleration region and the magnetic field, which are the parameters of the Hillas plot. We assume that photohadronic interactions between Fermi accelerated protons and synchrotron photons from co-accelerated electrons (or positrons) lead to charged pion production. We include synchrotron cooling of the charged pions and successively produced muons, which then further decay into neutrinos, as well as the helicity dependence of the muon decays. Our photohadronic interaction model includes direct production, higher resonances, and high energy processes in order to model the neutrino flavor and neutrino-antineutrino ratios with sufficient accuracy. Since we assume that the sources are optically thin to neutron interactions, we include the neutrino fluxes from neutron decays. We classify the Hillas plot into regions with different characteristic flavor ratios over 20x24 orders of magnitude in R and B. In some examples with sizable magnetic fields, we recover neutron beam, pion beam, and muon damped beam as a function of energy. However, we also find anomalous or new sources, such as muon beam sources with a flavor ratio nu_e:nu_mu:nu_tau of 1:1:0 in energy regions where the synchrotron-damped muons pile up. We also discuss the use of the Glashow resonance to identify pgamma optically thin sources with a strong imbalance between pi^+ and pi^- production. We find that these can, in principle, be identified in most cases in spite of the pi^- contamination from high energy photohadronic processes and the mixing parameter uncertainties.

We reanalyze the prompt muon neutrino flux from gamma-ray bursts (GRBs), at the example of the often used reference Waxman-Bahcall GRB flux, in terms of the particle physics involved. We first reproduce this reference flux treating synchrotron energy losses of the secondary pions explicitly. Then we include additional neutrino production modes, the neutrinos from muon decays, the magnetic field effects on all secondary species, and flavor mixing with the current parameter uncertainties. We demonstrate that the combination of these effects modifies the shape of the original Waxman-Bahcall GRB flux significantly, and changes the normalization by a factor of three to four. As a consequence, the gamma-ray burst search strategy of neutrino telescopes may be based on the wrong flux shape, and the constraints derived for the GRB neutrino flux, such as the baryonic loading, may in fact be already much stronger than anticipated.

We discuss simplified models for photo-meson production in cosmic accelerators, such as Active Galactic Nuclei and Gamma-Ray Bursts. Our self-consistent models are directly based on the underlying physics used in the SOPHIA software, and can be easily adapted if new data are included. They allow for the efficient computation of neutrino and photon spectra (from pi^0 decays), as a major requirement of modern time-dependent simulations of the astrophysical sources and parameter studies. In addition, the secondaries (pions and muons) are explicitely generated, a necessity if cooling processes are to be included. For the neutrino production, we include the helicity dependence of the muon decays which in fact leads to larger corrections than the details of the interaction model. The separate computation of the pi^0, pi^+, and pi^- fluxes allows, for instance, for flavor ratio predictions of the neutrinos at the source, which are a requirement of many tests of neutrino properties using astrophysical sources. We confirm that for charged pion generation, the often used production by the Delta(1232)-resonance is typically not the dominant process in Active Galactic Nuclei and Gamma-Ray Bursts, and we show, for arbitrary input spectra, that the number of neutrinos are underestimated by at least a factor of two if they are obtained from the neutral to charged pion ratio. We compare our results for several levels of simplification using isotropic synchrotron and thermal spectra, and we demonstrate that they are sufficiently close to the SOPHIA software.