Nuclei in high-energy neutrino sources: A multimessenger study of in-source propagation

Sintesi Tecnica: Nuclei nelle sorgenti di neutrini ad alta energia: uno studio multimessaggero della propagazione in situ

Problematica
La scoperta di neutrini astrofisici ad alta energia da parte di IceCube ha sollevato questioni critiche sulla natura delle loro sorgenti. Mentre i modelli standard assumono che questi neutrini originino da interazioni protone-protone ($pp$) o protone-fotone ($p\gamma$), la mancanza di controparti gamma corrispondenti nel range TeV–GeV suggerisce che queste sorgenti siano "nascoste" o opache ai fotoni ad alta energia. Gli studi precedenti hanno ampiamente assunto una composizione puramente protonica delle particelle accelerate. Tuttavia, la composizione chimica delle particelle iniettate (ad esempio, nuclei leggeri vs. pesanti) rimane una questione aperta. Se la sorgente contiene una frazione significativa di nuclei, la propagazione in situ comporta cascate nucleari complesse (fotodisintegrazione) e processi elettromagnetici che alterano le conseguenti rese di neutrini e raggi gamma. Il problema centrale affrontato è come il rilassare l'ipotesi del protone puro e considerare la composizione nucleare influenzi la fenomenologia multimessaggera, specificamente il rapporto tra i flussi di neutrini e raggi gamma e la luminosità della sorgente richiesta.

Metodologia
Gli autori eseguono simulazioni dettagliate Monte Carlo di cascate nucleari ed elettromagnetiche all'interno di una sorgente astrofisica ad alta energia, utilizzando la galassia Seyfert NGC 1068 come caso di studio rappresentativo. La simulazione è divisa in due fasi:

1. Interazioni fotonucleari e cascata nucleare: Il codice traccia l'evoluzione dei nuclei primari (da $Z=1$ a $Z=90$) attraverso un campo di fotoni di fondo. Include:
   • Fotodisintegrazione (PD): Modellata utilizzando le sezioni d'urto TALYS/TENDL-2023 per nuclei pesanti e parametrizzazioni dedicate per nuclei leggeri ($A \le 9$). Questa include canali di rottura esclusivi (ad esempio, $(\gamma, n)$, $(\gamma, p)$, $(\gamma, \alpha)$) e fluorescenza di risonanza nucleare.
   • Produzione di fotomesoni (PM): Modellata utilizzando il generatore di eventi SOPHIA con un'approssimazione di superposizione per i nuclei.
   • Produzione di coppie Bethe-Heitler (BH): Trattata come un processo di perdita di energia continua per i nuclei carichi.
   • Produzione di coppie con cattura (PPC): Una nuova inclusione della produzione di coppie in cui l'elettrone viene catturato nel guscio L (PPCL), seguito da deeccitazione atomica immediata con emissione di fotoni di linea.
   • Decadimento nucleare: Incluso per completezza, ma ritenuto trascurabile per gli intervalli temporali considerati.
2. Cascata elettromagnetica: I fotoni secondari, gli elettroni e i positroni derivanti dalla fase (i) sono propagati attraverso il campo di radiazione. Ciò include lo scattering Compton inverso (ICS) standard e la produzione di coppie elettroniche (EPP), nonché processi rari come la doppia produzione di coppie (DPP) e la produzione di coppie muoniche (MPP), sebbene quest'ultima sia annotata come trascurabile negli intervalli di energia studiati.

La sorgente è modellata come una sfera omogenea di raggio $R$ immersa in un campo di fotoni disco-corona (disco UV + corona a raggi X). Lo studio varia la dimensione della sorgente ($R$) e la composizione nucleare iniettata (protoni, neutroni e vari nuclei come Carbonio, Ferro e Zirconio) per mappare le profondità ottiche e la ripartizione dell'energia.

Contributi Chiave

• Framework Numerico: Il documento presenta un nuovo trattamento numerico delle cascate indotte dai nuclei che estende i processi fotonucleari fino agli isotopi $Z=90$ e tiene esplicitamente conto della produzione di coppie con cattura nel guscio L, un processo precedentemente trascurato in tali contesti.
• Analisi Sistematica della Composizione: Fornisce un'analisi completa di come la composizione nucleare iniettata alteri il rapporto di resa neutrini-raggi gamma, dimostrando che questa relazione è non monotona e altamente dipendente dalla profondità ottica della sorgente (densità di fotoni e dimensione).
• Rianalisi COMPTEL: Gli autori eseguono una rianalisi dedicata dei dati d'archivio COMPTEL per NGC 1068 utilizzando framework moderni (GammaLib/ctools), derivando limiti aggiornati di confidenza al 95% nella banda MeV che sono circa un fattore 2 più stretti rispetto ai limiti precedenti nei bin di energia inferiori.

Risultati
Lo studio identifica due regimi distinti basati sulla compattezza della sorgente (densità di fotoni):

1. Regime Compatto/Denso ($R \lesssim 2 R_S$):

   • La sorgente è otticamente spessa alla fotodisintegrazione. I nuclei pesanti vengono rapidamente frammentati in nucleoni liberi (protoni e neutroni) prima di uscire.
   • Una frazione significativa dell'energia viene trasferita a neutroni liberi. Poiché i neutroni non subiscono perdite di energia Bethe-Heitler (BH), l'energia disponibile per la produzione di neutrini viene preservata.
   • Risultato: L'iniezione di nuclei pesanti è più efficiente nel produrre neutrini per unità di luminosità dei raggi cosmici rispetto ai protoni. Di conseguenza, far corrispondere il flusso di neutrini osservato da IceCube richiede una luminosità dei raggi cosmici totale minore per i nuclei pesanti. La cascata elettromagnetica risultante (raggi gamma MeV–GeV) è pi più debole (del ~28–35%) rispetto al caso protonico, riducendo così i vincoli dalle osservazioni dei raggi gamma.

2. Regime Esteso/Raffollato ($R \gtrsim 5 R_S$):

   • La fotodisintegrazione è incompleta; una frazione sostanziale dell'energia iniettata rimane bloccata nei frammenti nucleari carichi durante l'attraversamento della sorgente.
   • Questi frammenti pesanti sperimentano perdite di energia BH potenziate (che scalano come $Z^2/A$) rispetto ai protoni alla stessa energia per nucleone.
   • Risultato: L'iniezione di nuclei pesanti è meno efficiente nel produrre neutrini. Per far corrispondere il flusso IceCube, è necessaria una luminosità dei raggi cosmici maggiore (un aumento del 20–30%). Inoltre, le perdite BH potenziate e la maggiore luminosità portano a una cascata elettromagnetica significativamente più forte (flusso di raggi gamma un fattore ~2.0–2.5 superiore), causando un superamento dei limiti superiori esistenti per i raggi gamma più severo rispetto ai casi protonici.

Significatività e Rivendicazioni
Il documento sostiene che la composizione nucleare delle particelle accelerate è un parametro critico, sebbene spesso trascurato, nella modellazione delle sorgenti di neutrini ad alta energia. La significatività del lavoro risiede nel dimostrare che l'interpretazione dei dati multimessaggeri (neutrini vs. raggi gamma) è intrinsecamente legata alla compattezza e alla composizione iniettata della sorgente.

Nello specifico, gli autori argomentano che:

• La banda dei raggi gamma MeV è uno strumento diagnostico potente per le sorgenti di neutrini "nascoste" perché la cascata elettromagnetica rielabora i fotoni ad alta energia in questa banda.
• I dati MeV attuali (specificamente i limiti COMPTEL rianalizzati) forniscono già vincoli lievi ma significativi sulla composizione e sulle dimensioni della sorgente.
• In sorgenti compatte, i nuclei pesanti possono spiegare il flusso di neutrini con meno potenza totale e meno emissione di raggi gamma, potenzialmente risolvendo la tensione con le non rilevazioni di raggi gamma. Al contrario, in sorgenti meno compatte, i nuclei pesanti sono sfavoriti dai limiti dei raggi gamma.
• Le future strutture per i raggi gamma MeV (ad esempio, COSI, AMEGO-X, e-ASTROGAM) saranno essenziali per affinare queste sonde e distinguere tra scenari di iniezione protonica e nucleare nell'era emergente dell'astrofisica multimessaggera.

Gli autori rimangono modesti, notando che il loro modello è un "modello semplificato" (toy-model) del sistema disco-corona e che i lavori futuri dovranno incorporare le perdite per sincrotrone e geometrie più complesse (2D/3D) per catturare pienamente la fisica della sorgente.