Black Hole-Neutron Star Binaries near Neutron Star Disruption Limit in the Mass Regime of Event GW230529

Resumen Técnico: Binarias de Agujero Negro-Estrella de Neutrones cerca del Límite de Disrupción de la Estrella de Neutrones en el Régimen de Masa de GW230529

Planteamiento del Problema
La detección del evento de ondas gravitacionales GW230529 por parte de la Colaboración LIGO–Virgo–KAGRA (LVK) marcó la primera posible fusión de agujero negro–estrella de neutrones (BHNS) con un potencial significativo de producir contrapartidas electromagnéticas (EM). Este evento involucra un agujero negro primario de 2.5–4.5 $M_\odot$ y un objeto compacto secundario de 1.2–2.0 $M_\odot$, situándolo en un régimen de masa donde la estrella de neutrones (NS) puede ser dislocada por marea en lugar de caer directamente en el agujero negro (BH). La disrupción por marea es un requisito previo para generar señales EM observables, tales como kilonovas y brotes de rayos gamma de binarias compactas (cbGRBs), que requieren la formación de un disco de acreción y la eyección de materia rica en neutrones.

Sin embargo, predecir si una fusión BHNS específica resultará en una disrupción depende fuertemente de fórmulas de ajuste analíticas, como el modelo FNH18 (Foucart et al. 2018), que estima la masa del remanente del disco post-fusión. Estos modelos fueron calibrados principalmente en simulaciones con razones de masa $q \gtrsim 3$ y espines alejados del umbral crítico requerido para la disrupción. En consecuencia, existe una falta de datos de alta fidelidad de relatividad numérica (NR) para sistemas cerca del límite de disrupción, particularmente para los parámetros de masa y espín consistentes con GW230529. Esta brecha limita la capacidad de predecir con precisión la prevalencia de señales EM y de restringir la ecuación de estado (EOS) de la estrella de neutrones a partir de observaciones EM.

Metodología
Los autores utilizan el Código de Einstein Espectral (SpEC) para realizar simulaciones de alta resolución de fusiones BHNS. El estudio se centra en una razón de masa de $q=3$ y una masa de estrella de neutrones de $M_{NS} = 1.35 M_\odot$, parámetros consistentes con el evento GW230529. Se emplean dos ecuaciones de estado distintas para probar diferentes niveles de compacidad de la estrella de neutrones:

1. SFHo: Una EOS más blanda con un radio $R_{NS} \approx 11.9$ km.
2. DD2: Una EOS más rígida con un radio $R_{NS} \approx 13.2$ km.

Las simulaciones varían la magnitud del espín del agujero negro ($\chi_{BH}$) para sondear la región inmediatamente circundante al espín crítico ($\chi_{crit}$) predicho por el modelo FNH18 para la disrupción por marea. Para SFHo, $\chi_{crit} \approx 0.037$; para DD2, $\chi_{crit} \approx -0.23$. El estudio incluye simulaciones con espines tanto ligeramente por debajo como por encima de estos umbrales críticos.

Características metodológicas clave:

• Transporte de Neutrinos de Monte Carlo: Esto marca la primera aplicación del transporte de radiación de neutrinos de Monte Carlo en simulaciones de fusión BHNS. A diferencia de los modelos de fuga, este método rastrea las trayectorias de los neutrinos en el espacio de posición y momento para modelar con precisión la absorción, emisión y dispersión. Esto es crucial para determinar la composición (fracción electrónica, $Y_e$) del remanente post-fusión y los flujos de salida, lo que dicta el potencial para la nucleosíntesis de r-proceso.
• Estudios de Resolución: Se realizaron pruebas de convergencia para el caso SFHo en tres resoluciones para cuantificar las incertidumbres numéricas.
• Condiciones Iniciales: Los datos iniciales se generaron utilizando el solver Spells, con excentricidades orbitales reducidas a $e \lesssim 0.009$.

Resultados Clave

1. Masa del Disco vs. Predicciones Analíticas:

   • EOS SFHo: Los resultados numéricos para la masa total del remanente generalmente caen dentro de los errores de modelado esperados de la fórmula FNH18. Sin embargo, las simulaciones muestran consistentemente masas de remanente ligeramente superiores a las que predice el modelo.
   • EOS DD2: Se observa una discrepancia significativa. Para espines cerca del límite de disrupción, el modelo FNH18 subestima la masa del disco remanente por un factor que excede los márgenes de error declarados del modelo. Por ejemplo, a $\chi_{BH} = -0.030$, el modelo predice una masa remanente de $\sim 0.012 M_\odot$, mientras que la simulación arroja $\sim 0.049 M_\odot$.
   • Umbral de Espín Crítico: A pesar de las discrepancias de masa, el valor del espín crítico requerido para activar la disrupción (donde la masa del remanente transiciona de cero a no nula) es consistente entre el modelo analítico y las simulaciones. La diferencia principal radica en la pendiente de la curva de masa del remanente ($dM_{rem}/d\chi_{BH}$) cerca del umbral; las simulaciones muestran un aumento más pronunciado que el modelo.

2. Emergencia de Discos de Baja Masa:

   • Las simulaciones revelan la formación de discos de acreción no nulos incluso para espines ligeramente por debajo del límite de disrupción predicho (por ejemplo, el caso SFHo con $\chi_{BH} = -0.0130$). Aunque estos discos son de baja masa ($\sim 10^{-3} M_\odot$), son potencialmente suficientes para alimentar cortas cbGRBs, como sugieren estudios previos (por ejemplo, Gottlieb et al.).

3. Eyección de Masa y Composición:

   • Eyecta Dinámica: La cantidad de masa eyectada dinámicamente durante la fusión es muy baja, típicamente $\lesssim 0.004 M_\odot$. Para las configuraciones de menor espín, la masa no ligada está cerca de los límites de la resolución numérica, con incertidumbres que potencialmente alcanzan el $\sim 100\%$.
   • Vientos de Disco: El estudio sugiere que para estos sistemas, la masa total eyectada a través de vientos de disco tardíos (impulsados por viscosidad y magnetismo) probablemente dominará sobre la eyección dinámica inicial.
   • Fracción Electrónica ($Y_e$): El eyecta dinámico es altamente rico en neutrones ($Y_e < 0.1$). Sin embargo, los flujos de salida y el material del disco muestran un rango más amplio de $Y_e$ (hasta $\sim 0.3$) debido a interacciones débiles (calentamiento por neutrinos). La EOS DD2 produce eyecta más ricos en neutrones en comparación con SFHo.

4. Implicaciones Observacionales:

   • Cualquier señal de kilonova de estas fusiones específicas tipo GW230529 se predice como tenue, probablemente dominada por la emisión óptica de los vientos de disco en lugar de la emisión infrarroja del eyecta dinámico.
   • La masa total de eyecta disponible para la nucleosíntesis de r-proceso es $\lesssim 0.02 M_\odot$.

Significancia y Reivindicaciones
El artículo sostiene que, si bien el modelo analítico FNH18 sigue siendo ampliamente fiable para predecir la existencia de un umbral de disrupción (el espín crítico), subestima sistemáticamente la masa del disco remanente resultante, particularmente para ecuaciones de estado más rígidas (DD2) cerca del límite de disrupción.

Los autores enfatizan que sus resultados resaltan los límites de la calibración actual cuando se aplica a fusiones BHNS cerca del umbral de disrupción. Específicamente:

• La suposición del modelo de una masa remanente "cero" por debajo del espín crítico no es estrictamente exacta; pueden formarse discos de baja masa incluso por debajo del límite predicho.
• Estos discos de baja masa son físicamente significativos como progenitores potenciales para brotes de rayos gamma cortos, incluso si no producen kilonovas brillantes.
• Los hallazgos sugieren que puede ser necesaria una recalibración modesta de la fórmula de ajuste a medida que haya más simulaciones disponibles en el régimen de baja razón de masa y espín cercano al crítico.

El estudio no pretende alterar definitivamente la evaluación de probabilidad de la LVK sobre si GW230529 fue una fusión disruptiva, señalando que las diferencias observadas difícilmente cambiarán cualitativamente los resultados del análisis de la LVK. Sin embargo, proporciona una base numérica más rigurosa para comprender el comportamiento de las binarias BHNS en este espacio de parámetros específico y observacionalmente crítico.