No Arabic abstract
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.
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.
The spin angular momentum S of a supermassive black hole (SBH) precesses due to torques from orbiting stars, and the stellar orbits precess due to dragging of inertial frames by the spinning hole. We solve the coupled post-Newtonian equations describing the joint evolution of S and the stellar angular momenta Lj, j = 1...N in spherical, rotating nuclear star clusters. In the absence of gravitational interactions between the stars, two evolutionary modes are found: (1) nearly uniform precession of S about the total angular momentum vector of the system; (2) damped precession, leading, in less than one precessional period, to alignment of S with the angular momentum of the rotating cluster. Beyond a certain distance from the SBH, the time scale for angular momentum changes due to gravitational encounters between the stars is shorter than spin-orbit precession times. We present a model, based on the Ornstein-Uhlenbeck equation, for the stochastic evolution of star clusters due to gravitational encounters and use it to evaluate the evolution of S in nuclei where changes in the Lj are due to frame dragging close to the SBH and to encounters farther out. Long-term evolution in this case is well described as uniform precession of the SBH about the clusters rotational axis, with an increasingly important stochastic contribution when SBH masses are small. Spin precessional periods are predicted to be strongly dependent on nuclear properties, but typical values are 10-100 Myr for low-mass SBHs in dense nuclei, 100 Myr - 10 Gyr for intermediate mass SBHs, and > 10 Gyr for the most massive SBHs. We compare the evolution of SBH spins in stellar nuclei to the case of torquing by an inclined, gaseous accretion disk.
(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.
Interaction of a binary supermassive black hole with stars in a galactic nucleus can result in changes to all the elements of the binarys orbit, including the angles that define its orientation. If the nucleus is rotating, the orientation changes can be large, causing large changes in the binarys orbital eccentricity as well. We present a general treatment of this problem based on the Fokker-Planck equation for f, defined as the probability distribution for the binarys orbital elements. First- and second-order diffusion coefficients are derived for the orbital elements of the binary using numerical scattering experiments, and analytic approximations are presented for some of these coefficients. Solutions of the Fokker-Planck equation are then derived under various assumptions about the initial rotational state of the nucleus and the binary hardening rate. We find that the evolution of the orbital elements can become qualitatively different when we introduce nuclear rotation: 1) the orientation of the binarys orbit evolves toward alignment with the plane of rotation of the nucleus; 2) binary orbital eccentricity decreases for aligned binaries and increases for counter-aligned ones. We find that the diffusive (random-walk) component of a binarys evolution is small in nuclei with non-negligible rotation, and we derive the time-evolution equations for the semimajor axis, eccentricity and inclination in that approximation. The aforementioned effects could influence gravitational wave production as well as the relative orientation of host galaxies and radio jets.
An extraordinary recent development in astrophysics was the discovery of the fossil relationship between central black hole mass and the stellar mass of galactic bulges. The physical process underpinning this relationship has become known as feedback. The Chandra X-ray Observatory was instrumental in realizing the physical basis for feedback, by demonstrating a tight coupling between the energy released by supermassive black holes and the gaseous structures surrounding them. This white paper discusses how a great leap forward in X-ray collecting area and spectral resolution will allow a qualitatively new way of studying how feedback from black holes influenced the growth of structure.