• ### Production and propagation of ultra-high energy photons using CRPropa 3(1710.11406)

Jan. 5, 2018 astro-ph.IM, astro-ph.HE
In order to interpret cosmic ray observations, detailed modeling of propagation effects invoking all important messengers is necessary. We introduce a new photon production and propagation code as an inherent part of the CRPropa 3 software framework. By implementing additional photon production channels, which are important for energies below 10**18 eV, this code can be used for multi-messenger studies connecting the TeV and sub EeV energy regime and for interpreting models of ultra-high energy cosmic ray sources. We discuss the importance of the individual production channels and propagation effects and present example applications.
• ### A Deep Learning-based Reconstruction of Cosmic Ray-induced Air Showers(1708.00647)

Aug. 2, 2017 astro-ph.IM, astro-ph.HE
We describe a method of reconstructing air showers induced by cosmic rays using deep learning techniques. We simulate an observatory consisting of ground-based particle detectors with fixed locations on a regular grid. The detector's responses to traversing shower particles are signal amplitudes as a function of time, which provide information on transverse and longitudinal shower properties. In order to take advantage of convolutional network techniques specialized in local pattern recognition, we convert all information to the image-like grid of the detectors. In this way, multiple features, such as arrival times of the first particles and optimized characterizations of time traces, are processed by the network. The reconstruction quality of the cosmic ray arrival direction turns out to be competitive with a standard reconstruction algorithm. The reconstructed shower direction, energy and shower depth show the expected dependencies on the cosmic ray energy.
• 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.
• ### CRPropa 3 - a Public Astrophysical Simulation Framework for Propagating Extraterrestrial Ultra-High Energy Particles(1603.07142)

We present the simulation framework CRPropa version 3 designed for efficient development of astrophysical predictions for ultra-high energy particles. Users can assemble modules of the most relevant propagation effects in galactic and extragalactic space, include their own physics modules with new features, and receive on output primary and secondary cosmic messengers including nuclei, neutrinos and photons. In extension to the propagation physics contained in a previous CRPropa version, the new version facilitates high-performance computing and comprises new physical features such as an interface for galactic propagation using lensing techniques, an improved photonuclear interaction calculation, and propagation in time dependent environments to take into account cosmic evolution effects in anisotropy studies and variable sources. First applications using highlighted features are presented as well.
• ### Effects of uncertainties in simulations of extragalactic UHECR propagation, using CRPropa and SimProp(1508.01824)

Dec. 13, 2015 astro-ph.IM, astro-ph.HE
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.
• 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.
• ### Cosmic ray propagation with CRPropa 3(1410.5323)

Oct. 20, 2014 astro-ph.IM, astro-ph.HE
Solving the question of the origin of ultra-high energy cosmic rays (UHECRs) requires the development of detailed simulation tools in order to interpret the experimental data and draw conclusions on the UHECR universe. CRPropa is a public Monte Carlo code for the galactic and extragalactic propagation of cosmic ray nuclei above $\sim 10^{17}$ eV, as well as their photon and neutrino secondaries. In this contribution the new algorithms and features of CRPropa 3, the next major release, are presented. CRPropa 3 introduces time-dependent scenarios to include cosmic evolution in the presence of cosmic ray deflections in magnetic fields. The usage of high resolution magnetic fields is facilitated by shared memory parallelism, modulated fields and fields with heterogeneous resolution. Galactic propagation is enabled through the implementation of galactic magnetic field models, as well as an efficient forward propagation technique through transformation matrices. To make use of the large Python ecosystem in astrophysics CRPropa 3 can be steered and extended in Python.
• ### A field study of data analysis exercises in a bachelor physics course using the internet platform VISPA(1402.2836)

Feb. 12, 2014 hep-ex, physics.ed-ph, astro-ph.HE
Bachelor physics lectures on particle physics and astrophysics were complemented by exercises related to data analysis and data interpretation at the RWTH Aachen University recently. The students performed these exercises using the internet platform VISPA, which provides a development environment for physics data analyses. We describe the platform and its application within the physics course, and present the results of a student survey. The students acceptance of the learning project was positive. The level of acceptance was related to their individual preference for learning with a computer. Furthermore, students with good programming skills favor working individually, while students who attribute themselves having low programming abilities favor working in teams. The students appreciated approaching actual research through the data analysis tasks.
• ### PARSEC: A Parametrized Simulation Engine for Ultra-High Energy Cosmic Ray Protons(1302.3761)

Dec. 12, 2013 astro-ph.HE
We present a new simulation engine for fast generation of ultra-high energy cosmic ray data based on parametrizations of common assumptions of UHECR origin and propagation. Implemented are deflections in unstructured turbulent extragalactic fields, energy losses for protons due to photo-pion production and electron-pair production, as well as effects from the expansion of the universe. Additionally, a simple model to estimate propagation effects from iron nuclei is included. Deflections in galactic magnetic fields are included using a matrix approach with precalculated lenses generated from backtracked cosmic rays. The PARSEC program is based on object oriented programming paradigms enabling users to extend the implemented models and is steerable with a graphical user interface.
• Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
• ### CRPropa 3.0 - a Public Framework for Propagating UHE Cosmic Rays through Galactic and Extragalactic Space(1307.2643)

July 10, 2013 astro-ph.IM, astro-ph.HE
The interpretation of experimental data of ultra-high energy cosmic rays (UHECRs) above 10^17 eV is still under controversial debate. The development and improvement of numerical tools to propagate UHECRs in galactic and extragalactic space is a crucial ingredient to interpret data and to draw conclusions on astrophysical parameters. In this contribution the next major release of the publicly available code CRPropa (3.0) is presented. It reflects a complete redesign of the code structure to facilitate high performance computing and comprises new physical features such as an interface for galactic propagation using lensing techniques and inclusion of cosmological effects in a three-dimensional environment. The performance is benchmarked and first applications are presented.
• ### A new Monte Carlo Generator for Ultra-High Energy Cosmic Rays from the Local and Distant Universe(1202.3005)

Feb. 14, 2012 astro-ph.HE
For the understanding of the origin and propagation of ultra-high energy cosmic rays (UHECR) we developed a new approach to simulating UHECRs from an arbitrary number of sources based on Monte Carlo technique. The method consists of a combination of three steps. For distant sources we apply commonly accepted parameterizations to calculate the contribution to the observed cosmic ray flux. For sources of the local universe we use forward tracking through realistic matter distributions and magnetic fields resulting from explicit simulations of large-scale structure formation. From the calculations and the forward tracking we generate maps of the probability to observe a particle with a given energy from a discrete direction. To account for deflections in the galactic field, these probability maps are transformed by matrices calculated from backtracking of antiparticles through field parameterizations. Based on the combined probability maps, Monte Carlo production of individual UHECR data is performed which are then used in comparisons with UHECR measurements. The simulated UHECR data serves for optimizing the analysis techniques used in UHECR measurements as well as for constraining the parameter space of the underlying source and magnetic field models.