Search for: All records

Abstract The detection of gravitational waves (GWs) from a binary neutron star (BNS) merger by Advanced LIGO and Advanced Virgo (GW170817), together with its electromagnetic counterpart, the short gamma-ray burst GRB 170817A, heralded the birth of multimessenger astronomy. The detection of TeV emission from GRBs motivates follow-up observations with the Cherenkov Telescope Array Observatory (CTAO), which is ideal for detecting such signals due to its unprecedented sensitivity, rapid response, and wide-field survey capabilities. The aim of this work is to evaluate GeV–TeV GW follow-up strategies for CTAO using a multistep simulation pipeline and to estimate the expected rate of joint GW–GRB detections during observing run O5. Using a simulated sample of BNS systems with corresponding GW detections, gamma-ray emission is simulated through phenomenological prescriptions based on the observed population of short GRBs, including off-axis jet scenarios. CTAO observations are simulated to account for instrument response, sky tiling strategies, integration times, and varying observing conditions. Strategies with variable and constant integration times are investigated. We find that, via an optimized follow-up strategy, about 5% of simulated GW-associated short GRBs produce GeV–TeV radiation detectable by CTAO. Detectability is strongly influenced by the jet opening angle and viewing angle, suggesting that even rough estimates of the viewing angle in GW alerts could enhance targeting. This framework motivates future follow-ups of GW-detectable events, including neutron star–black hole mergers, and further supports the development of advanced strategies incorporating galaxy distributions and synergies with future detectors such as the Einstein Telescope. 

We present a novel approach for assessing the muon content of air showers with large zenith angles on a combined analysis of their radio emission and particle footprint. We use the radiation energy reconstructed by the Auger engineering radio array (AERA) as an energy estimator and determine the muon number independently with the water-Cherenkov detector array of the Pierre Auger Observatory, deployed on a 1500 m grid. We focus our analysis on air showers with primary energy above 4 EeV to ensure full detection efficiency. Over approximately ten years of accumulated data, we identify a set of 40 high-quality events that are used in the analysis. The estimated muon contents in data are compatible with those for iron primaries as predicted by current-generation hadronic interaction models. This result can be interpreted as a deficit of muons in simulations as a lighter mass composition has been established from $X_{\max }$measurements. This muon deficit was already observed in previous analyses of the Auger Collaboration and is confirmed using hybrid events that include radio measurements for the first time. 

The energy spectrum of cosmic rays above 2.5 EeV has been measured across the declination range −90° ≤ δ ≤ þ44.8° using ∼310 000 events accrued at the Pierre Auger Observatory from an exposure of ð104 900 3 100Þ km2 sr yr. No significant variations of energy spectra with declination are observed, after allowing or not for nonuniformities across the sky arising from the well-established dipolar anisotropies in the arrival directions of ultrahigh-energy cosmic rays. The instep feature in the spectrum at ≃10 EeV reported previously is now established at a significance above 5σ. Within the statistics, the energy spectra are indistinguishable across declinations so disfavoring an origin for the instep from a few distinctive sources. 

Abstract Deflections in the propagation of charged ultra-high-energy cosmic rays (UHECRs) caused by magnetic fields make the identification of their sources challenging. On the other hand, the arrival directions at Earth of neutrons point directly to their origin. The emission of UHECRs from a source is expected to be accompanied by the production of neutrons in its vicinity through interactions with ambient matter and radiation. Since free neutrons travel a mean distanced kpc⁻¹ = 9.2 (E EeV⁻¹) before decaying, a neutron flux in the EeV range could be detected on Earth from sources of UHECRs in our Galaxy. Using cosmic-ray data from Phase I of the Surface Detector of the Pierre Auger Observatory, we search for neutron fluxes from Galactic candidate sources. We select more than 1000 objects of astrophysical interest, stacking them into target sets. The targets all have decl. within the exposure of the Observatory, ranging from −90° up to +45° for energies above 1 EeV (and up to +20° for energies down to 0.1 EeV). Given that a neutron air shower is indistinguishable from a proton one, there is a significant background due to cosmic rays. A neutron flux from the direction of a candidate source would be identified by a celestial density of events that significantly exceeds the expected density of cosmic rays for that direction. No significant excess is found at any tested target direction, and an upper limit on the neutron flux is calculated for each candidate source. 

ABSTRACT The dwarf spheroidal galaxies (dSphs) orbiting the Milky Way are widely regarded as systems supported by velocity dispersion against self-gravity, and as prime targets for the search for indirect dark matter (DM) signatures in the GeV-to-TeV $$\gamma$$-ray range owing to their lack of astrophysical $$\gamma$$-ray background. We present forecasts of the sensitivity of the forthcoming Cherenkov Telescope Array Observatory (CTAO) to annihilating or decaying DM signals in these targets. An original selection of candidates is performed from the current catalogue of known objects, including both classical and ultrafaint dSphs. For each, the expected DM content is derived using the most comprehensive photometric and spectroscopic data available, within a consistent framework of analysis. This approach enables the derivation of novel astrophysical factor profiles for indirect DM searches, which are compared with results from the literature. From an initial sample of 64 dSphs, eight promising targets are identified – Draco I, Coma Berenices, Ursa Major II, Ursa Minor, and Willman 1 in the North, Reticulum II, Sculptor, and Sagittarius II in the South – for which different DM density models yield consistent expectations, leading to robust predictions. CTAO is expected to provide the strongest limits above $$\sim$$10 TeV, reaching velocity-averaged annihilation cross sections of $$\sim 5\times 10^{-25}$$ cm$^3$ s$$^{-1}$$ and decay lifetimes up to $$\sim 10^{26}$$ s for combined limits. The dominant uncertainties arise from the imprecise determination of the DM content, particularly for ultrafaint dSphs. Observation strategies are proposed that optimize either deep exposures of the best candidates or diversified target selections. 

Abstract The modulation of low-energy galactic cosmic rays reflects interplanetary magnetic field variations and can provide useful information on solar activity. An array of ground-surface detectors can reveal the secondary particles, which originate from the interaction of cosmic rays with the atmosphere. In this work, we present an investigation of the low-threshold rate (scaler) time series recorded in 16 yr of operation by the Pierre Auger Observatory surface detectors in Malargüe, Argentina. Through an advanced spectral analysis, we detected highly statistically significant variations in the time series with periods ranging from the decadal to the daily scale. We investigate their origin, revealing a direct connection with solar variability. Thanks to their intrinsic very low noise level, the Auger scalers allow a thorough and detailed investigation of the galactic cosmic-ray flux variations in the heliosphere at different timescales and can, therefore, be considered a new proxy of solar variability. 

A dedicated search for upward-going air showers at zenith angles exceeding 110° and energies $E>0.1\mathrm{EeV}$has been performed using the Fluorescence Detector of the Pierre Auger Observatory. The search is motivated by two “anomalous” radio pulses observed by the ANITA flights I and III that appear inconsistent with the standard model of particle physics. Using simulations of both regular cosmic-ray showers and upward-going events, a selection procedure has been defined to separate potential upward-going candidate events and the corresponding exposure has been calculated in the energy range [0.1–33] EeV. One event has been found in the search period between January 1, 2004, and December 31, 2018, consistent with an expected background of $0.27\pm 0.12$events from misreconstructed cosmic-ray showers. This translates to an upper bound on the integral flux of $(7.2\pm 0.2)\times {10}^{-21}{\mathrm{cm}}^{-2}{\mathrm{sr}}^{-1}{\mathrm{y}}^{-1}$and $(3.6\pm 0.2)\times {10}^{-20}{\mathrm{cm}}^{-2}{\mathrm{sr}}^{-1}{\mathrm{y}}^{-1}$for an $E^{-1}$and $E^{-2}$spectrum, respectively. An upward-going flux of showers normalized to the ANITA observations is shown to predict over 34 events for an $E^{-3}$spectrum and over 8.1 events for a conservative $E^{-5}$spectrum, in strong disagreement with the interpretation of the anomalous events as upward-going showers. Published by the American Physical Society2025 

We present measurements of the atmospheric depth of the shower maximum $X_{\max }$, inferred for the first time on an event-by-event level using the surface detector of the Pierre Auger Observatory. Using deep learning, we were able to extend measurements of the $X_{\max }$distributions up to energies of 100 EeV ( ${10}^{20}\mathrm{eV}$), not yet revealed by current measurements, providing new insights into the mass composition of cosmic rays at extreme energies. Gaining a 10-fold increase in statistics compared to the fluorescence detector data, we find evidence that the rate of change of the average $X_{\max }$with the logarithm of energy features three breaks at $6.5\pm 0.6\left(\mathrm{stat}\right)\pm 1\left(\mathrm{syst}\right)\mathrm{EeV}$, $11\pm 2\left(\mathrm{stat}\right)\pm 1\left(\mathrm{syst}\right)\mathrm{EeV}$, and $31\pm 5\left(\mathrm{stat}\right)\pm 3\left(\mathrm{syst}\right)\mathrm{EeV}$, in the vicinity to the three prominent features (ankle, instep, suppression) of the cosmic-ray flux. The energy evolution of the mean and standard deviation of the measured $X_{\max }$distributions indicates that the mass composition becomes increasingly heavier and purer, thus being incompatible with a large fraction of light nuclei between 50 and 100 EeV. Published by the American Physical Society2025 

We report an investigation of the mass composition of cosmic rays with energies from 3 to 100 EeV ( $1\mathrm{EeV}={10}^{18}\mathrm{eV}$) using the distributions of the depth of shower maximum $X_{\max }$. The analysis relies on $\sim 50,000$events recorded by the surface detector of the Pierre Auger Observatory and a deep-learning-based reconstruction algorithm. Above energies of 5 EeV, the dataset offers a 10-fold increase in statistics with respect to fluorescence measurements at the Observatory. After cross-calibration using the fluorescence detector, this enables the first measurement of the evolution of the mean and the standard deviation of the $X_{\max }$distributions up to 100 EeV. Our findings are threefold: (i) The evolution of the mean logarithmic mass toward a heavier composition with increasing energy can be confirmed and is extended to 100 EeV. (ii) The evolution of the fluctuations of $X_{\max }$toward a heavier and purer composition with increasing energy can be confirmed with high statistics. We report a rather heavy composition and small fluctuations in $X_{\max }$at the highest energies. (iii) We find indications for a characteristic structure beyond a constant change in the mean logarithmic mass, featuring three breaks that are observed in proximity to the ankle, instep, and suppression features in the energy spectrum. Published by the American Physical Society2025 

Abstract The Pierre Auger Collaboration has embraced the concept of open access to their research data since its foundation, with the aim of giving access to the widest possible community. A gradual process of release began as early as 2007 when 1% of the cosmic-ray data was made public, along with 100% of the space-weather information. In February 2021, a portal was released containing 10% of cosmic-ray data collected by the Pierre Auger Observatory from 2004 to 2018, during the first phase of operation of the Observatory. The Open Data Portal includes detailed documentation about the detection and reconstruction procedures, analysis codes that can be easily used and modified and, additionally, visualization tools. Since then, the Portal has been updated and extended. In 2023, a catalog of the highest-energy cosmic-ray events examined in depth has been included. A specific section dedicated to educational use has been developed with the expectation that these data will be explored by a wide and diverse community, including professional and citizen scientists, and used for educational and outreach initiatives. This paper describes the context, the spirit, and the technical implementation of the release of data by the largest cosmic-ray detector ever built and anticipates its future developments.