• 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.
• Neutrinos in the cosmic ray flux with energies near 1 EeV and above are detectable with the Surface Detector array of the Pierre Auger Observatory. We report here on searches through Auger data from 1 January 2004 until 20 June 2013. No neutrino candidates were found, yielding a limit to the diffuse flux of ultra-high energy neutrinos that challenges the Waxman-Bahcall bound predictions. Neutrino identification is attempted using the broad time-structure of the signals expected in the SD stations, and is efficiently done for neutrinos of all flavors interacting in the atmosphere at large zenith angles, as well as for "Earth-skimming" neutrino interactions in the case of tau neutrinos. In this paper the searches for downward-going neutrinos in the zenith angle bins $60^\circ-75^\circ$ and $75^\circ-90^\circ$ as well as for upward-going neutrinos, are combined to give a single limit. The $90\%$ C.L. single-flavor limit to the diffuse flux of ultra-high energy neutrinos with an $E^{-2}$ spectrum in the energy range $1.0 \times 10^{17}$ eV - $2.5 \times 10^{19}$ eV is $E_\nu^2 dN_\nu/dE_\nu < 6.4 \times 10^{-9}~ {\rm GeV~ cm^{-2}~ s^{-1}~ sr^{-1}}$.
• 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.
• ### Exploiting the geomagnetic distortion of the inclined atmospheric showers(1508.04354)

Aug. 18, 2015 astro-ph.IM, astro-ph.HE
We propose a novel approach for the determination of the nature of ultra-high energy cosmic rays by exploiting the geomagnetic deviation of muons in nearly horizontal showers. The distribution of the muons at ground level is well described by a simple parametrization providing a few shape parameters tightly correlated to $X^\mu_\mathrm{max}$, the depth of maximal muon production, which is a mass indicator tightly correlated to the usual parameter $X_\mathrm{max}$, the depth of maximal development of the shower. We show that some constraints can be set on the predictions of hadronic models, especially by combining the geomagnetic distortion with standard measurement of the longitudinal profile. We discuss the precision needed to obtain significant results and we propose a schematic layout of a detector.
• ### Layered water Cherenkov detector for the study of ultra high energy cosmic rays(1405.5699)

May 22, 2014 astro-ph.IM
We present a new design for the water Cherenkov detectors that are in use in various cosmic ray observatories. This novel design can provide a significant improvement in the independent measurement of the muonic and electromagnetic component of extensive air showers. From such multi-component data an event by event classification of the primary cosmic ray mass becomes possible. According to popular hadronic interaction models, such as EPOS-LHC or QGSJetII-04, the discriminating power between iron and hydrogen primaries reaches Fisher values of $\sim$ 2 or above for energies in excess of $10^{19}$ eV with a detector array layout similar to that of the Pierre Auger Observatory.
• Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
• ### A phenomenological model of the muon density profile on the ground of very inclined air showers(0904.2372)

June 28, 2010 astro-ph.IM
Ultra-high energy cosmic rays generate extensive air showers in Earth's atmosphere. A standard approach to reconstruct the energy of an ultra-high energy cosmic rays is to sample the lateral profile of the particle density on the ground of the air shower with an array of surface detectors. For cosmic rays with large inclinations, this reconstruction is based on a model of the lateral profile of the muon density observed on the ground, which is fitted to the observed muon densities in individual surface detectors. The best models for this task are derived from detailed Monte-Carlo simulations of the air shower development. We present a phenomenological parametrization scheme which allows to derive a model of the average lateral profile of the muon density directly from a fit to a set of individual Monte-Carlo simulated air showers. The model reproduces the detailed simulations with a high precision. As an example, we generate a muon density model which is valid in the energy range 1e18 eV < E < 1e20 eV and the zenith angle range 60 deg < theta < 90 deg. We will further demonstrate a way to speed up the simulation of such muon profiles by three orders of magnitude, if only the muons in the shower are of interest.
• ### Estimates of multipolar coefficients to search for cosmic ray anisotropies with non-uniform or partial sky coverage(0710.2290)

Feb. 13, 2008 astro-ph
We study the possibility to extract the multipolar moments of an underlying distribution from a set of cosmic rays observed with non-uniform or even partial sky coverage. We show that if the degree is assumed to be upper bounded by $L$, each multipolar moment can be recovered whatever the coverage, but with a variance increasing exponentially with the bound $L$ if the coverage is zero somewhere. Despite this limitation, we show the possibility to test predictions of a model without any assumption on $L$ by building an estimate of the covariance matrix seen through the exposure function.
• ### A coverage independent method to analyze large scale anisotropies(0706.3705)

June 25, 2007 astro-ph
The arrival time distribution of cosmic ray events is well suited to extract information regarding sky anisotropies. For an experiment with nearly constant exposure, the frequency resolution one can achieve is given by the inverse of the time $T$ during which the data was recorded. For $T$ larger than one calendar year the resolution becomes sufficient to resolve the sidereal and diurnal frequencies. Using a Fourier expansion on a modified time parameter, we show in this note that one can accurately extract sidereal modulations without knowledge of the experimental coverage. This procedure also gives the full frequency pattern of the event sample under studies which contains important information about possible systematics entering in the sidereal analysis. We also show how this method allows to correct for those systematics. Finally, we show that a two dimensional analysis, in the form of the spherical harmonic ($Y_l^m$) decomposition, can be performed under the same conditions for all $m\ne 0$.
• ### Evaluation of the Primary Energy of UHE Photon-induced Atmospheric Showers from Ground Array Measurements(astro-ph/0701583)

Jan. 22, 2007 astro-ph
A photon induced shower at $E_{prim}\ge 10^{18}$ eV exhibits very specific features and is different from a hadronic one. At such energies, the LPM effect delays in average the first interactions of the photon in the atmosphere and hence slows down the whole shower development. They also have a smaller muonic content than hadronic ones. The response of a surface detector such as that of the Auger Observatory to these specific showers is thus different and has to be accounted for in order to enable potential photon candidates reconstruction correctly. The energy reconstruction in particular has to be adapted to the late development of photon showers. We propose in this article a method for the reconstruction of the energy of photon showers with a surface detector. The key feature of this method is to rely explicitly on the development stage of the shower. This approach leads to very satisfactory results ($\simeq 20%$). At even higher energies ($5.10^{19}$ eV and above) the probability for the photon to convert into a pair of e$^+$e$^-$ in the geomagnetic field becomes non negligible and requires a different function to evaluate the energy with the proposed method. We propose several approaches to deal with this issue in the scope of the establishment of an upper bound on the photon fraction in UHECR.
• ### Note on the Origin of the Highest Energy Cosmic Rays(astro-ph/0001427)

Jan. 25, 2000 astro-ph
In this note we argue that the galactic model chosen by E.-J. Ahn, G. Medina-Tanco, P.L. Bierman and T. Stanev in their paper discussing the origin of the highest energy cosmic rays, is alone responsible for the focussing of positive particles towards the North galactic pole. We discuss the validity of this model, in particular in terms of field reversals and radial extensions. We conclude that with such a model one cannot retreive any directional information from the observed direction of the cosmic rays. In particular one cannot identify point sources at least up to energies of about 200 EeV. Therefore the apparent clustering of the back-traced highest energy cosmic rays observed to date cannot be interpreted as an evidence for a point source nor for the identification of M87, which happens to be close to the North pole, as being such a source.