2025MNRAS.536.3554X

Observations suggest a co-evolution of supermassive black holes (SMBHs) and host galaxies. In this paper, we consider the mass and energy flow in a near-equilibrium bulge suffused by gases of varying temperatures. By assuming the rate of energy flow independent of the distance r from the bulge center and the local virial equilibrium for permeated gases on scale r, a key parameter $\varepsilon _b$ was identified that quantifies the mass and energy flow in gases and the efficiency of gas cooling (or the "specific" cooling rate per unit mass), and thus regulates the co-evolution of both SMBHs and hosts. With the help of Illustris simulations and observations, we determined the redshift variation $\varepsilon _b\propto (1+z)^{5/2}$. A higher $\varepsilon _b$ in the early Universe means a higher specific cooling rate that allows rapid evolution of SMBHs and hosts. This simple theory, characterized by a single parameter $\varepsilon _b$, provides the dominant mean cosmic evolution of SMBHs and hosts. All other transient phenomena may only contribute to the dispersion around this mean evolution. Based on this theory and relevant assumptions, scaling laws involving $\varepsilon _b$ were identified for the evolution of SMBHs and hosts. For host galaxies, the mass-size relation $M_b\propto \varepsilon _b^{2/3}r_b^{5/3}G^{-1}$, the dispersion-size relation $\sigma _b^2\propto (\varepsilon _b r_b)^{2/3}\propto (1+z)$, or the mass-dispersion relation $M_b\propto \varepsilon _b^{-1}G^{-1}\sigma _b^5$ were identified, where $r_b\propto (1+z)^{-1}$ is the bulge size. For SMBHs, three evolution phases were found involving an initial rapid growth stage with a rising luminosity $L_B\propto (\varepsilon _b M_{\rm BH})^{4/5}$, a transition stage with a declining $L_B\propto \varepsilon _b^2 M_{\rm BH} \propto (1+z)^5$, and a dormant stage with $L_B\propto (\varepsilon _b M_{\rm BH})^{4/3}$. Our results suggest a rapid initial super-Eddington growth in a short period with a new redshift-dependent luminosity limit $L_X\propto\varepsilon _b^{4/5}M_{\rm BH}^{4/5}G^{-1/5}c$, in contrast to the Eddington limit. Analytical solutions were formulated for the BH mass function $\Phi _{\rm BH}$, active galactic nucleus (AGN) mass function $\Phi _{\rm AGN}$, and duty cycle U that predict $\Phi _L\propto L^{-1/5}$ for the faint-end luminosity function, $\Phi _{\rm AGN}\propto M^{-1/5}$ for small-mass-end AGN mass function $\Phi _L$, and $U\propto M^{-1/5}$ at high redshift.