Energy estimation of cosmic rays by air shower radio signals

Technische Samenvatting: Energie-inschatting van Kosmische Straling door Radio-signalen van Luchtstromen

Probleemstelling
Het bepalen van de primaire energie van kosmische straling is essentieel voor het begrijpen van hun oorsprong, versnelling en voortplantingsmechanismen. Hoewel Uitgebreide Luchtstromen (EAS) geïnitieerd door kosmische straling de directe informatie van het primaire deeltje vernietigen, bieden secundaire deeltjes en radio-emissies reconstructiepaden. Eerdere experimenten (bijv. AERA, Tunka-Rex, LOPES) hebben gebruikgemaakt van de Laterale Distributiefunctie (LDF) van radio-signalen, waarbij de correlatie tussen de signaalintensiteit op specifieke afstanden van de kern van de luchtstroom en de primaire energie wordt onderzocht. Deze methoden vereisen echter vaak het identificeren van een optimale afstand die specifiek is voor de geometrie van elke station, en kunnen gevoelig zijn voor de locatie van de kern van de luchtstroom. Het SURA-experiment, een array van vier radio-antennes in Semnan, Iran, heeft met succes de aankomstrichtingen en de locaties van de kern gereconstrueerd, maar vereist een robuuste methode om de primaire energie te schatten die de gevoeligheid voor variaties in de kernpositie minimaliseert.

Methodologie
De auteurs stellen een reconstructiemethode voor die gebaseerd is op het vergelijken van experimentele radio-signalen met hoogwaardige simulaties met de CoREAS-extensie van CORSIKA. De kern van de methodologie omvat de volgende stappen:

1. Simulatie-opstelling:

   • SURA-configuratie: Simulaties werden uitgevoerd om te voldoen aan de specifieke parameters van het SURA-experiment, inclus€nd de vier Log-Periodieke Dipool-Arrays (LPDA), lokale magnetische veldcomponenten ($B_x = 39,53 \, \mu T$, $B_z = 28,1 \, \mu T$), en een hoogte van 1130 m.
   • Dichte Referentie-array: Een virtuele dichte array van 12.342 antennes met een tussenruimte van 1 meter werd gesimuleerd om een hoogresolutie referentie voor radio-signaalintensiteiten te genereren. Deze array dient als basislijn voor vergelijking.
   • Event Generatie: Proton-geïnduceerde stromen werden gesimuleerd met primaire energieën variërend van $10^{17}$ eV tot $7 \times 10^{17}$ eV. De primaire aankomstrichting werd ingesteld op een zenithhoek van $35^\circ$ en een azimuth van $40^\circ$ (gebaseerd op de piekdetectiegraad van SURA), met aanvullende simulaties voor zenithhoeken van $60^\circ$ en $65^\circ$ om effecten van inclinatie te testen.

2. Data-kalibratie:

   • De studie maakt gebruik van een tweestaps-kalibratieproces: real-time FPGA-filtering (high-pass, band-pass, notch filters) en post-processing (reconstructie van de aankomstrichting, toepassing van frequentierespons, wavelet-ruisonderdrukking en signaal up-sampling naar een resolutie van 0,2 ns).

3. De Schaalfactor ($C_{ij}$) Benadering:

   • Een schaalfactor, $C_{ij}$, wordt gedefinieerd als de ratio van de gemiddelde piekintensiteit van gesimuleerde signalen uit de dichte array ($\bar{A}_{sim}$) tot de piekintensiteit van signalen van de SURA-configuratie ($\bar{A}_{exp}$):
     $$C_{ij} = \frac{\bar{A}_{sim}}{\bar{A}_{exp}}$$
   • Deze factor wordt afgeleid door een $\chi^2$-functie te minimaliseren die de experimentele (gesimuleerde SURA) en de dichte array-signalen over de antenne-array vergelijkt.
   • De methode stelt dat hoewel de absolute signaalintensiteit afhankelijk is van de locatie van de kern, de ratio (schaalfactor) tussen de ijle SURA-array en de dichte referentie-array een specifieke correlatie vertoont met de primaire energie.

Belangrijkste Bijdragen

• Onafhankelijkheid van Kernlocatie: Het artikel toont aan dat de schaalfactor $C_{ij}$ grotendeels onafhankelijk is van de locatie van de kern van de kosmische stralingsstroom. Zelfs voor kernposities ver van het centrum van de array (tot $\pm 30$ m en $\pm 125$ m in verschillende tests), blijft de schaalfactor effectief constant voor een gegeven energie.
• Energierelatie: De auteurs stellen een duidelijke inverse correlatie vast tussen de schaalfactor en de primaire energie. Naarmate de primaire energie toeneemt, neemt de schaalfactor af. Deze relatie wordt gemodelleerd door een inverse eerste-orde functie:
  $$C_{ij} = y_0 + \frac{a}{E}$$
  waarbij $E$ de primaire energie is, genormaliseerd door $10^{17}$ eV, en $y_0$ en $a$ aanpassing-constanten zijn die afhankelijk zijn van de aankomstrichting.
• Robuustheid tegen Aankomstrichting: De methode is getest op geneigde stromen ($60^\circ$ en $65^\circ$ zenithhoeken). De inverse correlatie tussen $C_{ij}$ en energie blijft standhouden bij deze hoeken, hoewel de aanpassingsparameters ($a$ en $y_0$) licht variëren met de zenithhoek.

Resultaten

• Reconstructienauwkeurigheid: Met behulp van de afgeleide schaalfactor en de gefitteerde inverse functie werd de primaire energie van gesimuleerde kosmische straling gereconstrueerd. De methode bereikte een maximale reconstructiefout van ongeveer 11% (geobserveerd bij $5 \times 10^{17}$ eV voor de $35^\circ$ zenithhoek).
• Foutanalyse: De relatieve fout ($\Delta E / E$) werd berekend en bleef zo laag als 2% voor bepaalde energieniveaus (bijv. $2 \times 10^{17}$ eV) en bleef over het geteste energiespectrum over het algemeen binnen de 1%–11%.
• Geneigde Stromen: Voor meer geneigde stromen ($60^\circ$ en $65^\circ$ zenithhoeken) bleef de methode effectief, hoewel de precisie varieerde. De laagste precisie werd waargenomen bij $1 \times 10^{17}$ eV en $5 \times 10^{17}$ eV voor de $65^\circ$ hoek, terwijl de maximale fout voor de $60^\circ$ hoek optrad bij $7 \times 10^{17}$ eV.

Betekenis en Claims
Het artikel beweert dat deze techniek een duidelijk voordeel biedt ten opzichte van eerdere benaderingen door de gevoeligheid voor de locatie van de kern van de luchtstroom te verminderen, een veelvoorkomende bron van onzekerheid bij radio-gebaseerde energie-reconstructie. Door een dichte referentie-array te gebruiken om de ijle experimentele data te normaliseren, biedt de methode een betrouwbare estimator voor de primaire energie zonder dat de exacte bepaling van de optimale afstand van de kern voor elk specifiek event vereist is.

De auteurs concluderen dat de methode geldig is voor zowel verticale als geneigde stromen en kan worden toegepast op de data van het SURA-experiment zodra de werkelijke dataverzameling voltooid is. Zij merken op dat toekomstig werk zal bestaan uit het toepassen van deze techniek op werkelijke SURA-experimentele data en het uitbreiden van simulaties om stromen met nog verder verwijderde kernlocaties en variërende aankomstrichtingen te omvatten om de robuustheid van de methode verder te valideren. De studie claimt niet bestaande methoden te vervangen, maar biedt een complementaire benadering die specifiek is afgestemd op de geometrie en beperkingen van de SURA-array.