Chasing the neutrino blazar candidates II: SED modeling with hadronic model

Sintesi Tecnica: Inseguimento dei candidati blazar neutrino II: Modellazione SED con modello adronico

Enunciato del Problema
Sebbene i blazar siano considerati candidati promettenti per sorgenti di neutrini ad alta energia, l'origine fisica della loro emissione di neutrini rimane incerta. Sebbene lo scenario leptonico (scattering Compton inverso) costituisca il quadro standard per la modellazione delle distribuzioni di energia spettrale (SED) dei blazar, il rilevamento di neutrini ad alta energia suggerisce un ruolo significativo degli adroni relativistici. Un lavoro precedente degli autori (Carta I) ha identificato 199 candidati blazar neutrino (NBC) utilizzando un approccio di apprendimento con pochi esempi. Questo studio mira ad estendere tale lavoro indagando le proprietà fisiche di questi candidati attraverso una modellazione SED a banda larga all'interno di un quadro adronico, specificamente per stimare la massima produzione possibile di neutrini e vincolare la composizione delle particelle dei getti dei blazar.

Metodologia
Gli autori hanno costruito SED a banda larga per 103 NBC (selezionati dai 199 originali sulla base di una copertura multi-lunghezza d'onda sufficiente) e per la sorgente confermata TXS 0506+056. La modellazione ha impiegato un codice adronico dipendente dal tempo a una zona (AM3) con le seguenti ipotesi e procedure chiave:

• Dominio Adronico: Per stimare il limite superiore dell'emissione di neutrini, il modello assume che i fotoni ad alta energia (HE) e a energia molto alta (VHE) siano prodotti principalmente tramite interazioni $p\gamma$ e cascate elettromagnetiche associate. Il contributo dallo scattering Compton inverso (IC) leptonico è fortemente soppresso.
• Configurazione Fisica: La regione di emissione è modellata come un "blob" sferico in movimento con un fattore di Lorentz di bulk $\Gamma$ in un campo magnetico $B$. Elettroni e protoni non termici sono iniettati con spettri di legge di potenza.
• Processi Modellati: Il modello tiene conto della radiazione di sincrotrone, della produzione di fotopioni ($p\gamma \to n\pi^+$ o $p\pi^0$), della produzione di coppie Bethe-Heitler ($p\gamma \to pe^+e^-$), della produzione di coppie fotone-fotone e del successivo decadimento di pioni e muoni per produrre neutrini. Viene applicata l'attenuazione della luce di fondo extragalattica (EBL) ai raggi gamma ad alta energia.
• Strategia di Adattamento: Gli autori hanno utilizzato l'algoritmo di ottimizzazione dello sciame di particelle (PSO) tramite il pacchetto PySwarms per adattare nove parametri liberi (ad esempio, intensità del campo magnetico, dimensione della regione, luminosità di iniezione delle particelle e indici spettrali) ai dati multi-lunghezza d'onda che spaziano dalle bande radio a quelle $\gamma$. I dati radio al di sotto di $10^{11}$ Hz sono stati esclusi a causa delle limitazioni dell'auto-assorbimento da sincrotrone nel modello a una zona. È stata applicata una raffinazione "ad occhio" per garantire rappresentazioni SED accettabili.
• Analisi di Correlazione: È stata eseguita un'analisi di correlazione parziale per investigare le relazioni tra la luminosità dei neutrini ($L_\nu$) e le luminosità elettromagnetiche (radio, ottica, X, $\gamma$), nonché le proprietà del motore centrale (massa del buco nero $M_{BH}$ e luminosità del disco di accrescimento $L_{disk}$), controllando al contempo la dipendenza dal redshift.

Contributi e Risultati Chiave

1. Modellazione SED e Vincoli sui Parametri: Lo studio ha vincolato con successo nove parametri chiave per 103 NBC. La modellazione rivela che l'intensità del campo magnetico tende ad essere elevata (il 66% delle sorgenti ha $\log B \ge 1.0$) per dissipare efficientemente l'energia degli elettroni tramite radiazione di sincrotrone, sopprimendo così l'emissione IC e permettendo ai protoni di dominare l'impulso ad alta energia. La potenza trasportata dai protoni relativistici ($P_p$) domina sia sulla potenza degli elettroni ($P_e$) che sulla potenza del campo magnetico ($P_B$), con un rapporto medio $\log(P_p/P_B) \simeq 1.67$.
2. Predizione nella Banda MeV: Una caratteristica distintiva del modello è la previsione di una prominente emissione di sincrotrone protonico che raggiunge il picco nella banda MeV. Per 99 dei 103 NBC, il picco di sincrotrone protonico cade entro 0,1–100 MeV. In 90 di queste sorgenti, il flusso di sincrotrone protonico nella banda MeV supera il flusso di produzione di coppie Bethe-Heitler nella banda GeV.
3. Correlazioni con l'Emissione Elettromagnetica:
   • Raggi Gamma e Raggi X: Sono state trovate forti correlazioni positive tra la luminosità dei neutrini e le luminosità nei raggi $\gamma$ ($r=0,78$) e nei raggi X ($r=0,52$). Ciò supporta un'origine adronica in cui protoni ad alta energia producono sia neutrini che fotoni ad alta energia.
   • Ottico (Banda R): Esiste una correlazione debole/moderata ($r=0,30$) tra la luminosità dei neutrini e quella ottica nella banda R. Gli autori suggeriscono che ciò possa derivare dal fatto che i fotoni di sincrotrone ottici dal getto fungono da fotoni seed per la produzione di fotopioni, o potenzialmente a causa di una connessione con il disco di accrescimento.
   • Radio: Non è stata trovata alcuna correlazione significativa tra la luminosità dei neutrini e quella radio a 1,4 GHz.
   • Motore Centrale: Non è stata trovata alcuna correlazione statisticamente significativa tra la luminosità dei neutrini e la massa del buco nero ($M_{BH}$) o la luminosità del disco di accrescimento ($L_{disk}$), sebbene fosse visivamente evidente una debole tendenza positiva con $L_{disk}$.
4. Prospettive di Rilevamento: Basandosi sui flussi massimi di neutrini previsti dal modello:
   • IceCube: 3 NBC sono potenzialmente rilevabili.
   • Futuri Strumenti: Fino a 22 sorgenti potrebbero essere rilevabili da KM3NeT, 45 da NEON e 62 da TRIDENT.
   • Il flusso di neutrini sommato del campione non supera i limiti attuali di stacking di IceCube, suggerendo che lo scenario "massimo neutrino" sia fisicamente realizzabile.

Significato e Affermazioni
Il paper afferma che i suoi risultati forniscono previsioni verificabili per future osservazioni multi-messaggero e offrono nuove intuizioni sulla composizione dei getti dei blazar. Nello specifico:

• La Finestra MeV: Gli autori evidenziano la banda MeV come una finestra osservativa critica per distinguere tra scenari leptonici e adronici, poiché il modello adronico prevede un distinto picco di sincrotrone protonico in questo intervallo, assente nei modelli leptonici standard.
• Composizione del Getto: La modellazione supporta uno scenario in cui i protoni relativistici dominano il bilancio energetico dei getti in questi candidati, coerentemente con le prove multi-messaggero provenienti da TXS 0506+056.
• Limiti Superiori: I risultati rappresentano un limite superiore all'emissione di neutrini; i flussi effettivi potrebbero essere inferiori, il che spiega perché molti candidati non sono ancora stati confermati come emettitori di neutrini nonostante si trovino al di sopra delle curve di sensibilità attuali nel modello.
• Osservazioni Future: Lo studio sottolinea che missioni di raggi $\gamma$ MeV di prossima generazione (ad esempio, e-ASTROGAM, AMEGO) e rivelatori di neutrini di prossima generazione (KM3NeT, NEON, TRIDENT) sono essenziali per verificare questi scenari adronici e vincolare la composizione delle particelle dei getti dei blazar.

Gli autori mantengono un tono modesto riguardo all'origine della correlazione ottico-neutrino, notando che potrebbe derivare da effetti di selezione del campione o da molteplici meccanismi (getto vs disco), e riconoscono che la necessità di potenze di getto super-Eddington in alcune sorgenti potrebbe richiedere modelli multi-zona per una piena riconciliazione.