Multiwavelength variability of the high-energy neutrino candidate PKS 0735+178 over three decades

Resumo Técnico: Variabilidade Multicomprimento de Onda do Candidato a Neutrino de Alta Energia PKS 0735+178 ao Longo de Três Décadas

Problema
PKS 0735+178 é um objeto BL Lac clássico associado ao evento de neutrino de alta energia IC211208A (171 TeV). Embora a fonte seja conhecida por sua extrema variabilidade multicomprimento de onda e um jato relativístico visualizado sob pequenos ângulos, a conexão física entre as regiões de emissão em diferentes frequências, a natureza de seu comportamento quase periódico de longo prazo e os mecanismos específicos que impulsionam a atividade em torno do evento de neutrino permanecem mal compreendidos. Estudos anteriores focaram primariamente em propriedades espectrais de curto prazo ou em flares isolados, carecendo de uma análise unificada da evolução temporal da fonte ao longo de uma base de várias décadas para restringir a geometria do jato e a propagação de perturbações.

Metodologia
Os autores apresentam uma análise abrangente de aproximadamente 30 anos (1995–2026) de dados multicomprimento de onda, integrando:

• Dados de Rádio: Observações quase simultâneas de 1–22 GHz do telescópio RATAN-600 e dados de 230 GHz do Submillimeter Array (SMA).
• Dados Ópticos: Fotometria na banda R proveniente dos telescópios Zeiss-1000, AS-500/2, KAIT, ZTF e CRTS.
• Dados de Raios Gama: Curvas de luz do Fermi-LAT (1.5–300 GeV) construídas usando agrupamento adaptativo (adaptive binning) para resolver estruturas finas durante estados de flare.
• Restrições de VLBI: Medições de core-shift a 15 e 43 GHz usando dados do MOJAVE e VLBA-BU-BLAZAR para derivar a cinemática do jato e parâmetros físicos.

A análise empregou diversas técnicas de modelagem estatística e física:

1. Correlação Cruzada: Funções de Correlação Discreta (DCF) foram usadas para medir atrasos temporais entre as bandas.
2. Análise Espectral: Espectros de rádio foram agrupados em três grandes episódios de flare para rastrear a evolução do índice espectral ($\alpha$).
3. Busca de Quasi-Periodicidade: Os periodogramas de Transformada de Wavelet Ponderada (WWZ) e de Lomb-Scargle Generalizado (GLS) foram aplicados para detectar escalas de tempo características, com a significância avaliada contra fundos de ruído vermelho via simulações de Monte Carlo.
4. Modelagem Física:
   • Propagação de Choque: Atrasos temporais foram ajustados a uma dependência de lei de potência ($\Delta t \propto \nu^{-1/k_r}$) para testar um modelo de jato estratificado por opacidade.
   • Precessão de Jato: Um modelo de otimização Bayesiana (UltraNest) foi usado para ajustar as curvas de luz a um cenário de precessão, testando se efeitos geométricos poderiam explicar a modulação de longo prazo.
   • Core Shift de VLBI: Usado para estimar o ângulo de visão do jato, velocidade aparente, força do campo magnético e fluxo magnético em relação ao limite de Disco Magneticamente Arrestado (MAD).

Principais Resultados

• Atrasos Temporais Dependentes da Frequência: As curvas de luz em radio, óptico e raios gama estão correlacionadas, mas a emissão de rádio atrasa em relação às frequências mais altas. Os atrasos aumentam sistematicamente em direção a frequências mais baixas, atingindo $\sim$1200 dias em 1–2 GHz em relação aos raios gama. Isso segue uma lei de potência com $k_r = 1.96 \pm 0.24$, consistente com um jato estratificado por opacidade onde as perturbações propagam-se rio abaixo.
• Flares do Evento de Neutrino: Flares curtos e brilhantes nas bandas óptica e de raios gama ocorreram dentro de dias do evento IC211208A. Um atraso robusto de 12 dias foi detectado entre os picos óptico e de raios gama (óptico liderando), sugerindo regiões de emissão não co-espaciais ou contribuições de processos como o espalhamento Compton externo. O terceiro grande flare de rádio (2019–2026) exibe espectros invertidos ou planos em baixas frequências, indicando uma região de emissão compacta e opticamente espessa com aceleração de partículas intensificada.
• Variabilidade Quasi-Periódica: Uma escala de tempo de variabilidade característica de $\sim$10–11 anos é robustamente detectada ($\ge 3\sigma$) nas bandas de rádio e milimétrica (5–230 GHz). Esta modulação é coerente através dessas bandas, mas não é estritamente periódica. Em contraste, os dados ópticos e de raios gama mostram sinais mais fracos e de períodos mais curtos ($\sim$5–9 anos) que são localizados no tempo e menos coerentes.
• Geometria e Física do Jato:
  • Medições de core-shift de VLBI resultam em um ângulo de visão do jato de $\theta \approx 3.6^\circ$ e uma velocidade aparente de $\beta_{app} \approx 0.8c$, consistentes com as velocidades de propagação derivadas dos atrasos temporais.
  • O campo magnético a 1 pc é estimado em $\sim$0.36 G. A razão do fluxo magnético em relação ao limite MAD sugere que o disco de acreção não está em um estado MAD nos níveis atuais de luminosidade.
• Limitações do Modelo de Precessão: Embora um modelo de precessão de jato possa reproduzir a forma geral das curvas de luz, os parâmetros derivados (ângulo de visão $\sim$27.4°) contradizem as restrições de VLBI ($\theta \sim 1^\circ-3^\circ$). Isso indica que a precessão sozinha não pode explicar a variabilidade observada.

Significância e Alegações
O artigo alega que a variabilidade de PKS 0735+178 é governada por dois processos físicos distintos operando em diferentes escalas de tempo:

1. Modulação de Longo Prazo: Uma variação lenta na liberação de energia na base do jato impulsiona uma modulação quasi-periódica ($\sim$11 anos), provavelmente ligada ao fluxo de acreção ou à dinâmica da base do jato, que regula a taxa de injeção de perturbações.
2. Flares de Curto Prazo: Flares individuais, incluindo o associado ao evento de neutrino, surgem de choques propagando-se rio abaixo através de um jato estratificado por opacidade. Os atrasos dependentes da frequência e a evolução espectral são naturalmente explicados pelo choque cruzando fotossferas de sincrotron em diferentes distâncias.

O estudo fornece a primeira medição confiável do atraso de 12 dias entre as emissões óptica e de raios gama nesta fonte, desafiando cenários simples de synchrotron self-Compton de zona única. Além disso, a consistência entre as restrições independentes de core-shifts de VLBI e os atrasos temporais de rádio apoia fortemente um cenário unificado de propagação de choque em um fluxo estratificado. O agrupamento temporal do evento de neutrino com o flare multicomprimento de onda reforça a conexão entre PKS 0735+178 e a produção de neutrinos de alta energia, provavelmente ocorrendo nas regiões internas e opticamente espessas do jato.