R. Alves Batista

Effects of uncertainties in simulations of extragalactic UHECR propagation, using CRPropa and SimProp

The results of simulations of 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 of 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.

Cosmic ray propagation with CRPropa 3

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 ~ 1017 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.

Using spherical wavelets to search for magnetically-induced alignment in the arrival directions of ultra-high energy cosmic rays

Abstract Due to the action of the intervening cosmic magnetic fields, ultra-high energy cosmic rays (UHECRs) can be deflected in such a way as to create clustered energy-ordered filamentary structures in the arrival direction of these particles, the so-called multiplets. In this work we propose a new method based on the spherical wavelet transform to identify multiplets in sky maps containing arrival directions of UHECRs. The method is illustrated in simulations with a multiplet embedded in isotropic backgrounds with different numbers of events. The efficiency of the algorithm is assessed through the calculation of Type I and II errors.

CRPropa: A public framework to propagate UHECRs in the universe

To answer the fundamental questions concerning the origin and nature of ultra-high energy cosmic rays (UHECRs), it is important to confront data with simulated astrophysical scenarios. These scenarios should include detailed information on particle interactions and astrophysical environments. To achieve this goal one should make use of computational tools to simulate the propagation of these particles. For this reason the CRPropa framework was developed. It allows the propagation of UHECRs with energies ≳1017  eV and secondary gamma rays and neutrinos. The newest version, CRPropa 3, reflects an efficient redesign of the code as well as several new features such as time dependent propagation in three dimensions, galactic magnetic field effects and improved treatment of interactions, among other enhancements.