Forbidden Formation Histories: The Binary Black Hole Merger Rate Disfavors Long Delay Times

Résumé technique : Histories de formation interdites : Le taux de fusion des trous noirs binaires défavorise les temps de retard longs

Énoncé du problème
Le taux de fusion des trous noirs binaires (BBH), $R_{\text{merge}}(t)$, est la convolution du taux de formation des progéniteurs, $R_{\text{form}}(t)$, et de la distribution des temps de retard entre la formation et la fusion, $p(\tau)$. Les approches traditionnelles visant à comprendre les origines des BBH reposent souvent sur une modélisation directe : on suppose des formes fonctionnelles spécifiques pour $R_{\text{form}}$ (par exemple, liées au taux de formation stellaire global) et pour $p(\tau)$ (dérivées de la synthèse des populations), puis on compare le taux de fusion prédit aux données d'ondes gravitationnelles (GW). Cette approche souffre de dégénérescences inhérentes ; différentes combinaisons d'histoires de formation et de temps de retard peuvent produire des taux de fusion identiques, et des hypothèses paramétriques rigides peuvent obscurcir des caractéristiques des données ou forcer des résultats physiquement incohérents. De plus, des modèles forts pour $R_{\text{form}}$ et $p(\tau)$ peuvent fondamentalement ne pas correspondre à l'évolution en décalage vers le rouge observée du taux de fusion.

Méthodologie
Les auteurs proposent un cadre de « modélisation inverse » pour déconvoluer le taux de fusion observé des BBH en ses composantes constitutives sans supposer de forme fonctionnelle spécifique pour le taux de formation des progéniteurs.

1. Cadre mathématique : En utilisant le théorème de convolution, les auteurs transforment l'équation intégrale $R_{\text{merge}}(t) = \int R_{\text{form}}(t-\tau)p(\tau)d\tau$ dans le domaine fréquentiel via la transformée de Fourier rapide (FFT). Cela permet la reconstruction directe de $R_{\text{form}}$ étant donné un $R_{\text{merge}}$ observé et un $p(\tau)$ connu :
   $$R_{\text{form}} = \mathcal{F}^{-1} \left\{ \frac{\mathcal{F}\{R_{\text{merge}}\}}{\mathcal{F}\{p(\tau)\}} \right\}$$
2. Implémentation numérique : Pour gérer le caractère non périodique des données et l'horizon d'observation fini (données GWTC-4.0 jusqu'à $z \approx 1.5$), les auteurs utilisent un remplissage par zéro pour éviter les artefacts de repliement et appliquent un écrêtage en cosinus pour atténuer le bruit du phénomène de Gibbs aux frontières. Ils supposent une condition aux limites conservative où $R_{\text{merge}}(z > 1.5) = 0$, bien qu'ils testent des valeurs asymptotiques jusqu'à 500 Gpc$^{-3}$yr$^{-1}$.
3. Régularisation : Pour assurer la stabilité numérique lors de la division par la transformée de Fourier de la distribution des temps de retard (qui peut s'approcher de zéro), ils utilisent la déconvolution de Wiener, en ajoutant un terme de régularisation $\lambda$ au dénominateur.
4. Données d'entrée : La méthode est appliquée en utilisant l'ajustement non paramétrique par B-splines du taux de fusion des BBH issu de GWTC-4.0 et une série de distributions de temps de retard générées par le code de synthèse des populations COMPAS, qui font varier des paramètres liés au transfert de masse, aux supernovae et aux masses des résidus.

Contributions et résultats clés

• Reconstruction des histoires de formation : En déconvoluant les taux de fusion de GWTC-4.0 avec les distributions de temps de retard de COMPAS, les auteurs reconstruisent le taux de formation des progéniteurs implicite. Ils constatent que le $R_{\text{form}}$ déduit décroît plus rapidement vers les faibles décalages vers le rouge que le taux de formation stellaire global (SFR), ce qui est cohérent avec une efficacité de formation des BBH dépendante de la métallicité.
• Identification d'« histoires interdites » : L'étude démontre que certaines combinaisons de taux de fusion et de distributions de temps de retard sont physiquement impossibles car elles nécessitent que $R_{\text{form}}$ devienne négatif à des époques tardives. Cela se produit lorsqu'une distribution de temps de retard avec une longue queue (indice de loi de puissance faible) prédit plus de fusions issues d'une formation ancienne à haut décalage vers le rouge que ce qui est effectivement observé à faible décalage vers le rouge.
• Contraintes sur les distributions de temps de retard : En utilisant l'exigence que $R_{\text{form}} > 0$, les auteurs contraignent l'indice de loi de puissance de la distribution des temps de retard $\alpha$ (où $p(\tau) \propto \tau^\alpha$).
  • Les distributions de temps de retard avec des indices faibles ($\alpha \gtrsim -0.7$), caractéristiques des canaux de transfert de masse stable, sont fortement défavorisées.
  • Plus précisément, pour les données GWTC-4.0, $\alpha < -1.08$ est préféré à $\alpha > -0.67$ par un facteur de 19.
  • Les simulations COMPAS avec des distributions raides (favorisant des retards courts, $\alpha < -1.2$) sont cohérentes avec presque tous les tirages a posteriori du taux de fusion de GWTC-4.0, tandis que les distributions faibles ne sont compatibles qu'avec environ 30 % des tirages.
• Contraintes sur la physique de l'évolution binaire : À titre de preuve de concept, les auteurs contraignent le paramètre d'efficacité du transfert de masse $\beta$ dans COMPAS. Des valeurs plus faibles de $\beta$ (transfert de masse non conservateur, conduisant à un rétrécissement orbital et à des retards courts) montrent une compatibilité plus élevée avec les données GW. Des valeurs élevées de $\beta$ (transfert de masse conservateur, retards longs) sont défavorisées par un facteur d'environ 3,7.
• Impact sur l'inférence du taux de fusion : Imposer la contrainte physique d'un taux de formation positif comme a priori sur le modèle B-spline de GWTC-4.0 réduit le taux de fusion médian déduit à $z=1.5$ d'environ 10 % (de 334 à 299 Gpc$^{-3}$yr$^{-1}$), favorisant une évolution du taux de fusion plus faible que l'inférence standard basée sur les données.
• Cartographie de l'espace des paramètres : Les auteurs cartographient la région interdite dans l'espace des paramètres $\alpha$–$\kappa$ (où $\kappa$ décrit l'évolution du taux de fusion $R_{\text{merge}} \propto (1+z)^\kappa$). Ils montrent que pour un $\kappa$ donné, il existe un $\alpha_{\text{crit}}$ critique au-dessus duquel le taux de formation devient non physique. Les données actuelles excluent les distributions de temps de retard en loi de puissance avec $\alpha \gtrsim -0.7$.

Signification
L'article présente un cadre indépendant du modèle pour exclure des régions de l'évolution binaire et de l'espace des paramètres du taux de fusion. En évaluant directement la compatibilité des paramètres de synthèse des populations avec les observations d'ondes gravitationnelles à travers le prisme de la causalité (c'est-à-dire que le taux de formation ne peut pas être négatif), la méthode contourne le besoin d'a priori de haute dimension sur la formation stellaire et l'évolution de la métallicité. Les auteurs soutiennent que cette approche expose des « histoires » de formation interdites que la modélisation directe pourrait manquer, fournissant une vérification de cohérence robuste pour les théories de l'évolution des binaires stellaires. Les résultats suggèrent que les canaux de formation prédisant des temps de retard longs (tels que le transfert de masse stable) sont en tension avec la forte décroissance observée du taux de fusion des BBH à faible décalage vers le rouge, favorisant les canaux produisant des retards plus courts.