Convective stability analysis of massive neutron stars formed in binary mergers

Resumen Técnico: Análisis de Estabilidad Convectiva de Estrellas de Neutrones Masivas Formadas en Fusiones Binarias

Planteamiento del Problema
La evolución post-fusión de las fusiones de estrellas de neutrones binarias (BNS) está gobernada por la dinámica de la estrella de neutrones (NS) masiva resultante. Un aspecto crítico de esta evolución es el potencial desarrollo de inestabilidades hidrodinámicas, particularmente la inestabilidad convectiva, que puede impulsar movimientos de fluido a gran escala, alterar la estructura térmica y excitar modos de ondas gravitacionales (GW). Estudios previos han arrojado conclusiones contradictorias: algunas simulaciones utilizando el criterio de Schwarzschild sugirieron que las NS masivas desarrollan una inestabilidad convectiva a gran escala, excitando posteriormente modos inerciales observables [74, 75]. Por el contrario, otros estudios que incorporaron radiación de neutrinos encontraron una estratificación estable [55]. Una limitación clave en los análisis previos es la frecuente negligencia de las fuerzas restauradoras de rotación. Dado que los remanentes de fusión exhiben una rotación rápida y diferencial, la contribución de la rotación a la estabilidad (a menudo descrita por el criterio de Solberg-Høiland) es esencial. Este artículo tiene como objetivo resolver estas discrepancias realizando simulaciones de hidrodinámica totalmente relativistas y aplicando un criterio de estabilidad integral que tenga en cuenta tanto la flotabilidad térmica como los efectos de rotación.

Metodología
Los autores emplean el código de relatividad numérica SACRA-MPI para simular la fusión de sistemas BNS de masas iguales durante una duración de aproximadamente 100 ms post-fusión.

• Ecuación de Estado (EOS): Se utilizan tres EOS politrópicas por tramos (APR4, SLy4, MPA1), aumentadas con un componente térmico de ley $\Gamma$ ($\Gamma_{\text{th}} = 1.8$) para modelar el calentamiento por choque.
• Configuración Numérica: Las simulaciones utilizan una rejilla anidada cartesiana con un refinamiento de 2:1 y simetría de espejo a través del plano orbital. Se prueban tres resoluciones para el modelo APR4 para asegurar la convergencia.
• Análisis de Estabilidad: Los autores derivan y aplican los criterios de estabilidad de Solberg-Høiland relativistas (Criterio I y Criterio II) para estrellas calientes y diferencialmente rotantes. Estos criterios incorporan la interacción entre la flotabilidad térmica (discriminante de Schwarzschild) y las fuerzas restauradoras de rotación (gradientes de momento angular específico). El análisis se realiza sobre los perfiles promediados angularmente de entropía específica y momento angular específico.
• Análisis de Modos: El estudio analiza la excitación de modos de fluido y la emisión de GW, centrándose específicamente en el modo de un solo brazo $m=1$ y la potencial excitación de modos inerciales $m=2$ tras el amortiguamiento del modo $f$ cuadrupolar.

Contribuciones Clave

1. Derivación de Criterios de Estabilidad Relativistas: El artículo proporciona una aplicación rigurosa de los criterios de Solberg-Høiland relativistas (Ecs. 46 y 47) a las NS masivas post-fusión, contabilizando explícitamente la fuerza restauradora de rotación que fue omitida en estudios que dependían únicamente del criterio de Schwarzschild.
2. Marco de Simulación Independiente: Los autores llevan a cabo un análisis independiente utilizando configuraciones numéricas (elecciones de EOS, resoluciones, esquemas hidrodinámicos) distintas de los estudios conflictivos previos, evitando comparaciones directas uno a uno para asegurar una investigación autosuficiente de la física.
3. Pruebas de Simetría y Conservación del Momento: El estudio introduce una prueba de preservación de simetría para verificar que el código numérico no rompa artificialmente la simetría $\pi$. Además, investiga la correlación entre el crecimiento de la amplitud del modo $m=1$ y la violación de la conservación del momento lineal, un control crítico para distinguir inestabilidades físicas de artefactos numéricos.

Resultados

• Estabilidad Convectiva: Las simulaciones no muestran evidencia de inestabilidad convectiva a gran escala. Tanto la entropía específica como el momento angular específico promediados angularmente aumentan hacia afuera dentro de las NS masivas. Aunque regiones locales pueden parecer inestables bajo el criterio de Schwarzschild por sí solo, la inclusión de efectos de rotación (específicamente el aumento hacia afuera del momento angular específico) estabiliza estas regiones. Los criterios de estabilidad general (Criterio I y II) se cumplen para la gran mayoría de la estrella.
• Modos Inerciales: Contrario a los hallazgos en las Refs. [74, 75], no se detecta la excitación de modos inerciales observables tras el amortiguamiento de los modos $f$ cuadrupolares. La señal de GW permanece dominada por el modo $f$ durante toda la simulación, sin que emerjan modos inerciales discretos a frecuencias de kilohertz.
• Inestabilidad $m=1$: Se observa que un modo de un solo brazo $m=1$ crece y se satura, de acuerdo con la literatura previa. Sin embargo, los autores encuentran una fuerte correlación entre el crecimiento de la amplitud de este modo y la violación de la conservación del momento lineal. Esto sugiere que la excitación del modo de un solo brazo en estas simulaciones específicas puede tener un origen numérico más que ser una inestabilidad puramente física.
• Perfiles Térmicos y de Rotación: La NS masiva se establece en un estado de cuasi-equilibrio dentro de los 20–30 ms. El núcleo es frío y rota lentamente, mientras que la envoltura es calentada por interacciones de brazos espirales y rota diferencialmente. El momento angular específico aumenta hacia afuera, proporcionando una fuerte fuerza centrífuga estabilizadora.

Significancia y Reivindicaciones
El artículo sostiene que la inclusión de las fuerzas restauradoras de rotación altera fundamentalmente la evaluación de la estabilidad de las NS masivas post-fusión. Los autores afirman que sus resultados, que muestran una configuración estable, son incompatibles con estudios previos que predecían convección a gran escala y la subsiguiente excitación de modos inerciales. Argumentan que la falta de inestabilidad convectiva en sus simulaciones explica la ausencia de modos inercales excitados.

Respecto al modo $m=1$, los autores mantienen una postura modesta, afirmando que "no pueden rechazar la posibilidad de que la excitación del modo de un solo brazo tenga un origen numérico" debido a su correlación con la no conservación del momento lineal. Enfatizan que, si bien el código pasa las pruebas de preservación de simetría, la conservación del momento lineal en simulaciones con refinamiento de malla sigue siendo un desafío que requiere mayor investigación para establecer definitivamente la naturaleza física de la inestabilidad de un solo brazo.

El estudio concluye que es necesaria una investigación futura que incorpore microfísica (gradientes de composición) y campos magnéticos para proporcionar una conclusión definitiva, pero dentro del alcance de sus simulaciones hidrodinámicas actuales, las NS masivas formadas en fusiones BNS parecen ser convectivamente estables.