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

Sintesi Tecnica: Sorgenti Mirage e Grandi Offset tra Halo TeV e Pulsar

Problematica
Recenti osservazioni effettuate dal Large High Altitude Air Shower Observatory (LHAASO) e da HAWC hanno identificato estese emissioni di raggi gamma TeV ("aloni di pulsar") attorno a pulsar di età intermedia. Una sfida significativa nell'interpretazione di queste sorgenti è la frequente osservazione di grandi offset spaziali tra la posizione della pulsar e il centroide dell'emissione TeV. I modelli tradizionali che assumono una diffusione isotropa di coppie elettrone-posizione ($e^\pm$) nel mezzo interstellare (ISM) non possono spiegare naturalmente tali offset, poiché la diffusione isotropa è un'approssimazione valida solo quando la scala fisica supera molteplici lunghezze di coerenza magnetica. Tuttavia, per elettroni a energia ultra-elevata ($>100$ TeV) con tassi di raffreddamento rapidi, la scala di propagazione può essere inferiore alla lunghezza di coerenza del campo magnetico turbolento dell'ISM. Di conseguenza, il trasporto delle particelle può essere altamente asimmetrico, seguendo le linee del campo magnetico locale piuttosto che diffondendosi uniformemente. Questo articolo investiga se tale propagazione asimmetrica possa generare sorgenti "mirage" (sorgenti specchio o illusorie)—picchi di emissione apparenti spostati dall'acceleratore—e significativi offset spaziali.

Metodologia
Gli autori utilizzano una simulazione di particelle di test accelerata su GPU, basata sui primi principi, per tracciare le traiettorie di elettroni da 100 TeV iniettati in un campo magnetico turbolento.

• Generazione del Campo Magnetico: Il campo magnetico turbolento è modellato come una sovrapposizione di quattro campi indipendenti con spettro di Kolmogorov ($n_{in} = -11/3$) generati tramite un metodo nello spazio di Fourier. Questo approccio multi-scala (utilizzando spaziature di griglia da $4 \times 10^{-9}$ a $5 \times 10^{-4}$ kpc) assicura un intervallo inerziale realistico e campi privi di divergenza ($\nabla \cdot B = 0$). Il campo totale include sia componenti regolari ($B_r$) che turbolente ($B_t$).
• Propagazione delle Particelle: Gli elettroni sono propagati utilizzando un algoritmo Boris pusher per risolvere l'equazione del moto relativistica sotto la forza di Lorentz. La simulazione tiene conto delle perdite di energia radiativa tramite radiazione di sincrotrone e scattering Compton inverso (IC) rispetto alla radiazione cosmica di fondo (CMB) e al campo di radiazione interstellare (ISRF).
• Campionamento Statistico: Per campionare efficientemente lo spettro energetico iniziale ($dN/dE \propto E^{-\alpha} \exp(-E/E_{cut})$), gli autori utilizzano un approccio a macro-particelle pesate. Vengono iniettate 8.192 macro-particelle, ciascuna rappresentante un bin energetico logaritmico con un peso proporzionale al numero fisico di elettroni in quel bin.
• Sintesi Osservativa: Le distribuzioni di elettroni simulate vengono convertite in mappe sintetiche di raggi gamma (banda 25–40 TeV) tramite una matrice di scattering IC. Queste mappe sono convolute con una funzione di risposta strumentale (PSF) gaussiana di $\sigma_{PSF} = 0.43^\circ$ (mimando LHAASO-KM2A) e sommate a un modello di rumore di fondo di tipo Poissoniano.
• Identificazione della Sorgente: Viene applicato un framework di confronto di modelli Bayesiani utilizzando il metodo Monte Carlo a catena di Markov (MCMC). I modelli con un numero crescente di sorgenti gaussiane vengono adattati ai dati sintetici. Una sorgente è considerata "identificata" se la statistica di test (TS) aumenta di oltre 25 ($\Delta \log L > 12.5$) aggiungendo una componente.

Contributi Chiave e Risultati
Lo studio esplora sistematicamente lo spazio dei parametri, variando la distanza della pulsar ($d$), l'indice spettrale di iniezione degli elettroni ($\alpha$), l'intensità del campo magnetico totale ($B$), la lunghezza di coerenza della turbolenza ($L_c$) e il rapporto tra campo regolare e turbolento ($B_r/B_t$).

1. Formazione di Sorgenti Mirage: Le simulazioni confermano che la propagazione asimmetrica può creare molteplici sorgenti di raggi gamma distinte ("sorgenti mirage") da un singolo acceleratore. Ciò accade quando gli elettroni si propagano lungo le linee del campo magnetico che, a causa di effetti di proiezione, sembrano originarsi da una posizione distante dal punto di iniezione.
2. Molteplicità della Sorgente ($N$):
   • Il numero di sorgenti identificate è influenzato più fortemente dall'intensità del campo magnetico ($B$). Un $B$ più elevato porta a un raffreddamento più rapido, confinando gli elettroni su scale minori e riducendo la probabilità di formare molteplici sorgenti distinte (N inferiore).
   • Il numero di sorgenti è correlato negativamente con la distanza ($d$) e la lunghezza di coerenza ($L_c$), e positivamente con $B^2$.
3. Separazione delle Sorgenti e Rapporto di Flusso:
   • Separazione: La massima separazione tra le sorgenti identificate è correlata positivamente con $L_c$ e negativamente con $B^2$.
   • Rapporto di Flusso ($\chi$): Il rapporto tra il flusso della sorgente principale e il flusso della sorgente mirage diminuisce all'aumentare di $B^2$. In campi magnetici forti, la sorgente mirage può talvolta essere più luminosa della sorgente principale.
   • Campi Regolari vs Turbolenti: La distanza di separazione mostra una dipendenza non monotona da $B_r/B_t$. La separazione massima si verifica quando $B_r/B_t \approx 1$, dove il campo regolare fornisce una guida direzionale mentre la turbolenza permette uno scattering sufficiente per creare offset. Campi puramente turbolenti o dominanti regolarmente risultano in separazioni minori.
4. Offset Spaziali:
   • Gli offset tra la pulsar e la sorgente identificata più vicina sono massimizzati quando $B_r/B_t \approx 1$ e $L_c$ è grande.
   • In regioni con grandi lunghezze di coerenza (ad esempio, regioni tra i bracci o l'alone galattico), gli offset possono raggiungere i 50 pc.

Significato e Rivendicazioni
L'articolo sostiene di fornire una spiegazione quantitativa per le enigmatiche sorgenti LHAASO in cui gli aloni TeV appaiono spazialmente disconnessi dalle loro pulsar associate. Nello specifico, gli autori suggeriscono che sorgenti quali 1LHAASO J0216+4237u, J0212+4252u e J0206+4302u potrebbero essere aloni mirage alimentati da una singola pulsar a millisecondi (PSR J0218+4232), piuttosto che da acceleratori distinti.

Lo studio sottolinea che l'approssimazione della diffusione isotropa è insufficiente per modellare il trasporto delle particelle vicino agli acceleratori nell'ISM, in particolare per i leptoni ad alta energia dove i tempi di raffreddamento sono più brevi del tempo di coerenza magnetica. Dimostrando che gli effetti di proiezione nei campi magnetici turbolenti possono naturalmente produrre grandi offset e molteplici morfologie di sorgente, il lavoro offre un meccanismo per riconciliare il cielo TeV osservato con i modelli standard di nebulosa di vento di pulsar. Gli autori notano che, sebbene la loro turbolenza sintetica assuma l'isotropia (a differenza della turbolenza alchenica anisotropa della MHD), essa cattura qualitativamente la natura filamentosa su larga scala della turbolenza dell'ISM responsabile di queste emissioni asimmetriche.