The sources of ultra-high energy (UHE) cosmic rays, which can have energies up to 10^20 eV, remain a mystery. UHE neutrinos may provide important clues to understanding the nature of cosmic-ray sources. ARIANNA aims to detect UHE neutrinos via radio (Askaryan) emission from particle showers when a neutrino interacts with ice, which is an efficient method for neutrinos with energies between 10^16 eV and 10^20 eV. The ARIANNA radio detectors are located in Antarctic ice just beneath the surface. Neutrino observation requires that radio pulses propagate to the antennas at the surface with minimum distortion by the ice and firn medium. Using the residual hole from the South Pole Ice Core Project, radio pulses were emitted from a transmitter located up to 1.7 km below the snow surface. By measuring these signals with an ARIANNA surface station, the angular and polarization reconstruction abilities are quantified, which are required to measure the direction of the neutrino. After deconvolving the raw signals for the detector response and attenuation from propagation through the ice, the signal pulses show no significant distortion and agree with a reference measurement of the emitter made in an anechoic chamber. Furthermore, the signal pulses reveal no significant birefringence for our tested geometry of mostly vertical ice propagation. The origin of the transmitted radio pulse was measured with an angular resolution of 0.37 degrees indicating that the neutrino direction can be determined with good precision if the polarization of the radio-pulse can be well determined. In the present study we obtained a resolution of the polarization vector of 2.7 degrees. Neither measurement show a significant offset relative to expectation.

The proposed ARIANNA-200 neutrino detector, located at sea-level on the Ross Ice Shelf, Antarctica, consists of 200 autonomous and independent detector stations separated by 1 kilometer in a uniform triangular mesh, and serves as a pathfinder mission for the future IceCube-Gen2 project. The primary science mission of ARIANNA-200 is to search for sources of neutrinos with energies greater than 10^17 eV, complementing the reach of IceCube. An ARIANNA observation of a neutrino source would provide strong insight into the enigmatic sources of cosmic rays. ARIANNA observes the radio emission from high energy neutrino interactions in the Antarctic ice. Among radio based concepts under current investigation, ARIANNA-200 would uniquely survey the vast majority of the southern sky at any instant in time, and an important region of the northern sky, by virtue of its location on the surface of the Ross Ice Shelf in Antarctica. The broad sky coverage is specific to the Moore's Bay site, and makes ARIANNA-200 ideally suited to contribute to the multi-messenger thrust by the US National Science Foundation, Windows on the Universe - Multi-Messenger Astrophysics, providing capabilities to observe explosive sources from unknown directions. The ARIANNA architecture is designed to measure the angular direction to within 3 degrees for every neutrino candidate, which too plays an important role in the pursuit of multi-messenger observations of astrophysical sources.

The detection of the radio emission following a neutrino interaction in ice is a promising technique to obtain significant sensitivities to neutrinos with energies above PeV. The detectable radio emission stems from particle showers in the ice. So far, detector simulations have considered only the radio emission from the primary interaction of the neutrino. For this study, existing simulation tools have been extended to cover secondary interactions from $\mu$- and $\tau$-leptons. We find that secondary interactions of both leptons add up to 25\% to the effective volume of neutrino detectors. Also, muon and tau neutrinos can create several detectable showers, with the result that double signatures do not constitute an exclusive signature for $\tau$-neutrinos. We also find that the background of atmospheric $\mu$-leptons from cosmic rays is non-negligible for in-ice arrays and that an air shower veto should be considered helpful for radio detectors.

The primary mission of the ARIANNA ultra-high energy neutrino telescope is to uncover astrophysical sources of neutrinos with energies greater than $10^{16}\mathrm{eV}$. A pilot array, consisting of seven ARIANNA stations located on the surface of the Ross Ice Shelf in Antarctica, was commissioned in November 2014. We report on the search for astrophysical neutrinos using data collected between November 2014 and February 2019. A straight-forward template matching analysis yielded no neutrino candidates, with a signal efficiency of 79%. We find a 90% confidence upper limit on the diffuse neutrino flux of $E^2\Phi=1.7\times 10^{-6}\mathrm{GeV cm^{-2}s^{-1}sr^{-1}}$ for a decade wide logarithmic bin centered at a neutrino energy of $10^{18}\mathrm{eV}$, which is an order of magnitude improvement compared to the previous limit reported by the ARIANNA collaboration. The ARIANNA stations, including purpose built cosmic-ray stations at the Moore's Bay site and demonstrator stations at the South Pole, have operated reliably. Sustained operation at two distinct sites confirms that the flexible and adaptable architecture can be deployed in any deep ice, radio quiet environment. We show that the scientific capabilities, technical innovations, and logistical requirements of ARIANNA are sufficiently well understood to serve as the basis for large area radio-based neutrino telescope with a wide field-of-view.

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

Ultra high energy neutrinos ($E_\nu > 10^{16.5}$eV$)$ are efficiently measured via radio signals following a neutrino interaction in ice. An antenna placed $\mathcal{O}$(15 m) below the ice surface will measure two signals for the vast majority of events (90% at $E_\nu$=$10^{18}$eV$)$: a direct pulse and a second delayed pulse from a reflection off the ice surface. This allows for a unique identification of neutrinos against backgrounds arriving from above. Furthermore, the time delay between the direct and reflected signal (D'n'R) correlates with the distance to the neutrino interaction vertex, a crucial quantity to determine the neutrino energy. In a simulation study, we derive the relation between time delay and distance and study the corresponding experimental uncertainties in estimating neutrino energies. We find that the resulting contribution to the energy resolution is well below the natural limit set by the unknown inelasticity in the initial neutrino interaction. We present an in-situ measurement that proves the experimental feasibility of this technique. Continuous monitoring of the local snow accumulation in the vicinity of the transmit and receive antennas using this technique provide a precision of $\mathcal{O}$(1 mm) in surface elevation, which is much better than that needed to apply the D'n'R technique to neutrinos.

While the radio detection of cosmic rays has advanced to a standard method in astroparticle physics, the radio detection of neutrinos is just about to start its full bloom. The successes of pilot-arrays have to be accompanied by the development of modern and flexible software tools to ensure rapid progress in reconstruction algorithms and data processing. We present NuRadioReco as such a modern Python-based data analysis tool. It includes a suitable data-structure, a database-implementation of a time-dependent detector, modern browser-based data visualization tools, and fully separated analysis modules. We describe the framework and examples, as well as new reconstruction algorithms to obtain the full three-dimensional electric field from distributed antennas which is needed for high-precision energy reconstruction of particle showers.

Since 2003, significant efforts have been devoted to the understanding of the radio emission of extensive air showers above 20 MHz. Despite some studies led until the early nineties, the band available below 20 MHz has remained unused for 20 years. However, it has been claimed by some pioneering experiments that extensive air showers emit a strong electric field in this band and that there is evidence of a large increase of the radio pulse amplitude with decreasing frequencies. The EXTASIS experiment, located within the Nan\c{c}ay Radioastronomy Observatory and supported by the scintillator array and the autonomous radio stations of the CODALEMA experiment, aims to re-investigate the low-frequency band, and especially to study the so-called "sudden death" contribution, the expected electric field radiated by the shower front when hitting ground level. In this work, we present the instrumental setup of the EXTASIS experiment composed of~$7$ low-frequency antennas operating in [1.7-3.7] MHz and covering approximately \SI{1}{\kilo\metre\squared}. We report the observation of~$18$ air showers detected in coincidence in the three instruments, and estimate a detection threshold of \SI{23\pm4}{\micro\volt\per\metre} from comparisons with detailed SELFAS3 simulations. We also report a strong correlation of the low-frequency signal observation with the atmospheric electric field.

Since 2002, the CODALEMA experiment located within the Nan\c{c}ay Radio Observatory studies the ultra-high energy cosmic rays (above 10^{17} eV) arriving in the Earth atmosphere. These cosmic rays interact with the component of the atmosphere, inducing an extensive air shower (EAS) composed mainly of charged particles (electrons and positrons). During the development of the shower in the atmosphere, these charged particles in movement generate a fast electric field transient (a few nanoseconds to a few tens of ns), detected at ground by CODALEMA with dedicated radio antennas over a wide frequency band (between 20 MHz and 200 MHz). The study of this electric field emitted during the shower development aims to determine the characteristics of the primary cosmic ray which has induced the particle shower: its nature, its arrival direction and its energy. After some theoretical considerations and a short description of the SELFAS simulation code, we will present the CODALEMA experiment, its performances and main results. At last, we will show how the EAS radio-detection technique could be used to observe very high energy gamma rays sources, with the NenuFAR radio telescope.

Since 2003, significant efforts have been devoted to the understanding of the radio emission of extensive air shower in the range [20-200] MHz. Despite some studies led until the early nineties, the [1-10] MHz band has remained unused for 20 years. However, it has been measured by some pioneering experiments that extensive air shower emit a strong electric field in this band and that there is evidence of a large increase in the amplitude of the radio pulse at lower frequencies. The EXTASIS experiment, located within the Nan\c{c}ay Radioastronomy Observatory and supported by the CODALEMA experiment, aims to reinvestigate the [1-10] MHz band, and especially to study the so-called "Sudden Death" contribution, the expected electric field emitted by shower front when hitting the ground level. Currently, EXTASIS has confirmed some results obtained by the pioneering experiments, and tends to bring explanations to the other ones, for instance the role of the underlying atmospheric electric field. Moreover, CODALEMA has demonstrated that in the most commonly used frequency band ([20-80] MHz) the electric field profile of EAS can be well sampled, and contains all the information needed for the reconstruction of EAS: an automatic comparison between the SELFAS3 simulations and data has been developed, allowing us to reconstruct in an almost real time the primary cosmic ray characteristics.

The radio detection of cosmic rays consists in the estimation of the properties of a primary cosmic ray by observing the electric field emitted by the extensive air shower (EAS) created when the primary cosmic ray enters the atmosphere. This technique is fully operative nowadays and presents a good degree of maturity. In addition, several projects intend to employ this technique for the detection of neutrinos. In order for the technique to be useful, accurate methods for computing the electric field created by a particle shower in the context of a particular experiment must exist. Although current ground-based radio experiments lie on the air-soil interface and some planned experiments on the South Pole envision antennas near the air-ice interface, most of the analytical approaches and Monte Carlo codes used for calculating the electric field either do not take into account the effect of the boundary or calculate the radiation fields only (direct, reflected and transmitted radiation fields). When the particle shower and the antenna are close to the boundary, compared to the observation wavelength, the far-field approximation breaks down, which is the case for the low-frequency EXTASIS experiment, for instance. We present in this work a new formula for calculating the exact field emitted by a particle track in two semi-infinite media separated by a planar boundary. We also explore the validity of the far-field approximation and make some predictions for EAS using a simple shower model.

The composition of ultra-high energy cosmic rays is still poorly known and constitutes a very important topic in the field of high-energy astrophysics. Detection of ultra-high energy cosmic rays is carried out via the extensive air showers they create after interacting with the atmosphere constituents. The secondary electrons and positrons within the showers emit a detectable electric field in the kHz-GHz range. It is possible to use this radio signal for the estimation of the atmospheric depth of maximal development of the showers \xmax, with a good accuracy and a duty cycle close to $100\%$. This value of \xmax\ is strongly correlated to the nature of the primary cosmic ray that initiated the shower. We show in this paper the importance of using a realistic atmospheric model in order to correct for systematic errors that can prevent a correct and unbiased estimation of~\xmax.

The radio technique for the detection of high energy cosmic rays consists in measuring the electric field created by the particle showers created inside a medium by the primary cosmic ray. The electric field is then used to infer the properties of the primary particle. Nowadays, the radio technique is a standard, well-established technique. While most current experiments measure the field at frequencies above $20$ MHz, several experiments have reported a large emission at low frequencies, below $10$ MHz. The EXTASIS experiment aims at measuring again and understand this low-frequency electric field. Since at low frequencies the standard far-field approximation for the calculation of the electric field does not necessarily hold, in order to comprehend the low-frequency emission we need to go beyond the far-field approximation. We present in this work a formula for the electric field created by a particle track inside a dielectric medium that is valid for all frequencies. We then implement this formula in the SELFAS Monte Carlo code and calculate the low-frequency electric field of EAS. We also introduce the sudden death pulse (SDP), an emission caused by the coherent deceleration of the shower front, and study its properties.

Contributions of the CODALEMA/EXTASIS experiment to the 35th International Cosmic Ray Conference, 12-20 July 2017, Busan, South Korea.

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.