Cascade equations and hadronic interactions at high energies

The new matrix form of cascade equations is convenient to discuss numerical properties and to investigate the role of hadrons and phase-space for atmospheric leptons. Precision is often limited by hadronic interactions models. I improved the event generator DPMJET-III using experimental data including the LHC. A specific curvature in charged-particle multiplicity distributions raised questions about elastic scattering amplitude and the black disk limit at LHC energies

Nuclear Physics Meets the Sources of the Ultra-High Energy Cosmic Rays

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.Comment: 11 pages, 6 figures. Publication is open acces

Cosmic-Ray and Neutrino Emission from Gamma-Ray Bursts with a Nuclear Cascade

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.Comment: 36 pages, 25 figure

A new view on Auger data and cosmogenic neutrinos in light of different nuclear disintegration and air-shower models

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.Comment: 21 pages, 11 figures. Accepted for publication in Ap

Phenomenology of atmospheric neutrinos

The detection of astrophysical neutrinos, certainly a break-through result, introduced new experimental challenges and fundamental questions about acceleration mechanisms of cosmic rays. On one hand IceCube succeeded in finding an unambiguous proof for the existence of a diffuse astrophysical neutrino flux, on the other hand the precise determination of its spectral index and normalization requires a better knowledge about the atmospheric background at hundreds of TeV and PeV energies. Atmospheric neutrinos in this energy range originate mostly from decays of heavy-flavor mesons, which production in the phase space relevant for prompt leptons is uncertain. Current accelerator-based experiments are limited by detector acceptance and not so much by the collision energy. This paper recaps phenomenological aspects of atmospheric leptons and calculation methods, linking recent progress in flux predictions with particle physics at colliders, in particular the Large Hadron Collider

daemonflux: DAta-drivEn MuOn-calibrated Neutrino Flux

In this paper, we present a refined calculation of the atmospheric neutrino flux spanning from GeV to PeV energies. Our method, Daemonflux, utilizes data-driven inputs and incorporates adjustable parameters to take their uncertainties into account. By optimizing these parameters using a combination of muon data and constraints from fixed-target experiments, we achieve uncertainties in the calculated neutrino fluxes of less than 10% up to 1 TeV, with neutrino ratios constrained to below 10%. Our model performs particularly well at energies below 100 GeV, where the smallest errors are obtained. We make our model available as a software package that provides access to predictions of fluxes, ratios, and errors, including the covariance matrix obtained from the fit.Comment: 21 pages, 23 figure

Neutrinos and Ultra-High-Energy Cosmic-Ray Nuclei from Blazars

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.Comment: 23 pages, 20 figure

IceCube Neutrinos from Hadronically Powered Gamma-Ray Galaxies

In this work we use a multi-messenger approach to determine if the high energy diffuse neutrino flux observed by the IceCube Observatory can originate from $\gamma$-ray sources powered by Cosmic Rays interactions with gas. Typical representatives of such sources are Starburst and Ultra-Luminous Infrared Galaxies. Using the three most recent calculations of the non-blazar contribution to the extragalactic $\gamma$-ray background measured by the Fermi-LAT collaboration, we find that a hard power-law spectrum with spectral index $\alpha \leq 2.12$ is compatible with all the estimations for the allowed contribution from non-blazar sources, within 1$\sigma$. Using such a spectrum we are able to interpret the IceCube results, showing that various classes of hadronically powered $\gamma$-ray galaxies can provide the dominant contribution to the astrophysical signal above 100 TeV and about half of the contribution to the energy flux between 10-100 TeV. With the addition of neutrinos from the Galactic plane, it is possible to saturate the IceCube signal at high energy. Our result shows that these sources are still well motivated candidates.Comment: Accepted for publication on JCA

Astrophysical neutrino production and impact of associated uncertainties in photo-hadronic interactions of UHECRs

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.Comment: invited talk presented at ISVHECRI 2018, submitted for publication to EPJ Web of Conference