Black holes with millions or even billions of solar masses are found in the centers of most galaxies. The most powerfully active nuclear black holes are in quasars and these have been spotted as far away as the epoch when the universe was less than a billion years old, suggesting that the galaxy-black hole symbiosis was already underway at this early time. Accreting black holes can emit powerful jets or winds that reverse the accretion and drive material outward, sometimes quenching the star formation. These and other lines of evidence help to clarify the co-evolution mechanisms between black holes and galaxies and reveal the joint evolution of the galaxy and the supermassive black hole populations. CfA astronomers Lars Hernquist and Rainer Weinberger and their colleagues used the large-scale hydrodynamic simulation called IllustrisTNG to trace the development of galaxies and their black holes. The code is able to model the evolution of a wide range of black hole and galaxy properties as the universe ages. They successfully reproduce the observed correlation between star formation rate and galaxy mass. They find, among numerous other trends, that quiescent galaxies (those no longer actively making stars) first go through a phase of shrinking in size before they undergo a quenching event; they also find that in the cosmic epoch of peak star formation (about ten billion years ago) as many as twenty percent of galaxies hosted an active supermassive black hole. Reference: “Linking Galaxy Structural Properties and Star Formation Activity to Black Hole Activity with IllustrisTNG” by Melanie Habouzit, Shy Genel, Rachel S. Somerville, Dale Kocevski, Michaela Hirschmann, Avishai Dekel, Ena Choi, Dylan Nelson, Annalisa Pillepich, Paul Torrey, Lars Hernquist, Mark Vogelsberger, Rainer Weinberger, and Volker Springel, 16 January, MNRAS.DOI: 10.1093/mnras/stz102