Co-evolution Of Supermassive Black Holes Research Articles

OUTFLOW AND METALLICITY IN THE BROAD-LINE REGION OF LOW-REDSHIFT ACTIVE GALACTIC NUCLEI

ABSTRACT Outflows in active galactic nuclei (AGNs) are crucial to understand in investigating the co-evolution of supermassive black holes (SMBHs) and their host galaxies since outflows may play an important role as an AGN feedback mechanism. Based on archival UV spectra obtained with the Hubble Space Telescope and IUE, we investigate outflows in the broad-line region (BLR) in low-redshift AGNs (z < 0.4) through detailed analysis of the velocity profile of the C iv emission line. We find a dependence of the outflow strength on the Eddington ratio and the BLR metallicity in our low-redshift AGN sample, which is consistent with earlier results obtained for high-redshift quasars. These results suggest that BLR outflows, gas accretion onto SMBHs, and past star formation activity in host galaxies are physically related in low-redshift AGNs as in powerful high-redshift quasars.

The Keck OSIRIS Nearby AGN (KONA) Survey: AGN Fueling and Feedback

In an effort to better constrain the relevant physical processes dictating the co-evolution of supermassive black holes and the galaxies in which they reside we turn to local Seyfert AGN. It is only with these local AGN that we can reach the spatial resolution needed to adequately characterize the inflow and outflow mechanisms thought to be the driving forces in establishing the relationship between black holes and their host galaxies at higher redshift. We present the first results from the KONA (Keck OSIRIS Nearby AGN) survey, which takes advantage of the integral field unit OSIRIS plus laser and natural guide star adaptive optics to probe down to scales of 5-30 parsecs in a sample of 40 local Seyfert galaxies. With these K-band data we measure the two-dimensional distribution and kinematics of the nuclear stars, molecular gas, and ionized gas within the central few hundred parsecs.

The dust masses of powerful radio galaxies: clues to the triggering of their activity

Abstract We use deep Herschel Space Observatory observations of a 90 per cent complete sample of 32 intermediate-redshift 2Jy radio galaxies (0.05 < z < 0.7) with strong emission lines to estimate the dust masses of their host galaxies and thereby investigate the triggering mechanisms for their quasar-like AGN. The dust masses derived for the radio galaxies (7.2 × 105 < Md < 2.6 × 108 M⊙) are intermediate between those of quiescent elliptical galaxies on the one hand, and ultraluminous infrared galaxies (ULIRGs) on the other. Consistent with simple models for the co-evolution of supermassive black holes and their host galaxies, these results suggest that most radio galaxies represent the late time re-triggering of AGN activity via mergers between the host giant elliptical galaxies and companion galaxies with relatively low gas masses. However, a minority of the radio galaxies in our sample (∼20 per cent) have high, ULIRG-like dust masses, along with evidence for prodigious star formation activity. The latter objects are more likely to have been triggered in major, gas-rich mergers that represent a rapid growth phase for both their host galaxies and their supermassive black holes.

THE COEVOLUTION OF SUPERMASSIVE BLACK HOLES AND MASSIVE GALAXIES AT HIGH REDSHIFT

We exploit the recent, wide samples of far-infrared (FIR) selected galaxies followed-up in X rays and of X-ray/optically selected active galactic nuclei (AGNs) followed-up in the FIR band, along with the classic data on AGN and stellar luminosity functions at high redshift z>1.5, to probe different stages in the coevolution of supermassive black holes (BHs) and host galaxies. The results of our analysis indicate the following scenario: (i) the star formation in the host galaxy proceeds within a heavily dust-enshrouded medium at an almost constant rate over a timescale ~0.5-1 Gyr, and then abruptly declines due to quasar feedback; over the same timescale, (ii) part of the interstellar medium loses angular momentum, reaches the circum-nuclear regions at a rate proportional to the star formation and is temporarily stored into a massive reservoir/proto-torus wherefrom it can be promptly accreted; (iii) the BH grows by accretion in a self-regulated regime with radiative power that can slightly exceed the Eddington limit L/L_Edd< 4, particularly at the highest redshifts; (iv) for massive BHs the ensuing energy feedback at its maximum exceeds the stellar one and removes the interstellar gas, thus stopping the star formation and the fueling of the reservoir; (v) afterwards, if the latter has retained enough gas, a phase of supply-limited accretion follows exponentially declining with a timescale of about 2 e-folding times. We show that the ratio of the FIR luminosity of the host galaxy to the bolometric luminosity of the AGN maps the various stages of the above sequence. Finally, we discuss how the detailed properties and the specific evolution of the reservoir can be investigated via coordinated, high-resolution observations of starforming, strongly-lensed galaxies in the (sub-)mm band with ALMA and in the X-ray band with Chandra and the next generation X-ray instruments.

REFINING THEMBH-VcSCALING RELATION WITH H I ROTATION CURVES OF WATER MEGAMASER GALAXIES

Black hole - galaxy scaling relations provide information about the coevolution of supermassive black holes and their host galaxies. We compare the black hole mass - circular velocity (MBH - Vc) relation with the black hole mass - bulge stellar velocity dispersion (MBH - sigma) relation, to see whether the scaling relations can passively emerge from a large number of mergers, or require a physical mechanism, such as feedback from an active nucleus. We present VLA H I observations of five galaxies, including three water megamaser galaxies, to measure the circular velocity. Using twenty-two galaxies with dynamical MBH measurements and Vc measurements extending to large radius, our best-fit MBH - Vc relation, log MBH = alpha + beta log(Vc /200 km s^-1), yields alpha = 7.43+/-0.13, beta = 3.68+1.23/-1.20, and intrinsic scatter epsilon_int = 0.51+0.11/-0.09. The intrinsic scatter may well be higher than 0.51, as we take great care to ascribe conservatively large observational errors. We find comparable scatter in the MBH - sigma relations, epsilon_int = 0.48+0.10/-0.08, while pure merging scenarios would likely result in a tighter scaling with the dark halo (as traced by Vc) than baryonic (sigma) properties. Instead, feedback from the active nucleus may act on bulge scales to tighten the MBH - sigma relation with respect to the MBH - Vc relation, as observed.

The AGN-Starburst Connection traced by the Nitrogen Abundance

AbstractThe connection between the active galactic nuclei (AGNs) and star formation activity is one of the most important issues in understanding the coevolution of supermassive black holes (SMBHs) and galaxies. In our recent study, by using SDSS quasar spectra we found that the emission-line flux rations involving a nitrogen line correlate with the Eddington ratio. This correlation suggests that the mass accretion in quasars is associated with a post-starburst phase, when AGB stars enrich the interstellar medium with the nitrogen. Moreover, we focused on nitrogen-loud quasars, that have prominent emission lines of the nitrogen, to investigate whether this argument is correct. In this symposium, we present our results regarding the relation between the star formation and feeding to SMBHs traced by the nitrogen abundance.

The Co-Evolution of Supermassive Black Holes and Galaxies: Observational Constraints

AbstractThe connection between the growth of supermassive black holes (SMBHs) and the assembly of their host galaxies is termed ‘co-evolution’. Understanding co-evolution is one of the most fundamental issues in modern astrophysics. In this contribution, we review recent progress in addressing how the growth of SMBHs is linked to the properties of their host galaxies in the context of galaxy evolution, from the observational point of view. Although a coherent picture has not yet emerged, multiple pathways of co-evolution appear to be favored with a probable dependence on AGN luminosity and redshift.

A TALE OF TWO POPULATIONS: THE CONTRIBUTION OF MERGER AND SECULAR PROCESSES TO THE EVOLUTION OF ACTIVE GALACTIC NUCLEI

Due to the co-evolution of supermassive black holes and their host galaxies, understanding the mechanisms that trigger active galactic nuclei (AGN) are imperative to understanding galaxy evolution and the formation of massive galaxies. It is observationally difficult to determine the trigger of a given AGN due to the difference between the AGN lifetime and triggering timescales. Here, we utilize AGN population synthesis modeling to determine the importance of different AGN triggering mechanisms. An AGN population model is computed by combining an observationally motivated AGN triggering rate and a theoretical AGN light curve. The free parameters of the AGN light curve are constrained by minimizing a \chi squared test with respect to the observed AGN hard X-ray luminosity function. The observed black hole space density, AGN number counts, and X-ray background spectrum are also considered as observational constraints. It is found that major mergers are not able to account for the entire AGN population. Therefore, non-merger processes, such as secular mechanisms, must also trigger AGN. Indeed, non-merger processes are the dominant AGN triggering mechanism at z \lesssim 1--1.5. Furthermore, the shape and evolution of the black hole mass function of AGN triggered by major mergers is intrinsically different from the shape and evolution of the black hole mass function of AGN triggered by secular processes.

A Method of Identifying AGNs Based on Emission-Line Excess and the Nature of Low-Luminosity AGNs in the Sloan Digital Sky Survey. II. The Nature of Low-Luminosity AGNs

Abstract We have developed a new method of identifying active galactic nuclei (AGNs) and studied the nature of low-luminosity AGNs in the Sloan Digital Sky Survey. This is the latter part of a series of papers in which we consider correlations between the AGN activities and the host-galaxy properties. Based on a sample of AGNs identified by a new method developed in the former part (2012, PASJ, 64, 36), we found that AGNs typically show extinction of $\tau_V$$=$ 1.2, and exhibit a wide range of ionization levels. The finding of ionization levels motivated us to use [O II] $+$ [O III] as an indicator of AGN power. We found that AGNs are preferentially located in massive, red, early-type galaxies. Taking into account a selection bias of the Oxygen-excess method, we showed that strong AGNs are located in active star-forming galaxies, and that rapidly growing super-massive black holes are located in rapidly growing galaxies, which clearly shows the coevolution of super-massive black holes and their host galaxies. This is a surprising phenomenon, given that the growths of black holes and host galaxies occur on their respective physical scales which are very different. Interestingly, the AGN power does not strongly correlate with the host-galaxy mass. It seems that the mass works as a ``switch'' for activating AGNs. The absence of AGNs in low-mass galaxies might be due to the absence of super-massive black holes there, but a dedicated observation of the nuclear region of nearby low-mass galaxies would be necessary to obtain a deeper insight into it.

Dynamics Induced by the Central Supermassive Black Holes in Galaxies

Dynamics Induced by the Central Supermassive Black Holes in Galaxies

One-zone models for spheroidal galaxies with a central supermassive black-hole

By means of a one-zone evolutionary model, we study the co-evolution of supermassive black holes and their host galaxies, as a function of the accretion radiative efficiency, dark matter content, and cosmological infall of gas. In particular, the radiation feedback is computed by using the self-regulated Bondi accretion. The models are characterized by strong oscillations when the galaxy is in the AGN state with a high accretion luminosity. We found that these one-zone models are able to reproduce two important phases of galaxy evolution, namely an obscured-cold phase when the bulk of star formation and black hole accretion occur, and the following quiescent hot phase in which accretion remains highly sub-Eddington. A Compton-thick phase is also found in almost all models, associated with the cold phase. An exploration of the parameter space reveals that the closest agreement with the present-day Magorrian relation is obtained, independently of the dark matter halo mass, for galaxies with a low-mass seed black hole, and the accretion radiative efficiency ~0.1.

On the co-evolution of supermassive black holes and their host galaxies since z= 3

[Abridged] To investigate the evolution in the relation between galaxy stellar and central black hole mass we construct a volume limited complete sample of 85 AGN with host galaxy stellar masses M_{*} > 10^{10.5} M_{sol}, and specific X-ray luminosities L_{X} > 2.35 x 10^{43} erg s^{-1} at 0.4 < z < 3. We calculate the Eddington limiting masses of the supermassive black holes residing at the centre of these galaxies, and observe an increase in the average Eddington limiting black hole mass with redshift. By assuming that there is no evolution in the Eddington ratio (\mu) and then that there is maximum possible evolution to the Eddington limit, we quantify the maximum possible evolution in the M_{*} / M_{BH} ratio as lying in the range 700 < M_{*}/M_{BH} < 10000, compared with the local value of M_{*}/M_{BH} ~ 1000. We furthermore find that the fraction of galaxies which are AGN (with L_{X} > 2.35 x 10^{43} erg s^{-1}) rises with redshift from 1.2 +/- 0.2 % at z = 0.7 to 7.4 +/- 2.0 % at z = 2.5. We use our results to calculate the maximum timescales for which our sample of AGN can continue to accrete at their observed rates before surpassing the local galaxy-black hole mass relation. We use these timescales to calculate the total fraction of massive galaxies which will be active (with L_{X} > 2.35 x 10^{43} erg s^{-1}) since z = 3, finding that at least ~ 40% of all massive galaxies will be Seyfert luminosity AGN or brighter during this epoch. Further, we calculate the energy density due to AGN activity in the Universe as 1.0 (+/- 0.3) x 10^{57} erg Mpc^{-3} Gyr^{-1}, potentially providing a significant source of energy for AGN feedback on star formation. We also use this method to compute the evolution in the X-ray luminosity density of AGN with redshift, finding that massive galaxy Seyfert luminosity AGN are the dominant source of X-ray emission in the Universe at z < 3.

Grand unification of AGN activity in the ΛCDM cosmology

We track the co-evolution of supermassive black holes (SMBHs) and their host galaxies. The calculation is embedded in the GALFORM semi-analytical model which simulates the formation and evolution of galaxies in a cold dark matter (CDM) universe. During the evolution of the host galaxy, hot and cold gas are added to the SMBH by flows triggered by halo gas cooling, disc instabilities and galaxy mergers. This builds up the mass and spin of the BH, and the resulting accretion power regulates the gas cooling and subsequent star formation. The accretion flow is assumed to form a geometrically thin cool disc when the accretion rate exceeds 0.01\dot{M}_Edd, and a geometrically thick, radiatively inefficient hot flow when the accretion rate falls below this value. The resulting quasar optical luminosity function matches observations very well, and the mass of the SMBH correlates with the mass of the galaxy bulge as observed. The BH spin distribution depends strongly on whether the gas in any given accretion episode remains in the same plane (prolonged accretion) or whether, due to self-gravity, it fragments into multiple, randomly aligned accretion episodes (chaotic accretion). In the chaotic accretion model there is a clear correlation of spin with SMBH mass. Massive BHs (M>5\times10^8\Msun) are hosted by giant elliptical galaxies and are rapidly spinning, while lower mass BHs are hosted in spiral galaxies and have much lower spin. Using the Blandford-Znajek mechanism for jet production to calculate the jet power, our model is able to reproduce the radio loudness of radio galaxies, LINERS and Seyferts. This is the first confirmation that a CDM galaxy formation model can reproduce the observed phenomenology of AGN.

Large Dynamic Range Simulations of Galaxies Hosting Supermassive Black Holes

AbstractThe co-evolution of supermassive black holes (SMBHs) and their host galaxies is a rich problem, spanning a large-dynamic range and depending on many physical processes. Simulating the transport of gas and angular momentum from super-galactic scales all the way down to the outer edge of the black hole's accretion disk requires sophisticated numerical techniques with extensive treatment of baryonic physics. We use a hydrodynamic adaptive mesh refinement simulation to follow the growth and evolution of a typical disk galaxy hosting an SMBH, in a cosmological context (covering a dynamical range of 10 million!). We have adopted a piecemeal approach, focusing our attention on the gas dynamics in the central few hundred parsecs of the simulated galaxy (with boundary conditions provided by the larger cosmological simulation), and beginning with a simplified picture (no mergers or feedback). In this scenario, we find that the circumnuclear disk remains marginally stable against catastrophic fragmentation, allowing stochastic fueling of gas into the vicinity of the SMBH. I will discuss the successes and the limitations of these simulations, and their future direction.

Co-Evolution of Supermassive Black Holes and Their Host Galaxies

AbstractWe study the evolution of the MBH/Mhost relation up to z = 3 for a sample of 96 quasars with known host galaxy luminosities. Black hole masses are estimated assuming virial equilibrium in the broad-line regions, while the host galaxy masses are inferred from their luminosities. With this data, we are able to pin down the evolution of the MBH/Mhost relation over 85% of the age of the universe. While the MBH/Lhost relation remains nearly unchanged, taking into account the aging of the stellar population, we find that the MBH/Mhost ratio (Γ) increases by a factor ~ 7 from z = 0 to z = 3. We show that the evolution of Γ is independent of radio loudness and quasar luminosity. We propose that the most massive black holes, in their quasar phase at high-redshift, become extremely rare objects in host galaxies of similar mass in the local universe.

The spatial distribution of X-ray selected AGN in theChandradeep fields: a theoretical perspective

We study the spatial distribution of X-ray selected AGN in the framework of hierarchical co-evolution of supermassive black holes and their host galaxies and dark matter haloes. To this end, we have applied the model developed by Croton et al.(2006), De Lucia & Blaizot(2007) and Marulli et al.(2008) to the output of the Millennium Run and obtained hundreds of realizations of past light-cones from which we have extracted realistic mock AGN catalogues that mimic the Chandra deep fields. We find that the model AGN number counts are in fair agreement with observations, except at fluxes <1e-15 erg/cm^2/s. The spatial two-point correlation function predicted by the model is well described by a power-law relation out to 20 Mpc/h, in close agreement with observations. Our model matches the correlation length r_0 of AGN in the Chandra Deep Field North but underestimates it in the Chandra Deep Field South. When fixing the slope to gamma = 1.4, as in Gilli et al. (2005), the statistical significance of the mismatch is 2-2.5 sigma, suggesting that the predicted cosmic variance, which dominates the error budget, may not account for the different correlation length of the AGN in the two fields. While our results are robust to changes in the model prescriptions for the AGN lightcurves, the luminosity dependence of the clustering is sensitive to the different lightcurve models adopted. However, irrespective of the model considered, the luminosity dependence of the AGN clustering in our mock fields seems to be weaker than in the real Chandra fields. The significance of this mismatch needs to be confirmed using larger datasets.

The accuracy of supermassive black hole masses determined by the single-epoch spectrum (Dibai) method

The first set of supermassive black hole mass estimates, published from 1980 to 1984 by E. A. Dibai, are shown to be in excellent agreement with recent reverberation-mapping estimates. Comparison of the masses of 17 AGNs covering a mass range from about 10^6 to 10^9 solar masses shows that the Dibai mass estimates agree with reverberation-mapping mass estimates to significantly better than 0.3 dex and were, on average, only 0.14 dex (~ 40%) systematically lower than masses obtained from reverberation mapping. This surprising agreement with the results of over a quarter of a century ago has important implication for the structure and kinematics of AGNs and implies that AGNs are very similar. Our results give strong support to the use of the single-epoch-spectrum (Dibai) method for investigating the co-evolution of supermassive black holes and their host galaxies.

UPPER LIMITS ON THE MASSES OF 105 SUPERMASSIVE BLACK HOLES FROMHUBBLE SPACE TELESCOPE/SPACE TELESCOPE IMAGING SPECTROGRAPH ARCHIVAL DATA

Based on the modeling of the central emission-line width measured over sub-arcsecond apertures with the Hubble Space Telescope, we present stringent upper bounds on the mass of the central supermassive black hole, MBH, for a sample of 105 nearby galaxies (D<100Mpc) spanning a wide range of Hubble types (E-Sc) and values of the central stellar velocity dispersion, sigma (58-419km/s). For the vast majority of the objects the derived MBH upper limits run parallel and above the well-known MBH-sigma relation independently of the galaxy distance, suggesting that our nebular line-width measurements trace rather well the nuclear gravitational potential. For values of sigma between 90 and 220km/s the 68% of our upper limits falls immediately above the MBH-sigma relation without exceeding the expected MBH values by more than a factor 4.1. No systematic trends or offsets are observed in this sigma range as a function of the galaxy Hubble type or with respect to the presence of a bar. For 6 of our 12 MBH upper limits with sigma<90km/s our line-width measurements are more sensitive to the stellar contribution to the gravitational potential, either due to the presence of a nuclear stellar cluster or because of a greater distance compared to the other galaxies at the low-sigma end of the MBH-sigma relation. Conversely, our MBH upper bounds appear to lie closer to the expected MBH in the most massive elliptical galaxies with values of sigma above 220km/s. Such a flattening of the MBH-sigma relation at its high-sigma end would appear consistent with a coevolution of supermassive black holes and galaxies driven by dry mergers, although better and more consistent measurements for sigma and K-band luminosity are needed for these kind of objects before systematic effects can be ruled out.

Modelling the cosmological co-evolution of supermassive black holes and galaxies - II. The clustering of quasars and their dark environment

We use semi-analytic modeling on top of the Millennium simulation to study the joint formation of galaxies and their embedded supermassive black holes. Our goal is to test scenarios in which black hole accretion and quasar activity are triggered by galaxy mergers, and to constrain different models for the lightcurves associated with individual quasar events. In the present work we focus on studying the spatial distribution of simulated quasars. At all luminosities, we find that the simulated quasar two-point correlation function is fit well by a single power-law in the range 0.5 < r < 20 h^{-1} Mpc, but its normalization is a strong function of redshift. When we select only quasars with luminosities within the range typically accessible by today's quasar surveys, their clustering strength depends only weakly on luminosity, in agreement with observations. This holds independently of the assumed lightcurve model, since bright quasars are black holes accreting close to the Eddington limit, and are hosted by dark matter haloes with a narrow mass range of a few 10^12 h^{-1} M_sun. Therefore the clustering of bright quasars cannot be used to disentangle lightcurve models, but such a discrimination would become possible if the observational samples can be pushed to significantly fainter limits. Overall, our clustering results for the simulated quasar population agree rather well with observations, lending support to the conjecture that galaxy mergers could be the main physical process responsible for triggering black hole accretion and quasar activity.

Coevolution of Supermassive Black Holes and Circumnuclear Disks

We propose a new evolutionary model of a supermassive black hole (SMBH) and a circumnuclear disk (CND), taking into account the mass-supply from a host galaxy and the physical states of CND. In the model, two distinct accretion modes depending on gravitational stability of the CND play a key role on accreting gas to a SMBH. (i) If the CMD is gravitationally unstable, energy feedback from supernovae (SNe) supports a geometrically thick, turbulent gas disk. The accretion in this mode is dominated by turbulent viscosity, and it is significantly larger than that in the mode (ii), i.e., the CMD is supported by gas pressure. Once the gas supply from the host is stopped, the high accretion phase ($\sim 0.01- 0.1 M_{\odot} {\rm yr}^{-1}$) changes to the low one (mode (ii), $\sim 10^{-4} M_{\odot} {\rm yr}^{-1}$), but there is a delay with $\sim 10^{8}$ yr. Through this evolution, the gas-rich CND turns into the gas poor stellar disk. We found that not all the gas supplied from the host galaxy accrete onto the SMBH even in the high accretion phase (mode (i)), because the part of gas is used to form stars. As a result, the final SMBH mass ($M_{\rm BH,final}$) is not proportional to the total gas mass supplied from the host galaxy ($M_{\rm sup}$); $M_{\rm BH,final}/M_{\rm sup}$ decreases with $M_{\rm sup}$.This would indicate that it is difficult to form a SMBH with $\sim 10^{9} M_{\odot}$ observed at high-$z$ QSOs. The evolution of the SMBH and CND would be related to the evolutionary tracks of different type of AGNs.