On the treatment of thermal effects in the equation of state on neutron star merger remnants

Título: Sobre el tratamiento de los efectos térmicos en la ecuación de estado de los remanentes de fusiones de estrellas de neutrones

1. El Problema

La detección de ondas gravitacionales (GW) provenientes de fusiones de estrellas de neutrones binarias (BNS), como GW170817, ha abierto una nueva ventana para estudiar la materia nuclear a densidades supranucleares. Sin embargo, la interpretación de la señal de ondas gravitacionales del post-fusión (especialmente la fase de larga duración del remanente hipermasivo, HMNS) depende críticamente de la Ecuación de Estado (EOS) de la materia nuclear.

El problema central abordado en este trabajo es la incertidumbre sobre cómo el tratamiento de los efectos térmicos en la EOS afecta la dinámica del remanente. La mayoría de las simulaciones numéricas de larga duración utilizan EOS híbridas, que combinan una parte fría (politrópica por tramos) con una presión térmica simplificada (tipo gas ideal). Aunque computacionalmente eficientes, estas aproximaciones pueden no capturar con precisión la física termodinámica compleja (interacciones nucleares, degeneración, composición) presente en las EOS tabuladas de física nuclear. La cuestión clave es si las simplificaciones de las EOS híbridas alteran significativamente la evolución del remanente, la estabilidad contra el colapso, la excitación de modos de oscilación y la señal de ondas gravitacionales resultante, especialmente en las fases tardías (>100 ms).

2. Metodología

Los autores realizaron una suite de simulaciones de relatividad numérica a largo plazo (hasta $\gtrsim 150$ ms) de fusiones de estrellas de neutrones binarias con masas iguales ($M_0 = 1.4 M_\odot$).

• Comparación de EOS: Se compararon dos enfoques para cuatro modelos de EOS micro-físicas (SLy4, DD2, HShen, LS220):
  1. EOS Tabuladas: Utilizan tablas completas de física nuclear que dependen de la densidad bariónica ($\rho$), la temperatura ($T$) y la fracción de electrones ($Y_e$).
  2. EOS Híbridas: Combinan una parte fría (aproximada por politropos por tramos de 7 y 10 piezas) y una parte térmica simplificada ($P_{th} = \kappa_{th} \rho^{\Gamma_{th}}$ con $\Gamma_{th} = 1.8$).
• Códigos Numéricos: Se utilizó el código ILLINOISGRMHD dentro de la infraestructura Einstein Toolkit para evolucionar las ecuaciones de Einstein (formulación BSSN) y la hidrodinámica relativista general (GRHD).
• Análisis de Estabilidad: Se evaluó la estabilidad convectiva del remanente utilizando dos criterios:
  • Criterio de Ledoux (gradientes de entropía y composición).
  • Criterio generalizado de Solberg-Høiland (combinando la frecuencia de Brunt-Väisälä $N^2$ y la frecuencia epicyclic $K^2$ para perturbaciones axisimétricas en relatividad general).
• Análisis de Señales: Se extrajeron las señales de ondas gravitacionales ($h_{22}$), se calcularon sus espectros de frecuencia y se identificaron modos de oscilación (modo $f_2$ fundamental y modos inerciales) mediante el método de Prony (ESPRIT).

3. Contribuciones Clave

• Validación de Modos Inerciales: Confirmación de que la excitación de modos inerciales en las fases tardías del post-fusión, previamente reportada solo en simulaciones con EOS híbridas, también ocurre y es robusta cuando se utilizan EOS tabuladas completas.
• Definición de Temperatura: Un análisis crítico sobre la definición de temperatura en modelos híbridos. Se demuestra que la fórmula estándar utilizada en híbridos ($T_{th} \propto \epsilon_{th}$) subestima sistemáticamente la temperatura real en regímenes degenerados (densidades altas), lo que afecta la presión térmica y la evolución del remanente.
• Nuevas Relaciones Cuasi-Universales: Derivación de nuevas relaciones empíricas para las frecuencias de las ondas gravitacionales post-fusión, mostrando dónde y cuándo las desviaciones entre modelos híbridos y tabulados se vuelven significativas.
• Impacto en la Estabilidad del Remanente: Demostración de que el tratamiento térmico influye directamente en si el remanente colapsa a un agujero negro o permanece estable.

4. Resultados Principales

• Evolución Dinámica y Colapso:

  • Las EOS híbridas tienden a producir remanentes más compactos y propensos al colapso prematuro. Por ejemplo, para la EOS blanda LS220, el modelo híbrido colapsa a un agujero negro a los $\sim 66$ ms, mientras que el modelo tabulado permanece estable durante toda la simulación (>140 ms).
  • Esto se debe a que las EOS tabuladas capturan mejor la presión térmica y los gradientes, proporcionando un soporte contra la gravedad más efectivo que la aproximación de gas ideal.

• Señales de Ondas Gravitacionales (GW):

  • Fase Temprana (<30 ms): Las frecuencias del modo fundamental cuadrupolar ($f_{2,i}$) son similares entre ambos enfoques, aunque existen diferencias en la amplitud y el decaimiento.
  • Fase Tardía (>30-40 ms): Surgen diferencias significativas. Los modelos tabulados muestran una evolución de frecuencia distinta y la aparición de modos inerciales (frecuencias menores a 2 kHz) excitados por convección.
  • Las desviaciones de las relaciones cuasi-universales establecidas (basadas en datos híbridos) son significativas en las fases tardías, lo que implica que usar modelos simplificados podría llevar a inferencias erróneas sobre las propiedades de la estrella de neutrones.

• Inestabilidades Convectivas:

  • Se identificaron patrones convectivos locales y sostenidos en el remanente que persisten más allá de los 100 ms.
  • Aunque la estructura convectiva difiere entre modelos tabulados e híbridos (los híbridos muestran convección más "aparente" debido a la definición de temperatura), ambos generan la excitación de modos inerciales.
  • El uso del criterio de Solberg-Høiland (incluyendo efectos rotacionales) es crucial, ya que el análisis basado solo en la frecuencia de Brunt-Väisälä o en promedios angulares puede enmascarar inestabilidades locales.

• Perfiles de Rotación:

  • Los remanentes con EOS tabuladas muestran gradientes de rotación menos extremos que los híbridos. Los modelos híbridos, al subestimar la presión térmica, desarrollan gradientes de velocidad angular más pronunciados que pueden inducir inestabilidades de cizalla.

5. Significado e Implicaciones

Este estudio tiene implicaciones profundas para la astrofísica de ondas gravitacionales y la física nuclear:

1. Fiabilidad de las Inferencias: Las relaciones cuasi-universales utilizadas para extraer propiedades de la EOS a partir de las señales de GW deben ser revisadas o calibradas considerando los efectos térmicos completos. Las desviaciones en las fases tardías podrían llevar a errores en la estimación de radios o masas si se usan modelos híbridos simplificados.
2. Detección Futura: La excitación de modos inerciales, con frecuencias más bajas que el modo $f_2$ fundamental, podría ser detectable por detectores de tercera generación (como el Einstein Telescope o Cosmic Explorer). Estos modos son una firma directa de la dinámica térmica y rotacional del remanente.
3. Necesidad de Modelado Preciso: El trabajo subraya que, para simulaciones de larga duración y estudios de estabilidad a largo plazo, es esencial utilizar EOS tabuladas con efectos térmicos consistentes, ya que las aproximaciones híbridas, aunque eficientes, pueden fallar en capturar la física crítica que determina el destino del remanente (estabilidad vs. colapso).
4. Termodinámica de Materia Densa: La discusión sobre la definición de temperatura en fluidos degenerados revela limitaciones en los modelos actuales de EOS híbrida, sugiriendo que la temperatura "efectiva" en estos modelos no es un proxy fiable de la temperatura termodinámica real en el núcleo de las estrellas de neutrones.

En conclusión, el artículo demuestra que el tratamiento de los efectos térmicos no es un detalle menor, sino un factor determinante en la evolución dinámica, la estabilidad y la firma observacional de los remanentes de fusiones de estrellas de neutrones.