Multi-epoch VLBI observations of the blazar 3C 66A: Spatial twisting and temporal oscillation of the parsec-scale jet

Enunciado do Problema
Os blazares, uma subclasse de núcleos ativos de galáxias (AGNs) com jatos alinhados à linha de visão da Terra, frequentemente exibem comportamentos cinemáticos complexos, incluindo variações no ângulo de posição (PA) do jato interno. Embora tais variações de PA tenham sido detectadas em numerosos blazares, suas origens físicas permanecem debatidas. Os mecanismos propostos incluem instabilidades no fluxo do jato, precessão do disco, efeitos de Lense-Thirring (LT) e sistemas de buracos negros supermassivos binários (SMBHB). O blazar 3C 66A é um objeto bem conhecido com uma estrutura de jato torcida e uma periodicidade conhecida no fluxo óptico (~2 anos), no entanto, a cinemática do jato em escala de parsec e a relação entre suas mudanças de orientação e a variabilidade do fluxo não foram totalmente caracterizadas ao longo de longos períodos de tempo. Estudos anteriores utilizando dados do VLBA cobriram intervalos de tempo limitados (tipicamente <15 anos) ou resoluções menores, potencialmente perdendo periodicidades de longo prazo ou oscilações sutis.

Metodologia
Os autores realizaram uma análise multi-época combinando novas observações de alta resolução com extensos dados de arquivo:

1. Novas Observações: Observações do KaVA (KVN e VERA Array) foram realizadas a 22 e 43 GHz em 2014 através de três épocas para mapear a morfologia do jato em escala de parsec.
2. Dados de Arquivo: O estudo utilizou 109 épocas de dados de arquivo de 43 GHz do VLBA abrangendo de 1996 a 2025, incluindo dados dos programas de monitoramento VLBA-BU-BLAZAR e BEAM-ME.
3. Redução de Dados: Técnicas padrão de calibração e imagem (AIPS e DIFMAP) foram empregadas. A estrutura do jato foi parametrizada usando ajuste de modelo Gaussiano circular em visibilidades UV.
4. Análise:
   • Morfologia: A estrutura do jato foi mapeada para identificar curvaturas e torções.
   • Cinemática: O PA do jato interno foi rastreado usando dois métodos: ajuste de componente Gaussiana e médias ponderadas pelo fluxo de componentes clean em bins anulares.
   • Busca de Periodicidade: Três métodos (Jurkevich, Lomb-Scargle e Transformada de Wavelet Ponderada de Z) foram aplicados à série temporal do PA para detectar periodicidade.
   • Correlação: Uma Função de Correlação Discreta (DCF) foi usada para testar correlações entre o PA do jato interno e a densidade de fluxo do núcleo.
   • Modelagem: Um modelo de jato precessante foi ajustado aos dados de PA usando análise Bayesiana (MCMC) para restringir parâmetros geométricos.
   • Cenários Físicos: Os autores avaliaram vários modelos físicos, incluindo movimento orbital de SMBHB, precessão spin-órbita, precessão disco-órbita, efeitos de LT e instabilidades de jato (CDI, KHI), para explicar as escalas de tempo observadas.

Principais Contribuições e Resultados

• Morfologia de Jato Torcido: Mapas de alta resolução do KaVA confirmaram uma estrutura de jato torcida com múltiplas curvaturas (sudoeste para sudeste, depois sul, depois sudeste) estendendo-se até ~16 mas.
• Oscilação e Deslocamento de PA: Pela primeira vez, uma oscilação periódica no PA do jato interno foi detectada com uma amplitude de $7,77 \pm 0,79^\circ$ e um período de $10,94 \pm 0,22$ anos. Sobreposta a essa oscilação, há um deslocamento contínuo no sentido horário de $-0,83 \pm 0,07^\circ$/ano.
• Propagação Jusante: Uma oscilação semelhante foi detectada no jato jusante (a ~0,7 mas), mas com amplitude reduzida ($4,25^\circ$) e um deslocamento de fase de ~1,5 anos, sugerindo uma propagação do tipo onda.
• Variabilidade de Fluxo: A densidade de fluxo do núcleo exibe forte variabilidade com possíveis periodicidades em ~10, ~7 e ~3-4 anos. Uma correlação fraca, mas significativa ($2\sigma$), foi encontrada entre o PA do jato interno e a densidade de fluxo do núcleo, com mudanças de PA precedendo as mudanças de fluxo em ~2 anos.
• Linha de Base de Longo Prazo: Uma diminuição contínua na densidade de fluxo do núcleo de rádio de $-9,40 \pm 1,43$ mJy/ano foi observada ao longo do intervalo de 29 anos.
• Interpretação SMBHB: Os autores propõem que as características observadas são melhor explicadas por um sistema SMBHB.
  • A oscilação de PA de ~11 anos é atribuída ao movimento orbital do binário (período intrínseco de ~8,12 anos) combinado com a precessão disco-órbita ("wobbling").
  • O deslocamento de PA de longo prazo e a diminuição da linha de base de fluxo são atribuídos à precessão disco-órbita impulsionada pelo desalinhamento entre o disco de acreção e o plano orbital do binário.
  • Parâmetros derivados para o hipotético SMBHB: Massa total $M \approx 1,42 \times 10^8 M_\odot$, separação $r \approx 1,65 \times 10^{-2}$ pc, e um período orbital de ~16,7 anos (referencial).
• Modelos Alternativos: O artigo descarta ou desfavorece sistematicamente outros mecanismos:
  • Precessão LT: Pode explicar tanto o período de curto quanto o de longo prazo, mas não ambos simultaneamente com um único conjunto de parâmetros físicos.
  • Instabilidades de Jato (CDI/KHI): Embora instabilidades de kink induzidas por corrente (CDI) ou Kelvin-Helmholtz (KHI) possam explicar o padrão ondulatório de ~11 anos, elas não explicam o desvio secular de longo prazo ou a periodicidade do fluxo de forma tão abrangente quanto o modelo SMBHB.
  • Instabilidades de Disco: Mecanismos como empenamento induzido por radiação ou instabilidades térmico-viscosas produzem escalas de tempo ordens de magnitude superiores para coincidir com as observações.

Significância
O artigo afirma que a combinação de um jato torcido, uma oscilação de PA de 11 anos, um deslocamento de PA secular de longo prazo e variabilidade de fluxo correlacionada fornece evidências fortes para um jato precessante impulsionado por um sistema SMBHB. Este cenário naturalmente unifica as escalas de tempo de curto e longo prazo, bem como a periodicidade do fluxo óptico (2 anos, interpretada como emissão de jato com Doppler-boost), sem exigir parâmetros ajustados para mecanismos de buraco negro único. A detecção da oscilação de PA ao longo de quase três ciclos (30 anos) em alta resolução (43 GHz) representa um avanço significativo na compreensão da cinemática de jatos em blazars. Os autores concluem que, embora cenários alternativos envolvendo instabilidades ou efeitos de LT não possam ser inteiramente descartados para características específicas, o modelo SMBHB oferece a explicação mais parcimoniosa para o conjunto completo de propriedades observadas em 3C 66A. O sistema é previsto para fundir-se em ~934 Myr, embora sua deformação de onda gravitacional esteja atualmente abaixo da sensibilidade dos arrays de temporização de pulsares.