Search for: All records

Extragalactic background light (EBL) plays an important role in cosmology since it traces the history of galaxy formation and evolution. Such diffuse radiation from near-UV to far-infrared wavelengths can interact with γ -rays from distant sources such as active galactic nuclei (AGNs), and is responsible for the high-energy absorption observed in their spectra. However, probing the EBL from γ -ray spectra of AGNs is not trivial due to internal processes that can mimic its effect. Such processes are usually taken into account in terms of curvature of the intrinsic spectrum. Hence, an improper choice of parametrization for the latter can seriously affect EBL reconstruction. In this paper, we propose a statistical approach that avoids a priori assumptions on the intrinsic spectral curvature and that, for each source, selects the best-fit model on a solid statistical basis. By combining the Fermi -LAT observations of 490 blazars, we determine the γ -ray-inferred level of EBL for various state-of-the-art EBL models. We discuss the EBL level obtained from the spectra of both BL Lacs and flat spectrum radio quasars (FSRQ) in order to investigate the impact of internal absorption in different classes of objects. We further scrutinize constraints on the EBL evolution from γ -ray observations by reconstructing the EBL level in four redshift ranges, up to z  ∼ 2.5. The approach implemented in this paper, carefully addressing the question of the modeling of the intrinsic emission at the source, can serve as a solid stepping stone for studies of hundreds of high-quality spectra acquired by next-generation γ -ray instruments. 

The IceCube Neutrino Observatory at the geographic South Pole, with its surface array IceTop, detects three different components of extensive air showers: the total signal at the surface, low energy muons on the periphery of the showers, and high energy muons in the deep In Ice array of IceCube. These measurements enable determination of the energy spectrum and composition of cosmic rays from PeV to EeV energies, the anisotropy in the distribution of cosmic ray arrival directions, the muon density of cosmic ray air showers, and the PeV gamma-ray flux. Furthermore, IceTop can be used as a veto for the neutrino measurements. The latest results from these IceTop analyses will be presented along with future plans. 

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. 

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 Results are presented for the measurement of large-scale anisotropies in the arrival directions of ultra–high-energy cosmic rays detected at the Pierre Auger Observatory during 19 yr of operation, prior to AugerPrime, the upgrade of the observatory. The 3D dipole amplitude and direction are reconstructed above 4 EeV in four energy bins. Besides the established dipolar anisotropy in R.A. above 8 EeV, the Fourier amplitude of the 8–16 EeV energy bin is now also above the 5σdiscovery level. No time variation of the dipole moment above 8 EeV is found, setting an upper limit to the rate of change of such variations of 0.3% yr⁻¹at the 95% confidence level. Additionally, the results for the angular power spectrum are shown, demonstrating no other statistically significant multipoles. The results for the equatorial dipole component down to 0.03 EeV are presented, using for the first time a data set obtained with a trigger that has been optimized for lower energies. Finally, model predictions are discussed and compared with observations, based on two source emission scenarios obtained in the combined fit of spectrum and composition above 0.6 EeV. 

The Pierre Auger Observatory is the most sensitive instrument to detect photons with energies above $10^{17}\mathrm{eV}$. It measures extensive air showers generated by ultrahigh energy cosmic rays using a hybrid technique that exploits the combination of a fluorescence detector with a ground array of particle detectors. The signatures of a photon-induced air shower are a larger atmospheric depth of the shower maximum ( $X_{\max }$) and a steeper lateral distribution function, along with a lower number of muons with respect to the bulk of hadron-induced cascades. In this work, a new analysis technique in the energy interval between 1 and 30 EeV ( $1\mathrm{EeV}=10^{18}\mathrm{eV}$) has been developed by combining the fluorescence detector-based measurement of $X_{\max }$with the specific features of the surface detector signal through a parameter related to the air shower muon content, derived from the universality of the air shower development. No evidence of a statistically significant signal due to photon primaries was found using data collected in about 12 years of operation. Thus, upper bounds to the integral photon flux have been set using a detailed calculation of the detector exposure, in combination with a data-driven background estimation. The derived 95% confidence level upper limits are 0.0403, 0.01113, 0.0035, 0.0023, and $0.0021{\mathrm{km}}^{-2}{\mathrm{sr}}^{-1}{\mathrm{yr}}^{-1}$above 1, 2, 3, 5, and 10 EeV, respectively, leading to the most stringent upper limits on the photon flux in the EeV range. Compared with past results, the upper limits were improved by about 40% for the lowest energy threshold and by a factor 3 above 3 EeV, where no candidates were found and the expected background is negligible. The presented limits can be used to probe the assumptions on chemical composition of ultrahigh energy cosmic rays and allow for the constraint of the mass and lifetime phase space of super-heavy dark matter particles. Published by the American Physical Society2024 

Abstract Galaxy clusters are expected to be both dark matter (DM) reservoirs and storage rooms for the cosmic-ray protons (CRp) that accumulate along the cluster's formation history. Accordingly, they are excellent targets to search for signals of DM annihilation and decay atγ-ray energies and are predicted to be sources of large-scaleγ-ray emission due to hadronic interactions in the intracluster medium (ICM).In this paper, we estimate the sensitivity of the Cherenkov Telescope Array (CTA) to detect diffuseγ-ray emission from the Perseus galaxy cluster.We first perform a detailed spatial and spectral modelling of the expected signal for both the DM and the CRp components. For each case, we compute the expected CTA sensitivity accounting for the CTA instrument response functions. The CTA observing strategy of the Perseus cluster is also discussed.In the absence of a diffuse signal (non-detection), CTA should constrain the CRp to thermal energy ratioX₅₀₀within the characteristic radiusR₅₀₀down to aboutX₅₀₀< 3 × 10^-3, for a spatial CRp distribution that follows the thermal gas and a CRp spectral index α_CRp= 2.3. Under the optimistic assumption of a pure hadronic origin of the Perseus radio mini-halo and depending on the assumed magnetic field profile, CTA should measure α_CRpdown to about Δα_CRp≃ 0.1 and the CRp spatial distribution with 10% precision, respectively. Regarding DM, CTA should improve the current ground-basedγ-ray DM limits from clusters observations on the velocity-averaged annihilation cross-section by a factor of up to ∼ 5, depending on the modelling of DM halo substructure. In the case of decay of DM particles, CTA will explore a new region of the parameter space, reaching models withτ_χ> 10²⁷s for DM masses above 1 TeV.These constraints will provide unprecedented sensitivity to the physics of both CRp acceleration and transport at cluster scale and to TeV DM particle models, especially in the decay scenario. 

Abstract The flux of ultra-high energy cosmic rays reaching Earth above the ankle energy (5 EeV) can be described as a mixture of nuclei injected by extragalactic sources with very hard spectra and a low rigidity cutoff.Extragalactic magnetic fields existing between the Earth and the closest sources can affect the observed CR spectrum by reducing the flux of low-rigidity particles reaching Earth. We perform a combined fit of the spectrum and distributions of depth of shower maximum measured with the Pierre Auger Observatory including the effect of this magnetic horizon in the propagation of UHECRs in the intergalactic space.We find that, within a specific range of the various experimental and phenomenological systematics, the magnetic horizon effect can be relevant for turbulent magnetic field strengths in the local neighbourhood in which the closest sources lieof order B_rms≃ (50–100) nG (20 Mpc/d_s)( 100 kpc/L_coh)^1/2, with d_sthe typical intersource separation and L_cohthe magnetic field coherence length. When this is the case,the inferred slope of the source spectrum becomes softer and can be closer to the expectations of diffusive shock acceleration, i.e., ∝ E^-2.An additional cosmic-ray population with higher source density and softer spectra, presumably also extragalactic and dominating the cosmic-ray flux at EeV energies, is also required to reproduce the overall spectrum and composition results for all energies down to 0.6 EeV. 

Abstract Approximately one hundred sources of very-high-energy (VHE) gamma rays are known in the Milky Way, detected with a combination of targeted observations and surveys. A survey of the entire Galactic Plane in the energy range from a few tens of GeV to a few hundred TeV has been proposed as a Key Science Project for the upcoming Cherenkov Telescope Array Observatory (CTAO). This article presents the status of the studies towards the Galactic Plane Survey (GPS). We build and make publicly available a sky model that combines data from recent observations of known gamma-ray emitters with state-of-the-art physically-driven models of synthetic populations of the three main classes of established Galactic VHE sources (pulsar wind nebulae, young and interacting supernova remnants, and compact binary systems), as well as of interstellar emission from cosmic-ray interactions in the Milky Way. We also perform an optimisation of the observation strategy (pointing pattern and scheduling) based on recent estimations of the instrument performance. We use the improved sky model and observation strategy to simulate GPS data corresponding to a total observation time of 1620 hours spread over ten years. Data are then analysed using the methods and software tools under development for real data. Under our model assumptions and for the realisation considered, we show that the GPS has the potential to increase the number of known Galactic VHE emitters by almost a factor of five. This corresponds to the detection of more than two hundred pulsar wind nebulae and a few tens of supernova remnants at average integral fluxes one order of magnitude lower than in the existing sample above 1 TeV, therefore opening the possibility to perform unprecedented population studies. The GPS also has the potential to provide new VHE detections of binary systems and pulsars, to confirm the existence of a hypothetical population of gamma-ray pulsars with an additional TeV emission component, and to detect bright sources capable of accelerating particles to PeV energies (PeVatrons). Furthermore, the GPS will constitute a pathfinder for deeper follow-up observations of these source classes. Finally, we show that we can extract from GPS data an estimate of the contribution to diffuse emission from unresolved sources, and that there are good prospects of detecting interstellar emission and statistically distinguishing different scenarios.Thus, a survey of the entire Galactic plane carried out from both hemispheres with CTAO will ensure a transformational advance in our knowledge of Galactic VHE source populations and interstellar emission.