The Near-Infrared Spectrometer and Photometer (NISP) on board the satellite provides multiband photometry and $R slitless grism spectroscopy in the 950--2020\,nm wavelength range. In this reference article, we illuminate the background of NISP's functional and calibration requirements, describe the instrument's integral components, and provide all its key properties. We also sketch the processes needed to understand how NISP operates and is calibrated as well as its technical potentials and limitations. Links to articles providing more details and the technical background are included. The NISP's 16 H2RG detectors with a plate scale of $ deliver a field of view of 0.57\,deg$^2$. In photometric mode, NISP reaches a limiting magnitude of sim \,24.5\,AB\,mag in three photometric exposures of about 100\,s in exposure time for point sources and with a S/N of five. For spectroscopy, NISP's point-source sensitivity is a SNR = 3.5 detection of an emission line with flux sim \,$2 $ integrated over two resolution elements of 13.4\ in 3times 560\,s grism exposures at 1.6\ (redshifted Halpha ). Our calibration includes on-ground and in-flight characterisation and monitoring of the pixel-based detector baseline, dark current, non-linearity, and sensitivity to guarantee a relative photometric accuracy better than 1.5 and a relative spectrophotometry better than 0.7. The wavelength calibration must be accurate to 5\ or better. The NISP is the state-of-the-art instrument in the near-infrared for all science beyond small areas available from HST and JWST -- and it represents an enormous advance from any existing instrumentation due to its combination of field size and high throughput of telescope and instrument. During six-year survey covering 14\,000\,deg$^2$ of extragalactic sky, NISP will be the backbone in determining distances of more than a billion galaxies. Its near-infrared data will become a rich reference imaging and spectroscopy data set for the coming decades.

The near-infrared calibration unit (NI-CU) on board NISP is the first astronomical calibration lamp based on LED to be operated in space. is a mission in ESA's Cosmic Vision 2015--2025 framework to explore the dark universe and provide a next-level characterisation of the nature of gravitation, dark matter, and dark energy. Calibrating photometric and spectrometric measurements of galaxies to better than 1.5 accuracy in a survey homogeneously mapping sim \,14\,000\,deg$^2$ of extragalactic sky requires a very detailed characterisation of NIR detector properties as well as constant monitoring of them in flight. To cover two of the main contributions -- relative pixel-to-pixel sensitivity and non-linearity characteristics -- and to support other calibration activities, NI-CU was designed to provide spatially approximately homogeneous ($<$\,12 variations) and temporally stable illumination (0.1--0.2 over 1200\,s) over the NISP detector plane with minimal power consumption and energy dissipation. NI-CU covers the spectral range sim \,nm -- at cryo-operating temperature -- at five fixed independent wavelengths to capture wavelength-dependent behaviour of the detectors, with fluence over a dynamic range of gtrsim \,100 from sim $. For this functionality, NI-CU is based on LED . We describe the rationale behind the decision and design process, the challenges in sourcing the right LED and the qualification process and lessons learned. We also provide a description of the completed NI-CU, its capabilities, and performance as well as its limits. NI-CU has been integrated into NISP and the satellite, and since launch in July 2023, it has started supporting survey operations.

The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance.

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}  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±0.6(stat)±1(syst)  EeV, 11±2(stat)±1(syst)  EeV, and 31±5(stat)±3(syst)  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.

Dark matter particles could be superheavy, provided their lifetime is much longer than the age of the Universe. Using the sensitivity of the Pierre Auger Observatory to ultrahigh energy neutrinos and photons, we constrain a specific extension of the Standard Model of particle physics that meets the lifetime requirement for a superheavy particle by coupling it to a sector of ultralight sterile neutrinos. Our results show that, for a typical dark coupling constant of 0.1, the mixing angle.m between active and sterile neutrinos must satisfy, roughly,theta(m)less than or similar to 1.5 x 10(-6)(MX=10(9) GeV)(- 2) for a mass MX of the dark-matter particle between 10(8) GeV and 10(11) GeV.

We report an investigation of the mass composition of cosmic rays with energies from 3 to 100 EeV (1 EeV=$10^{18}$ eV) using the distributions of the depth of shower maximum $X_\mathrm{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 data set 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_\mathrm{max}$ distributions up to 100 EeV. Our findings are threefold: (1.) The evolution of the mean logarithmic mass towards a heavier composition with increasing energy can be confirmed and is extended to 100 EeV. (2.) The evolution of the fluctuations of $X_\mathrm{max}$ towards 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_\mathrm{max}$ at the highest energies. (3.) 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.

We test the predictions of hadronic interaction models regarding the depth of maximum of air-shower profiles, X-max, and ground-particle signals in water-Cherenkov detectors at 1000 m from the shower core, S(1000), using the data from the fluorescence and surface detectors of the Pierre Auger Observatory. The test consists of fitting the measured two-dimensional (S(1000), X-max) distributions using templates for simulated air showers produced with hadronic interaction models Epos-LHC, QGSJet-II-04, SIBYLL 2.3d and leaving the scales of predicted X-max and the signals from hadronic component at ground as free-fit parameters. The method relies on the assumption that the mass composition remains the same at all zenith angles, while the longitudinal shower development and attenuation of ground signal depend on the mass composition in a correlated way. The analysis was applied to 2239 events detected by both the fluorescence and surface detectors of the Pierre Auger Observatory with energies between 10(18.5) eV to 10(19.0) eV and zenith angles below 60 degrees. We found, that within the assumptions of the method, the best description of the data is achieved if the predictions of the hadronic interaction models are shifted to deeper X-max values and larger hadronic signals at all zenith angles. Given the magnitude of the shifts and the data sample size, the statistical significance of the improvement of data description using the modifications considered in the paper is larger than 5 sigma even for any linear combination of experimental systematic uncertainties.

To interpret the mean depth of cosmic ray air shower maximum and its dispersion, we parametrize those two observables as functions of the first two moments of the ln A distribution. We examine the goodness of this simple method through simulations of test mass distributions. The application of the parameterization to Pierre Auger Observatory data allows one to study the energy dependence of the mean ln A and of its variance under the assumption of selected hadronic interaction models. We discuss possible implications of these dependences in term of interaction models and astrophysical cosmic ray sources.

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 5s 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(-1) 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.

Large-format infrared detectors are at the heart of major ground and space-based astronomical instruments, and the HgCdTe HxRG is the most widely used. The Near Infrared Spectrometer and Photometer (NISP) of the ESA's Euclid mission launched in July 2023 hosts 16 H2RG detectors in the focal plane. Their performance relies heavily on the effect of image persistence, which results in residual images that can remain in the detector for a long time contaminating any subsequent observations. Deriving a precise model of image persistence is challenging due to the sensitivity of this effect to observation history going back hours or even days. Nevertheless, persistence removal is a critical part of image processing because it limits the accuracy of the derived cosmological parameters. We will present the empirical model of image persistence derived from ground characterization data, adapted to the Euclid observation sequence and compared with the data obtained during the in-orbit calibrations of the satellite.

Euclid, the M2 mission of the ESA's Cosmic Vision 2015-2025 program, aims to explore the Dark Universe by conducting a survey of approximately 14 000 deg(2) and creating a 3D map of the observable Universe of around 1.5 billion galaxies up to redshift z similar to 2. This mission uses two main cosmological probes: weak gravitational lensing and galaxy clustering, leveraging the high-resolution imaging capabilities of the Visual Imaging (VIS) instrument and the photometric and spectroscopic measurements of the Near Infrared Spectrometer and Photometer (NISP) instrument. This paper details some of the activities performed during the commissioning phase of the NISP instrument, following the launch of Euclid on July 1, 2023. In particular, we focus on the calibration of the NISP detectors' baseline and on the performance of a parameter provided by the onboard data processing (called NISP Quality Factor, QF) in detecting the variability of the flux of cosmic rays hitting the NISP detectors. The NISP focal plane hosts sixteen Teledyne HAWAII-2RG (H2RG) detectors. The calibration of these detectors includes the baseline optimization, which optimizes the dynamic range and stability of the signal acquisition. Additionally, this paper investigates the impact of Solar proton flux on the NISP QF, particularly during periods of high Solar activity. Applying a selection criterion on the QF (called NISP QF Proxy), the excess counts are used to monitor the amount of charged particles hitting the NISP detectors. A good correlation was found between the Solar proton flux component above 30 MeV and the NISP QF Proxy, revealing that NISP detectors are not subject to the lower energy components, which are absorbed by the shielding provided by the spacecraft.

The Pierre Auger Observatory is the most sensitive instrument to detect photons with energies above 1017 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 EeV = 10(18) eV) has been developed by combining the fluorescence detector-based measurement of Xmax 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 km(-2) sr(-1) 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.

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 of order $B_{\rm rms}\simeq (50-100)\,{\rm nG}\,(20\rm{Mpc}/{d_{\rm s})( 100\,\rm{kpc}/L_{\rm coh}})^{1/2}$, with $d_{\rm s}$ the typical intersource separation and $L_{\rm coh}$ the 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., $\propto 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.

The Pierre Auger Observatory is the most sensitive instrument to detect photons with energies above $10^{17}$ eV. It measures extensive air showers generated by ultra high 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 EeV = $10^{18}$ 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 km$^{-2}$ sr$^{-1}$ 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 ultra-high energy cosmic rays and allow for the constraint of the mass and lifetime phase space of super-heavy dark matter particles.

We show, for the first time, radio measurements of the depth of shower maximum (X_{max}) of air showers induced by cosmic rays that are compared to measurements of the established fluorescence method at the same location. Using measurements at the Pierre Auger Observatory we show full compatibility between our radio and the previously published fluorescence dataset, and between a subset of air showers observed simultaneously with both radio and fluorescence techniques, a measurement setup unique to the Pierre Auger Observatory. Furthermore, we show radio X_{max} resolution as a function of energy and demonstrate the ability to make competitive high-resolution X_{max} measurements with even a sparse radio array. With this, we show that the radio technique is capable of cosmic-ray mass composition studies, both at Auger and at other experiments.

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.

Operating since 2004, the Pierre Auger Observatory has led to major advances in our understanding of the ultra-high-energy cosmic rays. The latest findings have revealed new insights that led to the upgrade of the Observatory, with the primary goal of obtaining information on the primary mass of the most energetic cosmic rays on a shower-by-shower basis. In the framework of the upgrade, called AugerPrime, the 1660 water-Cherenkov detectors of the surface array are equipped with plastic scintillators and radio antennas, allowing us to enhance the composition sensitivity. To accommodate new detectors and to increase experimental capabilities, the electronics is also upgraded. This includes better timing with up-to-date GPS receivers, higher sampling frequency, increased dynamic range, and more powerful local processing of the data. In this paper, the design characteristics of the new electronics and the enhanced dynamic range will be described. The manufacturing and test processes will be outlined and the test results will be discussed. The calibration of the SD detector and various performance parameters obtained from the analysis of the first commissioning data will also be presented.

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.