Revisiting the Radio Lateral Distribution Function: An amplitude dependence on $X_{\rm max}$ and primary composition

Probleemstelling

De detectie van uitgestrekte luchtschoepen (EAS) veroorzaakt door kosmische straling met radiotechnieken is een gevestigde methode om eigenschappen van de primaire deeltjes (energie, aankomstrichting en massa-samenstelling) te reconstrueren. Hoewel bekend is dat de vorm van de laterale verdelingsfunctie (LDF) van het radiosignaal afhangt van de diepte van het schoepenmaximum ($X_{max}$), was er tot nu toe geen volledig fysiek coherent model dat verklaarde waarom de amplitude van het elektrische veld op grondniveau ook sterk afhankelijk is van $X_{max}$, zelfs na correctie voor het elektromagnetische (EM) energieverschil tussen verschillende schoepen.

Traditioneel werden amplitudeverschillen tussen lichte (bijv. protonen) en zware (bijv. ijzer) primaire deeltjes vaak toegeschreven aan verschillen in de totale EM-energie van de schoep. Echter, recente machine learning-analyses suggereerden dat er een sterke, onafhankelijke amplitudeafhankelijkheid van $X_{max}$ bestaat die informatie bevat over de primaire samenstelling, maar deze werd historisch over het hoofd gezien of als een artefact beschouwd.

Methodologie

De auteurs hebben een uitgebreide simulatiestudie uitgevoerd om de onderliggende fysica van deze amplitudeafhankelijkheid te ontrafelen:

1. Simulaties: Er werden radio-emissiesimulaties uitgevoerd met de ZHAIRES-pakket (een combinatie van ZHS-formalisme en AIRES voor deeltjesontwikkeling).

   • Locaties: Simulaties zijn gedaan voor twee locaties: het Pierre Auger Observatory (Auger) en de Giant Radio Array for Neutrino Detection (GRAND).
   • Parameters: Er werden duizenden proton- en ijzer-geïnduceerde schoepen gesimuleerd met verschillende zenit-hoeken (van 42° tot 85°) en vaste primaire energieën.
   • Normalisatie: Om de invloed van de EM-energie uit te sluiten, werden de gesimuleerde elektrische veldamplitudes genormaliseerd per schoep op basis van hun berekende EM-energie. Hierdoor hebben alle schoepen in de analyse dezelfde EM-energie, en blijven amplitudeverschillen puur gerelateerd aan de geometrie en samenstelling ($X_{max}$).

2. Fysieke Analyse: De auteurs analyseerden de schalingen van het elektrische veld op grondniveau door twee concurrerende effecten te modelleren:

   • Afstandsschaling ($1/R$): Het elektrische veld neemt toe naarmate de afstand ($R$) van het schoepenmaximum ($X_{max}$) tot de grond kleiner wordt (diepere schoepen).
   • Dichtheidsschaling ($1/\rho$): Het geomagnetische emissieproces is omgekeerd evenredig met de luchtdichtheid ($\rho$) op het moment van emissie. Diepere schoepen (hogere $X_{max}$) ontwikkelen zich in dikkere lucht, wat de driftsnelheid van geladen deeltjes verlaagt en de emissie vermindert.
   • Coherentieverlies: Bij hoge zenit-hoeken (lage dichtheid) en sterke geomagnetische velden treden grotere Lorentz-krachten op, wat leidt tot grotere tijdsvertragingen tussen deeltjesbanen. Dit vermindert de coherentie van de emissie. De auteurs introduceerden een niet-lineariteitsfactor $J(\theta)$ om dit verlies van coherentie te kwantificeren, waardoor de dichtheidsschaling wordt: $(1/\rho)J(\theta)$.

Belangrijkste Bijdragen

• Fysisch Model: Het presenteren van een compleet fysiek model dat de amplitude van het radiosignaal verklaart als een competitie tussen de $1/R$-schaling (die diepere schoepen bevoordeelt) en de $(1/\rho)J(\theta)$-schaling (die ondiepere schoepen bevoordeelt, maar gereduceerd wordt door coherentieverlies).
• Identificatie van Coherentieverlies: Het kwantificeren van hoe coherentieverlies, veroorzaakt door de combinatie van lage luchtdichtheid en sterke geomagnetische velden, de dichtheidsschaling kan verzwakken of zelfs omkeren.
• Nieuwe Kader voor Samenstellingsbepaling: Het aantonen dat de amplitude van de LDF direct informatie bevat over de primaire samenstelling, zelfs zonder een expliciete reconstructie van $X_{max}$.

Resultaten

De resultaten tonen een sterk verband tussen de LDF-amplitude en $X_{max}$, afhankelijk van de zenit-hoek en de locatie:

1. Concurrerende Schalingen:

   • Bij lage zenit-hoeken domineert de $1/R$-schaling. Omdat protonen (lichte deeltjes) gemiddeld een dieper $X_{max}$ hebben dan ijzer, vertonen proton-schoepen bij lage hoeken hogere amplitudes.
   • Bij hoge zenit-hoeken wordt het beeld complexer. De $1/\rho$-schaling zou ijzer (die ondieper ontwikkelen) bevoordelen, maar coherentieverlies speelt een grote rol.
   • Op GRAND (sterk geomagnetisch veld) is het coherentieverlies zo groot dat de $1/R$-schaling dominant blijft, zelfs bij hoge hoeken. Protonen hebben hier dus overal hogere amplitudes.
   • Op Auger (zwakker veld) is er een overgangshoek (rond 72°). Beneden deze hoek domineren protonen (hoge amplitude), boven deze hoek domineert ijzer (hoge amplitude), omdat de dichtheidsschaling hier sterker is dan de afstandsschaling.

2. Magic Angle: Rond de "magische hoek" (waar de vorm van de LDF onafhankelijk wordt van $X_{max}$) blijft de amplitude sterk afhankelijk van $X_{max}$. Dit bevestigt dat amplitude-informatie onafhankelijk is van de vorm-informatie.

3. Validatie: De voorspellingen van het schalingsmodel kwamen nauwkeurig overeen met de resultaten van de volledige ZHAIRES-simulaties voor zowel protonen als ijzer over alle onderzochte hoeken en locaties.

Significantie en Implicaties

De bevindingen hebben belangrijke gevolgen voor de radioastronomie van kosmische straling:

• Event-by-Event Samenstellingsbepaling: Omdat de amplitude direct afhankelijk is van $X_{max}$ (en dus de massa), kan de primaire samenstelling van individuele gebeurtenissen worden bepaald zonder eerst $X_{max}$ te reconstrueren. Dit biedt een nieuwe, directe route voor massadiscriminatie, wat eerder werd aangetoond met machine learning maar nu fysisch wordt onderbouwd.
• Verbetering van $X_{max}$-reconstructie: Bestaande methoden (zoals die van LOFAR) die de vorm van de LDF gebruiken, kunnen worden verbeterd door ook de amplitudeafhankelijkheid mee te nemen. Vooral bij hoge zenit-hoeken, waar de vorminformatie afneemt, is de amplitudeinformatie cruciaal.
• Energie-reconstructie: Er bestaat een risico op bias in de schatting van de EM-energie als de amplitudeafhankelijkheid van $X_{max}$ niet correct wordt gecorrigeerd. De huidige empirische correcties voor afstand en dichtheid krijgen hiermee een fysische onderbouwing.
• Locatiekeuze voor Experimenten: Voor experimenten die gericht zijn op zeer schuine schoepen (bijna horizontaal), is de sterkte van het lokale geomagnetische veld minder kritiek dan eerder gedacht, omdat coherentieverlies het netto-effect van het veld op de amplitude bij hoge hoeken neutraliseert.

Kortom, dit werk onthult een fundamentele, maar eerder over het hoofd geziene afhankelijkheid in radio-detectie, die de nauwkeurigheid van reconstructie-algoritmen voor zowel energie als samenstelling aanzienlijk kan verbeteren.