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

Titolo

Vincoli sulla dimensione della zona a diffusione lenta e sull'indice spettrale di iniezione degli elettroni per l'alone della pulsar Geminga.

1. Il Problema Scientifico

Lo studio dei "halo" di pulsar (alone di raggi gamma) è fondamentale per comprendere la propagazione dei raggi cosmici nella Galassia. In particolare, le osservazioni recenti hanno rivelato che gli elettroni e i positroni ad alta energia che sfuggono dalle nebulose di vento di pulsar (PWNe) diffondono nel mezzo interstellare (ISM) molto più lentamente del previsto (coefficiente di diffusione $D$ circa due ordini di grandezza inferiore alla media galattica).
Questo fenomeno di "diffusione lenta" ha implicazioni cruciali per:

• L'interpretazione dell'eccesso di positroni nei raggi cosmici.
• L'origine dell'eccesso diffuso di raggi gamma TeV.
• La determinazione dei parametri fisici delle PWNe, come lo spettro di iniezione degli elettroni.

Tuttavia, esisteva un'incertezza significativa riguardo alla dimensione spaziale della zona a diffusione lenta ($r^*$) e alla forma esatta dello spettro di iniezione degli elettroni ($p$) dalla pulsar Geminga. Le misurazioni precedenti fornivano solo vincoli sul coefficiente di diffusione, lasciando aperti questi parametri chiave.

2. Metodologia

L'autore utilizza un approccio fenomenologico basato sul modello di diffusione a due zone:

• Zona interna: Una regione sferica di raggio $r^*$ attorno alla pulsar dove il coefficiente di diffusione è soppresso ($D_1$).
• Zona esterna: Il mezzo interstellare standard con un coefficiente di diffusione più alto ($D_{ISM}$).

Il modello risolve l'equazione di trasporto-diffusione con perdite di energia (sinchrotrone e scattering Compton inverso) per gli elettroni/positroni. I parametri chiave variati durante l'analisi sono:

• $r^*$: Dimensione della zona a diffusione lenta.
• $p$: Indice dello spettro di potenza di iniezione ($q(E) \propto E^{-p}$).
• $\delta$: Indice energetico del coefficiente di diffusione ($D \propto E^\delta$).
• $E_c$: Energia di taglio nello spettro di iniezione.

Dati Utilizzati:

• HAWC (High-Altitude Water Cherenkov): Nuove misurazioni aggiornate della morfologia spaziale (Surface Brightness Profile - SBP) e dello spettro energetico dei raggi gamma dell'alone di Geminga.
• AMS-02: Dati sullo spettro dei positroni sulla Terra per verificare la coerenza con le previsioni del modello.

Procedura di Fitting:

1. Vincolo Morfologico: Adattamento del modello alla SBP 1D di HAWC per determinare il limite inferiore di $r^*$.
2. Vincolo Spettrale: Adattamento dello spettro integrato dei raggi gamma (entro 10°) di HAWC per determinare il limite superiore di $p$ e il limite superiore di $r^*$, assumendo vincoli teorici su $p$ ($p \gtrsim 1$).
3. Verifica Positroni: Utilizzo dei parametri ottimali trovati per prevedere il flusso di positroni sulla Terra e confrontarlo con i dati AMS-02.

3. Risultati Chiave

A. Dimensione della Zona a Diffusione Lenta ($r^*$)

• Limite Inferiore: L'analisi della morfologia (SBP) esclude modelli con $r^* \lesssim 30$ pc a un livello di confidenza del 95%. Modelli con $r^*$ troppo piccoli producono un calo troppo rapido del flusso angolare, non compatibile con i dati HAWC.
• Limite Superiore: L'analisi dello spettro energetico, combinata con l'assunzione teorica che l'indice spettrale di iniezione non possa essere inferiore a 1 ($p \gtrsim 1$), impone un limite superiore a $r^* \approx 70$ pc. Per $r^* \ge 70$ pc, il limite superiore di $p$ scende sotto 1, rendendo i modelli non fisici.
• Intervallo Definito: La dimensione della zona a diffusione lenta è vincolata a $30 \text{ pc} \le r^* \le 70 \text{ pc}$.

B. Indice Spettrale di Iniezione ($p$)

• Esiste una correlazione inversa tra $r^*$ e $p$: una zona di diffusione più piccola richiede uno spettro di iniezione più duro (valore di $p$ più basso) per spiegare lo spettro gamma osservato.
• Assumendo il limite inferiore di $r^* = 30$ pc e un'ipotesi di diffusione indipendente dall'energia ($\delta = 0$), il limite superiore di $p$ a 3$\sigma$ è 2.17.
• Se $r^*$ aumenta, il limite superiore per $p$ diminuisce drasticamente (es. per $r^*=70$ pc, $p < 1$).

C. Coerenza con l'Eccesso di Positroni

• Utilizzando i parametri più "ottimistici" ma compatibili con i dati HAWC ($r^* = 30$ pc, $\delta = 0$, $p = 2.17$), il modello predice uno spettro di positroni sulla Terra che coincide perfettamente con le misurazioni di AMS-02 nell'intervallo 50–500 GeV, sia per forma spettrale che per flusso assoluto.
• Questo suggerisce che Geminga potrebbe essere la fonte dominante dell'eccesso di positroni in questo intervallo di energia, a condizione che la zona di diffusione lenta sia vicina al limite inferiore e lo spettro di iniezione sia sufficientemente duro.

4. Contributi e Significatività

1. Superamento della misura del solo coefficiente di diffusione: Il lavoro va oltre la semplice misurazione di $D$, fornendo vincoli diretti sulla geometria della zona di confinamento ($r^*$) e sulla fisica dell'accelerazione ($p$).
2. Raffinamento dei limiti inferiori: Rispetto a studi precedenti (es. Schroer et al. 2023), questo studio esclude con alta significatività ($>95\%$) zone di diffusione più piccole di 30 pc, basandosi su un adattamento diretto ai dati morfologici aggiornati di HAWC.
3. Risoluzione del problema dell'eccesso di positroni: Dimostra che, all'interno del modello a due zone, è possibile conciliare le osservazioni dei raggi gamma TeV (HAWC) con lo spettro dei positroni (AMS-02) senza ricorrere a sorgenti multiple o parametri non fisici, purché $r^*$ sia vicino a 30 pc.
4. Degenerazione dei parametri: L'analisi evidenzia la degenerazione tra $r^*$, $p$ e $\delta$. Mostra che un aumento di $\delta$ (dipendenza energetica della diffusione) renderebbe i vincoli su $p$ ancora più stringenti.
5. Prospettive Future: Il paper sottolinea la necessità di misure morfologiche dipendenti dall'energia (attese da LHAASO) per rompere la degenerazione tra i parametri e verificare le discrepanze con le recenti osservazioni di H.E.S.S. (che mostrano uno spettro più morbido a basse energie).

Conclusione

Il paper stabilisce che la zona a diffusione lenta attorno a Geminga ha una dimensione compresa tra 30 e 70 parsec. Questo vincolo geometrico, combinato con lo spettro di iniezione degli elettroni, permette di spiegare coerentemente sia l'alone di raggi gamma osservato da HAWC sia l'eccesso di positroni misurato da AMS-02, fornendo un quadro fisico più completo della propagazione dei raggi cosmici nelle vicinanze delle pulsar.