Is the coexistence of strange quark stars and hadronic stars favored by astrophysical data? A Bayesian analysis

Resumen Técnico: Coexistencia de Estrellas de Quarks Extrañas y Estrellas Hadónicas

Planteamiento del Problema
La composición microscópica de las estrellas compactas y la ecuación de estado (EOS) de la materia fuertemente interactuante densa siguen siendo problemas abiertos centrales en la física nuclear y la astrofísica. Existe una tensión específica entre las restricciones observacionales que requieren una EOS "suave" a densidades intermedias (para explicar los radios pequeños de objetos de masa ordinaria y la deformabilidad de marea de GW170817) y la necesidad de una EOS "rígida" a altas densidades (para soportar púlsares con masas que exceden los $2 M_\odot$). Además, la inclusión de hiperones y resonancias $\Delta$ en la materia hadónica típicamente suaviza la EOS, lo que conduce al "problema del hiperón", donde las masas máximas caen por debajo de los valores observados.

El escenario de las "dos familias" (2F) propone una resolución: la coexistencia de dos ramas desconectadas de estrellas compactas. En este marco, las estrellas hadónicas (HSs) son metaestables y pueblan una rama con radios más pequeños, mientras que la materia de quarks extraños (SQM) absolutamente estable forma una rama auto-ligada de estrellas de quarks (QSs) capaz de soportar masas más altas. Si bien este escenario ha sido explorado teóricamente, no había sido sometido previamente a un análisis bayesiano riguroso que lo comparara directamente contra la hipótesis estándar de "una familia" (1F) (donde todas las estrellas son hadónicas) utilizando un conjunto unificado de datos astrofísicos y de laboratorio.

Metodología
Los autores realizan el primer análisis bayesiano detallado para restringir los parámetros de la EOS tanto para los escenarios 1F como 2F simultáneamente.

• Modelos:
  • Sector Hadrónico: Un modelo de campo medio relativista no lineal (RMF) que incluye nucleones, hiperones ($\Lambda, \Sigma, \Xi$) y resonancias $\Delta$. El modelo utiliza parámetros nucleares empíricos en saturación y acoplamientos dependientes de la densidad para el mesón $\rho$ para controlar la pendiente de la energía de simetría ($L_{sym}$).
  • Sector de Quarks (solo 2F): Una EOS de modelo de bolsa de color-sabor bloqueado (CFL) que describe la SQM absolutamente estable. Este sector incluye parámetros para la constante de bolsa ($B$), la brecha de apareamiento CFL ($\Delta_{CFL}$) y correcciones de QCD perturbativa ($a_4$).
• Datos y Verosimilitudes: El análisis incorpora:
  • Restricciones de Masa-Radio (M-R): Posteriors conjuntos de NICER (PSR J0030+0451, PSR J0740+6620, PSR J0614–3329) y observaciones de rayos X de HESS J1731–347.
  • Restricciones de Solo Masa: PSR J0952–0607 y el componente secundario de GW190814 (tratado como un potencial objeto compacto).
  • Ondas Gravitacionales: La posterior conjunta completa de GW170817 (masas y deformabilidades de marea), tratada mediante la marginalización sobre las etiquetas de familia (HS-HS, HS-QS, QS-HS, QS-QS) en el escenario 2F.
  • Restricciones de Laboratorio/Teóricas: Bandas de la teoría de campo efectivo quiral ($\chi$EFT) para materia de neutrones pura y datos de flujo de colisiones de iones pesados (HIC) para materia nuclear simétrica.
• Marco Estadístico:
  • El estudio utiliza muestreo anidado dinámico (vía el paquete dynesty) para computar las distribuciones posteriores y las evidencias bayesianas ($Z$).
  • El factor de Bayes ($B_{2F,1F} = Z_{2F}/Z_{1F}$) se utiliza para comparar los modelos, teniendo en cuenta las diferentes dimensionalidades del espacio de parámetros (11 parámetros para 1F frente a 14 para 2F).
  • En el escenario 2F, la asignación de familia de cada fuente observada se trata como una variable latente y se marginaliza, con priors definidos basados en la estabilidad estructural y los canales de formación astrofísica.

Contribuciones Clave

1. Primera Comparación Bayesiana: Este trabajo proporciona la primera comparación bayesiana rigurosa entre los escenarios 1F y 2F utilizando un conjunto de datos exhaustivo que incluye tanto restricciones nucleares de baja densidad como observaciones astrofísicas de alta masa.
2. Marco Unificado: Restringe la EOS hadónica metaestable y la EOS de la SQM absolutamente estable simultáneamente dentro de la hipótesis 2F, en lugar de asumir la naturaleza de las fuentes observadas a priori.
3. Resolución de la Tensión: El análisis demuestra cómo el escenario 2F resuelve naturalmente la tensión entre la suavidad de densidad intermedia y la rigidez de alta densidad al separar estos requisitos en dos ramas estelares distintas.

Resultados

• Preferencia del Modelo: El escenario 2F es favorecido sobre el escenario 1F en todas las combinaciones de conjuntos de datos. Para el conjunto de datos completo (incluyendo HIC y todas las restricciones astrofísicas), el log-factor de Bayes es $\ln B_{2F,1F} = 11.9$, lo que indica un apoyo "decisivo" para el escenario de las dos familias.
• Características de la EOS:
  • Escenario 1F: Para satisfacer tanto las restricciones de radio pequeño como el límite de $2 M_\odot$, la EOS 1F requiere un compromiso finamente equilibrado. Esto se logra mediante un inicio temprano de las resonancias $\Delta$ ($\sim 1.8 n_{sat}$) que suaviza la EOS a densidades intermedias, seguido de un fuerte endurecimiento a altas densidades. La masa máxima se restringe a $M_{max} \approx 2.05 M_\odot$.
  • Escenario 2F: La tensión se alivia. La rama hadónica permanece suave (soportada por las resonancias $\Delta$ e hiperones) con una masa máxima de $M_{HS,max} \approx 1.54 M_\odot$ y un radio de $R_{1.4} \approx 11.1$ km. Los púlsares pesados son soportados por la rama QS más rígida, alcanzando $M_{QS,max} \approx 2.33 M_\odot$.
• Restricciones de Parámetros: El posterior 2F prefiere una masa efectiva de Dirac ($m^*/m$) mayor y un ajuste menos extremo de los acoplamientos $\Delta$ en comparación con el caso 1F. Los parámetros del sector de quarks ($B, \Delta_{CFL}, a_4$) están moderadamente restringidos, con una correlación positiva entre la constante de bolsa y la brecha de apareamiento.
• Interpretación de GW170817: En el escenario 2F, el evento GW170817 es más probablemente una fusión HS-HS, pero los canales mixtos (HS-QS) contribuyen significativamente (~35% de la verosimilitud), mientras que el canal QS-QS es desfavorecido para la mejor EOS.

Significancia y Reivindicaciones
El artículo afirma que los datos astrofísicos y de laboratorio actuales están organizados más naturalmente por dos familias de estrellas compactas desconectadas que por una única secuencia puramente hadónica. El escenario 2F ofrece una explicación natural para la coexistencia de objetos ordinarios muy compactos (como HESS J1731–347) y púlsares muy masivos, eliminando efectivamente el problema del hiperón.

Los autores enfatizan que su conclusión es robusta frente a variaciones en los priors (potenciales de hiperones, rangos de parámetros de SQM) y supuestos de masa de la fuente, siempre que se utilice el conjunto completo de datos. Sin embargo, señalan que la preferencia estadística es menos definitiva si se consideran conjuntos de datos restringidos (excluyendo HESS J1731–347, PSR J0614–3329 y datos de HIC).

El estudio excluye explícitamente las estrellas híbridas (estrellas con núcleos de quarks) de la comparación actual, señalando que una selección de modelos completa debería incluirlas como una tercera hipótesis en trabajos futuros. Los autores concluyen que detecciones futuras de objetos compactos más masivos, objetos de masa ordinaria muy compactos, o mediciones precisas de dos masas distintas con el mismo radio, proporcionarían indicaciones aún más fuertes a favor del escenario de las dos familias.