Study of an interaction between the jet and an interstellar medium of M87 with a spectral analysis by using Chandra

Resumo Técnico: Estudo de uma Interação entre o Jato e um Meio Interestelar de M87 com uma Análise Espectral utilizando o Chandra

Definição do Problema
A galáxia de rádio M87, localizada no centro do aglomerado de Virgem, exibe um jato proeminente que se estende por 20 segundos de arco. Embora a origem dos raios gama de TeV em M87 seja geralmente atribuída ao núcleo ou a nós específicos como o HST-1, os mecanismos físicos que impulsionam a interação entre o jato e o meio interestelar (ISM) circundante permanecem sob investigação. Análises espectrais anteriores do jato de M87 usando dados do Chandra (ex: Wilson & Yang 2002; Perlman & Wilson 2005; Sun et al. 2018) focaram primordialmente no ajuste de modelos de lei de potência aos espectros de energia de raios X, assumindo emissão síncrotron. No entanto, esses estudos enfrentaram limitações:

1. Estatísticas Insuficientes: Muitas análises basearam-se em tempos de exposição curtos (ex: 13–38 ks), levando a estatísticas pobres onde modelos complexos poderiam parecer aceitáveis simplesmente devido à baixa qualidade dos dados.
2. Erros de Subtração de Background: Trabalhos anteriores frequentemente utilizaram uma região de background comum para todo o jato. Isso é fisicamente incorreto porque o gás quente do aglomerado de Virgem, que contribui para o background, varia em intensidade com a distância do núcleo.
3. Negligência de Componentes Térmicos: Existe a possibilidade de que o jato interaja com o ISM via choques, aquecendo o gás para produzir um componente térmico (bremsstrahlung) além da emissão síncrotron não térmica. Este componente pode ter sido omitido em estudos de baixa estatística.

Metodologia
Para abordar essas questões, este estudo utiliza um conjunto de dados significativamente maior do Observatório de Raios X Chandra, cobrindo observações de 2000 a 2018 com um tempo de exposição total de aproximadamente 800 ks.

• Seleção de Dados e Tratamento de Pile-up: A análise foca em quatro regiões específicas: o núcleo, o nó HST-1, o nó D e o nó A. Devido ao alto brilho do núcleo e de HST-1, os dados foram cuidadosamente selecionados para gerenciar o "pile-up" (onde múltiplos fótons atingem um único pixel, distorcendo o espectro). Os dados foram filtrados para garantir que as frações de pile-up estivessem abaixo de 5% para a maioria das análises, embora 10% tenha sido aceito para o núcleo e HST-1 no conjunto de dados anterior para maximizar as estatísticas.
• Modelagem de Background: Diferente de estudos anteriores, as regiões de background foram selecionadas individualmente para cada nó (norte e sul da fonte) para considerar a variação espacial do gás quente do aglomerado de Virgem. Para o nó A, uma região de background específica foi escolhida para evitar quedas conhecidas de raios X moles.
• Análise Espectral: Os espectros de energia foram extraídos usando o software CIAO (v4.15) e analisados com o XSPEC. A faixa de energia foi definida de 0,5 keV a 10 keV para evitar problemas de calibração em energias mais baixas.
• Ajuste de Modelo:
  1. Um modelo simples de lei de potência (representando a emissão síncrotron) com absorção fotoelétrica foi testado.
  2. Se o ajuste da lei de potência fosse estatisticamente rejeitado (indicado por um $\chi^2$ alto e excesso mole), um modelo combinado de lei de potência e um modelo APEC (bremsstrahlung térmico) era aplicado.
  3. Parâmetros físicos, incluindo densidade neutra, densidade de plasma e densidade de elétrons acelerados, foram derivados dos melhores modelos de ajuste.
  4. O fluxo de bremsstrahlung não térmico resultante da interação de elétrons acelerados com o ISM foi calculado e comparado com observações de raios gama do Fermi.

Resultos Principais

1. Características Espectrais:

   • Núcleo e Nó D: Os espectros de energia de raios X tanto para o núcleo quanto para o nó D são bem descritos por um modelo simples de lei de potência, consistente com a emissão síncrotron.
   • Nó HST-1 e Nó A: Para estas regiões, um modelo de lei de potência pura é estatisticamente rejeitado.
     • HST-1: No conjunto de dados combinado (incluindo dados com 10% de pile-up), um modelo de lei de potência resulta em uma probabilidade de chance de $2,6 \times 10^{-4}$. Adicionar um componente térmico (APEC) fornece um ajuste razoável.
     • Nó A: O modelo de lei de potência é fortemente rejeitado ($\chi^2 = 628,76$ para 355 d.o.f, probabilidade $1,6 \times 10^{-17}$). Um modelo combinado de lei de potência e bremsstrahlung térmico ajusta bem os dados.

2. Propriedades do Componente Térmico:

   • O componente térmico detectado nos nós HST-1 e A é caracterizado por uma temperatura de 0,2 keV e uma abundância metálica de 0,00 (essencialmente zero).
   • A abundância metálica é significativamente menor que a do gás quente circundante do aglomerado de Virgem (que é de ~1 solar), sugerindo que o gás aquecido origina-se dentro do próprio jato, e não do meio ambiente.
   • A razão de fluxo térmico no nó A é de aproximadamente 22% do fluxo total.

3. Parâmetros Físicos:

   • Densidade de Plasma: Derivada da normalização APEC, a densidade eletrônica no nó A é de aproximadamente $11,0 \, \text{cm}^{-3}$.
   • Aquecimento por Choque: A presença de um componente térmico com baixa metalicidade e alta temperatura em relação ao meio ambiente apoia a hipótese de que o gás dentro do jato é aquecido por uma onda de choque resultante da interação com o ISM.
   • Densidade Neutra: Calculada a partir da diferença entre a coluna de densidade derivada de raios X e observações de rádio de 21 cm, indicando a presença de material frio ao redor do jato.

4. Bremsstrahlung Não Térmico e Comparação com Fermi:

   • O estudo calculou o fluxo esperado de bremsstrahlung não térmico de elétrons acelerados interagindo com o ISM.
   • O fluxo calculado do núcleo e do nó HST-1 é ordens de magnitude inferior ao fluxo observado em GeV pelo Fermi.
   • Conclusão: Não há contribuição do bremsstrahlung não térmico no núcleo ou no HST-1 para o fluxo observado de GeV de M87. O artigo observa que o nó A pode contribuir se seu campo magnético for significativamente menor do que o assumido (<10 $\mu$G), mas isso requer mais estudos.

Significância e Alegações
O artigo alega que, ao utilizar dados de alta estatística (800 ks) e uma subtração de background espacialmente dependente e correta, identificou com sucesso um componente térmico nos espectros de raios X dos nós do jato de M87, HST-1 e A. Esta descoberta desafia interpretações anteriores que dependiam exclusivamente de modelos de lei de potência.

Os autores concluem que:

1. A emissão térmica é provavelmente gás dentro do jato aquecido por uma interação de choque com o ISM, evidenciado pela baixa abundância metálica e alta temperatura.
2. O mecanismo de bremsstrahlung não térmico (interação elétron-ISM) não explica o fluxo observado de raios gama em GeV de M87, reforçando que a origem da emissão de TeV/GeV provavelmente reside em outro lugar (ex: o núcleo ou diferentes mecanismos de aceleração de partículas).
3. Estudos anteriores podem ter perdido essas características térmicas devido a estatísticas insuficientes e modelagem inadequada de background.

O estudo não propõe novos setups experimentais, mas enfatiza a necessidade de análise de dados de arquivo de alta qualidade para resolver o espectro de energia multicomprimento e as condições físicas de jatos relativísticos.