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

Sintesi Tecnica: Determinazione del Raggio della Stella di Neutroni tramite Modellazione del Profilo di Pulsazione con Verosimiglianza di Profilo

Definizione del Problema
Il Neutron Star Interior Composition Explorer (NICER) utilizza la modellazione del profilo di pulsazione dei raggi X per determinare la massa e il raggio delle stelle di neutroni, impiegando principalmente il pacchetto X-ray Pulse Simulation and Inference (X-PSI). La pratica standard attuale si basa sull'inferenza bayesiana utilizzando algoritmi di nested sampling (ad esempio, MultiNest) per stimare le distribuzioni a posteriori. Sono state osservate discrepanze nei raggi inferiti in diversi studi, potenzialmente attribuibili alla scelta del campionatore o a problemi di convergenza, in particolare su superfici di verosimiglianza piatte per raggi superiori a 11 km. Inoltre, la marginalizzazione bayesiana su ampi intervalli di prior può introdurre effetti di volume che introducono bias nelle stime dei parametri. Sebbene i metodi bayesiani siano lo standard astrofisico, l'inferenza frequentista rimane dominante nella fisica delle particelle. Questo lavoro affronta la necessità di confrontare gli approcci bayesiano e frequentista per garantire la robustezza dei vincoli di fisica fondamentale derivati da X-PSI.

Metodologia
Gli autori hanno implementato una pipeline di inferenza frequentista utilizzando la tecnica della verosimiglianza di profilo (profile likelihood) all'interno dell'esistente framework X-PSI, garantendo che il modello fisico sottostante e l'implementazione della verosimiglianza rimanessero identici all'approccio bayesiano standard.

• Fonte dei Dati: L'analisi ha utilizzato un dataset sintetico NICER + XMM-Newton (designato come syntX1), precedentemente utilizzato in H25 per testare i nested sampler.
• Configurazione del Modello: Lo studio ha impiegato il modello di profilo di pulsazione ST-U, caratterizzato da due punti caldi circolari a temperatura singola e semplicemente connessi. L'atmosfera della stella di neutroni è stata modellata utilizzando il modello di idrogeno completamente ionizzato NSX, con la propagazione dei fotoni calcolata tramite ray tracing in relatività generale in uno spaziotempo di Schwarzschild oblato.
• Implementazione Frequentista:
  • Il vettore dei parametri è stato diviso nel parametro di interesse (raggio equatoriale, $R_{eq}$) e nei parametri di disturbo ($\alpha$).
  • I parametri di disturbo (inclusi massa, distanza, inclinazione, geometria del punto caldo, temperature e densità di colonna dell'idrogeno interstellare) sono stati profilati massimizzando la verosimiglianza $L(R_{eq}, \alpha)$ rispetto ad $\alpha$ per valori fissi di $R_{eq}$.
  • L'ottimizzazione è stata eseguita utilizzando l'algoritmo simplex di Nelder-Mead (scipy.optimize.minimize).
  • Gli intervalli di confidenza sono stati costruiti utilizzando la statistica $\Delta\chi^2$ ($\Delta\chi^2 = -2 \ln L_{profile} + 2 \ln L_{max}$), assumendo una distribuzione $\chi^2$ per un grado di libertà.
• Base di Confronto: I risultati sono stati confrontati con un'analisi bayesiana standard utilizzando l'algoritmo di nested sampling MultiNest con impostazioni identiche e configurazioni X-PSI a bassa risoluzione.

Contributi Chiave

1. Applicazione Proof-of-Principle: Questo lavoro dimostra la prima applicazione dell'inferenza frequentista tramite verosimiglianza di profilo per determinare i raggi delle stelle di neutroni utilizzando il pacchetto X-PSI.
2. Confronto Metodologico: Fornisce un confronto diretto tra le pipeline di inferenza bayesiana e frequentista su dati sintetici identici, isolando la metodologia di inferenza come l'unica variabile.
3. Open Science: Gli autori hanno rilasciato pubblicamente le modifiche al codice di inferenza frequentista per X-PSI per facilitare ulteriori ricerche.

Risultati
Utilizzando il dataset sintetico syntX1 con un raggio iniettato vero di 12.176 km:

• Inferenza Bayesiana: L'intervallo di credibilità al 68% marginalizzato per il raggio è stato $R_{68,bayes} = 11.970^{+0.180}_{-0.181}$ km. Questo risultato concordava con il valore iniettato entro $1.1\sigma$.
• Inferenza Frequentista: La verosimiglianza di profilo ha prodotto un raggio di miglior adattamento (best-fit) di 12.096 km. L'intervallo di confidenza al 68% è stato $r_{68} = 12.096^{+0.064}_{-0.182}$ km, e l'intervallo al 95% è stato $r_{95} = 12.096^{+0.316}_{-0.397}$ km. La stima frequentista concordava con il valore iniettato entro $0.4\sigma$.
• Precisione: La precisione del metodo frequentista è risultata comparabile all'analisi bayesiana.
• Efficienza Computazionale: L'analisi frequentista ha richiesto circa 9 ore-core di CPU, mentre l'analisi bayesiana ha richiesto circa 3.840 ore-core di CPU (80 ore su un nodo a 48 core). L'approccio frequentista è dunque di ordini di grandezza più veloce.

Significatività e Rivendicazioni
L'articolo posiziona questo lavoro come uno studio di prova di principio (proof-of-principle). Gli autori affermano che l'inferenza frequentista tramite verosimiglianza di profilo può recuperare con successo il vero raggio della stella di neutroni con un'accuratezza $< 1\sigma$ utilizzando dati sintetici. La tecnica offre una precisione comparabile all'analisi bayesiana, ma con un costo computazionale significativamente ridotto. Gli autori concludono che questo metodo funge da prezioso complemento all'inferenza bayesiana, offrendo un'alternativa robusta per la stima dei parametri in X-PSI. Essi dichiarano esplicitamente che il lavoro futuro comporterà l'applicazione di questa tecnica ai dati reali delle pulsar del telescopio NICER, ma nessuna tale applicazione è presentata in questo manoscritto.