GW231123 Formation from Population III Stars: Isolated Binary Evolution

Resumen Técnico: Formación de GW231123 a partir de Estrellas de Población III

1. Planteamiento del Problema

El evento de ondas gravitacionales GW231123 representa un hito en la astronomía de ondas gravitacionales, ya que involucra la fusión de dos agujeros negros (AN) con masas estimadas superiores a $100 M_\odot$ ($137^{+22}_{-17} M_\odot$ y $103^{+20}_{-51} M_\odot$). Estas masas sitúan a los componentes dentro o por encima del "hueco de masa de inestabilidad de pares" (PI gap, $\sim 60-130 M_\odot$), un rango donde la formación de agujeros negros mediante colapso estelar estándar se considera teóricamente improbable debido a la explosión total de la estrella (supernova de inestabilidad de pares, PISN).

Además, ambos agujeros negros exhiben espines potencialmente altos. El origen de tales eventos masivos es objeto de debate, con canales propuestos que incluyen fusiones jerárquicas en cúmulos estelares, discos de núcleos galácticos activos (AGN) o evolución de binarias aisladas. Este estudio se centra en evaluar si la evolución de binarias aisladas de estrellas de Población III (estrellas de metalicidad cero formadas en el universo temprano) puede explicar la existencia de GW231123.

2. Metodología

Los autores emplearon cálculos de síntesis de poblaciones binarias utilizando el código BSEEMP, una derivación del código BSE original adaptada para metalicidades extremadamente bajas ($Z = 2 \times 10^{-10}$), simulando así la evolución de estrellas de Población III.

• Población Inicial: Se simularon 3 millones de sistemas binarios. Se asumió una función de masa inicial (IMF) "top-heavy" (pesada en masas altas) con un índice de ley de potencias de -1 en el rango de $10-150 M_\odot$. La fracción binaria se fijó en 0.5, con distribuciones de periodos orbitales, relaciones de masa y excentricidades basadas en observaciones de estrellas masivas de Población I/II (Sana et al. 2012).
• Modelos de Evolución Estelar: Se compararon dos modelos de evolución de estrellas individuales para Población III basados en simulaciones hidrodinámicas 1D (código HOSHI):
  • Modelo M: Sobrepasos convectivos ineficientes (radios máximos pequeños, $\lesssim 10^2 R_\odot$).
  • Modelo L: Sobrepasos convectivos eficientes (radios máximos grandes, $\gtrsim 10^3 R_\odot$).
• Física Nuclear y Supernovas: Se varió la tasa de la reacción nuclear $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$, un parámetro crítico que determina los umbrales de masa para las supernovas de inestabilidad de pares (PISN) y pulsacionales (PPISN). Se probaron tasas estándar y tasas reducidas en $1\sigma$, $2\sigma$ y $3\sigma$.
• Definición de Eventos: Se definieron "eventos tipo GW231123" como fusiones donde las masas de los componentes y la relación de masas caen dentro de los intervalos de credibilidad del 90% reportados por LIGO/Virgo para GW231123.

3. Contribuciones Clave y Resultados

El estudio identificó las condiciones físicas estrictas necesarias para que la evolución de binarias aisladas de Población III produzca fusiones masivas como GW231123:

• Condiciones Necesarias: La formación de eventos tipo GW231123 es posible solo si se cumplen simultáneamente tres condiciones:
  1. Las estrellas de Población III evolucionan con sobrepasos convectivos ineficientes (Modelo M).
  2. La tasa de reacción $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ es al menos $2\sigma$ inferior al valor estándar.
  3. Las binarias de Población III comparten parámetros orbitales iniciales similares a las de Población I/II (periodos cortos).
• Tasas de Fusión: Bajo las condiciones anteriores (Modelo M con tasa $2\sigma$ o $3\sigma$ más baja), la tasa de fusión local calculada es de $\sim 2.1 \times 10^{-2} \text{ Gpc}^{-3} \text{ yr}^{-1}$. Esto es consistente con la estimación de un solo evento reportada por LIGO ($0.08^{+0.19}_{-0.07} \text{ Gpc}^{-3} \text{ yr}^{-1}$).
• Contraste con el Modelo L: En el modelo con sobrepasos convectivos eficientes (Modelo L), incluso con la tasa de reacción reducida ($3\sigma$), no se producen eventos tipo GW231123. Las estrellas masivas en este modelo se expanden enormemente, sufren interacciones binarias intensas (envolvente común) que reducen drásticamente sus masas finales a $\sim 70 M_\odot$, o mantienen separaciones orbitales demasiado grandes para fusionarse en el tiempo de Hubble.
• Espines y Parámetros Orbitales:
  • El modelo predice que los agujeros negros primarios tienen espines bajos ($a_1 \sim 0$) debido a la pérdida de momento angular durante la transferencia de masa estable, mientras que los secundarios pueden tener espines altos ($a_2 \sim 1$).
  • El espín efectivo de inspiración ($\chi_{\text{eff}}$) y la relación de masas de los eventos sintetizados caen dentro de los intervalos de credibilidad observados para GW231123.
  • Las órbitas iniciales de los agujeros negros binarios tienen excentricidades bajas ($\lesssim 0.5$) y semiejes mayores de $\sim 100 R_\odot$, lo que implica que circularizarán completamente antes de entrar en la banda de frecuencias de observatorios espaciales como LISA.

4. Significado e Implicaciones

• Restricción de Modelos Estelares: El descubrimiento de GW231123 impone restricciones mucho más estrictas que eventos anteriores como GW190521. Mientras que GW190521 podría formarse con tasas de reacción estándar o ligeramente reducidas, GW231123 requiere ineficiencia convectiva y una tasa de reacción nuclear significativamente reducida. Esto descarta grandes regiones del espacio de parámetros de los modelos de evolución estelar de una sola estrella.
• Discriminación de Canales de Formación: La predicción de órbitas circulares y espines específicos en el canal de evolución de binarias aisladas contrasta con el canal de formación dinámica en cúmulos estelares (discutido en un artículo complementario), que podría producir eventos con alta excentricidad. La detección (o no) de excentricidad orbital en futuros observatorios (LISA, TianQin) podría distinguir entre estos orígenes.
• Validación de Población III: El estudio demuestra que, bajo condiciones físicas específicas (nucleares y de mezcla convectiva), las estrellas de Población III son candidatos viables y suficientes para explicar la población de fusiones de agujeros negros ultra-masivos observada, sin necesidad de recurrir exclusivamente a mecanismos dinámicos complejos.

En conclusión, el artículo establece que la detección de GW231123 actúa como un "laboratorio cósmico" que valida o refina parámetros fundamentales de la física estelar (como la tasa de reacción $^{12}\text{C}(\alpha, \gamma)^{16}\text{O}$ y la eficiencia del sobrepaso convectivo) bajo la hipótesis de evolución binaria aislada.