We report on the construction of the SiRO-SiPM ReadOut muon detector, a detection system based on plastic scintillator bars designed for muography applications. Using six 1 m(2) layers of active medium, grouped two by two into three rectangular matrices of pixels, each separated by a variable distance, the spatial coordinates of the muon's impact point on every matrice are obtained and used for trajectory reconstruction. Validation studies have been performed using Monte Carlo simulations and later confirmed by preliminary measurements in our laboratory and in underground, in the Slanic Prahova salt mine, in Romania.

In this work we present the interpretation of the energy spectrum and mass composition data as measured by the Pierre Auger Collaboration above 6 × 10^17 eV. We use an astrophysical model with two extragalactic source populations to model the hardening of the cosmic-ray flux at around 5× 10^18 eV (the so-called "ankle" feature) as a transition between these two components. We find our data to be well reproduced if sources above the ankle emit a mixed composition with a hard spectrum and a low rigidity cutoff. The component below the ankle is required to have a very soft spectrum and a mix of protons and intermediate-mass nuclei. The origin of this intermediate-mass component is not well constrained and it could originate from either Galactic or extragalactic sources. To the aim of evaluating our capability to constrain astrophysical models, we discuss the impact on the fit results of the main experimental systematic uncertainties and of the assumptions about quantities affecting the air shower development as well as the propagation and redshift distribution of injected ultra-high-energy cosmic rays (UHECRs).

A catalog containing details of the highest-energy cosmic rays recorded through the detection of extensive air-showers at the Pierre Auger Observatory is presented with the aim of opening the data to detailed examination. Descriptions of the 100 showers created by the highest-energy particles recorded between 1 January 2004 and 31 December 2020 are given for cosmic rays that have energies in the range 78 EeV to 166 EeV. Details are also given of a further nine very-energetic events that have been used in the calibration procedure adopted to determine the energy of each primary. A sky plot of the arrival directions of the most energetic particles is shown. No interpretations of the data are offered.

Instantons, which are nonperturbative solutions to Yang-Mills equations, provide a signal for the occurrence of quantum tunneling between distinct classes of vacua. They can give rise to decays of particles otherwise forbidden. Using data collected at the Pierre Auger Observatory, we search for signatures of such instanton-induced processes that would be suggestive of super-heavy particles decaying in the Galactic halo. These particles could have been produced during the post-inflationary epoch and match the relic abundance of dark matter inferred today. The nonobservation of the signatures searched for allows us to derive a bound on the reduced coupling constant of gauge interactions in the dark sector: α_{X}≲0.09, for 10^{9}≲M_{X}/GeV<10^{19}. Conversely, we obtain that, for instance, a reduced coupling constant α_{X}=0.09 excludes masses M_{X}≳3×10^{13}  GeV. In the context of dark matter production from gravitational interactions alone, we illustrate how these bounds are complementary to those obtained on the Hubble rate at the end of inflation from the nonobservation of tensor modes in the cosmological microwave background.

Using the data of the Pierre Auger Observatory, we report on a search for signatures that would be suggestive of super-heavy particles decaying in the Galactic halo. From the lack of signal, we present upper limits for different energy thresholds above ${\gtrsim}10^8$\,GeV on the secondary by-product fluxes expected from the decay of the particles. Assuming that the energy density of these super-heavy particles matches that of dark matter observed today, we translate the upper bounds on the particle fluxes into tight constraints on the couplings governing the decay process as a function of the particle mass. Instantons, which are non-perturbative solutions to Yang-Mills equations, can give rise to decay channels otherwise forbidden and transform stable particles into meta-stable ones. Assuming such instanton-induced decay processes, we derive a bound on the reduced coupling constant of gauge interactions in the dark sector: $\alpha_X \lesssim 0.09$, for $10^{9} \lesssim M_X/\text{GeV} < 10^{19}$. Conversely, we obtain that, for instance, a reduced coupling constant $\alpha_X = 0.09$ excludes masses $M_X \gtrsim 3\times 10^{13}~$GeV. In the context of dark matter production from gravitational interactions alone during the reheating epoch, we derive constraints on the parameter space that involves, in addition to $M_X$ and $\alpha_X$, the Hubble rate at the end of inflation, the reheating efficiency, and the non-minimal coupling of the Higgs with curvature.

We use the surface detector of the Pierre Auger Observatory to search for air showers initiated by photons with an energy above 1019 eV. Photons in the zenith angle range from 30∘ to 60∘ can be identified in the overwhelming background of showers initiated by charged cosmic rays through the broader time structure of the signals induced in the water-Cherenkov detectors of the array and the steeper lateral distribution of shower particles reaching ground. Applying the search method to data collected between January 2004 and June 2020, upper limits at 95% CL are set to an E-2 diffuse flux of ultra-high energy photons above 1019 eV, 2 × 1019 eV and 4 × 1019 eV amounting to 2.11 × 10-3, 3.12 × 10-4 and 1.72 × 10-4 km-2 sr-1 yr-1, respectively. While the sensitivity of the present search around 2 × 1019 eV approaches expectations of cosmogenic photon fluxes in the case of a pure-proton composition, it is one order of magnitude above those from more realistic mixed-composition models. The inferred limits have also implications for the search of super-heavy dark matter that are discussed and illustrated.

A search for time-directional coincidences of ultra-high-energy (UHE) photons above 10 EeV with gravitational wave (GW) events from the LIGO/Virgo runs O1 to O3 is conducted with the Pierre Auger Observatory. Due to the distinctive properties of photon interactions and to the background expected from hadronic showers, a subset of the most interesting GW events is selected based on their localization quality and distance. Time periods of 1000 s around and 1 day after the GW events are analyzed. No coincidences are observed. Upper limits on the UHE photon fluence from a GW event are derived that are typically at & SIM;7 MeV cm(-2) (time period 1000 s) and & SIM;35 MeV cm(-2) (time period 1 day). Due to the proximity of the binary neutron star merger GW170817, the energy of the source transferred into UHE photons above 40 EeV is constrained to be less than 20% of its total GW energy. These are the first limits on UHE photons from GW sources.

Ultra-high-energy photons with energies exceeding 10(17) 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 10(15) eV range further motivate searches for even higher-energy photons. In this paper, we present a search for photons with energies exceeding 2 x 10(17) 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 10(17) and 10(18) eV.

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 ≡ 10 18 eV). Despite the flux of these particles being extremely low, the area of ∼3000 km 2 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 km 2 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.

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 10^19 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 10^20 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.

We use the surface detector of the Pierre Auger Observatory to search for air showers initiated by photons with an energy above 10^19 eV. Photons in the zenith angle range from 30^∘ to 60^∘ can be identified in the overwhelming background of showers initiated by charged cosmic rays through the broader time structure of the signals induced in the water-Cherenkov detectors of the array and the steeper lateral distribution of shower particles reaching ground. Applying the search method to data collected between January 2004 and June 2020, upper limits at 95% CL are set to an E^-2 diffuse flux of ultra-high energy photons above 10^19 eV, 2×10^19 eV and 4×10^19 eV amounting to 2.11×10^-3, 3.12×10^-4 and 1.72×10^-4 km^-2 sr^-1 yr^-1, respectively. While the sensitivity of the present search around 2 × 10^19 eV approaches expectations of cosmogenic photon fluxes in the case of a pure-proton composition, it is one order of magnitude above those from more realistic mixed-composition models. The inferred limits have also implications for the search of super-heavy dark matter that are discussed and illustrated.

We investigate instanton-induced decay processes of super-heavy dark matter particles $X$ produced during the inflationary epoch. Using data collected at the Pierre Auger Observatory we derive a bound on the reduced coupling constant of gauge interactions in the dark sector: $\alpha_X^{\rm eff} \lesssim 0.09$, for $10^{10}<M_X/{\rm GeV}<10^{16}$. We show that this upper limit on $\alpha_X^{\rm eff}$ is complementary to that obtained from the non-observation of tensor modes in the cosmic microwave background.

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 is provided by the IceCube Neutrino Observatory and ANTARES, while the UHECR data with energies above $\sim$50 EeV is 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 neutrinos 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 has found a significant excess, and previously reported over-fluctuations are reduced in significance. Based on these results, we further constrain the neutrino flux spatially correlated with UHECRs.

The development of an extensive air shower depends not only on the nature of the primary ultra-high-energy cosmic ray but also on the properties of the hadronic interactions. For energies above those achievable in human-made accelerators, hadronic interactions are only accessible through the studies of extensive air showers, which can be measured at the Pierre Auger Observatory. With its hybrid detector design, the Pierre Auger Observatory measures both the longitudinal development of showers in the atmosphere and the lateral distribution of particles that arrive at the ground. This way, observables that are sensitive to hadronic interactions at ultra-high energies can be obtained. While the hadronic interaction cross-section can be assessed from the longitudinal profiles, the number of muons and their fluctuations measured with the ground detectors are linked to other physical properties. In addition to these direct studies, we discuss here how measurements of the atmospheric depth of the maximum of air-shower profiles and the characteristics of the muon signal at the ground can be used to test the self-consistency of the post-LHC hadronic models.

AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to complement the study of ultra-high-energy cosmic rays (UHECR) by measuring the muon content of extensive air showers (EAS). It consists of an array of 61 water Cherenkov detectors on a denser spacing in combination with underground scintillation detectors used for muon density measurement. Each detector is composed of three scintillation modules, with 10 m^2 detection area per module, buried at 2.3 m depth, resulting in a total detection area of 30 m^2. Silicon photomultiplier sensors (SiPM) measure the amount of scintillation light generated by charged particles traversing the modules. In this paper, the design of the front-end electronics to process the signals of those SiPMs and test results from the laboratory and from the Pierre Auger Observatory are described. Compared to our previous prototype, the new electronics shows a higher performance, higher efficiency and lower power consumption, and it has a new acquisition system with increased dynamic range that allows measurements closer to the shower core. The new acquisition system is based on the measurement of the total charge signal that the muonic component of the cosmic ray shower generates in the detector.

The Auger Muon Infill Ground Array (AMIGA) is part of the AugerPrime upgrade of the Pierre Auger Observatory. It consists of particle counters buried 2.3 m underground next to the water-Cherenkov stations that form the 23.5 km^2 large infilled array. The reduced distance between detectors in this denser area allows the lowering of the energy threshold for primary cosmic ray reconstruction down to about 10^17 eV. At the depth of 2.3 m the electromagnetic component of cosmic ray showers is almost entirely absorbed so that the buried scintillators provide an independent and direct measurement of the air showers muon content. This work describes the design and implementation of the AMIGA embedded system, which provides centralized control, data acquisition and environment monitoring to its detectors. The presented system was firstly tested in the engineering array phase ended in 2017, and lately selected as the final design to be installed in all new detectors of the production phase. The system was proven to be robust and reliable and has worked in a stable manner since its first deployment.

We report the results of the analyses of the cosmic ray data collected with a 4 tonne (3×1×1 m3) active mass (volume) Liquid Argon Time-Projection Chamber (TPC) operated in a dual-phase mode. We present a detailed study of the TPC's response, its main detector parameters and performance. The results are important for the understanding and further developments of the dual-phase technology, thanks to the verification of key aspects, such as the extraction of electrons from liquid to gas and their amplification through the entire one square metre readout plain, gain stability, purity and charge sharing between readout views.

The atmospheric depth of the air shower maximum $X_{\mathrm{max}}$ is an observable commonly used for the determination of the nuclear mass composition of ultra-high energy cosmic rays. Direct measurements of $X_{\mathrm{max}}$ are performed using observations of the longitudinal shower development with fluorescence telescopes. At the same time, several methods have been proposed for an indirect estimation of $X_{\mathrm{max}}$ from the characteristics of the shower particles registered with surface detector arrays. In this paper, we present a deep neural network (DNN) for the estimation of $X_{\mathrm{max}}$. The reconstruction relies on the signals induced by shower particles in the ground based water-Cherenkov detectors of the Pierre Auger Observatory. The network architecture features recurrent long short-term memory layers to process the temporal structure of signals and hexagonal convolutions to exploit the symmetry of the surface detector array. We evaluate the performance of the network using air showers simulated with three different hadronic interaction models. Thereafter, we account for long-term detector effects and calibrate the reconstructed $X_{\mathrm{max}}$ using fluorescence measurements. Finally, we show that the event-by-event resolution in the reconstruction of the shower maximum improves with increasing shower energy and reaches less than $25~\mathrm{g/cm^{2}}$ at energies above $2\times 10^{19}~\mathrm{eV}$.