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

Abstract

Most observations and clustering analyses suggest that quasars inhabit a narrow range of dark-matter halo masses ($10^{12-13}$ M$_{\odot}$) across cosmic time ($z\lesssim7$). Recent hydrodynamical simulations in gigaparsec-scale volumes now enable direct comparison of this picture with self-consistent galaxy-formation models. We quantify the relation between AGN bolometric luminosity and host halo mass before Cosmic Noon in TNG300, TNG-Cluster, FLAMINGO L1_m8 and L2p8_m9, and in smaller-volume simulations (Illustris, EAGLE, TNG100, and Simba). For AGNs with $L^{\mathrm{AGN}}_{\mathrm{bol}} \ge 10^{42}$ erg s$^{-1}$, more massive haloes host more luminous AGNs on average, but only up to a certain mass. The median luminosity-halo mass relation is highly non-linear, with large scatter, and flattens (FLAMINGO) or turns over (TNG300+TNG-Cluster) at halo mass, $M_{\mathrm{200,crit}} \gtrsim10^{12}$ M$_{\odot}$, at least at $z<5$-6. This high mass AGN quenching also manifests as a characteristic quasar host halo mass: in TNG300+TNG-Cluster, quasars ($L^{\mathrm{AGN}}_{\mathrm{bol}} \sim10^{45-47}$ erg s$^{-1}$) typically reside in haloes of mass $10^{12-12.5}$ M$_{\odot}$ at $z=3$-6, while FLAMINGO quasars extend to median masses of $\sim10^{12.8}$ M$_{\odot}$ at $z\sim3$-4. All simulations predict substantially larger scatter in AGN luminosity at fixed halo mass than in halo mass at fixed luminosity (up to 3 dex versus $\lesssim1$ dex between the 5th and 95th percentiles), implying weak coupling between halo growth and instantaneous SMBH accretion. Consequently, simulated quasar host masses broadly agree with observational estimates. The most luminous AGNs occupy increasingly rare haloes at earlier epochs but typically do not reside in the most massive haloes at any redshift up to $z\approx7$.