A promising energy range to look for angular correlations between cosmic rays of extragalactic origin and their sources is at the highest energies, above a few tens of EeV (1 EeV ≡ 1018 eV). Despite the flux of these particles being extremely low, the area of ∼3000 km2 covered at the Pierre Auger Observatory, and the 17 yr data-taking period of the Phase 1 of its operations, have enabled us to measure the arrival directions of more than 2600 ultra-high-energy cosmic rays above 32 EeV. We publish this data set, the largest available at such energies from an integrated exposure of 122,000 km2 sr yr, and search it for anisotropies over the 3.4π steradians covered with the Observatory. Evidence for a deviation in excess of isotropy at intermediate angular scales, with ∼15° Gaussian spread or ∼25° top-hat radius, is obtained at the 4σ significance level for cosmic-ray energies above ∼40 EeV.

Ultra-high-energy photons with energies exceeding 1017 eV offer a wealth of connections to different aspects of cosmic-ray astrophysics as well as to gamma-ray and neutrino astronomy. The recent observations of photons with energies in the 1015 eV range further motivate searches for even higher-energy photons. In this paper, we present a search for photons with energies exceeding 2 × 1017 eV using about 5.5 yr of hybrid data from the low-energy extensions of the Pierre Auger Observatory. The upper limits on the integral photon flux derived here are the most stringent ones to date in the energy region between 1017 and 1018 eV.

No abstract

For several decades, the origin of ultra-high-energy cosmic rays (UHECRs) has been an unsolved question of high-energy astrophysics. One approach for solving this puzzle is to correlate UHECRs with high-energy neutrinos, since neutrinos are a direct probe of hadronic interactions of cosmic rays and are not deflected by magnetic fields. In this paper, we present three different approaches for correlating the arrival directions of neutrinos with the arrival directions of UHECRs. The neutrino data are provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above ∼50 EeV are provided by the Pierre Auger Observatory and the Telescope Array. All experiments provide increased statistics and improved reconstructions with respect to our previous results reported in 2015. The first analysis uses a high-statistics neutrino sample optimized for point-source searches to search for excesses of neutrino clustering in the vicinity of UHECR directions. The second analysis searches for an excess of UHECRs in the direction of the highest-energy neutrinos. The third analysis searches for an excess of pairs of UHECRs and highest-energy neutrinos on different angular scales. None of the analyses have found a significant excess, and previously reported overfluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

Lorentz invariance violation (LIV) is often described by dispersion relations of the form E i 2 = m i 2+p i 2+δi,n E 2+n with delta different based on particle type i, with energy E, momentum p and rest mass m. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constraints on the LIV coefficients δi,n tend to come from the highest energies. At sufficiently high energies, photons produced by cosmic ray interactions as they propagate through the Universe could be subluminal and unattenuated over cosmological distances. Cosmic ray interactions can also be modified and lead to detectable fingerprints in the energy spectrum and mass composition observed on Earth. The data collected at the Pierre Auger Observatory are therefore possibly sensitive to both the electromagnetic and hadronic sectors of LIV. In this article, we explore these two sectors by comparing the energy spectrum and the composition of cosmic rays and the upper limits on the photon flux from the Pierre Auger Observatory with simulations including LIV. Constraints on LIV parameters depend strongly on the mass composition of cosmic rays at the highest energies. For the electromagnetic sector, while no constraints can be obtained in the absence of protons beyond 1019 eV, we obtain δγ,0 > -10-21, δγ,1 > -10-40 eV-1 and δγ,2 > -10-58 eV-2 in the case of a subdominant proton component up to 1020 eV. For the hadronic sector, we study the best description of the data as a function of LIV coefficients and we derive constraints in the hadronic sector such as δhad,0 < 10-19, δhad,1 < 10-38 eV-1 and δhad,2 < 10-57 eV-2 at 5σ CL.