Inferring Neutron Star Nuclear Properties from Gravitational-Wave and Gamma-Ray Burst Observations

Resumen Técnico: Inferencia de las Propiedades Nucleares de las Estrellas de Neutrones a partir de Observaciones de Ondas Gravitacionales y Brotes de Rayos Gamma

Planteamiento del Problema
Observaciones recientes de brotes de rayos gamma de larga duración (lGRBs) acompañados de emisión de kilonova (por ejemplo, GRB 211211A, GRB 230307A) desafían el paradigma tradicional de que los lGRBs surgen únicamente del colapso de estrellas masivas. Si estos eventos se originan en fusiones de estrellas de neutrones binarias (BNS) o de una estrella de neutrones con un agujero negro (NSBH), surge una pregunta crítica: ¿qué propiedad física gobierna la duración del brote? Los modelos teóricos sugieren que la masa total de la fusión determina la vida útil del remanente de la fusión. Las fusiones más ligeras pueden formar remanentes de estrellas de neutrones (NS) de larga duración que alimentan brotes de rayos gamma cortos (sGRBs), mientras que las fusiones más pesadas resultan en remanentes de corta duración que colapsan en agujeros negros (BHs), lo que potencialmente alimenta lGRBs a través de discos de acreción masivos. Sin embargo, los umbrales de masa específicos que separan estos resultados y su relación con las propiedades nucleares fundamentales, específicamente la masa de Tolman-Oppenheimer-Volkoff ($M_{\text{TOV}}$), permanecen sin restricciones observacionales.

Metodología
Los autores desarrollan un marco bayesiano para inferir escalas de masa características comparando las tasas de fusión de ondas gravitacionales (GW) de la colaboración LIGO-Virgo-Kagra (LVK) con las tasas observadas de sGRBs y lGRBs asociados a kilonovas.

1. Definiciones de Masa Característica: El estudio parametriza tres umbrales de masa crítica relativos a la masa TOV ($M_{\text{TOV}}$):

   • $M_{\text{VL}} = k_1 M_{\text{TOV}}$: El límite entre remanentes de muy larga duración y de larga duración.
   • $M_{\text{LS}} = k_2 M_{\text{TOV}}$: El límite entre remanentes de larga duración y de corta duración. Se hipotetiza que las fusiones con $M_{\text{VL}} < M_{\text{tot}} \le M_{\text{LS}}$ producen sGRBs.
   • $M_{\text{PC}} = k_3 M_{\text{TOV}}$: El límite entre remanentes de corta duración y el colapso inmediato (prompt collapse). Las fusiones con $M_{\text{LS}} < M_{\text{tot}} \le M_{\text{PC}}$ se hipotetiza que producen lGRBs. Las fusiones con $M_{\text{tot}} > M_{\text{PC}}$ experimentan un colapso inmediato; los lGRBs se producen solo si la relación de masa ($q$) permite la formación de un disco suficiente ($1.2 \le q \le 3$ para BNS, $q \le 3$ para NSBH).

2. Fuentes de Datos:

   • Tasas de GW: El catálogo GWTC-4 de LVK (Abac et al. 2025a) y el modelo de población FullPop-4.0 (Abac et al. 2025b) proporcionan tasas de fusión locales en función de las masas de los componentes.
   • Tasas de GRB: Las tasas intrínsecas para sGRBs ($R_{\text{sGRB}}$) y la razón de tasa $\eta = R_{\text{lGRB}}/R_{\text{sGRB}}$ se tratan como restricciones observacionales con incertidumbres significativas, modeladas mediante distribuciones normales.

3. Inferencia: Utilizando el teorema de Bayes, los autores infieren las distribuciones posteriores para el conjunto de parámetros $\vec{\Theta} = (k_1, k_2, k_3, M_{\text{TOV}})$. La función de verosimilitud compara las tasas de GRB inducidas por GW predichas (derivadas de los umbrales de masa) contra las tasas observadas de GRB. Los priors para $M_{\text{TOV}}$ y el radio de la NS en la masa TOV ($R_{\text{TOV}}$) se extraen de Legred et al. (2021), condicionados por observaciones de púlsares y de GW.

Resultados Clave

• Relación de Masa Característica ($k_2$): El análisis arroja un valor mediano de $k_2 = 1.36^{+0.08}_{-0.09}$ (intervalo de confianza del 68%), con $k_2 > 1.24$ con un 90% de confianza. Esto indica que la transición de remanentes de larga duración a corta duración ocurre aproximadamente a $1.36 M_{\text{TOV}}$.
• Implicación para la Vida Útil de los Remanentes: El alto valor de $k_2$ sugiere que las estrellas de neutrones masivas pueden sobrevivir durante períodos prolongados (cientos de milisegundos) después de la fusión, en lugar de colapsar inmediatamente.
• Correlaciones: El estudio identifica una fuerte correlación entre $k_2$ y $k_3$, y una anticorrelación entre $k_2$ y $M_{\text{TOV}}$. Esta anticorrelación permite que las restricciones sobre la relación de masa característica $k_2$ se mapeen directamente a límites superiores de la masa TOV.
• Robustez: La inferencia de $k_2$ es robusta ante una amplia gama de tasas de sGRB asumidas (30 a $2000 \text{ Gpc}^{-3}\text{yr}^{-1}$) y razones de tasa ($\eta$ de 0.01 a 3). Las variaciones en los umbrales de la relación de masa ($q_{\text{low}}$ y $q_{\text{high}}$) y la inclusión o exclusión del límite $M_{\text{VL}}$ tienen un impacto mínimo en la $k_2$ inferida.

Significancia y Reivindicaciones
El artículo establece un método novedoso e independiente para restringir la ecuación de estado (EoS) de las estrellas de neutrones vinculando las poblaciones de fusión de GW con la duración de los GRBs. Al identificar una escala de masa característica ($1.36 M_{\text{TOV}}$) que separa a los progenitores de sGRB y lGRB, los autores proporcionan una restricción observacional directa sobre la masa máxima de las estrellas de neutrones.

Los autores señalan que sus hallazgos son consistentes con simulaciones de relatividad numérica (por ejemplo, Perna et al. 2025) que predicen una transición entre $1.3$y$1.4 M_{\text{TOV}}$. Enfatizan que, si bien los resultados actuales son robustos frente a las incertidumbres observacionales, las mejoras futuras en las mediciones de las tasas de GW y GRB, particularmente las detecciones conjuntas de GRBs con eventos de GW, estrecharán la correlación entre las masas características y $M_{\text{TOV}}$, reduciendo potencialmente la incertidumbre sobre la masa TOV a $\sim 0.1 M_{\odot}$. Este enfoque complementa las técnicas existentes y ofrece una vía para establecer límites más estrictos en la EoS de las NS.