High-energy Multi-messenger Emission from Galaxy Clusters in the Local Universe

Sintesi Tecnica: Emissione Multi-messaggero ad Alta Energia da Ammassi di Galassie nell'Universo Locale

Enunciato del Problema
L'origine dei fondi diffusi di neutrini e raggi $\gamma$ ad alta energia rimane un problema irrisolto di primaria importanza in astrofisica. Sebbene i nuclei galattici attivi (AGN) e le galassie a starburst siano noti contributori, il ruolo specifico degli ammassi di galassie come sorgenti di questi messaggeri cosmici non è pienamente vincolato. Gli ammassi di galassie, caratterizzati da un gas caldo del mezzo intraammasso (ICM) ($\sim 10^8$ K) e da forti campi magnetici ($\sim 1$ $\mu$G), sono teoricamente capaci di accelerare i raggi cosmici (CR) tramite shock di fusione, shock di accrescimento e turbolenza. Questi protoni accelerati possono interagire con l'ICM per produrre raggi $\gamma$ secondari e neutrini attraverso processi adronici. Tuttavia, studi precedenti si sono spesso basati su modelli semplificati a simmetria sferica o hanno assunto l'iniezione di CR esclusivamente dai nuclei centrali degli ammassi, potenzialmente rappresentando in modo errato la distribuzione spaziale e il flusso totale dell'emissione diffusa. Inoltre, distinguere l'emissione diffusa degli ammassi dall'emissione di sorgenti puntiformi (ad esempio da AGN centrali come NGC 1275 o M87) rimane una sfida a causa della risoluzione angolare limitata nelle osservazioni attuali.

Metodologia
Questo lavoro impiega un rigoroso quadro numerico che combina simulazioni Magnetoidrodinamiche (MHD) ad alta risoluzione con metodi di trasporto Monte Carlo per modellare l'emissione multi-messaggero da ammassi di galassie nell'universo locale ($z \lesssim 0.1$).

1. Simulazioni MHD (SLOW): Lo studio utilizza la simulazione "Simulating the LOcal Web" (SLOW), che incorpora condizioni iniziali vincolate osservativamente basate sul catalogo CosmicFLow. A differenza dei modelli precedenti, SLOW fornisce una descrizione 3D auto-consistente e ad alta risoluzione delle distribuzioni di densità del gas, temperatura e campo magnetico all'interno degli ammassi (in particolare Virgo, Perseo e Coma), catturando morfologie non sferiche e sottostrutture derivanti dalle fusioni.
2. Iniezione di Raggi Cosmici: Viene adottato un modello di iniezione motivato fisicamente in cui l'iniezione di CR segue la distribuzione spaziale della densità del gas dell'ammasso. Questo approccio tiene conto dell'accelerazione dei CR nelle regioni ad alta densità associate a turbolenza e shock in tutto il volume dell'ammasso, piuttosto che limitare l'iniezione al nucleo dell'ammasso. Il bilancio energetico totale dei CR è parametrizzato dal rapporto $X_{CR} = E_{CR}/E_{cluster} \simeq 0.01$, coerente con i limiti superiori osservativi.
3. Trasporto dei CR e Produzione Secondaria: La propagazione dei CR (assunti essere protoni) all'interno dell'ICM e del mezzo intergalattico (IGM) è simulata utilizzando il codice Monte Carlo CRPropa. Il codice risolve l'equazione di trasporto di Parker, tenendo conto di avvezione, diffusione spaziale, diffusione di momento e perdite di energia. I canali di interazione chiave includono collisioni protone-protone ($pp$) (dominanti sotto i $10^{17}$ eV), produzione di fotopioni e cascata elettromagnetica (scattering Compton inverso sul CMB e sull'EBL, perdite per sincrotrone). Raggi $\gamma$ e neutrini secondari sono generati in modo auto-consistente da queste interazioni.

Contributi Chiave

• Modellizzazione Realistica dell'Ambiente: L'integrazione delle simulazioni MHD SLOW permette di modellare la propagazione dei CR in ambienti di ammassi realistici e non sferici, superando le ipotesi semplificate di simmetria sferica degli studi analitici o semi-analitici precedenti.
• Schema di Iniezione Esteso: L'adozione di un modello di iniezione di CR spazialmente esteso, che traccia il profilo di densità del gas, rappresenta un significativo miglioramento metodologico rispetto ai lavori precedenti che assumevano un'iniezione esclusivamente centrale. Ciò riduce il tempo di confinamento efficace dei CR e i tassi di interazione nel nucleo, portando a un profilo di emissione più diffuso.
• Stima Completa del Flusso: Lo studio calcola sia i flussi individuali degli ammassi (Virgo, Perseo, Coma) sia il flusso diffuso cumulativo dell'intera popolazione di ammassi massicci ($M \gtrsim 5 \times 10^{13} M_\odot$) entro una distanza comovibile di 500 Mpc.

Risultati

• Singoli Ammassi: I flussi di raggi $\gamma$ diffusi previsti per gli ammassi di Virgo, Perseo e Coma sono ben al di sotto dei limiti superiori attuali riportati dalla collaborazione LHAASO. I flussi previsti sono anche inferiori all'emissione di raggi $\gamma$ osservata da MAGIC dalla sorgente centrale NGC 1275 in Perseo, coerentemente con l'aspettativa che MAGIC rilevi una sorgente puntiforme mentre il modello prevede un'emissione diffusa.
• Flusso Cumulativo: Il flusso di raggi $\gamma$ integrato da tutti gli ammassi entro 500 Mpc è significativamente inferiore al fondo diffuso di raggi $\gamma$ misurato da Fermi-LAT. Lo studio nota che a energie inferiori a $\sim 100$ GeV, altre sorgenti (blazar, AGN disallineati, galassie a starburst) dominano il fondo.
• Emissione di Neutrini: Il flusso di neutrini di tutti i sapori previsto dalla popolazione locale di ammassi è ben al di sotto del flusso diffuso osservato da IceCube, ma è comparabile ai limiti superiori esistenti di IceCube specificamente per gli ammassi di galassie ($M > 10^{14} M_\odot$). Il flusso di neutrini previsto dall'ammasso di Virgo risulta comparabile alla sensibilità proiettata del rivelatore KM3NeT/ARCA.
• Caratteristiche Spettrali: Gli indici spettrali dei flussi previsti ($\alpha \approx 2.0 - 2.5$) sono coerenti con le aspettative per l'accelerazione tramite shock di fusione e accrescimento. I risultati indicano che ammassi di massa inferiore contribuiscono in modo efficiente all'emissione di raggi $\gamma$ a $\sim 10$ GeV, mentre i sistemi massicci dominano il flusso di neutrini ad alta energia ($\gtrsim 10$ TeV) grazie alle loro superiori capacità di confinamento dei CR.

Significato e Affermazioni
Il lavoro afferma che i suoi risultati sono coerenti con i vincoli osservativi attuali, in particolare la non rilevazione di un'emissione diffusa significativa da parte di LHAASO e i limiti superiori stabiliti da IceCube. Gli autori ipotizzano che la tensione tra le previsioni dei modelli precedenti (che spesso assumevano flussi più elevati) e i limiti superiori di IceCube sia sostanzialmente ridotta dall'adozione del modello di iniezione spazialmente esteso e dall'uso di simulazioni SLOW ad alta risoluzione.

Lo studio conclude che, sebbene sia improbabile che gli ammassi di galassie dominino il fondo totale di neutrini ad alta energia diffuso, essi potrebbero rappresentare una frazione non trascurabile dello stesso. Il lavoro fornisce un quadro teorico che motiva ricerche dedicate per neutrini diffusi dall'ammasso di Virgo utilizzando KM3NeT e IceCube-Gen2, e per raggi $\gamma$ diffusi utilizzando CTA e LHAASO. Gli autori sottolineano che i loro risultati si allineano con le prospettive dei futuri osservatori multi-messaggero, offrendo un test tempestivo dei processi di accelerazione adronica dei CR nelle strutture su larga scala.