The host halo masses of AGNs and quasars at $z \sim 3-7$ with TNG-Cluster, FLAMINGO and other cosmological galaxy simulations

En utilisant des simulations cosmologiques à grande échelle, cette étude révèle que bien que la luminosité des AGN augmente généralement avec la masse du halo hôte jusqu'à un certain seuil, la relation est hautement non linéaire avec une dispersion significative, indiquant que les quasars les plus lumineux à $z \sim 3-7$ habitent typiquement des halos de masse intermédiaire ($10^{12-12,5}$ M$_{\odot}$) plutôt que les plus massifs, une découverte qui s'aligne bien avec les estimations observationnelles.

L'essentiel

La vue d'ensemble : Trouver le « point idéal » pour les superstars cosmiques

Imaginez l'univers comme une ville géante en pleine croissance. Dans cette ville, les « halos de matière noire » sont les quartiers, et les « noyaux galactiques actifs » (AGN) ou « quasars » sont les phares massifs et aveuglants situés au centre des bâtiments les plus importants. Ces phares sont alimentés par des trous noirs supermassifs dévorant du gaz.

Pendant longtemps, les astronomes ont tenté de répondre à une question simple : Quelle taille de quartier (halo) est nécessaire pour abriter le phare le plus brillant (quasar) ?

La plupart des observations suggèrent que ces quasars ultra-brillants vivent dans des quartiers de taille sensiblement la même, peu importe l'époque (ou la distance) à laquelle on regarde en arrière. Ils semblent préférer une zone « Goldilocks » (ni trop chaud, ni trop froid) : ni trop petit, ni trop grand.

Ce document utilise des simulations informatiques massives pour tester cette idée. Les chercheurs ont construit des univers numériques (en utilisant des modèles comme TNG-Cluster et FLAMINGO) pour voir si leurs quasars virtuels se comportent comme les vrais.

La découverte principale : Ce n'est pas une ligne droite

Les chercheurs ont découvert que la relation entre la taille du quartier et la luminosité du phare n'est pas une ligne droite. C'est plutôt comme une colline avec un sommet.

1. L'ascension : À mesure que le quartier s'agrandit (devient plus massif), le phare devient généralement plus brillant. Cela est logique ; les quartiers plus grands disposent de plus de gaz pour nourrir le trou noir.
2. Le sommet : Cette tendance ne progresse que jusqu'à un certain point (environ $10^{12}$ fois la masse de notre Soleil).
3. La chute : Une fois que le quartier devient trop massif, le phare devient en réalité plus faible ou reste stable.

L'analogie : Pensez-y comme à une fête.

• Dans une petite maison (petit halo), il y a peu de monde, et la musique est calme.
• Dans une maison de taille moyenne (le point idéal), vous avez la foule parfaite, et la fête est bruyante et énergique.
• Dans un stade massif (un énorme halo), la fête devient en fait plus calme. Pourquoi ? Parce que l'« hôte » (le trou noir) devient si gourmand et puissant qu'il commence à expulser les invités de la pièce (c'est ce qu'on appelle la rétroaction des AGN). Il repousse le gaz, s'affamant lui-même et empêchant la fête de devenir plus bruyante.

Le facteur « Chaos » : Une dispersion énorme

L'une des découvertes les plus surprenantes est à quel point les données sont désordonnées. Le document compare deux façons d'observer les données :

• Scénario A : Si vous choisissez une taille de quartier spécifique, quelle est la luminosité du phare ?
  • Résultat : Totalement imprévisible. La luminosité peut varier d'un facteur 1 000 à 10 000 (3 à 4 « décades » de différence). Une maison de même taille pourrait avoir une petite veilleuse, tandis qu'une autre possède un projecteur aveuglant.
• Scénario B : Si vous choisissez une luminosité spécifique (par exemple, un quasar ultra-brillant), quelle est la taille du quartier ?
  • Résultat : Beaucoup plus prévisible. Si vous voyez un phare super brillant, il se trouve presque certainement dans un quartier d'une gamme de taille spécifique.

L'analogie : Imaginez essayer de deviner le poids d'une personne rien qu'en regardant sa taille.

• Si vous choisissez une taille spécifique (par exemple, 1m80), les gens peuvent peser n'importe quoi, de 60 kg à 120 kg. C'est une énorme dispersion.
• Mais si vous choisissez un poids spécifique (par exemple, 100 kg), la personne fait presque certainement environ 1m80. Le poids prédit mieux la taille que la taille ne prédit le poids.

Le document conclut que la taille du quartier est un mauvais prédicteur de la luminosité actuelle du trou noir. La luminosité du trou noir est chaotique et dépend de beaucoup d'autres choses (comme la quantité de gaz actuellement disponible), et pas seulement de la taille du quartier.

Le « Point Idéal » confirmé

Malgré le chaos, les simulations ont confirmé l'idée observationnelle : les quasars préfèrent une taille de quartier spécifique.

• Que l'univers ait 3 milliards d'années ou 7 milliards d'années, les quasars les plus brillants vivent généralement dans des halos ayant une masse d'environ $10^{12}$ à $10^{12,5}$ masses solaires.
• Ils vivent rarement dans les halos les plus massifs (les super-villes).
• Ils vivent rarement dans les plus petits halos (les petits villages).

Cela suggère qu'il existe un « point idéal » pour la croissance des trous noirs. Si le quartier est trop petit, le trou noir ne peut pas obtenir assez de nourriture. S'il est trop grand, la propre rétroaction du trou noir (en expulsant le gaz) l'empêche de croître.

Pourquoi cela importe

Le document soutient que nous ne pouvons pas simplement regarder la taille du quartier d'une galaxie pour deviner la luminosité de son trou noir. La connexion est trop lâche et trop chaotique. Cependant, si nous voyons un quasar super brillant, nous pouvons être assez certains de la taille de son foyer.

Les simulations montrent également que ce « point idéal » reste sensiblement de la même taille tout au long de l'histoire cosmique, même si l'univers est en expansion et en mutation. Cela implique que les règles régissant la façon dont les trous noirs se nourrissent et la façon dont ils s'arrêtent de le faire sont cohérentes à travers des milliards d'années.

Résumé en une phrase

Le document conclut que, bien que les trous noirs supermassifs vivent dans un type de quartier cosmique de taille « Goldilocks » spécifique, la luminosité du trou noir est si chaotique que connaître la taille du quartier ne vous apprend presque rien sur la luminosité du trou noir à un instant donné.

Résumé technique

Énoncé du Problème
Les études observationnelles et les analyses de regroupement ont longtemps suggéré que les quasars habitent une plage étroite de masses de halos de matière noire ($\sim 10^{12} - 10^{13} \, M_\odot$) à travers les époques cosmiques ($z \lesssim 7$). Cependant, ces inférences sont souvent indirectes, reposant sur des échantillons limités par le flux et des hypothèses concernant la monotonicité de la relation entre la luminosité des Noyaux Actifs de Galaxies (AGN) et la masse du halo hôte. Une lacune critique existe dans la compréhension de l'écart intrinsèque et de la nature précise de cette relation avant le "Midi Cosmique" ($z \gtrsim 3$), particulièrement concernant la question de savoir si les AGN les plus lumineux résident dans les halos les plus massifs et comment les mécanismes de rétroaction régulent cette connexion. Les simulations précédentes étaient souvent limitées par le volume (incapables d'échantillonner les halos massifs et rares) ou par la résolution (incapables de résoudre la dynamique des gaz à petite échelle pilotant les luminosités des AGN).

Méthodologie
Ce travail utilise une suite de simulations hydrodynamiques cosmologiques pour comparer directement les images observationnelles avec des modèles de formation de galaxies auto-cohérents. L'étude se concentre sur les AGN ayant des luminosités bolométriques $L_{\text{AGN}}^{\text{bol}} \geq 10^{42} \, \text{erg s}^{-1}$ aux redshifts $z = 3 - 7$.

L'analyse principale repose sur deux suites de simulations à grand volume :

1. TNG300 + TNG-Cluster : Partie du projet IllustrisTNG, offrant une haute résolution (masse cible baryonique $\sim 1.1 \times 10^7 \, M_\odot$) sur un volume de $300^3 \, \text{cMpc}^3$ et une extension de type "zoom-in" ciblant des halos massifs jusqu'à $10^{15} \, M_\odot$.
2. FLAMINGO : Englobant des volumes de $1 - 3 \, \text{Gpc}^3$ par côté (runs L1_m8 et L2p8_m9), offrant un échantillonnage statistique supérieur de la partie haute de la fonction de masse des halos, bien qu'à une résolution plus faible que TNG.

Pour l'étalonnage, les auteurs incluent des simulations de plus petits volumes ($\sim 100 \, \text{cMpc}^3$) issues de générations précédentes : Illustris, EAGLE, TNG100 et Simba.

Les luminosités bolométriques des AGN ne sont pas des sorties directes des simulations, mais sont calculées en post-traitement à partir des taux d'accrétion instantanés des trous noirs supermassifs (SMBH) en utilisant un modèle de conversion bimodal (Churazov et al. 2005) qui distingue les régimes d'accrétion radiativement efficaces et inefficaces. Une efficacité radiative uniforme de $\epsilon_r = 0.1$ est adoptée dans tous les modèles afin d'isoler les différences de la physique d'accrétion sous-jacente.

Contributions Clés et Résultats

• Relations Non-Linéaires et Non-Monotones : À travers tous les modèles, la relation médiane entre la luminosité bolométrique des AGN et la masse du halo hôte est hautement non linéaire. Bien que les halos plus massifs aient tendance à héberger des AGN plus lumineux en moyenne, cette tendance ne persiste que jusqu'à une masse de halo caractéristique de $M_{200,\text{crit}} \sim 10^{12} \, M_\odot$.

  • Dans TNG300+TNG-Cluster, la relation médiane s'aplatit et même s'inverse (chute) à $M_{200,\text{crit}} \gtrsim 10^{12} \, M_\odot$ pour $z < 5-6$.
  • Dans FLAMINGO, la relation s'aplatit à la limite de haute masse mais ne montre pas le même basculement prononcé.
  • Cette extinction à haute masse est attribuée à la rétroaction des AGN (cinétique dans TNG, thermique dans FLAMINGO) qui supprime l'accrétion de gaz sur les SMBH dans les halos massifs.

• Grand Écart Intrinsèque : Le papier quantifie une variation objet à objet significative. L'écart de la luminosité bolométrique des AGN à une masse de halo donnée est beaucoup plus grand (jusqu'à 3 dex pour les percentiles 5e–95e) que l'écart de la masse du halo hôte à une luminosité d'AGN donnée ($\lesssim 1$ dex). Cela implique un couplage faible entre la croissance du halo et l'accrétion instantanée du SMBH, suggérant que la masse du halo seule est un mauvais prédicteur de la luminosité instantanée de l'AGN.

• Masses des Halos Hôtes des Quasars :

  • TNG300+TNG-Cluster : Les quasars ($L_{\text{AGN}}^{\text{bol}} \sim 10^{45} - 10^{47} \, \text{erg s}^{-1}$) résident typiquement dans des halos de $10^{12} - 10^{12.5} \, M_\odot$ à $z=3-6$.
  • FLAMINGO : Les quasars s'étendent à des masses médianes de $\sim 10^{12.8} \, M_\odot$ à $z \sim 3-4$.
  • Crucialement, les AGN les plus lumineux n'habitent pas les halos les plus massifs disponibles à chaque époque donnée jusqu'à $z \approx 7$. Au lieu de cela, ils résident dans des halos progressivement plus rares aux époques précoces, mais évitent typiquement les systèmes de masse très élevée ($M_{200,\text{crit}} \gtrsim 10^{13} \, M_\odot$) en raison de la rétroaction d'extinction.

• Comparaison avec les Observations : Les masses de halos hôtes simulées pour les quasars correspondent largement aux estimations observationnelles dérivées des mesures de regroupement et de corrélation croisée, qui infèrent généralement des masses caractéristiques de $10^{12} - 10^{13} \, M_\odot$. Le grand écart prédit par les simulations aide à réconcilier les larges contraintes observationnelles.

Signification et Revendications
Le papier affirme que les simulations cosmologiques modernes, malgré leurs différences de résolution, de volume et de physique de sous-grille (amorçage, accrétion, rétroaction), convergent vers une image cohérente : la relation luminosité d'AGN–masse de halo est caractérisée par un "point idéal" de masse de halo caractéristique ($\sim 10^{12} \, M_\odot$) où les quasars sont les plus prévalents, flanqué par une extinction dans les halos plus massifs et une suppression dans les moins massifs.

Les auteurs soulignent que le grand écart intrinsèque de la luminosité des AGN à masse de halo fixe remet en question les hypothèses de nombreux modèles semi-empiriques d'occupation de halo, qui supposent souvent une relation plus étroite. Les résultats suggèrent que si le quasar typique réside dans une gamme de masse spécifique, les systèmes individuels présentent une immense diversité, et que les halos les plus massifs ne sont pas nécessairement les hôtes des quasars les plus lumineux à un instant donné. Ce travail fournit une base quantitative pour interpréter les observations de quasars à haut redshift et met en évidence le rôle de la rétroaction des AGN dans le façonnement de la démographie des trous noirs les plus lumineux avant le Midi Cosmique.