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

Resumo Técnico: Conectando o comportamento de variabilidade de longo prazo de núcleos galácticos ativos aos seus motores centrais

Declaração do Problema
Núcleos Galácticos Ativos (NGAs) exibem variabilidade multi-comprimento de onda impulsionada por jatos relativísticos lançados por buracos negros supermassivos. Embora a análise de Densidade Espectral de Potência (PSD) seja um método padrão para identificar escalas de tempo características de variabilidade (especificamente a quebra de baixa frequência no espectro de ruído vermelho), a interpretação física dessas escalas de tempo permanece debatida. Estudos anteriores em binárias de raios-X e NGAs sugerem que as dobras da PSD escalonam com a massa do buraco negro ($M_{BH}$) e a taxa de acreção ($\dot{M}$), mas a característica específica da curva de luz correspondente à dobra da PSD em rádio não está clara. É incerto se o inverso da frequência de dobra corresponde à duração de erupções individuais ou à separação temporal entre elas. Além disso, a relação entre essas escalas de tempo de variabilidade e o tamanho físico da região de emissão (frequentemente associada à região de transição do jato) e aos parâmetros do motor central requer investigação adicional utilizando uma amostra maior do que as anteriormente disponíveis.

Metodologia
Os autores analisaram curvas de luz de rádio de longo prazo de 54 NGAs monitorados no Observatório de Rádio Metsähovi da Universidade Aalto (MRO) a 37 GHz. A amostra incluiu 11 fontes com escalas de tempo de PSD previamente restringidas (de Kankkunen et al. 2025, "Papel I") e 43 fontes adicionais (incluindo 3 casos limítrofes) selecionadas por alta densidade de dados e erupções proeminentes.

A metodologia central envolveu:

1. Decomposição da Curva de Luz: O código Magnetron foi empregado para decompor curvas de luz em erupções individuais usando modelagem hierárquica bayesiana e amostragem aninhada difusiva. Um processo de Ornstein-Uhlenbeck (OU) foi utilizado para modelar a variabilidade de fundo subjacente, prevenindo a identificação equivocada de variações turbulentas menores como erupções.
2. Extração de Escalas de Tempo: Para cada fonte, os autores calcularam a duração média e a separação média de pico das erupções mais brilhantes (definidas como aquelas que excedem a amplitude média em uma extração posterior). Tempos de subida também foram extraídos para as erupções mais brilhantes (amplitude > 50% do máximo) para servir como proxies para o tamanho da região de emissão.
3. Estimativa de Parâmetros:
   • Fatores Doppler ($\delta_{var}$): Estimados usando temperaturas de brilho de variabilidade derivadas dos tempos de subida e amplitudes das erupções, comparados contra uma distribuição gaussiana de temperaturas de brilho intrínsecas.
   • Fatores de Lorentz ($\Gamma$): Calculados usando $\delta_{var}$ e velocidades aparentes do jato ($\beta_{app}$) da pesquisa MOJAVE.
   • Parâmetros do Motor Central: Massas de buracos negros ($M_{BH}$) e luminosidades do disco de acreção ($L_{acc}$) foram estimadas usando luminosidades de linhas largas (H$\beta$, MgII ou [OIII]) e valores de largura total a meia altura (FWHM). Taxas de acreção de massa normalizadas ($\dot{m}$) foram derivadas assumindo uma eficiência radiativa de $\eta = 0.4$.
4. Análise Estatística: Testes de correlação de posto de Kendall foram realizados para avaliar associações entre as escalas de tempo derivadas (durações de erupção, separações e tempos de subida) e parâmetros físicos ($M_{BH}$, $\dot{m}$, $L_{acc}$, $\Gamma$). Análises de correlação parcial foram conduzidas para contabilizar potenciais dependências do desvio para o vermelho ($z$).

Principais Resultados

• Escala de Tempo da PSD vs. Características da Curva de Luz: O inverso da frequência de dobra da PSD corresponde melhor à duração média das erupções mais brilhantes. Embora as separações médias das erupções também se correlacionem com a escala de tempo da PSD, as durações das erupções mostram um acordo um-para-um mais forte, particularmente para as 11 fontes restringidas.
• Conexão com o Motor Central: Usando durações e separações de erupção como proxies para a escala de tempo da PSD, os autores encontraram uma correlação positiva moderada com a razão entre a massa do buraco negro e a taxa de acreção de massa normalizada ($M_{BH}/\dot{m}$). Essa correlação é significativa no referencial do observador, mas enfraquece ou torna-se não significativa quando corrigida para o referencial de emissão. A correlação parece ser impulsionada principalmente por $M_{BH}$, embora o papel de $\dot{m}$ não possa ser definitivamente excluído devido às limitações do tamanho da amostra.
• Tamanho da Região de Emissão: Os tempos de subida das erupções mais brilhantes, usados como proxies para a largura da região de transição do jato, mostraram correlações tentativas fracas com $M_{BH}$, $L_{acc}$ e $\Gamma$. No entanto, essas associações foram encontradas como impulsionadas por uma dependência comum no desvio para o vermelho. Quando o desvio para o vermelho foi contabilizado via correlação parcial, a significância das relações entre tempos de subida e parâmetros do motor central desapareceu.
• Dependência do Desvio para o Vermelho: Uma correlação positiva significativa foi encontrada entre o desvio para o vermelho e as durações de erupção no referencial de emissão, enquanto uma correlação negativa existia entre o desvio para o vermelho e as separações de erupção no referencial de emissão (implicando uma taxa de ejeção maior para fontes de maior desvio para o vermelho). Essas tendências opostas podem explicar a falta de correlação entre durações/separações de erupção observadas e fatores Doppler.

Significado e Alegações
O artigo alega que as escalas de tempo características derivadas de dobras de PSD em rádio estão fisicamente associadas ao motor central de NGAs, escalonando especificamente com $M_{BH}/\dot{m}$. Isso apoia a hipótese de que o motor central modula a variabilidade do jato, semelhante a achados no domínio de raios-X. O estudo sugere que a duração média das erupções mais brilhantes é o proxy de curva de luz mais apropriado para a escala de tempo da dobra da PSD.

Os autores observam modestamente que, embora existam correlações com parâmetros do motor central, elas são frequentemente fracas e potencialmente confundidas por efeitos de desvio para o vermelho e vieses de amostra (por exemplo, a seleção de fontes com erupções discerníveis). Os resultados relativos ao tamanho da região de emissão e ao fator de Lorentz são descritos como "correlações tentativas fracas" que podem não se sustentar uma vez que as dependências do desvio para o vermelho forem removidas. Consequentemente, o artigo conclui que, embora uma conexão com o motor central seja sugerida, amostras maiores e estimativas de parâmetros mais precisas são necessárias para caracterizar definitivamente essas relações. O trabalho não propõe novos arranjos experimentais, mas sim valida o uso da decomposição de curvas de luz como um método para estender análises baseadas em PSD a uma amostra mais ampla de NGAs.