Multiwavelength variability of the high-energy neutrino candidate PKS 0735+178 over three decades

Sintesi Tecnica: Variabilità Multiwavelength del Candidato Neutrino ad Alta Energia PKS 0735+178 su Tre Decenni

Problema
PKS 0735+178 è un oggetto BL Lac classico associato all'evento di neutrino ad alta energia IC211208A (171 TeV). Sebbene la sorgente sia nota per la sua estrema variabilità multiwavelength e per un getto relativistico osservato con piccoli angoli, la connessione fisica tra le regioni di emissione a diverse frequenze, la natura del suo comportamento quasi-periodico a lungo termine e i meccanismi specifici che guidano l'attività intorno all'evento del neutrino rimangono scarsamente compresi. Studi precedenti si sono concentrati principalmente sulle proprietà spettrali a breve termine o su flare isolati, mancando di un'analisi unificata dell'evoluzione temporale della sorgente su una linea di base multi-decennale per vincolare la geometria del getto e la propagazione delle perturbazioni.

Metodologia
Gli autori presentano un'analisi completa di circa 30 anni (1995–2026) di dati multiwavelength, integrando:

• Dati Radio: Osservazioni quasi-simultanee a 1–22 GHz dal radiotelescopio RATAN-600 e dati a 230 GHz dal Submillimeter Array (SMA).
• Dati Ottici: Fotometria in banda R dai telescopi Zeiss-1000, AS-500/2, KAIT, ZTF e CRTS.
• Dati Gamma: Curve di luce Fermi-LAT (1.5–300 GeV) costruite mediante binning adattivo per risolvere le strutture fini durante gli stati di flare.
• Vincoli VLBI: Misurazioni del core-shift a 1 al 15 e 43 GHz utilizzando i dati MOJAVE e VLBA-BU-BLAZAR per derivare la cinematica del getto e i parametri fisici.

L'analisi ha impiegato diverse tecniche statistiche e di modellazione fisica:

1. Correlazione Incrociata: Le Funzioni di Correlazione Discreta (DCF) sono state utilizzate per misurare i ritardi temporali tra le bande.
2. Analisi Spettrale: Gli spettri radio sono stati raggruppati in tre grandi episodi di flare per tracciare l'evoluzione dell'indice spettrale ($\alpha$).
3. Ricerca della Quasi-Periodicità: Sono stati applicati i periodogrammi Weighted Wavelet Z-transform (WWZ) e Generalized Lomb-Scargle (GLS) per rilevare le scale temporali caratteristiche, valutando la significatività rispetto ai background di rumore rosso tramite simulazioni Monte Carlo.
4. Modellazione Fisica:
   • Propagazione d'Urto: I ritardi temporali sono stati adattati a una dipendenza per legge di potenza ($\Delta t \propto \nu^{-1/k_r}$) per testare un modello di getto stratificato per opacità.
   • Precessione del Getto: Un modello di ottimizzazione bayesiana (UltraNest) è stato usato per adattare le curve di luce a uno scenario di precessione, testando se gli effetti geometrici potessero spiegare la modulazione a lungo termine.
   • VLBI Core Shift: Utilizzato per stimare l'angolo di visione del getto, la velocità apparente, l'intensità del campo magnetico e il flusso magnetico rispetto al limite del Disco Magnetizzato (MAD).

Risultati Chiave

• Ritardi Temporali Dipendenti dalla Frequenza: Le curve di luce attraverso le bande radio, ottica e gamma sono correlate, ma l'emissione radio ritarda rispetto alle frequenze più alte. I ritardi aumentano sistematicamente verso le frequenze più basse, raggiungendo circa 1200 giorni a 1–2 GHz rispetto ai raggi gamma; ciò segue una legge di potenza con $k_r = 1.96 \pm 0.24$, coerente con un getto stratificato per opacità in cui le perturbazioni si propagano a valle.
• Flare dell'Evento del Neutrino: Flare brevi e brillanti nelle bande ottica e gamma si sono verificati entro pochi giorni dall'evento IC211208A. È stato rilevato un ritardo robusto di 12 giorni tra i picchi ottici e gamma (l'ottico precede), suggerendo regioni di emissione non co-spaziali o contributi da processi come lo scattering Compton esterno. Il terzo grande flare radio (2019–2026) mostra spettri invertiti o piatti a basse frequenze, indicando una regione di emissione compatta e otticamente spessa con accelerazione di particelle potenziata.
• Variabilità Quasi-Periodica: Una scala temporale di variabilità caratteristica di $\sim$10–11 anni è robustamente rilevata ($\ge 3\sigma$) nelle bande radio e millimetriche (5–230 GHz). Questa modulazione è coerente attraverso queste bande, ma non è strettamente periodica. Al contrario, i dati ottici e gamma mostrano segnali più deboli e a periodo più breve ($\sim$5–9 anni), che sono localizzati nel tempo e meno coerenti.
• Geometria e Fisica del Getto:
  • Le misurazioni del core-shift VLBI forniscono un angolo di visione del getto $\theta \approx 3.6^\circ$ e una velocità apparente $\beta_{app} \approx 0.8c$, coerenti con le velocità di propagazione derivate dai ritardi temporali.
  • Il campo magnetico a 1 pc è stimato in $\sim$0.36 G. Il rapporto tra il flusso magnetico e il limite MAD suggerisce che il disco di accrescimento non si trovi in uno stato MAD ai livelli di luminosità attuali.
• Limitazioni del Modello di Precessione: Sebbene un modello di precessione del getto possa riprodurre la forma generale delle curve di luce, i parametri derivati (angolo di visione $\sim$27.4°) contraddicono i vincoli VLBI ($\theta \sim 1^\circ-3^\circ$). Ciò indica che la precessione da sola non può spiegare la variabilità osservata.

Significatività e Rivendicazioni
L'articolo sostiene che la variabilità di PKS 0735+178 è governata da due distinti processi fisici che operano su diverse scale temporali:

1. Modulazione a Lungo Termine: Una variazione lenta nel rilascio di energia alla base del getto guida una modulazione quasi-periodica ($\sim$11 anni), probabilmente legata al flusso di accrescimento o alla dinamica della base del getto, che regola il tasso di iniezione delle perturbazioni.
2. Flare a Breve Termine: I singoli flare, incluso quello associato all'evento del neutrino, derivano da urti che si propagano a valle attraverso un getto stratificato per opacità; i ritardi dipendenti dalla frequenza e l'evoluzione spettrale sono naturalmente spiegati dagli urti che attraversano le fotosfere di sincrotrone a diverse distanze.

Lo studio fornisce la prima misurazione affidabile del ritardo di 12 giorni tra le emissioni ottiche e gamma in questa sorgente, sfidando i semplici scenari di sincrotrone Compton auto-Compton a zona singola. Inoltre, la coerenza tra i vincoli indipendenti dai core-shift VLBI e i ritardi temporali radio supporta fortemente uno scenario unificato di propagazione di urti in un flusso stratificato. Il raggruppamento temporale dell'evento del neutrino con il flare multiwavelength rafforza la connessione tra PKS 0735+178 e la produzione di neutrini ad alta energia, che probabilmente avviene nelle regioni interne e otticamente spesse del getto.