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

Resumen Técnico: Variabilidad Multilongitud de Onda del Candidato a Neutrino de Alta Energía PKS 0735+178 durante Tres Décadas

Planteamiento del Problema
PKS 0735+178 es un objeto BL Lac clásico asociado con el evento de neutrino de alta energía IC211208A (171 TeV). Aunque se sabe que la fuente es conocida por su extrema variabilidad multi-longitud de onda y un chorro relativista visto desde ángulos pequeños, la conexión física entre las regiones de emisión en diferentes frecuencias, la naturaleza de su comportamiento cuasi-periódico a largo plazo y los mecanismos específicos que impulsan la actividad en torno al evento de neutrino siguen siendo insuficientemente comprendidos. Estudios previos se han centrado principalmente en propiedades espectrales de corto plazo o en llamaradas aisladas, careciendo de un análisis unificado de la evolución temporal de la fuente sobre una base de varias décadas para restringir la geometría del chorro y la propagación de perturbaciones.

Metodología
Los autores presentan un análisis exhaustivo de aproximadamente 30 años (1995–2026) de datos multi-longitud de onda, integrando:

• Datos de Radio: Observaciones cuasi-simultáneas a 1–22 GHz del radiotelescopio RATAN-600 y datos de 230 GHz del Submillimeter Array (SMA).
• Datos Ópticos: Fotometría en banda R de los telescopios Zeiss-1000, AS-500/2, KAIT, ZTF y CRTS.
• Datos de Rayos Gamma: Curvas de luz de Fermi-LAT (1.5–300 GeV) construidas mediante un agrupamiento adaptativo para resolver estructuras finas durante los estados de llamarada.
• Restricciones de VLBI: Mediciones de desplazamiento del núcleo (core-shift) a 15 y 43 GHz utilizando datos de MOJAVE y VLBA-BU-BLAZAR para derivar la cinemática del chorro y los parámetros físicos.

El análisis empleó diversas técnicas de modelado estadístico y físico:

1. Correlación Cruzada: Se utilizaron Funciones de Correlación Discreta (DCF) para medir los retardos temporales entre bandas.
2. Análisis Espectral: Los espectros de radio se agruparon en tres episodios principales de llamarada para rastrear la evolución del índice espectral ($\alpha$).
3. Búsqueda de Cuasi-periodicidad: Se aplicaron los periodogramas de Transformada de Wavelet Ponderada (WWZ) y de Lomb-Scargle Generalizada (GLS) para detectar escalas de tiempo características, evaluando la significancia contra fondos de ruido rojo mediante simulaciones de Monte Carlo.
4. Modelado Físico:
   • Propagación de Choques: Los retardos temporales se ajustaron a una dependencia de ley de potencia ($\Delta t \propto \nu^{-1/k_r}$) para probar un modelo de chorro estratificado por opacidad.
   • Precesión del Chorro: Se utilizó un modelo de optimización Bayesiana (UltraNest) para ajustar las curvas de luz a un escenario de precesión, probando si los efectos geométricos podrían explicar la modulación a largo plazo.
   • Desplazamiento del Núcleo de VLBI: Utilizado para estimar el ángulo de visión del chorro, la velocidad aparente, la fuerza del campo magnético y el flujo magnético en relación con el límite del Disco Magnetizados por Arresto Magnético (MAD).

Resultos Clave

• Retardos Temporales Dependientes de la Frecuencia: Las curvas de luz en las bandas de radio, óptica y rayos gamma están correlacionadas, pero la emisión de radio presenta un retraso respecto a las frecuencias más altas. Los retardos aumentan sistemáticamente hacia las frecuencias más bajas, alcanzando $\sim$1200 días a 1–2 GHz respecto a los rayos gamma. Esto sigue una ley de potencia con $k_r = 1.96 \pm 0.24$, lo cual es consistente con un chorro estratificado por opacidad donde las perturbaciones se propagan río abajo.
• Llamaradas del Evento de Neutrino: Se detectaron llamaradas cortas y brillantes en las bandas óptica y de rayos gamma dentro de los días previos al evento IC211208A. Se detectó un retraso robusto de 12 días entre los picos óptico y de rayos gamma (el óptico liderando), lo que sugiere regiones de emisión no co-espaciales o contribuciones de procesos como la dispersión Compton externa. La tercera llamarada mayor de radio (2019–2026) exhibe espectros invertidos o planos en bajas frecuencias, indicando una región de emisión compacta y ópticamente gruesa con aceleración de partículas mejorada.
• Variabilidad Cuasi-periódica: Se detectó de forma robusta ($\ge 3\sigma$) una escala de tiempo de variabilidad característica de $\sim$10–11 años en las bandas de radio y milimétricas (5–230 GHz). Esta modulación es coherente en estas bandas, pero no es estrictamente periódica. En contraste, los datos ópticos y de rayos gamma muestran señales de periodos más cortos ($\sim$5–9 años) que están localizadas en el tiempo y son menos coherentes.
• Geometría y Física del Chorro:
  • Las mediciones del desplazamiento del núcleo por VLBI arrojan un ángulo de visión del chorro de $\theta \approx 3.6^\circ$ y una velocidad aparente de $\beta_{app} \approx 0.8c$, consistentes con las velocidades de propagación derivadas de los retardos temporales.
  • El campo magnético a 1 pc se estima en $\sim$0.36 G. La relación del flujo magnético con el límite MAD sugiere que el disco de acreción no se encuentra en un estado MAD a los niveles de luminosidad actuales.
• Limitaciones del Modelo de Precesión: Aunque un modelo de precesión del chorro puede reproducir la forma general de las curvas de luz, el ángulo de visión derivado ($\sim$27.4°) contradice las restricciones de VLBI ($\theta \sim 1^\circ-3^\circ$). Esto indica que la precesión por sí sola no puede explicar la variabilidad observada.

Significancia y Reivindicaciones
El artículo sostiene que la variabilidad de PKS 0735+178 está gobernada por dos procesos físicos distintos que operan en diferentes escalas de tiempo:

1. Modulación a Largo Plazo: Una variación lenta en la liberación de energía en la base del chorro impulsa una modulación cuasi-periódica ($\sim$11 años), probablemente vinculada al flujo de acreción o a la dinámica de la base del chorro, la cual regula la tasa de inyección de perturbaciones.
2. Llamaradas de Corto Plazo: Las llamaradas individuales, incluyendo la asociada con el evento de neutrino, surgen de choques que se propagan río abajo a través de un chorro estratificado por opacidad; los retardos dependientes de la frecuencia y la evolución espectral se explican naturalmente mediante el cruce de los choques a través de las fotosferas de sincrotrón a diferentes distancias.

El estudio proporciona la primera medición fiable del retraso de 12 días entre las emisiones óptica y de rayos gamma en esta fuente, desafiando los escenarios simples de auto-Compton de sincrotrón de una sola zona. Además, la consistencia entre las restricciones independientes de los desplazamientos de núcleo de VLBI y los retardos temporales de radio respalda fuertemente un escenario unificado de propagación de choques en un flujo estratificado. El agrupamiento temporal del evento de neutrino con la llamarada multi-longitud de onda refuerza la conexión entre PKS 0735+178 y la producción de neutrinos de alta energía, ocurriendo probablemente en las regiones internas y ópticamente gruesas del chorro.