Search

The main challenge in detecting ultra-high energy (UHE) neutrinos is discriminating a neutrino-induced shower in the background of showers initiated by ultra-high energy nuclei. The resulting shower development from neutrinos exhibits different characteristics from hadron-induced showers because neutrinos penetrate the atmosphere more deeply than hadrons. This study focuses on simulations of highly inclined neutrino-induced showers above $75^\circ$ zenith angles, exploring an extensive energy range from $1 \text{EeV}$ to $120 \text{EeV}$. These simulated showers have different ranges of interaction depths corresponding to each zenith angle, presenting diverse detection challenges. Our methodology utilises timing data from radio antennas for the shower front calculation for extensive air showers induced by neutrinos and nuclei. Furthermore, we incorporate signals obtained from Water Cherenkov detectors and the spatial distribution of stations registering signals in both Water Cherenkov detectors and radio antennas. We aim to classify neutrino-induced showers and background events stemming from nuclei by harnessing a decision tree classifier employing the Gini impurity method. Our framework yields excellent accuracy for separating the neutrinos from the background. The findings of this study offer significant advancements in the domain of UHE neutrino detection, shedding light on astrophysical phenomena associated with these elusive particles amidst the complex background of UHE nuclei.

Charged current interactions of neutrinos inside the Earth can result in secondary muons and {\tau} - leptons which are detectable by a large swath of 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 present RDSim, a fast and comprehensive framework for the simulation of the radio emission and detection of downgoing air showers. It can handle any downgoing shower that can be simulated with ZHAireS including those induced by CC and NC neutrino interactions and $\tau$ decays. RDSim is based on a superposition toymodel that disentangles the Askaryan and geomagnetic components of the shower emission. By using full ZHAireS simulations as input, it is able to estimate the full radio footprint on the ground. A single input simulation at a given energy and arrival direction can be scaled in energy and rotated in azimuth by taking into account all relevant effects. This makes it possible to simulate a huge number of geometries and energies using just a few ZHAireS input simulations. The framework takes into account the main characteristics of the detector, such as trigger setups, thresholds and antenna patterns. To accommodate arrays that use particle detectors for triggering, such as the Auger RD extension, it also features a second toymodel to estimate the muon density at ground level, which is used to perform simple particle trigger simulations. It's speed makes it possible to investigate in detail events with a very low trigger probability, as well as many geometrical effects due to the array layout. In case more detailed studies of the radio detection are needed, RDSim can also be used to sweep the phase-space for the efficient creation of dedicated full simulation sets. This is particularly important in the case of neutrino events, that have extra variables that greatly impact shower characteristics, such as interaction or $\tau$ decay depth as well as the type of interaction and it's fluctuations., Comment: arXiv admin note: text overlap with arXiv:2307.07351

RDSim is a fast, accurate and flexible framework for the simulation of the radio emission of downgoing air showers and its detection by an arbitrary array, including showers initiated by neutrino interactions or tau-lepton decays. RDSim was build around speed and is based on simple and fast, yet still accurate, toymodel-like approaches. It models the radio emission using a superposition emission model that disentangles the Askaryan and geomagnetic components of the shower radio emission. It uses full ZHAireS simulations as an input to estimate the electric field at any position on the ground. A single input simulation can be scaled in energy and rotated in azimuth, taking into account all relevant effects. This makes it possible to simulate a huge number of geometries and energies using just a few ZHAireS input simulations. RDSim takes into account the main characteristics of the detector, such as trigger setups, thresholds and antenna patterns. To accommodate arrays that use particle detectors for triggering, such as the Auger RD extension, it also features a second toymodel to estimate the muon density at ground level and perform simple particle trigger simulations. Owing to the large statistics made possible by its speed, it can be used to investigate in detail events with a very low trigger probability and geometrical effects due to the array layout, making it specially suited to be used as a fast and accurate aperture calculator. In case more detailed studies of the radio emission and detector response are desired, RDSim can also be used to sweep the phase-space for the efficient creation of dedicated full simulation sets. This is particularly important in the case of neutrino events, that have extra variables that greatly impact shower characteristics, such as interaction or $\tau$ decay depth as well as the type of interaction and it's fluctuations.

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. When time synchronisation below $\sim$1 ns between different receivers can be achieved, 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 proposed method can be incorporated in operational and future cosmic particle observatories and with its high accuracy it has the potential to play a crucial role in unravelling the composition of the ultra-high-energy cosmic-particle flux., Comment: 8 pages, 4 figures. Added figures/references with respect to v1

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

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)

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.

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

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

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

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

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., Comment: 9 page 9 figures