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.

AbstractThe first gravitational wave transient GW150914 was observed by Advanced LIGO on September 14th, 2015 at 09:50:45 Universal Time. In addition to follow-up electromagnetic observations, the detection of neutrinos will probe deeply and more on the nature of astrophysical sources, especially in the ultra-high energy regime. Neutrinos in the EeV energy range were searched in data collected at the surface detector of the Pierre Auger Observatory within ± 500 s and 1 day after the GW150914 event. No neutrino candidates were found. Based on this non-observation, we derive the first and only neutrino fluence upper limit at EeV energies for this event at 90% CL, and report constraints on existence of accretion disk around mergers.

We present a search for ultrarelativistic magnetic monopoles with the Pierre Auger observatory. Such particles, possibly a relic of phase transitions in the early Universe, would deposit a large amount of energy along their path through the atmosphere, comparable to that of ultrahigh-energy cosmic rays (UHECRs). The air-shower profile of a magnetic monopole can be effectively distinguished by the fluorescence detector from that of standard UHECRs. No candidate was found in the data collected between 2004 and 2012, with an expected background of less than 0.1 event from UHECRs. The corresponding 90% confidence level (C.L.) upper limits on the flux of ultrarelativistic magnetic monopoles range from 10(-1)9 (cm(2) sr s)(-1) for a Lorentz factor gamma = 10(9) to 2.5 x 10(-21) (cm(2) sr s)(-1) for gamma = 10(12). These results-the first obtained with a UHECR detector-improve previously published limits by up to an order of magnitude.

We report a first measurement for ultra-high energy cosmic rays of the correlation between the depth of shower maximum and the signal in the water Cherenkov stations of air-showers registered simultaneously by the fluorescence and the surface detectors of the Pierre Auger Observatory. Such a correlation measurement is a unique feature of a hybrid air-shower observatory with sensitivity to both the electromagnetic and muonic components. It allows an accurate determination of the spread of primary masses in the cosmic-ray flux. Up till now, constraints on the spread of primary masses have been dominated by systematic uncertainties. The present correlation measurement is not affected by systematics in the measurement of the depth of shower maximum or the signal in the water Cherenkov stations. The analysis relies on general characteristics of air showers and is thus robust also with respect to uncertainties in hadronic event generators. The observed correlation in the energy range around the `ankle' at $\lg(E/{\rm eV})=18.5-19.0$ differs significantly from expectations for pure primary cosmic-ray compositions. A light composition made up of proton and helium only is equally inconsistent with observations. The data are explained well by a mixed composition including nuclei with mass $A > 4$. Scenarios such as the proton dip model, with almost pure compositions, are thus disfavoured as the sole explanation of the ultrahigh-energy cosmic-ray flux at Earth.

To exploit the full potential of radio measurements of cosmic-ray air showers at MHz frequencies, a detector timing synchronization within 1 ns is needed. Large distributed radio detector arrays such as the Auger Engineering Radio Array (AERA) rely on timing via the Global Positioning System (GPS) for the synchronization of individual detector station clocks. Unfortunately, GPS timing is expected to have an accuracy no better than about 5 ns. In practice, in particular in AERA, the GPS clocks exhibit drifts on the order of tens of ns. We developed a technique to correct for the GPS drifts, and an independent method is used to cross-check that indeed we reach a nanosecond-scale timing accuracy by this correction. First, we operate a ``beacon transmitter'' which emits defined sine waves detected by AERA antennas recorded within the physics data. The relative phasing of these sine waves can be used to correct for GPS clock drifts. In addition to this, we observe radio pulses emitted by commercial airplanes, the position of which we determine in real time from Automatic Dependent Surveillance Broadcasts intercepted with a software-defined radio. From the known source location and the measured arrival times of the pulses we determine relative timing offsets between radio detector stations. We demonstrate with a combined analysis that the two methods give a consistent timing calibration with an accuracy of 2 ns or better. Consequently, the beacon method alone can be used in the future to continuously determine and correct for GPS clock drifts in each individual event measured by AERA.

Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to extend its range of detection and to directly measure the muon content of the particle showers. It consists of an infill of surface water-Cherenkov detectors accompanied by buried scintillator detectors used for muon counting. The main objectives of the AMIGA engineering array, referred to as the Unitary Cell, are to identify and resolve all engineering issues as well as to understand the muon-number counting uncertainties related to the design of the detector. The mechanical design, fabrication and deployment processes of the muon counters of the Unitary Cell are described in this document. These muon counters modules comprise sealed PVC casings containing plastic scintillation bars, wavelength-shifter optical fibers, 64 pixel photomultiplier tubes, and acquisition electronics. The modules are buried approximately 2.25 m below ground level in order to minimize contamination from electromagnetic shower particles. The mechanical setup, which allows access to the electronics for maintenance, is also described in addition to tests of the modules' response and integrity. The completed Unitary Cell has measured a number of air showers of which a first analysis of a sample event is included here.

Ultrahigh energy cosmic ray air showers probe particle physics at energies beyond the reach of accelerators. Here we introduce a new method to test hadronic interaction models without relying on the absolute energy calibration, and apply it to events with primary energy 6-16 EeV (E-CM = 110-170 TeV), whose longitudinal development and lateral distribution were simultaneously measured by the Pierre Auger Observatory. The average hadronic shower is 1.33 +/- 0.16 (1.61 +/- 0.21) times larger than predicted using the leading LHC-tuned models EPOS-LHC (QGSJetII-04), with a corresponding excess of muons.

The azimuthal asymmetry in the risetime of signals in Auger surface detector stations is a source of information on shower development. The azimuthal asymmetry is due to a combination of the longitudinal evolution of the shower and geometrical effects related to the angles of incidence of the particles into the detectors. The magnitude of the effect depends upon the zenith angle and state of development of the shower and thus provides a novel observable, $(\sec \theta)_\mathrm{max}$, sensitive to the mass composition of cosmic rays above $3 \times 10^{18}$ eV. By comparing measurements with predictions from shower simulations, we find for both of our adopted models of hadronic physics (QGSJETII-04 and EPOS-LHC) an indication that the mean cosmic-ray mass increases slowly with energy, as has been inferred from other studies. However, the mass estimates are dependent on the shower model and on the range of distance from the shower core selected. Thus the method has uncovered further deficiencies in our understanding of shower modelling that must be resolved before the mass composition can be inferred from $(\sec \theta)_\mathrm{max}$.

On September 14, 2015 the Advanced LIGO detectors observed their first gravitational wave (GW) transient GW150914. This was followed by a second GW event observed on December 26, 2015. Both events were inferred to have arisen from the merger of black holes in binary systems. Such a system may emit neutrinos if there are magnetic fields and disk debris remaining from the formation of the two black holes. With the surface detector array of the Pierre Auger Observatory we can search for neutrinos with energy E-nu above 100 PeV from pointlike sources across the sky with equatorial declination from about -65 degrees to +60 degrees, and, in particular, from a fraction of the 90% confidence-level inferred positions in the sky of GW150914 and GW151226. A targeted search for highly inclined extensive air showers, produced either by interactions of downward-going neutrinos of all flavors in the atmosphere or by the decays of tau leptons originating from tau-neutrino interactions in the Earth's crust (Earth-skimming neutrinos), yielded no candidates in the Auger data collected within +/- 500 s around or 1 day after the coordinated universal time (UTC) of GW150914 and GW151226, as well as in the same search periods relative to the UTC time of the GW candidate event LVT151012. From the nonobservation we constrain the amount of energy radiated in ultrahigh-energy neutrinos from such remarkable events.

The Pierre Auger Observatory has begun a major Upgrade of its already impressive capabilities, with an emphasis on improved mass composition determination using the surface detectors of the Observatory. Known as AugerPrime, the upgrade will include new 4 m$^2$ plastic scintillator detectors on top of all 1660 water-Cherenkov detectors, updated and more flexible surface detector electronics, a large array of buried muon detectors, and an extended duty cycle for operations of the fluorescence detectors. This Preliminary Design Report was produced by the Collaboration in April 2015 as an internal document and information for funding agencies. It outlines the scientific and technical case for AugerPrime. We now release it to the public via the arXiv server. We invite you to review the large number of fundamental results already achieved by the Observatory and our plans for the future.

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±0.7(stat)±6.7(syst)  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-principles 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.

The DRIFT-IId dark matter detector is a m3-scale low-pressure TPC with directional sensitivity to WIMP-induced nuclear recoils. Its primary backgrounds were due to alpha decays from contamination on the central cathode. Efforts to reduce these backgrounds led to replacing the 20μm wire central cathode with one constructed from 0.9μm aluminized mylar, which is almost totally transparent to alpha particles. Detailed modeling of the nature and origin of the remaining backgrounds led to an in-situ, ppt-sensitive assay of alpha decay backgrounds from the central cathode. This led to further improvements in the thin-film cathode resulting in over 2 orders of magnitude reduction in backgrounds compared to the wire cathode. Finally, the addition of O2 to CS2 gas was found to produce multiple species of electronegative charge carriers, providing a method to determine the absolute position of nuclear recoils and reject all known remaining backgrounds while retaining a high efficiency for nuclear recoil detection.

We present the first hybrid measurement of the average muon number in air showers at ultrahigh energies, initiated by cosmic rays with zenith angles between 62 degrees and 80 degrees. The measurement is based on 174 hybrid events recorded simultaneously with the surface detector array and the fluorescence detector of the Pierre Auger Observatory. The muon number for each shower is derived by scaling a simulated reference profile of the lateral muon density distribution at the ground until it fits the data. A 10(19) eV shower with a zenith angle of 67 degrees, which arrives at the surface detector array at an altitude of 1450 m above sea level, contains on average (2.68 +/- 0.04 +/- 0.48(sys)) x 10(7) muons with energies larger than 0.3 GeV. The logarithmic gain d lnN mu/d ln E of muons with increasing energy between 4 x 10(18) eV and 5 x 10(19) eV is measured to be (1.029 +/- 0.024) 0.030(sys)).

Background events in the DRIFT-Ild dark matter detector, mimicking potential WIMP signals, are predominantly caused by alpha decays on the central cathode in which the alpha particle is completely or partially absorbed by the cathode material. We installed a 0.9 mu m thick aluminized-mylar cathode as a way to reduce the probability of producing these backgrounds. We study three generations of cathode (wire, thin-film, and radiologically clean thin-film) with a focus on the ratio of background events to alpha decays. Two independent methods of measuring the absolute alpha decay rate are used to ensure an accurate result, and agree to within 10%. Using alpha range spectroscopy, we measure the radiologically cleanest cathode version to have a contamination of 3.3 +/- 0.1 ppt U-234 and 73 +/- 2 ppb U-238. This cathode reduces the probability of producing an RPR from an alpha decay by a factor of 70 +/- 20 compared to the original stainless steel wire cathode. First results are presented from a texturized version of the cathode, intended to be even more transparent to alpha particles. These efforts, along with other background reduction measures, have resulted in a drop in the observed background rate from 500/day to With the recent implementation of full volume fiducialization, these remaining background events are identified, allowing for background free operation. Published by Elsevier B.V.

The energy spectrum above 4×10 eV is presented, obtained from 5936 events with zenith angles exceeding 62◦ collected by the Surface Detector of the Pierre Auger Observatory from 1 January 2004 to 31 December 2010. Showers with such large zenith angles are muon-dominated at ground level and the radial symmetry around the shower axis is broken due to geomagnetic deflections. They are analysed separately from showers with smaller zenith angles using two-dimensional models of the muon density at ground, allowing one to reconstruct a global muon number for every event. The conversion of the muon number to energy is obtained using the sub-sample of events detected simultaneously with both the Surface and the Fluorescence Detector. The spectrum obtained displays suppression near 4 × 10 eV compatible with the analysis that uses less inclined events.

The addition of O2 to gas mixtures in time projection chambers containing CS2 has recently been shown to produce multiple negative ions that travel at slightly different velocities. This allows a measurement of the absolute position of ionising events in the z (drift) direction. In this work, we apply the z-fiducialisation technique to a directional dark matter search. In particular, we present results from a 46.3 live-day source-free exposure of the DRIFT-IId detector run in this completely new mode. With full-volume fiducialisation, we have achieved the first background-free operation of a directional detector. The resulting exclusion curve for spin-dependent WIMP-proton interactions reaches 1.1 pb at 100 GeV/c2, a factor of 2 better than our previous work. We describe the automated analysis used here, and argue that detector upgrades, implemented after the acquisition of these data, will bring an additional factor of >3 improvement in the near future.

The optical calibration of the fluorescence telescopes is a significant contribution to the overall uncertainty of energy measurements made by the Pierre Auger Observatory and Telescope Array Project. Some sources of uncertainty, such as the fluorescence yield in air, affect both experiments similarly. However, the optical calibration of the fluorescence telescopes is a source of independent uncertainty. The Pierre Auger and Telescope Array collaborations have taken initial steps to establish a relative end-to-end optical calibration of the fluorescence telescopes. An Octocopter carrying a portable light source has been flown in front of fluorescence telescopes at both Pierre Auger and Telescope Array sites. Laboratory calibration measurements of the light source before and after the flights provide a common baseline for the relative end-to-end calibration. We expect this system will lead to a common photonic calibration for both experiments. After giving a brief description of the UV light source and the Octocopter used for the measurements, the parameters and the calibration procedures for the light source will be discussed. First results on telescope images of the light source will be presented.

We present the results of an analysis of the large angular scale distribution of the arrival directions of cosmic rays with energy above 4 EeV detected at the Pierre Auger Observatory including for the first time events with zenith angle between 60° and 80°. We perform two Rayleigh analyses, one in the right ascension and one in the azimuth angle distributions, that are sensitive to modulations in right ascension and declination, respectively. The largest departure from isotropy appears in the E > 8 ?> EeV energy bin, with an amplitude for the first harmonic in right ascension r 1 &agr; = ( 4.4 ± 1.0 ) × 10 − 2 ?> , that has a chance probability P ( ≥ r 1 &agr; ) = 6.4 × 10 − 5 ?> , reinforcing the hint previously reported with vertical events alone.