Optical polarization variability and its relation to gamma-ray activity in blazars

Resumen Técnico: Variabilidad de la Polarización Óptica y su Relación con la Actividad de Rayos $\gamma$ en Blazares

Planteamiento del Problema
La polarización óptica sirve como una sonda crítica para comprender los procesos de aceleración de partículas y emisión de alta energía dentro de los chorros relativistas de agujeros negros supermasivos. Sin embargo, la naturaleza estocástica de la variabilidad de la polarización requiere programas de monitoreo de largo plazo y con un muestreo denso, los cuales consumen muchos recursos. Si bien investigaciones previas (por ejemplo, RoboPol, Kanata, Steward Observatory) han proporcionado datos valiosos, los programas individuales suelen sufrir de una cadencia o cobertura limitadas, lo que puede hacer que se pierdan variaciones rápidas o eventos de rotación específicos. Además, la conexión física entre el grado de polarización óptica y la actividad de rayos $\gamma$ sigue siendo objeto de debate, con evidencia conflictiva respecto a las correlaciones y los mecanismos de emisión subyacentes (por ejemplo, auto-Compton de sincrotrón frente a Compton externo).

Metodología
Los autores construyeron curvas de luz de polarización óptica con un muestreo denso para 15 blazares combinando datos públicamente disponibles de tres programas de monitoreo principales: RoboPol, Kanata y el Observatorio Steward. El estudio se centró en el grado y el ángulo de polarización en la banda R.

• Selección de la Muestra: De una lista inicial de 20 fuentes, se excluyeron cinco debido a una cobertura insuficiente, dejando una muestra final de 15 objetos (11 objetos de tipo BL Lacertae y 4 cuásares de radio de espectro plano [FSRQs]; 7 de pico de sincrotrón bajo [LSP], 4 de pico intermedio [ISP] y 4 de pico de sincrotrón alto [HSP]).
• Datos de Rayos $\gamma$: Se extrajeron las curvas de luz de rayos $\gamma$ correspondientes (0.1–100 GeV) con un binning de 3 días desde el repositorio del Telescopio Espacial Fermi.
• Análisis Estadístico:
  • Detección de Variabilidad: Se aplicaron Bloques Bayesianos (BB) a las curvas de luz del grado de polarización y de rayos $\gamma$ para identificar estados de variabilidad significativos. El grado de polarización fue normalizado por su mediana para comparar las fuentes de forma intrínseca.
  • Correlación: Se realizaron pruebas de correlación de rango de Spearman para investigar las relaciones entre el grado de polarización y el flujo de rayos $\gamma$ para fuentes individuales y submuestras (FSRQs, BLLs, LSP, ISP, HSP).
  • Detección de Rotación: Se empleó un algoritmo automatizado (Glykopoulou et al. 2026) utilizando BB, la prueba t de Student y pruebas binomiales para identificar rotaciones del ángulo de polarización (amplitudes $>90^\circ$).
  • Comparación de Distribución: Se utilizó la prueba de Anderson-Darling para comparar las distribuciones del grado de polarización durante los periodos de rotación frente a los periodos sin rotación.

Contribuciones Clave

1. Síntesis de Datos: La creación de las curvas de luz de polarización mejor muestreadas hasta la fecha para 15 blazares comunes mediante la fusión de datos de múltiples observatorios, disponibles vía Harvard Dataverse.
2. Catálogo de Rotación: La identificación de 64 rotaciones del ángulo de polarización en 12 fuentes, 39 de las cuales son reportadas por primera vez en la literatura.
3. Análisis de Correlación Sistemático: Un examen exhaustivo de la relación entre el grado de polarización y el flujo de rayos $\gamma$ a través de fuentes individuales y subclases de blazares, desafiando la suposición de una correlación universal.

Resultados

• Grado de Polarización vs. Flujo de Rayos $\gamma$: No se encontró una correlación fuerte o una anticorrelación entre el grado de polarización y el flujo de rayos $\gamma$ para las fuentes individuales o para ninguna submuestra (FSRQs, BLLs, LSP, ISP, HSP). Aunque dos fuentes (3C 66A y 3C 454.3) mostraron indicios de anticorrelación, la inspección visual sugirió que estos eran probablemente artefactos impulsados por valores atípicos o la mezcla de estados activos y no activos, más que relaciones físicas intrínsecas.
• Rotaciones del Ángulo de Polarización: El estudio identificó 64 rotaciones con una amplitud promedio de 13 grados y una duración promedio de 28 días. La amplitud máxima fue de 400 grados (AO 0235+164) y la duración máxima fue de 111 días (PKS 1510-089). Notablemente, no se detectaron rotaciones en las tres fuentes de pico de sincrotrón alto (HSP) (Mrk 501, 1ES 1959+650, 1ES 2344+514) dentro de la banda óptica, lo cual es consistente con hallazgos previos de que los blazares HSP rotan con menos frecuencia en longitudes de onda ópticas.
• Grado de Polarización Durante las Rotaciones: Para la muestra combinada, los autores confirmaron una tendencia donde los grados de polarización son menores durante los periodos de rotación en comparación con los periodos sin rotación (rechazando marginalmente la hipótesis nula al 4.5%). Sin embargo, para fuentes individuales, las distribuciones del grado de polarización durante los periodos de rotación y de no rotación fueron a menudo indistinguibles.

Significancia y Afirmaciones
El artículo concluye que la falta de una correlación universal entre el grado de polarización y el flujo de rayos $\gamma$, junto con las distribuciones de grado de polarización indistinguibles en fuentes individuales específicas, proporciona evidencia adicional para la hipótesis de múltiples mecanismos de emisión. Los autores sugieren que diferentes blazares, o incluso el mismo blazar en diferentes momentos, pueden operar bajo regímenes físicos distintos (por ejemplo, variando entre escenarios SSC y EC, o diferentes modelos de reconexión magnética).

Los resultados no favorecen un único mecanismo de emisión (como la reconexión impulsada por plásmidos o modelos de dos componentes) sobre otro. En su lugar, los datos sugieren que los impulsores físicos de la variabilidad de la polarización son complejos y probablemente varían a través de la población de blazares y a lo largo del tiempo. Los autores advierten que su tamaño de muestra es limitado y que se requieren muestras más grandes para confirmar estos hallazgos de manera definitiva. También señalan que trabajos futuros que involucren el modelado de flares dedicados y el análisis del espectro de energía (SED) podrían aislar mejor los procesos de emisión específicos.