• ### Simulation of the Radiation Energy Release in Air Showers(1609.05743)

Sept. 19, 2016 astro-ph.HE
A simulation study of the energy released by extensive air showers in the form of MHz radiation is performed using the CoREAS simulation code. We develop an efficient method to extract this radiation energy from air-shower simulations. We determine the longitudinal profile of the radiation energy release and compare it to the longitudinal profile of the energy deposit by the electromagnetic component of the air shower. We find that the radiation energy corrected for the geometric dependence of the geomagnetic emission scales quadratically with the energy in the electromagnetic component of the air shower with a second order dependency on the atmospheric density at the position of the maximum of the shower development $X_\mathrm{max}$. In a measurement where $X_\mathrm{max}$ is not accessible, this second order dependence can be approximated using the zenith angle of the incoming direction of the air shower with only a minor deterioration in accuracy. This method results in an intrinsic uncertainty of 4% with respect to the electromagnetic shower energy which is well below current experimental uncertainties.
• ### Simulation of Radiation Energy Release in Air Showers(1606.01641)

Sept. 19, 2016 hep-ph, astro-ph.HE
A simulation study of the energy released by extensive air showers in the form of MHz radiation is performed using the CoREAS simulation code. We develop an efficient method to extract this radiation energy from air-shower simulations. We determine the longitudinal profile of the radiation energy release and compare it to the longitudinal profile of the energy deposit by the electromagnetic component of the air shower. We find that the radiation energy corrected for the geometric dependence of the geomagnetic emission scales quadratically with the energy in the electromagnetic component of the air shower with a second-order dependence on the atmospheric density at the position of the maximum shower development $X_\mathrm{max}$. In a measurement where $X_\mathrm{max}$ is not accessible, this second order dependence can be approximated using the zenith angle of the incoming direction of the air shower with only a minor loss in accuracy. Our method results in an intrinsic uncertainty of 4% in the determination of the energy in the electromagnetic air-shower component, which is well below current experimental uncertainties.
• 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.
Stochastic processes driven by stationary fractional Gaussian noise, that is, fractional Brownian motion and fractional Langevin equation motion, are usually considered to be ergodic in the sense that, after an algebraic relaxation, time and ensemble averages of physical observables coincide. Recently it was demonstrated that fractional Brownian motion and fractional Langevin equation motion under external confinement are transiently non-ergodic---time and ensemble averages behave differently---from the moment when the particle starts to sense the confinement. Here we show that these processes also exhibit transient ageing, that is, physical observables such as the time averaged mean squared displacement depend on the time lag between the initiation of the system at time $t=0$ and the start of the measurement at the ageing time $t_a$. In particular, it turns out that for fractional Langevin equation motion the ageing dependence on $t_a$ is different between the cases of free and confined motion. We obtain explicit analytical expressions for the aged moments of the particle position as well as the time averaged mean squared displacement and present a numerical analysis of this transient ageing phenomenon.