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

Résumé Technique : Détermination du rayon d'une étoile à neutrons par modélisation du profil d'impulsion à l'aide de la vraisemblance profilée

Énoncé du problème
Le Neutron Star Interior Composition Explorer (NICER) utilise la modélisation du profil d'impulsion de rayons X pour déterminer la masse et le rayon des étoiles à neutrons, en employant principalement le package X-ray Pulse Simulation and Inference (X-PSI). La pratique standard actuelle repose sur l'inférence bayésienne utilisant des algorithmes d'échantillonnage imbriqué (par exemple, MultiNest) pour estimer les distributions a posteriori. Cependant, des écarts dans les rayons inférés ont été observés entre différentes études, ce qui pourrait être attribué au choix de l'échantillonneur ou à des problèmes de convergence, particulièrement sur des surfaces de vraisemblance plates pour des rayons supérieurs à 11 km. De plus, la marginalisation bayésienne sur de larges plages de priors peut introduire des effets de volume qui biaisent les estimations de paramètres. Bien que les méthodes bayésiennes soient la norme en astrophysique, l'inférence fréquentiste reste dominante en physique des particules. Ce travail répond à la nécessité de comparer les approches bayésienne et fréquentiste pour garantir la robustesse des contraintes de physique fondamentale dérivées de X-PSI.

Méthodologie
Les auteurs ont implémenté un pipeline d'inférence fréquentiste utilisant la technique de la vraisemblance profilée (profile likelihood) au sein de l'infrastructure existante de X-PSI, garantissant que le modèle physique sous-jacent et l'implémentation de la vraisemblance restent identiques à l'approche bayésienne standard.

• Source de données : L'analyse a utilisé un ensemble de données synthétiques NICER + XMM-Newton (désigné syntX1), précédemment utilisé dans H25 pour tester les échantillonneurs imbriqués.
• Configuration du modèle : L'étude a employé le modèle de profil d'impulsion ST-U, présentant deux points chauds circulaires, simplement connectés et à température unique. L'atmosphère de l'étoile à neutrons a été modélisée avec le modèle d'hydrogène entièrement ionisé NSX, la propagation des photons étant calculée via un lancer de rayons en relativité générale dans un espace-temps de Schwarzschild oblate.
• Implémentation fréquentiste :
  • Le vecteur de paramètres a été divisé en le paramètre d'intérêt (rayon équatorial, $R_{eq}$) et les paramètres de nuisance ($\alpha$).
  • Les paramètres de nuisance (incluant la masse, la distance, l'inclinaison, la géométrie du point chaud, les températures et la colonne de densité d'hydrogène interstellaire) ont été profilés en maximisant la vraisemblance $L(R_{eq}, \alpha)$ par rapport à $\alpha$ pour des valeurs fixes de $R_{eq}$.
  • L'optimisation a été réalisée à l'aide de l'algorithme de simplex Nelder-Mead (scipy.optimize.minimize).
  • Les intervalles de confiance ont été construits à l'aide de la statistique $\Delta\chi^2$ ($\Delta\chi^2 = -2 \ln L_{profile} + 2 \ln L_{max}$), en supposant une distribution $\chi^2$ pour un degré de liberté.
• Base de comparaison : Les résultats ont été comparés à une analyse bayésienne standard utilisant l'algorithme d'échantillonnage imbriqué MultiNest avec des configurations X-PSI identiques à basse résolution.

Principales contributions

1. Application de preuve de concept : Ce travail démontre la première application de l'inférence fréquentiste par vraisemblance profilée pour déterminer les rayons des étoiles à neutrons à l'aide du package X-PSI.
2. Comparaison méthodologique : Il fournit une comparaison directe entre les pipelines d'inférence bayésienne et fréquentiste sur des données synthétiques identiques, isolant la méthodologie d'inférence comme seule variable.
3. Science ouverte : Les auteurs ont rendu publics les modules de modification de leur code d'inférence fréquentiste pour X-PSI afin de faciliter la recherche ultérieure.

Résultats
En utilisant l'ensemble de données synthétiques syntX1 avec un rayon injecté de 12,176 km :

• Inférence bayésienne : L'intervalle de crédibilité à 68 % marginalisé pour le rayon était $R_{68,bayes} = 11,970^{+0,180}_{-0,181}$ km. Ce résultat concordait avec la valeur injectée à $1,1\sigma$.
• Inférence fréquentiste : La vraisemblance profilée a produit un rayon de meilleur ajustement de 12,096 km. L'intervalle de confiance à 68 % était $r_{68} = 12,096^{+0,064}_{-0,182}$ km, et l'intervalle à 95 % était $r_{95} = 12,096^{+0,316}_{-0,397}$ km. L'estimation fréquentiste concordait avec la valeur injectée à $0,4\sigma$.
• Précision : La précision de la méthode fréquentiste s'est avérée comparable à celle de l'analyse bayésienne.
• Efficacité computationnelle : L'analyse fréquentiste a nécessité environ 9 heures-cœur CPU, tandis que l'analyse bayésienne a nécessité environ 3 840 heures-cœur CPU (80 heures sur un nœud de 48 cœurs). L'approche fréquentiste est donc de plusieurs ordres de grandeur plus rapide.

Signification et affirmations
Le document positionne ce travail comme une étude de preuve de concept. Les auteurs affirment que l'inférence fréquentiste via la vraisemblance profilée peut récupérer avec succès le vrai rayon de l'étoile à neutrons avec une précision $< 1\sigma$ en utilisant des données synthétiques. La technique offre une précision comparable à l'analyse bayésienne, mais avec un coût computationnel nettement réduit. Les auteurs concluent que cette méthode constitue un complément précieux à l'inférence bayésienne, offrant une alternative robuste pour l'estimation de paramètres dans X-PSI. Ils déclarent explicitement que les travaux futurs consisteront à appliquer cette technique à des données de pulsars réels provenant du télescope NICER, mais aucune application de ce type n'est présentée dans ce manuscrit.