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 distribution of depth in which a cosmic ray air shower reaches its maximum number of particles ($X_{max}$) is studied and parametrized. Three functions are studied for proton, carbon, silicon, and iron primary particles with energies ranging from $10^{17}$ eV to $10^{20}$ eV for three hadronic interaction models: EPOS-LHC, QGSJetII.04, and Sibyll2.3c. The function which best describes the $X_{max}$ distribution of a mixed composition is also studied. A very large number of simulated showers and a detailed analysis procedure were used to guarantee negligible effects of undersampling and of fitting in the final results. For the first time, a comparison of several functions is presented under the same assumption with the intention of selecting the best functional form to describe the $X_{max}$ distribution. The Generalized Gumbel distribution is shown to be the best option for a general description of all cases. The Log-normal distribution is also a good choice for some cases while the Exponentially Modified Gaussian distribution has shown to be the worst choice in almost all cases studied. All three functions are parametrized as a function of energy and primary mass.

The understanding of the basic properties of the ultra - high energy extensive air showers is strongly dependent on the description of the hadronic interactions in a energy range beyond that probed by the LHC. One of the uncertainties present in the modeling of the air showers is the treatment of diffractive interactions, which are dominated by non - perturbative physics and usually described by phenomenological models. These interactions are expect to affect the development of the air showers, since they provide a way of transporting substantial amounts of energy deep in the atmosphere, modifying the global characteristics of the shower profile. In this paper we investigate the impact of the diffractive interactions in the observables that can be measured in hadronic collisions at high energies and ultra - high energy cosmic ray interactions. We consider three distinct phenomenological models for the treatment of diffractive physics and estimate the influence of these interactions on the elasticity, number of secondaries, longitudinal air shower profiles and muon densities for proton - air and iron - air collisions at different primary energies. Our results demonstrate that the diffractive events has a non - negligible effect on the observables and that the distinct approaches for these interactions, present in the phenomenological models, are an important source of theoretical uncertainty for the description of the extensive air showers.