• A Comprehensive Approach to Tau-Lepton Production by High-Energy Tau Neutrinos Propagating Through Earth(1707.00334)

Dec. 2, 2017 astro-ph.HE
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): Present and Perspectives(1708.05128)

Aug. 17, 2017 astro-ph.IM, astro-ph.HE
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.
• Overview of lunar detection of ultra-high energy particles and new plans for the SKA(1704.05336)

April 18, 2017 astro-ph.IM, astro-ph.HE
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.
• Can transition radiation explain the ANITA event 3985267?(1606.07059)

Jan. 30, 2017 astro-ph.HE
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.
• Radio pulses from ultra-high energy atmospheric showers as the superposition of Askaryan and geomagnetic mechanisms(1402.3504)

Sept. 20, 2016 astro-ph.HE
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 lunar Askaryan technique with the Square Kilometre Array(1608.02408)

Aug. 8, 2016 astro-ph.IM, astro-ph.HE
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.
• 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.

March 21, 2016 hep-ph, astro-ph.IM, astro-ph.HE
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.
• 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.
• Report from the Multi-Messenger Working Group at UHECR-2014 Conference(1510.02050)

Nov. 12, 2015 astro-ph.IM, astro-ph.HE
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.
• Simulations of reflected radio signals from cosmic ray induced air showers(1502.02117)

Feb. 7, 2015 astro-ph.IM, astro-ph.HE
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.
• Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
• Calculations of electric fields for radio detection of Ultra-High Energy particles(1210.1052)

Oct. 3, 2012 astro-ph.HE
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.
• Coherent Radiation from Extensive Air Showers in the Ultra-High Frequency Band(1208.0951)

Aug. 4, 2012 astro-ph.HE
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.
• Monte Carlo simulations of radio pulses in atmospheric showers using ZHAireS(1107.1189)

April 16, 2012 astro-ph.HE
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.

April 16, 2012 astro-ph.HE
• Coherent Cherenkov radio pulses from hadronic showers up to EeV energies(1005.0552)

April 16, 2012 astro-ph.HE
• Cherenkov radio pulses from electromagnetic showers in the time-domain(1002.3873)

Feb. 20, 2010 astro-ph.HE
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.
• Coherent radio pulses from showers in different media: A unified parameterization(astro-ph/0512337)

Dec. 13, 2005 astro-ph
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.
• The large scale structure of the Galactic magnetic field and High Energy Cosmic Ray anisotropy(astro-ph/0507273)

July 12, 2005 astro-ph
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.
• High energy neutrinos from radio-quiet AGNs(astro-ph/0409034)

Sept. 1, 2004 astro-ph, hep-ph
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.
• Influence of shower fluctuations and primary composition on studies of the shower longitudinal development(astro-ph/0402092)

Feb. 4, 2004 astro-ph
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.
• GRB941017: A Case Study of Neutrino Production in Gamma Ray Bursts(astro-ph/0310417)

Oct. 15, 2003 astro-ph, hep-ph
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.