We introduce the new version of SimProp, a Monte Carlo code for simulating the propagation of ultra-high energy cosmic rays in intergalactic space. This version, SimProp v2r4, together with an overall improvement of the code capabilities with a substantial reduction in the computation time, also computes secondary cosmogenic particles such as electron-positron pairs and gamma rays produced during the propagation of ultra-high energy cosmic rays. As recently pointed out by several authors, the flux of this secondary radiation and its products, within reach of the current observatories, provides useful information about models of ultra-high energy cosmic ray sources which would be hard to discriminate otherwise.

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.

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.

Interactions of ultra-high-energy cosmic rays with background photons set a limit to the distance cosmic rays reaching us above a certain energy can originate from, making measurements of their fluxes insensitive to properties of sources at high redshifts. On the other hand, the secondary PeV--EeV neutrinos produced in UHECR propagation can reach the Earth even from very high redshifts, and electromagnetic cascades initiated by secondary photons and electrons/positrons contribute to the diffuse gamma-ray background. Therefore, a multi-messenger analysis combining UHECR, neutrino, and gamma-ray data can still provide information about the cosmological evolution of UHECR sources. In this work, we compare predicted particle fluxes from Monte Carlo simulations in various scenarios with recent experimental data, and discuss the conclusions that can be drawn about UHECR sources and their cosmological evolution.

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.

Ultra high energy cosmic rays (UHECRs) are expected to be accelerated in astrophysical sources and to travel through extragalactic space before hitting the Earth atmosphere. They interact both with the environment in the source and with the intergalactic photon fields they encounter, causing different processes at various scales depending on the photon energy in the nucleus rest frame. UHECR interactions are sensitive to uncertainties in the extragalactic background spectrum and in the photo-disintegration models.

If ultra-high-energy cosmic rays (UHECRs) have extragalactic origins, as is widely assumed to be the case at least for the majority of cosmic rays with energies above a few EeV, secondary neutrinos and photons can be expected to be produced during the propagation of UHECRs through intergalactic space via interactions with cosmic background photons. The fluxes of such secondary particles are strongly dependent on the redshift evolution of the emissivity (number density times luminosity) of UHECR sources. We show how cosmic rays, neutrinos, and gamma rays can potentially provide complementary information about UHECR source evolution.

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.

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.

We introduce the new version of SimProp, a Monte Carlo code for simulating the propagation of ultra-high energy cosmic rays in intergalactic space. This version, SimProp v2r3, allows the choice of many more models for the extragalactic background light spectrum and evolution and photodisintegration cross sections and branching ratios than previous versions of SimProp.

The results of simulations of the extragalactic propagation of ultra-high energy cosmic rays (UHECRs) have intrinsic uncertainties due to poorly known physical quantities and approximations used in the codes. We quantify the uncertainties in the simulated UHECR spectrum and composition due to different models for the extragalactic background light (EBL), different photodisintegration setups, approximations concerning photopion production and the use of different simulation codes. We discuss the results for several representative source scenarios with proton, nitrogen or iron at injection. For this purpose we used SimProp and CRPropa, two publicly available codes for Monte Carlo simulations of UHECR propagation. CRPropa is a detailed and extensive simulation code, while SimProp aims to achieve acceptable results using a simpler code. We show that especially the choices for the EBL model and the photodisintegration setup can have a considerable impact on the simulated UHECR spectrum and composition.

Ultra-high energy cosmic ray experimental data are now of very good statistical significance even in the region of the expected GZK feature. The identification of their sources requires sophisticate analysis of their propagation in the extragalactic space. When looking at the details of this propagation some unforeseen features emerge. We will discuss some of these "surprises".

Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.

Relativistic Invariance might be modified by Quantum Gravity effects. The interesting point which emerged in the last fifteen years is that remnants of possible Lorentz Invariance Violations could be present at energies much lower than their natural scale, and possibly affect Ultra High Energy Cosmic Rays phenomena. We discuss their status in the view of recent data from the Pierre Auger Observatory.

Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013