Search

Atmosphere-skimming air showers are initiated by cosmic rays with incoming directions that allow the cascade to develop entirely within the atmosphere, without reaching the ground. Radio pulses induced by this type of showers have already been observed in balloon-borne experiments such as ANITA, but a detailed characterisation of their properties is lacking. The extreme range of densities in which these cascades can develop gives rise to a wide range of shower profiles, with radio emission characteristics that can differ significantly from those of regular downward-going showers. In this work, we have used the ZHAireS-RASPASS program to characterise the expected radio emission from atmosphere-skimming air showers and its properties. We have studied the interplay between the magnetic field and atmospheric density profile in the expected radio signal, focusing on its detection aboard balloon-borne experiments. The almost horizontal geometry of the events gives rise to a significant \textit{refractive} asymmetry in the spatial distribution of the electric field, due to the propagation of the radio signals across a gradient of index of refraction. In addition, a unique \textit{coherence} asymmetry appears in the intensity of the signals, as a consequence of the cumulative effect of the Earth's magnetic field over the very long distances that these particle cascades traverse. The implications of the peculiar characteristics of the emission are discussed regarding their impact both on the interpretation of collected data and in the exposure of balloon-borne experiments, Comment: 22 pages, 13 figures, to be submitted to JCAP

Atmosphere-skimming showers are initiated by cosmic rays with incoming directions such that the full development of the cascade occurs inside the atmosphere without ever reaching the ground. This new class of showers has been observed in balloon-borne experiments such as ANITA, but a characterisation of their properties is lacking. The interplay between the Earth's magnetic field, the long distances over which atmosphere-skimming showers develop, and the low density of the atmosphere they traverse gives rise to several effects that are not seen in downward-going cascades, and require detailed modeling. In this article, we used the latest version of the ZHAireS-RASPASS shower simulation program to tackle this problem, and dwell on the particular phenomena that arises from the peculiar environment on which these showers develop. We focus in particular on the properties of the longitudinal profile of the shower and its fluctuations as a function of cosmic-ray energy, direction and primary mass. We have also studied the phase-space of cosmic-ray arrival directions where detection in high-altitude experiments is more likely, and have found that only in a small range of directions the showers are sufficiently developed before reaching the altitude of the detector. Our results are relevant for the design of high-altitude and in particular balloon-borne experiments, and for the interpretation of the data they collect., Comment: 18 pages, 11 figures. Submited to JCAP

The Antarctic Impulsive Transient Antenna (ANITA) detector has observed several radio pulses coming from the surface of the ice cap at the South Pole. These pulses were attributed to upward-going atmospheric particle showers instead of the downward-going showers induced by cosmic rays that exhibit a characteristic polarity inversion of the radio signal due to reflection in the ice. Coherent transition radiation from cosmic-ray showers developing in the atmosphere and intercepting the ice surface has been suggested as a possible and alternative explanation of these so-called "anomalous" events. To test this interpretation, we have developed an extension of ZHS, a program to calculate coherent pulses from electromagnetic showers, to deal with showers that transit a planar interface between two homogeneous and dielectric media, including transition radiation. By considering different geometries, it is found that all pulses from air showers intercepting the ice surface and detected at the height of ANITA, display the same polarity as pulses emitted by ultra-high-energy cosmic-ray showers that fully develop in the atmosphere and are reflected on the ice. We find that transition radiation is disfavored as a possible explanation of the anomalous ANITA events.

Charged current interactions of neutrinos inside the Earth can result in secondary muons and $\tau$-leptons which are detectable by several existing and planned neutrino experiments through a wide variety of event topologies. Consideration of such events can improve detector performance and provide unique signatures which help with event reconstruction. In this work, we describe NuLeptonSim, a propagation tool for neutrinos and charged leptons that builds on the fast NuTauSim framework. NuLeptonSim considers energy losses of charged leptons, modelled both continuously for performance or stochastically for accuracy, as well as interaction models for all flavors of neutrinos, including the Glashow resonance. We demonstrate the results from including these effects on the Earth emergence probability of various charged leptons from different flavors of primary neutrino and their corresponding energy distributions. We find that the emergence probability of muons can be higher than that of taus for energies below 100 PeV, whether from a primary muon or $\tau$ neutrino, and that the Glashow resonance contributes to a surplus of emerging leptons near the resonant energy.

We have studied the effect of changing the density and magnetic field strength in the coherent pulses that are emitted as energetic showers develop in the atmosphere. For this purpose we have developed an extension of ZHS, a program to calculate coherent radio pulses from electromagnetic showers in homogeneous media, to account for the Lorentz force due to a magnetic field. This makes it possible to perform quite realistic simulations of radio pulses from air showers in a medium similar to the atmosphere but without variations of density with altitude. The effects of independently changing the density, the refractive index and the magnetic field strength are studied in the frequency domain for observers in the Cherenkov direction at far distances from the shower. This approach is particularly enlightening providing an explanation of the spectral behavior of the induced electric field in terms of shower development parameters. More importantly, it clearly displays the complex scaling properties of the pulses as density and magnetic field intensity are varied. The usually assumed linear behavior of electric field amplitude with magnetic field intensity is shown to hold up to a given magnetic field strength at which the extra time delays due to the deflection in the magnetic field break it. Scaling properties of the pulses are obtained as the density of air decreases relative to sea level. A remarkably accurate scaling law is obtained that relates the spectra of pulses obtained when reducing the density and increasing the magnetic field., Comment: 27 pages, 14 figures

We assess the capabilities of a ground-based gamma-ray observatory to detect astrophysical neutrinos with energies in the $100\,{\rm TeV}$ to $100\,{\rm PeV}$ range. The identification of these events would be done through the measurement of very inclined extensive air showers induced by downward-going and upward-going neutrinos. The discrimination of neutrino-induced showers in the overwhelming cosmic-ray background is achieved by analysing the balance of the total electromagnetic and muonic signals of the shower at the ground. We demonstrate that a ${\rm km^2}$-scale wide field-of-view ground-based gamma-ray observatory could detect a couple of Very-High to Ultra-High energy (VHE-UHE) neutrino events per year with a reasonable pointing accuracy, making it an interesting facility for multi-messenger studies with both photons and neutrinos., Comment: 12 pages, 14 figures, accepted in Phys.Rev.D

Astrophysical neutrinos are excellent probes of astroparticle physics and high-energy physics. With energies far beyond solar, supernovae, atmospheric, and accelerator neutrinos, high-energy and ultra-high-energy neutrinos probe fundamental physics from the TeV scale to the EeV scale and beyond. They are sensitive to physics both within and beyond the Standard Model through their production mechanisms and in their propagation over cosmological distances. They carry unique information about their extreme non-thermal sources by giving insight into regions that are opaque to electromagnetic radiation. This white paper describes the opportunities astrophysical neutrino observations offer for astrophysics and high-energy physics, today and in coming years., Comment: Contribution to Snowmass 2021, updated to include community feedback

Tau neutrinos are the least studied particle in the Standard Model. This whitepaper discusses the current and expected upcoming status of tau neutrino physics with attention to the broad experimental and theoretical landscape spanning long-baseline, beam-dump, collider, and astrophysical experiments. This whitepaper was prepared as a part of the NuTau2021 Workshop., Comment: 155 pages, 63 figures, 10 tables, 69 endorsers, comments welcome, Contribution to Snowmass 2021. v2: published version

Compilation of the six contributions to the ICRC conference 2021 by the CORSIKA 8 Collaboration. The status of the project is illustrated. In particular, the secondary hadron as well as the electromagnetic cascades are being validated individually, and current results are reviewed. A novel framework for radio emission simulations is presented, which is designed given the modular nature of CORSIKA 8 to support, both, the CoREAS as well as the ZHS formalism. At the same time, first Cherenkov emission calculations are shown which are based on CORSIKA 8 coupled with a GPU Cherenkov emission code. Finally, a new powerful feature of CORSIKA 8 is illustrated, where the entire genealogy of air shower particles can be studied in all details., Comment: ICRC 2021 conference, 6 proceedings, 51 pages

Ultra-high energy cosmic rays (UHECRs) are the most energetic particles observed and serve as a probe of the extreme universe. A key question to understanding the violent processes responsible for their acceleration is identifying which classes of astrophysical objects (active galactic nuclei or starburst galaxies, for example) correlate to their arrival directions. While source clustering is limited by deflections in the Galactic magnetic field, at the highest energies the scattering angles are sufficiently low to retain correlation with source catalogues. While there have been several studies attempting to identify source catalogue correlations with data from the Pierre Auger Observatory and the Telescope Array, the significance above an isotropic background has not yet reached the threshold for discovery. It has been known for several decades that a full-sky UHECR observatory would provide a substantial increase in sensitivity to the anisotropic component of UHECRs. There have been several concepts developed in that time targeting the identification of UHECR sources such as OWL, JEM-EUSO, and POEMMA, using fluorescence detection in the Earth's atmosphere from orbit. In this white paper, we present a concept called the Zettavolt Askaryan Polarimeter (ZAP), designed to identify the source of UHECRs using radio detection of the Askaryan radio emissions produced by UHECRs interacting in the Moon's regolith from low lunar orbit., Comment: 6 pages, 2 figures, letter of interest submitted to Snowmass 2021

Tau neutrinos are expected to comprise roughly one third of both the astrophysical and cosmogenic neutrino flux, but currently the flavor ratio is poorly constrained and the expected flux at energies above $10^{17}$ eV is low. We present a detector concept aimed at measuring the diffuse flux of tau neutrinos in this energy range via a high-elevation mountaintop detector using the radio technique. The detector searches for radio signals from upgoing air showers generated by Earth-skimming tau neutrinos. Signals from several antennas in a compact array are coherently summed at the trigger level, permitting not only directional masking of anthropogenic backgrounds, but also a low trigger threshold. This design takes advantage of both the large viewing area available at high-elevation sites and the nearly full duty cycle available to radio instruments. We present trade studies that consider the station elevation, frequency band, number of antennas in the array, and the trigger threshold to develop a highly efficient station design. Such a mountaintop detector can achieve a factor of ten improvement in acceptance over existing instruments with 100 independent stations. With 1000 stations and three years of observation, it can achieve a sensitivity to an integrated $\mathcal{E}^{-2}$ flux of $<10^{-9}$ GeV cm$^{-2}$ sr$^{-1}$ s$^{-1}$, in the range of the expected flux of all-flavor cosmogenic neutrinos assuming a pure iron cosmic-ray composition., Comment: 26 pages, 11 figures

We present a semi-analytical method for the calculation of coherent Askaryan radiation in showers induced by neutrinos of any flavor in ice. We compare our results with those of a full Monte Carlo simulation based on the ZHAireS code. This approach is able to reproduce the vector potential and hence electric field at any experimentally relevant observer position in the time domain. This work extends published results only valid for electron-induced showers. We establish the validity of the semi-analytical calculation of the radio signal produced by all types of neutrino-induced showers in ice. The method is computationally efficient and only requires as inputs the longitudinal charge excess profile of the showers and a parameterization of the vector potential in the far-field region of the shower at the Cherenkov angle that we also provide. Our methodology avoids tracking the contributions to the electric field from millions of particles every time the radio pulse has to be calculated at a given observer position. These results can be readily used in the interpretation of the data taken by experiments, and in the planning and design of future initiatives based on the radio technique in ice., Comment: 18 pages, 11 figures. Accepted in Phys. Rev. D

High-energy astrophysical neutrinos, recently discovered by IceCube up to energies of several PeV, opened a new window to the high-energy Universe. Yet much remains to be known. IceCube has excellent muon flavor identification, but tau flavor identification is challenging. This limits its ability to probe neutrino physics and astrophysics. To address this limitation, we present a concept for a large-scale observatory of astrophysical tau neutrinos in the 1-100 PeV range, where a flux is guaranteed to exist. Its detection would allow us to characterize the neutrino sources observed by IceCube, to discover new ones, and test neutrino physics at high energies. The deep-valley air-shower array concept that we present provides highly background-suppressed neutrino detection with pointing resolution <1 degree, allowing us to begin the era of high-energy tau-neutrino astronomy., Comment: 15 pages, 6 figures, white paper submitted to the Latin American Strategy Forum for Research Infrastructure (LASF4RI)

NuRadioMC is a Monte Carlo framework designed to simulate ultra-high energy neutrino detectors that rely on the radio detection method. This method exploits the radio emission generated in the electromagnetic component of a particle shower following a neutrino interaction. NuRadioMC simulates everything from the neutrino interaction in a medium, the subsequent Askaryan radio emission, the propagation of the radio signal to the detector and finally the detector response. NuRadioMC is designed as a modern, modular Python-based framework, combining flexibility in detector design with user-friendliness. It includes a state-of-the-art event generator, an improved modelling of the radio emission, a revisited approach to signal propagation and increased flexibility and precision in the detector simulation. This paper focuses on the implemented physics processes and their implications for detector design. A variety of models and parameterizations for the radio emission of neutrino-induced showers are compared and reviewed. Comprehensive examples are used to discuss the capabilities of the code and different aspects of instrumental design decisions., Comment: replaced with published version

We report an error found during independent review of the NuTauSim publicly available code \url{https://github.com/harmscho/NuTauSim} that forms the basis of this publication. The error in the code was in tracking the density of the medium during particle propagation. After the first interaction, the code was referencing the depth of penetration back to the surface of the Earth rather than the location of the last interaction. The results were obtained using densities that were systematically underestimated when the particle was traversing the inner layers of the Earth by assigning the density of either ice or bedrock, depending on the particle energy or ice thickness of the simulation. This error was fixed and the repository updated on September 29, 2018., Comment: erratum to arXiv:1707.00334

The Giant Radio Array for Neutrino Detection (GRAND) is a planned large-scale observatory of ultra-high-energy (UHE) cosmic particles, with energies exceeding 10^8 GeV. Its goal is to solve the long-standing mystery of the origin of UHE cosmic rays. To do this, GRAND will detect an unprecedented number of UHE cosmic rays and search for the undiscovered UHE neutrinos and gamma rays associated to them with unmatched sensitivity. GRAND will use large arrays of antennas to detect the radio emission coming from extensive air showers initiated by UHE particles in the atmosphere. Its design is modular: 20 separate, independent sub-arrays, each of 10 000 radio antennas deployed over 10 000 km^2. A staged construction plan will validate key detection techniques while achieving important science goals early. Here we present the science goals, detection strategy, preliminary design, performance goals, and construction plans for GRAND., Comment: 45 pages, 34 figures. Matches published version

We present a new methodology to discriminate between light and heavy ultra-high energy cosmic-ray primaries on an event-by-event basis using information from the radio detection of extensive air showers at MHz frequencies. Similarly to other methods to determine primary cosmic ray composition, the one presented here is based on comparisons between detected radio signals and Monte Carlo simulations for multiple primary cosmic ray compositions. Unlike other methods that first reconstruct the depth of maximum shower development $X_{\rm max}$ to relate it to the nature of the primaries, we instead infer the cosmic-ray composition directly. The method is most effective in the case of inclined showers that arrive at large zenith angles with respect to the vertical to the ground, where methods based on the determination of $X_{\rm max}$ lose accuracy. We show that a discrimination efficiency between 65% and 80% can be reached for zenith angles $\theta \gtrsim 60^{\circ}$, even when typical uncertainties in radio detection are taken into account, including shower energy uncertainty. Our methodology could in principle be applied in large and sparse radio arrays, designed with the large radio footprint of inclined showers in mind, to significantly increase the statistics of ultra-high energy cosmic-ray composition studies.

The Giant Radio Array for Neutrino Detection (GRAND) aims at detecting ultra-high energy extraterrestrial neutrinos via the extensive air showers induced by the decay of tau leptons created in the interaction of neutrinos under the Earth's surface. Consisting of an array of $\sim10^5$ radio antennas deployed over $\sim 2\times10^5\,\rm {km}^2$, GRAND plans to reach, for the first time, an all-flavor sensitivity of $\sim1.5\times10^{-10} \,\rm GeV\, cm^{-2} \,s^{-1}\, sr^{-1}$ above $5\times10^{17}$ eV and a sub-degree angular resolution, beyond the reach of other planned detectors. We describe here preliminary designs and simulation results, plans for the ongoing, staged approach to the construction of GRAND, and the rich research program made possible by GRAND's design sensitivity and angular resolution., Comment: Proceedings of the 35th International Cosmic Ray Conference (ICRC2017), Bexco, Busan, Korea

There has been a recent surge in interest in the detection of $\tau$ lepton-induced air showers from detectors at altitude. When a $\tau$ neutrino ($\nu_\tau$) enters the Earth it produces $\tau$ leptons as a result of nuclear charged current interactions. In some cases, this process results in a $\tau$ lepton exiting the surface of the Earth, which can subsequently decay in the atmosphere and produce an extensive air shower. These upward-going air showers can be detected via fluorescence, optical Cherenkov, or geomagnetic radio emission. Several experiments have been proposed to detect these signals. We present a comprehensive simulation of the production of $\tau$ leptons by $\nu_\tau$'s propagating through Earth to aid the design of future experiments. These simulations for $\nu_\tau$'s and leptons in the energy range from $10^{15}$ eV to $10^{21}$ eV treat the full range of incidence angles from Earth-skimming to diametrically-traversing. Propagation of $\nu_\tau$'s and leptons include the effects of rock and an ocean or ice layer of various thicknesses. The interaction models include $\nu_\tau$ regeneration and account for uncertainties in the Standard Model neutrino cross-section and in the photo-nuclear contribution to the $\tau$ energy loss rate., Comment: 13 pages, 14 figures

The lunar technique is a method for maximising the collection area for ultra-high-energy (UHE) cosmic ray and neutrino searches. The method uses either ground-based radio telescopes or lunar orbiters to search for Askaryan emission from particles cascading near the lunar surface. While experiments using the technique have made important advances in the detection of nanosecond-scale pulses, only at the very highest energies has the lunar technique achieved competitive limits. This is expected to change with the advent of the Square Kilometre Array (SKA), the low-frequency component of which (SKA-low) is predicted to be able to detect an unprecedented number of UHE cosmic rays. In this contribution, the status of lunar particle detection is reviewed, with particular attention paid to outstanding theoretical questions, and the technical challenges of using a giant radio array to search for nanosecond pulses. The activities of SKA's High Energy Cosmic Particles Focus Group are described, as is a roadmap by which this group plans to incorporate this detection mode into SKA-low observations. Estimates for the sensitivity of SKA-low phases 1 and 2 to UHE particles are given, along with the achievable science goals with each stage. Prospects for near-future observations with other instruments are also described., Comment: 6 pages, 3 figures, presented at the 7th International Conference on Acoustic and Radio EeV Neutrino Detection Activities (ARENA 2016), Groningen, The Netherlands

The lunar Askaryan technique is a method to study the highest-energy cosmic rays, and their predicted counterparts, the ultra-high-energy neutrinos. By observing the Moon with a radio telescope, and searching for the characteristic nanosecond-scale Askaryan pulses emitted when a high-energy particle interacts in the outer layers of the Moon, the visible lunar surface can be used as a detection area. Several previous experiments, at Parkes, Goldstone, Kalyazin, Westerbork, the ATCA, Lovell, LOFAR, and the VLA, have developed the necessary techniques to search for these pulses, but existing instruments have lacked the necessary sensitivity to detect the known flux of cosmic rays from such a distance. This will change with the advent of the SKA. The Square Kilometre Array (SKA) will be the world's most powerful radio telescope. To be built in southern Africa, Australia and New Zealand during the next decade, it will have an unsurpassed sensitivity over the key 100 MHz to few-GHZ band. We introduce a planned experiment to use the SKA to observe the highest-energy cosmic rays and, potentially, neutrinos. The estimated event rate will be presented, along with the predicted energy and directional resolution. Prospects for directional studies with phase 1 of the SKA will be discussed, as will the major technical challenges to be overcome to make full use of this powerful instrument. Finally, we show how phase 2 of the SKA could provide a vast increase in the number of detected cosmic rays at the highest energies, and thus to provide new insight into their spectrum and origin., Comment: 8 pages, 6 figures. Proceedings of the 34th International Cosmic Ray Conference, Den Haag, The Netherlands, 2015

We investigate whether transition radiation from a particle shower crossing the interface between Earth and air and induced by an Earth-skimming neutrino can explain the upward event announced recently by the ANITA Collaboration. While the properties of the observed signal can in principle be explained with transition radiation, the flux necessary for a successful explanation is in tension with the current best limits from the Pierre Auger Observatory, the IceCube neutrino detector and the ANITA balloon. We also show that in this scenario, the direction of the incoming neutrino is determined precisely to within a few degrees, from the polarization properties of the observed events and thanks to the Earth opacity to ultra high energy neutrinos., Comment: 8 pages, 4 figures, version accepted to PRD

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., Comment: Replaced with published version. Added journal reference and DOI. Supplemental material in the ancillary files

The IceCube, Pierre Auger and Telescope Array Collaborations have recently reported results on neutral particles (neutrons, photons and neutrinos) which complement the measurements on charged primary cosmic rays at ultra-high energy. The complementarity between these messengers and between their detections are outlined. The current status of their search is reviewed and a cross-correlation analysis between the available results is performed. The expectations for photon and neutrino detections in the near future are also presented., Comment: 13 pages, 5 figures, Proceeding of UHECR 2014 conference, Springdale, Utah, USA; accepted for publication in JPS Conference Proceedings

Coherent radiation at radio frequencies from high-energy showers fully contained in a dense radio-transparent medium - like ice, salt or regolith - has been extensively investigated as a promising technique to search for ultra-high energy (UHE) neutrinos. Additional emission in the form of transition radiation may occur when a neutrino-induced shower produced close to the Earth surface emerges from the ground into atmospheric air. We present the first detailed evaluation of transition radiation from high-energy showers crossing the boundary between two different media. We found that transition radiation is sizable over a wide solid angle and coherent up to $\sim$ 1 GHz. These properties encourage further work to evaluate the potential of a large-aperture UHE neutrino experiment based on detection of transition radiation., Comment: 13 pages, 14 figures, matches PRD accepted version, typo for the emission strength for EeV showers corrected

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal., Comment: Replaced with published version. Added journal reference and DOI

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}}$., Comment: Replaced with published version. Added journal reference and DOI

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., Comment: Replaced with published version. Added journal reference and DOI

We present the calculation of coherent radio pulses emitted by extensive air showers induced by ultra-high energy cosmic rays accounting for reflection on the Earth's surface. Results have been obtained with a simulation program that calculates the contributions from shower particles after reflection at a surface plane. The properties of the radiation are discussed in detail emphasizing the effects of reflection. The shape of the frequency spectrum is shown to be closely related to the angle of the observer with respect to shower axis, becoming hardest in the Cherenkov direction. The intensity of the flux at a fixed observation angle is shown to scale with the square of the primary particle energy to very good accuracy indicating the coherent aspect of the emission. The simulation methods of this paper provide the foundations for energy reconstruction of experiments looking at the Earth from balloons and satellites. They can also be used in dedicated studies of existing and future experimental proposals.

Radio emission in atmospheric showers is currently interpreted in terms of radiation due to the deviation of the charged particles in the magnetic field of the Earth and to the charge excess (Askaryan effect). Each of these mechanisms has a distinctive polarization. The complex signal patterns can be qualitatively explained as the interference (superposition) of the fields induced by each mechanism. In this work we explicitly and quantitatively test a simple phenomenological model based on this idea. The model is constructed by isolating each of the two components at the simulation level and by making use of approximate symmetries for each of the contributions separately. The results of the model are then checked against full ZHAireS Monte Carlo simulations of the electric field calculated from first principles. We show that the simple model describes radio emission at a few percent level in a wide range of shower-observer geometries and on a shower-by-shower basis. As a consequence, this approach provides a simple method to reduce the computing time needed to accurately predict the electric field of radio pulses emitted from air showers, with many practical applications in experimental situations of interest., Comment: Replaced with revised version published in Astroparticle Physics

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

The detection of electromagnetic pulses from high energy showers is used as a means to search for Ultra-High Energy cosmic ray and neutrino interactions. An approximate formula has been obtained to numerically evaluate the radio pulse emitted by a charged particle that instantaneously accelerates, moves at constant speed along a straight track and halts again instantaneously. The approximate solution is applied to the particle track after dividing it in smaller subintervals. The resulting algorithm (often referred to as the ZHS algorithm) is also the basis for most of the simulations of the electric field produced in high energy showers in dense media. In this work, the electromagnetic pulses as predicted with the ZHS algorithm are compared to those obtained with an exact solution of the electric field produced by a charged particle track. The precise conditions that must apply for the algorithm to be valid are discussed and its accuracy is addressed. This comparison is also made for electromagnetic showers in dense media. The ZHS algorithm is shown to describe Cherenkov radiation and to be valid for most situations of interest concerning detectors searching for Ultra-High Energy neutrinos. The results of this work are also relevant for the simulation of pulses emitted from air showers., Comment: 12 pages, 7 figures

Using detailed Monte Carlo simulations we have characterized the features of the radio emission of inclined air showers in the Ultra-High Frequency band (300 MHz - 3 GHz). The Fourier-spectrum of the radiation is shown to have a sizable intensity well into the GHz frequency range. The emission is mainly due to transverse currents induced by the geomagnetic field and to the excess charge produced by the Askaryan effect. At these frequencies only a significantly reduced volume of the shower around the axis contributes coherently to the signal observed on the ground. The size of the coherently emitting volume depends on frequency, shower geometry and observer position, and is interpreted in terms of the relative time delays. At ground level, the maximum emission at high frequencies is concentrated in an elliptical ring-like region around the intersection of a Cherenkov cone with its vertex at shower maximum and the ground. The frequency spectrum of inclined showers when observed at positions that view shower maximum in the Cherenkov direction, is shown to be in broad agreement with the pulses detected by the Antarctic Impulsive Transient Antenna (ANITA) experiment, making the interpretation that they are due to Ultra-High Energy Cosmic Ray atmospheric showers consistent with our simulations. These results are also of great importance for experiments aiming to detect molecular bremsstrahlung radiation in the GHz range as they present an important background for its detection., Comment: 8 pages, 8 figures

We present predictions for the radio pulses emitted by extensive air showers using ZHAireS, an AIRES-based Monte Carlo code that takes into account the full complexity of ultra-high energy cosmic-ray induced shower development in the atmosphere, and allows the calculation of the electric field in both the time and frequency domains. We do not presuppose any emission mechanism, but our results are compatible with a superposition of geomagnetic and charge excess radio emission effects. We investigate the polarization of the electric field as well as the effects of the refractive index n and shower geometry on the radio pulses. We show that geometry, coupled to the relativistic effects that appear when using a realistic n>1, play a prominent role on the radio emission of air showers., Comment: Replaced with version published in Astroparticle Physics

An in-depth characterization of coherent radio Cherenkov pulses from particle showers in dense dielectric media, referred to as the Askaryan effect, is presented. The time-domain calculation developed in this article is based on a form factor to account for the lateral dimensions of the shower. It is computationally efficient and able to reproduce the results of detailed particle shower simulations with high fidelity in most regions of practical interest including Fresnel effects due to the longitudinal development of the shower. In addition, an intuitive interpretation of the characteristics of the Askaryan pulse is provided. We expect our approach to benefit the analysis of radio pulses in experiments exploiting the radio technique., Comment: Replaced with version published Phys. Rev. D

The Cherenkov radio pulse emitted by hadronic showers in ice is calculated for showers of energies in the EeV range. This is obtained with three dimensional simulations of both shower development and the coherent radio pulse emitted as the excess charge develops in the shower. A Monte Carlo, ZHAireS, has been developed for this purpose combining the high energy hadronic interaction capabilities of AIRES, and the dense media propagation capabilities of TIERRAS, with the precise low energy tracking and specific algorithms developed to calculate the radio emission in ZHS. A thinning technique is implemented and optimized to allow the simulation of radio pulses induced by showers up to 10 EeV in ice. The code is validated comparing the results for electromagnetic and hadronic showers to those obtained with GEANT4 and ZHS codes. The contribution to the pulse of other shower particles in addition to electrons and positrons, mainly pions and muons, is found to be below 1%. The characteristics of hadronic showers and the corresponding Cherenkov frequency spectra are compared with those from purely electromagnetic showers. The dependence of the spectra on shower energy and high-energy hadronic model is addressed and parameterizations for the radio emission in hadronic showers in ice are given for practical applications., Comment: Replaced with version published in Astroparticle Physics

The electric field of the Cherenkov radio pulse produced by a single charged particle track in a dielectric medium is derived from first principles. An algorithm is developed to obtain the pulse in the time domain for numerical calculations. The algorithm is implemented in a Monte Carlo simulation of electromagnetic showers in dense media (specifically designed for coherent radio emission applications) as might be induced by interactions of ultra-high energy neutrinos. The coherent Cherenkov radio emission produced by such showers is obtained simultaneously both in the time and frequency domains. A consistency check performed by Fourier-transforming the pulse in time and comparing it to the frequency spectrum obtained directly in the simulations yields, as expected, fully consistent results. The reversal of the time structure inside the Cherenkov cone and the signs of the corresponding pulses are addressed in detail. The results, besides testing algorithms used for reference calculations in the frequency domain, shed new light into the properties of the radio pulse in the time domain. The shape of the pulse in the time domain is directly related to the depth development of the excess charge in the shower and its width to the observation angle with respect to the Cherenkov direction. This information can be of great practical importance for interpreting actual data., Comment: 10 pages, 4 figures

We study the frequency and angular dependences of Cherenkov radio pulses originated by the excess of electrons in electromagnetic showers in different dense media. We develop a simple model to relate the main characteristics of the electric field spectrum to the properties of the shower such as longitudinal and lateral development. This model allows us to establish the scaling of the electric field spectrum with the properties of the medium such as density, radiation length, Moliere radius, critical energy and refraction index. We normalize the predictions of the scaling relations to the numerical results obtained in our own developed GEANT4-based Monte Carlo simulation, and we give a unified parameterization of the frequency spectrum and angular distribution of the electric field in ice, salt, and the lunar regolith, in terms of the relevant properties of the media. Our parameterizations are valid for electromagnetic showers below the energy at which the Landau-Pomeranchuk-Migdal effect starts to be relevant in these media. They also provide an approximate estimate of radio emission in hadronic showers induced by high energy cosmic rays or neutrinos., Comment: 7 pages, 4 figures

The magnetic field in our Galaxy is not well known and is difficult to measure. A spiral regular field in the disk between the Galactic arms is favored by observations, however it is still controversial if the field reverses from arm to arm. The parity of the field across the Galactic plane is also not well established. In this letter we demonstrate that cosmic ray protons in the energy range 10^18 to 10^19 eV, if accelerated near the center of the Galaxy, can probe the large scale structure of the Galactic Magnetic Field. In particular if the field is of even parity, and the spiral field reverses direction from arm to arm, i.e. if it is bi-symmetric (BSS), ultra high energy protons will predominantly come from the Southern Galactic hemisphere, and predominantly from the Northern Galactic hemisphere if the field is of even parity and axi-symmetric (ASS). There is no sensitivity to the BSS or ASS configurations if the field is of odd parity., Comment: 6 pages, 3 tables, 2 figures

Most active galactic nuclei (AGN) lack prominent jets, and show modest radio emission and significant X-ray emission which arises mainly from the galactic core, very near from the central black hole. We use a quantitative scenario of such core-dominated radio-quiet AGN, which attributes a substantial fraction of the X-ray emission to the presence of abortive jets involving the collision of gas blobs in the core. Here we investigate the consequences of the acceleration of protons in the shocks from such collisions. We find that protons will be accelerated up to energies above the pion photoproduction threshold on both the X-rays and the UV photons from the accretion disk. The secondary charged pions decay, producing neutrinos. We predict significant fluxes of TeV-PeV neutrinos, and show that the AMANDA II detector is already constraining several important astrophysical parameters of these sources. Larger cubic kilometer detectors such as IceCube will be able to detect such neutrinos in less than one year of operation, or otherwise rule out this scenario., Comment: 10 pages, 5 Postscript figures

We study the influence of shower fluctuations, and the possible presence of different nuclear species in the primary cosmic ray spectrum, on the experimental determination of both shower energy and the proton air inelastic cross section from studies of the longitudinal development of atmospheric showers in fluorescence experiments. We investigate the potential of track length integral and shower size at maximum as estimators of shower energy. We find that at very high energy (~10^19-10^20 eV) the error of the total energy assignment is dominated by the dependence on the hadronic interaction model, and is of the order of 5%. At lower energy (~10^17-10^18 eV), the uncertainty of the energy determination due to the limited knowledge of the primary cosmic ray composition is more important. The distribution of depth of shower maximum is discussed as a measure of the proton-air cross section. Uncertainties in a possible experimental measurement of this cross section introduced by intrinsic shower fluctuations, the model of hadronic interactions, and the unknown mixture of primary nuclei in the cosmic radiation are numerically evaluated., Comment: 12 pages, 11 figures, 4 tables

GRB941017, a gamma-ray burst of exceptional fluence, has recently been shown to have a high-energy component which is not consistent with the standard fireball phenomenology. If this component is the result of photomeson interactions in the burst fireball, it provides new and compelling support for substantial high-energy neutrino fluxes from this and similar sources. In this letter, we consider what impact this new information has on the neutrino spectra of gamma-ray bursts and discuss how this new evidence impacts the prospects for detection of such events in next generation neutrino telescopes.

We explore the constant intensity cut method that is widely used for the derivation of the cosmic ray energy spectrum, for comparisons of data obtained at different atmospheric depths, for measuring average shower profiles, and for estimates of the proton-air cross section from extensive air shower data. The constant intensity cut method is based on the selection of air showers by charged particle or muon size and therefore is subject to intrinsic shower fluctuations. We demostrate that, depending on the selection method, shower fluctuations can strongly influence the characteristics of the selected showers. Furthermore, a mixture of different primaries in the cosmic ray flux complicates the interpretation of measurements based on the method of constant intensity cuts. As an example we consider data published by the Akeno Collaboration. The interpretation of the Akeno measurements suggests that more than 60-70% of cosmic ray primaries in the energy range 10^{16}-19^{17} eV are heavy nuclei. Our conclusions depend only weakly on the hadronic interaction model chosen to perform the simulations, namely SIBYLL and QGSjet., Comment: 10 pages, 10 figures, RevTeX4

The observation of TeV-gamma rays of neutral pion origin from the supernova remnant RX J1713.7-3946 might have revealed the first specific site where protons are accelerated to energies typical of the main component of the cosmic rays. In this letter we calculate the high-energy neutrino flux associated with this source to be at least 40 muon-type neutrinos per kilometer-squared per year. We perform the same calculations for other known sources of TeV-gamma rays and show how neutrino observations can establish whether the TeV-gamma rays emitted by blazars and supernova remnants are the decay products of neutral pions and thus unequivocally establish the sources as cosmic accelerators., Comment: Replaced with version published in Astrophys. J

We present a fast one dimensional hybrid method to efficiently simulate extensive air showers up to the highest observed energies. Based on precalculated pion showers and a bootstrap technique, our method predicts the average shower profile, the number of muons at detector level above several energy thresholds as well as the fluctuations of the electromagnetic and hadronic components of the shower. We study the main characteristics of proton-induced air showers up to ultra-high energy, comparing the predictions of three different hadronic interaction models: SIBYLL 1.7, SIBYLL 2.1 and QGSjet98. The influence of the hadronic interaction models on the shower evolution, in particular the elongation rate, is discussed and the applicability of analytical approximations is investigated., Comment: 17 pages (Latex), 15 postscript figures, 7 tables

If spacetime has more than four dimensions, ultra-high energy cosmic rays may create microscopic black holes. Black holes created by cosmic neutrinos in the Earth will evaporate, and the resulting hadronic showers, muons, and taus may be detected in neutrino telescopes below the Earth's surface. We simulate such events in detail and consider black hole cross sections with and without an exponential suppression factor. We find observable rates in both cases: for conservative cosmogenic neutrino fluxes, several black hole events per year are observable at the IceCube detector; for fluxes at the Waxman-Bahcall bound, tens of events per year are possible. We also present zenith angle and energy distributions for all three channels. The ability of neutrino telescopes to differentiate hadrons, muons, and possibly taus, and to measure these distributions provides a unique opportunity to identify black holes, to experimentally constrain the form of black hole production cross sections, and to study Hawking evaporation., Comment: 20 pages, 9 figures

Recently the AGASA Collaboration presented data suggesting a significant clustering of ultra-high energy cosmic rays coming from the outer Galaxy region. In this paper we calculate expected cosmic ray arrival distributions for several simple, limiting source location scenarios and investigate the possibility of clustering and correlation effects. The role of the Galactic magnetic field is discussed in detail., Comment: 33 pages, 21 figures (5 in color), 2 tables