The muon charge asymmetry and the directional distribution of thunderstorm events observed by the GRAPES-3 muon telescope

Titel: De muon-ladingasymmetrie en de directionele verdeling van onweersbuien waargenomen door de GRAPES-3 muontelescoop

1. Probleemstelling

De GRAPES-3 muontelescoop (G3MT) in Ooty, India, heeft in de periode april 2011 tot december 2020 duizenden onweersbuien gedetecteerd door variaties in de intensiteit van secundaire kosmische straling (muonen) aan de grond te meten. Een opvallende en tot nu toe onverklaarde observatie was een extreme directionele asymmetrie in de detectie: er werden bijna zes keer zoveel onweersbuien waargenomen vanuit het oosten dan vanuit het westen.

De vraag was of dit een artefact van de instrumentatie was, een lokaal meteorologisch fenomeen, of dat er een fundamenteel fysica-proces aan ten grondslag lag. De auteurs vermoedden dat de ladingssamenstelling van de muonen (de verhouding tussen positieve en negatieve muonen) een cruciale rol speelde, beïnvloed door het aardmagnetisch veld.

2. Methodologie

De studie combineert observatiegegevens met uitgebreide Monte Carlo-simulaties:

• Observaties: Analyse van 487 onweersbuien gedetecteerd door de G3MT. De telescoop bestaat uit 3712 proportionele tellers verdeeld over 16 modules, die muonen in 169 richtingen kunnen reconstrueren (later samengevoegd tot 9 hoofdrichtingen voor analyse). Elektrische veldmeters (EFM's) werden gebruikt om de onweersbuien te valideren en hun potentiaal te schatten.
• Simulaties (CORSIKA): Er werden gedetailleerde simulaties uitgevoerd met de CORSIKA-pakket (Cosmic Ray SImulations for KAscade) om de ontwikkeling van uitgestrekte luchtshowers (EAS) te modelleren.
  • Hadronische interacties: Er werden zeven verschillende combinaties van high- en low-energy hadronische interactiegenerators getest (o.a. SIBYLL, EPOS-LHC, QGSJETIII, DPMJETIII gekoppeld aan FLUKA en UrQMD) om de onzekerheid in de muonproductie te kwantificeren.
  • Variabele parameters: De simulaties varieerden de onweersbui-potentiaal (-3 GV tot +3 GV) en de richting van aankomst om de invloed op de muonintensiteit ($\Delta I_\mu$) en de ladingverhouding ($R_\mu = N_{\mu^+}/N_{\mu^-}$) te bestuderen.
  • Controle-experiment: Er werden simulaties uitgevoerd met een constante ladingverhouding ($R_\mu$) om te testen of de directionele asymmetrie verdween zonder de variatie in $R_\mu$.

3. Belangrijkste Bijdragen en Mechanisme

De kern van de bevindingen ligt in de koppeling tussen het aardmagnetisch veld, de lading van primaire kosmische straling (PCR) en de detectie van onweersbuien:

1. Geomagnetische Filtering (East-West Effect): Het aardmagnetisch veld filtert geladen deeltjes afhankelijk van hun lading en rigiditeit. Omdat PCR's voornamelijk positief geladen zijn, worden deeltjes uit het oosten sterker afgebogen dan die uit het westen. Dit resulteert in een hogere cutoff-rigiditeit in het oosten.
2. Variatie in Muon Ladingverhouding ($R_\mu$): Door deze filtering is de flux van positieve muonen in het oosten relatief hoger dan in het westen. De simulaties tonen aan dat $R_\mu$ systematisch toeneemt van west naar oost:
   • Westelijke richting (WG): $R_\mu \approx 1,14$
   • Oostelijke richting (EG): $R_\mu \approx 1,37$
3. Interactie met Onweersbuien: Onweersbuien hebben doorgaans een positieve potentiaal aan de top. Dit remt positieve muonen af en versnelt negatieve muonen. Omdat er meer positieve muonen zijn ($R_\mu > 1$), resulteert dit in een netto afname van de muonintensiteit ($\Delta I_\mu$).
4. Selectie-effect: Een hogere $R_\mu$ (zoals in het oosten) zorgt voor een grotere verandering in de intensiteit voor dezelfde onweersbui-potentiaal. Hierdoor vallen meer oostelijke gebeurtenissen boven de detectiedrempel ($|\Delta I_\mu| > 0,3\%$) dan westelijke gebeurtenissen.

4. Resultaten

• Directionele Asymmetrie: De data toonde een verhouding van 5,9 ± 0,8 tussen het aantal gebeurtenissen uit het oosten en het westen.
• Simulatiebevestiging: De Monte Carlo-simulaties (met de SIBYLL-FLUKA combinatie) reproduceerden deze asymmetrie met een verhouding van 5,09 ± 0,02.
• Controle met constante $R_\mu$: Wanneer de simulaties werden uitgevoerd met een constante $R_\mu$ (geen directionele variatie), verdween de asymmetrie volledig (verhouding EG/WG = 1,00). Dit bewijst dat de variatie in $R_\mu$ de dominante drijvende kracht is.
• Invloed van Hadronische Generatoren: Verschillende generatorcombinaties leidden tot een spreiding van $R_\mu$ van maximaal 7%. Dit resulteerde in een variatie van maximaal 14% in de afgeleide onweersbui-potentiaal. De SIBYLL-FLUKA combinatie gaf de hoogste $R_\mu$ en dus de meest conservatieve (laagste) schatting van de potentiaal.
• Verdeling: Ongeveer 82% van de gedetecteerde gebeurtenissen kwam uit de oostelijke richtingen (NE, E, SE), terwijl de westelijke richtingen slechts 14% vertegenwoordigden.

5. Significatie

Deze studie levert de eerste kwantitatieve link tussen de door het aardmagnetisch veld veroorzaakte muon-ladingasymmetrie en de directionele verdeling van waargenomen onweersbuien.

• Validatie van G3MT: Het resultaat bevestigt dat de GRAPES-3 telescoop een uniek en betrouwbaar instrument is voor het bestuderen van gigavolt-potentialen in de atmosfeer, en dat de waargenomen asymmetrie geen instrumentale bias is.
• Fysisch Inzicht: Het versterkt het begrip van hoe geomagnetische filtering de detectie van atmosferische elektrische velden beïnvloedt.
• Methodologische Vooruitgang: De studie onderstreept het belang van het rekening houden met de directionele variatie van $R_\mu$ bij het interpreteren van muondata voor atmosferische elektriciteit, en biedt een robuust kader voor het omzetten van intensiteitsvariaties naar elektrische potentialen, rekening houdend met onzekerheden in hadronische interactiemodellen.

Kortom, de schijnbaar vreemde concentratie van onweersbuien in het oosten is een direct gevolg van de interactie tussen het aardmagnetisch veld en de lading van kosmische straling, wat de gevoeligheid van de telescoop voor oostelijke onweersbuien aanzienlijk verhoogt.