Connecting the long-term variability behaviour of active galactic nuclei to their central engines

Resumen Técnico: Vinculando el comportamiento de variabilidad a largo plazo de los núcleos galácticos activos con sus motores centrales

Planteamiento del Problema
Los Núcleos Galácticos Activos (NGA) exhiben variabilidad multibanda impulsada por chorros relativistas lanzados por agujeros negros supermasivos. Si bien el análisis de la Densidad Espectral de Potencia (DEP) es un método estándar para identificar escalas temporales de variabilidad características (específicamente la inflexión de baja frecuencia en el espectro de ruido rojo), la interpretación física de estas escalas temporales sigue siendo debatida. Estudios previos sobre binarias de rayos X y NGA sugieren que las inflexiones de la DEP escalan con la masa del agujero negro ($M_{BH}$) y la tasa de acreción ($\dot{M}$), pero la característica específica de la curva de luz que corresponde a la inflexión de la DEP en radio no está clara. Es incierto si la inversa de la frecuencia de inflexión corresponde a la duración de destellos individuales o a la separación temporal entre ellos. Además, la relación entre estas escalas temporales de variabilidad y el tamaño físico de la región de emisión (a menudo asociada con la región de transición del chorro) y los parámetros del motor central requiere una mayor investigación utilizando una muestra más grande que la disponible previamente.

Metodología
Los autores analizaron curvas de luz de radio a largo plazo de 54 NGA monitoreados en el Observatorio de Radio Metsähovi de la Universidad Aalto (MRO) a 37 GHz. La muestra incluyó 11 fuentes con escalas temporales de DEP previamente restringidas (de Kankkunen et al. 2025, "Papel I") y 43 fuentes adicionales (incluyendo 3 casos limítrofes) seleccionadas por su alta densidad de datos y destellos prominentes.

La metodología central involucró:

1. Descomposición de Curvas de Luz: Se empleó el código Magnetron para descomponer las curvas de luz en destellos individuales utilizando modelado jerárquico bayesiano y muestreo anidado difusivo. Se utilizó un proceso de Ornstein-Uhlenbeck (OU) para modelar la variabilidad de fondo subyacente, evitando la identificación errónea de variaciones turbulentas menores como destellos.
2. Extracción de Escalas Temporales: Para cada fuente, los autores calcularon la duración media y la separación media de los picos de los destellos más brillantes (definidos como aquellos que superan la amplitud media en un dibujo posterior). También se extrajeron los tiempos de ascenso para los destellos más brillantes (amplitud > 50% del máximo) para servir como proxies del tamaño de la región de emisión.
3. Estimación de Parámetros:
   • Factores Doppler ($\delta_{var}$): Estimados utilizando temperaturas de brillo de variabilidad derivadas de los tiempos de ascenso y amplitudes de los destellos, comparadas contra una distribución gaussiana de temperaturas de brillo intrínsecas.
   • Factores de Lorentz ($\Gamma$): Calculados utilizando $\delta_{var}$ y velocidades aparentes del chorro ($\beta_{app}$) de la encuesta MOJAVE.
   • Parámetros del Motor Central: Las masas de los agujeros negros ($M_{BH}$) y las luminosidades del disco de acreción ($L_{acc}$) se estimaron utilizando luminosidades de líneas anchas (H$\beta$, MgII o [OIII]) y valores de ancho completo a la mitad del máximo (FWHM). Las tasas de acreción de masa normalizadas ($\dot{m}$) se derivaron asumiendo una eficiencia radiativa de $\eta = 0.4$.
4. Análisis Estadístico: Se realizaron pruebas de correlación de rango de Kendall para evaluar las asociaciones entre las escalas temporales derivadas (duraciones de destellos, separaciones y tiempos de ascenso) y los parámetros físicos ($M_{BH}$, $\dot{m}$, $L_{acc}$, $\Gamma$). Se llevaron a cabo análisis de correlación parcial para tener en cuenta posibles dependencias del corrimiento al rojo ($z$).

Resultados Clave

• Escala Temporal de la DEP vs. Características de la Curva de Luz: La inversa de la frecuencia de inflexión de la DEP corresponde mejor a la duración media de los destellos más brillantes. Si bien las separaciones medias de los destellos también se correlacionan con la escala temporal de la DEP, las duraciones de los destellos muestran un acuerdo uno a uno más fuerte, particularmente para las 11 fuentes restringidas.
• Conexión con el Motor Central: Utilizando las duraciones y separaciones de los destellos como proxies para la escala temporal de la DEP, los autores encontraron una correlación positiva moderada con la relación entre la masa del agujero negro y la tasa de acreción de masa normalizada ($M_{BH}/\dot{m}$). Esta correlación es significativa en el marco del observador, pero se debilita o se vuelve no significativa cuando se corrige al marco de emisión. La correlación parece estar impulsada principalmente por $M_{BH}$, aunque el papel de $\dot{m}$ no puede excluirse definitivamente debido a las limitaciones del tamaño de la muestra.
• Tamaño de la Región de Emisión: Los tiempos de ascenso de los destellos más brillantes, utilizados como proxies para el ancho de la región de transición del chorro, mostraron correlaciones tentativas débiles con $M_{BH}$, $L_{acc}$ y $\Gamma$. Sin embargo, se encontró que estas asociaciones estaban impulsadas por una dependencia común del corrimiento al rojo. Cuando se tuvo en cuenta el corrimiento al rojo mediante correlación parcial, la significancia de las relaciones entre los tiempos de ascenso y los parámetros del motor central desapareció.
• Dependencia del Corrimiento al Rojo: Se encontró una correlación positiva significativa entre el corrimiento al rojo y las duraciones de los destellos en el marco de emisión, mientras que existía una correlación negativa entre el corrimiento al rojo y las separaciones de los destellos en el marco de emisión (lo que implica una tasa de eyección más alta para fuentes de mayor corrimiento al rojo). Estas tendencias opuestas pueden explicar la falta de correlación entre las duraciones/separaciones de destellos observadas y los factores Doppler.

Significado y Afirmaciones
El artículo afirma que las escalas temporales características derivadas de las inflexiones de la DEP en radio están físicamente asociadas con el motor central de los NGA, escalando específicamente con $M_{BH}/\dot{m}$. Esto respalda la hipótesis de que el motor central modula la variabilidad del chorro, similar a los hallazgos en el dominio de rayos X. El estudio sugiere que la duración media de los destellos más brillantes es el proxy de curva de luz más apropiado para la escala temporal de la inflexión de la DEP.

Los autores notan modestamente que, si bien existen correlaciones con los parámetros del motor central, a menudo son débiles y potencialmente confundidas por efectos de corrimiento al rojo y sesgos de muestra (por ejemplo, la selección de fuentes con destellos discernibles). Los resultados relacionados con el tamaño de la región de emisión y el factor de Lorentz se describen como "correlaciones tentativas débiles" que podrían no sostenerse una vez que se eliminan las dependencias del corrimiento al rojo. En consecuencia, el artículo concluye que, aunque se sugiere una conexión con el motor central, se requieren muestras más grandes y estimaciones de parámetros más precisas para caracterizar definitivamente estas relaciones. El trabajo no propone nuevos diseños experimentales, sino que valida el uso de la descomposición de curvas de luz como un método para extender los análisis basados en DEP a una muestra más amplia de NGA.