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

Énoncé du problème
Avec la publication du quatrième catalogue de l'onde gravitationnelle LIGO–Virgo–KAGRA (GWTC-4.0), le nombre de fusions de trous noirs binaires (BBH) détectées avec certitude a dépassé les 150. Alors que les premières attentes suggéraient qu'un petit nombre de détections pourrait suffire à identifier un canal de formation primaire unique, les données révèlent désormais une population complexe qui ne peut être expliquée facilement par un seul mécanisme. Les tentatives précédentes visant à identifier des sous-populations à l'aide des inclinaisons de spin (spin tilts) ont été confrontées à des difficultés dues aux difficultés de mesure et aux incertitudes théoriques. De plus, bien que des corrélations globales entre les paramètres existent, leur interprétation astrophysique reste non triviale. Le problème central abordé est de savoir si des sous-populations distinctes de BBH en fusion existent, caractérisées par des frontières de transition claires dans la masse primaire ($m_1$) et des distributions divergentes des rapports de masse ($q$) et des grandeurs de spin ($\chi$).

Méthodologie
Les auteurs adoptent une approche phénoménologique pour modéliser la population de BBH en utilisant les données cumulatives de GWTC-4.0 (en excluant GW190814 en raison de l'ambiguïté de la composante secondaire et GW231123 en raison de préoccupations concernant la systématique des formes d'ondes et les spins extrêmes). Ils emploient l'inférence bayésienne hiérarchique pour comparer des modèles présentant différents niveaux de complexité.

Caractéristiques méthodologiques clés :

• Paramétrage : La distribution de masse primaire ($p(m_1)$) est modélisée comme une loi de puissance brisée avec deux pics gaussiens. Les grandeurs de spin ($p(\chi)$) et les inclinaisons ($p(\cos \theta)$) sont modélisées à l'aide de normales tronquées et de mélanges de composantes isotropes/alignées, respectivement. Le redshift ($p(z)$) suit une loi de puissance.
• Masses de transition : L'innovation centrale consiste à permettre aux distributions de $q$, $\chi$ et $\cos \theta$ de changer de forme à jusqu'à trois masses de transition, créant ainsi jusqu'à quatre sous-populations. La forme fonctionnelle de la distribution peut rester la même à travers une transition, mais les hyperparamètres contrôlant la forme sont autorisés à différer.
• Sélection de modèle : Les auteurs utilisent les facteurs de Bayes pour effectuer la sélection de modèles, en appliquant une pénalité d'Occam pour éviter le surapprentissage (overfitting). Ils comparent un modèle « multipop » (autorisant les transitions) à un modèle fiduciel fortement paramétré (gwtc4-fiducial) et à des variantes où des caractéristiques spécifiques sont supprimées.

Contributions clés et résultats
L'analyse fournit des preuves solides de l'existence d'au moins trois sous-populations distinctes de fusions de trous noirs binaires, séparées par deux masses de transition : $m_t^{AB} = 27,7^{+4,1}_{-3,4} M_\odot$ et $m_t^{BC} = 40,2^{+4,7}_{-3,2} M_\odot$. Le modèle « multipop » est préféré au modèle fiduciel avec un logarithme naturel du facteur de Bayes de 12,3.

Les trois sous-populations sont caractérisées comme suit :

1. Sous-population A ($m_1 \lesssim 27,7 M_\odot$) :

   • Rapport de masse : Présente une distribution presque plate ($\beta_0 \approx -0,8$).
   • Spin : Caractérisée par de faibles grandeurs de spin, avec un support principalement pour $\chi \leq 0,5^{+0,1}_{-0,1}$.
   • Interprétation : Inclut probablement des contributions de multiples canaux de formation, potentiellement incluant des binaires isolées avec un transport efficace du moment angulaire.

2. Sous-population B ($27,7 M_\odot \lesssim m_1 \lesssim 40,2 M_\odot$) :

   • Rapport de masse : Montre une préférence beaucoup plus marquée pour l'appariement de masses égales ($q \approx 1$), avec $\beta_1 \approx 3,6$. Cette distribution d'appariement distincte résout un pic précédemment rapporté dans la distribution de masse primaire à $\sim 35 M_\odot$, suggérant que ce pic était un artefact de spécification erronée du modèle.
   • Spin : Les distributions de grandeur de spin et d'inclinaison sont statistiquement indiscernables de la sous-population A.
   • Interprétation : Les auteurs suggèrent l'évolution homogène chimiquement (CHE) ou les étoiles de Population III comme origines potentielles, bien qu'ils notent des difficultés à réconcilier l'absence de différences de spin avec les prédictions standard de la CHE. La formation dynamique est considérée comme moins probable compte tenu de l'absence de distributions de tilt de spin distinctes.

3. Sous-population C ($m_1 \gtrsim 40,2 M_\odot$) :

   • Spin : Caractérisée par une distribution de grandeur de spin nettement plus large avec un support substantiel pour $\chi > 0,5$ (culminant autour de $\chi \sim 0,7$).
   • Rapport de masse : Montre un support ténu pour une préférence vers $q \approx 0,5$.
   • Interprétation : Les propriétés sont largement cohérentes avec les fusions hiérarchiques (fusions de trous noirs de générations précédentes), particulièrement la large distribution de spin et la préférence pour $q \approx 0,5$. Cependant, la présence de petits spins au sein de ce groupe et l'absence de transition dans les inclinaisons de spin suggèrent que les fusions hiérarchiques dans les amas denses n'expliquent pas entièrement cette sous-population.

Signification et affirmations
L'article affirme que les données préfèrent fortement l'interprétation selon laquelle ces trois sous-populations existent conjointement, plutôt que comme une population continue unique. L'identification de frontières de transition nettes dans les distributions de rapport de masse et de grandeur de spin implique que l'hypothèse de séparabilité dans la distribution de la population des paramètres de BBH se brise au-dessus de $\sim 25 M_\odot$.

Les auteurs soulignent que bien que leurs résultats reposent sur des modèles fortement paramétrés et soient sujets à une potentielle erreur de spécification de modèle, la force statistique du signal et la confirmation indépendante de ces masses de transition dans d'autres analyses de la littérature suggèrent que l'existence de ces sous-populations est réelle. Ils concluent que différents chemins de formation dominent probablement différentes plages de masse, la sous-population B étant potentiellement liée à la CHE ou aux étoiles de Population III, et la sous-population C liée aux fusions hiérarchiques, bien que des scénarios alternatifs (tels que l'accrétion dans les disques d'AGN) restent possibles. L'étude stipule explicitement qu'elle ne trouve aucune preuve de transitions autour du pic de $10 M_\odot$, laissant cette sous-structure pour des investigations futures.