• Energy Spectrum of Fast Second Order Fermi Accelerators as Sources of Ultra-High-Energy Cosmic Rays(1612.03675)

April 18, 2018 astro-ph.HE
Stochastic acceleration of cosmic rays in second order Fermi processes is usually considered too slow to reach ultra-high energies, except in specific cases. In this paper we present the energy spectrum obtained from second order Fermi acceleration in highly turbulent magnetic fields as e.g.\ found in the outskirts of AGN jets in situations where it can be sufficiently fast to accelerate particles to the highest observed energies. We parametrize the resulting non-power-law spectra and show that these can describe the cosmic ray energy spectrum and mass-composition data at the highest energies if propagation effects are taken into account.
• Overview of lunar detection of ultra-high energy particles and new plans for the SKA(1704.05336)

April 18, 2017 astro-ph.IM, astro-ph.HE
The lunar technique is a method for maximising the collection area for ultra-high-energy (UHE) cosmic ray and neutrino searches. The method uses either ground-based radio telescopes or lunar orbiters to search for Askaryan emission from particles cascading near the lunar surface. While experiments using the technique have made important advances in the detection of nanosecond-scale pulses, only at the very highest energies has the lunar technique achieved competitive limits. This is expected to change with the advent of the Square Kilometre Array (SKA), the low-frequency component of which (SKA-low) is predicted to be able to detect an unprecedented number of UHE cosmic rays. In this contribution, the status of lunar particle detection is reviewed, with particular attention paid to outstanding theoretical questions, and the technical challenges of using a giant radio array to search for nanosecond pulses. The activities of SKA's High Energy Cosmic Particles Focus Group are described, as is a roadmap by which this group plans to incorporate this detection mode into SKA-low observations. Estimates for the sensitivity of SKA-low phases 1 and 2 to UHE particles are given, along with the achievable science goals with each stage. Prospects for near-future observations with other instruments are also described.
• Ultimate precision in cosmic-ray radio detection --- the SKA(1608.08869)

Feb. 6, 2017 hep-ex, astro-ph.IM, astro-ph.HE
As of 2023, the low-frequency part of the Square Kilometre Array will go online in Australia. It will constitute the largest and most powerful low-frequency radio-astronomical observatory to date, and will facilitate a rich science programme in astronomy and astrophysics. With modest engineering changes, it will also be able to measure cosmic rays via the radio emission from extensive air showers. The extreme antenna density and the homogeneous coverage provided by more than 60,000 antennas within an area of one km$^2$ will push radio detection of cosmic rays in the energy range around 10$^{17}$ eV to ultimate precision, with superior capabilities in the reconstruction of arrival direction, energy, and an expected depth-of-shower-maximum resolution of 6~g/cm${^2}$.
• 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.
• 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.
• Prospects of GPGPU in the Auger Offline Software Framework(1507.07733)

July 28, 2015 astro-ph.IM
The Pierre Auger Observatory is the currently largest experiment dedicated to unveil the nature and origin of the highest energetic cosmic rays. The software framework 'Offline' has been developed by the Pierre Auger Collaboration for joint analysis of data from different independent detector systems used in one observatory. While reconstruction modules are specific to the Pierre Auger Observatory components of the Offline framework are also used by other experiments. The software framework has recently been extended to incorporate data from the Auger Engineering Radio Array (AERA), the radio extension of the Pierre Auger Observatory. The reconstruction of the data of such radio detectors requires the repeated evaluation of complex antenna gain patterns which significantly increases the required computing resources in the joint analysis. In this contribution we explore the usability of massive parallelization of parts of the Offline code on the GPU. We present the result of a systematic profiling of the joint analysis of the Offline software framework aiming for the identification of code areas suitable for parallelization on GPUs. Possible strategies and obstacles for the usage of GPGPU in an existing experiment framework are discussed.
• 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.
• Detecting Local Deflection Patterns of Ultra-high Energy Cosmic Rays using the Principal Axes of the Directional Energy Distribution(1307.8273)

July 31, 2013 astro-ph.IM, astro-ph.HE
From deflections in galactic and extragalactic magnetic fields energy dependent structures in the arrival directions of ultra-high energy cosmic rays (UHECR) are expected. We propose to characterize these structures by the strength of collimation of energy along the principal axes in selected regions in the sky. While the strength of collimation are indicators of anisotropy in the arrival distribution of UHECR, the orientation of the principal system holds information about the direction of the deflections of UHECR. We discuss the method and present expected limits on the strength of deflection and density of sources using simulated scenarios of UHECR proton propagation.
• 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.