Is GW231123 a hierarchical merger?

Título: ¿Es GW231123 una fusión jerárquica?

1. El Problema

El evento de ondas gravitacionales GW231123 representa un desafío significativo para los modelos astrofísicos actuales de formación de binarias de agujeros negros. Se distingue por dos características extremas:

• Masas sin precedentes: Es el evento más masivo observado hasta la fecha, con componentes de $m_1 = 137^{+22}_{-17} M_\odot$ y $m_2 = 103^{+20}_{-52} M_\odot$. Estas masas caen dentro o superan el "hueco de masa" (mass gap) predicho por la física estelar estándar (50–130 $M_\odot$), donde la inestabilidad por pares impide la formación directa de agujeros negros.
• Espines extremos: Mide los espines más altos registrados, con valores de dimensión $\chi_1 = 0.90^{+0.10}_{-0.19}$ y $\chi_2 = 0.80^{+0.20}_{-0.51}$ (intervalos de credibilidad del 90%).

La formación de agujeros negros tan masivos y con espines tan altos a través de la evolución estelar aislada es difícil de explicar debido a la pérdida de momento angular en estrellas masivas y los mecanismos de inestabilidad por pares. Aunque se han propuesto varios escenarios hipotéticos (acreción de gas, estrellas de población III, agujeros negros primordiales), la ensamblaje jerárquico (fusiones de agujeros negros que son ellos mismos productos de fusiones anteriores) es un candidato prometedor, ya que los remanentes de fusiones heredan momento angular orbital, resultando en espines altos ($\chi \approx 0.7$).

Sin embargo, la interpretación de GW231123 está complicada por:

1. Sistemáticas de modelos de onda: Diferentes modelos de waveform (NRSur7dq4 vs. IMRPhenomXO4a) arrojan distribuciones posteriores de espín contradictorias.
2. Ruido no estacionario: La señal es corta (pocos ciclos) y se vio afectada por "glitches" (artefactos de ruido) en los detectores LIGO Hanford y Livingston, lo que podría sesgar la inferencia de parámetros.

2. Metodología

Los autores investigan la hipótesis de que GW231123 es el resultado de una fusión de dos agujeros negros de segunda generación (2G+2G) formados en un cúmulo estelar denso.

• Modelo de Población: Utilizan la distribución de espines de remanentes ($\chi_f$) para agujeros negros de segunda generación obtenida por Borchers et al. (2025) mediante el código Cluster Monte Carlo (CMC). Asumen un espín de nacimiento conservador de $\chi_{birth} = 0.2$ para los agujeros negros de primera generación.
• Formalismo Estadístico: Aplican el método de Passenger et al. (2024) para calcular un valor-p ($p$). Este método compara la evidencia normalizada ($Z$) de un evento observado contra la distribución de evidencia de una población simulada.
  • Se define una estadística $Z$ que compara la probabilidad de que la señal provenga del modelo de población jerárquica ($\pi_{2G}$) frente a un prior uniforme ($\pi_U$).
  • Se simulan eventos inyectados en ruido gaussiano coloreado, analizados con los mismos modelos de waveform que el evento real.
• Análisis de Waveforms: Realizan dos análisis independientes para abordar la incertidumbre del modelo:
  1. IMRPhenomXO4a (XO4a): El modelo favorecido por el factor de Bayes en el análisis original de GW231123.
  2. NRSur7dq4 (NRSur): El modelo mejor calibrado con Relatividad Numérica, aunque menos favorecido por el factor de Bayes en este caso específico.
• Mitigación de Ruido: Utilizan el paquete BayesWave para mitigar el efecto del glitch detectado en el detector H1.

3. Contribuciones Clave

• Prueba cuantitativa de la hipótesis jerárquica: Proporcionan el primer cálculo formal de probabilidad (valor-p) para determinar si los espines extremos de GW231123 son consistentes con un origen de fusión 2G+2G en cúmulos globulares.
• Análisis de dependencia del modelo de waveform: Demuestran cómo la conclusión sobre la naturaleza jerárquica del evento cambia drásticamente dependiendo del modelo de waveform utilizado, resaltando la necesidad crítica de modelos mejor calibrados en el régimen de espines altos.
• Evaluación de la viabilidad en cúmulos: Muestran que, bajo ciertas condiciones de espín de nacimiento y modelos de waveform específicos, la formación jerárquica es una explicación plausible, aunque no la única.

4. Resultados

Los resultados son altamente dependientes del modelo de waveform utilizado:

• Caso IMRPhenomXO4a (XO4a):
  • El valor-p calculado es $p \ll 1\%$.
  • Conclusión: Se rechaza la hipótesis de que GW231123 es una fusión 2G+2G bajo este modelo. El evento cae fuera de la distribución de espines predicha por el modelo de cúmulos para fusiones jerárquicas.
• Caso NRSur7dq4 (NRSur):
  • El valor-p calculado es $p \approx 5\%$.
  • Conclusión: No se puede rechazar la hipótesis jerárquica. Existe una probabilidad razonable (aunque no abrumadora) de que GW231123 sea una fusión 2G+2G. El modelo NRSur permite espines ligeramente menores para el componente primario ($\chi_1$) en comparación con XO4a, lo que lo hace más compatible con la distribución de remanentes de cúmulos.

Interpretación de la discrepancia:
La diferencia se debe a que XO4a favorece fuertemente espines cercanos a la unidad ($\chi \approx 1$), mientras que NRSur permite un rango más amplio que incluye valores más bajos. Dado que los modelos de waveform actuales no están bien calibrados para espines $\chi > 0.8$, la incertidumbre sistemática es el factor dominante.

Además, los autores señalan que GW231123 podría ser una fusión de generaciones mixtas (2G+1G o 3G+2G), especialmente dado el espín extremadamente alto del componente primario, aunque las tasas de fusión para estos escenarios (especialmente 3G) se estiman muy bajas en cúmulos globulares estándar.

5. Significado e Implicaciones

• Validación de Escenarios Exóticos: El estudio sugiere que, a pesar de las sistemáticas, es posible que eventos como GW231123 se originen en ensamblajes jerárquicos en cúmulos estelares densos. Esto abre la puerta a interpretaciones alternativas más allá de la evolución estelar aislada.
• Alerta sobre Sistemáticas: El resultado subraya la fragilidad de las conclusiones astrofísicas cuando se basan en eventos con espines extremos y señales cortas. La dependencia del modelo de waveform indica que se necesitan modelos de waveform mejor calibrados con Relatividad Numérica en el régimen de espines altos para distinguir entre diferentes canales de formación.
• Ruido y Glitches: Reafirma que los artefactos de ruido no estacionario (glitches) pueden sesgar significativamente la inferencia de espines, favoreciendo valores extremos ($\chi \approx 1$). El desarrollo de modelos de ruido más flexibles es crucial para futuros análisis.
• Futuro de la Astrofísica: Si se confirma que GW231123 es un producto de múltiples generaciones de fusiones, proporcionaría una prueba directa de la dinámica de cúmulos estelares y la existencia de agujeros negros de generaciones superiores (2G, 3G), lo cual tendría implicaciones profundas para la comprensión de la dinámica de sistemas estelares densos y la cosmología.

En resumen, el papel no descarta la naturaleza jerárquica de GW231123, pero establece que la evidencia es condicional a la precisión de los modelos de waveform y la gestión del ruido, con una probabilidad de compatibilidad del ~5% bajo el modelo de Relatividad Numérica más robusto.