Search

The LOFAR radio telescope measures air showers in the energy range 1017 to 1018 eV. For each measured shower, the depth of shower maximum Xmax is reconstructed by simulating the radio signal for an ensemble of showers using Corsika and CoREAS. Fitting their radio ‘footprints’ on the ground to the measured radio data yields an Xmax estimate to a precision of about 20 g/cm2 . Compared to previous works, we have improved the method in several ways. Local atmospheric data and refractive index profiles are now included into the simulations. The energy estimate and the fitting procedure are now done using the radio signals only, thus limiting systematic uncertainties due to the particle detector array (LORA). Using selection criteria from a more elaborate characterisation of the radio and particle detection, we reduce a composition bias in the Xmax reconstruction. A possible residual bias has been bounded from above. Thus, the systematic uncertainties on hXmaxi have been lowered, reducing an important limiting factor for composition studies at any level of statistics.

The LOw-Frequency ARray (LOFAR) measures radio emission from extensive air showers. Precise knowledge of the electric field at each antenna in the 30 − 80 MHz range is obtained using a newly developed, frequency-dependent calibration built on knowledge of the Galactic emission and a detailed model of the signal chain. The energy fluence for each event is then determined, allowing for the calculation of the radiation energy of the air shower. The radiation energy, corrected for geometrical effects, scales quadratically with the energy contained in the electromagnetic component of the air shower. These measurements, combined with predictions that rely only on first-principle electrodynamics, provide an energy estimate for the primary particle. In this contribution we present the radio-based energy scale of cosmic rays detected with LOFAR, and compare it to particle-based energy measurements made using the scintillator array located at the LOFAR core.

Estimating the depth of shower maximum Xmax with high precision is of great interest for the study of primary particle composition. One of the systematic uncertainties in reconstructing Xmax from the radio emission of air showers is the limited knowledge of the atmospheric parameters like humidity, pressure, temperature and the index-of refraction. Using the Global Data Assimilation System (GDAS), a global atmospheric model, we have implemented time-dependent realistic atmospheric profiles in the air shower simulation codes CORSIKA and the radio plug-in CoREAS. This program is now available within CORSIKA and flexible to be adapted for different air shower experiments. We have analyzed the LOFAR cosmic ray data with dedicated simulations for each detected air shower with event specific GDAS atmospheres and investigated the effects of pressure and humidity on the reconstructed Xmax. This study shows that for bulk of the events, where the ground pressure is close to US standard atmosphere values, there is a small systematic shift in Xmax that is less than 2 g/cm2 and for very low pressure values the shift is up to 15 g/cm2 .

The detection of astrophysical neutrinos by IceCube has opened a new window for the observation of our universe. The use of antennas for detecting astroparticles has been proven useful and accurate after many years of radio cosmic-ray detections, and will be a valuable asset for increasing the statistics for the detection of neutrinos < 10PeV and reaching out to the highest energies. The radio technique requires realistic and precise simulations of the electric field created after the neutrino interacts with the detector volume, how it propagates through media, and what measurable signal it produces in the detector. With the goal of providing accurate, fast simulations of neutrino-induced electric fields, we have created the NuRadioMC code. It builds on the experience of software used for current radio neutrino and cosmic ray experiments.