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.

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.

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.

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 investigate non-thermal electron and nuclei energy losses within the binary neutron star merger remnant produced by the event GW170817. The lack of a cooling feature within the detected synchrotron emission from the source is used to constrain the magnetic field at the mG level, assuming that this emission is electron synchrotron in origin, and that the accelerated spectrum in the electrons follows the form $dN/dE_e \propto E_e^{-2}$. The level of subsequent gamma-ray emission from the source is demonstrated to provide a further constraint on the source magnetic field strength. We also put forward alternative strong ($\sim$G) magnetic field scenarios able to support this emission. For such stronger fields, the photo-disintegration of non-thermal nuclei within the source is considered, and a bottleneck period of $\sim$5-30 days is found when this process peaks. We find that this class of source is in principle able to support the population of cosmic rays detected at Earth below the "ankle".

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