A Bayesian Method for Air-Shower Reconstruction using Information Field Theory

Abstract

The radio detection of extensive air showers provides a powerful method for studying the origin of high-energy cosmic rays. The Low-Frequency Array (LOFAR) offers unprecedentedly detailed measurements of the radio emission footprint. However, fully exploiting this information requires advanced reconstruction techniques. In this paper, we introduce a novel framework for air shower reconstruction based on Bayesian inference and Information Field Theory (IFT). Our method is built on a fully differentiable forward model of the radio signal, which incorporates a physical emission parameterization and a precise wavefront model. Additionally, we augment this physical model with Gaussian processes to account for systematic uncertainties in both the signal fluence and arrival timing. By leveraging gradient information, our approach enables efficient (three orders of magnitude acceleration w.r.t.\ the legacy method) and robust inference of the underlying physical shower parameters, such as primary energy and the depth of shower maximum, $X_{\text{max}}$. This work provides not only point estimates but also a rigorous quantification of uncertainties. We achieve a resolution in $X_{\text{max}}$ of $25g/c{\mathrm{m}}^2$ and a radiation energy resolution of $12\%$ on simulations for LOFAR.

19 pages, 11 figures

References in corpus (37)

• LOFAR: The LOw-Frequency ARray
• The Pierre Auger Cosmic Ray Observatory
• Measurements of the Cosmic Ray Composition with Air Shower Experiments
• Simulating radio emission from air showers with CoREAS
• A parsec-scale Galactic 3D dust map out to 1.25 kpc from the Sun
• Information field theory for cosmological perturbation reconstruction and non-linear signal analysis
• Radio detection of cosmic ray air showers in the digital era
• A large light-mass component of cosmic rays at 10^{17} - 10^{17.5} eV from radio observations
• A method for high precision reconstruction of air shower Xmax using two-dimensional radio intensity profiles
• Measurement of the Radiation Energy in the Radio Signal of Extensive Air Showers as a Universal Estimator of Cosmic-Ray Energy
• Detecting cosmic rays with the LOFAR radio telescope
• Bayesian power-spectrum inference for Large Scale Structure data
• General description of electromagnetic radiation processes based on instantaneous charge acceleration in `endpoints'
• NIFTY - Numerical Information Field Theory - a versatile Python library for signal inference
• Coherent Cherenkov Radiation from Cosmic-Ray-Induced Air Showers
• Variable structures in M87* from space, time and frequency resolved interferometry
• Polarized radio emission from extensive air showers measured with LOFAR
• Simulation of Radiation Energy Release in Air Showers
• NuRadioReco: A reconstruction framework for radio neutrino detectors
• Measuring a Cherenkov ring in the radio emission from air showers at 110-190 MHz with LOFAR
• The shape of the radio wavefront of extensive air showers as measured with LOFAR
• Radio measurements of the energy and the depth of the shower maximum of cosmic-ray air showers by Tunka-Rex
• LORA: A scintillator array for LOFAR to measure extensive air showers
• A comparison of the cosmic-ray energy scales of Tunka-133 and KASCADE-Grande via their radio extensions Tunka-Rex and LOPES
• Depth of shower maximum and mass composition of cosmic rays from 50 PeV to 2 EeV measured with the LOFAR radio telescope
• Geometric variational inference
• The wavefront of the radio signal emitted by cosmic ray air showers
• Calibration of the LOFAR low-band antennas using the Galaxy and a model of the signal chain
• An analytic description of the radio emission of air showers based on its emission mechanisms
• Reconstructing air shower parameters with LOFAR using event specific GDAS atmospheres
• Re-Envisioning Numerical Information Field Theory (NIFTy.re): A Library for Gaussian Processes and Variational Inference
• Radio Measurements of the Depth of Air-Shower Maximum at the Pierre Auger Observatory
• Demonstrating Agreement between Radio and Fluorescence Measurements of the Depth of Maximum of Extensive Air Showers at the Pierre Auger Observatory
• Reconstructing non-repeating radio pulses with Information Field Theory
• On the cosmic-ray energy scale of the LOFAR radio telescope
• A high-precision interpolation method for pulsed radio signals from cosmic-ray air showers
• A LOFAR-style reconstruction of cosmic-ray air showers with SKA-Low