Constraining the cosmic-ray mass composition by measuring the shower length with SKA

The current generation of air shower radio arrays has demonstrated that the atmospheric depth of the shower maximum Xmax can be reconstructed with high accuracy. These experiments are now contributing to mass composition studies in the energy range where a transition from galactic to extragalactic cosmic-ray sources is expected. However, we are still far away from an unambiguous interpretation of the data. Here we propose to use radio measurements to derive a new type of constraint on the mass composition, by reconstructing the shower length L. The low-frequency part of the Square Kilometer Array will have an extremely high antenna density of roughly 60.000 antennas within one square kilometer, and is the perfect site for high-resolution studies of air showers. In this contribution, we discuss the impact of being able to reconstruct L, and the unique contribution that SKA can make to cosmic-ray science.Comment: Proceedings 9th International Workshop on Acoustic and Radio EeV Neutrino Detection Activities - ARENA2022, 7-10 June 2022, Santiago de Compostela, Spain (8 pages

Depth of shower maximum and mass composition of cosmic rays from 50 PeV to 2 EeV measured with the LOFAR radio telescope

We present an updated cosmic-ray mass composition analysis in the energy range $10^{16.8}$ to $10^{18.3}$ eV from 334 air showers measured with the LOFAR radio telescope, and selected for minimal bias. In this energy range, the origin of cosmic rays is expected to shift from galactic to extragalactic sources. The analysis is based on an improved method to infer the depth of maximum $X_{\rm max}$ of extensive air showers from radio measurements and air shower simulations. We show results of the average and standard deviation of $X_{\rm max}$ versus primary energy, and analyze the $X_{\rm max}$-dataset at distribution level to estimate the cosmic ray mass composition. Our approach uses an unbinned maximum likelihood analysis, making use of existing parametrizations of $X_{\rm max}$-distributions per element. The analysis has been repeated for three main models of hadronic interactions. Results are consistent with a significant light-mass fraction, at best fit $23$ to $39$ $\%$ protons plus helium, depending on the choice of hadronic interaction model. The fraction of intermediate-mass nuclei dominates. This confirms earlier results from LOFAR, with systematic uncertainties on $X_{\rm max}$ now lowered to 7 to $9$ $\mathrm{g/cm^2}$. We find agreement in mass composition compared to results from Pierre Auger Observatory, within statistical and systematic uncertainties. However, in line with earlier LOFAR results, we find a slightly lower average $X_{\rm max}$. The values are in tension with those found at Pierre Auger Observatory, but agree with results from other cosmic ray observatories based in the Northern hemisphere.Comment: 24 pages, 14 figures. Accepted for publication in Phys. Rev.

A high-precision interpolation method for pulsed radio signals from cosmic-ray air showers

Analysis of radio signals from cosmic-ray induced air showers has been shown to be a reliable method to extract shower parameters such as primary energy and depth of shower maximum. The required detailed air shower simulations take 1 to 3 days of CPU time per shower for a few hundred antennas. With nearly $60,000$ antennas envisioned to be used for air shower studies at the Square Kilometre Array (SKA), simulating all of these would come at unreasonable costs. We present an interpolation algorithm to reconstruct the full pulse time series at any position in the radio footprint, from a set of antennas simulated on a polar grid. Relying on Fourier series representations and cubic splines, it significantly improves on existing linear methods. We show that simulating about 200 antennas is sufficient for high-precision analysis in the SKA era, including e.g. interferometry which relies on accurate pulse shapes and timings. We therefore propose the interpolation algorithm and its implementation as a useful extension of radio simulation codes, to limit computational effort while retaining accuracy.Comment: 19 pages, 12 figures. Submitted for publication in JINST (Journal of Instrumentation

A high-precision interpolation method for pulsed radio signals from cosmic-ray air showers

Analysis of radio signals from cosmic-ray induced air showers has been shown to be a reliable method to extract shower parameters such as primary energy and depth of shower maximum. The required detailed air shower simulations take 1 to 3 days of CPU time per shower for a few hundred antennas. With nearly 60,000 antennas envisioned to be used for air shower studies at the Square Kilometre Array (SKA), simulating all of these would come at unreasonable costs. We present an interpolation algorithm to reconstruct the full pulse time series at any position in the radio footprint, from a set of antennas simulated on a polar grid. Relying on Fourier series representations and cubic splines, it significantly improves on existing linear methods. We show that simulating about 200 antennas is sufficient for high-precision analysis in the SKA era, including e.g. interferometry which relies on accurate pulse shapes and timings. We therefore propose the interpolation algorithm and its implementation as a useful extension of radio simulation codes, to limit computational effort while retaining accuracy

Constraining the cosmic-ray mass composition by measuring the shower length with SKA

The current generation of air shower radio arrays has demonstrated that the atmospheric depth of the shower maximum Xmax can be reconstructed with high accuracy. These experiments are now contributing to mass composition studies in the energy range where a transition from galactic to extragalactic cosmic-ray sources is expected. However, we are still far away from an unambiguous interpretation of the data. Here we propose to use radio measurements to derive a new type of constraint on the mass composition, by reconstructing the shower length L. The low-frequency part of the Square Kilometer Array will have an extremely high antenna density of roughly 60.000 antennas within one square kilometer, and is the perfect site for high-resolution studies of air showers. In this contribution, we discuss the impact of being able to reconstruct L, and the unique contribution that SKA can make to cosmic-ray science.</p

On the cosmic-ray energy scale of the LOFAR radio telescope

Cosmic rays are routinely measured at LOFAR, both with a dense array of antennas and with the LOFAR Radboud air shower Array (LORA) which is an array of plastic scintillators. In this paper, we present two results relating to the cosmic-ray energy scale of LOFAR. First, we present the reconstruction of cosmic-ray energy using radio and particle techniques along with a discussion of the event-by-event and absolute scale uncertainties. The resulting energies reconstructed with each method are shown to be in good agreement, and because the radio-based reconstructed energy has smaller uncertainty on an event-to-event basis, LOFAR analyses will use that technique in the future. Second, we present the radiation energy of air showers measured at LOFAR and demonstrate how radiation energy can be used to compare the energy scales of different experiments. The radiation energy scales quadratically with the electromagnetic energy in an air shower, which can in turn be related to the energy of the primary particle. Once the local magnetic field is accounted for, the radiation energy allows for a direct comparison between the LORA particle-based energy scale and that of the Pierre Auger Observatory. They are shown to agree to within (6$\pm$20)% for a radiation energy of 1 MeV, where the uncertainty on the comparison is dominated by the antenna calibrations of each experiment. This study motivates the development of a portable radio array which will be used to cross-calibrate the energy scales of different experiments using radiation energy and the same antennas, thereby significantly reducing the uncertainty on the comparison

Reconstructing air shower parameters with LOFAR using event specific GDAS atmospheres

The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. However, time-dependent and location-specific atmospheric models are needed for the cosmic ray analysis method used for LOFAR data. There, dedicated simulation sets are used for each detected cosmic ray, to take into account the actual atmospheric conditions at the time of the measurement. Using the Global Data Assimilation System (GDAS), a global atmospheric model, we have implemented time-dependent, realistic atmospheric profiles in CORSIKA and CoREAS. We have produced realistic event-specific atmospheres for all air showers measured with LOFAR, an event set spanning several years and many different weather conditions. A complete re-analysis of our data set shows that for the majority of data, our previous correction factor performed rather well; we found only a small systematic shift of 2 g/cm$^2$ in the reconstructed $X_{\rm max}$. However, under extreme weather conditions, for example, very low air pressure, the shift can be up to 15 g/cm$^2$. We provide a correction formula to determine the shift in $X_{\rm max}$ resulting from a comparison of simulations done using the US-Std atmosphere and the GDAS-based atmosphere.Comment: Accepted for publication in Astroparticle Physics. arXiv admin note: text overlap with arXiv:1911.0285