Evidence for Three Subpopulations of Merging Binary Black Holes at Different Primary Masses

Planteamiento del Problema
Con el lanzamiento del cuarto catálogo de ondas gravitacionales LIGO–Virgo–KAGRA (GWTC-4.0), el número de fusiones de agujeros negros binarios (BBH) detectadas con confianza ha superado los 150. Aunque las expectativas iniciales sugerían que un pequeño número de detecciones podría bastar para identificar un único canal de formación primario, los datos revelan ahora una población compleja que no puede explicarse fácilmente mediante un solo mecanismo. Los intentos previos de identificar subpoblaciones mediante las inclinaciones de espín (spin tilts) han enfrentado desafíos debido a las dificultades de medición y las incertidumbres teóricas. Además, aunque existen correlaciones globales entre los parámetros, su interpretación astrofísica sigue siendo no trivial. El problema central abordado es si existen subpoblaciones distintas de BBH en fusión, caracterizadas por límites de transición claros en la masa primaria ($m_1$) y distribuciones diferenciadas de las razones de masa ($q$) y las magnitudes de espín ($\chi$).

Metodología
Los autores adoptan un enfoque fenomenológico para modelar la población de BBH utilizando los datos acumulados de GWTC-4.0 (excluyendo GW190814 debido a la ambigüedad de la secundaria y GW231123 debido a preocupaciones sobre la sistemática de la forma de onda y espines extremos). Emplean la inferencia bayesiana jerárquica para comparar modelos con distintos niveles de complejidad.

Características metodológicas clave incluyen:

• Parametrización: La distribución de masa primaria ($p(m_1)$) se modela como una ley de potencia rota con dos picos gaussianos. Las magnitudes de espín ($p(\chi)$) y las inclinaciones ($p(\cos \theta)$) se modelan mediante normales truncadas y mezclas de componentes isotrópicos/alineados, respectivamente. El desplazamiento al rojo ($p(z)$) sigue una ley de potencia.
• Masas de transición: La innovación central es permitir que las distribuciones de $q$, $\chi$ y $\cos \theta$ cambien de forma en hasta tres masas de transición, creando efectivamente hasta cuatro subpoblaciones. La forma funcional de la distribución puede permanecer igual a través de una transición, pero los hiperparámetros que controlan la forma pueden diferir.
• Selección de modelos: Los autores utilizan factores de Bayes para realizar la selección de modelos, aplicando una penalización de Occam para evitar el sobreajuste (overfitting). Comparan un modelo de "múltiples poblaciones" (que permite transiciones) contra un modelo fiduciario fuertemente parametrizado (gwtc4-fiducial) y variantes donde se eliminan características específicas.

Contribuciones Clave y Resultados
El análisis proporciona evidencia sólida de al menos tres subpoblaciones distintas de fusiones de agujeros negros binarios, separadas por dos masas de transición: $m_t^{AB} = 27.7^{+4.1}_{-3.4} M_\odot$ y $m_t^{BC} = 40.2^{+4.7}_{-3.2} M_\odot$. El modelo de "múltiples poblaciones" es preferido sobre el modelo fiduciario con un factor de Bayes de logaritmo natural de 12.3.

Las tres subpoblaciones se caracterizan de la siguiente manera:

1. Subpoblación A ($m_1 \lesssim 27.7 M_\odot$):

   • Razón de masa: Exhibe una distribución casi plana ($\beta_0 \approx -0.8$).
   • Espín: Se caracteriza por magnitudes de espín pequeñas, con soporte principalmente para $\chi \leq 0.5^{+0.1}_{-0.1}$.
   • Interpretación: Probablemente incluye contribuciones de múltiples canales de formación, incluyendo potencialmente binarias aisladas con un transporte de momento angular eficiente.

2. Subpoblación B ($27.7 M_\odot \lesssim m_1 \lesssim 40.2 M_\odot$):

   • Razón de masa: Muestra una preferencia mucho más marcada por el emparejamiento de masas iguales ($q \approx 1$), con $\beta_1 \approx 3.6$. Este distintivo emparejamiento de distribución resuelve un pico reportado previamente en la distribución de masa primaria en $\sim 35 M_\odot$, sugiriendo que dicho pico fue un artefacto de la especificación errónea del modelo.
   • Espín: Las distribuciones de magnitud y de inclinación de espín son estadísticamente indistinguibles de la Subpoblación A.
   • Interpretación: Los autores sugieren la Evolución Homogénea Químicamente (CHE) o estrellas de Población III como posibles orígenes, aunque señalan desafíos para reconciliar la falta de diferencias de espín con las predicciones estándar de la CHE. La formación dinámica se considera menos probable dada la falta de distribuciones de inclinación de espín distintivas.

3. Subpoblación C ($m_1 \gtrsim 40.2 M_\odot$):

   • Espín: Se caracteriza por una distribución de magnitud de espín significativamente más amplia con un soporte sustancial para $\chi > 0.5$ (con un pico alrededor de $\chi \sim 0.7$).
   • Razón de masa: Muestra un apoyo tentativo hacia una preferencia por $q \approx 0.5$.
   • Interpretación: Las propiedades son mayormente consistentes con fusiones jerárquicas (fusiones de agujeros negros de generaciones anteriores), particularmente la amplia distribución de espín y la preferencia por $q \approx 0.5$. Sin embargo, la presencia de espines pequeños dentro de este grupo y la falta de una transición en las inclinaciones de espín sugieren que las fusiones jerárquicas en cúmulos densos podrían no explicar completamente esta subpoblación.

Significancia y Reivindicaciones
El artículo afirma que los datos prefieren fuertemente la interpretación de que estas tres subpoblaciones existen en conjunto, en lugar de como una única población continua. La identificación de límites de transición agudos en las distribuciones de razón de masa y magnitud de espín implica que el supuesto de separabilidad en la distribución de los parámetros de la población de BBH se rompe por encima de $\sim 25 M_\odot$.

Los autores enfatizan que, aunque sus resultados dependen de modelos fuertemente parametrizados y están sujetos a una posible especificación errónea del modelo, la fuerza estadística de la señal y la confirmación independiente de estas masas de transición en otros análisis de la literatura sugieren que la existencia de estas subpoblaciones es real. Concluyen que diferentes vías de formación probablemente dominan diferentes rangos de masa, con la Subpoblación B potencialmente vinculada a la CHE o a estrellas de Población III, y la Subpoblación C vinculada a fusiones jerárquicas, aunque otros escenarios (como la acreción en discos de AGN) siguen siendo posibles. El estudio establece explícitamente que no encuentra evidencia de transiciones alrededor del pico de $10 M_\odot$, dejando esa subestructura para investigaciones futuras.