• ### Studying cosmological $\gamma$-ray propagation with the Cherenkov Telescope Array(1709.04185)

Sept. 26, 2017 astro-ph.HE
The measurement of $\gamma$-rays originating from active galactic nuclei offers the unique opportunity to study the propagation of very-high-energy photons over cosmological distances. Most prominently, $\gamma$-rays interact with the extragalactic background light (EBL) to produce $e^+e^-$ pairs, imprinting an attenuation signature on $\gamma$-ray spectra. The $e^+e^-$ pairs can also induce electromagnetic cascades whose detectability in $\gamma$-rays depends on the intergalactic magnetic field (IGMF). Furthermore, physics beyond the Standard Model such as Lorentz invariance violation (LIV) or oscillations between photons and weakly interacting sub-eV particles (WISPs) could affect the propagation of $\gamma$-rays. The future Cherenkov Telescope Array (CTA), with its unprecedented $\gamma$-ray source sensitivity, as well as enhanced energy and spatial resolution at very high energies, is perfectly suited to study cosmological effects on $\gamma$-ray propagation. Here, we present first results of a study designed to realistically assess the capabilities of CTA to probe the EBL, IGMF, LIV, and WISPs.
• ### The Giant Radio Array for Neutrino Detection (GRAND): Present and Perspectives(1708.05128)

Aug. 17, 2017 astro-ph.IM, astro-ph.HE
The Giant Radio Array for Neutrino Detection (GRAND) aims at detecting ultra-high energy extraterrestrial neutrinos via the extensive air showers induced by the decay of tau leptons created in the interaction of neutrinos under the Earth's surface. Consisting of an array of $\sim10^5$ radio antennas deployed over $\sim 2\times10^5\,\rm {km}^2$, GRAND plans to reach, for the first time, an all-flavor sensitivity of $\sim1.5\times10^{-10} \,\rm GeV\, cm^{-2} \,s^{-1}\, sr^{-1}$ above $5\times10^{17}$ eV and a sub-degree angular resolution, beyond the reach of other planned detectors. We describe here preliminary designs and simulation results, plans for the ongoing, staged approach to the construction of GRAND, and the rich research program made possible by GRAND's design sensitivity and angular resolution.
• ### 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.
• ### Morphological properties of blazar-induced gamma-ray haloes(1707.02376)

July 8, 2017 astro-ph.CO, astro-ph.HE
At TeV energies and above gamma rays can induce electromagnetic cascades, whose charged component is sensitive to intervening intergalactic magnetic fields (IGMFs). When interpreting gamma-ray measurements in the energy range between a few GeV and hundreds of TeV, one has to carefully account for effects due to IGMFs, which depend on their strength and power spectrum. Therefore, gamma-ray-induced electromagnetic cascades can be used as probes of cosmic magnetism, since their arrival distribution as well as spectral and temporal properties can provide unique information about IGMFs, whose origin and properties are currently poorly understood. In this contribution we present an efficient three-dimensional Monte Carlo code for simulations of gamma-ray propagation. We focus on the effects of different configurations of IGMFs, in particular magnetic helicity and the power spectrum of stochastic fields, on the morphology of the arrival directions of gamma rays, and discuss the prospects for detecting pair haloes around distant blazars.
• ### Ultra-high-energy cosmic rays from tidally-ignited white dwarfs(1702.06978)

Ultra-high-energy cosmic rays (UHECRs) can be accelerated by tidal disruption events of stars by black holes. We suggest a novel mechanism for UHECR acceleration wherein white dwarfs (WDs) are tidally compressed by intermediate-mass black holes (IMBHs), leading to their ignition and subsequent explosion as a supernova. Cosmic rays accelerated by the supernova may receive an energy boost when crossing the accretion-powered jet. The rate of encounters between WDs and IMBHs can be relatively high, as the number of IMBHs may be substantially augmented once account is taken of their likely presence in dwarf galaxies. In this Letter, we show that this kind of tidal disruption event naturally provides an intermediate/heavy composition for the observed UHECRs, and suggest that dwarf galaxies and globular clusters are suitable sites for particle acceleration to ultra-high energies.
• ### Implications of Strong Intergalactic Magnetic Fields for Ultra-High-Energy Cosmic-Ray Astronomy(1704.05869)

June 14, 2017 astro-ph.CO, astro-ph.HE
We study the propagation of ultra-high-energy cosmic rays in the magnetised cosmic web. We focus on the particular case of highly magnetised voids ($B \sim \text{nG}$), using the upper bounds from the Planck satellite. The cosmic web was obtained from purely magnetohydrodynamical cosmological simulations of structure formation considering different power spectra for the seed magnetic field in order to account for theoretical uncertainties. We investigate the impact of these uncertainties on the propagation of cosmic rays, showing that they can affect the measured spectrum and composition by up to $\simeq 80\%$ and $\simeq 5\%$, respectivelly. In our scenarios, even if magnetic fields in voids are strong, deflections of 50 EeV protons from sources closer than $\sim\;$50 Mpc are less than $15^\circ$ in approximately 10-50% of the sky, depending on the distribution of sources and magnetic power spectrum. Therefore, UHECR astronomy might be possible in a significant portion of the sky depending on the primordial magnetic power spectrum, provided that protons constitute a sizeable fraction of the observed UHECR flux.
• ### 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.
• ### Electromagnetic Cascades as Probes of Cosmic Magnetism(1704.05137)

April 14, 2017 astro-ph.HE
The existence of intergalactic magnetic fields (IGMFs) is an open problem in cosmology and has never been unambiguously confirmed. High-energy gamma rays emitted by blazars are unique probes of cosmic magnetism, as their interactions with pervasive radiation fields generate a short-lived charged component sensitive to intervening magnetic fields. Spatial and temporal properties of the secondary gamma rays generated in the electromagnetic cascade can provide information about the strength, power spectrum, and topology of IGMFs. To probe these fields, detailed simulations of gamma-ray propagation in the intergalactic medium are necessary. In this work the effects of magnetic fields on the spectrum and arrival directions of gamma rays are studied using three-dimensional simulations, emphasising the particular case of helical IGMFs.
• ### Probing Intergalactic Magnetic Fields with Simulations of Electromagnetic Cascades(1607.00320)

July 1, 2016 astro-ph.HE
We determine the effect of intergalactic magnetic fields on the distribution of high energy gamma rays by performing three-dimensional Monte Carlo simulations of the development of gamma-ray-induced electromagnetic cascades in the magnetized intergalactic medium. We employ the so-called 'Large Sphere Observer' method to efficiently simulate blazar gamma ray halos. We study magnetic fields with a Batchelor spectrum and with maximal left- and right-handed helicities. We also consider the case of sources whose jets are tilted with respect to the line of sight. We verify the formation of extended gamma ray halos around the source direction, and observe spiral-like patterns if the magnetic field is helical. We apply the $Q$-statistics to the simulated halos to extract their spiral nature and also propose an alternative method, the $S$-statistics. Both methods provide a quantative way to infer the helicity of the intervening magnetic fields from the morphology of individual blazar halos for magnetic field strengths $B \gtrsim 10^{-15}\,{\rm G}$ and magnetic coherence lengths $L_{\rm c} \gtrsim 100\,{\rm Mpc}$. We show that the $S$-statistics has a better performance than the $Q$-statistics when assessing magnetic helicity from the simulated halos.
• 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.
• ### Diffusion of cosmic rays at EeV energies in inhomogeneous extragalactic magnetic fields(1407.6150)

Nov. 21, 2014 astro-ph.CO, astro-ph.HE
Ultra-high energy cosmic rays can propagate diffusively in cosmic magnetic fields. When their propagation time is comparable to the age of the universe, a suppression in the flux relative to the case in the absence of magnetic fields will occur. In this work we find an approximate parametrization for this suppression for energies below $\sim$ Z EeV using several magnetic field distributions obtained from cosmological simulations of the magnetized cosmic web. We assume that the magnetic fields have a Kolmogorov power spectrum with the field strengths distributed according to these simulations. We show that, if magnetic fields are coupled to the matter distribution, low field strengths will fill most of the volume, making the suppression milder compared to the case of a constant magnetic field with strength equal to the mean value of this distribution. We also derive upper limits for this suppression to occur for some models of extragalactic magnetic fields, as a function of the coherence length of these fields.
• ### Magnetic horizons of ultra-high energy cosmic rays(1410.7800)

Oct. 28, 2014 astro-ph.CO, astro-ph.HE
The propagation of ultra-high energy cosmic rays in extragalactic magnetic fields can be diffusive, depending on the strength and properties of the fields. In some cases the propagation time of the particles can be comparable to the age of the universe, causing a suppression in the flux measured on Earth. In this work we use magnetic field distributions from cosmological simulations to assess the existence of a magnetic horizon at energies around 10$^{18}$ 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.
• ### Propagation of UHECRs in the Universe(1308.1530)

March 17, 2014 astro-ph.CO, astro-ph.HE
The origin, propagation, and mechanisms of acceleration of the ultra-high energy cosmic rays (UHECRs) are not yet well understood. Aiming for a better interpretation of the available experimental data, these data have to be confronted with theoretical models. A realistic simulation of the propagation of UHECRs in the universe should take into account all the relevant energy loss processes due to the interaction with astrophysical backgrounds, as well as the intervening cosmic magnetic fields. Cosmological effects, such as the redshift dependence of the photon backgrounds and the adiabatic expansion of the universe can play an important role in the aforementioned processes. Here we present results of simulations of the propagation of UHECR through the large scale structure of the universe considering cosmological and magnetic field effects simultaneously.
• 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.
• ### Identifying Patterns on Cosmic Ray Maps with Wavelets on the Sphere(1201.2183)

Jan. 24, 2012 astro-ph.GA, astro-ph.IM
The deflection of ultra-high energy cosmic rays depends on the shape of the injection spectrum of the source and the pervasive cosmic magnetic fields. In this work it is applied the wavelet transform on the sphere to search for energy ordered filamentary structures arisen from magnetic bending. These structures, the so-called multiplets, can bring relevant information concerning the intervening magnetic fields.
• ### Detection of Point Sources in Cosmic Ray Maps using the Mexican Hat Wavelet Family(1201.2799)

Jan. 13, 2012 physics.data-an, astro-ph.IM
An analysis of the sensitivity of gaussian and mexican hat wavelet family filters to the detection of point sources of ultra-high energy cosmic rays was performed. A source embedded in a background was simulated and the number of events and amplitude of this source was varied aiming to check the sensitivity of the method to detect faint sources with low statistic of events.
• ### Amplification of the Signal-to-Noise Ratio in Cosmic Ray Maps Using the Mexican Hat Wavelet Family(1108.6253)

Aug. 31, 2011 physics.data-an, astro-ph.IM
In this work we analyze the effect of smoothing maps containing arrival directions of cosmic rays with a gaussian kernel and kernels of the mexican hat wavelets of orders 1, 2 and 3. The analysis is performed by calculating the amplification of the signal-to-noise ratio for several background patterns (noise) and different number of events coming from a simulated source (signal) for an ideal detector capable of observing the full sky with uniform coverage. We extend this analysis for a virtual observatory with two sites, one in the northern hemisphere, the other in the southern, considering an acceptance law.