The redshift evolution of the M_(rm BH)-M_* scaling relation: new insights from cosmological simulations and semi-analytic models

We study the co-evolution of black holes (BHs) and their host galaxies in the ASTRID and Illustris-TNG300 cosmological simulations and the Dark Sage Semi-Analytic Model (SAM), focusing on the evolution of the BH mass - stellar mass ($M_{\rm BH}-M_*$) relation. Due to differences in the adopted sub-grid modeling of BH seeding, dynamics, and feedback, the models differ in their predicted redshift evolution of the $M_{\rm BH}-M_*$ relation. We find that it is the interplay between the star formation rate (SFR) and the black hole accretion rate (BHAR) which drives the evolution of the mean relation. We define a quantity $\mathcal{R}$, the ratio between the specific BHAR and SFR (i.e. $\mathcal{R} \equiv\ $sBHAR/sSFR), and demonstrate that it is $\mathcal{R}$ that governs the evolution of individual sources in the $M_{\rm BH}-M_*$ plane. The efficiency of BH growth versus stellar mass growth in the sSFR-sBHAR plane reflects the partitioning of gas between fueling star formation versus BH accretion. This partitioning depends on the implementation of BH dynamics and the nature of how AGN feedback quenches galaxies. In the cosmological simulations (ASTRID and Illustris-TNG300), the BHAR and SFR are intrinsically linked, resulting in a tight $M_{\rm BH}-M_*$ correlation, while the Dark Sage SAM produces a significantly larger scatter. We discuss these results in the context of recently discovered over-massive BHs and massive quenched galaxies at high redshift by the James Webb Space Telescope.