Population of Binary Black Holes Inferred from One Hundred and Fifty Gravitational Wave Signals

Resumen Técnico: Población de Agujeros Negros Binarios Inferida a partir de Ciento Cincuenta Señales de Ondas Gravitacionales

Planteamiento del Problema
Con la publicación del cuarto catálogo de ondas gravitacionales (GWTC-4) por parte de las colaboraciones LIGO–Virgo–KAGRA (LVK), el número de observaciones de Sistemas Binarios de Agujeros Negros (BBH) ha superado los 150, duplicando con creces el conjunto de datos disponible en catálogos anteriores. Un objetivo central de la astronomía de ondas gravitacionales es reconciliar los procesos físicos de formación de binarias compactas con las distribuciones observadas de masa, espín y corrimiento al rojo (redshift). Si bien estudios previos han identificado estructuras emergentes en la distribución de masa de los BBH —específicamente picos alrededor de 8–10 $M_\odot$ (masa de chirp) y 14 $M_\odot$—, la significancia estadística de estas características y su interpretación siguen siendo objeto de debate. Además, la distribución de espín ha resultado difícil de constreñir debido a las grandes incertidumbres y las señales débiles. El desafío radica en desarrollar metodologías que infieran con precisión las propiedades de la población y las conecten con canales de formación específicos, tales como las fusiones jerárquicas, que predicen firmas distintivas de masa y espín que aún no se han confirmado plenamente a nivel de población.

Metodología
Este estudio utiliza el marco de modelo de mezcla Vamana para inferir las propiedades de la población de 153 observaciones de BBH (seleccionadas con una tasa de falsa alarma $\le$ 1 año$^{-1}$ y una masa secundaria media $> 3 M_\odot$).

• Parámetros del Modelo: A diferencia de versiones anteriores de Vamana que utilizaban la masa de chirp y la relación de masas, esta iteración infiere la masa primaria ($m_1$) y la masa secundaria ($m_2$) directamente. El modelo también infiere los componentes de espín alineados (asumidos como distribuidos de forma idéntica e independiente) y la evolución con el corrimiento al rojo de la tasa de fusión, modelada como $R(z) = R_0(1+z)^\kappa$.
• Enfoque Estadístico: La población se modela utilizando Gaussianas multivariantes para las distribuciones de masa y espín, combinadas con una ley de potencia para la evolución del corrimiento al rojo. Esta inferencia de cuatro dimensiones permite capturar las correlaciones entre parámetros.
• Manejo de Datos: El análisis se basa exclusivamente en eventos reportados por las colaboraciones LVK para asegurar estimaciones de sensibilidad precisas. El estudio también incorpora datos de observaciones recientemente publicadas (GW241110 124123 y GW241011 233834) para análisis específicos de alto espín.
• Subpoblación de Alto Espín: Para investigar las firmas de fusiones jerárquicas, los autores aíslan una subpoblación de BBH con un valor medio de magnitud de espín efectivo $|\langle\chi_{\text{eff}}\rangle| > 0.2$ e infieren sus distribuciones de masa de componentes de forma separada.

Contribuciones Clave y Resultados

1. Distribución de Masa y Picos:

   • La distribución de masa de chirp inferida exhibe tres picos distintos aproximadamente en 8 $M_\odot$, 14 $M_\odot$ y 27–28 $M_\odot$. El primer pico (que abarca de 6 a 10 $M_\odot$) contribuye aproximadamente con el 66% de la tasa de fusión astrofísica.
   • Existe un "vacío" notable entre los dos primeros picos (8.9–11.7 $M_\odot$), donde la densidad media cae significativamente.
   • Crucialmente, mientras que la distribución de masa de chirp muestra estos tres picos con alta confianza (superando el 99% para los dos primeros), las distribuciones de masa primaria y secundaria no muestran picos comparables con la misma confianza. Esto implica una correlación única entre las masas de los componentes que produce los picos en la masa de chirp sin correspondientes picos en las distribuciones de masa individuales.
   • Los picos están separados por un factor de aproximadamente 1.9 (contabilizando la pérdida de masa por ondas gravitacionales), lo cual es consistente con un factor de escala multiplicativo.

2. Correlaciones de Espín y Relación de Masas:

   • La distribución global de la relación de masas es mayormente uniforme entre 0.6 y 0.9, sin evidencia fuerte de un emparejamiento preferencial de binarias de masa igual.
   • Los espines alineados son generalmente pequeños (con un pico cerca de 0.05), siendo la mayoría de las observaciones consistentes con magnitudes de espín bajas.
   • Sin embargo, existe una correlación donde las binarias más pesadas tienden a tener espines alineados mayores. Para masas de chirp superiores a 35 $M_\odot$, el intervalo de credibilidad del 90% para los espines alineados asciende a 0.1–0.4.
   • La relación de masa para el segundo pico de masa de chirp (14 $M_\odot$) varía significamente, formando una correlación dentada con la masa primaria, mientras que el tercer pico (27 $M_\odot$) consiste principalmente de binarias con masas comparables.

3. Evidencia de Fusiones Jerárquicas:

   • El espaciamiento de los picos de masa de chirp (factor de ~1.9) se alinea con el escenario de fusión jerárquica, donde los remanentes de fusiones de menor masa forman binarias de mayor masa.
   • Aunque la población global carece de los altos espines y las relaciones de masa específicas (p. ej., 0.5, 0.25) predichas por modelos jerárquicos ingenuos, una subpoblación específica de alto espín ($|\langle\chi_{\text{eff}}\rangle| > 0.2$) exhibe estas firmas esperadas.
   • Las masas de los componentes de estas binarias de alto espín se agrupan alrededor de valores que son complementarios a los picos en la distribución de la población completa. Específicamente, los componentes de alto espín se alinean con la escala de masa teórica de fusiones 1G, 2G y 3G (p. ej., ~8.5 $M_\odot$, ~16.2 $M_\odot$, ~30.7 $M_\odot$).

4. Evolución con el Corrimiento al Rojo:

   • La evolución de la tasa de fusión es decididamente positiva, con $\kappa$ variando entre 0.4 y 4.7 con una credibilidad del 90%.
   • No hay evidencia fuerte de evolución dependiente de la masa; la evolución parece ser consistente a través de todo el espectro de masas, aunque la precisión está limitada por las incertidumbres en la distancia de luminosidad.

Significancia y Reivindicaciones
El artículo sostiene que la estructura observada en la población de BBH, particularmente los tres picos bien ubicados en la distribución de masa de chirp, proporciona evidencia convincente de que las fusiones jerárquicas dan forma a la distribución de masa.

• Los autores argumentan que la correlación única entre las masas de los componentes —que produce picos en la masa de chirp sin picos correspondientes en las masas individuales— es difícil de explicar si los picos surgieran de procesos no relacionados.
• El estudio sugiere que las fusiones jerárquicas pueden desempeñar un papel más dominante en la configuración de la distribución de masa de los agujeros negros binarios de lo que se pensaba anteriormente, extendiéndose incluso a las masas más bajas.
• La presencia de una subpoblación de alto espín cuyas masas de componentes se alinean con las predicciones jerárquicas ofrece evidencia directa de este escenario, incluso si la población general está dominada por binarias de bajo espín, probablemente de origen estelar.
• Los autores mantienen un tono modesto, reconociendo que la mayoría de las características siguen siendo tentativas y requieren la mayor estadística de futuros catálogos (como el resto de la cuarta campaña de observación) para establecer una confianza robusta. Enfatizan que, aunque las mediciones de espín son imprecisas, la distribución de masa por sí sola porta las firmas de los procesos jerárquicos.