Probing millisecond magnetar formation in binary neutron star mergers through X-ray follow-up of gravitational wave alerts

Resumen Técnico: Sondeo de la formación de magnetares de milisegundo en fusiones de estrellas de neutrones binarias mediante el seguimiento en rayos X de alertas de ondas gravitacionales

Planteamiento del Problema
La naturaleza del remanente formado tras una fusión de estrellas de neutrones binarias (BNS) sigue siendo incierta. Mientras que las fusiones más masivas se espera que colapsen rápidamente en agujeros negros (BH), las curvas de luz de rayos X de los afterglows de brotes gamma cortos (sGRB) sugieren que una estrella de neutrones (NS) de rotación rápida y altamente magnetizada —un magnetar de milisegundo— puede actuar como el motor central de una fracción significativa de estos transitorios. Distinguir entre una NS estable, una NS supramasiva (SMNS) que colapsa tras la pérdida de rotación (spin-down) y un BH inmediato es crítico para comprender la ecuación de estado (EoS) de la materia ultra-densa y la física de los brotes de rayos gamma. Sin embargo, identificar estos remanentes es un desafío porque las detecciones de ondas gravitacionales (GW) por sí solas no siempre restringen la naturaleza del remanente, y los contrapartes electromagnéticas (EM) suelen estar oscurecidos por el eyecta o ser observados fuera de eje. Este artículo tiene como objetivo evaluar las capacidades actuales y futuras de detección de remanentes de magnetares de milisegundo a través del seguimiento en rayos X de alertas de GW, evaluando específicamente cómo las configuraciones de los detectores de GW y las sensibilidades de los instrumentos de rayos X influyen en las perspectivas de detección.

Metodología
Los autores desarrollaron un marco de simulación integral que combina la generación de señales de GW con un modelado detallado de las curvas de luz de rayos X:

1. Población Sintética de BNS y Detectabilidad de GW:

   • Se generó una población sintética de sistemas BNS utilizando distribuciones de masa derivadas de la ejecución O3 de LIGO-Virgo-KAGRA (específicamente el prior "Peak").
   • Las señales de GW se simularon utilizando el software PyCBC para los interferómetros actuales (LVK O4 y O5) y de próxima generación (Einstein Telescope [ET] y Cosmic Explorer [CE]).
   • La precisión de la localización en el cielo se computó usando BAYESTAR, filtrando eventos bien localizados ($\Delta\Omega \le 50$ deg$^2$) para asegurar un seguimiento viable en rayos X.
   • El estudio consideró varios modelos de EoS (por ejemplo, Shen, BSk21, DDME2) para determinar la masa máxima de la NS no rotante ($M_{TOV}$) y el destino resultante del remanente (NS estable, SMNS o colapso inmediato).

2. Simulación de la Curva de Luz de Rayos X:

   • Un código numérico dedicado modeló la emisión de rayos X de un nuevo magnetar de milisegundo, contabilizando el spin-down por dipolo magnético, los efectos de propulsor (propeller effects) y la emisión de ondas gravitacionales.
   • Modelado del Eyecta: Se implementó un novedoso marco de eyecta multicomponente, tratando el eyecta como una superposición de componentes dinámicos y post-fusión. Las masas y opacidades se derivaron de ajustes polinómicos a simulaciones de relatividad numérica de vanguardia (Nedora et al. 2021a; Mösta et al. 2020a), vinculando las propiedades del eyecta directamente con los parámetros orbitales de la BNS y la EoS.
   • Geometría: La simulación distinguió entre dos ángulos de visión: la "zona libre" (cerca del eje del jet, baja opacidad) y la "zona atrapada" (ecuatorial, alta opacidad). El modelo contabilizó el reprocesamiento de la radiación del magnetar por el eyecta, calculando la profundidad óptica y la emisión de cuerpo negro.
   • Criterios de Colapso: La simulación rastreó la evolución del spin-down, determinando el tiempo de colapso en un BH basado en la EoS y el periodo de rotación crítico.

3. Restricciones Observacionales:

   • Los flujos predichos se compararon con los umbrales de sensibilidad de los instrumentos de rayos X actuales y futuros, específicamente SVOM/MXT, Einstein Probe (EP)/FXT, y EP/WXT.
   • Se estimaron las tasas de detección asumiendo una tasa de fusión BNS de $R = 100$ Gpc$^{-3}$ año$^{-1}$, con incertidumbres estadísticas derivadas mediante métodos bootstrap no paramétricos.

Contribuciones Clave

• Marco de Simulación Integrado: El artículo presenta una simulación autoconsistente que vincula las probabilidades de detección de GW con la evolución detallada de la curva de luz de rayos X, incorporando un modelado de eyecta realista y dependiente de la EoS, en lugar de parametrizaciones simplificadas.
• Análisis de Sesgo del Ángulo de Visión: El estudio cuantifica el compromiso entre la detectabilidad de GW y la visibilidad EM. Demuestra que, si bien los sistemas de la "zona libre" son más brillantes en rayos X, son más difíciles de localizar mediante GW a grandes distancias, mientras que los sistemas de la "zona atrapada" son más numerosos y fáciles de localizar a rangos cercanos, pero sufren de una mayor opacidad del eyecta.
• Ventanas de Detección Óptimas: Los autores identifican ventanas de tiempo específicas para la máxima detectabilidad, señalando que el pico de detectabilidad se desplaza de $\sim$2 horas post-fusión para los detectores actuales a $\sim$3–4 horas para los detectores de próxima generación debido a la dilatación del tiempo cosmológico.
• Estrategia de Instrumentación: El trabajo proporciona estrategias observacionales específicas, abogando por el uso de instrumentos de campo estrecho y sensibles (ej. MXT, FXT) para eventos de zona atrapada y monitores de campo amplio (ej. WXT) para eventos de zona libre, dependiendo de la generación del detector de GW.

Resultados

• Tasas de Formación: Las simulaciones indican que entre el 2% y el 16% de las fusiones BNS forman magnetares de milisegundo de larga duración, dependiendo fuertemente de la EoS (específicamente $M_{TOV}$).
• Tasas de Detección (Generación Actual):
  • Para la ejecución LVK O5, la tasa de detección predicha de magnetares de milisegundo es de $1.0^{+0.3}_{-0.3}$ año$^{-1}$ asumiendo una EoS rígida (DDME2). Si la EoS es blanda (Shen), no se esperan detecciones con los instrumentos actuales.
  • La mejor detectabilidad ocurre aproximadamente dos horas después de la fusión.
  • Para eventos bien localizados, la fracción de magnetares detectables oscila entre el 0% y el 6% (hasta ~12% para la EoS más rígida) una hora después de la fusión.
• Tasas de Detección (Próxima Generación):
  • Con ET y CE, la tasa de detección podría aumentar hasta tres órdenes de magnitud, alcanzando $\sim$750 año$^{-1}$ para escenarios de EoS rígida con ET+CE.
  • Para los detectores de próxima generación, el pico de detectabilidad se desplaza a 3–4 horas post-fusión.
• Dependencia del Campo Magnético:
  • Los magnetares con campos de $\sim 7.5 \times 10^{14}$ G ofrecen la mayor detectabilidad a una hora post-fusión (16% de los eventos bien localizados).
  • Campos más altos ($10^{16}$ G) conducen a un spin-down rápido y un colapso temprano (dentro de minutos), limitando la ventana de observación.
  • Campos más bajos ($10^{14}$ G) tienen vidas más largas pero luminosidades iniciales menores, volviéndose detectables solo horas después de la fusión una vez que el eyecta se vuelve ópticamente delgado.
• Restricciones de la EoS: Los resultados sugieren que la detección (o no detección) de magnetares de milisegundo puede imponer restricciones sobre la EoS. Específicamente, la ausencia de mesetas (plateaus) de rayos X de larga duración favorecería una EoS más blanda o campos magnéticos más bajos.

Significancia y Reivindicaciones
Los autores afirman que la detección de la emisión de spin-down de un magnetar de milisegundo está "al alcance" incluso con las redes de GW actuales, particularmente durante la ejecución O5, siempre que la EoS sea suficientemente rígida. Enfatizan que una sola detección tendría implicaciones profundas para la física de las fusiones, los motores de los brotes de rayos gamma y la física de las estrellas de neutrones.

El artículo concluye modestamente que, aunque la detección es desafiante debido a las incertidumbres en la tasa de fusión, la EoS y la respuesta de los instrumentos, las estrategias de observación informadas por la simulación propuesta —específicamente el tiempo de las observaciones (2–4 horas post-fusión) y la selección de instrumentos basada en el ángulo de visión (zona atrapada vs. zona libre)— aumentan significamente la probabilidad de éxito. El trabajo subraya la necesidad de esfuerzos teóricos y observacionales sostenidos para resolver la naturaleza de los remanentes de las fusiones BNS.