• ### Limits on the Lorentz Invariance Violation from UHECR astrophysics(1701.04865)

Nov. 21, 2017 astro-ph.HE
In this paper, Lorentz Invariance Violation (LIV) is introduced in the calculations of photon propagation in the Universe. LIV is considered in the photon sector and the mean free path of the $\gamma \gamma \rightarrow e^{+} e^{-}$ interaction is calculated. The corresponding photon horizon including LIV effects is used to predict major changes in the propagation of photons with energy above $10^{18}$ eV. The flux of GZK photons on Earth considering LIV is calculated for several source models of ultra-high energy cosmic ray (UHECR). The predicted flux of GZK gamma-rays is compared to the new upper limits on the photon flux obtained by the Pierre Auger Observatory in order to impose upper limits on the LIV coefficients of order $n =$ 0, 1 and 2. The limits on the LIV coefficients derived here are more realistic than previous works and in some cases more restrictive. The analysis resulted in LIV upper limits in the photon sector of $\delta_{\gamma,0}^{limit} \sim -10^{-20}$, $\delta_{\gamma,1}^{limit} \sim -10^{-38} \; \mathrm{eV^{-1}}$ and $\delta_{\gamma,2}^{limit} \sim -10^{-56} \; \mathrm{eV^{-2}}$ in the astrophysical scenario which best describes UHECR data.
• ### Autonomous RPCs for a Cosmic Ray ground array(1709.09619)

Sept. 27, 2017 physics.ins-det, astro-ph.IM
We report on the behaviour of Resistive Plate Chambers (RPC) developed for muon detection in ultra-high energy cosmic ray (UHECR) experiments. The RPCs were developed for the MARTA project and were tested on field conditions. These RPCs cover an area of $1.5 \times 1.2\,{m^2}$ and are instrumented with 64 pickup electrodes providing a segmentation better than $20\,$cm. By shielding the detector units with enough slant mass to absorb the electromagnetic component in the air showers, a clean measurement of the muon content is allowed, a concept to be implemented in a next generation of UHECR experiments. The operation of a ground array detector poses challenging demands, as the RPC must operate remotely under extreme environmental conditions, with limited budgets for power and minimal maintenance. The RPC, DAQ, High Voltage and monitoring systems are enclosed in an aluminium-sealed case, providing a compact and robust unit suited for outdoor environments, which can be easily deployed and connected. The RPCs developed at LIP-Coimbra are able to operate using a very low gas flux, which allows running them for few years with a small gas reservoir. Several prototypes have already been built and tested both in the laboratory and outdoors. We report on the most recent tests done in the field that show that the developed RPCs have operated in a stable way for more than 2 years in field conditions.
• ### A Monte Carlo study of the relevance of fluorescence radiation in VHE gamma ray observations with Cherenkov telescopes(1708.09737)

Aug. 31, 2017 astro-ph.IM
It is generally assumed that fluorescence radiation does not play a significant role in the performance of Cherenkov telescopes. However, this assumption is required to be verified using detailed Monte Carlo simulations. In order to do this, we have implemented the production and tracking of fluorescence radiation inside the CORSIKA code, and simulated gamma-ray induced showers in the very high energy range. The most accurate fluorescence-yield data available so far was used for this purpose. The distribution of both light components on the ground has been studied as a function of various parameters affecting the detection and reconstruction of gamma-ray showers such as the angular aperture. From these distributions, we determined the conditions under which fluorescence radiation becomes significant. These results will also be useful to estimate the corresponding systematic errors in Cherenkov telescope observations. The full simulation results have been cross-checked, on a small sample of events, against numerical calculations based on a one-dimension shower profile and found to be compatible with each other. Both tools can be used for further investigations, like studying the possibility to modify Cherenkov telescopes for the measurement of fluorescence induced by extensive air showers.
• ### Astroparticle Physics Tests of Lorentz Invariance Violation(1708.00266)

Aug. 1, 2017 astro-ph.HE
Testing Lorentz invariance is essential as it is one of the pillars of modern physics. Moreover, its violation is foreseen in several popular Quantum Gravity models. Several authors study the effects of Lorentz invariance violation (LIV) in the propagation of ultra-high energy cosmic rays. These particles are the most energetic events ever detected and therefore represent a promising framework to test LIV. In this work we present an analytic calculation of the in-elasticity for any $a+b \rightarrow c+d$ interaction using first order perturbation in the dispersion relation that violates Lorentz invariance. The inelasticity can be calculated by solving a third-order polynomial equation containing: a) the kinematics of the interaction, b) the LIV term for each particle and c) the geometry of the interaction. We use the inelasticity we calculate to investigate the proton propagation in the intergalactic media. The photopion production of the proton interaction with the CMB is taken into account using the inelasticity and the attenuation length in different LIV scenarios. We show how the allowed phase space for the photopion production changes when LIV is considered for the interaction. The calculations presented here are going to be extended in order to calculated the modified ultra-high energy cosmic rays spectrum and compare it to the data.
• ### A new air-shower observable to constrain hadronic interaction models(1701.06764)

July 25, 2017 hep-ph, hep-ex, astro-ph.HE
The energy spectrum of muons at ground level in air showers are studied and a new observable is proposed to constrain hadronic interaction models used in air shower simulations. An asymmetric Gaussian function is proposed to describe the muon ground energy spectrum and its parameters are studied regarding primary particle, energy and hadronic interaction models. Based on two realistic measurements of the muon density at a given distance from the shower axis, a new observable (\rmu) is defined. Considering realistic values of detector resolutions and number of measured events, it is also shown \rmu can be successfully used to constrain low and high hadronic interaction models. The study is focused in the energy range between $10^{17.5}$ and $10^{18.0}$ eV because of the importance of this interval for particle physics and astrophysical models. The constraining power of the new observable is shown to be large within current experimental capabilities.
• ### The influence of the observatory latitude on the study of ultra high energy cosmic rays(1703.01529)

March 4, 2017 astro-ph.HE
Recent precision measurements of the Ultra High Energy Cosmic Rays (UHECR) arrival directions, spectrum and parameters related to the mass of the primary particle have been done by the HiRes, Pierre Auger and Telescope Array (TA) Observatories. In this paper, distributions of arrival directions of events in the nearby Universe are assumed to correlate with sources in the 2MASS Redshift Survey (2MRS), IRAS 1.2 Jy Survey, Palermo Swift-BAT and Swift-BAT catalogs, and the effect of the latitude of the observatory on the measurement of the energy spectrum and on the capability of measuring anisotropy is studied. The differences between given latitudes on the northern and southern hemispheres are quantified. It is shown that the latitude of the observatory: a) has an influence on the total flux measured and b) imposes an important limitation on the capability of measuring an anisotropic sky.
• ### Measurements of the depth of maximum of air-shower profiles at the Pierre Auger Observatory and their composition implications(1701.06812)

Jan. 24, 2017 astro-ph.HE
Air-showers measured by the Pierre Auger Observatory were analyzed in order to extract the depth of maximum (Xmax).The results allow the analysis of the Xmax distributions as a function of energy ($> 10^{17.8}$ eV). The Xmax distributions, their mean and standard deviation are analyzed with the help of shower simulations with the aim of interpreting the mass composition. The mean and standard deviation were used to derive <ln A> and its variance as a function of energy. The fraction of four components (p, He, N and Fe) were fit to the Xmax distributions. Regardless of the hadronic model used the data is better described by a mix of light, intermediate and heavy primaries. Also, independent of the hadronic models, a decrease of the proton flux with energy is observed. No significant contribution of iron nuclei is derived in the entire energy range studied.
• ### Interpretation of measurements of the number of muons in extensive air shower experiments(1606.08810)

June 28, 2016 astro-ph.HE
In this paper we analyze the energy evolution of the muon content of air showers between $10^{18.4}$ and $10^{19.6}$ eV to be able to determine the most likely mass composition scenario from future number of muons measurements. The energy and primary mass evolution of the number of muons is studied based on the Heitler-Matthews model and Monte Carlo simulation of the air shower. A simple model to describe the evolution of the first and second moments of number of muons distributions is proposed and validated. An analysis approach based on the comparison between this model's predictions and data to discriminate among a set of composition scenarios is presented and tested with simulations. It is shown that the composition scenarios can be potentially discriminated under the conditions imposed by the method. The discrimination power of the proposed analysis is stable under systematic changes of the absolute number of muons from model predictions and on the scale of the reconstructed energy.
• 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.
• ### Ultra high energy cosmic rays and possible signature of black strings(1510.04992)

March 15, 2016 hep-th, hep-ph, astro-ph.HE
Ultra high energy cosmic rays (UHECRs) probably originate in extreme conditions in which extra dimension effects might be important. In this paper we calculate the correction in black hole accretion mechanisms due to extra dimension effects in the static and rotating cases. A parametrization of the external Kerr horizons in both cases is presented and analysed. We use previous calculations of upper limits on the UHECR flux to set limits on the UHECR production efficiency of nine sources. The upper limit on the UHECR luminosity calculation is based on GeV-TeV gamma-ray measurements. The total luminosity due to the accretion mechanism is compared to the upper limit on UHECRs. The dependence of the UHECR production efficiency upper limit on black hole mass is also presented and discussed
• 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.
• ### Cosmic rays: the spectrum and chemical composition from $10^{10}$ to $10^{20}$ eV(1502.00305)

Feb. 1, 2015 astro-ph.HE
The production of energetic particles in the universe remains one of the great mysteries of modern science. The mechanisms of acceleration in astrophysical sources and the details about the propagation through the galactic and extragalactic media are still to be defined. In recent years, the cosmic ray flux has been measured with high precision in the energy range from \energy{10} to \energyEV{20.5} by several experiments using different techniques. In some energy ranges, it has been possible to determine the flux of individual elements (hydrogen to iron nuclei). This paper explores an astrophysical scenario in which only our Galaxy and the radio galaxy Cen A produce all particles measured on Earth in the energy range from \energy{10} to \energyEV{20.5}. Data from AMS-02, CREAM, KASCADE, KASCADE-Grande and the Pierre Auger Observatories are considered. The model developed here is able to describe the total and individual particle flux of all experiments considered. It is shown that the theory used here is able to describe the smooth transition from space-based to ground-based measurements. The flux of each element as determined by KASCADE and KASCADE-Grande and the mass sensitivity parameter \xmax measured by the Pierre Auger Observatory above \energyEV{18} are also explored within the framework of the model. The transition from \energy{16} to \energyEV{18} is carefully analyzed. It is shown that the data measured in this energy range suggest the existence of an extra component of cosmic rays yet to be understood.
• Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
• ### Comparison of the moments of the Xmax distribution predicted by different cosmic ray shower simulation models(1301.5555)

March 15, 2013 hep-ph, astro-ph.HE
In this paper we study the depth at which a cosmic ray shower reaches its maximum (\Xmax) as predicted by Monte Carlo simulation. The use of \Xmax in the determination of the primary particle mass can only be done by comparing the measured values with simulation predictions. For this reason it is important to study the differences between the available simulation models. We have done a study of the first and second moments of the \Xmax distribution using the \Corsika and \Conex programs. The study was done with high statistics in the energy range from $10^{17}$ to \energy{20.4}. We focus our analysis in the different implementations of the hadronic interaction models \Sibyll and \Qgsjet in \Corsika and \Conex. We show that the predictions of the \meanXmax and \sigmaXmax depend slightly on the combination of simulation program and hadronic interaction model. Although these differences are small, they are not negligible in some cases (up to 5 g/cm$^2$ for the worse case) and they should be considered as a systematic uncertainty of the model predictions for \meanXmax and \sigmaXmax. We have included a table with the suggested systematic uncertainties for the model predictions. Finally, we present a parametrization of the \Xmax distribution as a function of mass and energy according to the models \Sibyll and \Qgsjet, and showed an example of its application to obtain the predicted \Xmax distributions from cosmic ray propagation models.
• ### Centaurus A: the one extragalactic source of cosmic rays with energies above the knee(1106.0625)

Jan. 30, 2012 astro-ph.HE
The origin of cosmic rays at all energies is still uncertain. In this paper we present and explore an astrophysical scenario to produce cosmic rays with energy ranging from below $10^{15}$ o $3 \times 10^{20}$ eV. We show here that just our Galaxy and the radio galaxy Cen A, each with their own galactic cosmic ray particles, but with those from the radio galaxy pushed up in energy by a relativistic shock in the jet emanating from the active black hole, are sufficient to describe the most recent data in the energy range PeV to near ZeV. Data are available over this entire energy range from the experiments KASCADE, KASCADE-Grande and Pierre Auger Observatory. The energy spectrum calculated here correctly reproduces the measured spectrum beyond the knee, and contrary to widely held expectations, no other extragalactic source population is required to explain the data, even at energies far below the general cutoff expected at $6 \times 10^{19}$ eV, the Greisen-Zatsepin-Kuzmin turn-off due to interaction with the cosmological microwave background. We present several predictions for the source population, the cosmic ray composition and the propagation to Earth which can be tested in the near future.
• ### Ultra-High Energy Cosmic Rays from Centaurus A: Jet Interaction with Gaseous Shells(1006.5022)

June 25, 2010 hep-ex, astro-ph.CO, astro-ph.HE
Ultra high energy cosmic rays (UHECRs), with energies above ~6 x 10^19 eV, seem to show a weak correlation with the distribution of matter relatively near to us in the universe. It has earlier been proposed that UHECRs could be accelerated in either the nucleus or the outer lobes of the nearby radio galaxy Cen A. We show that UHECR production at a spatially intermediate location about 15 kpc northeast from the nucleus, where the jet emerging from the nucleus is observed to strike a large star-forming shell of gas, is a plausible alternative. A relativistic jet is capable of accelerating lower-energy heavy seed cosmic rays (CRs) to UHECRs on timescales comparable to the time it takes the jet to pierce the large gaseous cloud. In this model many CRs arising from a starburst, with a composition enhanced in heavy elements near the knee region around PeV, are boosted to ultra-high energies by the relativistic shock of a newly oriented jet. This model matches the overall spectrum shown by the Auger data and also makes a prediction for the chemical composition as a function of particle energy. We thus predict an observable anisotropy in the composition at high energy in the sense that lighter nuclei should preferentially be seen toward the general direction of Cen A. Taking into consideration the magnetic field models for the Galactic disk and a Galactic magnetic wind, this scenario may resolve the discrepancy between HiRes and Auger results concerning the chemical composition of UHECRs.
• ### Effects of the Energy Error Distribution of Fluorescence Telescopes on the UHECR energy spectrum(astro-ph/0702123)

Feb. 5, 2007 astro-ph, hep-ph, hep-ex
The measurement of the ultra high energy cosmic ray (UHECR) spectrum is strongly affected by uncertainties on the reconstructed energy. The determination of the presence or absence of the GZK cutoff and its position in the energy spectrum depends not only on high statistics but also on the shape of the energy error distribution. Here we determine the energy error distribution for fluorescence telescopes, based on a Monte Carlo simulation. The HiRes and Auger fluorescence telescopes are simulated in detail. We analyze the UHECR spectrum convolved with this energy error distribution. We compare this spectrum with one convolved with a lognormal error distribution as well as with a Gaussian error distribution. We show that the energy error distribution for fluorescence detectors can not be represented by these known distributions. We conclude that the convolved energy spectrum will be smeared but not enough to affect the GZK cutoff detection. This conclusion stands for both HiRes and Auger fluorescence telescopes. This result differs from the effect of the energy error distribution obtained with ground detectors and reinforces the importance of the fluorescence energy measurement. We also investigate the effect of possible fluorescence yield measurement errors in the energy spectrum.
• ### The shower size parameter as estimator of extensive air shower energy in fluorescence telescopes(astro-ph/0511677)

Nov. 23, 2005 astro-ph
The fluorescence technique has been successfully used to detect ultrahigh energy cosmic rays by indirect measurements. The underlying idea is that the number of charged particles in the atmospheric shower, i.e, its longitudinal profile, can be extracted from the amount of emitted nitrogen fluorescence light. However the influence of shower fluctuations and the very possible presence of different nuclear species in the primary cosmic ray spectrum makes the estimate of the shower energy from the fluorescence data analysis a difficult task. We investigate the potential of shower size at maximum depth as estimator of shower energy. The detection of the fluorescence light is simulated in detail and the reconstruction biases are carefully analyzed. We extend our calculations to both HiRes and EUSO experiments. This approach has shown some advantages to the reconstruction of the energy when compared to the standard analysis procedure.
• ### Acceptance of fluorescence detectors and its implication in energy spectrum inference at the highest energies(astro-ph/0509597)

Sept. 20, 2005 astro-ph
Along the years HiRes and AGASA experiments have explored the fluorescence and the ground array experimental techniques to measure extensive air showers, being both essential to investigate the ultra-high energy cosmic rays. However, such Collaborations have published contradictory energy spectra for energies above the GZK cut-off. In this article, we investigate the acceptance of fluorescence telescopes to different primary particles at the highest energies. Using CORSIKA and CONEX shower simulations without and with the new pre-showering scheme, which allows photons to interact in the Earth magnetic field, we estimate the aperture of the HiRes-I telescope for gammas, iron nuclei and protons primaries as a function of the number of simulated events and primary energy. We also investigate the possibility that systematic differences in shower development for hadrons and gammas could mask or distort vital features of the cosmic ray energy spectrum at energies above the photo-pion production threshold. The impact of these effects on the true acceptance of a fluorescence detector is analyzed in the context of top-down production models.
• ### Acceptance of fluorescence detectors for photons and its implication in energy spectrum inference at the highest energies(astro-ph/0509499)

Sept. 16, 2005 astro-ph
The HiRes and AGASA Collaborations have published contradictory energy spectra for energies above the GZK cut-off. In this article, we investigate the acceptance of fluorescence telescopes to different primary particles at the highest energies. Using CORSIKA shower simulations without and with the new pre-showering scheme, which allows photons to interact in the Earth magnetic field, we estimate the aperture of the HiRes I telescope for gammas and protons primaries. We calculate the dependence of the telescope sensitivity to primary particle identity. We also investigate the possibility that systematic differences in shower development for hadrons and gammas could mask or distort vital features of the cosmic ray energy spectrum at energies above the photo-pion production threshold. The impact of these effects on the true acceptance of a fluorescence detector is analyzed in the context of top-down production models.
• ### Alternative energy estimation from the shower lateral distribution function(astro-ph/0509501)

Sept. 16, 2005 astro-ph
The surface detector technique has been successfully used to detect cosmic ray showers for several decades. Scintillators or Cerenkov water tanks can be used to measure the number of particles and/or the energy density at a given depth in the atmosphere and reconstruct the primary particle properties. It has been shown that the experiment configuration and the resolution in reconstructing the core position determine a distance to the shower axis in which the lateral distribution function (LDF) of particles shows the least variation with respect to different primary particles type, simulation models and specific shapes of the LDF. Therefore, the signal at this distance (600 m for Haverah Park and 1000 m for Auger Observatory) has shown to be a good estimator of the shower energy. Revisiting the above technique, we show that a range of distances to the shower axis, instead of one single point, can be used as estimator of the shower energy. A comparison is done for the Auger Observatory configuration and the new estimator proposed here is shown to be a good and robust alternative to the standard single point procedure.
• ### Analysis of shower size as estimator of extensive air shower energy(astro-ph/0509502)

Sept. 16, 2005 astro-ph
The fluorescence technique has been successfully used to detect ultrahigh energy cosmic rays by indirect measurements. The underlying idea is that the number of charged particles in the atmospheric shower, i.e, its longitudinal profile, can be extracted from the amount of emitted nitrogen fluorescence light. However the influence of shower fluctuations and the very possible presence of different nuclear species in the primary cosmic ray spectrum make the estimate of the shower energy from the fluorescence data analysis a difficult task. We investigate the potential of shower size at maximum depth as estimator of shower energy. The detection of the fluorescence light is simulated in detail and the reconstruction biases are carefully analyzed. We extend our calculations to both Auger and EUSO experiments. This kind of approach is of particular interest for showers that are not fully contained inside the field of view of the detector.
• ### Analysis of extensive air showers with the hybrid code SENECA(astro-ph/0509503)

Sept. 16, 2005 astro-ph
The ultrahigh energy tail of the cosmic ray spectrum has been explored with unprecedented detail. For this reason, new experiments are exerting a severe pressure on extensive air shower modeling. Detailed fast codes are in need in order to extract and understand the richness of information now available. In this sense we explore the potential of SENECA, an efficient hybrid tridimensional simulation code, as a valid practical alternative to full Monte Carlo simulations of extensive air showers generated by ultrahigh energy cosmic rays. We discuss the influence of this approach on the main longitudinal characteristics of proton, iron nucleus and gamma induced air showers for different hadronic interaction models. We also show the comparisons of our predictions with those of CORSIKA code.