No Arabic abstract
(Abridged) We present a method that tracks the growth of supermassive black holes (BHs) and the feedback from AGN in cosmological simulations. Our model is a substantially modified version of the one by Springel et al. (2005). Because cosmological simulations lack both the resolution and the physics to model the multiphase interstellar medium, they tend to strongly underestimate the Bondi-Hoyle accretion rate. To allow low-mass BHs to grow, it is therefore necessary to increase the predicted Bondi-Hoyle rates in dense gas by large, ad-hoc factors. We explore the physical regimes where the use of such factors is reasonable, and through this introduce a new prescription for gas accretion. Feedback from AGN is modeled by coupling a fraction of the rest-mass energy of the accreted gas thermally into the surrounding medium. We describe the implementation as well as the limitations of the model and motivate all the changes relative to previous work. We investigate the robustness of the predictions for the cosmic star formation history, the redshift zero cosmic BH density, BH scaling relations, and galaxy specific star formation rates. We find that the freedom introduced by the need to increase the predicted accretion rates, the standard procedure in the literature, is the most significant source of uncertainty. Our simulations demonstrate that supermassive BHs are able to regulate their growth by releasing a fixed amount of energy for a given halo mass, independent of the assumed efficiency of AGN feedback, which sets the normalization of the BH scaling relations. Regardless of whether BH seeds are initially placed above or below the BH scaling relations they grow onto the same relations. AGN feedback efficiently suppresses star formation in high-mass galaxies.
We develop a simple evolutionary scenario for the growth of supermassive black holes (BHs), assuming growth due to accretion only, to learn about the evolution of the BH mass function from $z=3$ to 0 and from it calculate the energy budgets of different modes of feedback. We tune the parameters of the model by matching the derived X-ray luminosity function (XLF) with the observed XLF of active galactic nuclei. We then calculate the amount of comoving kinetic and bolometric feedback as a function of redshift, derive a kinetic luminosity function and estimate the amount of kinetic feedback and $PdV$ work done by classical double Fanaroff-Riley II (FR II) radio sources. We also derive the radio luminosity function for FR IIs from our synthesized population and set constraints on jet duty cycles. Around 1/6 of the jet power from FR II sources goes into $PdV$ work done in the expanding lobes during the time the jet is on. Anti hierarchical growth of BHs is seen in our model due to addition of an amount of mass being accreted on to all BHs independent of the BH mass. The contribution to the total kinetic feedback by active galaxies in a low accretion, kinetically efficient mode is found to be the most significant at $z<1.5$. FR II feedback is found to be a significant mode of feedback above redshifts $zsim 1.5$, which has not been highlighted by previous studies.
We present results from a new set of 30 cosmological simulations of galaxy clusters, including the effects of radiative cooling, star formation, supernova feedback, black hole growth and AGN feedback. We first demonstrate that our AGN model is capable of reproducing the observed cluster pressure profile at redshift, z~0, once the AGN heating temperature of the targeted particles is made to scale with the final virial temperature of the halo. This allows the ejected gas to reach larger radii in higher-mass clusters than would be possible had a fixed heating temperature been used. Such a model also successfully reduces the star formation rate in brightest cluster galaxies and broadly reproduces a number of other observational properties at low redshift, including baryon, gas and star fractions; entropy profiles outside the core; and the X-ray luminosity-mass relation. Our results are consistent with the notion that the excess entropy is generated via selective removal of the densest material through radiative cooling; supernova and AGN feedback largely serve as regulation mechanisms, moving heated gas out of galaxies and away from cluster cores. However, our simulations fail to address a number of serious issues; for example, they are incapable of reproducing the shape and diversity of the observed entropy profiles within the core region. We also show that the stellar and black hole masses are sensitive to numerical resolution, particularly the gravitational softening length; a smaller value leads to more efficient black hole growth at early times and a smaller central galaxy.
This is the first in a series of papers in which we study the application of spectroastrometry in the context of gas kinematical studies aimed at measuring the mass of supermassive black holes. The spectroastrometrical method consists in measuring the photocenter of light emission in different wavelength or velocity channels. In particular we explore the potential of spectroastrometry of gas emission lines in galaxy nuclei to constrain the kinematics of rotating gas disks and to measure the mass of putative supermassive black holes. By means of detailed simulations and test cases, we show that the fundamental advantage of spectroastrometry is that it can provide information on the gravitational potential of a galaxy on scales significantly smaller (~ 1/10) than the limit imposed by the spatial resolution of the observations. We then describe a simple method to infer detailed kinematical informations from spectroastrometry in longslit spectra and to measure the mass of nuclear mass concentrations. Such method can be applied straightforwardly to integral field spectra, which do not have the complexities due to a partial spatial covering of the source in the case of longslit spectra.
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.
The masses of 68 supermassive black holes (SMBHs) in nearby (z<0.15) active galactic nuclei (AGNs) detected by the INTEGRAL observatory in the hard X-ray energy band (17-60 keV) outside the Galactic plane (|b| > 5 degrees) have been estimated. Well-known relations between the SMBH mass and (1) the infrared luminosity of the stellar bulge (from 2MASS data) and (2) the characteristics of broad emission lines (from RTT-150 data) have been used. A comparison with the more accurate SMBH mass estimates obtained by the reverberation-mapping technique and from direct dynamical measurements is also made for several objects. The SMBH masses derived from the correlation with the bulge luminosity turn out to be systematically higher than the estimates made by other methods. The ratio of the bolometric luminosity to the critical Eddington luminosity has been found for all AGNs. It ranges from 1 to 100% for the overwhelming majority of objects.