Gamma-ray bursts are short-lived, luminous explosions at cosmological distances, thought to originate from relativistic jets launched at the deaths of massive stars. They are among the prime candidates to produce the observed cosmic rays at the highest energies. Recent neutrino data have, however, started to constrain this possibility in the simplest models with only one emission zone. In the classical theory of gamma-ray bursts, it is expected that particles are accelerated at mildly relativistic shocks generated by the collisions of material ejected from a central engine. We consider neutrino and cosmic-ray emission from multiple emission regions since these internal collisions must occur at very different radii, from below the photosphere all the way out to the circumburst medium, as a consequence of the efficient dissipation of kinetic energy. We demonstrate that the different messengers originate from different collision radii, which means that multi-messenger observations open windows for revealing the evolving GRB outflows.

In light of the latest IceCube data, we discuss the implications of the cosmic ray energy input from hypernovae and supernovae into the Universe, and their propagation in the hosting galaxy and galaxy clusters or groups. The magnetic confinement in these environments may lead to efficient $pp$ collisions, resulting in a diffuse neutrino spectrum extending from PeV down to 10 TeV energies, with a spectrum and flux level compatible with that recently reported by IceCube. If the diffuse 10 TeV neutrino background largely comes from such the CR reservoirs, the corresponding diffuse gamma-ray background should be compatible with the recent \textit{Fermi} data. In this scenario, the CR energy input from hypernovae should be dominant over that of supernovae, implying that the starburst scenario does not work if the supernova energy budget is a factor of two larger than the hypernova energy budget. Thus, this strong case scenario can be supported or ruled out in near future.

We reconsider the possibility that gamma-ray bursts (GRBs) are the sources of the ultra-high energy cosmic rays (UHECRs) within the internal shock model, assuming a pure proton composition of the UHECRs. For the first time, we combine the information from gamma-rays, cosmic rays, prompt neutrinos, and cosmogenic neutrinos quantitatively in a joint cosmic ray production and propagation model, and we show that the information on the cosmic energy budget can be obtained as a consequence. In addition to the neutron model, we consider alternative scenarios for the cosmic ray escape from the GRBs, i.e., that cosmic rays can leak from the sources. We find that the dip model, which describes the ankle in UHECR observations by the pair production dip, is strongly disfavored in combination with the internal shock model because a) unrealistically high baryonic loadings (energy in protons versus energy in electrons/gamma-rays) are needed for the individual GRBs and b) the prompt neutrino flux easily overshoots the corresponding neutrino bound. On the other hand, GRBs may account for the UHECRs in the ankle transition model if cosmic rays leak out from the source at the highest energies. In that case, we demonstrate that future neutrino observations can efficiently test most of the parameter space -- unless the baryonic loading is much larger than previously anticipated.

Gamma-ray bursts (GRBs) have long been held as one of the most promising sources of ultra-high energy (UHE) neutrinos. The internal shock model of GRB emission posits the joint production of UHE cosmic ray (UHECRs, above 10^8 GeV), photons, and neutrinos, through photohadronic interactions between source photons and magnetically-confined energetic protons, that occur when relativistically-expanding matter shells loaded with baryons collide with one another. While neutrino observations by IceCube have now ruled out the simplest version of the internal shock model, we show that a revised calculation of the emission, together with the consideration of the full photohadronic cross section and other particle physics effects, results in a prediction of the prompt GRB neutrino flux that still lies one order of magnitude below the current upper bounds, as recently exemplified by the results from ANTARES. In addition, we show that by allowing protons to directly escape their magnetic confinement without interacting at the source, we are able to partially decouple the cosmic ray and prompt neutrino emission, which grants the freedom to fit the UHECR observations while respecting the neutrino upper bounds. Finally, we briefly present advances towards pinning down the precise relation between UHECRs and UHE neutrinos, including the baryonic loading required to fit UHECR observations, and we will assess the role that very large volume neutrino telescopes play in this.

A search for muon neutrinos in coincidence with gamma-ray bursts with the ANTARES neutrino detector using data from the end of 2007 to 2011 is performed. Expected neutrino fluxes are calculated for each burst individually. The most recent numerical calculations of the spectra using the NeuCosmA code are employed, which include Monte Carlo simulations of the full underlying photohadronic interaction processes. The discovery probability for a selection of 296 gamma-ray bursts in the given period is optimised using an extended maximum-likelihood strategy. No significant excess over background is found in the data, and 90% confidence level upper limits are placed on the total expected flux according to the model.

In this work we give an estimate of the neutrino flux that can be expected from the microquasar Cyg X-3. We calculate the muon neutrino flux expected here on Earth as well as the corresponding number of neutrino events in the IceCube telescope based on the so-called hypersoft X-ray state of Cyg~X-3. If the average emission from Cyg~X-3 over a period of 5 yr were as high as during the used X-ray state, a total of 0.8 events should be observed by the full IceCube telescope. We also show that this conclusion holds by a factor of a few when we consider the other measured X-ray states. Using the correlation of AGILE data on the flaring episodes in 2009 June and July to the hypersoft X-ray state we calculate that the upper limits on the neutrino flux given by IceCube are starting to constrain the hadronic models, which have been introduced to interpret the high-energy emission detected by AGILE.

We present a model of ultra-high-energy cosmic ray (UHECR) production in the shock-accelerated fireball of a gamma-ray burst. In addition to the standard UHECR origin from neutron escape and decay into protons, our model considers direct proton emission through leakage from the edges of the accelerated baryon-loaded shells that make up the fireball. Depending on the optical thickness of the shells to photohadronic interactions, the source falls in one of three scenarios: the usual, optically thin source dominated by neutron escape, an optically thick source to neutron escape, or a "direct escape" scenario, where the main contribution to UHECRs comes from the leaked protons. The associated neutrino production will be different for each scenario, and we see that the standard "one neutrino per cosmic ray" assumption is valid only in the optically thin case, while more than one neutrino per cosmic ray is expected in the optically thick scenario. In addition, the extra direct escape component enhances the high-energy part of the UHECR flux, thus improving the agreement between the predictions and the observed flux.

The paradigm that gamma-ray burst (GRB) fireballs are the sources of the ultra-high energy cosmic rays (UHECRs) is being probed by neutrino observations. Very stringent bounds can be obtained from the cosmic ray (proton)--neutrino connection, assuming that the UHECRs escape as neutrons. In this study, we identify three different regimes as a function of the fireball parameters: the standard "one neutrino per cosmic ray" case, the optically thick (to neutron escape) case, and the case where leakage of protons from the boundaries of the shells (direct escape) dominates. In the optically thick regime, photomeson production is very efficient, and more neutrinos will be emitted per cosmic ray than in the standard case, whereas in the direct escape-dominated regime, more cosmic rays than neutrinos will be emitted. We demonstrate that, for efficient proton acceleration, which is required to describe the observed UHECR spectrum, the standard case only applies to a very narrow region of the fireball parameter space. We illustrate with several observed examples that conclusions on the cosmic ray--neutrino connection will depend on the actual burst parameters. We show that the definition of the pion production efficiency currently used by the IceCube collaboration underestimates the neutrino production in the optically thick case. Finally, we point out that the direct escape component leads to a spectral break in the cosmic ray spectrum emitted from a single source. The resulting "two-component model" can be used to even more strongly pronounce the spectral features of the observed UHECR spectrum than the dip model.

We discuss decays of ultra-relativistic neutrinos over cosmological distances by solving the decay equation in terms of its redshift dependence. We demonstrate that there are significant conceptual differences compared to more simplified treatments of neutrino decay. For instance, the maximum distance the neutrinos have traveled is limited by the Hubble length, which means that the common belief that longer neutrino lifetimes can be probed by longer distances does not apply. As a consequence, the neutrino lifetime limit from supernova 1987A cannot be exceeded by high-energy astrophysical neutrinos. We discuss the implications for neutrino spectra and flavor ratios from gamma-ray bursts as one example of extragalactic sources, using up-to-date neutrino flux predictions. If the observation of SN 1987A implies that \nu_1 is stable and the other mass eigenstates decay with rates much smaller than their current bounds, the muon track rate can be substantially suppressed compared to the cascade rate in the region IceCube is most sensitive to. In this scenario, no gamma-ray burst neutrinos may be found using muon tracks even with the full scale experiment, whereas reliable information on high-energy astrophysical sources can only be obtained from cascade measurements. As another consequence, the recently observed two cascade event candidates at PeV energies will not be accompanied by corresponding muon tracks.

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 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.