Search for neutrino emission from blazar $γ$-ray flares accounting for possible neutrino time delays

Sintesi Tecnica: Ricerca dell'emissione di neutrini da flare di raggi $\gamma$ di blazar tenendo conto di possibili ritardi temporali dei neutrini

Definizione del Probleamento
Sebbene l'osservatorio per i neutrini IceCube abbia stabilito l'esistenza di neutrini astrofisici ad alta energia ($E_\nu \gtrsim 1$ TeV), le sorgenti specifiche che contribuiscono al flusso diffuso rimangono ampiamente non identificate, nonostante le prove che collegano eventi specifici ai blazar (ad esempio, TXS 0506+056) e alle galassie attive (ad esempio, NGC 1068). Le ricerche precedenti di correlazioni tra blazar e neutrini di IceCube tipicamente assumevano la coincidenza temporale tra l'arrivo del neutrino e i principali flare elettromagnetici (specificamente di raggi $\gamma$). Tuttavia, considerazioni teoriche suggeriscono che tale coincidenza potrebbe non sussistere. Se i neutrini sono prodotti tramite interazioni $p\gamma$, le scale temporali per l'accelerazione dei protoni ($t'_{acc}$) e per le perdite di energia dei protoni ($t'_{p\gamma}$) nel frame del getto possono essere di ordini di grandezza superiori rispetto a quelle degli elettroni responsabili dei flare elettromagnetici. Di conseguenza, i neutrini ad alta energia ($\gtrsim 100$ TeV) possono arrivare mesi o anni dopo il picco del flare di raggi $\gamma$ associato. Questo articolo affronta la mancanza di ricerche che tengano conto di tali potenziali ritardi temporali.

Metodologia
Gli autori hanno eseguito un'analisi di cross-matching tra gli alert di neutrini ad alta energia dal catalogo pubblico IceCat-1 e due distinti campioni di sorgenti di blazar, scansionando la ricerca di correlazioni su un ampio intervallo di ritardi temporali nel frame del getto ($t'_{delay}$).

1. Dataset:

   • Neutrini: 348 alert dal catalogo IceCat-1 ($E_\nu \gtrsim 100$ TeV), escludendo quelli con segnali IceTop coincidenti.
   • Liste di Sorgenti:
     • H21+: Un campione di 294 AGN dal database MOJAVE (Homan et al. 2021) con fattori Doppler ($D$), redshift ($z$) e curve di luce Fermi-LAT noti.
     • R24+: Un campione di 295 blazar dal catalogo CGRaBS (Rodrigues et al. 2024b) con fattori Doppler derivati dalla modellazione dell'energia spettrale (SED), incrociati anche con le curve di luce Fermi-LAT.
   • Curve di Luce: Curve di luce del flusso fotonico con binning mensile dal Fermi-LAT Light Curve Repository (LCR).

2. Framework Statistico:

   • Associazione Spaziale: I neutrini sono stati associati ai blazar se la sorgente ricadeva entro una regione di incertezza angolare ampliata (incorporando un'incertezza sistematica $\Delta = 0.78^\circ$).
   • Schema di Ponderazione: Per ogni alert, è stato calcolato un peso $w_\nu(t'_{delay})$ basato sulla significatività dell'alert, un fattore di penalità per le grandi incertezze angolari e un fattore temporale $G_{\nu a}$.
   • Fattore Temporale: Il fattore temporale ha collegato il tempo di arrivo del neutrino ($t_\nu$) al massimo globale della curva di luce $\gamma$ del blazar precedente l'arrivo ($t^{\gamma \max}_a$). Il fattore assegnava il peso massimo se la differenza temporale corrispondeva al ritardo previsto nel frame dell'osservatore ($t_{delay} = t'_{delay}(1+z)/D$) entro una dimensione di bin, applicando una penalità gaussiana per le deviazioni.
   • Statistica di Test (TS): Il TS globale è definito come la somma dei singoli pesi su tutti gli alert.
   • Test di Significatività: Il TS per i dati reali è stato confrontato con una distribuzione di valori di TS generati da 49.999 dataset mock dove le ascisse di ampiezza e i tempi di arrivo dei neutrini sono stati randomizzati. I $p$-values pre-trial sono stati calcolati e successivamente corretti per i trial (scansione su $t'_{delay}$).

Contributi Chiave

• Strategia di Ricerca Novella: Questo lavoro rappresenta il primo tentativo di cercare associazioni tra i neutrini di IceCube e i flare di raggi $\gamma$ di blazar tenendo specificamente conto di un ritardo temporale universale nel frame del getto ($t'_{delay}$) causato da una lenta accelerazione dei protoni o dalle perdite di energia.
• Doppio Campione di Sorgenti: Lo studio utilizza due cataloghi indipendenti di blazar con diversi metodi per determinare i fattori Doppler (argomenti di temperatura di brillamento rispetto alla modellazione SED) per testare la robustezza di potenziali correlazioni contro le incertezze nella cinematica dei getti.
• Scansione Sistematica del Ritardo: L'analisi scansiona $t'_{delay}$ da 1 mese a $\sim 3.5 \times 10^4$ giorni, coprendo l'intervallo teorico in cui la produzione di neutrini potrebbe ritardare rispetto all'emissione elettromagnetica.

Risultati

• Significatività Statistica: La ricerca non ha prodotto alcuna evidenza statisticamente significativa di una correlazione tra i neutrini IceCat-1 e gli AGN nelle liste di sorgenti H21+ o R24+.
  • Per la lista H21+, il minimo $p$-value pre-trial è stato 0.037 (a $t'_{delay} \approx 7.7 \times 10^3$ d).
  • Per la lista R24+, il minimo $p$-value pre-trial è stato 0.026 (a $t'_{delay} \approx 1.9 \times 10^3$ d).
  • Dopo la correzione per i trial, questi minimi corrispondono rispettivamente a $p$-values post-trial di 0.14 (1.06$\sigma$) e 0.11 (1.21$\sigma$), coerenti con l'ipotesi nulla.
• Top Candidati: Le associazioni più significative sono state guidate dagli alert IC120523B (collegato al flare del 2010 di 3C 454.3) e IC150926A (collegato a un flare di 3C 279). Sebbene questi eventi specifici mostrassero un allineamento temporale con i ritardi previsti, essi sono stati insufficienti per stabilire una correlazione globale attraverso l'intero dataset.
• Discrepanza nei Fattori Doppler: I tempi di ritardo ottimali differivano significativamente tra i due campioni di sorgenti (un fattore di 4), in gran parte a causa delle discrepanze nei fattori Doppler stimati per 3C 454.3 e 3C 279 tra le metodologie H21 e R24.

Significato e Rivendicazioni
Gli autori concludono con moderazione che i loro risultati non supportano una correlazione significativa tra i neutrini IceCat-1 e le specifiche popolazioni di blazar studiate, anche tenendo conto dei ritardi temporali. Suggeriscono diverse ragioni, non mutuamente esclusive, per questi risultati nulli:

1. Inefficienza: I blazar potrebbero essere sorgenti inefficienti di neutrini a $\sim 100$ TeV, contribuendo solo una piccola frazione ($\sim 0.5\%$) al flusso di IceCube.
2. Scala di Energia: Se i blazar sono sorgenti efficienti di neutrini, il picco della loro emissione potrebbe trovarsi a energie molto più elevate ($\gtrsim 10$ PeV), oltre la sensibilità primaria del campione di alert attuale di IceCube.
3. Assunzioni del Modello: L'assunzione di un ritardo temporale universale nel frame del getto ($t'_{delay}$) potrebbe essere invalida a causa di variazioni sorgente-a-sorgente nei meccanismi di accelerazione e nei campi magnetici.
4. Disaccoppiamento: L'accelerazione dei protoni e la dissipazione dell'energia potrebbero essere spazialmente o temporalmente disaccoppiate dai processi elettronici, il che significa che la produzione di neutrini non è causalmente legata ai flare di raggi $\gamma$ osservati.

Il paper sostiene che i miglioramenti futivi nella sensibilità richiederanno probabilmente dataset più ampi (ad esempio, i prossimi cataloghi IceCat-2 e CAZ) e l'inclusione di curve di luce radio, che sono note per ritardare rispetto ai flare di raggi $\gamma$ e potrebbero correlarsi meglio con le scale temporali di produzione dei neutrini.