Abstract The coevolution of supermassive black holes (SMBHs) and their host galaxies remains one of the central open questions in cosmology, rooted in the coupling between accretion, feedback, and the multiscale physics that links the event horizon to the circumgalactic medium. Here we bridge these scales by embedding a first-principles, GRMHD-informed prescription for black hole accretion and feedback—derived from multizone simulations that self-consistently connect inflows and outflows from the horizon to the Bondi radius—within cosmological magnetohydrodynamic zoom-in simulations of ∼10 14 M ⊙ halos. These GRMHD results predict a “suppressed Bondi” regime in which magnetic stresses and relativistic winds strongly reduce effective accretion rates in a spin-dependent manner. We find that black holes cannot grow efficiently by accretion until they exceed ∼10 7 M ⊙ , regardless of the feedback strength. Beyond this threshold, systems bifurcate: low-spin ( η ∼ 0.02) black holes continue to accrete without quenching star formation, while high-spin ( η ≳ 0.3) black holes quench effectively but become starved of further growth. Early, massive seeding partially alleviates this tension through merger-driven assembly, yet an additional cold or super-Eddington accretion mode appears essential to reproduce the observed SMBH population and the empirical black hole–galaxy scaling relations. Our results demonstrate that GRMHD-informed feedback models can account for the maintenance-mode behavior of low-luminosity active galactic nuclei like M87*, but cannot by themselves explain the full buildup of SMBH mass across cosmic time. A unified, multiregime framework is required to capture the evolving interplay between spin-dependent feedback, cold inflows, and mergers in driving coevolution.
Su et al. (2026) studied this question.