Search

Active galactic nuclei (AGN) are among the main candidates for ultra-high-energy cosmic ray (UHECR) sources. However, while theoretical and some phenomenological works favor AGNs as the main sources, recent works have shown that using the very-high-energy (VHE) $\gamma$-ray flux as a proxy for the UHECR flux leads to a bad agreement with data. In this context, the energy spectrum and composition data are hardly fitted. At the same time, the arrival directions map is badly described and a spurious dipole direction is produced. In this work, we propose a possible solution to these contradictions. Using the observed $\gamma$-ray flux as a proxy may carry the implicit assumption of beamed UHECR emission and, consequently, its beam will remain collimated up to its detection on Earth. We show that assuming an isotropic UHECR emission and correcting the $\gamma$-ray emission proxy by Doppler boosting can overcome the problem. The combined fit of the spectrum and composition is improved by $3.56\sigma$, while the predicted arrival directions agree much better with the data. In particular, a spurious direction of the dipole can be reduced from $10.3 \ (5.4)\sigma$ away from the data to $2.2 \ (1.5)\sigma$ for $E > 8$ EeV ($E > 32$ EeV). We also show that this effect is particularly important when including AGNs of different classes in the same analysis, such as radio galaxies and blazars.

The predictions of hadronic interaction models for cosmic-ray induced air showers contain inherent uncertainties due to limitations of available accelerator data. This leads to differences in shower simulations using each of those models. Many studies have been carried out to track those differences by investigating the shower development or the particle content. In this work, we propose a new approach to search for discrepancies and similarities between the models, via the IACT images resulting from the observations of hadronic air showers. We use simulations of H.E.S.S. as a show-case scenario and, by investigating variables of the camera images, we find potential indicators to highlight differences between models. Number of pixels, Hillas image size, and density showed the largest difference between the models. We then further explore the (in)compatibility of the models by combining all the variables and using Boosted Decision Trees. For protons, a significant difference in the classifier output is found for EPOS-LHC when compared to both QGSJET-II04 and Sybill 2.3d. For helium and nitrogen, QGSJET-II04 is shown to be the outlier case. No significant differences are found for silicon and iron. The distribution of (in)compatibility between the models in the phase space of shower parameters shows that a targeted search can be fruitful, with showers with energies of a few TeV and core closer to the large telescope presenting the largest power of separation.

The Telescope Array experiment has recently reported the most energetic event detected in the hybrid technique era, with a reconstructed energy of 240 EeV, which has been named ``Amaterasu'' after the Shinto deity. Its origin is intriguing since no powerful enough candidate sources are located within the region consistent with its propagation horizon and arrival direction. In this work, we investigate the possibility of describing its origin in a scenario of new physics, specifically under a Lorentz Invariance Violation (LIV) assumption. The kinematics of UHECR propagation under a phenomenological LIV approach is investigated. The total mean free path for a particle with Amaterasu's energy increases from a few Mpc to hundreds of Mpc for $-\delta_{\rm{had},0} > 10^{-22}$, expanding significantly the region from which it could have originated. A combined fit of the spectrum and composition data of Telescope Array under different LIV assumptions was also performed. The data is best fitted with some level of LIV both with and without Amaterasu. Robustness with data from the Pierre Auger Observatory is investigated by exploring an intermediate composition scenario. Similar improvements in the description of the data with LIV are found for that. New physics in the form of LIV could, thus, provide a plausible and robust explanation for the Amaterasu particle.

Recent data from the Pierre Auger Observatory has revealed the presence of a large-scale dipole in the arrival direction distribution of ultra-high energy cosmic rays (UHECR). In this work, we build up an understanding of the diffusive origin of such a dipolar behavior as well as its dependency on energy and astrophysical source assumptions such as extragalactic magnetic field strength and cosmic ray composition. We present a novel analytical approach for calculating the angular distribution of CR coming from a single source and discuss the regimes in which the steady-state dipole result is expected. We also present a semianalytical method for calculating the evolution with energy of the resultant dipole for an ensemble of sources. We show that a local source allows for a strong growth of the dipole with energy over a large energy range. The possibility of a transition from a dipolar to non-dipolar regime at the highest energies and its implications for the source density, magnetic field intensity, and cosmic ray composition are discussed.

We update the constraints on the location of the nearest UHECR source. By analyzing recent data from the Pierre Auger Observatory using state-of-the-art CR propagation models, we reaffirm the need of local sources with a distance less than 25-100 Mpc, depending on mass composition. A new fast semi-analytical method for the propagation of UHECR in environments with turbulent magnetic fields is developed. The onset of an enhancement and a low-energy magnetic horizon of cosmic rays from sources located within a particular distance range is demonstrated. We investigate the distance to the nearest source, taking into account these magnetic field effects. The results obtained highlight the robustness of our constrained distances to the nearest source.

Lorentz invariance (LI) has a central role in science and its violation (LIV) at some high-energy scale has been related to possible solutions for several of the most intriguing puzzles in nature such as dark matter, dark energy, cosmic rays generation in extreme astrophysical objects and quantum gravity. We report on a search for LIV signal based on the propagation of gamma rays from astrophysical sources to Earth. An innovative data analysis is presented which allowed us to extract unprecedented information from the most updated data set composed of 111 energy spectra of 38 different sources measured by current gamma-ray observatories. No LIV signal was found, and we show that the data are best described by LI assumption. We derived limits for the LIV energy scale at least 3 times better than the ones currently available in the literature for subluminal signatures of LIV in high-energy gamma rays., Comment: Published in PRD

Testing Lorentz invariance is essential as it is one of the pillars of modern physics. Moreover, its violation is foreseen in several popular Quantum Gravity models. Several authors study the effects of Lorentz invariance violation (LIV) in the propagation of ultra-high energy cosmic rays. These particles are the most energetic events ever detected and therefore represent a promising framework to test LIV. In this work we present an analytic calculation of the in-elasticity for any $a+b \rightarrow c+d$ interaction using first order perturbation in the dispersion relation that violates Lorentz invariance. The inelasticity can be calculated by solving a third-order polynomial equation containing: a) the kinematics of the interaction, b) the LIV term for each particle and c) the geometry of the interaction. We use the inelasticity we calculate to investigate the proton propagation in the intergalactic media. The photopion production of the proton interaction with the CMB is taken into account using the inelasticity and the attenuation length in different LIV scenarios. We show how the allowed phase space for the photopion production changes when LIV is considered for the interaction. The calculations presented here are going to be extended in order to calculated the modified ultra-high energy cosmic rays spectrum and compare it to the data.

In this paper, Lorentz Invariance Violation (LIV) is introduced in the calculations of photon propagation in the Universe. LIV is considered in the photon sector and the mean free path of the $\gamma \gamma \rightarrow e^{+} e^{-}$ interaction is calculated. The corresponding photon horizon including LIV effects is used to predict major changes in the propagation of photons with energy above $10^{18}$ eV. The flux of GZK photons on Earth considering LIV is calculated for several source models of ultra-high energy cosmic ray (UHECR). The predicted flux of GZK gamma-rays is compared to the new upper limits on the photon flux obtained by the Pierre Auger Observatory in order to impose upper limits on the LIV coefficients of order $n =$ 0, 1 and 2. The limits on the LIV coefficients derived here are more realistic than previous works and in some cases more restrictive. The analysis resulted in LIV upper limits in the photon sector of $\delta_{\gamma,0}^{limit} \sim -10^{-20}$, $\delta_{\gamma,1}^{limit} \sim -10^{-38} \; \mathrm{eV^{-1}}$ and $\delta_{\gamma,2}^{limit} \sim -10^{-56} \; \mathrm{eV^{-2}}$ in the astrophysical scenario which best describes UHECR data.