We review the observed demographics and inferred evolution of supermassive black holes (BHs) found by dynamical modeling of spatially resolved kinematics. Most influential was the discovery of a tight correlation between BH mass and the velocity dispersion of the host-galaxy bulge. It and other correlations led to the belief that BHs and bulges coevolve by regulating each other's growth. New results are now replacing this simple story with a richer and more plausible picture in which BHs correlate differently with different galaxy components. BHs are found in pure-disk galaxies, so classical (elliptical-galaxy-like) bulges are not necessary to grow BHs. But BHs do not correlate with galaxy disks. And any correlations with disk-grown pseudobulges or halo dark matter are so weak as to imply no close coevolution. We suggest that there are four regimes of BH feedback. 1- Local, stochastic feeding of small BHs in mainly bulgeless galaxies involves too little energy to result in coevolution. 2- Global feeding in major, wet galaxy mergers grows giant BHs in short, quasar-like "AGN" events whose feedback does affect galaxies. This makes classical bulges and coreless-rotating ellipticals. 3- At the highest BH masses, maintenance-mode feedback into X-ray gas has the negative effect of helping to keep baryons locked up in hot gas. This happens in giant, core-nonrotating ellipticals. They inherit coevolution magic from smaller progenitors. 4- Independent of any feedback physics, the averaging that results from successive mergers helps to engineer tight BH correlations.Comment: 136 pages, 38 postscript figures, 4 tables; requires cittable.tex, psfig.tex, annrev4K-E.tex; accepted for publication in Volume 51 (2013) of Annual Review of Astronomy and Astrophysics; Supplementary Information will be submitted to arXiv separately in approximately 2013 Jun

Observations of nearby galaxies reveal a strong correlation between the mass of the central dark object M BH and the velocity dispersion of the host galaxy, of the form logðM BH =M Þ ¼ þ logð = 0 Þ; however, published estimates of the slope span a wide range (3.75-5.3). Merritt & Ferrarese have argued that low slopes (d4) arise because of neglect of random measurement errors in the dispersions and an incorrect choice for the dispersion of the Milky Way Galaxy. We show that these explanations and several others account for at most a small part of the slope range. Instead, the range of slopes arises mostly because of systematic differences in the velocity dispersions used by different groups for the same galaxies. The origin of these differences remains unclear, but we suggest that one significant component of the difference results from Ferrarese & Merritt's extrapolation of central velocity dispersions to r e =8 (r e is the effective radius) using an empirical formula. Another component may arise from dispersion-dependent systematic errors in the measurements. A new determination of the slope using 31 galaxies yields ¼ 4:02 AE 0:32, ¼ 8:13 AE 0:06 for 0 ¼ 200 km s À1 . The M BH -relation has an intrinsic dispersion in log M BH that is no larger than 0.25-0.3 dex and may be smaller if observational errors have been underestimated. In an appendix, we present a simple kinematic model for the velocity-dispersion profile of the Galactic bulge.

We describe a correlation between the mass of a galaxy's central black hole and the luminosity-weighted M bh line-of-sight velocity dispersion within the half-light radius. The result is based on a sample of 26 galaxies, j e including 13 galaxies with new determinations of black hole masses from errors. The -relation is of interest not only for its strong predictive power but also because it implies that M j bh e central black hole mass is constrained by and closely related to properties of the host galaxy's bulge.

We present a two-dimensional (2-D) fitting algorithm (GALFIT) designed to extract structural components from galaxy images, with emphasis on closely modeling light profiles of spatially well-resolved, nearby galaxies observed with the Hubble Space Telescope. Our algorithm improves on previous techniques in two areas, by being able to simultaneously fit a galaxy with an arbitrary number of components, and with optimization in computation speed, suited for working on large galaxy images. We use 2-D models such as the ``Nuker'' law, the Sersic (de Vaucouleurs) profile, an exponential disk, and Gaussian or Moffat functions. The azimuthal shapes are generalized ellipses that can fit disky and boxy components. Many galaxies with complex isophotes, ellipticity changes, and position-angle twists can be modeled accurately in 2-D. When examined in detail, we find that even simple-looking galaxies generally require at least three components to be modeled accurately, rather than the one or two components more often employed. We illustrate this by way of 7 case studies, which include regular and barred spiral galaxies, highly disky lenticular galaxies, and elliptical galaxies displaying various levels of complexities. A useful extension of this algorithm is to accurately extract nuclear point sources in galaxies. We compare 2-D and 1-D extraction techniques on simulated images of galaxies having nuclear slopes with different degrees of cuspiness, and we then illustrate the application of the program to several examples of nearby galaxies with weak nuclei.Comment: 29 pages, 14 figures, abridged version. Full version AJ accepted. For full version with high resolution figures, go to: http://zwicky.as.arizona.edu/~cyp/work/galfit.ps.g

We present a two-dimensional (2-D) fitting algorithm (Galfit, Version 3) with new capabilities to study the structural components of galaxies and other astronomical objects in digital images. Our technique improves on previous 2-D fitting algorithms by allowing for irregular, curved, logarithmic and power-law spirals, ring and truncated shapes in otherwise traditional parametric functions like the Sérsic, Moffat, King, Ferrer, etc., profiles. One can mix and match these new shape features freely, with or without constraints, apply them to an arbitrary number of model components and of numerous profile types, so as to produce realistic-looking galaxy model images. Yet, despite the potential for extreme complexity, the meaning of the key parameters like the Sérsic index, effective radius or luminosity remain intuitive and essentially unchanged. The new features have an interesting potential for use to quantify the degree of asymmetry of galaxies, to quantify low surface brightness tidal features beneath and beyond luminous galaxies, to allow more realistic decompositions of galaxy subcomponents in the presence of strong rings and spiral arms, and to enable ways to gauge the uncertainties when decomposing galaxy subcomponents. We illustrate these new features by way of several case studies that display various levels of complexity.

When surrounded by a transparent emission region, black holes are expected to reveal a dark shadow caused by gravitational light bending and photon capture at the event horizon. To image and study this phenomenon, we have assembled the Event Horizon Telescope, a global very long baseline interferometry array observing at a wavelength of 1.3 mm. This allows us to reconstruct event-horizon-scale images of the supermassive black hole candidate in the center of the giant elliptical galaxy M87. We have resolved the central compact radio source as an asymmetric bright emission ring with a diameter of 42±3 μas, which is circular and encompasses a central depression in brightness with a flux ratio 10:1. The emission ring is recovered using different calibration and imaging schemes, with its diameter and width remaining stable over four different observations carried out in different days. Overall, the observed image is consistent with expectations for the shadow of a Kerr black hole as predicted by general relativity. The asymmetry in brightness in the ring can be explained in terms of relativistic beaming of the emission from a plasma rotating close to the speed of light around a black hole. We compare our images to an extensive library of ray-traced general-relativistic magnetohydrodynamic simulations of black holes and derive a central mass of M=(6.5±0.7)×10 9 M e . Our radiowave observations thus provide powerful evidence for the presence of supermassive black holes in centers of galaxies and as the central engines of active galactic nuclei. They also present a new tool to explore gravity in its most extreme limit and on a mass scale that was so far not accessible.

A significant fraction of nearby galaxies show evidence of weak nuclear activity unrelated to normal stellar processes. Recent high-resolution, multiwavelength observations indicate that the bulk of this activity derives from black hole accretion with a wide range of accretion rates. The low accretion rates that typify most low-luminosity active galactic nuclei induce significant modifications to their central engine. The broad-line region and obscuring torus disappear in some of the faintest sources, and the optically thick accretion disk transforms into a three-component structure consisting of an inner radiatively inefficient accretion flow, a truncated outer thin disk, and a jet or outflow. The local census of nuclear activity supports the notion that most, perhaps all, bulges host a central supermassive black hole, although the existence of active nuclei in at least some late-type galaxies suggests that a classical bulge is not a prerequisite to seed a nuclear black hole.

It has been established that virial masses for black holes in low-redshift active galaxies can be estimated from measurements of the optical continuum strength and the width of the broad H line. Under various circumstances, however, both of these quantities can be challenging to measure or can be subject to large systematic uncertainties. To mitigate these difficulties, we present a new method for estimating black hole masses. From analysis of a new sample of broad-line active galactic nuclei, we find that H luminosity scales almost linearly with optical continuum luminosity and that a strong correlation exists between H and H line widths. These two empirical correlations allow us to translate the standard virial mass system to a new one based solely on observations of the broad H emission line.