We scrutinize the paradigm that conventional long-duration Gamma-Ray Bursts (GRBs) are the dominant source of the ultra-high energy cosmic rays (UHECRs) within the internal shock scenario by describing UHECR spectrum and composition and by studying the predicted (source and cosmogenic) neutrino fluxes. Since it has been demonstrated that the stacking searches for astrophysical GRB neutrinos strongly constrain the parameter space in single-zone models, we focus on the dynamics of multiple collisions for which different messengers are expected to come from different regions of the same object. We propose a model which can describe both stochastic and deterministic engines, which we study in a systematic way. We find that GRBs can indeed describe the UHECRs for a wide range of different model assumptions with comparable quality; the required heavy mass fraction at injection is found to be larger than 70% (95% CL). We demonstrate that the post-dicted (from UHECR data) neutrino fluxes from sources and UHECR propagation are indeed below the current sensitivities but will be reached by the next generation of experiments. We finally critically review the required source energetics with the specific examples found in this study.

We introduce a phenomenological concordance scenario with a relativistic jet for the Tidal Disruption Event (TDE) AT2019dsg, which has been proposed as source of the astrophysical neutrino event IceCube-191001A. Noting that AT2019dsg is one of the brightest (and few) TDEs observed in X-rays, we connect the neutrino production with the X-rays: an expanding cocoon causes the progressive obscuration of the X-rays emitted by the accretion disk, while at the same time it provides a sufficiently intense external target of back-scattered X-rays for photo-pion production off protons. We also describe the late-term emission of the neutrino (about 150 days after the peak), by scaling the production radius with the black body radius. Our energetics and assumptions for the jet and the cocoon are compatible with expectations from numerical simulations of TDEs. We predict 0.26 neutrino events in the right energy range in IceCube.

We discuss the production of multiple astrophysical messengers (neutrinos, cosmic rays, gamma-rays) in the Gamma-Ray Burst (GRB) internal shock scenario, focusing on the impact of the collision dynamics between two shells on the fireball evolution. In addition to the inelastic case, in which plasma shells merge when they collide, we study the Ultra Efficient Shock scenario, in which a fraction of the internal energy is re-converted into kinetic energy and, consequently, the two shells survive and remain in the system. We find that in all cases a quasi-diffuse neutrino flux from GRBs at the level of $10^{-11}$ to $10^{-10} \, \mathrm{GeV \, cm^{-2} \, s^{-1} \, sr^{-1}}$ (per flavor) is expected for protons and a baryonic loading of ten, which is potentially within the reach of IceCube-Gen2. The highest impact of the collision model for multi-messenger production is observed for the Ultra Efficient Shock scenario, that promises high conversion efficiencies from kinetic to radiated energy. However, the assumption that the plasma shells separate after a collision and survive as separate shells within the fireball is found to be justified too rarely in a multi-collision model that uses hydrodynamical simulations with the PLUTO code for individual shell collisions.

Since astrophysical neutrinos are produced in the interactions of cosmic rays, identifying the origin of cosmic rays using directional correlations with neutrinos is one of the most interesting possibilities of the field. For that purpose, especially the Ultra-High Energy Cosmic Rays (UHECRs) are promising, as they are deflected less by extragalactic and Galactic magnetic fields than cosmic rays at lower energies. However, photo-hadronic interactions of the UHECRs limit their horizon, while neutrinos do not interact over cosmological distances. We study the possibility to search for anisotropies by investigating neutrino-UHECR correlations from the theoretical perspective, taking into account the UHECR horizon, magnetic-field deflections, and the cosmological source evolution. Under the assumption that the neutrinos and UHECRs all come from the same source class, we demonstrate that the non-observation of neutrino multiplets strongly constrains the possibility to find neutrino-UHECR correlations.

We demonstrate that a population of blazars can describe the observed spectrum and composition of Ultra-High-Energy Cosmic Rays (UHECRs), and that the dominant contribution comes from low-luminosity BL Lacs. However, a sub-dominant contribution from high-luminosity blazars can lead to a substantial neutrino flux which will peak at EeV energies. We find that it is plausible that this neutrino flux from the sources outshines the cosmogenic neutrino flux produced during the propagation of UHECRs. This has profound implications for the high-energy (EeV) neutrino detection experiments: their main target may be source neutrinos, for which additional search strategies can be used, such as stacking searches, flare analyses, or multi-messenger follow-ups.

Photon-hadronic interactions are important for the sources and the transport of Ultra-High Energy Cosmic Rays (UHECRs). Current state-of-the-art cosmic ray transport simulations handle nuclear disintegration at energies of the Giant Dipole Resonance at a more sophisticated level, as well as the photohadronic interactions of nucleons in the high-energy regime above the pion production threshold. However, the interactions of nuclei above the pion production threshold are commonly modeled by treating the nucleus as a superposition of free nucleons -- ignoring the effect of the nuclear medium. We construct an improved, inclusive model for the photomeson regime for nuclei with $A \leq 56$ by employing more accurate, data-driven parametrizations of the interaction cross section, the fragmentation of the primary nucleus and the inclusive pion production cross section that directly affects the production of astrophysical neutrinos. We apply our results to two multi-messenger scenarios (Tidal Disruption Events and Gamma-Ray Bursts) in which photonuclear interactions in the photomeson regime are the dominant cooling process for the highest energy cosmic rays. While we find moderate changes to the mass composition of UHECRs, the astrophysical neutrino fluxes exhibit a significant (factor of a few) reduction compared to the na\"ive superposition of free nucleons for sources of UHECR nuclei with a populated cascade. The numerical code implementing the model has been made publicly available, which facilitates the integration of our results in similar frameworks.

We discuss possible interpretations of the neutrinos observed from the AGN blazar TXS 0506+056 in the multi-messenger and multi-wavelength context, including both the 2014-15 and 2017 neutrino flares. While the neutrino observed in September 2017 has to describe contemporary data in e.g. the X-ray and VHE gamma-ray ranges, data at the 2014-15 excess are much sparser. We demonstrate that in both cases the simplest possible one-zone AGN blazar models face challenges. While the 2017 flare can be well interpreted by considering more sophisticated source geometries, the 2014-15 flare is much harder to describe with conventional models. One challenge is the energy injected into the electromagnetic cascade coming together with the neutrino production, which cannot be reconciled with the 13 observed neutrino events. We also speculate if a common interpretation of both flares is feasible.

Tidal Disruption Events (TDEs) are processes where stars are torn apart by the strong gravitational force near to a massive or supermassive black hole. If a jet is launched in such a process, particle acceleration may take place in internal shocks. We demonstrate that jetted TDEs can simultaneously describe the observed neutrino and cosmic ray fluxes at the highest energies if stars with heavier compositions, such as carbon-oxygen white dwarfs, are tidally disrupted and these events are sufficiently abundant. We simulate the photo-hadronic interactions both in the TDE jet and in the propagation through the extragalactic space and we show that the simultaneous description of Ultra-High Energy Cosmic Ray (UHECR) and PeV neutrino data implies that a nuclear cascade in the jet develops by photo-hadronic interactions.

In September 2017, the IceCube Neutrino Observatory recorded a very-high-energy neutrino in directional coincidence with a blazar in an unusually bright gamma-ray state, TXS0506+056. Blazars are prominent photon sources in the universe because they harbor a relativistic jet whose radiation is strongly collimated and amplified. High-energy atomic nuclei known as cosmic rays can produce neutrinos; thus the recent detection may help identifying the sources of the diffuse neutrino flux and the energetic cosmic rays. Here we report on a self-consistent analysis of the physical relation between the observed neutrino and the blazar, in particular the time evolution and spectral behavior of neutrino and photon emission. We demonstrate that a moderate enhancement in the number of cosmic rays during the flare can yield a very strong increase of the neutrino flux which is limited by co-produced hard X-rays and TeV gamma rays. We also test typical radiation models for compatibility and identify several model classes as incompatible with the observations. We investigate to what degree the findings can be generalized to the entire population of blazars, to determine the relation between their output in photons, neutrinos, and cosmic rays, and suggest how to optimize the strategy of future observations.

We investigate whether the emission of neutrinos observed in 2014-15 from the direction of the blazar TXS 0506+056 can be accommodated with leptohadronic multi-wavelength models of the source commonly adopted for the 2017 flare. While multi-wavelength data during the neutrino flare are sparse, the large number of neutrino events ($13\pm5$) challenges the missing activity in gamma rays. We illustrate that two to five neutrino events during the flare can be explained with leptohadronic models of different categories: a one-zone model, a compact core model, and an external radiation field model. If, however, significantly more events were to be accommodated, the predicted multi-wavelength emission levels would be in conflict with observational X-ray constraints, or with the high-energy gamma ray fluxes observed by the Fermi LAT, depending on the model. For example, while the external radiation field model can predict up to five neutrino events without violating X-ray constraints, the absorption of high-energy gamma rays is in minor tension with data. We therefore do not find any model that can simultaneously explain the high event number quoted by IceCube and the (sparse) electromagnetic data during the neutrino flare.

We study the implications of Ultra-High Energy Cosmic Ray (UHECR) data from the Pierre Auger Observatory for potential accelerator candidates and cosmogenic neutrino fluxes for different combinations of nuclear disintegration and air-shower models. We exploit the most recent spectral and mass composition data (2017) with a new, computationally very efficient simulation code PriNCe. We extend the systematic framework originally developed by the Pierre Auger Collaboration with the cosmological source evolution as an additional free parameter. In this framework, an ensemble of generalized UHECR accelerators is characterized by a universal spectral index (equal for all injection species), a maximal rigidity, and the normalizations for five nuclear element groups. We find that the 2017 data favor a small but constrained contribution of heavy elements (iron) at the source. We demonstrate that the results moderately depend on the nuclear disintegration (PSB, Peanut, or Talys) model, and more strongly on the air-shower (EPOS-LHC, Sibyll-2.3, or QGSjet-II-04) model. Variations of these models result in different source evolutions and spectral indices, limiting the interpretation in terms of a particular class of cosmic accelerators. Better constrained parameters include the maximal rigidity and the mass composition at the source. Hence, the cosmogenic neutrino flux can be robustly predicted, since it originates from interactions with the cosmic infrared background and peaks at $10^8 \, \mathrm{GeV}$. Depending on the source evolution at high redshifts the flux is likely out of reach of future neutrino observatories in most cases, and a minimal cosmogenic neutrino flux cannot be claimed from data without assuming a cosmological distribution of the sources.

We demonstrate that the UHECRs produced in the nuclear cascade in the jet of Low-Luminosity Gamma-Ray Bursts (LL-GRBs) can describe the UHECR spectrum and composition and, at the same time, the diffuse neutrino flux at the highest energies. The radiation density in the source simultaneously controls the neutrino production and the development of the nuclear cascade, leading to a flux of nucleons and light nuclei describing even the cosmic-ray ankle at $5 \times 10^{18}$ eV. The derived source parameters are consistent with population studies, indicating a baryonic loading factor of about ten. Our results motivate the continued experimental search of LL-GRBs as a unique GRB population.

We study if the diffuse astrophysical neutrino flux can come from blazar jets -- a subclass of Active Galactic Nuclei (AGNs) -- while it, at the same time, respects the blazar stacking limit based on source catalogs and is consistent with the observation from TXS 0506+056. We compute the neutrino flux from resolved and unresolved sources using an averaged, empirical relationship between electromagnetic spectrum and luminosity, known as the {\em blazar sequence}, for two populations of blazars (BL Lacs and FSRQs). Using a source model with realistic neutrino flux computations, we demonstrate that blazars can indeed power the diffuse neutrino flux at the highest energies and obey the stacking limit at the same time, and we derive the conditions for the baryonic loading (proton versus $\gamma$-ray luminosity) evolving over the blazar sequence. Under the hypothesis that low-luminosity blazars power the diffuse astrophysical neutrino flux, we find that the dominant contribution of the diffuse flux up to PeV energies must come from unresolved BL Lacs with baryonic loadings larger than about $10^5$ -- while only a very small contribution may come from resolved high-luminosity BL Lacs or FSRQs, which can be directly tested by the stacking limit. We find that the blazar TXS 0506+056 is on the verge of these populations in our baseline scenario, at a relatively high luminosity and redshift; as a consequence we predict about 0.3 $\gamma$-ray-neutrino associations per year from the whole population, dominated by BL Lacs with $L_\gamma \simeq 10^{45} \, \mathrm{erg/s}$ and $z \sim 0.1$.

High energy neutrinos can be produced by interactions of ultra-high energy cosmic rays (UHECRs) in the dense radiation fields of their sources as well as off the cosmic backgrounds when they propagate through the universe. Multi-messenger interpretations of current measurements deeply rely on the understanding of these interactions. In order to efficiently produce neutrinos in the sources of UHECRs, at least a moderate level of interactions is needed, which means that a nuclear cascade develops if nuclei are involved. On the other hand, the available cross-section data and interaction models turn out to make poor predictions for most nuclei heavier than protons. We show the impact of these uncertainties in state-of-the-art photo-disintegration models and motivate nuclear cross-section measurements. Further, we discuss extensions for photo-meson models currently used in astrophysics and demonstrate the importance of understanding the details of UHECR interaction with the Glashow resonance.

We discuss the production of ultra-high-energy cosmic ray (UHECR) nuclei and neutrinos from blazars. We compute the nuclear cascade in the jet for both BL Lac objects and flat-spectrum radio quasars (FSRQs), and in the ambient radiation zones for FSRQs as well. By modeling representative spectral energy distributions along the blazar sequence, two distinct regimes are identified, which we call "nuclear survival" -- typically found in low-luminosity BL Lacs, and "nuclear cascade" -- typically found in high-luminosity FSRQs. We quantify how the neutrino and cosmic-ray (CR) emission efficiencies evolve over the blazar sequence, and demonstrate that neutrinos and CRs come from very different object classes. For example, high-frequency peaked BL Lacs (HBLs) tend to produce CRs, and HL-FSRQs are the more efficient neutrino emitters. This conclusion does not depend on the CR escape mechanism, for which we discuss two alternatives (diffusive and advective escape). Finally, the neutrino spectrum from blazars is shown to significantly depend on the injection composition into the jet, especially in the nuclear cascade case: Injection compositions heavier than protons lead to reduced neutrino production at the peak, which moves at the same time to lower energies. Thus, these sources will exhibit better compatibility with the observed IceCube and UHECR data.

We propose a multi-component model for the observed diffuse neutrino flux, including the residual atmospheric backgrounds, a Galactic contribution (such as from cosmic ray interactions with gas), an extra-galactic contribution from pp interactions (such as from starburst galaxies) and a hard extragalatic contribution from photo-hadronic interactions at the highest energies (such as from Tidal Disruption Events or Active Galactic Nuclei). We demonstrate that this model can address the key problems of astrophysical neutrino data, such as the different observed spectral indices in the high-energy starting and through-going muon samples, a possible anisotropy due to Galactic events, the non-observation of point sources, and the constraint from the extragalatic diffuse gamma-ray background. Furthermore, the recently observed muon track with a deposited energy of 4.5 PeV might be interpreted as evidence for the extragalactic photo-hadronic contribution. We perform the analysis based on the observed events instead of the unfolded fluxes by computing the probability distributions for the event type and reconstructed neutrino energy. As a consequence, we give the probability to belong to each of these astrophysical components on an event-to-event basis.

We compute the expected neutrino fluence from SGRB 170817A, associated with the gravitational wave event GW 170817, directly based on Fermi observations in two scenarios: structured jet and off-axis (observed) top-hat jet. While the expected neutrino fluence for the structured jet case is very small, large off-axis angles imply high radiation densities in the jet, which can enhance the neutrino production efficiency. In the most optimistic allowed scenario, the neutrino fluence can reach only $10^{-4}$ of the sensitivity of the neutrino telescopes. We furthermore demonstrate that the fact that gamma-rays can escape limits the baryonic loading (energy in protons versus photons) and the off-axis angle for the internal shock scenario. In particular, for a baryonic loading of ten, the off-axis angle is more strongly constrained by the baryonic loading than by the time delay between the gravitational wave event and the onset of the gamma-ray emission.

We systematically study the allowed parameter space for the flavor composition of astrophysical neutrinos measured at Earth, including beyond the Standard Model theories at production, during propagation, and at detection. One motivation is to illustrate the discrimination power of the next-generation neutrino telescopes such as IceCube-Gen2. We identify several examples that lead to potential deviations from the standard neutrino mixing expectation such as significant sterile neutrino production at the source, effective operators modifying the neutrino propagation at high energies, dark matter interactions in neutrino propagation, or non-standard interactions in Earth matter. IceCube-Gen2 can exclude about 90% of the allowed parameter space in these cases, and hence will allow to efficiently test and discriminate models. More detailed information can be obtained from additional observables such as the energy-dependence of the effect, fraction of electron antineutrinos at the Glashow resonance, or number of tau neutrino events.

The determination of the injection composition of cosmic ray nuclei within astrophysical sources requires sufficiently accurate descriptions of the source physics and the propagation - apart from controlling astrophysical uncertainties. We therefore study the implications of nuclear data and models for cosmic ray astrophysics, which involves the photo-disintegration of nuclei up to iron in astrophysical environments. We demonstrate that the impact of nuclear model uncertainties is potentially larger in environments with non-thermal radiation fields than in the cosmic microwave background. We also study the impact of nuclear models on the nuclear cascade in a gamma-ray burst radiation field, simulated at a level of complexity comparable to the most precise cosmic ray propagation code. We conclude with an isotope chart describing which information is in principle necessary to describe nuclear interactions in cosmic ray sources and propagation.

We study the frequently used assumption in multi-messenger astrophysics that the gamma-ray and neutrino fluxes are directly connected because they are assumed to be produced by the same photohadronic production chain. An interesting candidate source for this test is the flat-spectrum radio quasar PKS B1424-418, which recently called attention of a potential correlation between an IceCube PeV-neutrino event and its burst phase. We simulate both the multi-waveband photon and the neutrino emission from this source using a self-consistent radiation model. We demonstrate that a simple hadronic model cannot adequately describe the spectral energy distribution for this source, but a lepto-hadronic model with sub-dominant hadronic component can reproduce the multi-waveband photon spectrum observed during various activity phases of the blazar. As a conclusion, up to about 0.3 neutrino events may coincide with the burst, which implies that the leptonic contribution dominates in the relevant energy band. We also demonstrate that the time-wise correlation between the neutrino event and burst phase is weak.

We discuss neutrino and cosmic-ray emission from Gamma-Ray Bursts (GRBs) with the injection of nuclei, where we take into account that a nuclear cascade from photo-disintegration can fully develop in the source. One of our main objectives is to test if recent results from the IceCube and the Pierre Auger Observatory can be accommodated with the paradigm that GRBs are the sources of Ultra-High Energy Cosmic Rays (UHECRs). While our key results are obtained using an internal shock model, we discuss how the secondary emission from a GRB shell can be interpreted in terms of other astrophysical models. It is demonstrated that the expected neutrino flux from GRBs weakly depends on the injection composition, which implies that prompt neutrinos from GRBs can efficiently test the GRB-UHECR paradigm even if the UHECRs are nuclei. We show that the UHECR spectrum and composition, as measured by the Pierre Auger Observatory, can be self-consistently reproduced in a combined source-propagation model. In an attempt to describe the energy range including the ankle, we find tension with the IceCube bounds from the GRB stacking analyses. In an alternative scenario, where only the UHECRs beyond the ankle originate from GRBs, the requirement for a joint description of cosmic-ray and neutrino observations favors lower luminosities, which does not correspond to the typical expectation from {\gamma}-ray observations.

We study the production of high energy neutrinos in jets from the tidal disruption of stars by supermassive black holes. The diffuse neutrino flux expected from these tidal disruption events (TDEs) is calculated both analytically and numerically, taking account the dependence of the rate of TDEs on the redshift and black hole mass. We find that ~10% of the observed diffuse flux at IceCube at an energy of about 1 PeV can come from TDEs if the characteristics of known jetted tidal disruption events are assumed to apply to the whole population of these sources. If, however, plausible scalings of the jet Lorentz factor or variability timescale with the black hole mass are taken into account, the contribution of the lowest mass black holes to the neutrino flux is enhanced. In this case, TDEs can account for most of the neutrino flux detected at IceCube, describing both the neutrino flux normalization and spectral shape with moderate baryonic loadings. While the uncertainties on our assumptions are large, a possible signature of TDEs as the origin of the IceCube signal is the transition of the flux flavor composition from a pion beam to a muon damped source at the highest energies, which will also result in a suppression of Glashow resonance events.

We study the Glashow resonance $\bar{\nu}_e + e^- \rightarrow W^- \rightarrow$ hadrons at 6.3 PeV as diagnostic of the production processes of ultra-high energy neutrinos. The focus lies on describing the physics of neutrino production from pion decay as accurate as possible by including the kinematics of weak decays and Monte Carlo simulations of pp and p$\gamma$ interactions. We discuss optically thick (to photohadronic interactions) sources, sources of cosmic ray nuclei and muon damped sources. Even in the proposed upgrade IceCube-Gen2, a discrimination of scenarios such as pp versus p$\gamma$ is extremely challenging under realistic assumptions. Nonetheless, the Glashow resonance can serve as a smoking gun signature of neutrino production from photohadronic (A$\gamma$) interactions of heavier nuclei, as the expected Glashow event rate exceeds that of pp interactions. We finally quantify the exposures for which the non-observation of Glashow events exerts pressure on certain scenarios.

Gamma-ray bursts (GRBs) are promising as sources of neutrinos and cosmic rays. In the internal shock scenario, blobs of plasma emitted from a central engine collide within a relativistic jet and form shocks, leading to particle acceleration and emission. Motivated by present experimental constraints and sensitivities, we improve the predictions of particle emission by investigating time-dependent effects from multiple shocks. We produce synthetic light curves with different variability timescales that stem from properties of the central engine. For individual GRBs, qualitative conclusions about model parameters, neutrino production efficiency, and delays in high-energy gamma rays can be deduced from inspection of the gamma-ray light curves. GRBs with fast time variability without additional prominent pulse structure tend to be efficient neutrino emitters, whereas GRBs with fast variability modulated by a broad pulse structure can be inefficient neutrino emitters and produce delayed high-energy gamma-ray signals. Our results can be applied to quantitative tests of the GRB origin of ultra-high-energy cosmic rays, and have the potential to impact current and future multi-messenger searches.

The origin and composition of ultra-high-energy cosmic rays (UHECRs) remain a mystery. The proton dip model describes their spectral shape in the energy range above $10^9$ GeV by pair production and photohadronic interactions with the cosmic microwave background. The photohadronic interactions also produce cosmogenic neutrinos peaking around $10^9$ GeV. We test whether this model is still viable in light of recent UHECR spectrum measurements from the Telescope Array experiment, and upper limits on the cosmogenic neutrino flux from IceCube. While two-parameter fits have been already presented, we perform a full scan of the three main physical model parameters: source redshift evolution, injected proton maximal energy, and spectral index. We find qualitatively different conclusions compared to earlier two-parameter fits in the literature: a mild preference for a maximal energy cutoff at the sources instead of the Greisen--Zatsepin--Kuzmin (GZK) cutoff, hard injection spectra, and strong source evolution. The predicted cosmogenic neutrino flux exceeds the IceCube limit for any parameter combination. As a result, the proton dip model is challenged at more than 95\% C.L. This is strong evidence against this model independent of mass composition measurements.

The flavor composition of high-energy astrophysical neutrinos can reveal the physics governing their production, propagation, and interaction. The IceCube Collaboration has published the first experimental determination of the ratio of the flux in each flavor to the total. We present, as a theoretical counterpart, new results for the allowed ranges of flavor ratios at Earth for arbitrary flavor ratios in the sources. Our results will allow IceCube to more quickly identify when their data imply standard physics, a general class of new physics with arbitrary (incoherent) combinations of mass eigenstates, or new physics that goes beyond that, e.g., with terms that dominate the Hamiltonian at high energy.

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.

IceCube have observed neutrinos which are presumably of extra-galactic origin. Since specific sources have not yet been identified, we discuss what could be learned from the conceptual point of view. We use a simple model for neutrino production from the interactions between nuclei and matter, and we focus on the description of the spectral shape and flavor composition observed by IceCube. Our main parameters are spectral index, maximal energy, magnetic field, and composition of the accelerated nuclei. We show that a cutoff at PeV energies can be achieved by soft enough spectra, a cutoff of the primary energy, or strong enough magnetic fields. These options, however, are difficult to reconcile with the hypothesis that these neutrinos originate from the same sources as the ultra-high energy cosmic rays. We demonstrate that heavier nuclei accelerated in the sources may be a possible way out if the maximal energy scales appropriately with the mass number of the nuclei. In this scenario, neutrino observations can actually be used to test the UHECR acceleration mechanism. We also emphasize the need for a volume upgrade of the IceCube detector for future precision physics, for which the flavor information becomes a statistical meaningful model discriminator as qualitatively new ingredient.

Cosmic rays diffuse through the interstellar medium and interact with matter and radiations as long as they are trapped in the Galactic magnetic field. The IceCube experiment has detected some TeV-PeV neutrino events whose origin is yet unknown. We study if all or a fraction of these events can be described by the interactions of cosmic rays with matter. We consider the average target density needed to explain them for different halo sizes and shapes, the effect of the chemical composition of the cosmic rays, the impact of the directional information of the neutrino events, and the constraints from gamma ray bounds and their direction. We do not require knowledge of the cosmic ray escape time or injection for our approach. We find that, given all constraints, at most 0.1 of the observed neutrino events in Ice-Cube can be described by cosmic ray interactions with matter. In addition, we demonstrate that the currently established chemical composition of the cosmic rays contradicts a peak of the neutrino spectrum at PeV energies.

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.

We discuss the acceleration of secondary muons, pions, and kaons in gamma-ray bursts within the internal shock scenario, and their impact on the neutrino fluxes. We introduce a two-zone model consisting of an acceleration zone (the shocks) and a radiation zone (the plasma downstream the shocks). The acceleration in the shocks, which is an unavoidable consequence of the efficient proton acceleration, requires efficient transport from the radiation back to the acceleration zone. On the other hand, stochastic acceleration in the radiation zone can enhance the secondary spectra of muons and kaons significantly if there is a sufficiently large turbulent region. Overall, it is plausible that neutrino spectra can be enhanced by up to a factor of two at the peak by stochastic acceleration, that an additional spectral peaks appears from shock acceleration of the secondary muons and pions, and that the neutrino production from kaon decays is enhanced. Depending on the GRB parameters, the general conclusions concerning the limits to the internal shock scenario obtained by recent IceCube and ANTARES analyses may be affected by up to a factor of two by secondary acceleration. Most of the changes occur at energies above 10^7 GeV, so the effects for next-generation radio-detection experiments will be more pronounced. In the future, however, if GRBs are detected as high-energy neutrino sources, the detection of one or several pronounced peaks around 10^6 GeV or higher energies could help to derive the basic properties of the magnetic field strength in the GRB.

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.

We perform an unbiased search of the origin of the recently observed 28 events above ~30 TeV in the IceCube neutrino observatory, assuming that these are (apart from the atmospheric background) of astrophysical origin produced by photohadronic interactions. Instead of relying on the normalization of the neutrino flux, we demonstrate that spectral shape and flavor composition can be used to constrain or identify the source class. In order to quantify our observations, we use a model where the target photons are produced by the synchrotron emission of co-accelerated electrons, and we include magnetic field effects on the secondary muons, pions, and kaons. We find that the lack of observed events with energies much larger than PeV points towards sources with strong magnetic fields, which do not exhibit a direct correlation between highest cosmic ray and neutrino energies. While the simplest AGN models with efficient proton acceleration plausibly describe the current data at about the 3sigma confidence level, we show that IceCube can rule out that the observed neutrinos stem from the sources of the ultra-high energy cosmic rays with a factor of ten increased statistics at more than 5sigma if the current observations are confirmed. A possible caveat are sources with strong magnetic fields and high Lorentz factors, such as magnetic energy dominated GRBs.

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.

We discuss the neutrino mass hierarchy determination with atmospheric neutrinos in PINGU (Precision IceCube Next Generation Upgrade), based on a simulation with the GLoBES software including the full three flavor framework and parameter degeneracy, and we compare it to long-baseline experiment options. We demonstrate that the atmospheric mass hierarchy sensitivity depends on the achievable experiment properties and we identify the main targets for optimization, whereas the impact of a large number of tested systematical errors turns out to be small. Depending on the values of theta_23, delta, and the true hierarchy, a 90% CL to 3sigma discovery after three years of operation seems conceivable. We also emphasize the synergy with existing beam and reactor experiments, driven by NOvA, such as the complementary coverage of the parameter space. Finally, we point out that a low intensity neutrino beam with a relatively short decay pipe could be used to determine the mass hierarchy with a sensitivity comparable to the LBNE experiment irrespective of the directional resolution of the detector.

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.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

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 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, kaons), and flavor mixing, where we focus on p-gamma 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 p-gamma from p-p 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 co-accelerated electrons and 2) the target photon spectrum follows the observed photon spectrum of gamma-ray bursts. In the latter case, the multi-messenger extrapolation from the gamma-ray fluence to the expected neutrino flux is highlighted.

There are recent considerations to increase the photomultiplier density in the IceCube detector array beyond that of DeepCore, which will lead to a lower detection threshold and a huge fiducial mass for the neutrino detection. This initiative is known as "Phased IceCube Next Generation Upgrade" (PINGU). We discuss the possibility to send a neutrino beam from one of the major accelerator laboratories in the Northern hemisphere to such a detector. Such an experiment would be unique in the sense that it would be the only neutrino beam where the baseline crosses the Earth's core. We study the detector requirements for a beta beam, a neutrino factory beam, and a superbeam, where we consider both the cases of small theta_13 and large theta_13, as suggested by the recent T2K and Double Chooz results. We illustrate that a flavor-clean beta beam best suits the requirements of such a detector, in particular, that PINGU may replace a magic baseline detector for small values of theta_13 -- even in the absence of any energy resolution capability. For large theta_13, however, a single-baseline beta beam experiment cannot compete if it is constrained by the CERN-SPS. For a neutrino factory, because of the missing charge identification possibility in the detector, a very good energy resolution is required. If this can be achieved, especially a low energy neutrino factory, which does not suffer from the tau contamination, may be an interesting option for large theta_13. For the superbeam, where we use the LBNE beam as a reference, electron neutrino flavor identification and statistics are two of the main limitations. Finally, we demonstrate that, at least in principle, neutrino factory and superbeam can measure the density of the Earth's core to the sub-percent level for sin^2 2theta_13 larger than 0.01.

We discuss the interplay between spectral shape and detector response beyond a simple E^-2 neutrino flux at neutrino telescopes, at the example of time-integrated point source searches using IceCube-40 data. We use a self-consistent model for the neutrino production, in which protons interact with synchrotron photons from co-accelerated electrons, and we fully take into account the relevant pion and kaon production modes, the flavor composition at the source, flavor mixing, and magnetic field effects on the secondaries (pions, muon, and kaons). Since some of the model parameters can be related to the Hillas parameters R (size of the acceleration region) and B (magnetic field), we relate the detector response to the Hillas plane. In order to compare the response to different spectral shapes, we use the energy flux density as a measure for the pion production efficiency times luminosity of the source. We demonstrate that IceCube has a very good reach in this quantity for AGN nuclei and jets for all source declinations, while the spectra of sources with strong magnetic fields are found outside the optimal reach. We also demonstrate where neutrinos from kaon decays and muon tracks from tau decays can be relevant for the detector response. Finally, we point out the complementarity between IceCube and other experiments sensitive to high-energy neutrinos, at the example of 2004-2008 Earth-skimming neutrino data from Auger. We illustrate that Auger, in principle, is better sensitive to the parameter region in the Hillas plane from which the highest-energetic cosmic rays may be expected in this model.

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 discuss the importance of flavor ratio measurements in neutrino telescopes, such as by measuring the ratio between muon tracks to cascades, for the purpose of extracting new physics signals encountered by astrophysical neutrinos during propagation from the source to the detector. The detected flavor ratios not only carry the energy information of specific new physics scenarios which alter the transition probabilities in distinctive ways, but also the energy dependent flavor composition at the source. In the present work, we discuss the interplay of these two energy dependent effects and identify which new physics scenarios can be distinguished from the detected flavor ratios as a function of astrophysical parameters. We use a recently developed self-consistent neutrino production model as our toy model to generate energy dependent source flavor ratios and discuss (invisible) neutrino decay and quantum decoherence as specific new physics examples. Furthermore, we identify potentially interesting classes of sources on the Hillas plot for the purpose of new physics searches. We find that sources with substantial magnetic fields 10^3 Gauss <= B <= 10^6 Gauss, such as Active Galactic Nuclei (AGN) cores, white dwarfs, or maybe gamma-ray bursts, have, in principle, the best discrimination power for the considered new physics scenarios, whereas AGN jets, which typically perform as pion beam sources, can only discriminate few sub cases in the new physics effects. The optimal parameter region somewhat depends on the class of new physics effect considered.

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.

We discuss the interplay of neutrino oscillation and decay properties at neutrino telescopes. Motivated by recent unparticle scenarios, which open the possibility of new neutrino decay modes over astrophysical distances, we perform a complete classification of possible decay schemes, and we illustrate how different scenarios can be identified. Moreover, we show that the sensitivity of neutrino telescopes to standard neutrino properties, such as the mass hierarchy or delta_CP, is greatly enhanced in specific decay scenarios. In particular we discuss the impact of an astrophysical neutrino detection on terrestrial experiments, such as on the mass hierarchy measurement at NOvA. For example, we find that the scenario where only nu_1 is stable can be uniquely identified against all the other decay schemes, and that in this case CP violation can be established (for large theta_13) by the combination of Double Chooz with the track-to-shower ratio at a neutrino telescope, even if the flavor composition at the source is unknown. Our statements are based on a complete analysis of all the present solar, atmospheric, reactor and accelerator neutrino data, as well as on realistic simulation of future terrestrial neutrino oscillation experiments.

We discuss the possible impact of astrophysical neutrino flux measurements at neutrino telescopes on the neutrino oscillation program of reactor experiments and neutrino beams. We consider neutrino fluxes from neutron sources, muon damped sources, and pion sources, where we parameterize the input from these sources in terms of the flux ratio $R=\phi_\mu/(\phi_e+\phi_\tau)$ which can be extracted from the muon track to shower ratio in a neutrino telescope. While it is difficult to obtain any information from this ratio alone, we demonstrate that the dependence on the oscillation parameters is very complementary to the one of reactor experiments and neutrino beams. We find that for large values of $\sin^2 2 \theta_{13}$, a measurement of R with a precision of about 20% or better may not only improve the measurement of the leptonic CP phase, but also help the determination of the mass hierarchy. In some cases, early information on $\delta_{CP}$ may even be obtained from Double Chooz and an astrophysical flux alone without the help of superbeams. For small values of $\sin^2 2 \theta_{13}$, we find that using the information from an astrophysical neutrino flux could eliminate the octant degeneracy better than reactor experiments and beams alone. Finally, we demonstrate that implementing an additional observable based on the electromagnetic to hadronic shower ratio at a neutrino telescope (such as at higher energies) could be especially beneficial for pion beam sources.

We investigate the possibility to use the neutrinos coming from a future galactic supernova explosion to perform neutrino oscillation tomography of the Earth's core. We propose to use existing or planned detectors, resulting in an additional payoff. Provided that all of the discussed uncertainties can be reduced as expected, we find that the average matter densities of the Earth's inner and outer cores could be measured with a precision competitive with geophysics. However, since seismic waves are more sensitive to matter density jumps than average matter densities, neutrino physics would give partly complementary information.

Supernova neutrinos could be well-suited for probing neutrino decay, since decay may be observed even for very small decay rates or coupling constants. We will introduce an effective operator framework for the combined description of neutrino decay and neutrino oscillations for supernova neutrinos, which can especially take into account two properties: One is the radially symmetric neutrino flux, allowing a decay product to be re-directed towards the observer even if the parent neutrino had a different original direction of propagation. The other is decoherence because of the long baselines for coherently produced neutrinos. We will demonstrate how to use this effective theory to calculate the time-dependent fluxes at the detector. In addition, we will show the implications of a Majoron-like decay model. As a result, we will demonstrate that for certain parameter values one may observe some effects which could also mimic signals similar to the ones expected from supernova models, making it in general harder to separate neutrino and supernova properties.

Neutrino decay in vacuum has often been considered as an alternative to neutrino oscillations. Because non-zero neutrino masses imply the possibility of both neutrino decay and neutrino oscillations, we present a model-independent formal treatment of these combined scenarios. For that, we show for the example of Majoron decay that in many cases decay products are observable and may even oscillate. Furthermore, we construct a minimal scenario where we study the physical implications of neutrino oscillations with intermediate decays.