El análisis de la distribución del espín efectivo (χeff) de agujeros negros binarios en el catálogo GWTC-5.0 mediante inferencia bayesiana jerárquica reveló una transición en la masa del agujero negro principal. La población se modela como una mezcla de dos componentes de espín con una masa frontera extraída de los datos: 15,2(+4,3/−3,6) M⊙. El análisis de catálogos artificiales mostró que es poco probable que esta transición sea causada por fluctuaciones de la muestra finita. Por debajo de la masa de transición, la distribución de χeff es estrecha, con un pico en un valor positivo pequeño, pero con un soporte significativo para espines negativos; por encima, la distribución es más ancha, el pico se desplaza hacia χeff≈0 y el soporte para valores positivos y negativos es aproximadamente igual. Estos resultados apuntan a poblaciones estadísticamente diferentes: el grupo dominante con masas de ~10 M⊙ probablemente se forma de otra manera. La población de baja masa es compatible con la formación a partir de estrellas múltiples masivas en el campo, ya sea mediante la evolución de binarias aisladas con grandes patadas natales de agujeros negros (≳100 km/s) por debajo de la masa frontera, o mediante la evolución dinámica de sistemas triples jerárquicos; los modelos estándar de evolución binaria aislada no reproducen el soporte de χeff negativos.
Las fusiones de agujeros negros binarios, predichos por Karl Schwarzschild, son la principal fuente de ondas gravitacionales registradas por detectores como LIGO. El origen de estos sistemas sigue siendo uno de los problemas centrales de la astrofísica. El espín efectivo χeff es un observable clave que contiene información sobre las masas, la rotación y la orientación de los ejes de los componentes. Diferentes canales de formación predicen distintos comportamientos del espín: desde ordenados en la evolución aislada hasta caóticos en cúmulos densos. Comprender estas diferencias es esencial para reconstruir la evolución de estrellas masivas y la dinámica de poblaciones estelares en las galaxias.
En el trabajo se utilizaron datos de 259 eventos confiables del catálogo GWTC-5.0, registrados por observatorios de ondas gravitacionales. El análisis poblacional se realizó mediante inferencia jerárquica bayesiana, modelando la distribución de χeff como una mezcla de dos componentes separadas por una masa de transición ˜m, derivada directamente de los datos. Los componentes se describieron de manera no paramétrica mediante procesos gaussianos, lo que evitó suposiciones previas sobre la forma. Para verificar la robustez, también se utilizó un modelo paramétrico. Para descartar artefactos, se simularon 100 catálogos artificiales con una distribución suave del espín independiente de la masa.
Se detectó una transición marcada a una masa ˜m = 15,2+4,3–3,6 M⊙. Por debajo de este valor, la distribución de χeff es estrecha (anchura aproximada de 0,05) con un pico alrededor de +0,03, aunque con una fracción notable de valores negativos. Por encima de la transición, la distribución se vuelve significativamente más ancha (∼0,12), el pico se desplaza hacia cero y los espines positivos y negativos aparecen casi con la misma frecuencia. La diferencia entre las dos poblaciones es estadísticamente significativa a un nivel >99% para las medianas. Las simulaciones mostraron que tal transición no puede surgir de fluctuaciones aleatorias si la distribución del espín es independiente de la masa. Es importante destacar que el resultado no depende de la inclusión de eventos inusuales particulares, como GW241011 y GW241110.
Los resultados muestran que la población dominante de fusiones con masas alrededor de 10 M⊙ es estadísticamente distinta de las demás. Esto pone en duda la universalidad de los escenarios de evolución aislada de binarias, que en su forma estándar no pueden reproducir la fracción de χeff negativos. Se requieren, para explicarlo, ya sea grandes 'patadas' natales durante el colapso a agujeros negros (≳100 km/s), lo que está relacionado con la física de las supernovas, o efectos dinámicos en sistemas triples. Así, el origen de los agujeros negros de baja masa resulta estar estrechamente ligado a las etapas finales de la evolución de estrellas masivas, donde desempeñan un papel importante el hidrógeno y el helio en las envolturas, así como la síntesis de elementos hasta el carbono.
En el futuro, con la acumulación de datos de los detectores actuales y planificados (LIGO, Virgo, KAGRA, así como los de tercera generación como Einstein Telescope y Cosmic Explorer), será posible precisar la estructura fina de la transición y estudiar su dependencia con el corrimiento al rojo. Un mejor modelado del colapso estelar y la hidrodinámica de las supernovas ayudará a esclarecer los mecanismos de formación de agujeros negros con grandes impulsos. También se esperan avances en simulaciones numéricas directas de sistemas triples y múltiples, lo que permitirá contrastar las predicciones con los datos.
El trabajo abarca un amplio espectro de campos: desde la física de ondas gravitacionales y la astrofísica de objetos compactos hasta la teoría de evolución de estrellas masivas y la astronomía dinámica. También impulsa investigaciones en astrofísica nuclear, donde los procesos durante el colapso del núcleo son cruciales.
Los próximos pasos incluyen una comparación detallada de las predicciones de síntesis poblacional de sistemas triples con una muestra ampliada de eventos, así como la búsqueda de marcadores adicionales, como la distribución de las relaciones de masa y la excentricidad. Se planea una investigación más exhaustiva de la dependencia con la metalicidad y el entorno.
El límite de espín descubierto está directamente relacionado con preguntas no resueltas: ¿cuáles son las velocidades típicas de las patadas natales de los agujeros negros? ¿Con qué frecuencia se forman estrellas masivas en sistemas triples o más complejos? ¿Con qué eficacia la interacción de marea sincroniza la rotación en pares cercanos? Las respuestas a estas preguntas llenarán los vacíos en nuestra comprensión de las etapas finales de la evolución estelar y la naturaleza de la población de fuentes de ondas gravitacionales.
🎯 La masa del Sol es aproximadamente 2×10³⁰ kg, y un agujero negro típico de 15 M⊙ tiene un radio del horizonte de eventos de solo unos 45 kilómetros, menos que la distancia de Madrid a Toledo.
🎬 En la novela 'Contacto' de Carl Sagan y la película homónima, la señal proveniente del sistema Vega contenía patrones de ondas gravitacionales encriptados. Hoy nosotros mismos desciframos las 'señales' desde las profundidades del Universo para entender cómo funcionan las fábricas de agujeros negros.