Assuming that the short gamma-ray burst detected by the Fermi Gamma-Ray Space Telescope about 0.4 seconds after the gravitational waves observed by the LIGO and VIRGO Collaborations originated from the same black hole merger event, we perform a model-independent analysis of different quantum gravity scenarios based on (modified) dispersion relations (typical of quantum gravity models) for the graviton and the photon. We find that only scenarios where at least one of the two particles is luminal (the other being sub- or super-luminal) are allowed, while scenarios where none of the two particles is luminal are ruled out. Moreover, the physical request of having acceptable values for the quantum gravity scale imposes stringent bounds on the difference between the velocities of electromagnetic and gravitational waves, much more stringent than any previously known bound.

Spherical harmonic moments are well-suited for capturing anisotropy at any scale in the flux of cosmic rays. An unambiguous measurement of the full set of spherical harmonic coefficients requires full-sky coverage. This can be achieved by combining data from observatories located in both the northern and southern hemispheres. To this end, a joint analysis using data recorded at the Telescope Array and the Pierre Auger Observatory above $10^{19}$ eV is presented in this work. The resulting multipolar expansion of the flux of cosmic rays allows us to perform a series of anisotropy searches, and in particular to report on the angular power spectrum of cosmic rays above $10^{19}$ eV. No significant deviation from isotropic expectations is found throughout the analyses performed. Upper limits on the amplitudes of the dipole and quadrupole moments are derived as a function of the direction in the sky, varying between 7% and 13% for the dipole and between 7% and 10% for a symmetric quadrupole.

Ultra-high energy (UHE) photons play an important role as an independent probe of the photo-pion production mechanism by UHE cosmic rays. Their observation, or non-observation, may constrain astrophysical scenarios for the origin of UHECRs and help to understand the nature of the flux suppression observed by several experiments at energies above $10^{19.5}$ eV. Whereas the interaction length of UHE photons above $10^{17}$ eV ranges from a few hundred kpc up to tenths of Mpc, photons can interact with the extragalactic background radiation initiating the development of electromagnetic cascades which affect the fluxes of photons observed at Earth. The interpretation of the current experimental results rely on the simulations of the UHE photon propagation. In this paper, we present the novel Monte Carlo code EleCa to simulate the $Ele$ctromagnetic $Ca$scading initiated by high-energy photons and electrons. We provide an estimation of the surviving probability for photons inducing electromagnetic cascades as a function of their distance from the observer and we calculate the distances within which we expect to observe UHE photons with energy between $10^{17}$ and $10^{19}$ eV. Furthermore, the flux of GZK photons at Earth is investigated in several astrophysical scenarios where we vary both injection spectrum and composition at the source and the intensity of the intervening extragalactic magnetic field. Although the photon propagation depends on several astrophysical factors, our numerical predictions combined with future experimental observations (or non-observations) of UHE photons -- in the energy range between $10^{17.5}$ eV and $10^{20}$ eV -- can help to constrain these scenarios.

We describe the method devised to reconstruct inclined cosmic-ray air showers with zenith angles greater than $60^\circ$ detected with the surface array of the Pierre Auger Observatory. The measured signals at the ground level are fitted to muon density distributions predicted with atmospheric cascade models to obtain the relative shower size as an overall normalization parameter. The method is evaluated using simulated showers to test its performance. The energy of the cosmic rays is calibrated using a sub-sample of events reconstructed with both the fluorescence and surface array techniques. The reconstruction method described here provides the basis of complementary analyses including an independent measurement of the energy spectrum of ultra-high energy cosmic rays using very inclined events collected by the Pierre Auger Observatory.

A flux of neutrons from an astrophysical source in the Galaxy can be detected in the Pierre Auger Observatory as an excess of cosmic-ray air showers arriving from the direction of the source. To avoid the statistical penalty for making many trials, classes of objects are tested in combinations as nine "target sets", in addition to the search for a neutron flux from the Galactic Center or from the Galactic Plane. Within a target set, each candidate source is weighted in proportion to its electromagnetic flux, its exposure to the Auger Observatory, and its flux attenuation factor due to neutron decay. These searches do not find evidence for a neutron flux from any class of candidate sources. Tabulated results give the combined p-value for each class, with and without the weights, and also the flux upper limit for the most significant candidate source within each class. These limits on fluxes of neutrons significantly constrain models of EeV proton emission from non-transient discrete sources in the Galaxy.

Measurements of air showers made using the hybrid technique developed with the fluorescence and surface detectors of the Pierre Auger Observatory allow a sensitive search for point sources of EeV photons anywhere in the exposed sky. A multivariate analysis reduces the background of hadronic cosmic rays. The search is sensitive to a declination band from -85{\deg} to +20{\deg}, in an energy range from 10^17.3 eV to 10^18.5 eV. No photon point source has been detected. An upper limit on the photon flux has been derived for every direction. The mean value of the energy flux limit that results from this, assuming a photon spectral index of -2, is 0.06 eV cm^-2 s^-1, and no celestial direction exceeds 0.25 eV cm^-2 s^-1. These upper limits constrain scenarios in which EeV cosmic ray protons are emitted by non-transient sources in the Galaxy.

The Pierre Auger Observatory is the world's largest cosmic ray observatory. Our current exposure reaches nearly 40,000 km$^2$ str and provides us with an unprecedented quality data set. The performance and stability of the detectors and their enhancements are described. Data analyses have led to a number of major breakthroughs. Among these we discuss the energy spectrum and the searches for large-scale anisotropies. We present analyses of our X$_{max}$ data and show how it can be interpreted in terms of mass composition. We also describe some new analyses that extract mass sensitive parameters from the 100% duty cycle SD data. A coherent interpretation of all these recent results opens new directions. The consequences regarding the cosmic ray composition and the properties of UHECR sources are briefly discussed.

Joint contributions of the Pierre Auger and Telescope Array Collaborations to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013: cross-calibration of the fluorescence telescopes, large scale anisotropies and mass composition.

Ultra-high energy cosmic rays (UHECRs) interacting with the atmosphere generate extensive air showers (EAS) of secondary particles. The depth corresponding to the maximum development of the shower, $\Xmax$, is a well-known observable for determining the nature of the primary cosmic ray which initiated the cascade process. In this paper, we present an empirical model to describe the distribution of $\Xmax$ for EAS initiated by nuclei, in the energy range from $10^{17}$ eV up to $10^{21}$ eV, and by photons, in the energy range from $10^{17}$ eV up to $10^{19.6}$ eV. Our model adopts the generalized Gumbel distribution motivated by the relationship between the generalized Gumbel statistics and the distribution of the sum of non-identically distributed variables in dissipative stochastic systems. We provide an analytical expression for describing the $\Xmax$ distribution for photons and for nuclei, and for their first two statistical moments, namely $\langle \Xmax\rangle$ and $\sigma^{2}(\Xmax)$. The impact of the hadronic interaction model is investigated in detail, even in the case of the most up-to-date models accounting for LHC observations. We also briefly discuss the differences with a more classical approach and an application to the experimental data based on information theory.

Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013

Although the recent experimental efforts to improve the observation of Ultra-High Energy Cosmic Rays (UHECRs) above $10^{18}$ eV, the origin and the composition of such particles is still unknown. In this work, we present the novel Monte Carlo code (HERMES) simulating the propagation of UHE nuclei, in the energy range between $10^{16}$ and $10^{22}$ eV, accounting for propagation in the intervening extragalactic and Galactic magnetic fields and nuclear interactions with relic photons of the extragalactic background radiation. In order to show the potential applications of HERMES for astroparticle studies, we estimate the expected flux of UHE nuclei in different astrophysical scenarios, the GZK horizons and we show the expected arrival direction distributions in the presence of turbulent extragalactic magnetic fields. A stable version of HERMES will be released in the next future for public use together with libraries of already propagated nuclei to allow the community to perform mass composition and energy spectrum analysis with our simulator.

Ultra high energy photons play an important role as an independent probe of the photo-pion production mechanism by UHE cosmic rays. Their observation, or non-observation, may constrain astrophysical scenarios for the origin of UHECRs and help to understand the nature of the flux suppression observed by several experiments at energies above 10$^{19.5}$ eV. Whereas the interaction length of UHE photons above 10$^{17}$ eV is only of a few hundred kpc up to tenths of Mpc, photons can interact with the extragalactic background radiation leading to the development of electromagnetic cascades which affect the fluxes of photons observed at Earth. The interpretation of the current experimental results rely on the simulations of the UHE photon propagation. In this contribution, we present the novel Monte Carlo code "EleCa" to simulate the \emph{Ele}ctromagnetic \emph{Ca}scading initiated by high-energy photons and electrons. The distance within which we expect to observe UHE photons is discussed and the flux of GZK photons at Earth is investigated in several astrophysical scenarios.

The study of ultra-high energy cosmic rays (UHECR) at Earth cannot prescind from the study of their propagation in the Universe. In this paper, we present HERMES, the \emph{ad hoc} Monte Carlo code we have developed for the realistic simulation of UHECR propagation. We discuss the modeling adopted to simulate the cosmology, the magnetic fields, the interactions with relic photons and the production of secondary particles. In order to show the potential applications of HERMES for astroparticle studies, we provide an estimation of the surviving probability of UHE protons, the GZK horizons of nuclei and the all-particle spectrum observed at Earth in different astrophysical scenarios. Finally, we show the expected arrival direction distribution of UHECR produced from nearby candidate sources. A stable version of HERMES will be released in the next future for public use together with libraries of already propagated nuclei to allow the community to perform mass composition and energy spectrum analysis with our simulator.

We investigate how the density of baryonic and cold dark matter, the density of dark energy and the value of the Hubble parameter at the present time influence the propagation of ultrahigh energy protons in the nearby Universe. We take into account energy losses in the cosmic microwave radiation, the only one relevant for protons above $10^{18}$ eV, and we explore the dependence of Greisen-Zatsepin-Kuz'min (GZK) horizon on the cosmology. We investigate several cosmological scenarios, from matter dominated to energy dominated ones, and we consider the impact of uncertainties in the Hubble parameter in a $\Lambda-$Cold Dark Matter (CDM) Universe, estimated from recent observations, on the GZK horizon. The impact of the (unknown) extragalactic magnetic field on our study is discussed, as well as possible probes of the Hubble parameter attainable by current and future experiments.

The Pierre Auger Observatory is exploring the potential of the radio detection technique to study extensive air showers induced by ultra-high energy cosmic rays. The Auger Engineering Radio Array (AERA) addresses both technological and scientific aspects of the radio technique. A first phase of AERA has been operating since September 2010 with detector stations observing radio signals at frequencies between 30 and 80 MHz. In this paper we present comparative studies to identify and optimize the antenna design for the final configuration of AERA consisting of 160 individual radio detector stations. The transient nature of the air shower signal requires a detailed description of the antenna sensor. As the ultra-wideband reception of pulses is not widely discussed in antenna literature, we review the relevant antenna characteristics and enhance theoretical considerations towards the impulse response of antennas including polarization effects and multiple signal reflections. On the basis of the vector effective length we study the transient response characteristics of three candidate antennas in the time domain. Observing the variation of the continuous galactic background intensity we rank the antennas with respect to the noise level added to the galactic signal.

The Pierre Auger Observatory is a facility built to detect air showers produced by cosmic rays above 10^17 eV. During clear nights with a low illuminated moon fraction, the UV fluorescence light produced by air showers is recorded by optical telescopes at the Observatory. To correct the observations for variations in atmospheric conditions, atmospheric monitoring is performed at regular intervals ranging from several minutes (for cloud identification) to several hours (for aerosol conditions) to several days (for vertical profiles of temperature, pressure, and humidity). In 2009, the monitoring program was upgraded to allow for additional targeted measurements of atmospheric conditions shortly after the detection of air showers of special interest, e.g., showers produced by very high-energy cosmic rays or showers with atypical longitudinal profiles. The former events are of particular importance for the determination of the energy scale of the Observatory, and the latter are characteristic of unusual air shower physics or exotic primary particle types. The purpose of targeted (or "rapid") monitoring is to improve the resolution of the atmospheric measurements for such events. In this paper, we report on the implementation of the rapid monitoring program and its current status. The rapid monitoring data have been analyzed and applied to the reconstruction of air showers of high interest, and indicate that the air fluorescence measurements affected by clouds and aerosols are effectively corrected using measurements from the regular atmospheric monitoring program. We find that the rapid monitoring program has potential for supporting dedicated physics analyses beyond the standard event reconstruction.

We present a novel method to generate a synthetic distribution of objects (mock) on a spherical surface (i.e. a sky), by using a real distribution. The resulting surrogate map mimics the clustering features of the real data, including the effects of non-uniform exposure, if any. The method is model-independent, also preserving the angular correlation function, as well as the angular power spectrum, of the original data. It can be reliably adopted to mimic the angular clustering of objects distributed on a spherical surface and it can be easily extended to include further information, as the spatial clustering of objects distributed inside a sphere. Applications to real data are presented and discussed. In particular, we consider the distribution of galaxies recently presented in the 2MASS Redshift Survey (2MRS) catalog.

The Surface Detector of the Pierre Auger Observatory is sensitive to neutrinos of all flavours above 0.1 EeV. These interact through charged and neutral currents in the atmosphere giving rise to extensive air showers. When interacting deeply in the atmosphere at nearly horizontal incidence, neutrinos can be distinguished from regular hadronic cosmic rays by the broad time structure of their shower signals in the water-Cherenkov detectors. In this paper we present for the first time an analysis based on down-going neutrinos. We describe the search procedure, the possible sources of background, the method to compute the exposure and the associated systematic uncertainties. No candidate neutrinos have been found in data collected from 1 January 2004 to 31 May 2010. Assuming an E^-2 differential energy spectrum the limit on the single flavour neutrino is (E^2 * dN/dE) < 1.74x10^-7 GeV cm^-2 s^-1 sr^-1 at 90% C.L. in the energy range 1x10^17 eV < E < 1x10^20 eV.

Atmospheric conditions at the site of a cosmic ray observatory must be known for reconstructing observed extensive air showers. The Global Data Assimilation System (GDAS) is a global atmospheric model predicated on meteorological measurements and numerical weather predictions. GDAS provides altitude-dependent profiles of the main state variables of the atmosphere like temperature, pressure, and humidity. The original data and their application to the air shower reconstruction of the Pierre Auger Observatory are described. By comparisons with radiosonde and weather station measurements obtained on-site in Malarg\"ue and averaged monthly models, the utility of the GDAS data is shown.

We present a comprehensive study of the influence of the geomagnetic field on the energy estimation of extensive air showers with a zenith angle smaller than $60^\circ$, detected at the Pierre Auger Observatory. The geomagnetic field induces an azimuthal modulation of the estimated energy of cosmic rays up to the ~2% level at large zenith angles. We present a method to account for this modulation of the reconstructed energy. We analyse the effect of the modulation on large scale anisotropy searches in the arrival direction distributions of cosmic rays. At a given energy, the geomagnetic effect is shown to induce a pseudo-dipolar pattern at the percent level in the declination distribution that needs to be accounted for.

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an extensive air shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10^{17} and 10^{19} eV and zenith angles up to 65 degs. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte-Carlo results showing how LTP functions from data are in good agreement with simulations.