NICER Perspective on TeV Blazar Mrk~421: X-ray Variability and Particle Acceleration

Resumo Técnico: Perspectiva do NICER sobre o Blazar de TeV Mrk 421: Variabilidade de Raios-X e Aceleração de Partículas

Problema e Motivação
Blazares, particularmente objetos BL Lac de Pico de Síncrotron Alto (HSP) como Markarian 421 (Mrk 421), exibem variabilidade significativa em todo o espectro eletromagnético, impulsionada por jatos relativísticos orientados ao longo da linha de visada do observador. Embora o Mrk 421 seja um emissor de TeV prototípico e um dos blazars mais próximos ($z=0,031$), a compreensão detalhada da evolução espectral de raios-X e dos mecanismos físicos que impulsionam sua rápida variabilidade permanece uma área ativa de pesquisa. Estudos anteriores utilizaram vários instrumentos (ex: Swift, XMM-Newton, NuSTAR), mas uma análise abrangente do comportamento de raios-X de Mrk 421 utilizando as capacidades de alta sensibilidade e alta resolução temporal do Neutron Star Interior Composition Explorer (NICER) ao longo de um período sustentado era carente. Este estudo aborda a necessidade de caracterizar a variabilidade de longo prazo e intradia, a curvatura espectral e a correlação entre parâmetros espectrais para restringir os mecanismos de aceleração de partículas dentro do jato.

Metodologia
Os autores analisaram 45 observações do NICER de Mrk 421 coletadas entre 2022 e 2024. A redução de dados utilizou HEASoft v6.34 e XSPEC v12.14.1, restringindo a faixa de energia para 0,4–10 keV para otimizar a calibração e minimizar o ruído de fundo.

• Análise de Tempo: As curvas de luz foram extraídas com um tamanho de bin de 60 segundos. O estudo computou a variabilidade de Raiz do Quadrado Médio (RMS), a variabilidade fracionária ($F_{var}$) e examinou a distribuição de fluxo para determinar se seguia perfis Gaussianos ou log-normais.
• Análise Espectral: Os espectros foram ajustados usando três modelos: Lei de Potência (PL), Lei de Potência Quebrada (BPL) e Log-Parábola (LP). A densidade de coluna de hidrogênio galáctico foi fixada em $N_H = 0,0134 \times 10^{22} \text{ cm}^{-2}$. A seleção de modelo foi realizada usando o qui-quadrado reduzido ($\chi^2_r$) e o teste F para determinar a significância estatística.
• Estudos de Correlação: Os autores investigaram as relações entre parâmetros espectrais (índice de fótons $\Gamma$ ou $\alpha$, curvatura $\beta$, energia de pico de síncrotron $E_p$) e o fluxo. Razões de dureza (HR) foram calculadas para avaliar mudanças na forma espectral.
• Modelagem Física: Para interpretar a curvatura espectral, os autores simularam espectros de síncrotron usando uma distribuição de energia de elétrons log-parabólica ($N(\gamma) \propto \gamma^{-p-r \log_{10}(\gamma)}$) para reproduzir a anti-correlação observada entre $E_p$ e o parâmetro de curvatura.

Principais Resultados

1. Variabilidade: O Mrk 421 exibiu forte variabilidade de raios-X ao longo do período de dois anos, com níveis de fluxo variando por um fator de $\sim28$ (de $\sim50$ para $\sim1380$ contagens s$^{-1}$). A variabilidade fracionária total em todas as observações foi de aproximadamente 48%. A variabilidade intradia também foi significativa, com mudanças ocorrendo em escalas de tempo de algumas horas.
2. Distribuição de Fluxo: A distribuição das taxas de contagem de raios-X foi multimodal, sendo melhor ajustada por uma superposição de três componentes Gaussianas correspondentes a estados de fluxo baixo, intermediário e alto. Isso sugere que a variabilidade não é impulsionada por um único processo estocástico, mas provavelmente envolve múltiplos estados ou regiões de emissão.
3. Modelagem Espectral: O modelo Log-Parábola (LP) forneceu a descrição mais precisa para 42 das 45 observações, superando significativamente o modelo de Lei de Potência. As três observações restantes foram melhor descritas por uma Lei de Potência Quebrada. Os parâmetros espectrais médios para o modelo LP foram $\alpha \approx 2,32$ e $\beta \approx 0,32$.
4. Tendências Espectrais:
   • Mais duro quando mais brilhante (Harder-when-brighter): Uma forte correlação positiva foi encontrada entre a Razão de Dureza (HR) e a taxa de contagem ($r=0,736$), indicando que o espectro endurece conforme a fonte brilha.
   • Correlações de Parâmetros: Uma correlação positiva foi observada entre o índice de fótons $\alpha$ e o parâmetro de curvatura $\beta$ ($r=0,72$). Inversamente, $\beta$ mostrou uma correlação negativa com a energia de pico de síncrotron $E_p$ ($r=-0,52$), e $E_p$ mostrou uma correlação positiva com o fluxo ($r=0,82$).
   • Energia de Pico vs. Curvatura: A análise confirmou que, à medida que o pico de síncrotron se desloca para energias mais altas (estados mais brilhantes), a curvatura espectral ($\beta$) diminui.
5. Restrições da Região de Emissão: Com base nas escalas de tempo de variabilidade rápida ($\sim40$ min) e assumindo o resfriamento por síncrotron, a força do campo magnético foi estimada em $B \sim 0,5$ G. Argumentos de causalidade restringiram o tamanho da região de emissão a $R \lesssim 7,0 \times 10^{14}$ cm, sugerindo que os raios-X se originam de uma região compacta próxima ao motor central (algumas centenas de raios de Schwarzschild).

Significância e Alegações
O artigo afirma apresentar a primeira análise detalhada de tempo e espectro de Mrk 421 baseada em dados do NICER. A principal significância reside em confirmar que os espectros de raios-X deste blazar HSP são consistentemente melhor descritos por um modelo log-parabólico, apoiando cenários de aceleração de partículas estocástica ou dependente de energia.

A anti-correlação observada entre a energia de pico de síncrotron ($E_p$) e o parâmetro de curvatura ($\beta$) é interpretada dentro do arcabouço da aceleração de partículas dependente de energia (EDAP) e aceleração estocástica (SA). Nesses modelos, a probabilidade de aceleração de partículas diminui conforme a energia da partícula aumenta, produzindo naturalmente uma distribuição de energia de elétrons log-parabólica. O estudo conclui que a evolução espectral e a variabilidade de Mrk 421 são moldadas por processos não térmicos complexos envolvendo turbulência, campos magnéticos fortes e fluxos relativísticos, onde a curvatura do espectro serve como um diagnóstico para o mecanismo de aceleração subjacente. A distribuição de fluxo multimodal implica ainda que o jato pode hospedar múltiplas zonas de emissão distintas, em vez de uma única zona homogênea.