Automating the detection of polarization angle rotations in blazars. Re-analysis of RoboPol data reveals 27 new rotations

Resumen Técnico: Automatización de la Detección de Rotaciones del Ángulo de Polarización en Blázares

Enunciado del Problema
Los blázares exhiben una variabilidad significativa en su polarización óptica, específicamente en el Ángulo de Posición del Vector Eléctrico (EVPA). Las rotaciones grandes y sistemáticas del EVPA sirven como sondas únicas de las estructuras del campo magnético y la aceleración de partículas dentro de los chorros relativistas. Si bien estudios anteriores, particularmente aquellos dentro del programa RoboPol, establecieron catálogos estadísticos de estas rotaciones, los métodos de identificación dependían de procedimientos sensibles a la incertidumbre de la medición y a la variabilidad intrínseca. Además, la segmentación manual o semiautomática podía introducir sesgos subjetivos y potencialmente pasar por alto eventos que no se ajustaban estrictamente a las restricciones de suavidad predefinidas. Existía la necesidad de una tubería completamente automatizada y reproducible para identificar robustamente las rotaciones del EVPA en grandes muestras y extender el análisis a datos previamente no examinados, como la temporada de observación final de RoboPol de 2016–2017.

Metodología
Los autores desarrollaron una tubería completamente automatizada diseñada para procesar datos de monitoreo de RoboPol. El flujo de trabajo consta de cuatro etapas principales:

1. Preprocesamiento y Selección de Datos: La tubería ingiere mediciones de polarización desde archivos CSV, convirtiendo las marcas de tiempo a Fechas Julianas Modificadas (MJD). Para garantizar la fiabilidad estadística, solo se retienen las mediciones con una significancia del grado de polarización de al menos tres ($PD/errPD \geq 3$).
2. Resolución de la Ambigüedad de 180°: Dado que el EVPA se define módulo 180°, pueden ocurrir discontinuidades artificiales. La tubería implementa un procedimiento de ajuste ponderado por el error. Para cada nuevo punto de datos, el ángulo medido se compara con el valor ajustado anterior. Si la diferencia supera los 90° (después de ponderar por las incertidumbres), se suman múltiplos enteros de 180° para minimizar la diferencia ponderada. Los huecos mayores de 30 días se tratan como límites de segmento naturales donde la desambiguación es poco fiable.
3. Segmentación con Bloques Bayesianos: Las series temporales de EVPA ajustadas se segmentan utilizando el algoritmo de Bloques Bayesianos (Scargle et al. 2013) con una probabilidad de falsa alarma ($p_0$) de 0.001. Este agrupamiento adaptativo crea bloques estadísticamente homogéneos sin imponer restricciones de suavidad sobre la curva, permitiendo la detección de eventos que podrían dividirse o pasarse por alto mediante métodos basados en derivadas.
4. Detección y Validación de Rotaciones:
   • Segmentación: Las series temporales se dividen en segmentos contiguos que contienen al menos cuatro puntos de datos sin huecos que superen los 30 días.
   • Identificación de Extremos: Se extraen máximos y mínimos locales de las medias de los Bloques Bayesianos.
   • Selección de Candidatos: Los pares de extremos adyacentes (máx→mín o mín→máx) se tratan como rotaciones candidatas si la amplitud del EVPA supera los 90°.
   • Criterios Estructurales y Estadísticos: Un intervalo candidato se retiene solo si comprende al menos tres Bloques Bayesianos y contiene al menos cuatro puntos de datos. Cada intervalo se somete a una prueba t de una muestra (desviación de una media de referencia) y a una prueba binomial (asimetría en la distribución del EVPA).
   • Salida: La tubería calcula la amplitud de la rotación, la duración ($\Delta$MJD) y la velocidad angular, propagando las incertidumbres desde los errores de medición.

Contribuciones Clave

• Generación Automatizada de Catálogos: La tubería identificó 48 rotaciones del EVPA a través de 25 fuentes únicas de RoboPol. Varias fuentes exhibieron múltiples eventos (por ejemplo, RBPLJ2232+1143 con 6 eventos).
• Descubrimiento de Nuevos Eventos: El análisis reveló 27 rotaciones previamente no reportadas. Esto incluye 16 eventos de las temporadas 2013–2016 que fueron omitidos por catálogos anteriores y 11 eventos nuevos de la temporada 2016–2017, que no habían sido analizados previamente para rotaciones.
• Comparación Metodológica: El estudio proporciona una comparación sistemática entre el nuevo enfoque de Bloques Bayesianos y los métodos de los catálogos anteriores de RoboPol.
• Análisis de Correlación con Rayos Gamma: Los autores realizaron un análisis de correlación entre las rotaciones detectadas y las curvas de luz contemporáneas de rayos $\gamma$ de Fermi–LAT utilizando coeficientes de correlación de rango de Spearman, validados mediante un análisis de Monte Carlo (50,000 realizaciones).

Resultados

• Estadísticas de Rotación: Las rotaciones detectadas abarcan amplitudes desde 90.8° hasta 359.7° (media 159.7°), duraciones desde 7.0 hasta 111.3 días (media 44.4 días) y tasas de rotación que promedian 5.0°/día.
• Diferencias Sistemáticas: En comparación con el catálogo anterior de RoboPol, el método de Bloques Bayesianos produce rotaciones que son, en promedio, ~10% mayores en amplitud, ~2 veces más largas en duración y ~2/3 más lentas en velocidad angular. Los autores atribuyen estas desviaciones a la naturaleza de agrupamiento adaptativo de los Bloques Bayesianos, que tiende a fusionar oscilaciones adyacentes en eventos más largos, en contraste con la segmentación manual de los catálogos anteriores.
• Correlación con Rayos Gamma:
  • Duración vs. Actividad: Se encontró una correlación positiva significativa entre la duración de la rotación y el flujo de energía de rayos $\gamma$ pico-promedio ($\rho = 0.601$, $p_{med} = 6.4 \times 10^{-6}$). Las rotaciones más largas coinciden consistentemente con una actividad de rayos $\gamma$ aumentada.
  • Amplitud vs. Actividad: No se encontró una correlación significativa entre la amplitud de la rotación y el brillo de rayos $\gamma$ ($\rho = 0.146$, $p_{med} = 0.322$). La amplitud por sí sola no es predictiva de la mejora en rayos $\gamma$.

Significado y Afirmaciones
El artículo afirma que la tubería desarrollada minimiza los sesgos subjetivos y proporciona un marco reproducible para el catalogado de la variabilidad de polarización. Al recuperar eventos pasados por alto y extender la cobertura temporal de la muestra, el estudio enriquece el inventario estadístico de las rotaciones del EVPA.

Físicamente, los resultados sugieren que la emisión de rayos $\gamma$ de alta energía está vinculada a reconfiguraciones sostenidas del campo magnético del chorro (indicado por la correlación con la duración), en lugar de la magnitud de la rotación por sí sola. La falta de correlación con la amplitud implica un componente estocástico en los mecanismos de emisión. Los autores concluyen que el monitoreo de la polarización sirve como una sonda crítica para la dinámica del chorro y la aceleración de partículas, y que la tubería automatizada ofrece una base rigurosa para futuros estudios multibanda, incluidas las aplicaciones potenciales a datos de polarización en rayos X.