Usando un modelo flexible de mezcla de poblaciones para 259 fusiones binarias de agujeros negros de GWTC-5, se muestra que la función de masa de la subpoblación con espines altos reproduce, pico a pico, la distribución de masas predicha para los remanentes de fusiones de la subpoblación con espines bajos (de origen estelar) hasta ~80 M☉. La coincidencia morfológica, cuantificada por el coeficiente de Bhattacharyya (hasta ~0.95), se explica de forma natural si los agujeros negros de giro rápido son ellos mismos productos de fusiones previas; escenarios alternativos requerirían un ajuste fino de parámetros. Esto constituye una prueba decisiva de las fusiones jerárquicas. Adicionalmente, el corte abrupto superior de masas de la subpoblación de espín bajo en m_max,1 = 54.2^{+7.7}_{-7.2} M☉ proporciona un factor astrofísico S para la reacción 12C(α,γ)16O: S_300 = 151^{+30}_{-26} keV·b (intervalo de confianza del 68%), consistente con el valor teórico de referencia. Toda la población observada de agujeros negros se explica por el colapso estelar seguido de un ensamblaje jerárquico dinámico, sin necesidad de agujeros negros primordiales.
Cómo crecen los agujeros negros de masa estelar después de nacer es una de las preguntas centrales de la astrofísica. Según la relatividad general de Einstein, las masas aceleradas emiten ondas gravitacionales, y Schwarzschild fue el primero en encontrar una solución exacta para un agujero negro sin rotación. Las observaciones con los detectores LIGO, Virgo y KAGRA (creados con la participación clave de Rainer Weiss) revelaron una población de agujeros negros fusionándose con grandes masas y altos espines. Existen dos hipótesis competitivas sobre su origen: crecimiento a través de repetidas fusiones de agujeros negros en cúmulos estelares densos (fusiones jerárquicas) o acreción de gas en agujeros aislados. Distinguirlas significa entender la evolución de las poblaciones estelares y la dinámica de los cúmulos, así como trazar la cadena desde la evolución estelar hasta la formación de objetos compactos masivos.
Los investigadores aplicaron un análisis bayesiano jerárquico a 259 eventos de fusiones de agujeros negros del catálogo de ondas gravitacionales GWTC-5. Construyeron un modelo flexible de mezcla de dos subpoblaciones que describe la distribución de masas, espines y corrimiento al rojo, en el marco de las ecuaciones de Einstein, que predicen las características de la señal. Las funciones de masa se modelaron con splines teniendo en cuenta los efectos de selección y el ruido de los detectores. Se recuperaron distribuciones separadas para la primera subpoblación (bajos espines, interpretados como resultado de la evolución estelar y el colapso del núcleo) y la segunda (altos espines, candidatos a remanentes jerárquicos). La emisión de ondas gravitacionales durante la fusión se describe mediante la fórmula cuadrupolar, lo que permite extraer los parámetros de las fuentes.
El resultado principal es la forma casi idéntica de la función de masa de la segunda subpoblación y la distribución de masas finales después de la fusión de agujeros negros de la primera. Esta similitud se muestra en la figura: en el rango de hasta 80 masas solares, las curvas coinciden prácticamente. El coeficiente de Bhattacharyya, que mide el solapamiento de distribuciones normalizadas, alcanza 0.95 en el intervalo [40,80] M⊙, lo que indica una coincidencia estadísticamente significativa. Los picos de masa de la primera subpoblación a ~10 M⊙ y ~35 M⊙ se reproducen en la segunda como picos a ~17 M⊙ y ~58 M⊙, lo que corresponde a la acumulación esperada de masa durante la fusión. Además, el corte abrupto de la población de bajo espín a una masa de 54.2+7.7-7.2 M⊙ se interpreta como el límite inferior de la brecha de masa causada por supernovas de inestabilidad de pares. Esto permitió calcular el factor S de la reacción 12C(α,γ)16O, clave para la nucleosíntesis en estrellas masivas: S300 = 151+30-26 keV·barn, lo que concuerda excelentemente con las estimaciones teóricas.
La conexión morfológica descubierta es una prueba directa de que los agujeros negros masivos de alto espín son productos de fusiones de agujeros negros anteriores de la primera subpoblación. Un escenario alternativo, como el crecimiento por acreción desde un disco de acreción, requeriría un ajuste fino de parámetros para reproducir las mismas características espectrales. Por lo tanto, el ensamblaje jerárquico domina la formación de la población observada de agujeros negros binarios, y no se requiere la introducción de agujeros negros primordiales para explicarla. La restricción obtenida sobre la masa máxima de los agujeros de bajo espín también desafía los modelos de evolución estelar aislada, que predicen la formación de agujeros negros de rotación lenta con masas de hasta ~93 M⊙.
Con la acumulación de datos de nuevas temporadas de observación de LIGO–Virgo–KAGRA y futuros detectores como el Telescopio Einstein y el Explorador Cósmico, las funciones de masa se medirán con una precisión aún mayor. Esto permitirá no solo reforzar las pruebas de las fusiones jerárquicas de agujeros negros, sino también investigar en detalle la contribución de generaciones superiores (tercera y posteriores) y vincularlas con el nacimiento de agujeros negros de masa intermedia. El modelado de la dinámica de los cúmulos y los flujos de acreción también recibirá nuevas restricciones observacionales.
Los resultados impactarán en la astrofísica de objetos compactos, la teoría de la evolución estelar, la astrofísica nuclear (a través del refinamiento de la nucleosíntesis) y la cosmología, al precisar la tasa de fusiones de agujeros negros y la historia de formación estelar.
Los próximos pasos incluyen la inclusión de los efectos de retraso entre generaciones de fusiones y el análisis de la contribución de diferentes entornos dinámicos (cúmulos globulares, discos de núcleos galácticos activos) mediante simulación numérica detallada y futuras observaciones de ondas gravitacionales.
El trabajo aborda directamente problemas no resueltos: el origen de la función de masa observada de los agujeros negros, la naturaleza de la brecha de masa causada por supernovas de inestabilidad de pares, la velocidad de la reacción nuclear 12C(α,γ)16O y los mecanismos de formación de sistemas binarios en entornos densos.
🎯 La brecha de masa de las supernovas de inestabilidad de pares, donde casi no nacen agujeros negros, se consideró durante mucho tiempo una 'zona prohibida', pero ahora vemos cómo las fusiones jerárquicas la llenan con remanentes de generaciones anteriores, como una matrioska cósmica.
🎬 La idea de fusiones jerárquicas recuerda al 'universo matrioska' de la novela 'El problema de los tres cuerpos' de Liu Cixin, donde las civilizaciones anidan mundos unos dentro de otros, pero aquí la propia naturaleza construye agujeros negros cada vez más masivos a través de una cascada de colisiones.