A Multiwavelength Interpretation of HESS J1857+026 Emission Using the Fermi-LAT, VERITAS, and HAWC Observatories

Sintesi Tecnica: Un'interpretazione multi-lunghezza d'onda dell'emissione di HESS J1857+026

Problema
HESS J1857+026 è una sorgente di raggi $\gamma$ a energia molto elevata (VHE) estesa e non identificata, coincidente con l'energetica pulsar PSR J1856+0245. Sebbene la sorgente sia un candidato per una nebulosa di vento di pulsar (PWN), la sua natura multi-lunghezza d'onda rimane complessa. Le osservazioni precedenti hanno rivelato un disallineamento spaziale tra l'estensione dell'emissione nei GeV e nei TeV e una mancanza di controparti X confermate, sollevando interrogativi sul meccanismo di emissione dominante. Nello specifico, non è chiaro se l'emissione sia puramente leptonica (scattering Compton inverso da una PWN), puramente adronica (decadimento del pione da una supernova che interagisce con gas molecolare), o una combinazione di entrambe. Inoltre, le proprietà di diffusione delle particelle nell'ambiente immediato delle PWN evolute rimangono scarsamente vincolate rispetto al mezzo interstellare (ISM) medio.

Metodologia
Questo studio presenta un'analisi completa multi-strumentale di HESS J1857+026 utilizzando 15 anni di dati Fermi-LAT (300 MeV–2 TeV), dati VERITAS (0.3–10 TeV) e nuovi dati HAWC (1–316 TeV). L'analisi procede attraverso diverse fasi:

1. Analisi dei Dati:

   • Fermi-LAT: È stata eseguita un'analisi di verosimiglianza binata (binned likelihood analysis) su dati della classe Pass 8 SOURCE. La regione di interesse (ROI) è stata modellata utilizzando il catalogo 4FGL-DR3 e template diffusi. I test di estensione spaziale hanno confrontato template di punto, disco radiale e Gaussiana radiale.
   • VERITAS: È stata condotta un'analisi di massima verosimiglianza 3D su 30 ore di esposizione utilizzando il framework Gammapy, modellando le sorgenti come leggi di potenza con forme spaziali Gaussiane simmetriche.
   • HAWC: È stata applicata un'analisi di massima verosimiglianza a folding diretto (forward-folding) sul dataset Pass 5 (2860 giorni). Lo studio ha testato molteplici template spaziali (Gaussiano vs. basato sulla diffusione) e forme spettrali (legge di potenza, log-parabola, legge di potenza con cutoff esponenziale) utilizzando i criteri di informazione di Bayes e di Akaike (BIC/AIC) per selezionare i modelli di miglior adattamento.

2. Modellazione Radiativa:

   • Indipendente dal Tempo: Le distribuzioni spettrali di energia (SED) a banda larga sono state adattate utilizzando il pacchetto NAIMA per testare scenari puramente adronici (SNR) e lepto-adronici (PWN/SNR).
   • Dipendente dal Tempo: È stato impiegato un modello evolutivo (simile a Gelfand et al. 2009) per simulare l'espansione e il raffreddamento della PWN all'interno di un guscio di SNR. Una procedura Monte Carlo a catena di Markov (MCMC) ha vincolato i parametri fisici quali l'età del sistema, l'intensità del campo magnetico e gli spettri di iniezione delle particelle.

3. Analisi della Diffusione:

   • Sono stati estratti profili di luminosità superficiale radiale da tutti e tre gli strumenti. Tali profili sono stati adattati con un modello analitico basato sulla diffusione (Equazione 16) per vincolare il coefficiente di diffusione delle particelle, $D(E)$, e la soppressione della diffusione rispetto all'ISM.

Risultati Chiave

• Caratteristiche Spettrali e Spaziali:

  • I dati Fermi-LAT rivelano una sorgente estesa che si adatta meglio a una Gaussiana radiale con un raggio di $0.38^\circ \pm 0.04^\circ$ e un indice spettrale $\Gamma \approx 2.07$. Si osserva un break spettrale, con un indice più morbido ($\Gamma \approx 2.63$) al di sotto di 5 GeV e un indice più duro ($\Gamma \approx 1.99$) al di sopra di 5 GeV.
  • VERITAS rileva la sorgente con una significatività di $15.2\sigma$, misurando un indice spettrale di $\Gamma = 2.48 \pm 0.06$ e una larghezza Gaussiana di $0.236^\circ \pm 0.017^\circ$.
  • I dati HAWC mostrano che la sorgente si estende fino a $\sim 37$ TeV. Il modello spaziale di miglior adattamento favorisce un template di diffusione rispetto a una semplice Gaussiana per un campo magnetico di $1 \mu$G, con un'energia di cutoff $E_c \approx 14$ TeV.
  • La dimensione della sorgente diminuisce all'aumentare dell'energia (da $0.38^\circ$ alle energie dei GeV a $0.26^\circ$ alle energie di HAWC), indicando una morfologia dipendente dall'energia coerente con una PWN evoluta.

• Origine Radiativa:

  • Lepto-adronica vs. Adronica: Sebbene un modello puramente adronico non possa essere interamente escluso per la componente a bassa energia ($<10$ GeV), il modello lepto-adronico fornisce un adattamento superiore ai dati a banda larga. La mancanza di un guscio di SNR rilevato e la bassa densità di gas ambientale ($n_0 \lesssim 0.13$ cm$^{-3}$) derivata dal modello evolutivo scoraggiano un'origine prevalentemente adronica per l'emissione TeV.
  • Evoluzione della PWN: Il modello evolutivo suggerisce che il sistema sia una PWN evoluta con un'età di $\tau = [16, 21]$ kyr e un basso valore di intensità del campo magnetico di $B = [0.4, 1.6] \mu$G. Questi parametri spiegano l'emissione X debole (dovuta al raffreddamento per sincrotrone) e la morfologia dipendente dall'energia. Il basso campo magnetico supporta inoltre la dominanza dello scattering Compton inverso (ICS) rispetto ai campi fotonici locali per l'emissione TeV.

• Diffusione delle Particelle:

  • Il coefficiente di diffusione attorno a HESS J1857+026 risulta essere soppresso rispetto alla media dell'ISM galattico. A 50 TeV, il coefficiente di diffusione è $D_{PWN} \sim 2.14 \times 10^{28}$ cm$^2$ s$^{-1}$ (per $B=1 \mu$G), che è circa due ordini di grandezza inferiore al valore dell'ISM ($D_{ISM} \sim 1.2 \times 10^{30}$ cm$^2$ s$^{-1}$).
  • La soppressione riguarda principalmente la normalizzazione ($D_0$) piuttosto che la dipendenza dalla rigidità (indice spettrale $\delta$), che rimane coerente con i modelli di turbolenza di Kolmogorov ($\delta=1/3$) o Kraichnan ($\delta=1/2$). Ciò supporta lo scenario di una turbolenza magnetica localmente potenziata vicino alla sorgente.

Significatività e Rivendicazioni
L'articolo afferma che i dati combinati multi-lunghezza d'onda forniscono uno strumento potente per vincolare le proprietà di diffusione delle particelle nelle vicinanze delle PWN evolute. La conclusione primaria è che HESS J1857+026 è molto probabilmente una PWN alimentata da PSR J1856+0245, dove l'osservata emissione $\gamma$ è dominata da elettroni relativistici che subiscono scattering Compton inverso.

Lo studio evidenzia che la componente spettrale a bassa energia ($E < 10$ GeV) potrebbe ancora ospitare un contributo adronico da una SNR ospite o vicina, o alternativamente, una popolazione distinta di elettroni all'interno della stessa PWN. Tuttavia, l'emissione ad alta energia è fermamente leptonica. Il fattore di soppressione della diffusione derivato è coerente con altre PWN TeV (es. Geminga, Monogem), rafforzando l'ipotesi che il trasporto di particelle vicino a queste sorgenti sia significativamente ostacolato dalla turbolenza magnetica locale rispetto al generale ISM. Gli autori notano che, sebbene sia stato riportato un eccesso di neutrini a bassa significatività, il grande limite superiore sull'emissione di neutrini rispetto alla regione di emissione $\gamma$ suggerisce un'origine complessa o non correlata per qualsiasi potenziale segnale di neutrini.