Measurement of the cosmic ray spectrum above 4×10184{\times}10^{18} eV using inclined events detected with the Pierre Auger Observatory

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.Comment: Replaced with published version. Added journal reference and DO

Energy Estimation of Cosmic Rays with the Engineering Radio Array of the Pierre Auger Observatory

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.Comment: Replaced with published version. Added journal reference and DO

Measurement of the Radiation Energy in the Radio Signal of Extensive Air Showers as a Universal Estimator of Cosmic-Ray Energy

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.Comment: Replaced with published version. Added journal reference and DOI. Supplemental material in the ancillary file

Multi-messenger observations of a binary neutron star merger

On 2017 August 17 a binary neutron star coalescence candidate (later designated GW170817) with merger time 12:41:04 UTC was observed through gravitational waves by the Advanced LIGO and Advanced Virgo detectors. The Fermi Gamma-ray Burst Monitor independently detected a gamma-ray burst (GRB 170817A) with a time delay of ~1.7 s with respect to the merger time. From the gravitational-wave signal, the source was initially localized to a sky region of 31 deg2 at a luminosity distance of 40+8-8 Mpc and with component masses consistent with neutron stars. The component masses were later measured to be in the range 0.86 to 2.26 Mo. An extensive observing campaign was launched across the electromagnetic spectrum leading to the discovery of a bright optical transient (SSS17a, now with the IAU identification of AT 2017gfo) in NGC 4993 (at ~40 Mpc) less than 11 hours after the merger by the One- Meter, Two Hemisphere (1M2H) team using the 1 m Swope Telescope. The optical transient was independently detected by multiple teams within an hour. Subsequent observations targeted the object and its environment. Early ultraviolet observations revealed a blue transient that faded within 48 hours. Optical and infrared observations showed a redward evolution over ~10 days. Following early non-detections, X-ray and radio emission were discovered at the transient’s position ~9 and ~16 days, respectively, after the merger. Both the X-ray and radio emission likely arise from a physical process that is distinct from the one that generates the UV/optical/near-infrared emission. No ultra-high-energy gamma-rays and no neutrino candidates consistent with the source were found in follow-up searches. These observations support the hypothesis that GW170817 was produced by the merger of two neutron stars in NGC4993 followed by a short gamma-ray burst (GRB 170817A) and a kilonova/macronova powered by the radioactive decay of r-process nuclei synthesized in the ejecta

Combined fit to the spectrum and composition data measured by the Pierre Auger Observatory including magnetic horizon effects

The measurements by the Pierre Auger Observatory of the energy spectrum and mass composition of cosmic rays can be interpreted assuming the presence of two extragalactic source populations, one dominating the flux at energies above a few EeV and the other below. To fit the data ignoring magnetic field effects, the high-energy population needs to accelerate a mixture of nuclei with very hard spectra, at odds with the approximate E$^{-2}$ shape expected from diffusive shock acceleration. The presence of turbulent extragalactic magnetic fields in the region between the closest sources and the Earth can significantly modify the observed CR spectrum with respect to that emitted by the sources, reducing the flux of low-rigidity particles that reach the Earth. We here take into account this magnetic horizon effect in the combined fit of the spectrum and shower depth distributions, exploring the possibility that a spectrum for the high-energy population sources with a shape closer to E$^{-2}$ be able to explain the observations

Studies of the mass composition of cosmic rays and proton-proton interaction cross-sections at ultra-high energies with the Pierre Auger Observatory

In this work, we present an estimate of the cosmic-ray mass composition from the distributions of the depth of the shower maximum (Xmax) measured by the fluorescence detector of the Pierre Auger Observatory. We discuss the sensitivity of the mass composition measurements to the uncertainties in the properties of the hadronic interactions, particularly in the predictions of the particle interaction cross-sections. For this purpose, we adjust the fractions of cosmic-ray mass groups to fit the data with Xmax distributions from air shower simulations. We modify the proton-proton cross-sections at ultra-high energies, and the corresponding air shower simulations with rescaled nucleus-air cross-sections are obtained via Glauber theory. We compare the energy-dependent composition of ultra-high-energy cosmic rays obtained for the different extrapolations of the proton-proton cross-sections from low-energy accelerator data

Study of downward Terrestrial Gamma-ray Flashes with the surface detector of the Pierre Auger Observatory

The surface detector (SD) of the Pierre Auger Observatory, consisting of 1660 water-Cherenkov detectors (WCDs), covers 3000 km2 in the Argentinian pampa. Thanks to the high efficiency of WCDs in detecting gamma rays, it represents a unique instrument for studying downward Terrestrial Gamma-ray Flashes (TGFs) over a large area. Peculiar events, likely related to downward TGFs, were detected at the Auger Observatory. Their experimental signature and time evolution are very different from those of a shower produced by an ultrahigh-energy cosmic ray. They happen in coincidence with low thunderclouds and lightning, and their large deposited energy at the ground is compatible with that of a standard downward TGF with the source a few kilometers above the ground. A new trigger algorithm to increase the TGF-like event statistics was installed in the whole array. The study of the performance of the new trigger system during the lightning season is ongoing and will provide a handle to develop improved algorithms to implement in the Auger upgraded electronic boards. The available data sample, even if small, can give important clues about the TGF production models, in particular, the shape of WCD signals. Moreover, the SD allows us to observe more than one point in the TGF beam, providing information on the emission angle

The dynamic range of the upgraded surface-detector stations of AugerPrime

The detection of ultra-high-energy cosmic rays by means of giant detector arrays is often limited by the saturation of the recorded signals near the impact point of the shower core at the ground, where the particle density dramatically increases. The saturation affects in particular the highest energy events, worsening the systematic uncertainties in the reconstruction of the shower characteristics. The upgrade of the Pierre Auger Observatory, called AugerPrime, includes the installation of an 1-inch Small PhotoMultiplier Tube (SPMT) inside each water-Cherenkov station (WCD) of the surface detector array. The SPMT allows an unambiguous measurement of signals down to about 250m from the shower core, thus reducing the number of events featuring a saturated station to a negligible level. In addition, a 3.8m2 plastic scintillator (Scintillator Surface Detector, SSD) is installed on top of each WCD. The SSD is designed to match the WCD (with SPMT) dynamic range, providing a complementary measurement of the shower components up to the highest energies. In this work, the design and performances of the upgraded AugerPrime surface-detector stations in the extended dynamic range are described, highlighting the accuracy of the measurements. A first analysis employing the unsaturated signals in the event reconstruction is also presented