Implications of the nanoHertz Gravitational-Wave Background for Galactic Feedback and Massive Black Hole Growth

Résumé technique : Implications du fond d'ondes gravitationnelles nanoHertz pour le retour galactique et la croissance des trous noirs massifs

Énoncé du problème
Les mesures récentes effectuées par les réseaux de chronométrage de pulsars (PTA) ont fourni des preuves solides d'un fond stochastique d'ondes gravitationnelles (GWB) dans le domaine nanoHertz, probablement généré par la superposition incohérente de binaires de trous noirs supermassifs (SMBHB). L'amplitude de ce fond est sensible à la fonction de masse des trous noirs supermassifs (SMBH). Cependant, une tension existe : les simulations hydrodynamiques cosmologiques (telles que IllustrisTNG, MillenniumTNG et Simba) sous-estiment systématiquement l'amplitude du GWB déduite des données PTA. Cette divergence suggère que les modèles actuels de croissance des SMBH, de leur formation initiale (seeding) et des mécanismes de retour qui les régulent (retour AGN et retour stellaire) pourraient être incomplets. Plus précisément, il est incertain comment les variations dans les prescriptions de retour galactique modifient l'extrémité haute de la fonction de masse des trous noirs (BHMF) et, par conséquent, le signal GWB prédit.

Méthodologie
Les auteurs emploient un cadre semi-analytique basé sur Casey-Clyde et al. (2022) pour relier les BHMF prédits par les simulations à l'amplitude du GWB. La méthodologie se déroule comme suit :

1. Ensembles de simulations : L'étude analyse quatre grands ensembles de simulations : IllustrisTNG, MillenniumTNG, Simba et le projet CAMELS (Cosmology and Astrophysics with MachinE Learning Simulations).
2. Taille de la boîte et convergence : Pour répondre à la sensibilité du GWB à l'extrémité haute de la BHMF, les auteurs mènent des études de convergence en utilisant plusieurs tailles de boîtes (allant de 25 Mpc/h à 500 Mpc/h). Ils déterminent que les boîtes inférieures à ~35 Mpc/h (pour TNG) et ~50 Mpc/h (pour Simba) souffrent de bruit de Poisson et d'un échantillonnage insuffisant des halos massifs ($M_{BH} > 10^9 M_\odot$).
3. Variations du retour : L'étude explore l'impact du retour galactique en :
   • Comparant les modèles de référence (fiduciaux) entre différents codes de simulation.
   • Analysant les variations « 1P » (un paramètre) dans les ensembles CAMELS, où 23 paramètres astrophysiques contrôlant le retour stellaire et AGN sont variés individuellement.
   • Examinant des variantes spécifiques du retour dans Simba où les modes de jets AGN sont supprimés ou modifiés (par exemple, « SW only », « SW + Rad »).
4. Calcul du GWB : Les BHMF issus des simulations sont convertis en une population de binaires à l'aide d'une fonction de distribution du rapport d'Eddington (ERDF) pour identifier les quasars actifs. Cette population est mise à l'échelle par un facteur de normalisation ($N_{cal} \approx 4$) pour tenir compte du rapport binaire/quasar, et les masses secondaires sont attribuées via une distribution du rapport de masses. La population binaire résultante est intégrée pour calculer le spectre de déformation caractéristique et l'amplitude du GWB à une fréquence de référence de $f = 1 \text{ an}^{-1}$.

Contributions clés

• Quantification de l'impact du retour : L'article démontre que le choix des modèles de retour stellaire et AGN modifie significativement l'extrémité haute de la BHMF. Un retour AGN inefficace ou absent conduit à une surabondance de SMBH massifs, augmentant l'amplitude GWB prédite par des facteurs allant de 2 (dans les variations fiduciales) à 10 (dans les variations extrêmes de CAMELS).
• Analyse de convergence : Les auteurs établissent les tailles minimales de boîtes requises pour des prédictions GWB robustes, identifiant que les statistiques pour $M_{BH} > 10^9 M_\odot$ ne sont pas convergées dans les boîtes CAMELS standards de 25 Mpc/h.
• Sous-estimation systématique : L'étude confirme que les modèles de référence d'IllustrisTNG, MillenniumTNG et Simba sous-estiment l'amplitude du GWB déduite des PTA d'un facteur d'environ 2. Pour correspondre aux observations, ces modèles nécessiteraient 1,5 à 5,5 fois plus de SMBH dans la plage $M_{BH} > 10^8 M_\odot$.

Résultats

• Modèles de référence : Tous les modèles de simulation standards (fiduciaux) sous-estiment l'amplitude du GWB par rapport aux données NANOGrav de 15 ans. Cela implique que les prescriptions de retour actuelles suppriment trop efficacement la croissance des SMBH ou que la formation initiale (seeding) se produit trop tard.
• Variantes du retour : Les simulations avec un retour AGN réduit ou absent (notamment les exécutions Simba sans modes de jets) produisent des amplitudes GWB cohérentes avec les données PTA. Cependant, ces modèles échouent à reproduire des populations galactiques réalistes (par exemple, ils ne parviennent pas à éteindre efficacement la formation d'étoiles), les rendant physiquement irréalistes.
• Sensibilité aux paramètres : Dans les ensembles 1P de CAMELS, des paramètres spécifiques pilotent l'amplitude du GWB :
  • Pour CAMELS-TNG, le paramètre $\beta$ (contrôlant la transition entre les états AGN à faible et à fort taux d'accrétion) a l'impact le plus important. Une valeur de $\beta=0$ correspond aux données PTA mais sous-estime la fonction de masse stellaire à basse masse.
  • Pour CAMELS-Simba, le paramètre du rapport de masse de la graine ($A_{BHseed}$) est le plus critique. Abaisser le seuil de la graine permet une croissance plus précoce des SMBH, correspondant aux données PTA, mais pourrait perturber la formation d'étoiles précoce.
• Retour stellaire : Bien qu'on suppose souvent qu'il n'affecte que les trous noirs de faible masse, l'étude révèle que certains paramètres de retour stellaire (par exemple, la normalisation de la vitesse du vent et la fraction thermique) influencent également la BHMF haute masse et l'amplitude du GWB en régulant le carburant disponible pour la croissance des trous noirs.

Signification et affirmations
L'article affirme que les mesures PTA fournissent une sonde puissante et indépendante de la physique du retour galactique, complétant les contraintes électromagnétiques traditionnelles. Le décalage entre les prédictions des simulations et les données PTA suggère que la croissance des SMBH dans l'univers primordial est plus efficace ou se produit plus tôt que ce que capturent les modèles actuels.

Les auteurs concluent que, bien que la modification des paramètres de retour ou des prescriptions de formation initiale puisse aligner artificiellement les simulations sur les données PTA, de telles modifications rompent souvent l'accord avec les statistiques galactiques observées (telles que la fonction de masse stellaire). Par conséquent, la tension entre l'amplitude du GWB et les prédictions des simulations indique un besoin fondamental de réévaluer les prescriptions de formation initiale et de croissance des SMBH dans les simulations hydrodynamiques cosmologiques. L'étude souligne que les modèles futurs doivent simultanément reproduire des populations galactiques réalistes et l'amplitude observée du GWB, suggérant que les modèles de retour actuels pourraient supprimer excessivement la croissance des SMBH. Les auteurs notent que leur cadre repose sur des hypothèses simplificatrices (par exemple, négligeant le durcissement environnemental et les détails temporels des fusions) et qu'un traitement auto-cohérent de l'évolution binaire est requis pour la précision future.