The detection of the first electromagnetic counterpart to the binary neutron star (BNS) merger remnant GW170817 established the connection between short $\gamma$-ray bursts and BNS mergers. It also confirmed the forging of heavy elements in the ejecta (a so-called kilonova) via the r-process nucleosynthesis. The appearance of non-thermal radio and X-ray emission, as well as the brightening, which lasted more than 100 days, were somewhat unexpected. Current theoretical models attempt to explain this temporal behavior as either originating from a relativistic off-axis jet or a kilonova-like outflow. In either scenario, there is some ambiguity regarding how much energy is transported in the non-thermal electrons versus the magnetic field of the emission region. Combining the VLA (radio) and Chandra (X-ray) measurements with observations in the GeV-TeV domain can help break this ambiguity, almost independently of the assumed origin of the emission. Here we report for the first time on deep H.E.S.S. observations of GW170817 / GRB 170817A between 124 and 272 days after the BNS merger with the full H.E.S.S. array of telescopes, as well as on an updated analysis of the prompt (<5 days) observations with the upgraded H.E.S.S. phase-I telescopes. We discuss implications of the H.E.S.S. measurement for the magnetic field in the context of different source scenarios.

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.

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.

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

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

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.