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

This is an index of the contributions by the Giant Radio Array for Neutrino Detection (GRAND) Collaboration to the 10th International Workshop on Acoustic and Radio EeV Neutrino Detection Activities (ARENA 2024, University of Chicago, June 11-14, 2024). The contributions include an overview of GRAND in its present and future incarnations, methods of radio-detection that are being developed for them, and ongoing joint work between the GRAND and BEACON experiments., Comment: Note: To access the list of contributions, please follow the "HTML" link that can be found on the arXiv page

When performing radio detection of ultra-high energy astroparticles, the planar wavefront model is often used as a first step to evaluate the arrival direction of primary particles, by adjusting the wavefront orientation based on the peak of signal traces of individual antennas. The benefits of this approach are however limited by the lack of a good evaluation of its robustness. In order to mitigate this limitation, we derive in this work an analytical solution for the arrival direction as well as the corresponding analytical reconstruction uncertainty. Because it is fast and robust, this reconstruction method can be used as a proxy for more complex estimators, or implemented online for low-level triggering., Comment: Submitted to NIMA

The Giant Radio Array for Neutrino Detection (GRAND) aims to detect highly inclined extensive air showers (EAS) with down-going and up-going trajectories. Several working groups in the GRAND collaboration are developing methods to reconstruct the incoming direction, core position, primary energy, and composition of the showers. The reconstruction pipeline -- currently under development in the France/IAP working group -- relies on a model of spherical wavefront emission for arrival times, which is possible because the radio signals are generated far away from the antenna stations. The amplitude distribution of the signals at the antenna level is described by an Angular Distribution Function that considers various asymmetries in the data, including geomagnetic effects. In this contribution, we present preliminary results from testing our EAS reconstruction procedure using realistic mock observations., Comment: V2: Fixed typo in axis labels of Fig.5 and collaboration author list updated

The operation of upcoming ultra-high-energy cosmic-ray, gamma-ray, and neutrino radio-detection experiments, like the Giant Radio Array for Neutrino Detection (GRAND), poses significant computational challenges involving the production of numerous simulations of particle showers and their detection, and a high data throughput. GRANDlib is an open-source software tool designed to meet these challenges. Its primary goal is to perform end-to-end simulations of the detector operation, from the interaction of ultra-high-energy particles, through -- by interfacing with external air-shower simulations -- the ensuing particle shower development and its radio emission, to its detection by antenna arrays and its processing by data-acquisition systems. Additionally, GRANDlib manages the visualization, storage, and retrieval of experimental and simulated data. We present an overview of GRANDlib to serve as the basis of future GRAND analyses., Comment: 11 pages, 9 figures, plus appendices

Next-generation radio experiments such as the Radio Detector of the upgraded Pierre Auger Observatory and the planned GRAND and BEACON arrays target the detection of ultra-high-energy particle air showers arriving at low elevation angles. These inclined cosmic-ray air showers develop higher in the atmosphere than vertical ones, enhancing magnetic deflections of electrons and positrons inside the cascade. We evidence two novel features in their radio emission: a new polarization pattern, consistent with a geo-synchrotron emission model and a coherence loss of the radio emission, both for showers with zenith angle $\theta \gtrsim 65^{\circ}$ and strong enough magnetic field amplitude (typical strength of $B\sim 50\, \rm \mu T$). Our model is compared with both ZHAireS and CoREAS Monte-Carlo simulations. Our results break the cannonical description of a radio signal made of Askaryan and transverse current emission only, and provide guidelines for the detection and reconstruction strategies of next-generation experiments, including cosmic-ray/neutrino discrimination., Comment: 6 pages, 2 figures. Accepted in Physical Review Letter

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 Giant Radio Array for Neutrino Detection (GRAND) is an envisioned observatory of ultra-high-energy particles of cosmic origin, with energies in excess of 100 PeV. GRAND uses large surface arrays of autonomous radio-detection units to look for the radio emission from extensive air showers that are triggered by the interaction of ultra-high-energy cosmic rays, gamma rays, and neutrinos in the atmosphere or underground. In particular, for ultra-high-energy neutrinos, the future final phase of GRAND aims to be sensitive enough to discover them in spite of their plausibly tiny flux. Presently, three prototype GRAND radio arrays are in operation: GRANDProto300, in China, GRAND@Auger, in Argentina, and GRAND@Nancay, in France. Their goals are to field-test the design of the radio-detection units, understand the radio background to which they are exposed, and develop tools for diagnostic, data gathering, and data analysis. This list of contributions to the 38th International Cosmic Ray Conference (ICRC 2023) presents an overview of GRAND, in its present and future incarnations, and a look at the first data collected by GRANDProto13, the first phase of GRANDProto300., Comment: Note: To access the list of contributions, please follow the "HTML" link that can be found on the arXiv page (v2: fixed name spelling of one author)

Autonomous radio-detection, i.e., detection of air-showers with standalone radio arrays, is one of the major technical challenges to overcome for the next generation astroparticle detectors. In this context, we study polarisation signatures of simulated radio signals to perform an identification of the associated air-showers initiated by cosmic-rays and neutrinos. We compare the two sources of radio emission (the charge excess and geomagnetic) and show that the former is almost negligible for inclined (zenith angle $>65^{\circ}$) cosmic-ray air-showers. This provides an efficient background rejection criterion at the DAQ level, based on the projection of the total electric field along the direction of the local magnetic field. This relevant quantity can be easily computed -- even in an online treatment -- for antennas measuring three orthogonal polarisations. Independently of the experimental antenna layout, we estimate that assuming a random polarisation of noise events, a rejection from $\approx 72\%$ (for a non favorable detector location) to $\approx 93\%$ (for a favorable location) of the noise induced events and a trigger efficiency of 87% (93%) with a $3\sigma$ ($5\sigma$) trigger threshold level should be achievable. We also show that neutrino-induced showers present a charge excess to geomagnetic signal ratio up to $\sim 10$ times higher than for cosmic ray showers. Although this characteristic makes the identification of neutrino-induced showers challenging via the method developed here, it provides an efficient criterion to perform an offline discrimination between cosmic-ray and neutrino primaries. The stronger charge excess emission will also help the reconstruction of air-shower parameters, such as the core position., Comment: 19 pages, 26 figures, 3 tables. Accepted in Astroparticle Physics; caption formatting corrected in Fig. 1 and Fig. 7

Radio detection of air-showers is a mature technique that has gained momentum over the past decades. With increasingly large-scale experiments, massive air-shower simulations are needed to evaluate the radio signal at each antenna position. Radio Morphing was developed for this purpose. It is a semi-analytical tool that enables a fast computation of the radio signal emitted by any air-shower at any location, from the simulation data of one single reference shower at given positions. It relies on simple electromagnetic scaling laws of the radio emission (i.e., electric field) at the antenna level and then an interpolation of the radio pulse at the desired positions. We present here major improvements on the Radio Morphing method that have been implemented recently. The upgraded version is based on revised and refined scaling laws, derived from physical principles. It also includes shower-to-shower fluctuations and a new spatial interpolation technique, thanks to which an excellent signal timing accuracy of a fraction of nanosecond can be reached. This new implementation, provides simulated signals with relative differences on the peak-to-peak amplitude of ZHAireS simulations below 10\% (respectively 25\%) for 91\% (99\%) of antennas while the computation time was reduced by more than 2 orders of magnitude compared to standard simulations. This makes Radio Morphing an efficient tool that allows for a fast and accurate computation of air-shower radio signals. Further implementation of Askaryan emission or enabling to use an input value of the geomagnetic field should reduce relative differences with ZHAireS by few percents and make the method more universal., Comment: 8 pages, 4 figures. ICRC proceeding

Radio-detection is becoming an established technique for the detection of air showers induced by cosmic particles. This is in particular true at the highest energies, where very large detection areas are required. A proper description of the shape of the radio wavefront emitted by air showers may allow to reconstruct the properties of its parent particle. In this article, we show that for showers with zenith angles larger than 60{\deg} --those targeted by giant radio arrays detecting extensive air showers induced by cosmic particles--, a point-source-like description of the radio wavefront allows to constraint the lateral position of the shower axis within a few meters. Following, we show that the reconstructed longitudinal position of this point source is correlated with the nature of the cosmic rays initiating the shower. Further systematic studies are pending to determine the robustness of this parameter and its validity as a proxy for cosmic ray composition studies., Comment: 24 pages, 13 figures

Very inclined extensive air showers (EAS), with both down-going and up-going trajectories, are particularly targeted by the next generation of extended radio arrays, such as GRAND. Methods to reconstruct the incoming direction, core position, primary energy and composition of showers with these specific geometries, remain to be developed. Towards that goal, we present a new reconstruction procedure based on the arrival times and the amplitudes of the radio signal, measured at each antenna station. This hybrid reconstruction method, harnesses the fact that the emission is observed, at the antenna level, far away from the emission region, thus allowing for a point-like emission description. Thanks to this assumption, the arrival times are modelled following a spherical wavefront emission, which offers the possibility to reconstruct the radio emission zone as a fixed point along the shower axis. From that point the amplitude distribution at the antenna level is described through an Angular Distribution Function (ADF) taking into account at once all geo-magnetic asymmetries and early late effects as well as additional signal asymmetries featured by very inclined EAS. This method shows promising results in terms of arrival direction reconstruction, within the 0.1{\deg} range, even when taking into account experimental uncertainties, and interesting potential for the energy reconstruction and primary composition identification., Comment: 8 pages, ICRC2021 proceeding

The detection of extensive air showers (EAS) through their radio signal is becoming one of the most promising techniques for the study of Neutrinos and Cosmic rays at the highest energies. For the design, optimization and characterization of radio arrays, and of their associated reconstruction algorithms, tens of thousands of Monte Carlo simulations are needed. Current available simulation codes can take several days to compute the signals produced by a single shower, making it impossible to produce the required simulations in a reasonable amount of time, in a cost-effective and environmental-conscious way. In this article we present a method to synthesize the expected signals (the full time trace, not just the peak amplitude) at any point around the shower core, given a set of signals simulated in a finite number of antennas strategically located in a pattern that exploits the signature features of the radio wavefront. The method can be applied indistinctly to the electric field or to the antenna response to the electric field, in the three polarization directions. The synthesized signal can be used to evaluate trigger conditions, compute the fluence or reconstruct the shower incoming direction, allowing for the production of one single library of simulations that can be used and re-used for the characterization and optimization of radio arrays and their associated reconstruction methods, for a thousandth part of the otherwise required CPU time., Comment: 22 pages, 13 figures

Neutrinos of astrophysical origin could be detected through the electromagnetic radiation of the particle showers induced in the atmosphere by their interaction in the Earth. This applies in particular for tau neutrinos of energies $E>10^{16}$eV following Earth-skimming trajectories. The $\sim1^{\circ}$ beaming of the radio emission in the forward direction however implies that the radio signal will likely fly above a detector deployed over a flat site and would therefore not be detected. We study here how a non-flat detector topography can improve the detection probability of these neutrino-induced air showers. We do this by computing with three distinct tools the neutrino detection efficiency for a radio array deployed over a toy-model mountainous terrain, also taking into account experimental and topographic constraints. We show in particular that ground topographies inclined by few degrees only induce detection efficiencies typically three times larger than those obtained for flat areas for favorable trajectories. We conclude that the topography of the area where the detector is deployed will be a key factor for an experiment like GRAND., Comment: 29 pages, 13 figures

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., Comment: published as Astroparticle Physics, Volume 114, 2020, pp. 10-21, 14 pages, 13 figures

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

TREND50 is a radio detection setup of 50 self-triggered antennas working in the 50-100MHz frequency range and deployed in a radio-quiet valley of the Tianshan mountains (China). TREND50 achieved its goal: the autonomous radiodetection and identification of air showers. Thanks to a dedicated offine selection algorithm, 564 air shower candidates were indeed selected out of $7\cdot10^8$ transient radio signals recorded during the 314 live days of data taken during the first two years of operation of this setup (2011 and 2012). This event rate, as well as the distribution of the candidate directions of arrival, is consistent with what is expected from cosmic-ray-induced air showers according to simulations, assuming an additional $\sim$20% contamination of the final sample by background events. This result is obtained at the cost of a reduced air shower detection efficiency, estimated to be $\sim$3%. This low efficiency is mostly due to the large amount of dead time of the setup. This result paves the way for the GRANDProto35 experiment, the first stage of the GRAND project., Comment: 41 pages, 19 figures, published in Astropart. Phys

Over the last decades, radio detection of air showers has been established as a detection technique for ultra-high-energy cosmic-rays impinging on the Earth's atmosphere with energies far beyond LHC energies. Today's second-generation of digital radio-detection experiments, as e.g. AERA or LOFAR, are becoming competitive in comparison to already standard techniques e.g. fluorescence light detection. Thanks to a detailed understanding of the physics of the radio emission in extensive air showers, simulations of the radio signal are already successfully tested and applied in the reconstruction of cosmic rays. However, the limits of the computational power resources are easily reached when it comes to computing electric fields at the numerous positions requested by large or dense antenna arrays. In the case of mountainous areas as e.g. for the GRAND array, where 3D shower simulations are necessary, the problem arises with even stronger acuity. Therefore we developed a full parametrisation of the emitted radio signal on the basis of generic shower simulations which will reduce the simulation time by orders of magnitudes. In this talk we will present this concept after a short introduction to the concept of the radio detection of air-shower induced by cosmic rays., Comment: published in the proceedings of the SF2A 2017 conference

The detection of air-shower events via radio signals requires to develop a trigger algorithm for a clean discrimination between signal and background events in order to reduce the data stream coming from false triggers. In this contribution we will describe an approach to trigger air-shower events on a single-antenna level as well as performing an online reconstruction of the shower parameters using neural networks., Comment: To be published in the proceedings of the ARENA2018 conference, implemented referee's comments

Contributions of the CODALEMA/EXTASIS experiment to the 35th International Cosmic Ray Conference, 12-20 July 2017, Busan, South Korea., Comment: 7 proceedings; will appear in PoS(ICRC2017)

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

After the so-called cosmic recombination, the expanding universe entered into a period of darkness since most of the matter was in a neutral state. About a billion years later, however, the intergalactic space was once again ionized. The process, known as the cosmic reionization, required the operation of mechanisms that are not well understood. Among other ionizing sources, Population III stars, mini-quasars, and X-ray emitting microquasars have been invoked. In this article we propose that primordial cosmic rays, accelerated at the termination points of the jets of the first microquasars, may have contributed to the reionization of the intergalactic space as well. For this we quantify the ionization power of cosmic rays (electrons and protons) in the primordial intergalactic medium using extensive particle cascade simulations. We establish that, depending on the fraction of electrons to protons accelerated in the microquasar jets, cosmic rays should have contributed to the reionization of the primordial intergalactic medium as much as X-rays from microquasar accretion disks. If the primordial magnetic field was of the order of $10^{-17}$ G, as some models suggest, cosmic rays had an important role in ionizing the neutral material far beyond the birth places of the first stars., Comment: 4 pages, 2 figures, accepted for publication as a letter to the editor in Astronomy and Astrophysics

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

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

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