Extensive air showers, originating from ultra-high energy cosmic rays, have been successfully measured through the use of arrays of water-Cherenkov detectors (WCDs). Sophisticated analyses exploiting WCD data have made it possible to demonstrate that shower simulations, based on different hadronic-interaction models, cannot reproduce the observed number of muons at the ground. The accurate knowledge of the WCD response to muons is paramount in establishing the exact level of this discrepancy. In this work, we report on a study of the response of a WCD of the Pierre Auger Observatory to atmospheric muons performed with a hodoscope made of resistive plate chambers (RPCs), enabling us to select and reconstruct nearly 600 thousand single muon trajectories with zenith angles ranging from 0$^\circ$ to 55$^\circ$. Comparison of distributions of key observables between the hodoscope data and the predictions of dedicated simulations allows us to demonstrate the accuracy of the latter at a level of 2%. As the WCD calibration is based on its response to atmospheric muons, the hodoscope data are also exploited to show the long-term stability of the procedure.

We search for signals of magnetically-induced effects in the arrival directions of ultra-high-energy cosmic rays detected at the Pierre Auger Observatory. We apply two different methods. One is a search for sets of events that show a correlation between their arrival direction and the inverse of their energy, which would be expected if they come from the same point-like source, they have the same electric charge and their deflection is relatively small and coherent. We refer to these sets of events as "multiplets". The second method, called "thrust", is a principal axis analysis aimed to detect the elongated patterns in a region of interest. We study the sensitivity of both methods using a benchmark simulation and we apply them to data in two different searches. The first search is done assuming as source candidates a list of nearby active galactic nuclei and starburst galaxies. The second is an all-sky blind search. We report the results and we find no statistically significant features. We discuss the compatibility of these results with the indications on the mass composition inferred from data of the Pierre Auger Observatory.

We present measurements of the large-scale cosmic-ray anisotropies in right ascension, using data collected by the surface detector array of the Pierre Auger Observatory over more than 14 years. We determine the equatorial dipole component, $\vec{d}_\perp$, through a Fourier analysis in right ascension that includes weights for each event so as to account for the main detector-induced systematic effects. For the energies at which the trigger efficiency of the array is small, the ``East-West'' method is employed. Besides using the data from the array with detectors separated by 1500 m, we also include data from the smaller but denser sub-array of detectors with 750 m separation, which allows us to extend the analysis down to $\sim 0.03$ EeV. The most significant equatorial dipole amplitude obtained is that in the cumulative bin above 8~EeV, $d_\perp=6.0^{+1.0}_{-0.9}$%, which is inconsistent with isotropy at the 6$\sigma$ level. In the bins below 8 EeV, we obtain 99% CL upper-bounds on $d_\perp$ at the level of 1 to 3 percent. At energies below 1 EeV, even though the amplitudes are not significant, the phases determined in most of the bins are not far from the right ascension of the Galactic center, at $\alpha_{\rm GC}=-94^\circ$, suggesting a predominantly Galactic origin for anisotropies at these energies. The reconstructed dipole phases in the energy bins above 4 EeV point instead to right ascensions that are almost opposite to the Galactic center one, indicative of an extragalactic cosmic ray origin.

The Fluorescence detector Array of Single-pixel Telescopes (FAST) is a design concept for a next-generation UHECR observatory, addressing the requirements for a large-area, low-cost detector suitable for measuring the properties of ultra-high energy cosmic rays (UHECRs), having energies exceeding $10^{19.5}$\,eV, with an unprecedented aperture. We have developed a full-scale prototype consisting of four 200\,mm diameter photo-multiplier tubes at the focus of a segmented mirror of 1.6\,m in diameter. In October 2016, September 2017, and September 2018 we installed three such prototypes at the Black Rock Mesa site of the Telescope Array experiment in central Utah, USA. All three telescopes have been steadily taking data since installation. We report on the design and installation of these prototypes, and present some preliminary results, including measurements of artificial light sources, distant ultraviolet lasers, and UHECRs. Furthermore, we discuss some additional uses for these simplified low-cost fluorescence telescopes, including the facilitation of a systematic comparison of the transparency of the atmosphere above the Telescope Array experiment and the Pierre Auger Observatory, a study of the systematic uncertainty associated with the existing fluorescence detectors of these two experiments, and a cross-calibration of their energy and $X_{\text{max}}$ scales.

With the Surface Detector array (SD) of the Pierre Auger Observatory we can detect neutrinos with energy between $10^{17}\,$eV and $10^{20}\,$eV from point-like sources across the sky, from close to the Southern Celestial Pole up to $60^\circ$ in declination, with peak sensitivities at declinations around $\sim -53^\circ$ and $\sim+55^\circ$, and an unmatched sensitivity for arrival directions in the Northern hemisphere. A search has been performed for highly-inclined air showers induced by neutrinos of all flavours with no candidate events found in data taken between 1 Jan 2004 and 31 Aug 2018. Upper limits on the neutrino flux from point-like steady sources have been derived as a function of source declination. An unrivaled sensitivity is achieved in searches for transient sources with emission lasting over an hour or less, if they occur within the field of view corresponding to the zenith angle range between $60^\circ$ and $~95^\circ$ where the SD of the Pierre Auger Observatory is most sensitive to neutrinos.

The determination of the primary energy of extensive air showers using the fluorescence detection technique requires an estimation of the energy carried away by particles that do not deposit all their energy in the atmosphere. This estimation is typically made using Monte Carlo simulations and thus depends on the assumed primary particle mass and on model predictions for neutrino and muon production. In this work we present a new method to obtain the invisible energy from events detected by the Pierre Auger Observatory. The method uses measurements of the muon number at ground level, and it allows us to reduce significantly the systematic uncertainties related to the mass composition and the high energy hadronic interaction models, and consequently to improve the estimation of the energy scale of the Observatory.

Neutrinos with energies above $10^{17}$ eV are detectable with the Surface Detector Array of the Pierre Auger Observatory. The identification is efficiently performed for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for Earth-skimming $\tau$ neutrinos with nearly tangential trajectories relative to the earth. No neutrino candidates were found in $\sim\,14.7$ years of data taken up to 31 August 2018. This leads to restrictive upper bounds on their flux. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high-energy neutrinos with an $E_\nu^{-2}$ spectrum in the energy range $1.0 \times 10^{17}~{\rm eV} - 2.5 \times 10^{19}~{\rm eV}$ is $E^2 {\rm d}N_\nu/{\rm d}E_\nu < 4.4 \times 10^{-9}~{\rm GeV~cm^{-2}~s^{-1}~sr^{-1}}$, placing strong constraints on several models of neutrino production at EeV energies and on the properties of the sources of ultra-high-energy cosmic rays.

Contributions of the Pierre Auger Collaboration to the 36th International Cosmic Ray Conference (ICRC 2019), 24 July - 1 August 2019, Madison, Wisconsin, USA.

Advanced knowledge of the detailed atmospheric properties of both the future sites of the Cherenkov Telescope Array is essential in preparation of the arrival of the first scientific data. Meteorological variables are studied using a dedicated characterization station installed at the southern site in Chile and a wealth of data from existing observatories around the northern site on the La Palma island. Campaigns using radiosondes launched on balloons are foreseen to complement these data in the near future. Cloudiness during the night has been continuously monitored at both sites for several years using All-sky Cameras which assess the presence of clouds based on detection of stars. The integrated aerosol optical depth over the southern site has been measured using a Sun/Moon Photometer since 2016 and the small robotic FRAM telescope since 2017; identical instruments have been deployed at the northern site in autumn 2018. Also in October 2018, the ARCADE Raman lidar (RL) has started to take measurements on routine basis at the northern site, providing data on the vertical profile of the aerosol optical properties (i.e., extinction and scattering) and of the water vapour mixing ratio. We present the data currently available from these instruments from both sites with emphasis on characteristics important for the (future) operation of Imaging Atmospheric Cherenkov Telescopes.

The atmospheric monitoring devices for the planned calibration system of the Cherenkov Telescope Array (CTA) are undergoing intensive development, prototyping and testing. The All-Sky Cameras, the Sun/Moon Photometers and the FRAM telescopes have been gradually deployed at the future CTA sites with the primary goal of site characterization, simultaneously allowing the assessment of their operational reliability in realistic environmental conditions. All three devices have shown the ability to work smoothly in both the extreme dryness and the large temperature variations of the southern site as well as in the occasional adverse weather during winter months at the northern site. The target availability of 95\% of time has not yet been reached mostly due to minor hardware failures that have proven difficult to fix due to the remoteness of the installation in the absence of the future CTA infrastructure. The experience gathered during the prototype operations will contribute to the improved reliability of the final instruments. The Raman LIDARs, described in separate proceedings of this conference, and the infrared Ceilometer, ready for testing in Prague, will complement the set of atmospheric calibration devices in near future. The final operational procedures for the atmospheric calibration of the CTA during its operation are being finalized foreseeing the use of the All-sky Cameras and the Ceilometer for the monitoring of clouds over the whole sky and the LIDARs and FRAMs for precision measurements of the atmospheric transmission as a function of altitude and position within the field-of-view of the CTA array.

The next generation of Cherenkov telescope cameras feature Silicon Photo Multipliers (SiPM), which can guarantee excellent performance and allow for observation also under moonlight, increasing duty-cycle and therefore the physics reach. A 4 m-diameter Davies-Cotton prototype telescope with a 9-degree optical FoV and a 1296-pixel SiPM camera, has been designed to meet the requirements of the next generation of ground-based gamma-ray observatories at the highest energies. The large-scale production of the telescopes for array deployment has required the development of a fully automated calibration strategy which relies on a dedicated hardware, the Camera Test Setup (CTS). For each camera pixel, the CTS is equipped with two LEDs, one operated in pulsed mode to reproduce signal and one in continuous mode to reproduce night-sky background. In this contribution we will present the camera calibration strategy, from the laboratory measurement to the on-site monitoring with emphasis on the results obtained with the first camera prototype. In addition, key performances such as charge resolution, time resolution and trigger efficiencies and their degradation with increasing night-sky background level will be presented.

The origin and nature of ultra-high energy cosmic rays (UHECRs) are hot topics in the astroparticle physics community. The Fluorescence detector Array of Single-pixel Telescopes (FAST) is a design for a next-generation ground-based UHECR observatory, addressing the requirements for a large-area, low-cost detector suitable for measuring the properties of the highest energy cosmic rays with an unprecedented aperture. We have developed a full-scale prototype consisting of four 200 mm photomultiplier tubes at the focus of a segmented mirror of 1.6 m in diameter. Over the last three years, we have installed three prototypes at the Telescope Array Experiment in Utah, USA. These telescopes have been steadily taking data since installation. We report on preliminary results of the full-scale FAST prototypes, including measurements of UHECRs, and distant ultra-violet lasers used to study the atmospheric transparency. Furthermore, we discuss the installation of an additional identical FAST prototype at the Pierre Auger Observatory in Argentina. Possible benefits to the Telescope Array Experiment and the Pierre Auger Observatory include a comparison of the transparency of the atmosphere above both experiments, a study of the systematic uncertainty associated with their existing fluorescence detectors, and a cross-calibration of their energy and Xmax scales.

The SST-1M telescope was developed as a prototype of a Small-Size-Telescope for the Cherenkov Telescope Array observatory and it has been extensively tested in Krakow since 2017. In this contribution we present validation of the Monte Carlo model of the prototype and expected performance in Krakow conditions. We focus on gamma/hadron separation and mono reconstruction of energy and gamma photon arrival direction using Machine learning methods.

The profile of the longitudinal development of showers produced by ultra-high energy cosmic rays carries information related to the interaction properties of the primary particles with atmospheric nuclei. In this work, we present the first measurement of the average shower profile in traversed atmospheric depth at the Pierre Auger Observatory. The shapes of profiles are well reproduced by the Gaisser-Hillas parametrization within the range studied, for E > 10^{17.8} eV. A detailed analysis of the systematic uncertainties is performed using 10 years of data and a full detector simulation. The average shape is quantified using two variables related to the width and asymmetry of the profile, and the results are compared with predictions of hadronic interaction models for different primary particles.

An overview of the multi-messenger capabilities of the Pierre Auger Observatory is presented. The techniques and performance of searching for Ultra-High Energy neutrinos, photons and neutrons are described. Some of the most relevant results are reviewed, such as stringent upper bounds that were placed to a flux of diffuse cosmogenic neutrinos and photons, bounds placed on neutrinos emitted from compact binary mergers that were detected by LIGO and Virgo during their first and second observing runs, as well as searches for high energy photons and neutrons from the Galactic center that constrain the properties of the putative Galactic PeVatron, observed by the H.E.S.S.\ collaboration. The observation of directional correlations between ultra-high energy cosmic rays and either high energy astrophysical neutrinos or specific source populations, weighted by their electromagnetic radiation, are also discussed. They constitute additional multi-messenger approaches aimed at identifying the sources of high energy cosmic rays.

We present a detailed study of the large-scale anisotropies of cosmic rays with energies above 4 EeV measured using the Pierre Auger Observatory. For the energy bins [4,8] EeV and $E\geq 8$ EeV, the most significant signal is a dipolar modulation in right ascension at energies above 8 EeV, as previously reported. In this paper we further scrutinize the highest-energy bin by splitting it into three energy ranges. We find that the amplitude of the dipole increases with energy above 4 EeV. The growth can be fitted with a power law with index $\beta=0.79\pm 0.19$. The directions of the dipoles are consistent with an extragalactic origin of these anisotropies at all the energies considered. Additionally we have estimated the quadrupolar components of the anisotropy: they are not statistically significant. We discuss the results in the context of the predictions from different models for the distribution of ultrahigh-energy sources and cosmic magnetic fields.

With the Auger Engineering Radio Array (AERA) of the Pierre Auger Observatory, we have observed the radio emission from 561 extensive air showers with zenith angles between 60$^\circ$ and 84$^\circ$. In contrast to air showers with more vertical incidence, these inclined air showers illuminate large ground areas of several km$^2$ with radio signals detectable in the 30 to 80\,MHz band. A comparison of the measured radio-signal amplitudes with Monte Carlo simulations of a subset of 50 events for which we reconstruct the energy using the Auger surface detector shows agreement within the uncertainties of the current analysis. As expected for forward-beamed radio emission undergoing no significant absorption or scattering in the atmosphere, the area illuminated by radio signals grows with the zenith angle of the air shower. Inclined air showers with EeV energies are thus measurable with sparse radio-antenna arrays with grid sizes of a km or more. This is particularly attractive as radio detection provides direct access to the energy in the electromagnetic cascade of an air shower, which in case of inclined air showers is not accessible by arrays of particle detectors on the ground.