No Arabic abstract
We model the cosmological co-evolution of galaxies and their central supermassive black holes (BHs) within a semi-analytical framework developed on the outputs of the Millennium Simulation (Croton et al., 2006; De Lucia & Blaizot, 2007). In this work, we analyze the model BH scaling relations, fundamental plane and mass function, and compare them with the most recent observational data. Furthermore, we extend the original code developed by Croton et al. (2006) to follow the evolution of the BH mass accretion and its conversion into radiation, and compare the derived AGN bolometric luminosity function with the observed one. We find, for the most part, a very good agreement between predicted and observed BH properties. Moreover, the model is in good agreement with the observed AGN number density in 0<z<5, provided it is assumed that the cold gas fraction accreted by BHs at high redshifts is larger than at low redshifts (Marulli et al., 2008).
The star formation rate (SFR) and black hole accretion rate (BHAR) functions are measured to be proportional to each other at z < ~3. This close correspondence between SF and BHA would naturally yield a BH mass-galaxy mass correlation, whereas a BH mass-bulge mass correlation is observed. To explore this apparent contradiction we study the SF in spheroid-dominated galaxies between z=1 and the present day. We use 903 galaxies from the COMBO-17 survey with M* >2x10^10M_sun, ultraviolet and infrared-derived SFRs from Spitzer and GALEX, and morphologies from GEMS HST/ACS imaging. Using stacking techniques, we find that <25% of all SF occurs in spheroid-dominated galaxies (Sersic index n>2.5), while the BHAR that we would expect if the global scalings held is three times higher. This rules out the simplest picture of co-evolution, in which SF and BHA trace each other at all times. These results could be explained if SF and BHA occur in the same events, but offset in time, for example at different stages of a merger event. However, one would then expect to see the corresponding star formation activity in early-stage mergers, in conflict with observations. We conclude that the major episodes of SF and BHA occur in different events, with the bulk of SF happening in isolated disks and most BHA occurring in major mergers. The apparent global co-evolution results from the regulation of the BH growth by the potential well of the galactic spheroid, which includes a major contribution from disrupted disk stars.
Feedback from accreting SMBHs is often identified as the main mechanism responsible for regulating star-formation in AGN host galaxies. However, the relationships between AGN activity, radiation, winds, and star-formation are complex and still far from being understood. We study scaling relations between AGN properties, host galaxy properties and AGN winds. We then evaluate the wind mean impact on the global star-formation history, taking into account the short AGN duty cycle with respect to that of star-formation. We first collect AGN wind observations for 94 AGN with detected massive winds at sub-pc to kpc spatial scales. We then fold AGN wind scaling relations with AGN luminosity functions, to evaluate the average AGN wind mass-loading factor as a function of cosmic time. We find strong correlations between the AGN molecular and ionised wind mass outflow rates and the AGN bolometric luminosity. The power law scaling is steeper for ionised winds (slope 1.29+/-0.38) than for molecular winds (0.76+/-0.06), meaning that the two rates converge at high bolometric luminosities. The molecular gas depletion timescale and the molecular gas fraction of galaxies hosting powerful AGN winds are 3-10 times shorter and smaller than those of main-sequence galaxies with similar SFR, stellar mass and redshift. These findings suggest that, at high AGN bolometric luminosity, the reduced molecular gas fraction may be due to the destruction of molecules by the wind, leading to a larger fraction of gas in the atomic ionised phase. The AGN wind mass-loading factor $eta=dot M_{OF}/SFR$ is systematically higher than that of starburst driven winds. Our analysis shows that AGN winds are, on average, powerful enough to clean galaxies from their molecular gas only in massive systems at z<=2, i.e. a strong form of co-evolution between SMBHs and galaxies appears to break down for the least massive galaxies.
We develop a formalism to model the luminosity functions (LFs) of radio-loud Active Galactic Nuclei (AGN) at GHz frequencies by the cosmological evolution of the supermassive black hole (SMBH). The mass function and Eddington ratio distributions of SMBHs have been taken as starting point for this analysis. Our approach is based on physical and phenomenological relations that allow to statistically calculate the radio luminosity of AGN cores, corrected for beaming effects, by linking it with the SMBH at their center, through the Fundamental Plane of black hole activity. Moreover, radio luminosity from extended jets and lobes is also computed through a power-law relationship that reflects the expected correlation between the inner radio core and the extended jets/lobes. Radio-loud AGN are further divided in two classes, characterized by different accretion modes onto the central BH. If the Eddington ratio is <0.01 they are called low-kinetic (LK-mode) AGN; if >0.01, they are called high-kinetic (HK-mode) AGN, being this critical value roughly corresponding to the transition between radiatively inefficient and efficient accretion flows. The few free parameters used in the present model are determined by fitting two different types of observational data sets: local LFs of radio-loud AGN at 1.4 GHz and differential number counts of extragalactic radio sources at 1.4 and 5 GHz. Our present model fits well almost all published data on LFs of LK-mode AGN and of the total AGN population up to redshifts z<=1.5 and, moreover, in the full range of luminosities currently probed by data. On the other hand, it tends to underestimate some recent measures of the LF of HK-mode AGN at low redshifts, but only at low radio luminosities. The present model is also able to provide a very good fit to almost all data on number counts of radio-loud sources at 1.4 and 5 GHz. (abridged)
We summarize what large surveys of the contemporary universe have taught us about the physics and phenomenology of the processes that link the formation and evolution of galaxies and their central supermassive black holes. We present a picture in which the population of AGN can be divided into two distinct populations. The Radiative-Mode AGN are associated with black holes that produce radiant energy powered by accretion at rates in excess of ~1% of the Eddington Limit. They are primarily associated with less massive black holes growing in high-density pseudo-bulges at a rate sufficient to produce the total mass budget in these black holes in ~10 Gyr. The circum-nuclear environment contains high density cold gas and associated star-formation. Major mergers are not the primary mechanism for transporting this gas inward; secular processes appear dominant. Stellar feedback will be generic in these objects and strong AGN feedback is seen only in the most powerful AGN. In Jet-Mode AGN the bulk of energetic output takes the form of collimated outflows (jets). These AGN are associated with the more massive black holes in more massive (classical) bulges and elliptical galaxies. Neither the accretion onto these black holes nor star-formation in their host bulge is significant today. These AGN are probably fueled by the accretion of slowly cooling hot gas that is limited by the feedback/heating provided by AGN radio sources. Surveys of the high-redshift universe are painting a similar picture. (Abridged).
It is well established that the properties of supermassive black holes and their host galaxies are correlated through scaling relations. While hydrodynamical cosmological simulations have begun to account for the co-evolution of BHs and galaxies, they typically have neglected the BH spin, even though it may play an important role in modulating the growth and feedback of BHs. Here we introduce a new sub-grid model for the BH spin evolution in the moving-mesh code {small AREPO} in order to improve the physical faithfulness of the BH modelling in galaxy formation simulations. We account for several different channels of spin evolution, in particular gas accretion through a Shakura-Sunyaev $alpha$-disc, chaotic accretion, and BH mergers. For BH feedback, we extend the IllustrisTNG model, which considers two different BH feedback modes, a thermal quasar mode for high accretion states and a kinetic mode for low Eddington ratios, with a self-consistent accounting of spin-dependent radiative efficiencies and thus feedback strength. We find that BHs with mass $M_{rm{bh}}lesssim 10^{8}, {rm M}_odot$ reach high spin values as they typically evolve in the coherent gas accretion regime. On the other hand, BHs with mass $M_{rm{bh}}gtrsim 10^{8}, {rm M}_odot$ have lower spins as BH mergers become more frequent, and their accretion discs fragment due to self-gravity, inducing chaotic accretion. We also explore the hypothesis that the transition between the quasar and kinetic feedback modes is mediated by the accretion mode of the BH disc itself, i.e.~the kinetic feedback mode is activated when the disc enters the self-gravity regime. We find excellent agreement between the galaxy and BH populations for this approach and the fiducial TNG model with no spin evolution. Furthermore, our new approach alleviates a tension in the galaxy morphology -- colour relation of the original TNG model.