Black Hole -- Galaxy Correlations in Simba

We examine the co-evolution of galaxies and supermassive black holes in the Simba cosmological hydrodynamic simulation. Simba grows black holes via gravitational torque-limited accretion from cold gas and Bondi accretion from hot gas, while feedback from black holes is modeled in radiative and jet modes depending on the Eddington ratio ($f_{Edd}$). Simba shows generally good agreement with local studies of black hole properties, such as the black hole mass--stellar velocity dispersion ($M_{BH}-\sigma$) relation, 2 the black hole accretion rate vs. star formation rate (BHAR--SFR), and the black hole mass function. $M_{BH}-\sigma$ evolves such that galaxies at a given $M_{BH}$ have higher $\sigma$ at higher redshift, consistent with no evolution in $M_{BH}-M_*$. For $M_{BH}<\sim 10^8 M_\odot$, $f_{Edd}$ is anti-correlated with $M_{BH}$ since the BHAR is approximately independent of $M_{BH}$, while at higher masses $f_{Edd}-M_{BH}$ flattens and has a larger scatter. BHAR vs. SFR is invariant with redshift, but $f_{Edd}$ drops steadily with time at a given $M_{BH}$, such that all but the most massive black holes are accreting in a radiatively efficient mode at $z>\sim 2$. The black hole mass function amplitude decreases with redshift and is locally dominated by quiescent galaxies for $M_{BH}>10^{8}M_{\odot}$, but for $z>\sim 1$ star forming galaxies dominate at all $M_{BH}$. The $z=0$ $f_{Edd}$ distribution is roughly lognormal with a peak at $f_{Edd}<\sim 0.01$ as observed, shifting to higher $f_{Edd}$ at higher redshifts. Finally, we study the dependence of black hole properties with \HI\ content and find that the correlation between gas content and star formation rate is modulated by black hole properties, such that higher SFR galaxies at a given gas content have smaller black holes with higher $f_{Edd}$