Investigating Twin Star Equation of States in Light of Recent Astrophysical Observations

Resumen Técnico: Investigación de las Ecuaciones de Estado de Estrellas Gemelas a la Luz de Recientes Observaciones Astrofísicas

Planteamiento del Problema
La composición interna de las estrellas de neutrones (EN) y la naturaleza de la transición de fase (TF) de la materia hadrónica (MH) a la materia de quarks (MQ) siguen sin resolverse. Aunque la Cromodinámica Cuántica (QCD) predice una transición a altas densidades, las condiciones específicas —como la densidad de inicio, la presión de transición y la fuerza de la discontinuidad— son desconocidas. Una transición de fase de primer orden fuerte podría dar lugar a la existencia de "estrellas gemelas" (ET): un tercer brazo estable de objetos compactos en el diagrama masa-radio ($M-R$), separado del brazo estándar de las estrellas de neutrones por una región inestable. Estas estrellas gemelas poseerían una masa igual pero radios menores que sus contrapartes de neutrones. El desafío radica en identificar las ecuaciones de estado (EoS) híbridas que satisfagan tanto el equilibrio hidrostático de la relatividad general como las restricciones cada vez más precisas de la astrofísica multi-mensajero, incluyendo las mediciones de masa-radio de NICER y las observaciones de ondas gravitacionales (GW170817).

Metodología
Los autores adoptan un enfoque agnóstico al modelo para construir un conjunto exhaustivo de EoS híbridas capaces de producir soluciones de estrellas gemelas. La metodología se basa en las ecuaciones de Tolman-Oppenheimer-Volkoff (TOV) para resolver el equilibrio hidrostático. La construcción de la EoS involucra cuatro regiones distintas:

1. Región de la Corteza: Una EoS politrópica segmentada fija por debajo de $\sim 0.4\rho_{sat}$.
2. Región Hadrónica: Modelada mediante una politrope con un índice adiabático $\gamma$ libre, determinado por las condiciones de contorno en la transición.
3. Región de la Transición de Fase: Modelada mediante una construcción de Maxwell, caracterizada por una presión de transición constante ($P_{tr}$) y un salto discontinuo en la densidad de energía ($\Delta e$).
4. Región de Quarks: Modelada mediante una parametrización de velocidad del sonido constante (CSS, por sus siglas en inglés) ($c^2_{s,QM}$).

El espacio de parámetros se define por cuatro variables independientes: la densidad de energía de transición ($e_{tr}$), la presión de transición ($P_{tr}$), la fuerza de la discontinuidad de la densidad de energía ($\Delta e$) y la velocidad del sonido en la materia de quarks ($c^2_{s,QM}$). El estudio varía sistemáticamente estos parámetros dentro de los límites de causalidad para identificar las condiciones que produzcan una tercera rama estable en la secuencia $M-R$.

Contribuciones Clave

• Definición de Curvas de "Sombrero de Bruja": El artículo introduce una representación geométrica en el espacio de parámetros $(P_{tr}, \Delta e)$ para valores fijos de $e_{tr}$ y $c^2_{s,QM}$. Estas curvas delinean los límites de las fuerzas de discontinuidad permitidas que permiten la formación de estrellas gemelas. Los autores identifican tres regiones distintas en este espacio:
  • Región sin Solución: No se forman estrellas gemelas independientemente de $\Delta e$.
  • Región Limitada: Las estrellas gemelas se forman para un rango específico de $\Delta e$ por debajo de una presión de transición máxima.
  • Región No Limitada: Las estrellas gemelas se forman para un $\Delta e$ por encima de un límite mínimo, sin un límite superior para $\Delta e$ (dentro de la configuración numérica).
• Análisis de Causalidad y Completitud: El estudio revela que las curvas de "sombrero de bruja" pueden ser "perforadas" o incompletas. Esto ocurre cuando la presión de transición máxima permitida ($P_{tr}^{max}$) requerida para formar una estrella gemela excede el límite de causalidad ($c_s = 1$) o la presión central de la configuración de masa máxima hadrónica.
• Integración de Restricciones Observacionales Recientes: Los autores aplican intervalos de confianza de 2$\sigma$ (y 1$\sigma$/3$\sigma$) de mediciones recientes de NICER (PSR J0614–3329, PSR J0437–4715, PSR J0740+6620, PSR J0030+0451), el púlsar masivo PSR J0348+0432, y la deformabilidad de marea de GW170817 para restringir el espacio de parámetros.

Resultados
La aplicación de datos observacionales restringe significativamente el espacio de parámetros viable para las estrellas gemelas:

• Restricciones sobre la Densidad de Transición: Mediciones recientes de PSR J0614–3329 y PSR J0437–4715 sugieren una EoS más blanda alrededor de la masa canónica, lo que restringe fuertemente la densidad de energía de transición. El estudio establece un límite superior en la densidad de energía de transición de $e_{tr}^{max} \approx 1.40 \times 10^{15}$ g/cm$^3$ ($\approx 4.03\rho_{sat}$) para la EoS de quarks más rígida ($c^2_{s,QM}=1.0$).
• Restricciones sobre la Fuerza de la Discontinuidad: El máximo de la discontinuidad permitida en la densidad de masa en reposo se encuentra en $\Delta\rho_{max} \approx 7.76\rho_{sat}$ ($\Delta e_{max} \approx 2.94 \times 10^{15}$ g/cm$^3$). Este límite está gobernado por las EoS que producen estrellas gemelas de Categoría II (donde el segundo brazo soporta $>2.0 M_\odot$ pero el tercero no).
• Masa Máxima de las Estrellas Gemelas: La masa máxima alcanzable por las estrellas gemelas se calcula como $M_{TS}^{max} \approx 2.05 M_\odot$ (para $c^2_{s,QM}=1.0$ y 2$\sigma$ de confianza). Este límite es establecido por las estrellas gemelas de Categoría I (donde ambos brazos soportan $>2.0 M_\odot$).
• Impacto de los Niveles de Confianza: Los resultados son sensibles al nivel de confianza de las inferencias observacionales. En un nivel de confianza estricto de 1$\sigma$, las estrellas gemelas de Categoría I y II quedan totalmente descartadas, dejando solo soluciones de Categoría III y IV con masas máximas significativamente menores ($\approx 1.37 M_\odot$). Por el contrario, en 3$\sigma$, las restricciones se relajan, permitiendo un rango más amplio de parámetros, incluyendo casos donde la presión de transición se encuentra en la región "no limitada".
• Dependencia de la Categoría: El estudio categoriza las estrellas gemelas basándose en las masas máximas de sus segundos ($M_{TOV,2}$) y terceros ($M_{TOV,3}$) brazos. Encuentra que las estrellas gemelas de Categoría I (masivas en ambos brazos) solo están permitidas en densidades de transición extremadamente altas ($\approx 4\rho_{sat}$), mientras que las densidades de transición más blandas favorecen las Categorías II, III y IV.

Significancia
El artículo afirma que la observación de estrellas gemelas, particularmente las masivas, serviría como una firma de "pistola humeante" (smoking gun) de una fuerte transición de fase de primer orden de hadrón a quark en los núcleos de las estrellas de neutrones. Al mapear las curvas de "sombrero de bruja" y aplicar rigurosas restricciones observacionales, el estudio proporciona límites termodinámicos superiores específicos para la densidad de transición ($e_{tr}$), la fuerza de la discontinuidad ($\Delta e$) y la masa máxima de las estrellas gemelas. Estos límites ofrecen un marco crítico para interpretar futuros datos astrofísicos y reducir la posible naturaleza de la materia densa. Los autores señalan que, aunque su análisis se basa en el equilibrio hidrostático, la existencia de las estrellas gemelas de Categoría II es particularmente frágil; estudios dinámicos previos sugieren que pueden ser inestables ante perturbaciones, migrando a los brazos hadrónicos o colapsando en agujeros negros. Por lo tanto, la identificación de EoS de estrellas gemelas viables es un paso necesario, aunque potencialmente insuficiente, hacia la confirmación de la existencia de materia de quarks en las estrellas de neutrones.