Modelling Galactic neutrino emission: contributions from massive star clusters and interstellar cosmic rays

Sintesi Tecnica: Modellazione dell'emissione di neutrini galattici da ammassi stellari massicci e raggi cosmici interstellari

Problema
L'origine dei raggi cosmici (CR) più energetici dell'Universo, in particolare quelli che raggiungono energie nell'ordine dei PeV (PeVatroni), rimane incerta. Sebbene i resti di supernova (SNR) siano tradizionalmente considerati i principali acceleratori, le osservazioni di raggi $\gamma$ suggeriscono che essi fatichino a raggiungere energie di tipo PeV, rendendo necessari altri tipi di classi di sorgenti. Gli Ammassi Stellari Massicci Giovani (YMSC) sono stati proposti come candidati PeVatroni validi grazie all'aumento della turbolenza magnetica generata dallo scontro tra venti stellari e da molteplici eventi di SNR. Parallelamente, l'Osservatorio IceCube ha rilevato un flusso di neutrini ad alta energia dal Piano Galattico. Tuttavia, gli attuali studi si basano su template che tengono conto solo dell'emissione diffusa pura di neutrini generata dall'interazione del "mare" di CR con il gas interstellare. Questi modelli non risolvono esplicitamente il contributo degli acceleratori olistici non risolti, come la popolazione di YMSC, il che potrebbe mascherare le firme specifiche di queste sorgenti.

Metodologia
Gli autori presentano un nuovo quadro teorico per calcolare l'emissione di neutrini attesa dalla Via Lattea, integrando due componenti distinte:

1. Popolazione Sintetica di YMSC: Per superare il limite di osservare gli YMSC solo nel nostro vicinato solare ($\lesssim 2$ kpc), gli autori hanno generato 100 realizzazioni indipendenti della popolazione galattica di YMSC. Questa popolazione sintetica campiona masse degli ammassi ($10^3 - 6,3 \times 10^4 M_\odot$), età ($<30$ Myr) e distribuzioni spaziali (seguendo le braccia a spirale, la Barra Galattica e lo Spur Locale) basandosi su vincoli osservativi.
2. Accelerazione e Trasporto di Particelle: Per ogni ammasso sintetico, il modello simula la popolazione stellare (Sequenza Principale, Wolf-Rayet e progenitori di Supernovae) per calcolare la luminosità meccanica collettiva da venti e supernovae. Gli autori modellano le conseguenti bolle generate dai venti (wind-blown bubbles), calcolando i profili di densità del gas e del campo magnetico.
   • Meccanismi di Accelerazione: L'accelerazione delle particelle è modellata in due siti: lo shock di terminazione del vento collettivo (WTS) e gli shock dei resti di supernova (SNR).
   • Scenari di Diffusione: La propagazione dei CR all'interno delle bolle degli ammassi è computata sotto tre diverse assunzioni di turbolenza: diffusione di Kolmogorov, Kraichnan e Bohm.
   • Produzione di Neutrini: I flussi di neutrini sono calcolati dalle collisioni 있는지 (adroniche) di protoni accelerati con il gas bersaglio (principalmente nella densa calotta accumulata delle bolle).
3. Componente Galattica Diffusa: Il modello calcola anche l'emissione diffusa di neutrini derivante dall'interazione del mare di CR galattici con il gas interstellare. Questo utilizza modelli di propagazione dei CR (incorporando l'indurimento spettrale nell'interno della Galassia) e mappe di gas (GALPROP e un modello alternativo 3D).
4. Generazione di Template: Il flusso totale predetto è la somma del contributo YMSC e della componente diffusa indotta dai CR. Questi vengono sintetizzati in template spaziali e spettrali per il confronto con i dati IceCube.

Contributi Chiave

• Nuovi Template di Neutrini: Il documento fornisce i primi template completi di neutrini che integrano esplicitamente l'emissività della popolazione YMSC insieme al componente diffuso standard indotto dai CR.
• Modellazione della Popolazione Sintetica: Offre una stima statisticamente robusta della produzione cumulativa di neutrini dall'intera popolazione galattica di YMSC, piuttosto che affidarsi solo ai pochi ammassi vicini e risolti.
• Analisi Multi-componente: Il framework distingue tra la produzione "in situ" di neutrini all'interno degli ammassi e l'emissione diffusa di fondo, permettendo una separazione dei contributi delle sorgenti.
• Sensibilità dei Parametri: Lo studio esplora la dipendenza dei flussi predetti dai regimi di diffusione delle particelle (Kolmogorov, Kraichnan, Bohm) e dall'indurimento spettrale dei CR nell'interno della Galassia.

Risultati

• Rilevabilità delle Singole Sorgenti: Si trova che il flusso predetto dalle singole YMSC è marginalmente rilevabile anche da rivelatori su scala chilometro cubico (es. IceCube, KM3NeT). Di conseguenza, l'emissione cumulativa della popolazione YMSC dovrebbe apparire come una componente diffusa non risolta.
• Magnitudo del Flusso: La normalizzazione dei modelli è coerente con il best-fit di IceCube per i template galattici esistenti. Gli autori riscontrano che il contributo non risolto dall'emissione degli ammassi potrebbe essere non trascurabile.
• Dipendenza dalla Diffusione:
  • Nello scenario di diffusione di Bohm (confinamento più efficiente), gli YMSC dominano il flusso di neutrini sotto i $\sim 50$ TeV, potenzialmente sovrappredicendo il flusso rispetto ai dati IceCube negli scenari di Kraichnan e Kolmogorov se non viene assunto alcun indurimento spettrale.
  • Negli scenari Kolmogorov e Kraichnan, la componente diffusa indotta dai CR generalmente domina o fornisce contributi comparabili alla componente YMSC, particolarmente ad energie più elevate.
• Distribuzione Spaziale: Il contributo YMSC è concentrato nell'intervallo di longitudine galattica interna ($|l| \le 30^\circ$, la Cresta Galattica), mentre l'emissione indotta dai CR domina a longitudini maggiori.
• Indurimento Spettrale: L'inclusione dell'indurimento spettrale del mare di CR nell'interno della Galassia riduce il contributo relativo degli YMSC negli scenari Kolmogorov e Kraichnan, portando lo spettro totale predetto in un migliore accordo con le osservazioni di IceCube a $\sim 100$ TeV.
• Coerenza: Le previsioni del modello per la controparte di raggi $\gamma$ delle collisioni adroniche rimangono coerenti con le attuali osservazioni diffuse di raggi $\gamma$ nell'intervallo TeV-PeV.

Significato e Rivendicazioni
Gli autori affermano che il loro lavoro fornisce un'evoluzione necessaria nella modellazione teorica dell'emissione di neutrini galattici. Passando da template puramente diffusi a modelli che includono l'emissività specifica della popolazione YMSC, essi offrono uno strumento per "disaggregare" le due componenti dell'emissione di neutrini galattici. Il documento sostiene che l'applicazione di questi nuovi template ai futuri dataset di neutrini permetterà alla comunità di vincolare il ruolo degli ammassi stellari nell'accelerazione estrema dei CR e nella produzione di neutrini. Gli autori sottolineano che i loro risultati suggeriscono che il contributo non risolto degli ammassi sia probabilmente non trascurabile e che un fitting congiunto dei contributi delle sorgenti e della componente diffusa sia cruciale per identificare gli acceleratori galattici e risolvere l'origine dei raggi cosmici. Il lavoro non pretende di aver identificato definitivamente gli YMSC come la sorgente del segnale di IceCube, ma fornisce il quadro teorico necessario per testare questa ipotesi contro i dati osservativi.