Are NICER and GW170817 constraints suggesting a compactified scenario for Neutron stars?

Resumen Técnico: ¿Sugieren las restricciones de NICER y GW170817 un escenario de compactación para las estrellas de neutrones?

Planteamiento del problema
La ecuación de estado (EoS, por sus siglas en inglés) de la materia a densidades que exceden la saturación nuclear sigue siendo un problema fundamental no resuelto en la física. Mientras que la teoría de campo efectivo quiral ($\chi$EFT) modela de manera fiable las densidades subnucleares y la QCD perturbativa describe los regímenes de alta densidad asintótica, la región de densidad intermedia dentro de los núcleos de las estrellas de neutrones (NS) es poco comprendida. Observaciones astrofísicas recientes del Neutron Star Interior Composition Explorer (NICER) y el evento de ondas gravitacionales GW170817 han proporcionado restricciones precisas sobre las masas, radios y deformabilidades de marea de las NS. Sin embargo, sigue sin estar claro si estas restricciones favorecen una clase específica de EoS (monotónica, no monotónica o discontinua) o si revelan un patrón estructural universal en la materia densa. Específicamente, este trabajo investiga si los datos actuales sugieren un escenario de "compactación" donde la gravedad favorece naturalmente las configuraciones que maximizan la compacidad estelar por encima de aquellas que simplemente maximizan la masa.

Metodología
Los autores emplean un marco de inferencia bayesiana integral para analizar tres conjuntos de EoS distintos y ajenos a clases específicas:

1. Monotónica: Caracterizada por una velocidad del sonido ($c_s^2$) que aumenta monotonicamente con la densidad (típica de las EoS hadrónicas).
2. No monotónica: Caracterizada por un pico en $c_s^2$ seguido de un descenso (asociada con transiciones de fase suaves).
3. Discontinua: Caracterizada por un salto de densidad ($\Delta n$) en una transición de fase (asociada con transiciones de fase de primer orden fuertes).

Construcción de la EoS:
La EoS se construye mediante un enfoque híbrido. La materia de baja densidad ($< 0.5 n_0$) se define mediante la EoS tabulada BPS. La región de $0.5 n_0$ a $1.1 n_0$ se modela mediante una politropa con un índice adiabático $\Gamma$ restringido por las bandas de incertidumbre de $\chi$EFT. Para densidades superiores a $1.1 n_0$, la velocidad del sonido al cuadrado $c_s^2$ se parametriza como una función lineal por tramos del potencial químico $\mu$ a través de cuatro segmentos aleatorizados, extendiéndose hasta el régimen de pQCD ($\mu \approx 2.6$ GeV).

Marco de inferencia:
A diferencia de los métodos previos de muestreo por rechazo que generan millones de muestras de EoS y las filtran, este estudio utiliza el algoritmo de muestreo anidado UltraNest para explorar directamente la distribución posterior. El modelo contiene 10 parámetros libres (potenciales químicos y valores de la velocidad del sonido al cuadrado para los segmentos, además de $\Gamma_{CET}$ y un parámetro de salto de densidad $\delta n$ para la clase discontinua).
El análisis incorpora:

• Datos de NICER: Distribuciones posteriores conjuntas para masa y radio de PSR J0030+0451, PSR J0740+6620, PSR J0437−4715 y el recientemente reportado PSR J0614−3329.
• Datos de GW170817: Muestras posteriores de masas de componentes y deformabilidades de marea de la fusión de estrellas de neutrones binarias.
• Restricciones de $\chi$EFT: Verosimilitudes gaussianas para la presión a bajas densidades.

El estudio compara las distribuciones posteriores de los parámetros de la EoS y las propiedades estelares derivadas (secuencias Masa-Radio) bajo diferentes combinaciones de restricciones observacionales.

Resultados clave

1. Capacidad de restricción: El análisis bayesiano confirma que los datos observacionales restringen eficazmente la región de baja densidad de la EoS (específicamente $\Gamma_{CET}$ y el primer segmento de la velocidad del sonido $c_{s,1}^2$), la cual permanece robusta en todas las clases de EoS. En contraste, los parámetros de alta densidad ($c_{s,2}^2, c_{s,3}^2$) exhiben distribuciones posteriores más amplias, lo que indica una incertidumbre significativa restante.
2. Impacto de los nuevos datos: La inclusión de PSR J0614−3329 y GW170817 desplaza sistemáticamente los parámetros del potencial químico hacia valores más altos y reduce los límites de la masa máxima. Sin embargo, la evidencia bayesiana ($\ln Z$) no favorece fuertemente una clase de EoS sobre otra; las diferencias en la evidencia caen en el rango de "inconclusivo a débil".
3. Masa máxima vs. Compacidad máxima: Se establece una distincción crítica entre la EoS que produce la masa posible máxima ($M_{TOV}$) y la Eosa que produce la compacidad máxima ($C = GM/Rc^2$).
   • La EoS de Masa Máxima es rígida desde densidades muy bajas, proporcionando un soporte de presión uniforme para sostener las masas más grandes. Estas configuraciones ocupan el borde de radio grande de la posterior M-R.
   • La EoS de Compacidad Máxima presenta una EoS de baja a intermedia densidad más blanda (minimizando el radio) seguida de un endurecimiento pronunciado a altas densidades (sosteniendo la masa).
4. Preferencia observacional: La EoS más favorecida por los datos combinados de NICER y GW170817 (mayor $\ln Z$) se asemeja estructuralmente a la EoS de Compacidad Máxima, no a la de Masa Máxima. La EoS favorecida por la observación exhibe una relativa blandura a densidades intermedias, transicionando hacia la rigidez a altas densidades.
5. Transiciones de fase: Para la clase discontinua, los datos prefieren discontinuidades de densidad pequeñas ($\delta n \approx 0.09-0.10$ fm$^{-3}$) y transiciones de fase tardías (que ocurren a densidades más altas). Aunque una transición de fase de baja densidad podría teóricamente producir una estrella altamente compacta con un núcleo de quarks grande, las predicciones bayesianas favorecen configuraciones donde la transición de fase ocurre más tarde, confinando la materia de quarks a los núcleos de estrellas suficientemente masivas.

Significación y afirmaciones
El artículo afirma que las restricciones astrofísicas actuales están seleccionando configuraciones de EoS que producen estrellas compactas en lugar de las máximamente masivas. La convergencia de la EoS favorecida por la observación con la secuencia de máxima compacidad sugiere una propiedad fundamental de la materia densa: la EoS probablemente experimenta un endurecimiento significativo a densidades intermedias (señalando potencialmente el inicio de nuevos grados de libertad) antes de volver a endurecerse para sostener las estrellas de neutrones observadas de $\sim 2 M_\odot$.

Los autores concluyen que esta estructura de "blando-a-rígido" es lo que la gravedad favorece naturalmente para lograr la compactación. Si bien los datos no descartan definitivamente clases de EoS específicas, apuntan hacia un patrón universal donde la ecuación de estado transiciona de la blandura a la rigidez, impulsando a las estrellas de neutrones hacia sus límites de compacidad. El artículo señala modestamente que será necesario contar con futuras restricciones de púlsares de baja masa (para definir mejor la EoS de baja densidad) y de púlsares de alta masa (para restringir la Eos de alta densidad) para refinar aún más estas conclusiones y potencialmente descartar escenarios específicos, como las transiciones de fase tempranas.