ForewordThe 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 1 . 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.The Pierre Auger Collaboration 1 As a result of continuing R&D, slight changes have been implemented in the baseline design since this Report was written. These changes will be documented in a forthcoming Technical Design Report. ix x Executive Summary Present Results from the Pierre Auger ObservatoryMeasurements of the Auger Observatory have dramatically advanced our understanding of ultra-high energy cosmic rays. The suppression of the flux around 5×10 19 eV is now confirmed without any doubt. Strong limits have been placed on the photon and neutrino components of the flux indicating that "top-down" source processes, such as the decay of superheavy particles, cannot account for a significant part of the observed particle flux. A largescale dipole anisotropy of ∼7% amplitude has been found for energies above 8×10 18 eV. In addition there is also an indication of the presence of a large scale anisotropy below the ankle. Particularly exciting is the observed behavior of the depth of shower maximum with energy, which changes in an unexpected, non-trivial way. Around 3×10 18 eV it shows a distinct change of slope with energy, and the shower-to-shower variance decreases. Interpreted with the leading LHC-tuned shower models, this implies a gradual shift to a heavier composition. A number of fundamentally different astrophysical model scenarios have been developed to describe this evolution. The high degree of isotropy observed in numerous tests of the small-scale angular distribution of UHECR above 4×10 19 eV is remarkable, challenging original expectations that assumed only a few cosmic ray sources with a light composition at the highest energies. Interestingly, the largest departures from isotropy are observed for cosmic rays with E > 5.8×10 19 eV in ∼20 • sky-windows. Due to a duty cycle of ∼15% of the fluorescence telescopes, the data on the depth of shower maximum extend only up to the flux suppression region, i.e. 4×10 19 eV. Obtaining more information on the composition of cosmic rays at higher energies will provide crucial means to discriminate between the model classes and to understand the origin of the observed flux suppre...

We present the simulation framework CRPropa version 3 designed for efficient development of astrophysical predictions for ultra-high energy particles. Users can assemble modules of the most relevant propagation effects in galactic and extragalactic space, include their own physics modules with new features, and receive on output primary and secondary cosmic messengers including nuclei, neutrinos and photons. In extension to the propagation physics contained in a previous CRPropa version, the new version facilitates high-performance computing and comprises new physical features such as an interface for galactic propagation using lensing techniques, an improved photonuclear interaction calculation, and propagation in time dependent environments to take into account cosmic evolution effects in anisotropy studies and variable sources. First applications using highlighted features are presented as well.Keywords: ultra high energy cosmic rays, magnetic fields, ultra high energy photons and neutrinos 1 Corresponding author.

We report a study of the distributions of the depth of maximum, Xmax, of extensive air-shower profiles with energies above 10 17.8 eV as observed with the fluorescence telescopes of the Pierre Auger Observatory. The analysis method for selecting a data sample with minimal sampling bias is described in detail as well as the experimental cross-checks and systematic uncertainties. Furthermore, we discuss the detector acceptance and the resolution of the Xmax measurement and provide parameterizations thereof as a function of energy. The energy dependence of the mean and standard 4 deviation of the Xmax-distributions are compared to air-shower simulations for different nuclear primaries and interpreted in terms of the mean and variance of the logarithmic mass distribution at the top of the atmosphere.

A new analysis of the data set from the Pierre Auger Observatory provides evidence for anisotropy in the arrival directions of ultra-high-energy cosmic rays on an intermediate angular scale, which is indicative of excess arrivals from strong, nearby sources. The data consist of 5514 events above with zenith angles up to 80° recorded before 2017 April 30. Sky models have been created for two distinct populations of extragalactic gamma-ray emitters: active galactic nuclei from the second catalog of hard Fermi-LAT sources (2FHL) and starburst galaxies from a sample that was examined with Fermi-LAT. Flux-limited samples, which include all types of galaxies from the Swift-BAT and 2MASS surveys, have been investigated for comparison. The sky model of cosmic-ray density constructed using each catalog has two free parameters, the fraction of events correlating with astrophysical objects, and an angular scale characterizing the clustering of cosmic rays around extragalactic sources. A maximum-likelihood ratio test is used to evaluate the best values of these parameters and to quantify the strength of each model by contrast with isotropy. It is found that the starburst model fits the data better than the hypothesis of isotropy with a statistical significance of 4.0σ, the highest value of the test statistic being for energies above . The three alternative models are favored against isotropy with 2.7σ–3.2σ significance. The origin of the indicated deviation from isotropy is examined and prospects for more sensitive future studies are discussed.

A. AAB et al. we have examined the implications of the distributions of depths of atmospheric shower maximum (X max ), using a hybrid technique, for composition and hadronic interaction models. We do this by fitting the distributions with predictions from a variety of hadronic interaction models for variations in the composition of the primary cosmic rays and examining the quality of the fit. Regardless of what interaction model is assumed, we find that our data are not well described by a mix of protons and iron nuclei over most of the energy range. Acceptable fits can be obtained when intermediate masses are included, and when this is done consistent results for the proton and iron-nuclei contributions can be found using the available models. We observe a strong energy dependence of the resulting proton fractions, and find no support from any of the models for a significant contribution from iron nuclei. However, we also observe a significant disagreement between the models with respect to the relative contributions of the intermediate components.