• ### IMAGINE: A comprehensive view of the interstellar medium, Galactic magnetic fields and cosmic rays(1805.02496)

May 7, 2018 astro-ph.GA
In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fields and their interactions with cosmic rays. To achieve a higher level of self-consistency, depth and rigour can only be achieved by the coordinated efforts of experts in diverse areas of astrophysics involved in observational, theoretical and numerical work. We present our view of the present status of this topic, identify its key unsolved problems and suggest a strategy that will underpin our work. The backbone of the consortium is the Interstellar MAGnetic field INference Engine, a publicly available Bayesian platform that employs robust statistical methods to explore the multi-dimensional likelihood space using any number of modular inputs. It provides an interpretation and modelling framework that has the power and flexibility to include a variety of observational, theoretical and numerical lines of evidence into a self-consistent and comprehensive picture of the thermal and non-thermal interstellar media. An important innovation is that a consistent understanding of the phenomena that are directly or indirectly influenced by the Galactic magnetic field, such as the deflection of ultra-high energy cosmic rays or extragalactic backgrounds, is made an integral part of the modelling. The IMAGINE Consortium, which is informal by nature and open to new participants, hereby presents a methodological framework for the modelling and understanding of Galactic magnetic fields that is available to all communities whose research relies on a state-of-the-art solution to this problem. (Abridged.)
• ### Ultra-High-Energy Cosmic Rays from Radio Galaxies(1701.06792)

March 8, 2018 astro-ph.HE
Radio galaxies are intensively discussed as the sources of cosmic rays observed above about $3\,{\times}\,10^{18}\,\text{eV}$, called ultra-high energy cosmic rays (UHECRs). We present a first, systematic approach that takes the individual characteristics of these sources into account, as well as the impact of the extragalactic magnetic-field structures up to a distance of 120 Mpc. We use a mixed simulation setup, based on 3D simulations of UHECRs ejected by observed, individual radio galaxies taken out to a distance of 120 Mpc, and on 1D simulations over a continuous source distribution contributing from beyond 120 Mpc. Additionally, we include the ultra-luminous radio galaxy Cygnus A at a distance of about $250\,$Mpc, as its contribution is so strong that it must be considered as an individual point source. The implementation of the UHECR ejection in our simulation setup is based on a detailed consideration of the physics of radio jets and standard first-order Fermi acceleration. We show that the average contribution of radio galaxies taken over a very large volume cannot explain the observed features of UHECRs measured at Earth. However, we obtain excellent agreement with the spectrum, composition, and arrival-direction distribution of UHECRs measured by the Pierre Auger Observatory, if we assume that most UHECRs observed arise from only two sources: The ultra-luminous radio galaxy Cygnus A, providing a mostly light composition of nuclear species dominating up to about $6\,{\times}\,10^{19}\,$eV, and the nearest radio galaxy Centaurus A, providing a heavy composition dominating above $6\,{\times}\,10^{19}\,$eV. Here we have to assume that extragalactic magnetic fields out to 250 Mpc, which we did not include in the simulation, are able to isotropize the UHECR events at about 8 EeV arriving from Cygnus A.
• ### Cosmogenic gamma-rays and neutrinos constrain UHECR source models(1707.04511)

July 14, 2017 astro-ph.HE
We use CRPropa 3 to show how the expected cosmogenic neutrino and gamma-ray spectra depend on the maximum energy of ultra-high energy cosmic rays (UHECRs) at their sources, on the spectral index at injection and on the chemical composition of UHECRs. The isotropic diffuse gamma-ray background measured by Fermi/LAT is already close to touching upon a model with co-moving source evolution and with the chemical composition, spectral index and maximum acceleration energy optimized to provide the best fit to the UHECR spectrum and composition measured by the Pierre Auger Collaboration. Additionally, the detectable fraction of protons present at the highest energies in UHECRs, for experiments with sensitivities to the single-flavor neutrino flux at $\sim1$ EeV in the range of $\sim 10^{-8}$ - $10^{-10}$ GeV cm$^{-2}$ s$^{-1}$ sr$^{-1}$, is shown as a function of the evolution of UHECR sources. Experiments that reach this sensitivity will be able to significantly constrain the proton fraction for realistic source evolution models.
• ### Extragalactic Sources and Propagation of UHECRs(1705.02293)

May 5, 2017 astro-ph.HE
With the publicly available astrophysical simulation framework for propagating extraterrestrial UHE particles, CRPropa 3, it is now possible to study realistic UHECR source scenarios including deflections in Galactic and extragalactic magnetic fields in an efficient way. Here we discuss three recent studies that have already been done in that direction. The first one investigates what can be expected in the case of maximum allowed intergalactic magnetic fields. Here is shown that, even if voids contain strong magnetic fields, deflections of protons with energies $\gtrsim 60 \; \text{EeV}$ from nearby sources might be small enough to allow for UHECR astronomy. The second study looks into several scenarios with a smaller magnetization focusing on large-scale anisotropies. Here is shown that the local source distribution can have a more significant effect on the large-scale anisotropy than the EGMF model. A significant dipole component could, for instance, be explained by a dominant source within 5 Mpc distance. The third study looks into whether UHECRs can come from local radio galaxies. If this is the case it is difficult to reproduce the observed low level of anisotropy. Therefore is concluded that the magnetic field strength in voids in the EGMF model used here is too low and/or there are additional sources of UHECRs that were not taken into account in these simulations.
• ### Cosmogenic photons strongly constrain UHECR source models(1609.03336)

Sept. 12, 2016 astro-ph.HE
With the newest version of our Monte Carlo code for ultra-high-energy cosmic ray (UHECR) propagation, CRPropa 3, the flux of neutrinos and photons due to interactions of UHECRs with extragalactic background light can be predicted. Together with the recently updated data for the isotropic diffuse gamma-ray background (IGRB) by Fermi LAT, it is now possible to severely constrain UHECR source models. The evolution of the UHECR sources especially plays an important role in the determination of the expected secondary photon spectrum. Pure proton UHECR models are already strongly constrained, primarily by the highest energy bins of Fermi LAT's IGRB, as long as their number density is not strongly peaked at recent times.
• 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.
• ### PeV Neutrinos from Ultra-High-Energy Cosmic Rays(1407.6577)

July 24, 2014 astro-ph.HE
We investigate the possibility that the recently detected TeV-PeV neutrino events by IceCube can originate from extragalactic ultra-high-energy cosmic ray interactions with the cosmic microwave background or the UV/optical/IR background. This is done by simulating the propagation of the cosmic rays from their sources to the observer, including the production and propagation of secondary neutrinos and gamma rays. For this purpose we use the publicly available simulation package CRPropa 2.0. We find that in all the scenarios considered here the simulated neutrino flux level remains at least one order of magnitude below the flux level indicated by the IceCube events, thus showing it is difficult to interpret the IceCube events in terms of a cosmogenic neutrino flux.
• 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.
• ### PeV neutrinos from the propagation of ultra-high energy cosmic rays(1209.4033)

Jan. 23, 2013 hep-ph, astro-ph.HE
We discuss the possibility that the PeV neutrinos recently observed by IceCube are produced by the interactions of extragalactic cosmic rays during their propagation through the radiation backgrounds. We show that the fluxes resulting from the decays of neutrons produced in the interactions of cosmic ray protons with the CMB background are suppressed ($E_\nu^2$d$\Phi_\nu/$d$E< 10^{-10}$ GeV/cm$^2$ s sr), with those resulting from the decays of pions produced in the interactions with the UV/optical/IR backgrounds being the dominant ones at PeV energies. The anti-neutrino fluxes produced by the decay of neutrons resulting from the photodisintegration of heavy nuclei with CMB photons are also shown to be quite suppressed ($E_\nu^2$d$\Phi_\nu/$d$E< 10^{-11}$ GeV/cm$^2$ s sr), while those produced by photo-pion processes with UV/optical/IR backgrounds may be larger, although they are not expected to be above those achievable in the pure proton case. Scenarios with mixed composition and low cutoff rigidities can lead to PeV neutrino fluxes enhanced with respect to those in the pure Fe scenarios. We also discuss the possible impact of the Glashow resonance for the detection of these scenarios, showing that it plays a moderate role.
• ### Determining neutrino absorption spectra at Ultra-High Energies(0801.3342)

April 16, 2008 astro-ph
A very efficient method to measure the flux of Ultra-high energy (UHE) neutrinos is through the detection of radio waves which are emitted by the particle shower in the lunar regolith. The highest acceptance is reached for radio waves in the frequency band of 100-200 MHz which can be measured with modern radio telescopes. In this work we investigate the sensitivity of this detection method to structures in the UHE neutrino spectrum caused by their absorption on the low-energy relic anti-neutrino background through the Z-boson resonance. The position of the absorption peak is sensitive to the neutrino mass and the redshift of the source. A new generation of low-frequency digital radio telescopes will provide excellent detection capabilities for measuring these radio pulses, thus making our consideration here very timely.