El origen de los agujeros negros masivos de giro rápido es una cuestión abierta: fusiones jerárquicas en cúmulos densos o acreción en agujeros aislados. El análisis de 259 eventos de GWTC-5 mediante un modelo de mezcla muestra que la distribución de masas de la submuestra de alto espín reproduce con gran precisión la distribución esperada de masas residuales para la submuestra de bajo espín (origen estelar) hasta ~80 M☉. La cercanía de estas funciones (coeficiente de Bhattacharyya ~0.95) es una evidencia decisiva de fusiones jerárquicas, mientras que otros mecanismos requerirían un ajuste fino. El corte superior de las masas de los agujeros lentos en 54.2 (+7.7/-7.2) M☉ permite extraer el factor astrofísico S_300 = 151 (+30/-26) keV·b para la reacción 12C(α,γ)16O, lo que coincide con las predicciones teóricas. Así, toda la población de agujeros negros observados se explica por el colapso estelar y el posterior ensamblaje jerárquico, sin participación de agujeros negros primordiales.
La cuestión de cómo los agujeros negros de masa estelar ganan peso tras su formación sigue siendo uno de los misterios centrales de la astrofísica. En los últimos años, los observatorios de ondas gravitacionales LIGO, Virgo y KAGRA han registrado cientos de fusiones, entre las que destaca un grupo de objetos masivos con alta rotación propia. Se suponía que estos agujeros negros podían nacer en cúmulos densos, fusionándose sucesivamente entre sí, o crecer por acreción en sistemas aislados. Hasta ahora no se había logrado una separación clara de estos escenarios: las señales son demasiado difusas. Un nuevo estudio, basado en el catálogo GWTC-5, ofrece por primera vez una prueba directa a favor de las fusiones jerárquicas, comparando las funciones de masa de distintas subpoblaciones.
Los autores aplicaron un análisis bayesiano jerárquico a 259 eventos de fusión de agujeros negros binarios del catálogo GWTC-5. Se utilizó un modelo mixto flexible (dos subpoblaciones) con descripciones separadas de las distribuciones de masas, relaciones de masa y parámetros de espín. Las funciones de masa se modelaron con splines, lo que permitió revelar estructuras finas: picos y valles. Para cada subpoblación se reconstruyeron independientemente las distribuciones del espín efectivo χeff, el espín precesante χp y la tasa de fusiones. Luego, para la subpoblación de bajo espín (interpretada como resultado del colapso estelar) se calculó la distribución de masas residuales tras las fusiones (aproximación m_final = 0.95(m1+m2)). La similitud entre la forma de la distribución obtenida y el espectro de masas de los agujeros negros de alto espín se evaluó mediante el coeficiente de Bhattacharyya, las distancias de Jensen-Shannon y Wasserstein.
El resultado principal es una coincidencia casi perfecta entre la forma de la distribución de masas de la subpoblación de alto espín y la distribución de masas residuales de las fusiones de objetos de la subpoblación de bajo espín en el rango de ~20 a ~80 M⊙. El coeficiente de Bhattacharyya alcanzó un valor de ~0.95, y la similitud se mantiene para diferentes umbrales. Esto significa que los picos y valles de un espectro reproducen uno a uno la estructura del otro. Si los agujeros negros de alto espín crecieran por acreción, su función de masa no estaría tan morfológicamente ligada a la primera. Por tanto, esto es una prueba directa del origen jerárquico: los objetos masivos de rápida rotación son descendientes de fusiones previas. Adicionalmente, el análisis del corte de masas de la subpoblación de bajo espín (mmax,1 = 54.2+7.7−7.2 M⊙) proporcionó una restricción sobre el factor S astrofísico de la reacción 12C(α,γ)16O: S300 keV = 151+30−26 keV·b. Este valor concuerda magníficamente con los cálculos teóricos y confirma que el corte está relacionado con el límite superior de los colapsos estelares debido a la supernova de pares inestable (pulsacional). La fracción de fusiones jerárquicas obtenida es de aproximadamente el 23%, y aumenta para masas >30 M⊙.
El trabajo zanja un largo debate: los agujeros negros de alto espín realmente nacen de fusiones múltiples, no de acreción. Esto cambia fundamentalmente la comprensión del crecimiento de los agujeros negros en cúmulos galácticos y núcleos activos. Además, la identificación fiable de la subpoblación de origen estelar proporciona una «regla» precisa para la astrofísica nuclear, permitiendo calibrar la tasa de la reacción de síntesis del carbono — la famosa reacción del «Santo Grial» que determina la evolución de las estrellas masivas y la composición química del Universo. Finalmente, del análisis se deduce que no son necesarios agujeros negros primordiales como candidatos a materia oscura para interpretar los datos de GWTC-5.
A medida que se acumulen eventos de ondas gravitacionales, especialmente con detectores de nueva generación (como el Einstein Telescope), será posible la detección directa de descendientes en segunda e incluso tercera generación de fusiones. Surgirán mapas detallados de distribución de masas de hasta ~100 M⊙ y más, lo que permitirá estudiar el ensamblaje jerárquico de múltiples pasos en cúmulos y buscar semillas de agujeros negros de masa intermedia. Una medición más precisa de mmax,1 reducirá la incertidumbre del factor S a un nivel comparable a la precisión de experimentos de laboratorio, uniendo la física nuclear y la astronomía de ondas gravitacionales.
Los resultados impactarán la astrofísica de poblaciones estelares, la física de reacciones nucleares en estrellas y la cosmología, limitando la contribución de los agujeros negros primordiales a la materia oscura.
Es necesario medir con mayor precisión la cola de alta masa de la distribución para confirmar la existencia de agujeros negros de tercera generación, y refinar la relación entre la forma del «valle» en las masas y los parámetros de las supernovas de pares inestables con diferente metalicidad.
El trabajo conecta problemas no resueltos del crecimiento de los agujeros negros, el origen del corte superior de masas (brecha de pares inestables) y la determinación de la tasa de la reacción clave de combustión del carbono, uno de los principales «escollos» en la teoría de evolución de estrellas masivas.
🎯 El coeficiente de Bhattacharyya es una medida de solapamiento entre dos distribuciones, habitualmente utilizado en procesamiento de señales. ¡Aquí se ha aplicado por primera vez para demostrar que dos tipos de agujeros negros son «antepasado» y «descendiente», alcanzando un valor récord de 0.95!