The origin of the very-high-energy radiation along the jet of Centaurus A

Declaração do Problema
Centaurus A (Cen A), o núcleo galáctico ativo (AGN) mais próximo, exibe emissão de raios gama de muito alta energia (VHE > 240 GeV) detectada pelo observatório HESS, indicando a presença de elétrons ultra-relativísticos ao longo de seu jato. No entanto, o mecanismo de aceleração específico responsável por essas partículas permanece incerto. Embora vários autores tenham proposto que as subestruturas do jato, conhecidas como "nós" (knots), atuam como aceleradores eficientes de partículas, carecia-se de um modelo quantitativo que ligasse a hidrodinâmica desses nós ao espectro VHE observado. Especificamente, não era claro se a colisão do jato com ventos estelares (SW) poderia gerar as condições de choque necessárias para acelerar elétrons às energias requeridas e reproduzir a morfologia multi-comprimento de onda.

Metodologia
Os autores investigam a hipótese de que os nós estacionários no jato de Cen A resultam de interações entre o jato relativístico e os ventos estelares de estrelas poderosas (especificamente estrelas Wolf-Rayet). O estudo emprega uma abordagem de modelagem de múltiplas etapas:

1. Simulações Hidrodinâmicas Relativísticas (RHD): Usando o código PLUTO, os autores simulam a interação entre o jato (assumindo velocidade $v_j = 0,6c$, potência $L_j \sim 2 \times 10^{43}$ erg s$^{-1}$) e um vento estelar (caracterizado por $V_\infty = 3000$ km s$^{-1}$ e $\dot{M} = 10^{-4} M_\odot$ yr$^{-1}$). As simulações são realizadas em um referencial cilíndrico 2D axisimétrico para resolver a estrutura de choque duplo (choque de proa e choque de vento) e a descontinuidade de contato.
2. Aceleração e Resfriamento de Partículas: Com base na teoria da aceleração de choque difusiva (Fermi-I), os autores assumem que os elétrons são acelerados nas ondas de choque com um espectro de lei de potência e um corte exponencial em uma energia máxima ($E_{max}$). A energia máxima é determinada pelo equilíbrio entre o tempo de escala de aceleração e os tempos de escala de resfriamento por sincrotron e Compton Inverso (IC). O resfriamento depende dos campos magnéticos locais (derivados da pressão térmica via um parâmetro de equipartição $\beta_m$) e dos campos de radiação (SSC do núcleo, luz estelar, trilha de poeira, CMB e radiação de estrela WR).
3. Modelagem Radiativa: A Distribuição Espectral de Energia (SED) não térmica, do rádio aos raios gama, é calculada usando a biblioteca Naima. O modelo é restringido pelas medições de fluxo de rádio (VLA) e raios X (Chandra) e índices espectrais para nós específicos (AX1A, AX1C, AX2, BX2).
4. Comparação Morfológica: Imagens sintéticas de raios X são geradas projetando a emissividade simulada e aplicando um filtro Gaussiano para corresponder à resolução do Chandra, permitindo uma comparação direta com os tamanhos e formas dos nós observados.

Principais Contribuições

• Primeira Análise Combinada: Este trabalho representa a primeira tentativa de combinar simulações de fluidos, teoria de aceleração de choque e dados observacionais de múltiplos comprimentos de onda (rádio, raios X e raios gama VHE) para testar a hipótese de colisão jato-vento estelar para os nós de Cen A.
• Estrutura Hidrodinâmica: As simulações reproduzem com sucesso a estrutura de choque duplo esperada das interações jato-SW, fornecendo uma base física para a natureza estacionária dos nós observados.
• Modelagem de Emissão Estendida: Os autores vão além de modelar apenas os centroides dos nós para incluir a emissão estendida de toda a região perturbada/pós-choque, o que é crucial para explicar o fluxo VHE integrado.
• Restrições de Parâmetros: O estudo restringe parâmetros físicos fundamentais, incluindo a eficiência de aceleração ($\chi$) e a força do campo magnético ($\beta_m$), através do ajuste do modelo aos dados de rádio e raios X.

Resultos

• Energias de Elétrons: As simulações sugerem que os elétrons podem ser acelerados a energias de $\sim 3-4$ PeV nos nós AX2 e BX2, e até $\sim 6$ TeV nos nós internos AX1A/AX1C (assumindo menor eficiência de aceleração). Essas energias são suficientes para produzir os raios gama VHE observados via espalhamento IC.
• Distribuição Espectral de Energia: O modelo reproduz satisfatoriamente os fluxos de rádio e raios X para os nós AX2 e BX2. No entanto, para AX1A e AX1C, o modelo padrão superestima a emissão de raios X; uma eficiência de aceleração significativamente menor ($\chi \approx 2 \times 10^{-3}$) é necessária para coincidir com os dados, sugerindo que esses nós internos podem ter condições físicas diferentes (ex: choques mais relativísticos ou magnetização diferente).
• Origem de Raios Gama: A emissão dos centroides dos nós isoladamente é insuficiente para explicar os dados do HESS. Contudo, a emissão estendida de toda a região de choque (jusante do choque de proa) fornece uma explicação satisfatória para o fluxo de raios gama VHE detectado pelo HESS.
• Excesso do Fermi-LAT: O modelo sugere que o sinal do Fermi-LAT (faixa de GeV) pode ser explicado por elétrons que esfriaram conforme propagam ao longo do jato (até $\sim 3,6$ kpc) ou que estão presos nos grandes lobos, em vez de serem produzidos diretamente nos centroides dos nós.
• Morfologia: Os tamanhos de nós de raios X previstos para AX1C, AX2 e BX2 concordam com as medições do Chandra dentro de 5–25%. O nó AX1A é superestimado em tamanho, possivelmente indicando um vento estelar mais fraco ou mecanismos de formação alternativos.
• Potencial de UHECR: Embora os nós possam acelerar núcleos para $\sim 0,2 Z$ EeV em um cenário otimista de difusão de Bohm, os autores consideram isso provavelmente insuficiente para explicar diretamente os UHECRs de altíssima energia ($\sim 10^{19}$ eV). No entanto, os nós poderiam injetar núcleos pesados (de estrelas WR) que sofrem aceleração subsequente nos lobos.

Significância e Alegações
O artigo afirma que a aceleração de partículas nos locais de interação jato-vento estelar fornece uma explicação plausível para a emissão de rádio, raios X e raios gama VHE observada em Cen A. Os autores enfatizam que a emissão estendida das regiões de choque, e não apenas os centroides dos nós, é o contribuinte dominante para o fluxo VHE.

O estudo destaca que:

1. Os nós estacionários são provavelmente formados por colisões com estrelas Wolf-Rayet.
2. O espectro de raios gama VHE é consistente com elétrons acelerados nesses choques e resfriados via espalhamento IC em campos de radiação de poeira e luz estelar.
3. Observações futuras com o Cherenkov Telescope Array (CTAO) poderão resolver espacialmente as contribuições de diferentes complexos de nós (A, B, C, E vs. F, G) na faixa de 10–100 TeV, oferecendo um teste direto da eficiência de aceleração e dos nós dominantes na banda VHE.

Os autores mantêm-se modestos quanto aos Raios Cósmicos de Ultra-Alta Energia (UHECRs), observando que, embora os nós possam injetar núcleos pesados, a atividade atual do jato provavelmente não pode acelerá-los às energias mais altas observadas sem uma re-aceleração adicional nos lobos. Da mesma forma, reconhecem que efeitos de difusão, não totalmente capturados em seu modelo dominado por advecção, podem ser necessários para explicar plenamente a contagem de fótons e os raios dos nós.