• ### Modelling uncertainty of the radiation energy emitted by extensive air showers(1708.04830)

Aug. 16, 2017 astro-ph.HE
• ### Towards a fast and precise forward model for air shower radio simulation(1708.02481)

Aug. 8, 2017 astro-ph.IM, astro-ph.HE
The radio detection method for cosmic rays relies on coherent emission from electrons and positrons which is beamed in a narrow cone along the shower axis. Currently the only mod- els to reproduce this emission with sufficient accuracy are Monte Carlo based simulations of the particle and radio emission physics, which require large investments of computation time. The work presented here focuses on condensing the simulation results into a semi-analytical model. This relies on building a framework based on theoretical predictions of radio emission, but instead of calculating the radio signal directly these models are used to map template simu- lations to the specifications of a given radio event. Our current approach slices the radio signal based on atmospheric depth of origin and weights these slices based on a shower parameter such as electron number or an effective dipole moment. One significant gain over the existing Monte Carlo codes lies in the fact this makes the depth of the shower maximum a direct input to the simulation where currently one has to pre-select showers based on their random number seed. Such a model has great potential for heavily simulation-based analysis methods, for example the LOFAR air shower reconstruction. These techniques are severely limited by the available computation time but have the lowest errors in real measurement applications.
• ### 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.
• ### Radio Detection of Cosmic Rays -- Achievements and Future Potential(1704.05348)

April 18, 2017 astro-ph.IM
When modern efforts for radio detection of cosmic rays started about a decade ago, hopes were high but the true potential was unknown. Since then, we have achieved a detailed understanding of the radio emission physics and have consequently succeeded in developing sophisticated detection schemes and analysis approaches. In particular, we have demonstrated that the important air-shower parameters arrival direction, particle energy and depth of shower maximum can be reconstructed reliably from radio measurements, with a precision that is comparable with that of other detection techniques. At the same time, limitations inherent to the radio-emission mechanisms have become apparent. In this article, I shortly review the capabilities of radio detection in the very high-frequency band, and discuss the potential for future application in existing and new experiments for cosmic-ray detection.
• ### Ultimate precision in cosmic-ray radio detection --- the SKA(1608.08869)

Feb. 6, 2017 hep-ex, astro-ph.IM, astro-ph.HE
As of 2023, the low-frequency part of the Square Kilometre Array will go online in Australia. It will constitute the largest and most powerful low-frequency radio-astronomical observatory to date, and will facilitate a rich science programme in astronomy and astrophysics. With modest engineering changes, it will also be able to measure cosmic rays via the radio emission from extensive air showers. The extreme antenna density and the homogeneous coverage provided by more than 60,000 antennas within an area of one km$^2$ will push radio detection of cosmic rays in the energy range around 10$^{17}$ eV to ultimate precision, with superior capabilities in the reconstruction of arrival direction, energy, and an expected depth-of-shower-maximum resolution of 6~g/cm${^2}$.
• ### Initial simulation study on high-precision radio measurements of the depth of shower maximum with SKA1-low(1702.00283)

Feb. 1, 2017 astro-ph.IM
As LOFAR has shown, using a dense array of radio antennas for detecting extensive air showers initiated by cosmic rays in the Earth's atmosphere makes it possible to measure the depth of shower maximum for individual showers with a statistical uncertainty less than $20\,g/cm^2$. This allows detailed studies of the mass composition in the energy region around $10^{17}\,eV$ where the transition from a Galactic to an Extragalactic origin could occur. Since SKA1-low will provide a much denser and very homogeneous antenna array with a large bandwidth of $50-350\,MHz$ it is expected to reach an uncertainty on the $X_{\max}$ reconstruction of less than $10\,g/cm^2$. We present first results of a simulation study with focus on the potential to reconstruct the depth of shower maximum for individual showers to be measured with SKA1-low. In addition, possible influences of various parameters such as the numbers of antennas included in the analysis or the considered frequency bandwidth will be discussed.
• ### Radiowave Detection of Ultra-High Energy Neutrinos and Cosmic Rays(1701.02987)

Jan. 11, 2017 astro-ph.IM
• ### Radio detection of extensive air showers(1701.02995)

Jan. 11, 2017 hep-ex, astro-ph.IM
Radio detection of extensive air showers initiated in the Earth's atmosphere has made tremendous progress in the last decade. Today, radio detection is routinely used in several cosmic-ray observatories. The physics of the radio emission in air showers is well-understood, and analysis techniques have been developed to determine the arrival direction, the energy and an estimate for the mass of the primary particle from the radio measurements. The achieved resolutions are competitive with those of more traditional techniques. In this article, I shortly review the most important achievements and discuss the potential for future applications.
• ### Simulation of the Radiation Energy Release in Air Showers(1609.05743)

Sept. 19, 2016 astro-ph.HE
A simulation study of the energy released by extensive air showers in the form of MHz radiation is performed using the CoREAS simulation code. We develop an efficient method to extract this radiation energy from air-shower simulations. We determine the longitudinal profile of the radiation energy release and compare it to the longitudinal profile of the energy deposit by the electromagnetic component of the air shower. We find that the radiation energy corrected for the geometric dependence of the geomagnetic emission scales quadratically with the energy in the electromagnetic component of the air shower with a second order dependency on the atmospheric density at the position of the maximum of the shower development $X_\mathrm{max}$. In a measurement where $X_\mathrm{max}$ is not accessible, this second order dependence can be approximated using the zenith angle of the incoming direction of the air shower with only a minor deterioration in accuracy. This method results in an intrinsic uncertainty of 4% with respect to the electromagnetic shower energy which is well below current experimental uncertainties.
• ### Simulation of Radiation Energy Release in Air Showers(1606.01641)

Sept. 19, 2016 hep-ph, astro-ph.HE
A simulation study of the energy released by extensive air showers in the form of MHz radiation is performed using the CoREAS simulation code. We develop an efficient method to extract this radiation energy from air-shower simulations. We determine the longitudinal profile of the radiation energy release and compare it to the longitudinal profile of the energy deposit by the electromagnetic component of the air shower. We find that the radiation energy corrected for the geometric dependence of the geomagnetic emission scales quadratically with the energy in the electromagnetic component of the air shower with a second-order dependence on the atmospheric density at the position of the maximum shower development $X_\mathrm{max}$. In a measurement where $X_\mathrm{max}$ is not accessible, this second order dependence can be approximated using the zenith angle of the incoming direction of the air shower with only a minor loss in accuracy. Our method results in an intrinsic uncertainty of 4% in the determination of the energy in the electromagnetic air-shower component, which is well below current experimental uncertainties.
• ### Nanosecond-level time synchronization of AERA using a beacon reference transmitter and commercial airplanes(1609.00630)

Sept. 2, 2016 physics.ins-det
Radio detection of cosmic-ray air showers requires time synchronization of detectors on a nanosecond level, especially for advanced reconstruction algorithms based on the wavefront curvature and for interferometric analysis approaches. At the Auger Engineering Radio Array, the distributed, autonomous detector stations are time-synchronized via the Global Positioning System which, however, does not provide sufficient timing accuracy. We thus employ a dedicated beacon reference transmitter to correct for event-by-event clock drifts in our offline data analysis. In an independent cross-check of this "beacon correction" using radio pulses emitted by commercial airplanes, we have shown that the combined timing accuracy of the two methods is better than 2 nanoseconds.
• ### 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.
• 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.
• ### A parameterization for the radio emission of air showers as predicted by CoREAS simulations and applied to LOFAR measurements(1402.2872)

June 19, 2014 astro-ph.IM, astro-ph.HE
Measuring radio emission from air showers provides excellent opportunities to directly measure all air shower properties, including the shower development. To exploit this in large-scale experiments, a simple and analytic parameterization of the distribution of the radio signal at ground level is needed. Data taken with the Low-Frequency Array (LOFAR) show a complex two-dimensional pattern of pulse powers, which is sensitive to the shower geometry. Earlier parameterizations of the lateral signal distribution have proven insufficient to describe these data. In this article, we present a parameterization derived from air-shower simulations. We are able to fit the two-dimensional distribution with a double Gaussian, requiring five fit parameters. All parameters show strong correlations with air shower properties, such as the energy of the shower, the arrival direction, and the shower maximum. We successfully apply the parameterization to data taken with LOFAR and discuss implications for air shower experiments.
• ### The renaissance of radio detection of cosmic rays(1310.6927)

Oct. 25, 2013 astro-ph.IM, astro-ph.HE
Nearly 50 years ago, the first radio signals from cosmic ray air showers were detected. After many successful studies, however, research ceased not even 10 years later. Only a decade ago, the field was revived with the application of powerful digital signal processing techniques. Since then, the detection technique has matured, and we are now in a phase of transition from small-scale experiments accessing energies below 1018 eV to experiments with a reach for energies beyond 1019 eV. We have demonstrated that air shower radio signals carry information on both the energy and the mass of the primary particle, and current experiments are in the process of quantifying the precision with which this information can be accessed. All of this rests on a solid understanding of the radio emission processes which can be interpreted as a coherent superposition of geomagnetic emission, Askaryan charge-excess radiation, and Cherenkov-like coherence effects arising in the density gradient of the atmosphere. In this article, I highlight the "state of the art" of radio detection of cosmic rays and briefly discuss its perspectives for the next few years.
• ### Full Monte Carlo simulations of radio emission from extensive air showers with CoREAS(1307.7566)

July 29, 2013 astro-ph.IM, astro-ph.HE
CoREAS is a Monte Carlo simulation code for the calculation of radio emission from extensive air showers. It is based on the "endpoint formalism" for radiation from moving charges implemented directly in CORSIKA. Consequently, the full complexity of the air-shower physics is taken into account without the need for approximations or assumptions on the emission mechanism. We present results of simulations for an unthinned shower performed with CoREAS for both MHz and GHz frequencies. At MHz frequencies, the simulations predict the well-known mixture of geomagnetic and charge excess radiation. At GHz frequencies, the emission is strongly influenced by Cherenkov effects arising from the varying refractive index in the atmosphere. In addition, a qualitative difference in the symmetry of the GHz radiation pattern is observed when compared to the ones at lower frequencies. We also discuss the strong increase in the ground area subtended by the radio emission when going from near-vertical to very inclined geometries, making very inclined air showers the most promising ones for cosmic ray radio detection.
• Contributions of the Pierre Auger Collaboration to the 33rd International Cosmic Ray Conference, Rio de Janeiro, Brazil, July 2013
• ### A detailed comparison of REAS3 and MGMR simulations for radio emission from EAS(1202.3352)

Feb. 15, 2012 astro-ph.HE
In the last years, several models and simulations calculating the radio emission from cosmic ray air showers have been developed. However, a number of those made conflicting predictions on the pulse shapes and the amplitudes of the radio signal. In the scope of this paper, we discuss a detailed comparison of two independent and complementary theoretical approaches, namely MGMR and REAS3. Furthermore, we study the influence of the underlying air shower models on the predicted pulse shapes and amplitudes and show that remaining discrepancies between MGMR and REAS3 are mostly determined by the air shower models. With this general agreement, a breakthrough in the understanding of the modelling of radio emission from air showers has been achieved.
• ### The endpoint formalism for the calculation of electromagnetic radiation and its applications in astroparticle physics(1112.2126)

Dec. 9, 2011 astro-ph.HE
Over the past years, the freely available Monte Carlo-code REAS which simulates radio emission from air showers based on the geosynchrotron model, was used regularly for comparisons with data. However, it emerged that in the previous version of the code, emission due to the variation of the number of charged particles within an air shower was not taken into account. In the following article, we show the implementation of these emission contributions in REAS3 by the inclusion of end-point contributions'' and discuss the changes on the predictions of REAS obtained by this revision. The basis for describing radiation processes is an universal description which is gained by the use of the end-point formulation. Hence, not only pure geomagnetic radiation is simulated with REAS3 but also radiation due to the variation of the net charge excess in the air shower, independent of the Earth's magnetic field. Furthermore, we present a comparison of lateral distributions of LOPES data with REAS3-simulated distributions. The comparison shows a good argeement between both, data and REAS3 simulations.