Determination of neutron star radius from pulse profile modeling using profile likelihood

Resumo Técnico: Determinação do Raio de Estrela de Neutrons através da Modelagem de Perfil de Pulso Usando Verossimilhança de Perfil

Definição do Problema
O Neutron Star Interior Composition Explorer (NICER) utiliza a modelagem de perfil de pulso de raios-X para determinar a massa e o raio de estrelas de nêutrons, empregando principalmente o pacote X-ray Pulse Simulation and Inference (X-PSI). A prática padrão atual baseia-se na inferência Bayesiana utilizando algoritmos de amostragem aninhada (nested sampling, como o MultiNest) para estimar distribuições posteriores. No entanto, discrepâncias nos raios inferidos foram observadas em diferentes estudos, possivelmente atribuídas à escolha do amostrador ou a problemas de convergência, particularmente em superfícies de verossimilhança planas para raios superiores a 11 km. Além disso, a marginalização Bayesiana sobre amplos intervalos de prior pode introduzir efeitos de volume que enviesam as estimativas dos parâmetros. Embora os métodos Bayesianos sejam o padrão astrofísico, a inferência frequentista permanece dominante na física de partículas. Este trabalho aborda a necessidade de comparar as abordagens Bayesiana e Frequentista para garantir a robustez das restrições de física fundamental derivadas do X-PSI.

Metodologia
Os autores implementaram um pipeline de inferência frequentista utilizando a técnica de verossimilhança de perfil (profile likelihood) dentro do atual quadro do X-PSI, garantindo que o modelo físico subjacente e a implementação da verossimilhança permanecessem idênticos à abordagem Bayesiana padrão.

• Fonte de Dados: A análise utilizou um conjunto de dados sintéticos NICER + XMM-Newton (designado como syntX1), anteriormente utilizado em H25 para testar amostradores aninhados.
• Configuração do Modelo: O estudo empregou o modelo de perfil de pulso ST-U, apresentando dois pontos quentes circulares, de temperatura única e simplesmente conectados. A atmosfera da estrela de nêutrons foi modelada usando o modelo de hidrogênio totalmente ionizado NSX, com a propagação de fótons calculada via traçado de raios em um espaço-tempo de Schwarzschild oblato de relatividade geral.
• Implementação Frequentista:
  • O vetor de parâmetros foi dividido entre o parâmetro de interesse (raio equatorial, $R_{eq}$) e os parâmetros de confusão/nuance ($\alpha$).
  • Os parâmetros de nuance (incluindo massa, distância, inclinação, geometria do ponto quente, temperaturas e densidade de coluna de hidrogênio interestelar) foram perfilados através da maximização da verossimilhança $L(R_{eq}, \alpha)$ em relação a $\alpha$ para valores fixos de $R_{eq}$.
  • A otimização foi realizada utilizando o algoritmo simplex Nelder-Mead (scipy.optimize.minimize).
  • Os intervalos de confiança foram construídos utilizando a estatística $\Delta\chi^2$ ($\Delta\chi^2 = -2 \ln L_{profile} + 2 \ln L_{max}$), assumindo uma distribuição $\chi^2$ para um grau de liberdade.
• Linha de Base de Comparação: Os resultados foram comparados contra uma análise Bayesiana padrão utilizando o algoritmo de amostragem aninhada MultiNest com configurações idênticas e configurações de baixa resolução do X-PSI.

Principais Contribuições

1. Aplicação de Prova de Conceito: Este trabalho demonstra a primeira aplicação da inferência de verossimilhança de perfil frequentista para determinar raios de estrelas de nêutrons utilizando o pacote X-PSI.
2. Comparação Metodológica: Fornece uma comparação direta entre os pipelines de inferência Bayesiana e frequentista em dados sintéticos idênticos, isolando a metodologia de inferência como a única variável.
3. Ciência Aberta: Os autores liberaram publicamente as modificações de código de inferência frequentista para o X-PSI para facilitar pesquisas futuras.

Resultados
Utilizando o conjunto de dados sintéticos syntX1 com um raio injetado verdadeiro de 12,176 km:

• Inferência Bayesiana: O intervalo de credibilidade de 68% marginalizado para o raio foi $R_{68,bayes} = 11,970^{+0,180}_{-0,181}$ km. Este resultado concordou com o valor injetado dentro de $1,1\sigma$.
• Inferência Frequentista: A verossimilhança de perfil resultou em um raio de melhor ajuste de 12,096 km. O intervalo de confiança de 68% foi $r_{68} = 12,096^{+0,064}_{-0,182}$ km, e o intervalo de 95% foi $r_{95} = 12,096^{+0,316}_{-0,397}$ km. A estimativa frequentista concordou com o valor injetado dentro de $0,4\sigma$.
• Precisão: A precisão do método frequentista mostrou-se comparável à análise Bayesiana.
• Eficiência Computacional: A análise frequentista exigiu aproximadamente 9 horas de núcleo de CPU, enquanto a análise Bayesiana exigiu aproximadamente 3.840 horas de núcleo de CPU (80 horas em um nó de 48 núcleos). O método frequentista é, portanto, ordens de magnitude mais rápido.

Significância e Alegações
O artigo posiciona este trabalho como um estudo de prova de conceito. Os autores alegam que a inferência frequentista via verossimilhança de perfil consegue recuperar o raio verdadeiro de uma estrela de nêutrons com uma precisão de $< 1\sigma$ usando dados sintéticos. A técnica oferece precisão comparável à análise Bayesiana, mas com um custo computacional significativamente reduzido. Os autores concluem que este método serve como um complemento valioso à inferência Bayesiana, oferecendo uma alternativa robusta para a estimativa de parâmetros no X-PSI. Eles declaram explicitamente que trabalhos futuros envolverão a aplicação desta técnica a dados reais de pulsares do telescópio NICER, mas nenhuma tal aplicação é apresentada neste manuscrito.