Co-evolution Of Massive Black Holes Research Articles

Co-evolution of massive black holes and their host galaxies at high redshift: discrepancies from six cosmological simulations and the key role of JWST

ABSTRACT The James Webb Space Telescope will have the power to characterize high-redshift quasars at z ≥ 6 with an unprecedented depth and spatial resolution. While the brightest quasars at such redshift (i.e. with bolometric luminosity $L_{\rm bol}\geqslant 10^{46}\, \rm erg/s$) provide us with key information on the most extreme objects in the Universe, measuring the black hole (BH) mass and Eddington ratios of fainter quasars with $L_{\rm bol}= 10^{45}-10^{46}\, \rm erg\,s^{ -1}$ opens a path to understand the build-up of more normal BHs at z ≥ 6. In this paper, we show that the Illustris, TNG100, TNG300, Horizon-AGN, EAGLE, and SIMBA large-scale cosmological simulations do not agree on whether BHs at z ≥ 4 are overmassive or undermassive at fixed galaxy stellar mass with respect to the MBH − M⋆ scaling relation at z = 0 (BH mass offsets). Our conclusions are unchanged when using the local scaling relation produced by each simulation or empirical relations. We find that the BH mass offsets of the simulated faint quasar population at z ≥ 4, unlike those of bright quasars, represent the BH mass offsets of the entire BH population, for all the simulations. Thus, a population of faint quasars with $L_{\rm bol}= 10^{45}-10^{46}\, \rm erg\,s^{ -1}$ observed by JWST can provide key constraints on the assembly of BHs at high redshift. Moreover, this will help constraining the high-redshift regime of cosmological simulations, including BH seeding, early growth, and co-evolution with the host galaxies. Our results also motivate the need for simulations of larger cosmological volumes down to z ∼ 6, with the same diversity of subgrid physics, in order to gain statistics on the most extreme objects at high redshift.

Investigating the co-evolution of massive black holes in dual active galactic nuclei and their host galaxies via galaxy merger simulations

Major galaxy mergers can trigger nuclear activities and are responsible for high-luminosity quasi-stellar objects/active galactic nuclei (QSOs/AGNs). In certain circumstances, such mergers may cause dual active galactic nuclei (dAGN) phenomenon. This study investigates dAGN triggering and evolution of massive black holes (MBHs) during the merging processes using hydrodynamic code GADGET-2 to simulate several gas-rich major mergers at redshift $z=2$ and $3$, respectively. Results reveal that gas-rich major mergers can trigger significant nuclear activities after the second and third pericentric passages and the formation of dAGN with significant time duration ($\\sim$ 10-390 Myr). During the merging processes, galactic bulge evolves with time because of the rapid star formation in each (or both) galactic centers and initial mixing of stars in galactic disks due to violent relaxation. MBHs grow substantially due to accretion and finally merge into a bigger black hole. The growth of galactic bulges and corresponding increases of its velocity dispersions predate the growth of MBHs in the dAGN stages. The MBHs in these stages deviate below the relation between MBH mass and bulge mass (or velocity dispersion), and they revert to the relation after the final mergers due to the significant accretion that occurs mostly at a separation less than a few kpc. Then, the two MBHs merge with each other.

Testing the Evolutionary Link between Type 1 and Type 2 Quasars with Measurements of the Interstellar Medium

Abstract In a popular scenario for the coevolution of massive black holes and galaxies, major mergers of gas-rich galaxies fuel vigorous star formation and obscured (type 2) quasar activity until energy feedback from the active galactic nucleus (AGN) clears away the gas and dust to reveal an unobscured (type 1) quasar. Under this scenario, the precursor type 2 quasars should be more gas-rich than their type 1 counterparts, and both types of quasars are expected to be gas-deficient relative to normal, star-forming galaxies of similar stellar mass. We test this evolutionary hypothesis by investigating the infrared (∼1–500 μm) spectral energy distribution of 86 optically selected z < 0.5 type 2 quasars, matched in redshift and [O iii] luminosity to a comparison sample of type 1 quasars. Contrary to expectations, the gas content of the host galaxies of type 2 quasars is nearly indistinguishable from that of type 1 quasar hosts, and neither type exhibits the predicted deficit in gas relative to normal galaxies. The gas mass fraction of quasar hosts appears unaffected by the bolometric luminosity of the active nucleus, although their interstellar radiation field is preferentially higher than that of normal galaxies, potentially implicating AGN heating of the large-scale galactic dust.

A relationship between halo mass, cooling, active galactic nuclei heating and the co-evolution of massive black holes

We derive X-ray mass, luminosity, and temperature profiles for 45 galaxy clusters to explore relationships between halo mass, AGN feedback, and central cooling time. We find that radio--mechanical feedback power (referred to here as "AGN power") in central cluster galaxies correlates with halo mass as P$_{\rm mech}$ $\propto$ M$^{1.55\pm0.26}$, but only in halos with central atmospheric cooling times shorter than 1 Gyr. The trend of AGN power with halo mass is consistent with the scaling expected from a self-regulating AGN feedback loop, as well as with galaxy and central black hole co-evolution along the $M_{\rm BH} - \sigma$ relation. AGN power in clusters with central atmospheric cooling times longer than $\sim 1$ Gyr typically lies two orders of magnitude below those with shorter central cooling times. Galaxies centred in clusters with long central cooling times nevertheless experience ongoing and occasionally powerful AGN outbursts. We further investigate the impact of feedback on cluster scaling relations. We find $L-T$, and $M-T$ relations in clusters with direct evidence of feedback which are steeper than self-similar, but not atypical compared to previous studies of the full cluster population. While the gas mass rises, the stellar mass remains nearly constant with rising total mass, consistent with earlier studies. This trend is found regardless of central cooling time, implying tight regulation of star formation in central galaxies as their halos grew, and long-term balance between AGN heating and atmospheric cooling. Our scaling relations are presented in forms that can be incorporated easily into galaxy evolution models.

LINKING THE SPIN EVOLUTION OF MASSIVE BLACK HOLES TO GALAXY KINEMATICS

We present the results of a semianalytical model that evolves the masses and spins of massive black holes together with the properties of their host galaxies along the cosmic history. As a consistency check, our model broadly reproduces a number of observations, e.g. the cosmic star formation history, the black hole mass and luminosity function and the galaxy mass function at low redshift, the black hole to bulge mass relation, and the morphological distribution at low redshift. For the first time in a semianalytical investigation, we relax the simplifying assumptions of perfect coherency or perfect isotropy of the gas fueling the black holes. The dynamics of gas is instead linked to the morphological properties of the host galaxies, resulting in different spin distributions for black holes hosted in different galaxy types. We compare our results with the observed sample of spin measurements obtained through broad K$\alpha$ iron line fitting. The observational data disfavor both accretion along a fixed direction and isotropic fueling. Conversely, when the properties of the accretion flow are anchored to the kinematics of the host galaxy, we obtain a good match between theoretical expectations and observations. A mixture of coherent accretion and phases of activity in which the gas dynamics is similar to that of the stars in bulges (i.e., with a significant velocity dispersion superimposed to a net rotation) best describes the data, adding further evidence in support to the coevolution of massive black holes and their hosts.

MASSIVE BLACK HOLES IN CENTRAL CLUSTER GALAXIES

We explore how the co-evolution of massive black holes (MBHs) and galaxies is affected by environmental effects, addressing in particular MBHs hosted in the central galaxies of clusters (we will refer to these galaxies in general as 'CGs'). Recently the sample of MBHs in CGs with dynamically measured masses has increased, and it has been suggested that these MBH masses (M_BH) deviate from the expected correlations with velocity dispersion (sigma) and mass of the bulge (M_bulge) of the host galaxy: MBHs in CGs appear to be `over-massive'. This discrepancy is more pronounced when considering the M_BH-sigma relation than the M_BH-M_bulge one. We show that this behavior stems from a combination of two natural factors, (i) that CGs experience more mergers involving spheroidal galaxies and their MBHs, and (ii) that such mergers are preferentially gas-poor. We use a combination of analytical and semi-analytical models to investigate the MBH-galaxy co-evolution in different environments and find that the combination of these two factors explains the trends observed in current data-sets.

Dark energy, co-evolution of massive black holes with galaxies, and ASTROD-GW

The detection of low frequency band (100nHz–100mHz) and very low frequency band (300pHz–100nHz) gravitational waves (GWs) is important for exploration of the equation of state of dark energy and the co-evolution of massive black holes (MBHs) with galaxies. Most galaxies are believed to have a massive black hole in the galactic core. In the formation of these black holes, merging and accretion are the two main processes. Merging of massive black holes generate GWs which could be detected by space GW detectors and Pulsar Timing Arrays (PTAs) to cosmological distances. LISA (Laser-Interferometric Space Antenna) is most sensitive to the frequency band 1mHz–100mHz, ASTROD-GW (ASTROD [Astrodynamical Space Test of Relativity using Optical Devices] optimized for Gravitational Wave detection) is most sensitive to the frequency band 100nHz–1mHz and PTAs are most sensitive to the frequency band 300pHz–100nHz. In this paper, we discuss the sensitivities and outlooks of detection of GWs from binary massive black holes in these frequency bands with an emphasis on ASTROD-GW. The GWs generated by the inspirals, merging and subsequent ringdowns of binary black holes are standard sirens to the cosmological distance. Using GW observations, we discuss the methods for determining the equation of state of dark energy and for testing the co-evolution models of massive black holes. ASTROD-GW is an optimization of ASTROD to focus on the goal of detection of GWs. The mission orbits of the 3 spacecraft forming a nearly equilateral triangular array are chosen to be near the Sun-Earth Lagrange points L3, L4 and L5. The 3 spacecraft range interferometrically with one another with arm length about 260million kilometers. With 52 times longer in arm length compared to that of LISA, the strain detection sensitivity is 52 times better toward larger wavelength. The scientific aim is focused for gravitational wave detection at low frequency. The science goals include detection of GWs from MBHs, and Extreme-Mass-Ratio Black Hole Inspirals (EMRI), and using these observations to find the evolution of the equation of state of dark energy and to explore the co-evolution of massive black holes with galaxies.

Mass Function of Binary Massive Black Holes in Active Galactic Nuclei

Abstract If the activity of active galactic nuclei (AGNs) is predominantly induced by major galaxy mergers, then a significant fraction of AGNs should harbor binary massive black holes in their centers. We have studied the mass function of binary massive black holes in nearby AGNs based on the observed AGN black-hole mass function and the theory of evolution of binary massive black holes interacting with a massive circumbinary disk within the framework of the coevolution of massive black holes and their host galaxies. The circumbinary disk is assumed to be steady, axisymmetric, geometrically thin, self-regulated, self-gravitating but nonfragmenting with a fixed fraction of the Eddington accretion rate, which is typically one tenth of the Eddington value. The timescale of orbital decay is then estimated to be $\sim\ $10$^{8}\ $yr for equal mass black holes, being independent of the black-hole mass, semimajor axis, and viscosity parameter, but dependent on the black-hole mass ratio, Eddington ratio, and mass-to-energy conversion efficiency. This makes it possible for any binary massive black holes to merge within the Hubble time under the binary–disk interaction. We find that (1.5%$\ \pm\ $0.6%) of the total number of nearby AGNs for the equal-mass ratio and (1.3%$\ \pm\ $0.5%) for the one-to-ten mass ratio have close binary massive black holes with an orbital period of less than 10 yr in their centers, detectable with ongoing highly sensitive X-ray monitors, such as Monitor of All-sky X-ray Image and/or Swift/Burst Alert Telescope. Assuming that all binary massive black holes have the equal-mass ratio, about 10% of AGNs with black-hole masses of 10$^{6.5-7}\ M_{\odot}$ have close binaries, and thus provide the best chance to detect them.

Radiative Feedback from Massive Black Holes in Elliptical Galaxies: AGN Flaring and Central Starburst Fueled by Recycled Gas

The importance of the radiative feedback from massive black holes at the centers of elliptical galaxies is not in doubt, given the well established relations among electromagnetic output, black hole mass and galaxy optical luminosity. We show how this AGN radiative output affects the hot ISM of an isolated elliptical galaxy with the aid of a high-resolution hydrodynamical code, where the cooling and heating functions include photoionization plus Compton heating. We find that radiative heating is a key factor in the self-regulated coevolution of massive black holes and their host galaxies and that 1) the mass accumulated by the central black hole is limited by feedback to the range observed today, and 2) relaxation instabilities occur so that duty cycles are small enough (~0.03) to account for the very small fraction of massive ellipticals observed to be in the "on" -QSO- phase, when the accretion luminosity approaches the Eddington luminosity. The duty cycle of the hot bubbles inflated at the galaxy center during major accretion episodes is of the order of 0.1-0.4. Major accretion episodes caused by cooling flows in the recycled gas produced by normal stellar evolution trigger nuclear starbursts coincident with AGN flaring. During such episodes the central sources are often obscured; but overall, in the bursting phase (1<z<3), the duty cycle of the black hole in its "on" phase is of the order of percents and it is unobscured approximately one-third of the time. Mechanical energy output from non-relativistic gas winds integrates to 2.3 10^{59} erg, with most of it caused by broadline AGN outflows. [abridged]

Radiative feedback from quasars and the growth of massive black holes in stellar spheroids

We discuss the importance of feedback via photoionization and Compton heating on the co-evolution of massive black holes (MBHs) at the center of spheroidal galaxies and their stellar and gaseous components. We first assess the energetics of the radiative feedback from a typical quasar on the ambient interstellar medium. We then demonstrate that the observed Mbh-sigma relation could be established following the conversion of most of the gas of an elliptical progenitor into stars, specifically when the gas-to-stars mass ratio in the central regions has dropped to a low level ~0.01 or less, so that gas cooling is no longer able to keep up with the radiative heating by the growing central MBH. A considerable amount of the remaining gas will be expelled and both MBH accretion and star formation will proceed at significantly reduced rates thereafter, in agreement with observations of present day ellipticals. We find further support for this scenario by evolving over an equivalent Hubble time a simple, physically based toy model that additionally takes into account the mass and energy return for the spheroid evolving stellar population, a physical ingredient often neglected in similar approaches.

X-raying cosmic star formation history

Based on our previous work about a co-evolution of massive black holes (MBH) and their host spheroids, we estimate the cosmic star formation history associated with AGN accretion by ROSAT X-ray All sky surveys. We show: 1) the total amount of star formation associated with MBH growth is at least half of the net star formation at high redshift, which probably totally missed by the current UV/optical deep surveys; 2) the FIR emission from the dust heated by star formation on-going during the MBH growth could be a significant energy source for SCUBA populations; 3)the peak redshift of the massive spheroid formation in this case is around 1.5, not necessary to be at much higher redshift.