Analysis and implications of the spatio-spectral morphology of the Fermi Bubbles

Sintesi Tecnica: Analisi e Implicazioni della Morfologia Spazio-Spettrale delle Bolle di Fermi

Problema
L'origine e la fonte di energia delle Bolle di Fermi (FBs)—gigantesche strutture di raggi $\gamma$ che si estendono per circa $50^\circ$ sopra e sotto il piano galattico—rimangono incerte. Sebbene la loro morfologia suggerisca un'origine nel Centro Galattico (GC), il meccanismo microfisico che guida l'emissione di raggi $\gamma$ è oggetto di dibattito tra due canali primari: il canale adronico (collisioni inelastiche di protoni di raggi cosmici con il gas ambiente) e il canale leptone (diffusione Compton inversa di fotoni del campo di radiazione interstellare da parte di elettroni di raggi cosmici). Le analisi precedenti si sono basate su rappresentazioni spaziali grossolane (ad esempio, sezioni medie in latitudine), potenzialmente oscurando variazioni spettrali dipendenti dalla posizione. Questo articolo mira a risolvere tali incertezze eseguendo un'analisi spettrale ad alta risoluzione, pixel per pixel, per vincolare le popolazioni di raggi cosmici (CR) genitrici necessarie per riprodurre gli spettri di raggi $\gamma$ osservati.

Metodologia
Gli autori utilizzano una ricostruzione priva di template di dieci anni di dati all-sky di Fermi/LAT forniti da Platz et al. (2023). Nello specifico, impiegano la componente "non-dust" ($I_{nd}$) del modello M2, che isola l'emissione non correlata con le mappe di polvere termica, rivelando chiaramente le FBs rispetto al foreground galattico. Il cielo è pixelizzato utilizzando lo schema HEALPix ($N_{side}=128$), ottenendo circa 20.000 pixel che coprono le FBs.

Per modellare l'emissione, gli autori:

1. Costruiscono una Geometria 3D: Approssimano le FBs come figure di rotazione generate ruotando un confine tracciato a occhio attorno all'asse galattico $z$. Ciò consente il calcolo delle profondità lungo la linea di vista e l'integrazione delle quantità fisiche lungo ogni raggio.
2. Modellano il Campo di Radiazione Interstellare (ISRF): Utilizzano il modello di trasferimento radiativo assiale di Popescu et al. (2017) per determinare l'ISRF dipendente dalla posizione (componenti stellare, della polvere e CMB) lungo ogni linea di vista. Per le regioni oltre la griglia del modello ($z > 10$ kpc), effettuano un'estrapolazione utilizzando un decadimento a legge di potenza o scalatura $z^{-1}$.
3. Forward Modeling: Utilizzando il pacchetto Python Naima, adattano gli spettri di raggi $\gamma$ osservati in ogni pixel contro tre modelli di popolazione genitrice: Power-Law (PL), Exponential Cut-off Power-Law (EPL) e Broken Power-Law (BPL).
4. Adronico vs. Leptone: Per i modelli adronici, assumono una densità numerica uniforme di idrogeno fiduciale di $n_H = 10^{-3} \text{ cm}^{-3}$ e normalizzano i risultati di conseguenza. Per i modelli leptone, calcolano la diffusione Compton inversa (IC) sulla media dell'ISRF derivata per ogni pixel.
5. Goodness of Fit: Valutano i modelli utilizzando il $\chi^2$ ridotto, notando che, a causa delle incertezze correlate nella decomposizione dei dati, questa metrica serve a confrontare le prestazioni relative dei modelli piuttosto che la significatività statistica assoluta.

Risultati Chiave

• Forma Spettrale: I modelli a legge di potenza semplice (PL) sia per protoni che per elettroni forniscono adattamenti scarsi e sono esclusi. Sia i modelli EPL che i modelli BPL forniscono adattamenti ugualmente buoni per entrambi gli scenari adronici e leptone.
• Indici Spettrali e Indurimento: Gli indici spettrali migliori alle basse energie mostrano un evidente indurimento verso la punta meridionale delle bolle sia per i protoni dei raggi cosmici (CRp) che per gli elettroni dei raggi cosmici (CRe). Questa tendenza è coerente con i risultati precedenti ma si estende a latitudini più elevate rispetto alle analisi precedenti.
• Energie di Break/Cutoff:
  • Per i modelli EPL, le energie di cutoff sono dell'ordine di pochi TeV e mostrano poca variazione con la latitudine.
  • Per i modelli BPL, le energie di break sono $\sim 100$ GeV con una variazione latitudinale minima.
  • La variazione dell'indice spettrale attraverso il break nei modelli BPL è $\lesssim 0.5$. Questo è significativamente inferiore al $\Delta \alpha \approx 1$ atteso per un break di raffreddamento dovuto a perdite di sincrotrone o IC, suggerendo che i break spettrali non siano caratteristiche di raffreddamento.
• Contributi del Campo di Radiazione: Nei modelli leptone, la componente della polvere dell'ISRF domina la produzione di raggi $\gamma$ attraverso tutte le energie di Fermi/LAT, seguita dalla componente stellare (che è soppressa dagli effetti di Klein-Nishina ad alte energie) e dalla CMB. Di conseguenza, le perdite IC possono essere approssimate come avvengano nel regime di Thomson in tutte le bolle.
• Densità di Energia e Pressione:
  • Adronico: L'energia totale dei CR è stimata in $\sim 1.7 \times 10^{57} / n_{H,-3}$ erg. La pressione dei CR è probabilmente dominante sulla pressione del gas termico, ma solo di un fattore di poche unità, fortemente dipendente dalla densità di gas target assunta.
  • Leptone: La densità di energia degli CRe aumenta con la distanza dal piano galattico e raggiunge il picco verso i bordi delle bolle, qualitativamente coerente con il brillante bordo (limb brightening).
• Tempi di Raffreddamento (Cooling Timescales):
  • Adronico: Il tempo di raffreddamento è estremamente lungo ($\sim 100$ Gyr per $n_H = 10^{-3} \text{ cm}^{-3}$), implicando che le bolle potrebbero essere molto vecchie o lontane da uno stato stazionario.
  • Leptone: I tempi di raffreddamento per gli elettroni ad alta energia sono estremamente brevi ($\lesssim 1$ Myr ai cappi delle bolle, scendendo a $\sim 0.1$ Myr vicino al GC).

Significatività e Implicazioni
Il lavoro conclude che né i modelli adronici né quelli leptone possono essere definitivamente esclusi basandosi solo sulla forma spettrale, poiché sia i modelli EPL che i BPL si adattano bene ai dati per entrambi i tipi di particelle. Tuttavia, i risultati impongono vincoli rigorosi agli scenari fisici:

1. Trasporto Leptone vs. Accelerazione In-Situ: La combinazione di tempi di raffreddamento brevi e l'osservato aumento della densità di energia degli CRe verso i bordi delle bolle pone una sfida severa ai modelli in cui gli elettroni vengono iniettati nel Centro Galattico e trasportati verso l'esterno. Trasportare gli elettroni ai cappi delle bolle ($\sim 8$ kpc) entro il loro tempo di raffreddamento di $\sim 1$ Myr richiederebbe velocità di espansione ($\sim 8000$ km/s) incoerenti con la fenomenologia dei raggi X. Gli autori suggeriscono che la ri-accelerazione in-situ dei leptoni alla periferia della bolla sia una spiegazione più naturale per la morfologia e la distribuzione della densità di energia osservate.
2. Vincoli di Età Adronica: Il lungo tempo di raffreddamento adronico implica che, se le bolle sono giovani (ad esempio, $< 100$ Myr), il budget energetico richiesto necessiterebbe di un'origine da AGN piuttosto che da formazione stellare. Viceversa, se le bolle sono vecchie (multi-Gyr), un'origine da formazione stellare rimane plausibile, a condizione che la popolazione di CR abbia avuto tempo sufficiente per accumularsi.
3. Indurimento Spettrale: L'indurimento spettrale osservato verso la punta meridionale e la mancanza di un break di raffreddamento ($\Delta \alpha < 0.5$) suggeriscono che le caratteristiche spettrali siano guidate dal meccanismo di iniezione o accelerazione piuttosto che dalle perdite di raffreddamento radiativo.

In sintesi, questo lavoro fornisce i primi vincoli spettrali pixel per pixel sulle Bolle di Fermi, dimostrando che, sebbene le leggi di potenza semplici siano insufficienti, i modelli con cutoff esponenziali o break si adattano bene ai dati. I risultati favoriscono fortemente scenari che implicano o il trasporto relativistico dal nucleo o, più probabilmente, l'accelerazione distribuita in-situ dei raggi cosmici ai bordi delle bolle, particolarmente per il canale leptone.