We developed a radio interferometric technique for the observation of extensive air showers initiated by cosmic particles. In this proof-of-principle study we show that properties of extensive air showers can be derived with high accuracy in a straightforward manner. Direction reconstruction resolution of $< 0.2^\circ$ and resolution on the depth of shower maximum of $<10$\,g/cm$^2$ are obtained over the full parameter range studied, with even higher accuracy for inclined incoming directions. In addition, by applying the developed method to dense arrays of radio antennas, the energy threshold for the radio detection of extensive air showers can be significantly lowered. The method can be applied to several operational experiments and offers good prospects for planned cosmic particle observatories.

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.

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.

Over the last years, radio detection has matured to become a competitive method for the detection of air showers. Arrays of thousands of antennas are now envisioned for the detection of cosmic rays of ultra high energy or neutrinos of astrophysical origin. The data exploitation of such detectors requires to run massive air-shower simulations to evaluate the radio signal at each antenna position. In order to reduce the associated computational cost, we have developed a semi-analytical method for the computation of the emitted radio signal called Radio Morphing. The method consists in computing the radio signal of any air-shower at any location from the simulation of one single reference shower at given positions by i) a scaling of the electric-field amplitude of this reference shower, ii) an isometry on the simulated positions and iii) an interpolation of the radio pulse at the desired position. This technique enables one to compute electric field time traces with characteristics very similar to those obtained with standard computation methods, but with computation times reduced by several orders of magnitude. In this paper, we present this novel tool, explain its methodology, and discuss its limitations. Furthermore, we validate the method on a typical event set for the future GRAND experiment showing that the calculated peak amplitudes are consistent with the results from ZHAireS simulations with a mean offset of +8.5% and a standard deviation of 27.2% in this specific case. This overestimation of the signal strength by Radio Morphing arises mainly from the choice of the underlying reference shower.

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.

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.

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.

Recently, the ANITA collaboration reported on two upward-going extensive air shower events consistent with a primary particle that emerges from the surface of the ice. These events may be of $\nu_\tau$ origin, in which the neutrino interacts within the Earth to produce a $\tau$ lepton that emerges from the Earth, decays in the atmosphere, and initiates an extensive air shower. In this paper we estimate an upper bound on the ANITA acceptance to a diffuse $\nu_\tau$ flux detected via $\tau$-lepton-induced air showers within the bounds of Standard Model (SM) uncertainties. By comparing this estimate with the acceptance of Pierre Auger Observatory and IceCube and assuming SM interactions, we conclude that a $\nu_\tau$ origin of these events would imply a neutrino flux at least two orders of magnitude above current bounds.

Over the last decades, radio detection of air showers has been established as a promising detection technique for ultrahigh-energy cosmic rays and neutrinos. Very large or dense antenna arrays are necessary to be proficient at collecting information about these particles and understanding their properties accurately. The exploitation of such arrays requires to run massive air-shower simulations to evaluate the radio signal at each antenna position, taking into account features such as the ground topology. In order to reduce computational costs, we have developed a fast computation of the emitted radio signal on the basis of generic shower simulations, called Radio Morphing. The method consists in the calculation of the radio signal of any air-shower by i) a scaling of the electric-field amplitude of a reference air shower to the target shower, ii) an isometry on the simulated positions and iii) an interpolation of the radio pulse at the desired position. This technique enables one to gain many orders of magnitude in CPU time compared to a standard computation. In this contribution, we present this novel tool and explain its methodology. In particular, Radio Morphing will be a key element for the simulation chain of the Giant Radio Array for Neutrino Detection (GRAND) project, that aims at detecting ultra-high-energy neutrinos with an array of 200 000 radio antennas in mountainous regions.os with an array of 200 000 radio antennas in mountainous regions.

These proceedings address a recent publication by the ANITA collaboration of four upward- pointing cosmic-ray-like events observed in the first flight of ANITA. Three of these events were consistent with stratospheric cosmic-ray air showers where the axis of propagation does not inter- sect the surface of the Earth. The fourth event was consistent with a primary particle that emerges from the surface of the ice suggesting a possible {\tau}-lepton decay as the origin of this event. These proceedings follow-up on the modeling and testing of the hypothesis that this event was of {\tau} neutrino origin.

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.

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.

The Giant Radio Array for Neutrino Detection (GRAND) is a planned array of ~200 000 radio antennas deployed over ~200 000 km2 in a mountainous site. It aims primarly at detecting high-energy neutrinos via the observation of extensive air showers induced by the decay in the atmosphere of taus produced by the interaction of cosmic neutrinos under the Earth surface. GRAND aims at reaching a neutrino sensitivity of 5.10$^{11}$ E$^{-2}$ GeV$^{-1}$cm$^{-2}$s$^{-1}$sr$^{-1}$ above 3.10$^{16}$ eV. This ensures the detection of cosmogenic neutrinos in the most pessimistic source models, and ~50 events per year are expected for the standard models. The instrument will also detect UHECRs and possibly FRBs. Here we show how our preliminary design should enable us to reach our sensitivity goals, and discuss the steps to be taken to achieve GRAND.

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.

The first flight of the Antarctic Impulsive Transient Antenna (ANITA) experiment recorded 16 radio signals that were emitted by cosmic-ray induced air showers. For 14 of these events, this radiation was reflected from the ice. The dominant contribution to the radiation from the deflection of positrons and electrons in the geomagnetic field, which is beamed in the direction of motion of the air shower. This radiation is reflected from the ice and subsequently detected by the ANITA experiment at a flight altitude of 36km. In this paper, we estimate the energy of the 14 individual events and find that the mean energy of the cosmic-ray sample is 2.9 EeV. By simulating the ANITA flight, we calculate its exposure for ultra-high energy cosmic rays. We estimate for the first time the cosmic-ray flux derived only from radio observations. In addition, we find that the Monte Carlo simulation of the ANITA data set is in agreement with the total number of observed events and with the properties of those events.

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.

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.

We present a search for microwave emission from air showers induced by ultrahigh energy cosmic rays with the microwave detection of air showers experiment. No events were found, ruling out a wide range of power flux and coherence of the putative emission, including those suggested by recent laboratory measurements.

We present the design, implementation and data taking performance of the MIcrowave Detection of Air Showers (MIDAS) experiment, a large field of view imaging telescope designed to detect microwave radiation from extensive air showers induced by ultra-high energy cosmic rays. This novel technique may bring a tenfold increase in detector duty cycle when compared to the standard fluorescence technique based on detection of ultraviolet photons. The MIDAS telescope consists of a 4.5 m diameter dish with a 53-pixel receiver camera, instrumented with feed horns operating in the commercial extended C-Band (3.4 -- 4.2 GHz). A self-trigger capability is implemented in the digital electronics. The main objectives of this first prototype of the MIDAS telescope - to validate the telescope design, and to demonstrate a large detector duty cycle - were successfully accomplished in a dedicated data taking run at the University of Chicago campus prior to installation at the Pierre Auger Observatory.

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.

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.

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.

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.

Proceedings of the 13th International Conference on Elastic and Diffractive Scattering (Blois Workshop) - Moving Forward into the LHC Era

Recent measurements suggest that extensive air showers initiated by ultra-high energy cosmic rays (UHECR) emit signals in the microwave band of the electromagnetic spectrum caused by the collisions of the free-electrons with the atmospheric neutral molecules in the plasma produced by the passage of the shower. Such emission is isotropic and could allow the detection of air showers with 100% duty cycle and a calorimetric-like energy measurement, a significant improvement over current detection techniques. We have built MIDAS (MIcrowave Detection of Air Showers), a prototype of microwave detector, which consists of a 4.5 m diameter antenna with a cluster of 53 feed-horns in the 4 GHz range. The details of the prototype and first results will be presented.

Recent measurements suggest free electrons created in ultra-high energy cosmic ray extensive air showers (EAS) can interact with neutral air molecules producing Bremsstrahlung radiation in the microwave regime. The microwave radiation produced is expected to scale with the number of free electrons in the shower, which itself is a function of the energy of the primary particle and atmospheric depth. Using these properties a calorimetric measurement of the EAS is possible. This technique is analogous to fluorescence detection with the added benefit of a nearly 100% duty cycle and practically no atmospheric attenuation. The Microwave Detection of Air Showers (MIDAS) prototype is currently being developed at the University of Chicago. MIDAS consists of a 53 feed receiver operating in the 3.4 to 4.2 GHz band. The camera is deployed on a 4.5 meter parabolic reflector and is instrumented with high speed power detectors and autonomous FPGA trigger electronics. We present the current status of the MIDAS instrument and an outlook for future development.

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.

Technical reports on operations and features of the Pierre Auger Observatory, including ongoing and planned enhancements and the status of the future northern hemisphere portion of the Observatory. Contributions to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009.

Studies of the correlations of ultra-high energy cosmic ray directions with extra-Galactic objects, of general anisotropy, of photons and neutrinos, and of other astrophysical effects, with the Pierre Auger Observatory. Contributions to the 31st ICRC, Lodz, Poland, July 2009.

Studies of the cosmic ray energy spectrum at the highest energies with the Pierre Auger Observatory.

Reports on the atmospheric monitoring, calibration, and other operating systems of the Pierre Auger Observatory. Contributions to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009.

Studies of the composition of the highest energy cosmic rays with the Pierre Auger Observatory, including examination of hadronic physics effects on the structure of extensive air showers. Submissions to the 31st ICRC, Lodz, Poland (July 2009).

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV (1 EeV = 10^18 eV) respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top-down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section.

We investigate the detection of exotic massive strongly interacting hadrons (uhecrons) in ultra high energy cosmic ray telescopes. The conclusion is that experiments such as the Pierre Auger Observatory have the potential to detect these particles. It is shown that uhecron showers have clear distinctive features when compared to proton and nuclear showers. The simulation of uhecron air showers, and its detection and reconstruction by fluorescence telescopes is described. We determine basic cuts in observables that will separate uhecrons from the cosmic ray bulk, assuming this is composed by protons. If these are composed by heavier nucleus the separation will be much improved. We also discuss photon induced showers. The complementarity between uhecron detection in accelerator experiments is discussed.

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.