Constraining the slow-diffusion zone size and electron injection spectral index for the Geminga pulsar halo

Probleemstelling

De studie van gammastraalhalos rond pulsars, zoals die van Geminga, is cruciaal voor het begrijpen van de propagatie van kosmische straling in de Melkweg. Hoewel eerdere metingen (o.a. door HAWC) de diffusiecoëfficiënt van elektronen rondom pulsars hebben bepaald als zijnde ongeveer twee orde van grootte lager dan het galactisch gemiddelde ("slow diffusion"), bleven andere fysieke parameters onzeker. Specifiek zijn de grootte van het traag-diffusiegebied ($r_*$) en het spectrum van de injectie van elektronen uit de pulsar-windnevel (PWN) slecht beperkt. Deze onzekerheden belemmeren een nauwkeurige interpretatie van gerelateerde fenomenen, zoals het overschot aan positronen in de kosmische straling en de diffuse TeV-gammastraal-excess. De vraag is of de waargenomen morfologie en het spectrum van de Geminga-halo kunnen worden verklaard binnen een fysisch realistisch model, en wat de implicaties zijn voor de injectieparameters.

Methodologie

De auteur hanteert een twee-zone diffusiemodel om de propagatie van elektronen rondom de Geminga-pulsar te simuleren.

• Het Model: Het model onderscheidt twee zones:
  1. Een binnenste zone (straal $r_*$) rond de pulsar met een lage diffusiecoëfficiënt ($D_1$), veroorzaakt door turbulente magnetische velden (bijv. van een supernova-remnant of door streaming-instabiliteit).
  2. Een buitenste zone met de standaard interstellaire diffusiecoëfficiënt ($D_{ISM}$).
• Vergelijkingen: De elektronenpropagatie wordt beschreven door de diffusie-verliesvergelijking, inclusief brontermen (injectiespectrum), energieverlies (synchrotron en inverse Compton) en diffusie.
• Injectiespectrum: Het bron-spectrum wordt gemodelleerd als een machtswet met een super-exponentiële cutoff: $q(E) \propto E^{-p} \exp[-(E/E_c)^2]$.
• Data-analyse: De berekeningen worden vergeleken met de nieuwste gegevens van het HAWC-experiment (High-Altitude Water Cherenkov), specifiek:
  • De Surface Brightness Profile (SBP) (ruimtelijke morfologie) om de grootte van het traag-diffusiegebied ($r_*$) te beperken.
  • Het gammastraalspectrum (1–100 TeV) om de spectrale index van de injectie ($p$) en de maximale grootte van $r_*$ te beperken.
• Fitting: Er wordt gebruik gemaakt van een $\chi^2$-test voor de morfologie en een likelihood-ratiotest voor het spectrum om de parameters te fitten en betrouwbaarheidsintervallen (2$\sigma$ en 3$\sigma$) te bepalen.

Belangrijkste Bijdragen

1. Beperking van de traag-diffusiezone: Voor het eerst wordt de grootte van het traag-diffusiegebied ($r_*$) rond Geminga direct beperkt door een combinatie van ruimtelijke morfologie en spectrale fitting, in plaats van alleen op theoretische schattingen te vertrouwen.
2. Correlatie tussen parameters: De studie onthult een sterke anticorrelatie tussen de grootte van het traag-diffusiegebied ($r_*$) en de spectrale index van de injectie ($p$). Een grotere $r_*$ vereist een zachter (steiler) injectiespectrum om het waargenomen harde gammastraalspectrum te verklaren.
3. Koppeling aan positronen: De paper koppelt de uit HAWC-gegevens afgeleide parameters direct aan het voorspelde positronenspectrum op Aarde, specifiek in het kader van het AMS-02-positronen-overschot.

Resultaten

• Grootte van de traag-diffusiezone ($r_*$):
  • Ondergrens: $r_*$ moet groter zijn dan 30 pc. Modellen met $r_* \lesssim 30$ pc worden uitgesloten met een betrouwbaarheid van >95% omdat ze de SBP te snel laten afnemen vergeleken met de HAWC-data.
  • Bovengrens: $r_*$ moet kleiner zijn dan 70 pc. Als $r_* \ge 70$ pc, zou de vereiste spectrale index $p$ kleiner zijn dan 1 om het gammastraalspectrum te fitten. Dit is fysisch onwaarschijnlijk gezien simulaties en waarnemingen van PWN's suggereren dat $p \gtrsim 1$.
  • Conclusie: De traag-diffusiezone heeft een straal van 30–70 pc.
• Injectiespectrale index ($p$):
  • Bij de ondergrens van de zone ($r_* = 30$ pc) is de maximale toegestane spectrale index $p$ 2.17 (bij een 3$\sigma$ betrouwbaarheid).
  • De beste fit voor het spectrum geeft een $p$ van ongeveer 1.19, maar de onzekerheid laat waarden tot 2.17 toe.
  • De cutoff-energie ($E_c$) wordt bepaald op ongeveer 100 TeV.
• Implicaties voor positronen:
  • Wanneer $r_* = 30$ pc en $p = 2.17$ wordt aangenomen, komt het voorspelde positronenspectrum van Geminga op Aarde uitstekend overeen met de AMS-02-metingen in het energiebereik van 50–500 GeV (zowel in flux als in spectrale vorm).
  • Dit suggereert dat Geminga een belangrijke bijdrage kan leveren aan het positronen-overschot, mits de diffusieparameters binnen deze specifieke beperkingen vallen.

Betekenis en Conclusie

Deze studie verscherpt de fysieke beperkingen op de omgeving van de Geminga-pulsar aanzienlijk. Door de combinatie van ruimtelijke en spectrale data van HAWC binnen een tweezonen-model, kan men een reeks van fysisch plausibele parameters identificeren die eerder onbekend waren.

• De bevinding dat $r_*$ tussen 30 en 70 pc ligt, biedt een steviger basis voor het modelleren van de propagatie van kosmische straling in de lokale omgeving van de Melkweg.
• De resultaten tonen aan dat de harde gammastraalspectrum van Geminga niet noodzakelijk een extreem hard injectiespectrum vereist ($p \approx 1$), maar verklaard kan worden met een redelijker $p \approx 2.2$ als de traag-diffusiezone groot genoeg is.
• Cruciaal is dat de parameters die nodig zijn om het gammastraalspectrum te verklaren, ook leiden tot een positronenspectrum dat consistent is met AMS-02. Dit versterkt het argument dat nabije, middel-oude pulsars zoals Geminga de bron zijn van het geobserveerde positronen-overschot.
• De studie benadrukt ook de noodzaak van toekomstige metingen (bijv. door LHAASO) om de degeneratie tussen de diffusie-index ($\delta$), de zonegrootte ($r_*$) en de injectie-index ($p$) volledig te doorbreken.