Mirage Sources and Large TeV Halo-Pulsar Offsets: Exploring the Parameter Space

Technische Samenvatting: Mirage-bronnen en Grote TeV Halo-Pulsar Offsets

Probleemstelling
Recente observaties door het Large High Altitude Air Shower Observatory (LHAASO) en HAWC hebben uitgebreide TeV gamma-straal emissies ("pulsar halos") geïdentificeerd rond middeloude pulsars. Een aanzienlijke uitdaging bij het interpreteren van deze bronnen is de frequente observatie van grote ruimtelijke offsets tussen de pulsarpositie en het zwaartepunt van de TeV-emissie. Traditionele modellen die uitgaan van isotrope diffusie van elektron-positron ($e^\pm$) paren in het interstellaire medium (ISM), kunnen deze offsets niet op natuurlijke wijze verklaren, aangezien isotrope diffusie alleen een geldende benadering is wanneer de fysieke schaal groter is dan meerdere magnetische coherentielengtes. Voor ultra-hoogenergetische elektronen ($>100$ TeV) met snelle koelingssnelheden kan de propagatieschaal echter kleiner zijn dan de coherentielengte van het turbulente magnetische veld van het ISM. Bijgevolg kan het transport van deeltjes sterk asymmetrisch zijn, waarbij ze de lokale magnetische veldlijnen volgen in plaats van uniform te diffunderen. Dit artikel onderzoekt of dergelijke asymmetrische propagatie "mirage" bronnen (schijnbronnen) kan genereren—schijnbare emissiepieken die verschoven zijn ten opzichte van de versneller—en significante ruimtelijke offsets.

Methodologie
De auteurs maken gebruik van een first-principles, GPU-versnelde testdeeltjessimulatie om de trajecten van 100 TeV elektronen te volgen die worden geïnjecteerd in een turbulent magnetisch veld.

• Generatie van het Magnetisch Veld: Het turbulente magnetische veld wordt gemodelleerd als een superpositie van vier onafhankelijke Kolmogorov-spectrum ($n_{in} = -11/3$) velden gegenereerd via een Fourier-ruimte methode. Deze multi-schaal aanpak (met roosterafstanden van $4 \times 10^{-9}$ tot $5 \times 10^{-4}$ kpc) waarborgt een realistisch inertiebereik en divergentievrije velden ($\nabla \cdot B = 0$). Het totale veld bevat zowel reguliere ($B_r$) als turbulente ($B_t$) componenten.
• Deeltjespropagatie: Elektronen worden gepropageerd met behulp van een Boris pusher algoritme om de relativistische bewegingsvergelijking onder de Lorentz-kracht op te lossen. De simulatie houdt radiatieve energieverliezen rekening via synchrotronstraling en inverse Compton (IC) verstrooiing tegen de kosmische microgolfachtergrond (CMB) en het interstellaire stralingsveld (ISRF).
• Statistische Sampling: Om de initiële energiespectrum ($dN/dE \propto E^{-\alpha} \exp(-E/E_{cut})$) efficiënt te samplen, gebruiken de auteurs een gewogen macro-deeltjesbenadering. 8.192 macro-deeltjes worden geïnjecteerd, waarbij elk een logaritmische energiebin vertegenwoordigt met een gewicht dat proportioneel is aan het fysieke aantal elektronen in die bin.
• Observationele Synthese: De gesimuleerde elektronverdelingen worden omgezet in synthetische gamma-straal kaarten (25–40 TeV band) via een IC-verstrooiingsmatrix. Deze kaarten worden geconvolueerd met een Gaussische Point Spread Function (PSF) van $\sigma_{PSF} = 0.43^\circ$ (om LHAASO-KM2A na te bootsen) en toegevoegd aan een Poissoniaans achtergrondruismodel.
• Bronidentificatie: Een Bayesiaans modelvergelijkingskader met behulp van Markov Chain Monte Carlo (MCMC) wordt toegepast. Modellen met een toenemend aantal Gaussische bronnen worden gefit aan de synthetische data. Een bron wordt als "geïdentificeerd" beschouwd als de Test Statistic (TS) met meer dan 25 ($\Delta \log L > 12.5$) toeneemt bij het toevoegen van een component.

Belangrijkste Bijdragen en Resultaten
De studie verkent systematisch de parameterruimte door de afstand van de pulsar ($d$), de spectrale index van de elektroninjectie ($\alpha$), de totale magnetische veldsterkte ($B$), de turbulentie-coherentielengte ($L_c$) en de ratio van het reguliere naar het turbulente veld ($B_r/B_t$) te variëren.

1. Vorming van Mirage-bronnen: De simulaties bevestigen dat asymmetrische propagatie meerdere onderscheidende gamma-straal bronnen ("mirage sources") kan creëren vanuit een enkele versneller. Dit gebeurt wanneer elektronen langs magnetische veldlijnen propageren die, door projectie-effecten, lijken voort te komen uit een locatie die ver verwijderd is van het injectiepunt.
2. Bron-multipliciteit ($N$):
   • Het aantal geïdentificeerde bronnen wordt het sterkst beïnvloed door de magnetische veldsterkte ($B$). Een hogere $B$ leidt tot snellere koeling, waardoor elektronen beperkt worden tot kleinere schalen en de waarschijnlijkheid van het vormen van meerdere onderscheidende bronnen afneemt (lagere $N$).
   • Het aantal bronnen is negatief gecorreleerd met afstand ($d$) en coherentielengte ($L_c$), en positief gecorreleerd met $B^2$.
3. Bron-scheiding en Fluxratio:
   • Scheiding: De maximale scheiding tussen geïdentificeerde bronnen is positief gecorreleerd met $L_c$ en negatief gecorreleerd met $B^2$.
   • Fluxratio ($\chi$): De ratio van de flux van de hoofdbron naar de flux van de mirage-bron neemt af met een toenemende $B^2$. In sterke magnetische velden kan de mirage-bron soms helderder zijn dan de hoofdbron.
   • Reguliere versus Turbulente Velden: De scheidingsafstand vertoont een niet-monotone afhankelijkheid van $B_r/B_t$. Maximale scheiding treedt op bij $B_r/B_t \approx 1$, waar het reguliere veld directionele begeleiding biedt terwijl turbulentie voldoende verstrooiing mogelijk maakt om offsets te creëren. Zuivere turbulentie of dominante reguliere velden resulteren in kleinere scheidingen.
4. Ruimtelijke Offsets:
   • Offsets tussen de pulsar en de dichtstbijzijnde geïdentificeerde bron zijn maximaal wanneer $B_r/B_t \approx 1$ en $L_c$ groot is.
   • In regio's met grote coherentielengtes (bijv. tussenarmen-regio's of de Galactische halo) kunnen offsets oplopen tot 50 pc.

Betekenis en Claims
Het artikel beweert een kwantitatieve verklaring te bieden voor puzzelende LHAASO-bronnen waarbij TeV-halo's ruimtelijk losgekoppeld lijken van hun geassocieerde pulsars. Specifiek suggereren de auteurs dat bronnen zoals 1LHAASO J0216+4237u, J0212+4254u en J0206+4302u mirage-halo's zouden kunnen zijn die worden aangedreven door een enkele milliseconde pulsar (PSR J0218+4232), in plaats van door afzonderlijke versnellers.

De studie benadrukt dat de isotrope diffusie-benadering onvoldoende is voor het modelleren van deeltjestransport nabij versnellers in het ISM, met name voor hoogenergetische leptonen waarvan de koelingstijden korter zijn dan de magnetische coherentietijd. Door aan te tonen dat projectie-effecten in turbulente magnetische velden van nature grote offsets en meerdere bron-morfologieën kunnen produceren, biedt het papier een mechanisme om de geobserveerde TeV-hemel te verzoenen met standaard pulsar wind nebula-modellen. De auteurs merken op dat hoewel hun synthetische turbulentie isotroop aanneemt (in tegen tegenstelling tot de anisotrope Alfvénische turbulentie in MHD), het kwalitatief gezien de grootschalige filamentaire aard van de ISM-turbulentie die verantwoordelijk is voor deze asymmetrische emissies goed weergeeft.