Existing models of galaxy formation have not yet explained striking correlations between structure and star-formation activity in galaxies, notably the sloped and moving boundaries that divide star-forming from quenched galaxies in key structural diagrams. This paper uses these and other relations to ``reverse-engineer'' the quenching process for central galaxies. The basic idea is that star-forming galaxies with larger radii (at a given stellar mass) have lower black-hole masses due to lower central densities. Galaxies cross into the green valley when the cumulative effective energy radiated by their black hole equals $\sim4\times$ their halo-gas binding energy. Since larger-radii galaxies have smaller black holes, one finds they must evolve to higher stellar masses in order to meet this halo-energy criterion, which explains the sloping boundaries. A possible cause of radii differences among star-forming galaxies is halo concentration. The evolutionary tracks of star-forming galaxies are nearly parallel to the green-valley boundaries, and it is mainly the sideways motions of these boundaries with cosmic time that cause galaxies to quench. BH-scaling laws for star-forming, quenched, and green-valley galaxies are different, and most BH mass growth takes place in the green valley. Implications include: the radii of star-forming galaxies are an important second parameter in shaping their black holes; black holes are connected to their halos but in different ways for star-forming, quenched, and green-valley galaxies; and the same BH-halo quenching mechanism has been in place since $z \sim 3$. We conclude with a discussion of black hole-galaxy co-evolution, the origin and interpretation of BH scaling laws.

We study black hole - host galaxy correlations, and the relation between the over-massiveness (the distance from the average $M_{BH}-\sigma$ relation) of super-massive black holes (SMBHs) and star formation histories of their host galaxies in the Illustris and TNG100 simulations. We find that both simulations are able to produce black hole scaling relations in general agreement with observations at $z=0$, but with noticeable discrepancies. Both simulations show an offset from the observations for the $M_{BH}-\sigma$ relation, and the relation between $M_{BH}$ and the Sersic index. The relation between $M_{BH}$ and stellar mass $M_*$ is tighter than the observations, especially for TNG100. For massive galaxies in both simulations, the hosts of over-massive SMBHs (those above the mean $M_{BH}-\sigma$ relation) tend to have larger Sersic indices and lower baryon conversion efficiency, suggesting a multidimensional link between SMBHs and properties of their hosts. In Illustris, the hosts of over-massive SMBHs have formed earlier and have lower present-day star formation rates, in qualitative agreement with the observations for massive galaxies with $\sigma>100 \rm km/s$. For low-mass galaxies, such a correlation still holds in Illustris but does not exist in the observed data. For TNG100, the correlation between SMBH over-massiveness and star formation history is much weaker. The hosts of over-massive SMBHs generally have consistently larger star formation rates throughout history. These galaxies have higher stellar mass as well, due to the strong $M_{BH}-M_*$ correlation. Our findings show that simulated SMBH scaling relations and correlations are sensitive to features in the modeling of SMBHs.

Supermassive black hole feedback is thought to be responsible for the lack of star formation, or quiescence, in a significant fraction of galaxies. We explore how observable correlations between the specific star formation rate (sSFR), stellar mass (M$_{\rm{star}}$), and black hole mass (M$_{\rm{BH}}$) are sensitive to the physics of black hole feedback in a galaxy formation model. We use the IllustrisTNG simulation suite, specifically the TNG100 simulation and ten model variations that alter the parameters of the black hole model. Focusing on central galaxies at $z = 0$ with M$_{\rm{star}} > 10^{10}$ M$_{\odot}$, we find that the sSFR of galaxies in IllustrisTNG decreases once the energy from black hole kinetic winds at low accretion rates becomes larger than the gravitational binding energy of gas within the galaxy stellar radius. This occurs at a particular M$_{\rm{BH}}$ threshold above which galaxies are found to sharply transition from being mostly star-forming to mostly quiescent. As a result of this behavior, the fraction of quiescent galaxies as a function of M$_{\rm{star}}$ is sensitive to both the normalization of the M$_{\rm{BH}}$-M$_{\rm{star}}$ relation and the M$_{\rm{BH}}$ threshold for quiescence in IllustrisTNG. Finally, we compare these model results to observations of 91 central galaxies with dynamical M$_{\rm{BH}}$ measurements with the caveat that this sample is not representative of the whole galaxy population. While IllustrisTNG reproduces the observed trend that quiescent galaxies host more massive black holes, the observations exhibit a broader scatter in M$_{\rm{BH}}$ at a given M$_{\rm{star}}$ and show a smoother decline in sSFR with M$_{\rm{BH}}$.

We use the Henriques et al. (2015) version of the Munich galaxy formation model (L-GALAXIES) to investigate why the halo and stellar mass scales above which galaxies are quenched are constant with redshift and coincide with the scale where baryons are most efficiently converted into stars. This model assumes that central galaxies are quenched by AGN feedback when hot halo gas accretes onto a supermassive black hole. Nevertheless, we find that supernova (SN) feedback sets both mass scales. As haloes grow above a threshold mass, SNe can no longer eject material so their hot gas content increases, enhancing the cooling rate onto the central galaxy, its cold gas content, its star formation rate and the growth rate of its central black hole. Strong AGN feedback terminates this short-lived phase by suppressing the fuel supply for star formation. Despite strong evolution of the halo mass - temperature relation, quenching occurs at a redshift-independent halo and stellar mass which coincides with the mass where baryons have been converted into stars with maximal efficiency. These regularities and coincidences are a result of the specific parameters selected by MCMC tuning of the model to fit the observed abundance and passive fraction of galaxies over the redshift range 0<z<3. Thus they are required by the observed evolution of the galaxy population, at least in the context of models of this type.

We present a relationship between the black hole mass, stellar mass, and star formation rate of a diverse group of 91 galaxies with dynamically-measured black hole masses. For our sample of galaxies with a variety of morphologies and other galactic properties, we find that the specific star formation rate is a smoothly decreasing function of the ratio between black hole mass and stellar mass, or what we call the specific black hole mass. In order to explain this relation, we propose a physical framework where the gradual suppression of a galaxy's star formation activity results from the adjustment to an increase in specific black hole mass and, accordingly, an increase in the amount of heating. From this framework, it follows that at least some galaxies with intermediate specific black hole masses are in a steady state of partial quiescence with intermediate specific star formation rates, implying that both transitioning and steady-state galaxies live within this region known as the "green valley." With respect to galaxy formation models, our results present an important diagnostic with which to test various prescriptions of black hole feedback and its effects on star formation activity.

Roughly half of all stars reside in galaxies without significant ongoing star formation. However, galaxy formation models indicate that it is energetically challenging to suppress the cooling of gas and the formation of stars in galaxies that lie at the centers of their dark matter halos. In this Letter, we show that the dependence of quiescence on black hole and stellar mass is a powerful discriminant between differing models for the mechanisms that suppress star formation. Using observations of 91 star-forming and quiescent central galaxies with directly-measured black hole masses, we find that quiescent galaxies host more massive black holes than star-forming galaxies with similar stellar masses. This observational result is in qualitative agreement with models that assume that effective, more-or-less continuous AGN feedback suppresses star formation, strongly suggesting the importance of the black hole in producing quiescence in central galaxies.

We use the semi-analytic model developed by Henriques et al. (2015) to explore the origin of star formation history diversity for galaxies that lie at the centre of their dark matter haloes and have present-day stellar masses in the range 5-8 $\times$ 10$^{10}$ M$_{\odot}$, similar to that of the Milky Way. In this model, quenching is the dominant physical mechanism for introducing scatter in the growth histories of these galaxies. We find that present-day quiescent galaxies have a larger variety of growth histories than star-formers since they underwent 'staggered quenching' - a term describing the correlation between the time of quenching and present-day halo mass. While halo mass correlates broadly with quiescence, we find that quiescence is primarily a function of black hole mass, where galaxies quench when heating from their active galactic nuclei becomes sufficient to offset the redshift-dependent cooling rate. In this model, the emergence of a prominent quiescent population is the main process that flattens the stellar mass-halo mass relation at mass scales at or above that of the Milky Way.

The first generation of stars produces a background of Lyman-Werner (LW) radiation which can photo-dissociate molecular hydrogen, increasing the mass of dark matter halos required to host star formation. Previous studies have determined the critical mass required for efficient molecular cooling with a constant LW background. However, the true background is expected to increase rapidly at early times. Neglecting this evolution could underestimate star formation in small halos that may have started to cool in the past when the LW intensity was much lower. Background evolution is a large source of uncertainty in pre-reionization predictions of the cosmological 21cm signal, which can be observed with future radio telescopes. To address this, we perform zero-dimentional one-zone calculations that follow the density, chemical abundances, and temperature of gas in the central regions of dark matter halos, including hierarchical growth and an evolving LW background. We begin by studying the physics of halos subjected to a background that increases exponentially with redshift. We find that when the intensity increases more slowly than $J_{\rm LW}(z) \propto 10^{-z/5}$, cooling in the past is a relatively small effect. We then self-consistently compute the cosmological LW background over $z=15-50$ and find that cooling in the past due to an evolving background has a modest impact. Finally, we compare these results to three-dimensional hydrodynamical cosmological simulations with varying LW histories. While only a small number of halos were simulated, the results are consistent with our one-zone calculations.