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.

Deep Learning methods are among the state-of-art of several computer vision tasks, intelligent control systems, fast and reliable signal processing and inference in big data regimes.It is also a promising tool for scientific analysis, such as gamma/hadron discrimination.We present an approach based on Deep Learning for the regression of shower parameters, namely the its core position and ground energy, using water Cherenkov detectors.We design our method using simulations.We evaluate the limits of such estimation near the borders of the arrays, including when the center is outside the detector's range.We used Bayesian Neural Networks and derived and quantified systematic errors arising from Deep Learning models and in an EfficientNet model design.The method could be easily adapted to estimate other parameters.

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.

The Southern Wide-field Gamma-ray Observatory (SWGO) is the next-generation gamma-ray observatory, currently in an R&D phase. The experiment is expected to have a large array of water Cherenkov detectors (WCD) placed at a high elevation (> 4.4 km a.s.l.) in South America. Here we present a WCD concept with reduced surface area and height of stations comprising four PMTs at the bottom. We show that it is possible to reach an excellent gamma/hadron discrimination by analysing the data gathered by this station with machine learning techniques. Such performance can be achieved by analysing the shower patterns at the ground or through the PMTs signal time structure to tag muons. Moreover, it is shown that the station's performance does not depend on the array configuration (dense or sparse) nor on the shower inclination (theta < 40 degrees). Such a concept reduces the cost associated with the transport of massive amounts of water to high elevation sites while keeping a high physics performance. Therefore, it could be a good candidate station for SWGO, enabling to reach good sensitivities from low energies (similar to 100 GeV) up to the PeV region, covering large ground surface areas (few square km).

The concept of a small, single-layer water Cherenkov detector, with three photomultiplier tubes (PMTs), placed at its bottom in a 120^∘ star configuration ( Mercedes Water Cherenkov Detector) is presented. The PMTs are placed near the lateral walls of the stations with an adjustable inclination and may be installed inside or outside the water volume. To illustrate the technical viability of this concept and obtain a first-order estimation of its cost, an engineering design was elaborated. The sensitivity of these stations to low energy Extensive Air Shower (EAS) electrons, photons and muons is discussed, both in compact and sparse array configurations. It is shown that the analysis of the intensity and time patterns of the PMT signals, using machine learning techniques, enables the tagging of muons, achieving an excellent gamma/hadron discrimination for TeV showers. This concept minimises the station production and maintenance costs, allowing for a highly flexible and fast installation. Mercedes Water Cherenkov Detectors (WCDs) are thus well-suited for use in high-altitude large gamma-ray observatories covering an extended energy range from the low energies, closing the gap between satellite and ground-based measurements, to very high energy regions, beyond the PeV scale.

Authors: C. Boisson (LUTH), A.M. Brown (Durham U.), A. Burtovoi (CISAS, UNIPD, INAF), M. Cerruti (APC, CNRS/IN2P3), M. Chernyakova (Dublin City U.), T. Hassan (CIEMAT), J.-P. Lenain (LPNHE, CNRS/IN2P3), M. Manganaro (Univ. Rijeka), P. Romano (INAF), H. Sol (LUTH), F. Tavecchio (INAF), S. Vercellone (INAF), L. Zampieri (INAF), R. Zanin (CTAO), A. Zech (LUTH), (continues on the next page). Contact points: jlenain@in2p3.fr, patrizia.romano@inaf.it

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.

The measurement of the mass composition of ultra-high energy cosmic rays constitutes a prime challenge in astroparticle physics. Most detailed information on the composition can be obtained from measurements of the depth of maximum of air showers, $X_{\mathrm{max}}$, with the use of fluorescence telescopes, which can be operated only during clear and moonless nights. Using deep neural networks, it is now possible for the first time to perform an event-by-event reconstruction of $X_{\mathrm{max}}$ with the Surface Detector (SD) of the Pierre Auger Observatory. Therefore, previously recorded data can be analyzed for information on $X_{\mathrm{max}}$, and thus, the cosmic-ray composition. Since the SD operates with a duty cycle of almost $100\%$ and its event selection is less strict than for the Fluorescence Detector (FD), the gain in statistics with respect to the FD is almost a factor of 15 for energies above $10^{19.5}~\mathrm{eV}$. In this contribution, we introduce the neural network particularly designed for the SD of the Pierre Auger Observatory. We evaluate its performance using three different hadronic interaction models, verify its functionality using Auger hybrid measurements, and find that the method can extract mass information on an event level.

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}$.

The physics of the non-thermal Universe provides information on the acceleration mechanisms in extreme environments, such as black holes and relativistic jets, neutron stars, supernovae or clusters of galaxies. In the presence of magnetic fields, particles can be accelerated towards relativistic energies. As a consequence, radiation along the entire electromagnetic spectrum can be observed, and extreme environments are also the most likely sources of multi-messenger emission. The most energetic part of the electromagnetic spectrum corresponds to the very-high-energy (VHE, E>100 GeV) gamma-ray regime, which can be extensively studied with ground based Imaging Atmospheric Cherenkov Telescopes (IACTs). The results obtained by the current generation of IACTs, such as H.E.S.S., MAGIC, and VERITAS, demonstrate the crucial importance of the VHE band in understanding the non-thermal emission of extreme environments in our Universe. In some objects, the energy output in gamma rays can even outshine the rest of the broadband spectrum. The Cherenkov Telescope Array (CTA) is the next generation of IACTs, which, with cutting edge technology and a strategic configuration of ~100 telescopes distributed in two observing sites, in the northern and southern hemispheres, will reach better sensitivity, angular and energy resolution, and broader energy coverage than currently operational IACTs. With CTA we can probe the most extreme environments and considerably boost our knowledge of the non-thermal Universe.

The first Large-Sized Telescope (LST-1) proposed for the forthcoming Cherenkov Telescope Array (CTA) has started to operate in 2019 in La Palma. The large structure of LST-1 - with a 23 m mirror dish diameter - imposes a strict control of its deformations that could affect the pointing accuracy and its overall performance. According to CTA specifications that are conceived to resolve e.g. the fine structure of galactic sources, the LST post-calibration pointing accuracy should be better than 14 arcseconds. To fulfill this requirement, the telescope pointing precision is monitored with two dedicated CCD cameras located at the dish center. The analysis of their images allows us to disentangle different systematic deformations of the structure. In this work, we investigate a complementary approach that offers the possibility to monitor the pointing of the telescope during the acquisition of sky data. After properly cleaning the events from the Cherenkov showers, the reconstructed positions of the stars imaged in the camera field of view are compared to their nominal expected positions in catalogues. This provides a direct measurement of the telescope pointing, that can be used to cross-check the other methods and as a real-time monitoring of the optical properties of the telescope and of the pointing corrections applied by the bending models. Additionally, this method benefits from not relying on specific hardware or dedicated observations. In this contribution we will illustrate this analysis and show results based on simulations of LST-1.

We present a measurement of the cosmic-ray spectrum above 100 PeV using the part of the surface detector of the Pierre Auger Observatory that has a spacing of 750 m. An inflection of the spectrum is observed, confirming the presence of the so-called second-knee feature. The spectrum is then combined with that of the 1500 m array to produce a single measurement of the flux, linking this spectral feature with the three additional breaks at the highest energies. The combined spectrum, with an energy scale set calorimetrically via fluorescence telescopes and using a single detector type, results in the most statistically and systematically precise measurement of spectral breaks yet obtained. These measurements are critical for furthering our understanding of the highest energy cosmic rays.

Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

Elves are a class of transient luminous events, with a radial extent typically greater than 250km, that occur in the lower ionosphere above strong electrical storms. We report the observation of 1,598 elves, from 2014 to 2016, recorded with unprecedented time resolution (100 ns) using the fluorescence detector (FD) of the Pierre Auger Cosmic-Ray Observatory. The Auger Observatory is located in the Mendoza province of Argentina with a viewing footprint for elve observations of 3106km 2, reaching areas above the Pacific and Atlantic Oceans, as well as the Cordoba region, which is known for severe convective thunderstorms. Primarily designed for ultrahigh energy cosmic-ray observations, the Auger FD turns out to be very sensitive to the ultraviolet emission in elves. The detector features modified Schmidt optics with large apertures resulting in a field of view that spans the horizon, and year-round operation on dark nights with low moonlight background, when the local weather is favorable. The measured light profiles of 18% of the elve events have more than one peak, compatible with intracloud activity. Within the 3-year sample, 72% of the elves correlate with the far-field radiation measurements of the World Wide Lightning Location Network. The Auger Observatory plans to continue operations until at least 2025, including elve observations and analysis. To the best of our knowledge, this observatory is the only facility on Earth that measures elves with year-round operation and full horizon coverage.

The hybrid design of the Pierre Auger Observatory allows for the measurement of the properties of extensive air showers initiated by ultra-high energy cosmic rays with unprecedented precision. By using an array of prototype underground muon detectors, we have performed the first direct measurement, by the Auger Collaboration, of the muon content of air showers between 2× 10^17 and 2× 10^18 eV. We have studied the energy evolution of the attenuation-corrected muon density, and compared it to predictions from air shower simulations. The observed densities are found to be larger than those predicted by models. We quantify this discrepancy by combining the measurements from the muon detector with those from the Auger fluorescence detector at 10^17.5 eV and 10^18 eV . We find that, for the models to explain the data, an increase in the muon density of 38% ± 4% (12%) ±^21%_18% for EPOS-LHC , and of 50% (53%) ± 4% (13%) ±^23%_20% for QGSJetII-04 , is respectively needed.