Evolution Of Supermassive Black Holes Research Articles

Supermassive black holes (SMBH) at work: M87, a case study of the effects of SMBH outbursts

AbstractSupermassive black holes (SMBHs) play key roles in galaxy and cluster evolution. This is most clearly seen in the “fossil record” that is imprinted in the gas rich atmospheres of early type galaxies, groups, and clusters by powerful SMBH outbursts. From a detailed X-ray study of M87, we present the properties of a typical SMBH outburst, its evolution, and the energy partition between shocks and the enthalpy of the gas cavities inflated by the SMBH. About 12 Myr ago, the SMBH in M87 inflated a cavity of relativistic plasma which is still centered near the galaxy nucleus. This outburst drove a shock into the surrounding gas. For M87, we show that the outburst duration is a few Myr and that about 50% of the total energy (5 × 1057 ergs) resides in the bubble inflated by the jet from the SMBH, that 25% of the outburst energy is deposited directly into the ambient atmosphere by the shock, and that 25% of the outburst energy is lost from the radiatively bright core as the weak shock moves to large radii. We conclude by describing a future X-ray mission, SMART-X, with < 1” angular resolution that would allow us to study the evolution of SMBHs and the hot, X-ray emitting atmospheres from high redshifts to the present for M87-like systems.

Resonant motions of supermassive black hole triples

AbstractTriple supermassive black holes (SMBH) can form during the hierarchical mergers of massive galaxies with an existing binary. Perturbations by a third black hole may accelerate the merging process of an inner binary, for example through the Kozai mechanism. We analyze the evolution of simulated hierarchical triple SMBHs in galactic centers, and find resonances in the evolution of the semi-major axis, the eccentricity and the inclination, for both the inner and the outer orbits of the triple system, which are not only Kozai like. Through resonant oscillations, SMBH can trigger a significant increase of the inner SMBH binary eccentricity shortening the merger timescale expected from gravitational wave (GW) emission. As hierarchical triple SMBHs may be frequent in massive galaxies, the influence of orbital resonances is of great importance to our understanding of black hole coalescence and gravitational wave detection. Although Kozai mechanism is believed to play an important role in this process, detailed studies on the pattern of these resonances is necessary.

Evolution of binary supermassive black holes and the final-parsec problem

AbstractI review the evolution of binary supermassive black holes and focus on the stellar-dynamical mechanisms that may help to overcome the final-parsec problem – the possible stalling of the binary at a separation much larger than is required for an efficient gravitational wave emission. RecentN-body simulations have suggested that a departure from spherical symmetry in the nucleus of the galaxy may keep the rate of interaction of stars with the binary at a high enough level so that the binary continues to shrink rather rapidly. However, a major problem of all these simulations is that they do not probe the regime where collisionless effects are dominant – in other words, the number of particles in the simulation is still not sufficient to reach the asymptotic behavior of the system. I present a novel Monte Carlo method for simulating both collisional and collisionless evolution of non-spherical stellar systems, and apply it for the problem of binary supermassive black hole evolution. I show that in triaxial galaxies the final-parsec problem is largely non-existent, while in the axisymmetric case it seems to still exist in the limit of purely collisionless regime relevant for real galaxies, but disappears in theN-body simulations where the feasible values ofNare still too low to get rid of collisional effects.

THE BLACK HOLE MASS FUNCTION DERIVED FROM LOCAL SPIRAL GALAXIES

We present our determination of the nuclear supermassive black hole (SMBH) mass function for spiral galaxies in the Local Universe, established from a volume-limited sample consisting of a statistically complete collection of the brightest spiral galaxies in the Southern Hemisphere. Our SMBH mass function agrees well at the high-mass end with previous values given in the literature. At the low-mass end, inconsistencies exist in previous works that still need to be resolved, but our work is more in line with expectations based on modeling of SMBH evolution. This low-mass end of the spectrum is critical to our understanding of the mass function and evolution of SMBHs since the epoch of maximum quasar activity. A luminosity distance $\leq$ 25.4 $Mpc$ and an absolute B-band magnitude $\leq$ -19.12 define the sample. These limits define a sample of 140 spiral galaxies, with 128 measurable pitch angles to establish the pitch angle distribution for this sample. This pitch angle distribution function may be useful in the study of the morphology of late-type galaxies. We then use an established relationship between the pitch angle and the mass of the central SMBH in a host galaxy in order to estimate the mass of the 128 respective SMBHs in this sample. This result effectively gives us the distribution of mass for SMBHs residing in spiral galaxies over a lookback time $\leq$ 82.1 $h_{67.77}^{-1}$ $Myr$ and contained within a comoving volume of 3.37 $\times$ $10^4$ $h_{67.77}^{-3}$ $Mpc^3$. We estimate the density of SMBHs residing in spiral galaxies in the Local Universe is $5.54_{-2.73}^{+6.55}$ $\times$ $10^4$ $h_{67.77}^3$ $M_{\odot}$ $Mpc^{-3}$. Thus, our derived cosmological SMBH mass density for spiral galaxies is $\Omega_{BH} = 4.35_{-2.15}^{+5.14}$ $\times$ $10^{-7}$ $h_{67.77}$.

Binary supermassive black hole environments diminish the gravitational wave signal in the pulsar timing band

We assess the effects of super-massive black hole (SMBH) environments on the gravitational-wave (GW) signal from binary SMBHs. To date, searches with pulsar timing arrays for GWs from binary SMBHs, in the frequency band $\sim1-100$\,nHz, include the assumptions that all binaries are circular and evolve only through GW emission. However, dynamical studies have shown that the only way that binary SMBH orbits can decay to separations where GW emission dominates the evolution is through interactions with their environments. We augment an existing galaxy and SMBH formation and evolution model with calculations of binary SMBH evolution in stellar environments, accounting for non-zero binary eccentricities. We find that coupling between binaries and their environments causes the expected GW spectral energy distribution to be reduced with respect to the standard assumption of circular, GW-driven binaries, for frequencies up to $\sim20$\,nHz. Larger eccentricities at binary formation further reduce the signal in this regime. We also find that GW bursts from individual eccentric binary SMBHs are unlikely to be detectable with current pulsar timing arrays. The uncertainties in these predictions are large, owing to observational uncertainty in SMBH-galaxy scaling relations and the galaxy stellar mass function, uncertainty in the nature of binary-environment coupling, and uncertainty in the numbers of the most massive binary SMBHs. We conclude, however, that low-frequency GWs from binary SMBHs may be more difficult to detect with pulsar timing arrays than currently thought.

EROSITA – Mapping the X‐ray universe

AbstracteROSITA (extended ROentgen Survey with an Imaging Telescope Array) is the core instrument on the Russian/German Spektrum‐Roentgen‐Gamma (SRG) mission which is currently scheduled for launch late 2015/early 2016. eROSITA will perform a deep survey of the entire X‐ray sky. In the soft band (0.5–2 keV), it will be about 30 times more sensitive than ROSAT, while in the hard band (2–8 keV) it will provide the first ever true imaging survey of the sky. The design driving science is the detection of large samples of galaxy clusters to redshifts z > 1 in order to study the large scale structure in the universe and test cosmological models including dark energy. In addition, eROSITA is expected to yield a sample of a few million active galactic nuclei, including obscured objects, revolutionizing our view of the evolution of supermassive black holes. The survey will also provide new insights into a wide range of astro‐physical phenomena, including X‐ray binaries, active stars, and diffuse emission within the Galaxy. eROSITA is currently (Jan 2014) in its flight model and calibration phase. All seven flight mirror modules (plus one spare) have been delivered and measured in X‐rays. The first camera including the complete electronics has been extensively tested. So far, all subsystems and components are well within the expected performance. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)

THE INTRINSIC QUASAR LUMINOSITY FUNCTION: ACCOUNTING FOR ACCRETION DISK ANISOTROPY

Quasar luminosity functions are a fundamental probe of the growth and evolution of supermassive black holes. Measuring the intrinsic luminosity function is difficult in practice, due to a multitude of observational and systematic effects. As sample sizes increase and measurement errors drop, characterizing the systematic effects is becoming more important. It is well known that the continuum emission from the accretion disk of quasars is anisotropic --- in part due to its disk-like structure --- but current luminosity function calculations effectively assume isotropy over the range of unobscured lines of sight. Here, we provide the first steps in characterizing the effect of random quasar orientations and simple models of anisotropy on observed luminosity functions. We find that the effect of orientation is not insignificant and exceeds other potential corrections such as those from gravitational lensing of foreground structures. We argue that current observational constraints may overestimate the intrinsic luminosity function by as much as a factor of ~2 on the bright end. This has implications for models of quasars and their role in the Universe, such as quasars' contribution to cosmological backgrounds.

PROSPECTS FOR MEASURING SUPERMASSIVE BLACK HOLE MASSES WITH FUTURE EXTREMELY LARGE TELESCOPES

The next generation of giant-segmented mirror telescopes ($>$ 20 m) will enable us to observe galactic nuclei at much higher angular resolution and sensitivity than ever before. These capabilities will introduce a revolutionary shift in our understanding of the origin and evolution of supermassive black holes by enabling more precise black hole mass measurements in a mass range that is unreachable today. We present simulations and predictions of the observations of nuclei that will be made with the Thirty Meter Telescope (TMT) and the adaptive optics assisted integral-field spectrograph IRIS, which is capable of diffraction-limited spectroscopy from $Z$ band (0.9 $\mu$m) to $K$ band (2.2 $\mu$m). These simulations, for the first time, use realistic values for the sky, telescope, adaptive optics system, and instrument, to determine the expected signal-to-noise ratio of a range of possible targets spanning intermediate mass black holes of $\sim10^4$ \msun to the most massive black holes known today of $>10^{10}$ $M_\odot$. We find that IRIS will be able to observe Milky Way-mass black holes out the distance of the Virgo cluster, and will allow us to observe many more brightest cluster galaxies where the most massive black holes are thought to reside. We also evaluate how well the kinematic moments of the velocity distributions can be constrained at the different spectral resolutions and plate scales designed for IRIS. We find that a spectral resolution of $\sim8000$ will be necessary to measure the masses of intermediate mass black holes. By simulating the observations of galaxies found in SDSS DR7, we find that over $10^5$ massive black holes will be observable at distances between $0.005 < z < 0.18$ with the estimated sensitivity and angular resolution provided by access to $Z$-band (0.9 $\mu$m) spectroscopy from IRIS and the TMT adaptive optics system. (Abridged)

Nuclear coups: dynamics of black holes in galaxy mergers

We study the dynamical evolution of supermassive black holes (BHs) in merging galaxies on scales of hundreds of kpc to 10 pc, to identify the physical processes that aid or hinder the orbital decay of BHs. We present hydrodynamical simulations of galaxy mergers with a resolution of $\leq$20 pc, chosen to accurately track the motion of the nuclei and provide a realistic environment for the evolution of the BHs. We find that, during the late stages of the merger, tidal shocks inject energy in the nuclei, causing one or both nuclei to be disrupted and leaving their BH `naked', without any bound gas or stars. In many cases, the nucleus that is ultimately disrupted is that of the larger galaxy (`nuclear coup'), as star formation grows a denser nuclear cusp in the smaller galaxy. We supplement our simulations with an analytical estimate of the orbital-decay time required for the BHs to form a binary at unresolved scales, due to dynamical friction. We find that, when a nuclear coup occurs, the time-scale is much shorter than when the secondary's nucleus is disrupted, as the infalling BH is more massive, and it also finds itself in a denser stellar environment.

Star formation and accretion in the circumnuclear disks of active galaxies

We explore the evolution of supermassive black holes (SMBH) centered in a circumnuclear disk (CND) as a function of the mass supply from the host galaxy and considering different star formation laws, which may give rise to a self-regulation via the injection of supernova-driven turbulence. A system of equations describing star formation, black hole accretion and angular momentum transport was solved for an axisymmetric disk in which the gravitational potential includes contributions from the black hole, the disk and the hosting galaxy. Our model extends the framework provided by Kawakatu et al. (2008) by separately considering the inner and outer part of the disk, and by introducing a potentially non-linear dependence of the star formation rate on the gas surface density and the turbulent velocity. The star formation recipes are calibrated using observational data for NGC 1097, while the accretion model is based on turbulent viscosity as a source of angular momentum transport in a thin viscous accretion disk. We find that current data provide no strong constraint on the star formation recipe, and can in particular not distinguish between models entirely regulated by the surface density, and models including a dependence on the turbulent velocity. The evolution of the black hole mass, on the other hand, strongly depends on the applied star formation law, as well as the mass supply from the host galaxy. We suggest to explore the star formation process in local AGN with high-resolution ALMA observations to break the degeneracy between different star formation models.

The discovery of eight z ~ 6 quasars from Pan-STARRS1

AbstractHigh-redshift quasars are unique probes of the evolution of supermassive black holes and the intergalactic medium at the end of the epoch of reionization. We present the optical spectra of eight new z ~ 6 quasars selected from the Panoramic Survey Telescope &amp; Rapid Response System 1 (Pan-STARRS1). Details of the selection strategy can be found in Bañados et al. (2014). With this work we increase the number of known quasars at z &lt; 5.7 by more than 10%. The quasars discovered here span a large range of luminosities (19.6 ≤ zP1 ≤ 21.2) and are remarkably heterogeneous in their spectral features: half of them show bright emission lines whereas the other half show weak or no Lyα emission line. We find a larger fraction of weak–line emission quasars than in lower redshift studies, although still based on low number statistics, this may imply that the quasar population could be more diverse than previously thought.

Obscured quasars at redshift z ~ 2

AbstractCompleting the census of AGN in the Universe is the key to understanding the cosmic evolution of supermassive black holes (SMBH) and galaxies, and to resolving the spectrum of the X-ray background (XRB). However, a large population of AGN, especially the heavily obscured, Compton-thick AGN, are still missing from even the deepest X-ray surveys. The infrared spectra energy distribution (SED) of distant star-forming galaxies can reveal the presence of bright AGN activity. Using some of the deepest infrared, X-ray and radio data available in the GOODS fields, we identify a population of infrared bright quasars at redshift z ~ 2, which are often missed in the X-ray band. Amongst these sources the number of obscured and heavily-obscured quasars is much higher than those previously found in several X-ray and optical selected samples. A unique view on these heavily-obscured quasars is now given at high energies by NuSTAR. I will present the first NuSTAR detection of a heavily obscured quasar at z 2. This source is a potential archetype of the heavily-obscured high-z AGN in which most of the black hole growth is happening, that can explain the mysterious missing fraction of the XRB.

Mapping the active Universe with eROSITA

AbstracteROSITA (extended Röntgen Survey with an Imaging Telescope Array) is the core instrument on the Russian Spektrum-Röntgen-Gamma (SRG) mission which is current scheduled for launch in Q4 2014. eROSITA will perform a deep survey of the entire X-ray sky. In the soft band (0.5–2 keV), it will be about 30 times more sensitive than ROSAT, while in the hard band (2–8 keV) it will provide the first ever true imaging survey of the sky. The design driving science is the detection of large samples of galaxy clusters up to redshifts z ~ 1, in order to study the large scale structure in the Universe and test cosmological models including Dark Energy. In addition, eROSITA is expected to yield a sample of about 3 million active galactic nuclei, which is bound to revolutionize our view of the evolution of supermassive black holes and their impact on the process of structure formation in the Universe. The survey will also provide new insights into a wide range of astrophysical phenomena, including isolated Neutron Stars and Black Holes, X-ray binaries, active stars and diffuse emission within the Galaxy, as well as more exotic ones such as gamma-ray bursts, tidal disruption of stars in galactic nuclei and binary black holes. In this talk I presented the main characteristics of the mission and focus on the scientific drivers for extragalactic all-sky surveys of AGN. All what was presented at the Symposium (plots, simulations, expected numbers of various kind of sources –QSO, obscured and CT AGN– their properties and evolution with redshift) can be found in the official eROSITA Science Book (Merloni et al., 2012).

SUPERMASSIVE BLACK HOLE BINARY EVOLUTION IN AXISYMMETRIC GALAXIES: THE FINAL PARSEC PROBLEM IS NOT A PROBLEM

During a galaxy merger, the supermassive black hole (SMBH) in each galaxy is thought to sink to the center of the potential and form a supermassive black hole binary; this binary can eject stars via 3-body scattering, bringing the SMBHs ever closer. In a static spherical galaxy model, the binary stalls at a separation of about a parsec after ejecting all the stars in its loss cone -- this is the well-known final parsec problem. Earlier work has shown that the centrophilic orbits in triaxial galaxy models are key in refilling the loss cone at a high enough rate to prevent the black holes from stalling. However, the evolution of binary SMBHs has never been explored in axisymmetric galaxies, so it is not clear if the final parsec problem persists in these systems. Here we use a suite of direct N-body simulations to follow SMBH binary evolution in galaxy models with a range of ellipticity. For the first time, we show that mere axisymmetry can solve the final parsec problem; we find the the SMBH evolution is independent of N for an axis ratio of c/a=0.8, and that the SMBH binary separation reaches the gravitational radiation regime for c/a=0.75.

DISK ASSEMBLY AND THEMBH-σeRELATION OF SUPERMASSIVE BLACK HOLES

Recent Hubble Space Telescope (HST) observations have revealed that a majority of active galactic nuclei (AGN) at z ~ 1-3 are resident in isolated disk galaxies, contrary to the usual expectation that AGN are triggered by mergers. Here we develop a new test of the cosmic evolution of supermassive black holes (SMBHs) in disk galaxies by considering the local population of SMBHs. We show that substantial SMBH growth in spiral galaxies is required as disks assemble. SMBHs exhibit a tight relation between their mass and the velocity dispersion of the spheroid within which they reside, the M_BH-sigma relation. In disk galaxies the bulge is the spheroid of interest. We explore the evolution of the M_BH-sigma relation when bulges form together with SMBHs on the M_BH-sigma relation and then slowly reform a disk around them. The formation of the disk compresses the bulge raising its sigma. We present evidence for such compression in the form of larger velocity dispersion of classical bulges compared with elliptical galaxies at the same mass. This compression leads to an offset in the M_BH-sigma relation if it is not accompanied by an increased M_BH. We quantify the expected offset based on photometric data and show that, on average, SMBHs must grow by ~ 50-65% just to remain on the M_BH-sigma relation. We find no significant offset in the M_BH-sigma relations of classical bulges and of ellipticals, implying that SMBHs have been growing along with disks. Our simulations demonstrate that SMBH growth is necessary for the local population of disk galaxies to have remained on the M_BH-sigma relation.

Spin evolution of supermassive black holes and galactic nuclei

The spin angular momentum $\mathbit{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 $\mathbit{S}$ and the stellar angular momenta ${\mathbit{L}}_{j}$, $j=1\dots{}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 $\mathbit{S}$ about the total angular momentum vector of the system and (2) damped precession, leading, in less than one precessional period, to alignment of $\mathbit{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 $\mathbit{S}$ in nuclei where changes in the ${\mathbit{L}}_{j}$ 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 cluster's 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 $\ensuremath{\sim}{10}^{7}--{10}^{8}\text{ }\text{ }\mathrm{yr}$ for low-mass SBHs in dense nuclei, $\ensuremath{\sim}{10}^{8}--{10}^{10}\text{ }\text{ }\mathrm{yr}$ for SBH masses $\ensuremath{\sim}{10}^{8}{\mathcal{M}}_{\ensuremath{\bigodot}}$, and $\ensuremath{\sim}{10}^{10}--{10}^{11}\text{ }\text{ }\mathrm{yr}$ 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.

Observing the Dynamics of Supermassive Black Hole Binaries with Pulsar Timing Arrays

Pulsar timing arrays are a prime tool to study unexplored astrophysical regimes with gravitational waves. Here, we show that the detection of gravitational radiation from individually resolvable supermassive black hole binary systems can yield direct information about the masses and spins of the black holes, provided that the gravitational-wave-induced timing fluctuations both at the pulsar and at Earth are detected. This in turn provides a map of the nonlinear dynamics of the gravitational field and a new avenue to tackle open problems in astrophysics connected to the formation and evolution of supermassive black holes. We discuss the potential, the challenges, and the limitations of these observations.

Cosmological Evolution of Supermassive Black Holes: Mass Functions and Spins

AbstractWe derive the mass function of supermassive black holes (SMBHs) over the redshift range 0 &gt; z ≲ 2, using the latest deep luminosity and mass functions of field galaxies. Applying this mass function, combined with the bolometric luminosity function of active galactic nuclei (AGNs), into the the continuity equation of SMBH number density, we explicitly obtain the mass-dependent cosmological evolution of the radiative efficiency for accretion. We suggest that the accretion history of SMBHs and their spins evolve in two distinct regimes: an early phase of prolonged accretion, plausibly driven by major mergers, during which the black hole spins up, then switching to a period of random, episodic accretion, governed by minor mergers and internal secular processes, during which the hole spins down. The transition epoch depends on mass, mirroring other evidence for “cosmic downsizing” in the AGN population.

The Cosmogony of Super-Massive Black Holes

We summarize our recent results on the evolution of super-massive black holes in the cores of galaxies. Our models are based on the viscous, Eddington-limited evolution of self-gravitating accretion disks. We find that, within this framework, one can explain the growth time scales and the luminosities of individual objects as well as the variation of these quantities for AGN with different black hole masses.

CLUSTERING OF OBSCURED AND UNOBSCURED QUASARS IN THE BOÖTES FIELD: PLACING RAPIDLY GROWING BLACK HOLES IN THE COSMIC WEB

We present the first measurement of the spatial clustering of mid-infrared selected obscured and unobscured quasars, using a sample in the redshift range 0.7 < z < 1.8 selected from the 9 deg^2 Bootes multiwavelength survey. Recently the Spitzer Space Telescope and X-ray observations have revealed large populations of obscured quasars that have been inferred from models of the X-ray background and supermassive black hole evolution. To date, little is known about obscured quasar clustering, which allows us to measure the masses of their host dark matter halos and explore their role in the cosmic evolution of black holes and galaxies. In this study we use a sample of 806 mid-infrared selected quasars and ~250,000 galaxies to calculate the projected quasar-galaxy cross-correlation function w_p(R). The observed clustering yields characteristic dark matter halo masses of log (M_halo [h^-1 M_sun]) = 12.7^+0.4_-0.6 and 13.3^+0.3_-0.4 for unobscured quasars (QSO-1s) and obscured quasars (Obs-QSOs), respectively. The results for QSO-1s are in excellent agreement with previous measurements for optically-selected quasars, while we conclude that the Obs-QSOs are at least as strongly clustered as the QSO-1s. We test for the effects of photometric redshift errors on the optically-faint Obs-QSOs, and find that our method yields a robust lower limit on the clustering; photo-z errors may cause us to underestimate the clustering amplitude of the Obs-QSOs by at most ~20%. We compare our results to previous studies, and speculate on physical implications of stronger clustering for obscured quasars.