First eccentric inspiral-merger-ringdown analysis of neutron star-black hole mergers

Resumen Técnico: Primer Análisis de Inspiración-Fusión-Relajación Excéntrica de Fusiones de Estrella de Neutrones-Agujero Negro

Problema y Contexto
La detección de ondas gravitacionales (GW) de coalescencias de binarias compactas ha permitido el estudio de canales de formación astrofísica. Mientras que las binarias de masa estelar formadas mediante evolución aislada se espera que se circularicen antes de entrar en la banda del detector, las interacciones dinámicas en entornos densos (p. ej., cúmulos estelares) o sistemas triples pueden retener excentricidad orbital. La detección de esta excentricidad es crucial para distinguir los canales de formación. Análisis previos de la evento NSBH GW200105 utilizando un modelo de forma de onda excéntrico, de espín precesante y solo de inspiración (Morras et al. 2025a), proporcionaron evidencia sólida de excentricidad orbital. Sin embargo, aún no se había realizado un análisis completo que incorporara las fases completas de inspiración, fusión y relajación (IMR), incluyendo armónicos subdominantes, para sistemas NSBH. Además, los modelos IMR excéntricos existentes suelen carecer de calibración con relatividad numérica (NR) en el sector excéntrico o son computacionalmente prohibitivos para la estimación de parámetros (PE) de sistemas de baja masa.

Metodología
Los autores realizan el primer análisis de estimación de parámetros IMR completo de tres eventos NSBH: GW200105, GW200115 y GW230529. El análisis utiliza el modelo de forma de onda IMRPhenomTEHM, un modelo excéntrico de espín alineado que extiende al modelo IMRPhenomTHM cuasi-circular (QC). Las características clave de la metodología incluyen:

• Modelo de forma de onda: IMRPhenomTEHM incluye el modo dominante $(2, \pm 2)$ y armónicos subdominantes $(l, m) = \{(2, \pm 1), (3, \pm 3), (4, \pm 4), (5, \pm 5)\}$. Está calibrado con simulaciones de NR cuasi-circulares hasta una relación de masa de 1/18 y asume que la binaria se circulariza en la fusión. No incluye efectos de materia (disrupción de marea), los cuales se consideran despreciables para las relaciones de masa y las relaciones señal-ruido (SNR) estudiadas.
• Marco de Inferencia: La inferencia bayesiana se lleva a cabo utilizando el marco bilby con el algoritmo de muestreo anidado dynesty. El análisis muestrea un espacio de parámetros de 13 dimensiones que incluye masas de los componentes, espines alineados, excentricidad de referencia ($e_{20\text{Hz}}$), anomalía media y parámetros extrínsecos.
• Priors (Distribuciones a priori): Se prueban dos priors de excentricidad: un prior uniforme ($e \in [0, 0.5]$) y un prior log-uniforme ($e \in [10^{-4}, 0.5]$).
• Control de Sistemáticas: Los autores abordan los desafíos computacionales específicos de los modelos excéntricos en el dominio del tiempo, particularmente la necesidad de comenzar la generación de la forma de onda lo suficientemente temprano como para capturar todos los armónicos relevantes en la frecuencia de inicio de la integración de la verosimilitud (20 Hz). Prueban frecuencias de inicio de 17 Hz, 20 Hz y 13.3 Hz, y varían las duraciones de la señal (32s frente a 64s) para evaluar la resolución de baja frecuencia y la contaminación por ruido.

Resultados Clave

1. GW200105 (Evidencia de Excentricidad): El análisis confirma evidencia sólida de excentricidad orbital no nula.

   • Usando un prior uniforme, el pico de la posterior se sitúa en $e_{20\text{Hz}} = 0.12^{+0.02}_{-0.03}$, con un factor de Bayes de $\log_{10} B_{E/QC} = 1.22^{+0.12}_{-0.12}$ que favorece la hipótesis excéntrica sobre la hipótesis QC.
   • Usando un prior log-uniforme, el modo primario se desplaza hacia cero, pero un modo subdominante permanece en $e_{20\text{Hz}} \sim 0.11$, con un factor de Bayes modesto de $\log_{10} B_{E/QC} = 0.11^{+0.11}_{-0.11}$.
   • La posterior excéntrica exhibe bimodalidad, con picos en $e \sim 0.11$ y $e \sim 0.14$. Esto se atribuye a las degeneraciones entre la excentricidad, las masas de los componentes y el espín alineado, que influyen en la duración de la señal. El modo de mayor excentricidad se correlaciona con una masa primaria más baja ($m_1 \sim 8 M_\odot$) y un espín más anti-alineado.
   • Las pruebas variando las frecuencias de inicio y las duraciones confirman la robustez de la bimodalidad y la señal de excentricidad no nula, aunque los autores señalan que los límites de resolución de baja frecuencia (debido a la ventana de datos de 32s) pueden influir en la estructura de la posterior.

2. GW200115 y GW230529 (Consistencia con QC):

   • Ambos eventos son consistentes con binarias cuasi-circulares.
   • GW200115 arroja $e_{20\text{Hz}} = 0.02^{+0.03}_{-0.02}$ con $\log_{10} B_{E/QC} = -0.91$.
   • GW230529 arroja $e_{20\text{Hz}} = 0.01^{+0.02}_{-0.01}$ con $\log_{10} B_{E/QC} = -1.36$.
   • Estos resultados indican que no hay evidencia de excentricidad en estos eventos específicos.

3. Degeneraciones de Parámetros: La inclusión de la excentricidad en el modelo induce cambios en los parámetros QC recuperados (masas y espines) en comparación con los análisis solo de QC. Para GW200105, el modelo excéntrico favorece masas más bajas para compensar el acortamiento de la señal causado por la excentricidad.

Significancia y Reivindicaciones
El artículo afirma presentar el primer análisis de estimación de parámetros de espín alineado excéntrico IMR completo de fusiones NSBH. Su principal significancia radica en:

• Confirmación: Validar la afirmación de excentricidad para GW200105 utilizando un modelo IMR completo con multipolos superiores, reforzando los hallazgos del análisis de solo inspiración (Morras et al. 2025a).
• Avance Metodológico: Demostrar que modelos excéntricos IMR computacionalmente eficientes pueden aplicarse a sistemas de baja masa (NSBH) para medir la excentricidad, a pesar de la falta de calibración directa de NR en el sector excéntrico.
• Información Poblacional: Proporcionar evidencia de que al menos una fusión NSBH confirmada (GW200105) probablemente se formó a través de un canal dinámico en lugar de la evolución binaria aislada, mientras que otras (GW200115, GW230529) parecen consistentes con canales de formación estándar.
• Conciencia de Sistemáticas: Destacar los desafíos computacionales y sistemáticos de la PE de tiempo de la excentricidad, específicamente respecto a la resolución de baja frecuencia, la duración de la señal y la interacción entre las degeneraciones de excentricidad y espín/masa. Los autores señalan que los problemas de calidad de datos (glitches de ruido) en GW200105 y GW200129 (otro candidato excéntrico) requieren una evaluación cuidadosa en futuros análisis.

Los autores concluyen que, si bien las incertidumbres sistemáticas (como la bimodalidad y los límites de resolución de baja frecuencia) persisten, los resultados marcan un paso hacia mediciones de excentricidad robustas en la astronomía de ondas gravitacionales.