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

Resumen Técnico: Determinación del Radio de una Estrella de Neutrones mediante el Modelado de Perfil de Pulso usando Verosimilitud de Perfil

Planteamiento del Problema
El Explorador de Composición del Interior de Estrellas de Neutrones (NICER, por sus siglas en inglés) utiliza el modelado de perfiles de pulso de rayos X para determinar la masa y el radio de estrellas de neutrones, empleando principalmente el paquete de Inferencia y Simulación de Pulsos de Rayos X (X-PSI). La práctica estándar actual se basa en la inferencia bayesiana utilizando algoritmos de muestreo anidado (por ejemplo, MultiNest) para estimar las distribuciones a posteriori. Sin embargo, se han observado discrepancias en los radios inferidos en diferentes estudios, lo que podría atribuirse a la elección del muestreador o a problemas de convergencia, particularmente en superficies de verosimilitud planas para radios superiores a 11 km. Además, la marginalización bayesiana sobre rangos amplios de la distribución a priori puede introducir efectos de volumen que sesgan las estimaciones de los parámetros. Si bien los métodos bayesianos son el estándar astrofísico, la inferencia frecuentista sigue siendo dominante en la física de partículas. Este trabajo aborda la necesidad de comparar los enfoques bayesiano y frecuentista para garantizar la robustez de las restricciones de la física fundamental derivadas de X-PSI.

Metodología
Los autores implementaron un flujo de trabajo de inferencia frecuentista utilizando la técnica de verosimilitud de perfil dentro del marco existente de X-PSI, asegurando que el modelo físico subyacente y la implementación de la verosimilitud permanecieran idénticos al enfoque bayesiano estándar.

• Fuente de Datos: El análisis utilizó un conjunto de datos sintéticos de NICER + XMM-Newton (designado como syntX1), utilizado previamente en H25 para probar muestreadores anidados.
• Configuración del Modelo: El estudio empleó el modelo de perfil de pulso ST-U, que presenta dos puntos calientes circulares, de temperatura única y simplemente conectados. La atmósfera de la estrella de neutrones se modeló utilizando el modelo de hidrógeno totalmente ionizado NSX, con la propagación de fotones calculada mediante trazado de rayos en un espacio de Schwarzschild oblato con efectos de relatividad general.
• Implementación Frecuentista:
  • El vector de parámetros se dividió en el parámetro de interés (radio ecuatorial, $R_{eq}$) y los parámetros molestos ($\alpha$).
  • Los parámetros molestos (incluyendo masa, distancia, inclinación, geometría del punto caliente, temperaturas y densidad de columna de hidrógeno interestelar) se perfilaron maximizando la verosimilitud $L(R_{eq}, \alpha)$ con respecto a $\alpha$ para valores fijos de $R_{eq}$.
  • La optimización se realizó utilizando el algoritmo de simplex Nelder-Mead (scipy.optimize.minimize).
  • Los intervalos de confianza se construyeron utilizando el estadístico $\Delta\chi^2$ ($\Delta\chi^2 = -2 \ln L_{profile} + 2 \ln L_{max}$), asumiendo una distribución $\chi^2$ para un grado de libertad.
• Línea Base de Comparación: Los resultados se compararon contra un análisis bayesiano estándar utilizando el algoritmo de muestreo anidado MultiNest con configuraciones de X-PSI de baja resolución e idénticas.

Contribuciones Clave

1. Aplicación de Prueba de Principio: Este trabajo demuestra la primera aplicación de la inferencia de verosimilitud de perfil frecuentista para determinar radios de estrellas de neutrones utilizando el paquete X-PSI.
2. Comparación Metodológica: Proporciona una comparación directa entre flujos de trabajo de inferencia bayesiana y frecuentista sobre datos sintéticos idénticos, aislando la metodología de inferencia como la única variable.
3. Ciencia Abierta: Los autores han publicado sus modificaciones de código de inferencia frecuentista para X-PSI para facilitar investigaciones posteriores.

Resultados
Utilizando el conjunto de datos sintéticos syntX1 con un radio real inyectado de 12.176 km:

• Inferencia Bayesiana: El intervalo de credibilidad del 68% marginalizado para el radio fue $R_{68,bayes} = 11.970^{+0.180}_{-0.181}$ km. Este resultado concordó con el valor inyectado dentro de $1.1\sigma$.
• Inferencia Frecuentista: La verosimilitud de perfil arrojó un radio de mejor ajuste de 12.096 km. El intervalo de confianza del 68% fue $r_{68} = 12.096^{+0.064}_{-0.182}$ km, y el intervalo del 95% fue $r_{95} = 12.096^{+0.316}_{-0.397}$ km. La estimación frecuentista concordó con el valor inyectado dentro de $0.4\sigma$.
• Precisión: Se encontró que la precisión del método frecuentista es comparable a la del análisis bayesiano.
• Eficiencia Computacional: El análisis frecuentista requirió aproximadamente 9 horas de núcleos de CPU, mientras que el análisis bayesiano requirió aproximadamente 3,840 horas de núcleos de CPU (80 horas en un nodo de 48 núcleos). El enfoque frecuentista es, por tanto, órdenes de magnitud más rápido.

Significancia y Reivindicaciones
El artículo posiciona este trabajo como un estudio de prueba de principio. Los autores afirman que la inferencia frecuentista mediante verosimilitud de perfil puede recuperar con éxito el radio real de una estrella de neutrones con una precisión de $< 1\sigma$ utilizando datos sintéticos. La técnica ofrece una precisión comparable al análisis bayesiano pero con un costo computacional significativamente reducido. Los autores concluyen que este método sirve como un complemento valioso para la inferencia bayesiana, ofreciendo una alternativa robusta para la estimación de parámetros en X-PSI. Expresan explícitamente que el trabajo futuro implicará la aplicación de esta técnica a datos reales de púlsares del telescopio NICER, pero tal aplicación no se presenta en este manuscrito.