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

Technische Samenvatting: Bepaling van de Neutronensterstraal via Pulsatieprofielmodellering met behulp van Profile Likelihood

Probleemstelling
De Neutron Star Interior Composition Explorer (NICER) maakt gebruik van röntgenpulsatieprofielmodellering om de massa en straal van neutronensterren te bepalen, waarbij primair het X-ray Pulse Simulation and Inference (X-PSI) pakket wordt ingezet. De huidige standaardpraktijk vertrouwt op Bayesiaanse inferentie met behulp van nested sampling-algoritmen (bijv. MultiNest) om posterieure verdelingen te schatten. Er zijn echter discrepanties waargenomen in de afgeleide radii tussen verschillende studies, wat mogelijk toe te schrijven is aan de keuze van de sampler of convergentieproblemen, met name op vlakke likelihood-oppervlakken voor radii groter dan 11 km. Hoewel Bayesiaanse methoden de astrofysische standaard zijn, blijft frequentistische inferentie dominant in de deeltjeskunde. Dit werk adresseert de noodzaak om Bayesiaanse en frequentistische benaderingen te vergelijken om de robuustheid van fundamentele natuurkundige beperkingen afgeleid uit X-PSI te waarborgen.

Methodologie
De auteurs hebben een frequentistische inferentie-pipeline geïmplementeerd met behulp van de profile likelihood techniek binnen het bestaande X-PSI-framework, waarbij ervoor is gezorgd dat het onderliggende fysieke model en de likelihood-implementatie identiek bleven aan de standaard Bayesiaanse benadering.

• Databron: De analyse maakte gebruik van een synthetische NICER + XMM-Newton dataset (aangeduid als syntX1), die eerder in H25 werd gebruikt om nested samplers te testen.
• Modelconfiguratie: De studie maakte gebruik van het ST-U pulsatieprofielmodel, met twee enkelvoudige temperatuur, simpel verbonden circulaire hotspots. De atmosfeer van de neutronenster werd gemodelleerd met het NSX volledig geïoniseerde waterstofmodel, waarbij de fotonpropagatie werd berekend via algemene-relativistische straaltracking in een afgeplatte Schwarzschild-ruimtetijd.
• Frequentistische Implementatie:
  • De parametervector werd verdeeld in de parameter van belang (equatoriale straal, $R_{eq}$) en de nuisance-parameters ($\alpha$).
  • De nuisance-parameters (inclusief massa, afstand, inclinatie, hotspot-geometrie, temperaturen en interstellaire waterstofkolomdichtheid) werden geprofileerd door de likelihood $L(R_{eq}, \alpha)$ te maximaliseren met betrekking tot $\alpha$ voor vaste waarden van $R_{eq}$.
  • Optimalisatie werd uitgevoerd met behulp van het Nelder-Mead simplex algoritme (scipy.optimize.minimize).
  • Betrouwbaarheidsintervallen werden geconstrueerd met behulp van de $\Delta\chi^2$-statistiek ($\Delta\chi^2 = -2 \ln L_{profile} + 2 \ln L_{max}$), uitgaande van een $\chi^2$-verdeling voor één vrijheidsgraad.
• Vergelijkingsbasis: De resultaten werden vergeleken met een standaard Bayesiaanse analyse met behulp van het MultiNest nested sampling algoritme met identieke instellingen en lage-resolutie X-PSI configuraties.

Belangrijkste Bijdragen

1. Proof-of-Principle Toepassing: Dit werk vormt de eerste toepassing van frequentistische profile likelihood inferentie voor het bepalen van neutronensterradii met behulp van het X-PSI pakket.
2. Methodologische Vergelijking: Het biedt een directe vergelijking tussen Bayesiaanse en frequentistische inferentie-pipelines op identieke synthetische data, waarbij de inferentiemethodologie de enige variabele is.
3. Open Science: De auteurs hebben hun frequentistische inferentie-code-modificaties voor X-PSI publiekelijk vrijgegeven om verder onderzoek te faciliteren.

Resultaten
Gebruikmakend van de synthetische dataset syntX1 met een geïnjecteerde ware straal van 12,176 km:

• Bayesiaanse Inferentie: Het gemarginaliseerde 68% geloofwaardigheidsinterval voor de straal was $R_{68,bayes} = 11,970^{+0,180}_{-0,181}$ km. Dit resultaat kwam binnen $1,1\sigma$ overeen met de geïnjecteerde waarde.
• Frequentistische Inferentie: De profile likelihood leverde een best-fit straal van 12,096 km. Het 68% betrouwbaarheidsinterval was $r_{68} = 12,096^{+0,064}_{-0,182}$ km, en het 95% interval was $r_{95} = 12,096^{+0,316}_{-0,397}$ km. De frequentistische schatting kwam binnen $0,4\sigma$ overeen met de geïnjecteerde waarde.
• Precisie: De precisie van de frequentistische methode bleek vergelijkbaar met de Bayesiaanse analyse.
• Computationele Efficiëntie: De frequentistische analyse vereiste ongeveer 9 CPU-core uren, terwijl de Bayesiaanse analyse ongeveer 3.840 CPU-core uren vereiste (80 uur op een 48-core node). De frequentistische benadering is hiermee ordes van grootte sneller.

Betekenis en Claims
Het artikel positioneert dit werk als een proof-of-principle studie. De auteurs claimen dat frequentistische inferentie via profile likelihood succesvol de ware neutronensterstraal kan herstellen met een nauwkeurigheid van $< 1\sigma$ met behulp van synthetische data. De techniek biedt een precisie die vergelijkbaar is met de Bayesiaanse analyse, maar met aanzienlijk lagere computationele kosten. De auteurs concluderen dat deze methode dient als een waardevolle aanvulling op Bayesiaanse inferentie en een robuust alternatief biedt voor parameterestimatie in X-PSI. Zij stellen expliciet dat toekomstig werk zal bestaan uit het toepassen van deze techniek op echte pulsar-data van de NICER-telescoop, maar geen dergelijke toepassing wordt gepresenteerd in dit manuscript.