Photon-Count Statistics of Crab X-ray Pulses: Skellam Behavior and Excess Variance in the Main Pulse

Resumo Técnico: Estatísticas de Contagem de Fótons de Pulsos de Raios-X do Crabe

1. Problema e Contexto

O pulsar do Crabe (PSR B0531+21) é uma fonte extremamente brilhante que permite o estudo direto da variabilidade pulso-a-pulso em raios-X. Embora a distribuição de Poisson seja o padrão para contagem de fótons, a análise de dados de raios-X astronômicos frequentemente envolve a subtração de fundo (background). A diferença entre duas variáveis de Poisson independentes (sinal do pulso menos fundo) segue uma Distribuição de Skellam, e não uma distribuição Gaussiana ou Poissoniana simples.

Apesar de a distribuição de Skellam ser conhecida há décadas, sua aplicação explícita na análise de dados de raios-X é rara na literatura astrofísica, onde aproximações Gaussianas ou métodos baseados em Poisson para baixas contagens são comumente utilizados. O objetivo deste trabalho é caracterizar estatisticamente as distribuições de contagem de fótons dos componentes de pulso principal (MP - Main Pulse) e interpulso (IP - Interpulse) do Crabe, testando a adequação do modelo de Skellam e investigando a existência de variabilidade excessiva e correlações de curto prazo ("memória") entre pulsos consecutivos.

2. Metodologia

• Dados: Foram utilizados dados de raios-X do telescópio NICER (Neutron star Interior Composition Explorer), cobrindo um período de agosto de 2017 a janeiro de 2021, totalizando 78,8 ks de tempo de observação. O conjunto de dados contém mais de 2 milhões de pulsos individuais.
• Processamento de Dados:
  • Os tempos de chegada dos fótons foram corrigidos para o baricentro e sincronizados com o modelo de tempo do pulsar (usando ephemeris do Jodrell Bank e o software PINT).
  • Foram definidas janelas de fase para o MP (0,9 a 1,1), IP (0,3 a 0,45) e janelas de fundo ("OFF") de largura equivalente.
  • A contagem de fótons do pulso ($M$ e $I$) foi calculada subtraindo-se a contagem do fundo das contagens nas janelas "ON": $M = N_{MP} - N_{MP,OFF}$ e $I = N_{IP} - N_{IP,OFF}$.
• Análise Estatística:
  • As distribuições de contagem de fótons foram ajustadas a modelos Gaussianos e de Skellam.
  • A qualidade do ajuste foi avaliada usando o Qui-quadrado reduzido ($\chi^2_\nu$), o Critério de Informação de Akaike (AIC) e o Critério de Informação Bayesiano (BIC).
  • Foram testadas correlações de curto prazo ("efeito de memória") somando pares de pulsos consecutivos e calculando coeficientes de autocorrelação de lag-1.
  • Foram realizadas análises de sensibilidade dividindo as janelas de fase para investigar a influência de Pulsos de Rádio Gigantes (GRPs).

3. Principais Contribuições e Resultados

A. Comportamento de Skellam no Interpulso (IP)

• A distribuição de contagem de fótons do Interpulso (IP) é perfeitamente descrita pela distribuição de Skellam.
• O ajuste de Skellam apresentou um $\chi^2_\nu$ de 1,037 (98 graus de liberdade), indicando um ajuste excelente.
• Em contraste, o modelo Gaussiano falhou significativamente ($\chi^2_\nu = 1,716$).
• Significado: Este é um raro e robusto exemplo empírico de alta estatística demonstrando o comportamento de Skellam em um contexto astrofísico de contagem de fótons, validando o uso deste modelo para sinais subtraídos de fundo quando as contagens são suficientes.

B. Variabilidade Excessiva no Pulso Principal (MP)

• O Pulso Principal (MP) também é melhor descrito por Skellam do que por uma Gaussiana, mas apresenta um excesso de variância significativo.
• O ajuste de Skellam para o MP resultou em um $\chi^2_\nu$ alto de 1,57, indicando que uma distribuição de Skellam pura não captura toda a variabilidade observada.
• A análise de resíduos revela que o desvio ocorre principalmente na cauda de alta contagem (eventos com muitos fótons).
• Averaging (Média): Quando as contagens de fótons de dois pulsos consecutivos são somadas, a distribuição do MP torna-se consistente com o modelo de Skellam ($\chi^2_\nu \approx 1,006$). Isso indica que a variabilidade excessiva não persiste entre rotações, mas é um fenômeno localizado no pulso individual que se "suaviza" ao ser somado.

C. Ausência de Correlação de Memória

• A busca por correlações de curto prazo entre pulsos consecutivos (efeito de memória) não encontrou resultados estatisticamente significativos.
• O coeficiente de autocorrelação de lag-1 foi $r_1 = 0,002$, com um limite superior de confiança de 95% de $\rho_1 < 3,4 \times 10^{-3}$.
• Isso sugere que a variabilidade excessiva no MP não é causada por processos que persistem de um pulso para o seguinte.

D. Relação com Pulsos de Rádio Gigantes (GRPs)

• Embora os GRPs ocorram na janela de fase do MP e causem um aumento conhecido no fluxo de raios-X (~3,8%), simulações indicam que essa contribuição é insuficiente para explicar a magnitude da variância excessiva observada.
• A divisão da janela do MP em metades (antes e depois) sugeriu que a variabilidade excessiva pode estar ligeiramente concentrada na metade traseira do pulso, mas a análise foi limitada pela baixa contagem de fótons nessas sub-janelas.

4. Significado e Conclusões

• Distinção Estatística Clara: O trabalho estabelece uma distinção estatística clara entre o MP e o IP do Crabe. O IP segue um modelo puramente Poissoniano (diferença de Poissons = Skellam), enquanto o MP exibe uma física adicional de variabilidade de curta duração.
• Validação de Modelos: O estudo reforça a importância de utilizar modelos estatisticamente apropriados (Skellam) em vez de aproximações Gaussianas para dados de contagem de fótons subtraídos de fundo, especialmente em regimes de alta estatística.
• Implicações para Modelos de Emissão: A variabilidade excessiva no MP, que desaparece ao somar pulsos e não mostra correlação lag-1, sugere processos de emissão intrínsecos ao pulso, possivelmente relacionados a instabilidades rápidas ou fenômenos não estacionários dentro da janela de fase do MP, mas não ligados a GRPs de forma direta e dominante.
• Ferramenta para Astrofísica: Os resultados fornecem restrições quantitativas para modelos de emissão de pulsares e demonstram a utilidade de abordagens baseadas em Skellam para analisar dados de alta energia com grandes estatísticas de fótons.

Em suma, o artigo utiliza dados de alta precisão do NICER para demonstrar que, embora o comportamento estatístico fundamental do Crabe seja de Skellam, o pulso principal carrega uma assinatura de variabilidade intrínseca que desafia modelos puramente estocásticos simples, exigindo atenção cuidadosa na modelagem de dados de raios-X de alta energia.