Mrk 1216 is a nearby, early-type galaxy with a small effective radius of 2.8 kpc and a large stellar velocity dispersion of 308 km/s for its K-band luminosity of 1.4x10^11 L_sun. Using integral-field spectroscopy assisted by adaptive optics from Gemini North, we measure spatially resolved stellar kinematics within ~450 pc of the galaxy nucleus. The galaxy exhibits regular rotation with velocities of \pm 180 km/s and a sharply peaked velocity dispersion profile that reaches 425 km/s at the center. We fit axisymmetric, orbit-based dynamical models to the combination of these high angular resolution kinematics, large-scale kinematics extending to roughly three effective radii, and Hubble Space Telescope imaging, resulting in a constraint of the mass of the central black hole in Mrk 1216. After exploring several possible sources of systematics that commonly affect stellar-dynamical black hole mass measurements, we find a black hole mass of (4.9\pm1.7)x10^9 M_sun and a H-band stellar mass-to-light ratio of 1.3\pm0.4 M_sun/L_sun (1-sigma uncertainties). Mrk 1216 is consistent with the local black hole mass - stellar velocity dispersion relation, but is a factor of ~5-10 larger than expectations from the black hole mass - bulge luminosity and black hole mass - bulge mass correlations when conservatively using the galaxy's total luminosity or stellar mass. This behavior is quite similar to the extensively studied compact galaxy NGC 1277. Resembling the z~2 quiescent galaxies, Mrk 1216 may be a passively evolved descendant, and perhaps reflects a previous era when galaxies contained over-massive black holes relative to their bulge luminosities/masses, and the growth of host galaxies had yet to catch up.

The nearby lenticular galaxy NGC 1277 is thought to host one of the largest black holes known, however the black hole mass measurement is based on low spatial resolution spectroscopy. In this paper, we present Gemini Near-infrared Integral Field Spectrometer observations assisted by adaptive optics. We map out the galaxy's stellar kinematics within ~440 pc of the nucleus with an angular resolution that allows us to probe well within the region where the potential from the black hole dominates. We find that the stellar velocity dispersion rises dramatically, reaching ~550 km/s at the center. Through orbit-based, stellar-dynamical models we obtain a black hole mass of (4.9 \pm 1.6) x 10^9 Msun (1-sigma uncertainties). Although the black hole mass measurement is smaller by a factor of ~3 compared to previous claims based on large-scale kinematics, NGC 1277 does indeed contain one of the most massive black holes detected to date, and the black hole mass is an order of magnitude larger than expectations from the empirical relation between black hole mass and galaxy luminosity. Given the galaxy's similarities to the higher redshift (z~2) massive quiescent galaxies, NGC 1277 could be a relic, passively evolving since that period. A population of local analogs to the higher redshift quiescent galaxies that also contain over-massive black holes may suggest that black hole growth precedes that of the host galaxy.

Located in the Perseus cluster, NGC 1271 is an early-type galaxy with a small effective radius of 2.2 kpc and a large stellar velocity dispersion of 276 km/s for its K-band luminosity of 8.9x10^{10} L_sun. We present a mass measurement for the black hole in this compact, high-dispersion galaxy using observations from the integral field spectrograph NIFS on the Gemini North telescope assisted by laser guide star adaptive optics, large-scale integral field unit observations with PPAK at the Calar Alto Observatory, and Hubble Space Telescope WFC3 imaging observations. We are able to map out the stellar kinematics on small spatial scales, within the black hole sphere of influence, and on large scales that extend out to four times the galaxy's effective radius. We find that the galaxy is rapidly rotating and exhibits a sharp rise in the velocity dispersion. Through the use of orbit-based stellar dynamical models, we determine that the black hole has a mass of (3.0^{+1.0}_{-1.1}) x 10^9 M_sun and the H-band stellar mass-to-light ratio is 1.40^{+0.13}_{-0.11} M_sun/L_sun (1-sigma uncertainties). NGC 1271 occupies the sparsely-populated upper end of the black hole mass distribution, but is very different from the Brightest Cluster Galaxies (BCGs) and giant elliptical galaxies that are expected to host the most massive black holes. Interestingly, the black hole mass is an order of magnitude larger than expectations based on the galaxy's bulge luminosity, but is consistent with the mass predicted using the galaxy's bulge stellar velocity dispersion. More compact, high-dispersion galaxies need to be studied using high spatial resolution observations to securely determine black hole masses, as there could be systematic differences in the black hole scaling relations between these types of galaxies and the BCGs/giant ellipticals, thereby implying different pathways for black hole and galaxy growth.

We determine the mass of the nuclear black hole ($M$) in NGC 3706, an early type galaxy with a central surface brightness minimum arising from an apparent stellar ring, which is misaligned with respect to the galaxy's major axis at larger radii. We fit new HST/STIS and archival data with axisymmetric orbit models to determine $M$, mass-to-light ratio ($\Upsilon_V$), and dark matter halo profile. The best-fit model parameters with 1$\sigma$ uncertainties are $M = (6.0^{+0.7}_{-0.9}) \times 10^8\ M_{\scriptscriptstyle \odot}$ and $\Upsilon_V = 6.0 \pm 0.2\ M_{\scriptscriptstyle \odot}\ L_{{\scriptscriptstyle \odot},V}^{-1}$ at an assumed distance of 46 Mpc. The models are inconsistent with no black hole at a significance of $\Delta\chi^2 = 15.4$ and require a dark matter halo to adequately fit the kinematic data, but the fits are consistent with a large range of plausible dark matter halo parameters. The ring is inconsistent with a population of co-rotating stars on circular orbits, which would produce a narrow line-of-sight velocity distribution (LOSVD). Instead, the ring's LOSVD has a small value of $|V|/\sigma$, the ratio of mean velocity to velocity dispersion. Based on the observed low $|V|/\sigma$, our orbit modeling, and a kinematic decomposition of the ring from the bulge, we conclude that the stellar ring contains stars that orbit in both directions. We consider potential origins for this unique feature, including multiple tidal disruptions of stellar clusters, a change in the gravitational potential from triaxial to axisymmetric, resonant capture and inclining of orbits by a binary black hole, and multiple mergers leading to gas being funneled to the center of the galaxy.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

We present stellar kinematics and orbit superposition models for the central regions of four Brightest Cluster Galaxies (BCGs), based upon integral-field spectroscopy at Gemini, Keck, and McDonald Observatories. Our integral-field data span radii from < 100 pc to tens of kpc. We report black hole masses, M_BH, of 2.1 +/- 1.6 x 10^10 M_Sun for NGC 4889, 9.7 + 3.0 - 2.6 x 10^9 M_Sun for NGC 3842, and 1.3 + 0.5 - 0.4 x 10^9 M_Sun for NGC 7768. For NGC 2832 we report an upper limit of M_BH < 9 x 10^9 M_Sun. Stellar orbits near the center of each galaxy are tangentially biased, on comparable spatial scales to the galaxies' photometric cores. We find possible photometric and kinematic evidence for an eccentric torus of stars in NGC 4889, with a radius of nearly 1 kpc. We compare our measurements of M_BH to the predicted black hole masses from various fits to the relations between M_BH and stellar velocity dispersion, luminosity, or stellar mass. The black holes in NGC 4889 and NGC 3842 are significantly more massive than all dispersion-based predictions and most luminosity-based predictions. The black hole in NGC 7768 is consistent with a broader range of predictions.

We present Chandra observations of 12 galaxies that contain supermassive black holes with dynamical mass measurements. Each galaxy was observed for 30 ksec and resulted in a total of 68 point source detections in the target galaxies including supermassive black hole sources, ultraluminous X-ray sources, and extragalactic X-ray binaries. Based on our fits of the X-ray spectra, we report fluxes, luminosities, Eddington ratios, and slope of the power-law spectrum. Normalized to the Eddington luminosity, the 2--10 keV band X-ray luminosities of the SMBH sources range from $10^{-8}$ to $10^{-6}$, and the power-law slopes are centered at $\sim2$ with a slight trend towards steeper (softer) slopes at smaller Eddington fractions, implying a change in the physical processes responsible for their emission at low accretion rates. We find 20 ULX candidates, of which six are likely ($>90%$ chance) to be true ULXs. The most promising ULX candidate has an isotropic luminosity in the 0.3--10 keV band of $1.0_{-0.3}^{+0.6} \times 10^{40}$ erg/s.

Observational work conducted over the last few decades indicates that all massive galaxies have supermassive black holes at their centres. Although the luminosities and brightness fluctuations of quasars in the early Universe suggest that some are powered by black holes with masses greater than 10 billion solar masses, the remnants of these objects have not been found in the nearby Universe. The giant elliptical galaxy Messier 87 hosts the hitherto most massive known black hole, which has a mass of 6.3 billion solar masses. Here we report that NGC 3842, the brightest galaxy in a cluster at a distance from Earth of 98 megaparsecs, has a central black hole with a mass of 9.7 billion solar masses, and that a black hole of comparable or greater mass is present in NGC 4889, the brightest galaxy in the Coma cluster (at a distance of 103 megaparsecs). These two black holes are significantly more massive than predicted by linearly extrapolating the widely-used correlations between black hole mass and the stellar velocity dispersion or bulge luminosity of the host galaxy. Although these correlations remain useful for predicting black hole masses in less massive elliptical galaxies, our measurements suggest that different evolutionary processes influence the growth of the largest galaxies and their black holes.

We present HST STIS observations of the galaxy NGC 4382 (M85) and axisymmetric models of the galaxy to determine mass-to-light ration (M/L, V-band) and central black hole mass (M_BH). We find M/L = 3.74 +/- 0.1 (solar units) and M_BH = 1.3 (+5.2, -1.2) \times 10^7 M_sun at an assumed distance of 17.9 Mpc, consistent with no black hole. The upper limit, M_BH < 9.6 \times 10^7 M_sun (2{\sigma}) or M_BH < 1.4 \times 10^8 M_sun (3{\sigma}) is consistent with the current M-{\sigma} relation, which predicts M_BH = 8.8 \times 10^7 M_sun at {\sigma}_e = 182 km/s, but low for the current M-L relation, which predicts M_BH = 7.8 \times 10^8 M_sun at L_V = 8.9 \times 10^10 L_sun,V. HST images show the nucleus to be double, suggesting the presence of a nuclear eccentric stellar disk, in analogy to the Tremaine disk in M31. This conclusion is supported by the HST velocity dispersion profile. Despite the presence of this non-axisymmetric feature and evidence of a recent merger, we conclude that the reliability of our black hole mass determination is not hindered. The inferred low black hole mass may explain the lack of nuclear activity.

We present axisymmetric, orbit-based models to study the central black hole, stellar mass-to-light ratio, and dark matter halo of NGC 4594 (M104, the Sombrero Galaxy). For stellar kinematics, we use published high-resolution kinematics of the central region taken with the Hubble Space Telescope, newly obtained Gemini long-slit spectra of the major axis, and integral field kinematics from the SAURON instrument. At large radii, we use globular cluster kinematics to trace the mass profile and apply extra leverage to recovering the dark matter halo parameters. We find a black hole of mass M_{\bullet}=(6.6 +/- 0.4) x 10^8 M_{\odot}, and determine the stellar M/L_I=3.4 +/- 0.05 (uncertainties are the 68% confidence band marginalized over the other parameters). Our best fit dark matter halo is a cored logarithmic model with asymptotic circular speed V_c=376 +/- 12 km/s and core radius r_c= 4.7 +/- 0.6 kpc. The fraction of dark to total mass contained within the half-light radius is 0.52. Taking the bulge and disk components into account in our calculation of \sigma_e puts NGC 4594 squarely on the M-\sigma relation. We also determine that NGC 4594 lies directly on the M-L relation.

We examine the possibility that the observed relation between black-hole mass and host-galaxy stellar velocity dispersion (the M-sigma relation) is biased by an observational selection effect, the difficulty of detecting a black hole whose sphere of influence is smaller than the telescope resolution. In particular, we critically investigate recent claims that the M-sigma relation only represents the upper limit to a broad distribution of black-hole masses in galaxies of a given velocity dispersion. We find that this hypothesis can be rejected at a high confidence level, at least for the early-type galaxies with relatively high velocity dispersions (median 268 km/s) that comprise most of our sample. We also describe a general procedure for incorporating observational selection effects in estimates of the properties of the M-sigma relation. Applying this procedure we find results that are consistent with earlier estimates that did not account for selection effects, although with larger error bars. In particular, (i) the width of the M-sigma relation is not significantly increased; (ii) the slope and normalization of the M-sigma relation are not significantly changed; (iii) most or all luminous early-type galaxies contain central black holes at zero redshift. Our results may not apply to late-type or small galaxies, which are not well-represented in our sample.

We present the stellar kinematics in the central 2" of the luminous elliptical galaxy M87 (NGC 4486), using laser adaptive optics to feed the Gemini telescope integral-field spectrograph, NIFS. The velocity dispersion rises to 480 km/s at 0.2". We combine these data with extensive stellar kinematics out to large radii to derive a black-hole mass equal to (6.6+-0.4)x10^9 Msun, using orbit-based axisymmetric models and including only the NIFS data in the central region. Including previously-reported ground-based data in the central region drops the uncertainty to 0.25x10^9 Msun with no change in the best-fit mass; however, we rely on the values derived from the NIFS-only data in the central region in order to limit systematic differences. The best-fit model shows a significant increase in the tangential velocity anisotropy of stars orbiting in the central region with decreasing radius; similar to that seen in the centers of other core galaxies. The black-hole mass is insensitive to the inclusion of a dark halo in the models --- the high angular-resolution provided by the adaptive optics breaks the degeneracy between black-hole mass and stellar mass-to-light ratio. The present black-hole mass is in excellent agreement with the Gebhardt & Thomas value, implying that the dark halo must be included when the kinematic influence of the black hole is poorly resolved. This degeneracy implies that the black-hole masses of luminous core galaxies, where this effect is important, may need to be re-evaluated. The present value exceeds the prediction of the black hole-dispersion and black hole-luminosity relations, both of which predict about 1x10^9 Msun for M87, by close to twice the intrinsic scatter in the relations. The high-end of the black hole correlations may be poorly determined at present.

We present the first direct measurement of the central black hole mass, M_BH, in NGC 6086, the Brightest Cluster Galaxy (BCG) in Abell 2162. Our investigation demonstrates for the first time that stellar dynamical measurements of M_BH in BCGs are possible beyond the nearest few galaxy clusters. We observed NGC 6086 with laser guide star adaptive optics and the integral-field spectrograph (IFS) OSIRIS at the W.M. Keck Observatory, and with the seeing-limited IFS GMOS-N at Gemini Observatory North. We combined the two IFS data sets with existing major-axis kinematics, and used axisymmetric stellar orbit models to determine M_BH and the R-band stellar mass-to-light ratio, M*/L_R. We find M_BH = 3.6(+1.7)(-1.1) x 10^9 M_Sun and M*/L_R = 4.6(+0.3)(-0.7) M_Sun/L_Sun (68% confidence), from models using the most massive dark matter halo allowed within the gravitational potential of the host cluster. Models fitting only IFS data confirm M_BH ~ 3 x 10^9 M_Sun and M*/L_R ~ 4 M_Sun/L_Sun, with weak dependence on the dark matter halo structure. When data out to 19 kpc are included, the unrealistic omission of dark matter causes the best-fit black hole mass to decrease dramatically, to 0.6 x 10^9 M_Sun, and the best-fit stellar mass-to-light ratio to increase to 6.7 M_Sun/L_Sun. The latter value is at further odds with stellar population studies favoring M*/L ~ 2 M_Sun/L_Sun,R. Biases from dark matter omission could extend to dynamical models of other galaxies with central stellar cores, and new measurements of M_BH from models with dark matter could steepen the empirical scaling relationships between black holes and their host galaxies.

Black hole accretion and jet production are areas of intensive study in astrophysics. Recent work has found a relation between radio luminosity, X-ray luminosity, and black hole mass. With the assumption that radio and X-ray luminosity are suitable proxies for jet power and accretion power, respectively, a broad fundamental connection between accretion and jet production is implied. In an effort to refine these links and enhance their power, we have explored the above relations exclusively among black holes with direct, dynamical mass-measurements. This approach not only eliminates systematic errors incurred through the use of secondary mass measurements, but also effectively restricts the range of distances considered to a volume-limited sample. Further, we have exclusively used archival data from the Chandra X-ray Observatory to best isolate nuclear sources. We find log(L_R) = (4.80 +/- 0.24) + (0.78 +/- 0.27) log(M_BH) + (0.67 +/- 0.12) log(L_X), in broad agreement with prior efforts. Owing to the nature of our sample, the plane can be turned into an effective mass predictor. When the full sample is considered, masses are predicted less accurately than with the well-known M-sigma relation. If obscured AGN are excluded, the plane is potentially a better predictor than other scaling measures.

We report five new measurements of central black hole masses based on STIS and WFPC2 observations with the Hubble Space Telescope and on axisymmetric, three-integral, Schwarzschild orbit-library kinematic models. We selected a sample of galaxies within a narrow range in velocity dispersion that cover a range of galaxy parameters (including Hubble type and core/power-law surface density profile) where we expected to be able to resolve the galaxy's sphere of influence based on the predicted value of the black hole mass from the M-sigma relation. We find masses in units of 10^8 solar masses for the following galaxies: NGC 3585, M_BH = 3.4 (+1.5, -0.6); NGC 3607, M_BH = 1.2 (+0.4, -0.4); NGC 4026, M_BH = 2.1 (+0.7, -0.4); and NGC 5576, M_BH = 1.8 (+0.3, -0.4), all significantly excluding M_BH = 0. For NGC 3945, M_BH = 0.09 (+0.17, -0.21), which is significantly below predictions from M-sigma and M-L relations and consistent with M_BH = 0, though the presence of a double bar in this galaxy may present problems for our axisymmetric code.

We derive improved versions of the relations between supermassive black hole mass (M_BH) and host-galaxy bulge velocity dispersion (sigma) and luminosity (L) (the M-sigma and M-L relations), based on 49 M_BH measurements and 19 upper limits. Particular attention is paid to recovery of the intrinsic scatter (epsilon_0) in both relations. We find log(M_BH / M_sun) = alpha + beta * log(sigma / 200 km/s) with (alpha, beta, epsilon_0) = (8.12 +/- 0.08, 4.24 +/- 0.41, 0.44 +/- 0.06) for all galaxies and (alpha, beta, epsilon_0) = (8.23 +/- 0.08, 3.96 +/- 0.42, 0.31 +/- 0.06) for ellipticals. The results for ellipticals are consistent with previous studies, but the intrinsic scatter recovered for spirals is significantly larger. The scatter inferred reinforces the need for its consideration when calculating local black hole mass function based on the M-sigma relation, and further implies that there may be substantial selection bias in studies of the evolution of the M-sigma relation. We estimate the M-L relationship as log(M_BH / M_sun) = alpha + beta * log(L_V / 10^11 L_sun,V) of (alpha, beta, epsilon_0) = (8.95 +/- 0.11, 1.11 +/- 0.18, 0.38 +/- 0.09); using only early-type galaxies. These results appear to be insensitive to a wide range of assumptions about the measurement errors and the distribution of intrinsic scatter. We show that culling the sample according to the resolution of the black hole's sphere of influence biases the relations to larger mean masses, larger slopes, and incorrect intrinsic residuals.

Electromagnetic observations over the last 15 years have yielded a growing appreciation for the importance of supermassive black holes (SMBH) to the evolution of galaxies, and for the intricacies of dynamical interactions in our own Galactic center. Here we show that future low-frequency gravitational wave observations, alone or in combination with electromagnetic data, will open up unique windows to these processes. In particular, gravitational wave detections in the 10^{-5}-10^{-1} Hz range will yield SMBH masses and spins to unprecedented precision and will provide clues to the properties of the otherwise undetectable stellar remnants expected to populate the centers of galaxies. Such observations are therefore keys to understanding the interplay between SMBHs and their environments.

This White Paper to the National Academy of Sciences Astro2010 Decadal Review Committee outlines some of the outstanding questions regarding the assembly history of Massive Black Holes in the nuclei of galaxies and the revolutionary contributions anticipated in this field from low-frequency gravitational wave astronomy.

We determine the mass of the black hole at the center of the spiral galaxy NGC 4258 by constructing axisymmetric dynamical models of the galaxy. These models are constrained by high spatial resolution imaging and long-slit spectroscopy of the nuclear region obtained with the {\em Hubble Space Telescope}, complemented by ground-based observations extending to larger radii. Our best mass estimate is $\MBH = (3.3 \pm 0.2) \times 10^7 \MSun $ for a distance of 7.28 Mpc (statistical errors only). This is within 15% of $ (3.82\pm 0.01) \times 10^7 \MSun$, the mass determined from the kinematics of water masers (rescaled to the same distance) assuming they are in Keplerian rotation in a warped disk. The construction of accurate dynamical models of NGC 4258 is somewhat compromised by an unresolved active nucleus and color gradients, the latter caused by variations in the stellar population and/or obscuring dust. These problems are not present in the $\sim 30$ other black hole mass determinations from stellar dynamics that have been published by us and other groups; thus, the relatively close agreement between the stellar dynamical mass and the maser mass in NGC 4258 enhances our confidence in the black hole masses determined in other galaxies from stellar dynamics using similar methods and data of comparable quality.

The largest galaxies, and in particular central galaxies in clusters, offer unique insight into understanding the mechanism for the growth of nuclear black holes. We present Hubble Space Telescope kinematics for NGC1399, the central galaxy in Fornax. We find the best-fit model contains a black hole of 5.1 +-0.7 x 10^8 Msun (at a distance of 21.1 Mpc), a factor of over 2 below the correlation of black hole mass and velocity dispersion. We also find a dramatic signature for central tangential anisotropy. The velocity profiles on adjacent sides 0.5" away from the nucleus show strong bimodality, and the central spectrum shows a large drop in the dispersion. Both of these observations point to an orbital distribution that is tangentially biased. The best-fit orbital model suggests a ratio of the tangential to radial internal velocity dispersions of three. This ratio is the largest seen in any galaxy to date and will provide an important measure for the mode by which the central black hole has grown.

Programs to observe evolution in the Mbh-sigma or Mbh-L relations typically compare black-hole masses, Mbh, in high-redshift galaxies selected by nuclear activity to Mbh in local galaxies selected by luminosity L, or stellar velocity dispersion sigma. Because AGN luminosity is likely to depend on Mbh, selection effects are different for high-redshift and local samples, potentially producing a false signal of evolution. This bias arises because cosmic scatter in the Mbh-sigma and Mbh-L relations means that the mean log(L) or log(sigma) among galaxies that host a black hole of given Mbh, may be substantially different than the log(L) or log(sigma) obtained from inverting the Mbh-L or Mbh-sigma relations for the same nominal Mbh. The bias is particularly strong at high Mbh, where the luminosity and dispersion functions of galaxies are falling rapidly. The most massive black holes occur more often as rare outliers in galaxies of modest mass than in the even rarer high-mass galaxies, which would otherwise be the sole location of such black holes in the absence of cosmic scatter. Because of this bias, Mbh will typically appear to be too large in the distant sample for a given L or sigma. For the largest black holes and the largest plausible cosmic scatter, the bias can reach a factor of 3 in Mbh for the Mbh-sigma relation and a factor of 9 for the Mbh-L relation. Unfortunately, the actual cosmic scatter is not known well enough to correct for the bias. Measuring evolution of the Mbh and galaxy property relations requires object selection to be precisely defined and exactly the same at all redshifts.

A variety of host galaxy (bulge) parameters are examined in order to determine their predictive power in ascertaining the masses of the supermassive black holes (SMBH) at the centers of the galaxies. Based on a sample of 23 nearby galaxies, comprised of both elliptical galaxies and spiral/lenticular bulges, we identify a strong correlation between the bulge gravitational binding energy ($E_g$), as traced by the stellar light profile, and the SMBH mass ($M_{\bullet}$), such that $M_{\bullet} \propto E_g^{0.6}$. The scatter about the relationship indicates that this is as strong a predictor of $M_{\bullet}$ as the velocity dispersion ($\sigma$), for the elliptical galaxy subsample. Improved mass-to-light ratios, obtained with IFU spectroscopy and I-band photometry by the SAURON group, were used for those sample galaxies where available, resulting in an energy predictor with the same slope, but with reduced scatter. Alternative $M_{\bullet}$ predictors such as the gravitational potential and the bulge mass are also explored, but these are found to be inferior when compared with both the bulge gravitational binding energy and bulge velocity dispersion predictors, for the full galaxy sample.

We describe results from a new ground-based monitoring campaign on NGC 5548, the best studied reverberation-mapped AGN. We find that it was in the lowest luminosity state yet recorded during a monitoring program, namely L(5100) = 4.7 x 10^42 ergs s^-1. We determine a rest-frame time lag between flux variations in the continuum and the Hbeta line of 6.3 (+2.6/-2.3) days. Combining our measurements with those of previous campaigns, we determine a weighted black hole mass of M_BH = 6.54 (+0.26/-0.25) x 10^7 M_sun based on all broad emission lines with suitable variability data. We confirm the previously-discovered virial relationship between the time lag of emission lines relative to the continuum and the width of the emission lines in NGC 5548, which is the expected signature of a gravity-dominated broad-line region. Using this lowest luminosity state, we extend the range of the relationship between the luminosity and the time lag in NGC 5548 and measure a slope that is consistent with alpha = 0.5, the naive expectation for the broad line region for an assumed form of r ~ L^alpha. This value is also consistent with the slope recently determined by Bentz et al. for the population of reverberation-mapped AGNs as a whole.

Black hole masses predicted from the Mbh-sigma relationship conflict with those predicted from the Mbh-L relationship for the most luminous galaxies, such as brightest cluster galaxies (BCGs). This is because stellar velocity dispersion, sigma, increases only weakly with L for BCGs and other giant ellipticals. The Mbh-L relationship predicts that the most luminous BCGs may have Mbh approaching 10^{10}M_sol, while the M-sigma relationship always predicts Mbh<3X10^9M_sol. We argue that the Mbh-L relationship is a plausible or even preferred description for BCGs and other galaxies of similar luminosity. If cores in central stellar density are formed by binary BHs, then the inner-core cusp radius, r_gamma, may be an independent witness of Mbh. Using structural parameters derived from a large sample of early-type galaxies observed by HST, we argue that L is superior to sigma as an indicator of r_gamma in luminous galaxies. The observed r_gamma-Mbh relationship for 11 core galaxies with measured Mbh appears to be consistent with the Mbh-L relationship for BCGs. BCGs have large cores appropriate for their large luminosities that may be difficult to generate with the modest BH masses inferred from the Mbh-sigma relationship. Mbh~L may be expected to hold for BCGs, if they were formed in dissipationless mergers, which should preserve ratio of BH to stellar mass. This picture appears to be consistent with the slow increase in sigma with L and the more rapid increase in effective radii with L seen in BCGs. If BCGs have large BHs commensurate with their luminosities, then the local BH mass function for Mbh>3X10^9M_sol may be nearly an order of magnitude richer than that inferred from the Mbh-sigma relationship. The volume density of QSOs at earlier epochs may favor the predictions from the Mbh-L relationship.

We combine the results from several HST investigations of the central structure of early-type galaxies to generate a large sample of parameterized surface photometry. The studies included were those that used the "Nuker law" to characterize the inner light distributions of the galaxies. The sample comprises WFPC1 and WFPC2 V band observations published earlier by our group, R band WFPC2 photometry of Rest et al., NICMOS H band photometry by Ravindranath et al. and Quillen et al., and the BCG WFPC2 I band photometry of Laine et al. The distribution of the logarithmic slopes of the central profiles strongly affirms that the central structure of elliptical galaxies with Mv < -19 is bimodal, based on both parametric and non-parametric analysis. At the HST resolution limit, most galaxies are either power-law systems, which have steep cusps in surface brightness, or core systems, which have shallow cusps interior to a steeper envelope brightness distribution. A rapid transition between the two forms occurs over the luminosity range -22 < Mv < -20, with cores dominating at the highest luminosities, and power-laws at the lowest. There are a few "intermediate" systems that have both cusp slopes and total luminosities that fall within the core/power-law transition, but they are rare and do not fill in the overall bimodal distribution of cusp slopes. These results are inconsistent with the Ferrarese et al. Virgo Cluster Survey (VCS) analysis. However, using galaxies common to the VCS samples, we demonstrate that the VCS models of the cusps are either a poor match to the observations or consist of forms fitted to the galaxy envelopes and extrapolated inward to the HST resolution limit.

We present new observations leading to an improved black hole mass estimate for the Seyfert 1 galaxy NGC 4593 as part of a reverberation-mapping campaign conducted at the MDM Observatory. Cross-correlation analysis of the H_beta emission-line light curve with the optical continuum light curve reveals an emission-line time delay of 3.73 (+-0.75) days. By combining this time delay with the H_beta line width, we derive a central black hole mass of M_BH = 9.8(+-2.1)x10^6 M_sun, an improvement in precision of a factor of several over past results.

We have undertaken a new ground-based monitoring campaign to improve the estimates of the mass of the central black hole in NGC 4151. We measure the lag time of the broad H beta line response compared to the optical continuum at 5100 A and find a lag of 6.6 (+1.1/-0.8) days. We combine our data with the recent reanalysis of UV emission lines by Metzroth et al. to calculate a weighted mean of the black hole mass, M_BH = 4.57 (+0.57/-0.47) x 10^7 M_sun. The absolute calibration of the black hole mass is based on normalization of the AGN black hole mass - stellar velocity dispersion (M_BH - sigma_*) relationship to that of quiescent galaxies by Onken et al. The scatter in the M_BH - sigma_* relationship suggests that reverberation-mapping based mass measurements are typically uncertain by a factor of 3-4.

We present HST spectroscopy of the nucleus of M31 obtained with STIS. Spectra taken around the CaT lines at 8500 see only the red giants in the double bright- ness peaks P1 and P2. In contrast, spectra taken at 3600-5100 A are sensitive to the tiny blue nucleus embedded in P2, the lower surface brightness red nucleus. P2 has a K-type spectrum, but the embedded blue nucleus has an A-type spectrum with strong Balmer absorption lines. Given the small likelihood for stellar collisions, a 200 Myr old starburst appears to be the most plausible origin of the blue nucleus. In stellar population, size, and velocity dispersion, the blue nucleus is so different from P1 and P2 that we call it P3. The line-of-sight velocity distributions of the red stars in P1+P2 strengthen the support for Tremaine s eccentric disk model. The kinematics of P3 is consistent with a circular stellar disk in Keplerian rotation around a super-massive black hole with M_bh = 1.4 x 10^8 M_sun. The P3 and the P1+P2 disks rotate in the same sense and are almost coplanar. The observed velocity dispersion of P3 is due to blurred rotation and has a maximum value of sigma = 1183+-201 km/s. The observed peak rotation velocity of P3 is V = 618+-81 km/s at radius 0.05" = 0.19 pc corresponding to a circular rotation velocity at this radius of ~1700 km/s. Any dark star cluster alternative to a black hole must have a half-mass radius <= 0.03" = 0.11 pc. We show that this excludes clusters of brown dwarfs or dead stars on astrophysical grounds.

The infrared spectrograph GNIRS on Gemini South unlocks new possibilities to study the central black holes in dusty galaxies that have been inaccessible to previous black hole studies. We exploit good near-infrared seeing to measure the central black hole (BH) of Centaurus A (NGC 5128). We measure the stellar kinematics of NGC 5128 using the region around the CO bandheads at 2.3 microns and determine the black hole mass using orbit-based models. Black holes are believed to be essential components of galaxies, and their evolutionary states appear to be closely linked to those of their hosts. Our current knowledge does not go much beyond this; galaxies such as NGC 5128 (an AGN and recent merger) can further develop this knowledge. However, NGC 5128 and galaxies like it contain large amounts of dust which hamper optical spectroscopy, making near-infrared measurements an attractive alternative. We find a BH mass of 2.4e8 solar masses for an edge-on model, 1.8e8 solar masses for a model with inclination of 45 degrees, and 1.5e8 solar masses for a model with inclination of 20 degrees. We adopt the value for the edge-on model, which has significantly lower chi^2. These estimates are consistent with a previous gas dynamical estimate and are five to ten times higher than that predicted by the correlation between BH mass and velocity dispersion. If NGC 5128 will eventually follow the trend for quiescent galaxies, this result suggests that its BH assembled first before its host component. NGC 5128 thus provides an important example for our knowledge of central black holes; this technique can be applied to other such galaxies to further explore this question.

We analyze HST+WFPC2 images of 77 early-type galaxies. Brightness profiles are classed into "core" or "power-law" forms. Cores are typically rounder than power-law galaxies. Nearly all power-laws with central ellipticity >=0.3 have stellar disks, implying that disks are present in power-laws with epsilon <0.3, but are not visible due to unfavorable geometry. A few low-luminosity core galaxies also have disks; these may be transition forms from power-laws. Cores and power-laws both have twisting isophotes at small radii. Core galaxies have somewhat weaker color gradients than power-laws. Nuclei are found in 29% of the cores and 60% of the power-laws. Nuclei are typically bluer than the surrounding galaxy. NGC 4073 and 4382 have central minima in their intrinsic starlight distributions; NGC 4382 resembles the double nucleus of M31. In general, the peak brightness location is coincident with the photocenter of the core to <1 pc. Five galaxies, however, have centers significantly displaced from their cores; these may be unresolved asymmetric double nuclei. Central dust is visible in half of the galaxies. The presence and strength of dust correlates with nuclear emission. The prevalence of dust and its morphology suggests that dust clouds form, settle to the center, and disappear repeatedly on ~10^8 yr timescales. We discuss the hypothesis that cores are created by the decay of a massive black hole binary. Apart from their brightness profiles, there are no strong differences between cores and power-laws that demand this scenario; however, the rounder shapes of cores, their lack of disks, and their reduced color gradients may be consistent with it.

This is the first of a series of papers dedicated to unveil the mass composition and dynamical structure of a sample of flattened early type galaxies in the Coma cluster. We describe our modifications to the Schwarzschild code of Richstone et al. (in preparation). Applying a Voronoi tessellation in the surface of section we are able to assign accurate phase-space volumes to individual orbits and to reconstruct the full three-integral phase-space distribution function (DF) of any axisymmetric orbit library. Two types of tests have been performed to check the accuracy with which DFs can be represented by appropriate orbit libraries. First, by mapping DFs of spherical gamma-models and flattened Plummer models on to the library we show that the resulting line-of-sight velocity distributions and internal velocity moments of the library match those directly derived from the DF to a precision better than that of present day observational errors. Second, by fitting libraries to the projected kinematics of the same DFs we show that the distribution function reconstructed from the fitted library matches the input DF to a rms of about 15 per cent over a region in phase-space covering 90 per cent of the mass of the library. The achieved accuracy allows us to implement effective entropy-based regularisation to fit real, noisy and spatially incomplete data.

The quasar HE 0047-1756, at z=1.67, is found to be split into two images 1.44" apart by an intervening galaxy acting as a gravitational lens. The flux ratio for the two components is roughly 3.5:1, depending slightly upon wavelength. The lensing galaxy is seen on images obtained at 800 nm and 2.1 \mu; there is also a nearby faint object which may be responsible for some shear. The spectra of the two quasar images are nearly identical, but the emission line ratio between the two components scale differently from the continuum. Moreover, the fainter component has a bluer continuum slope than the brighter one. We argue that these small differences are probably due to microlensing. There are hints of an Einstein ring emanating from the brighter image toward the fainter one.

We show that orbit-superposition dynamical models (Schwarzschild's method) provide reliable estimates of nuclear black hole masses and errors when constructed from adequate orbit libraries and kinematic data. We thus rebut two recent papers that argue that BH masses obtained from this method are unreliable. These papers claim to demonstrate that the range of allowable BH masses derived from a given dataset is artificially too narrow as a result of an inadequate number of orbits in the library used to construct dynamical models. This is an elementary error that is easily avoided. We describe a method to estimate the number and nature of orbits needed for the library. We provide an example that shows that this prescription is adequate, in the sense that the range of allowable BH masses is not artificially narrowed by use of too few orbits. A second point raised by critics is that kinematic data are generally obtained with insufficient spatial resolution to obtain a reliable mass. We make the distinction between unreliable determinations and imprecise ones. We show that there are several different properties of a kinematic dataset that can lead to imprecise BH determinations, but none of the attributes we have investigated leads to an unreliable determination. In short, the degree to which the BH radius of influence is resolved by spectroscopic observations is already reflected in the BH-mass error envelope, and is not a hidden source of error. The BH masses published by our group and the Leiden group are reliable.

We present stellar kinematics for a sample of 10 early-type galaxies observed using the STIS aboard the Hubble Space Telescope, and the Modular Spectrograph on the MDM Observatory 2.4-m telescope. The spectra are used to derive line-of-sight velocity distributions (LOSVDs) of the stars using a Maximum Penalized Likelihood method. We use Gauss-Hermite polynomials to parameterize the LOSVDs and find predominantly negative h4 values (boxy distributions) in the central regions of our galaxies. One galaxy, NGC 4697, has significantly positive central h4 (high tail weight). The majority of galaxies have a central velocity dispersion excess in the STIS kinematics over ground-based velocity dispersions. The galaxies with the strongest rotational support, as quantified with v_MAX/sigma_STIS, have the smallest dispersion excess at STIS resolution. The best-fitting, general, axisymmetric dynamical models (described in a companion paper) require black holes in all cases, with masses ranging from 10^6.5 to 10^9.3 Msun. We replot these updated masses on the BH/sigma relation, and show that the fit to only these 10 galaxies has a slope consistent with the fits to larger samples. The greatest outlier is NGC 2778, a dwarf elliptical with relatively poorly constrained black hole mass. The two best candidates for pseudobulges, NGC 3384 and 7457, do not deviate significantly from the established relation between black hole and sigma. Neither do the three galaxies which show the most evidence of a recent merger, NGC 3608, 4473, and 4697.

Motivated by the observed correlation between black hole masses $\mbh$ and the velocity dispersion $\sigma$ of host galaxies, we develop a theoretical model of black hole formation in galactic bulges (this paper generalizes an earlier ApJ Letter). The model assumes an initial state specified by a a uniform rotation rate $\Omega$ and a density distribution of the form $\rho = \aeff^2 / 2 \pi G r^2$ (so that $\aeff$ is an effective transport speed). The black hole mass is determined when the centrifugal radius of the collapse flow exceeds the capture radius of the central black hole (for Schwarzschild geometry). This model reproduces the observed correlation between the estimated black hole masses and the velocity dispersions of galactic bulges, i.e., $\mbh \approx 10^8 M_\odot (\sigma/200 {\rm km s^{-1}})^4$, where $\sigma = \sqrt{2} \aeff$. To obtain this normalization, the rotation rate $\Omega \approx 2 \times 10^{15}$ rad/s. The model also defines a bulge mass scale $M_B$. If we identify the scale $M_B$ with the bulge mass, the model determines the ratio $\mrat$ of black hole mass to the host mass: $\mrat$ $\approx$ 0.0024 $(\sigma/200 {\rm km s^{-1}})$, again in reasonable agreement with observed values. In this scenario, supermassive black holes form quickly (in $\sim10^5$ yr) and are born rapidly rotating (with $a/M \sim 0.9$). This paper also shows how these results depend on the assumed initial conditions; the most important quantity is the initial distribution of specific angular momentum in the pre-collapse state.

Supermassive black holes are thought to be relics of quasars, and their numbers and masses are therefore related to the quasar luminosity function and its evolution with redshift. We have used the relationship between black hole mass and bulge velocity dispersion (the M_bullet - sigma relation) to make an improved estimate of the mass density and mass spectrum of supermassive black holes. Uncertainties in the M_bullet - sigma relation have little effect on the mass density. We find a mass density of (4.8 +/- 1.6) h^2 x 10^5 M_sun Mpc^-3. Some of the variance in published density estimates comes from the use of different values of the Hubble constant.

We present axisymmetric, orbit superposition models for 12 galaxies using data taken with the Hubble Space Telescope (HST) and ground-based observatories. In each galaxy, we detect a central black hole (BH) and measure its mass to accuracies ranging from 10% to 70%. We demonstrate that in most cases the BH detection requires both the HST and ground-based data. Using the ground-based data alone does provide an unbiased measure of the BH mass (provided they are fit with fully general models), but at a greatly reduced significance. The most significant correlation with host galaxy properties is the relation between the BH mass and the velocity dispersion of the host galaxy; we find no other equally strong correlation, and no second parameter that improves the quality of the mass-dispersion relation. We are also able to measure the stellar orbital properties from these general models. The most massive galaxies are strongly biased to tangential orbits near the BH, consistent with binary BH models, while lower-mass galaxies have a range of anisotropies, consistent with an adiabatic growth of the BH.

We used HST WFPC2 images to identify six early-type galaxies with surface- brightness profiles that decrease inward over a limited range of radii near their centers. The implied luminosity density profiles of these galaxies have local minima interior to their core break radii. NGC 3706 harbors a high surface brightness ring of starlight with radius ~20 pc. Its central structure may be related to that in the double-nucleus galaxies M31 and NGC 4486B. NGC 4406 and NGC 6876 have nearly flat cores that on close inspection are centrally depressed. Colors for both galaxies imply that this is not due to dust absorption. The surface brightness distributions of both galaxies are consistent with stellar tori that are more diffuse than the sharply defined system in NGC 3706. The remaining three galaxies are the brightest cluster galaxies in A260, A347, and A3574. Color information is not available for these objects, but they strongly resemble NGC 4406 and NGC 6876 in their cores. The thin ring in NGC 3706 may have formed dissipatively. The five other galaxies resemble the endpoints of some simulations of the merging of two gas-free stellar systems, each harboring a massive nuclear black hole. In one version of this scenario, diffuse stellar tori are produced when stars initially bound to one black hole are tidally stripped away by the second black hole. Alternatively, some inward-decreasing surface-brightness profiles may reflect the ejection of stars from a core during the hardening of the binary black hole created during the merger.

Observations of nearby galaxies reveal a strong correlation between the mass of the central dark object M and the velocity dispersion sigma of the host galaxy, of the form log(M/M_sun) = a + b*log(sigma/sigma_0); however, published estimates of the slope b span a wide range (3.75 to 5.3). Merritt & Ferrarese have argued that low slopes (<4) 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 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 b=4.02 +/- 0.32, a=8.13 +/- 0.06, for sigma_0=200 km/s. The M-sigma relation has an intrinsic dispersion in log M that is no larger than 0.3 dex. In an Appendix, we present a simple model for the velocity-dispersion profile of the Galactic bulge.

Galaxies that contain bulges appear to contain central black holes whose masses correlate with the velocity dispersion of the bulge. We show that no corresponding relationship applies in the pure disk galaxy M33. Three-integral dynamical models fit Hubble Space Telescope WFPC2 photometry and STIS spectroscopy best if the central black hole mass is zero. The upper limit is 1500 M_sun. This is significantly below the mass expected from the velocity dispersion of the nucleus and far below any mass predicted from the disk kinematics. Our results suggest that supermassive black holes are associated only with galaxy bulges and not with their disks.

We analyze the nuclear stellar dynamics of the SB0 galaxy NGC 1023, utilizing observational data both from the Space Telescope Imaging Spectrograph aboard the Hubble Space Telescope and from the ground. The stellar kinematics measured from these long-slit spectra show rapid rotation (V = 70 km/s at a distance of 0.1 arcsec = 4.9 pc from the nucleus) and increasing velocity dispersion toward the nucleus (where sigma = 295 +/- 30 km/s). We model the observed stellar kinematics assuming an axisymmetric mass distribution with both two and three integrals of motion. Both modeling techniques point to the presence of a central dark compact mass (which presumably is a supermassive black hole) with confidence > 99%. The isotropic two-integral models yield a best-fitting black hole mass of (6.0 +/- 1.4) x 10^7 M_sun and mass-to-light ratio (M/L_V) of 5.38 +/- 0.08, and the goodness-of-fit (chi^2) is insensitive to reasonable values for the galaxy's inclination. The three-integral models, which non-parametrically fit the observed line-of-sight velocity distribution as a function of position in the galaxy, suggest a black hole mass of (3.9 +/- 0.4) x 10^7 M_sun and M/L_V of 5.56 +/- 0.02 (internal errors), and the edge-on models are vastly superior fits over models at other inclinations. The internal dynamics in NGC 1023 as suggested by our best-fit three-integral model shows that the velocity distribution function at the nucleus is tangentially anisotropic, suggesting the presence of a nuclear stellar disk. The nuclear line of sight velocity distribution has enhanced wings at velocities >= 600 km/s from systemic, suggesting that perhaps we have detected a group of stars very close to the central dark mass.

We construct a model for the formation of black holes within galactic bulges. The initial state is a slowly rotating isothermal sphere, characterized by effective transport speed $\aeff$ and rotation rate $\Omega$. The black hole mass is determined when the centrifugal radius of the collapse flow exceeds the capture radius of the central black hole. This model reproduces the observed correlation between black hole masses and galactic velocity dispersions, $\mbh \approx 10^8 M_\odot (\sigma/200 \kms)^4$, where $\sigma = \sqrt{2} \aeff$. This model also predicts the ratio $\mrat$ of black hole mass to host mass: $\mrat$ $\approx$ 0.004 $(\sigma/200 \kms)$.

Black hole (BH) masses that have been measured by reverberation mapping in active galaxies fall significantly below the correlation between bulge luminosity and BH mass determined from spatially resolved kinematics of nearby normal galaxies. This discrepancy has created concern that one or both techniques suffer from systematic errors. We show that BH masses from reverberation mapping are consistent with the recently discovered relationship between BH mass and galaxy velocity dispersion. Therefore the bulge luminosities are the probable source of the disagreement, not problems with either method of mass measurement. This result underscores the utility of the BH mass -- velocity dispersion relationship. Reverberation mapping can now be applied with increased confidence to galaxies whose active nuclei are too bright or whose distances are too large for BH searches based on spatially resolved kinematics.

We describe a correlation between the mass M_BH of a galaxy's central black hole and the luminosity-weighted line-of-sight velocity dispersion sigma_e within the half-light radius. The result is based on a sample of 26 galaxies, including 13 galaxies with new determinations of black hole masses from Hubble Space Telescope measurements of stellar kinematics. The best-fit correlation is M_BH = 1.2 (+-0.2) x 10^8 M_sun (sigma_e/200 km/s)^(3.75 (+-0.3))over almost three orders of magnitude in M_BH; the scatter in M_BH at fixed sigma_e is only 0.30 dex and most of this is due to observational errors. The M_BH-sigma_e relation is of interest not only for its strong predictive power but also because it implies that central black hole mass is constrained by and closely related to properties of the host galaxy's bulge.

We fit axisymmetric 3-integral dynamical models to NGC3379 using the line-of-sight velocity distribution obtained from HST/FOS spectra of the galaxy center and ground-based long-slit spectroscopy along four position angles, with the light distribution constrained by WFPC2 and ground-based images. We have fitted models with inclinations from 29 (intrinsic galaxy type E5) to 90 degrees (intrinsic E1) and black hole masses from 0 to 1e9 M_solar. The best-fit black hole masses range from 6e7 to 2e8 M_solar, depending on inclination. The velocity ellipsoid of the best model is not consistent with either isotropy or a two-integral distribution function. Along the major axis, the velocity ellipsoid becomes tangential at the innermost bin, radial in the mid-range radii, and tangential again at the outermost bins. For the acceptable models, the radial to tangential dispersion in the mid-range radii ranges from 1.1 < sigma_r / sigma_t < 1.7. Compared with these 3-integral models, 2-integral isotropic models overestimate the black hole mass since they cannot provide adequate radial motion. However, the models presented in this paper still contain restrictive assumptions-namely assumptions of constant M/L and spheroidal symmetry-requiring yet more models to study black hole properties in complete generality.

In order for the core Fundamental Plane (cFP) to endure, small ellipticals must not survive mergers with giant ellipticals, despite the fact the small secondary galaxy can be as much as a million times more dense than the primary. However, our previous set of experiments has shown that, for purely stellar galaxies, the secondary does in fact survive mergers with primaries up to 100 times more massive. In this paper, we demonstrate the effect that a massive central black hole has on mergers of cFP galaxies. Our results indicate that the addition of a massive central singularity inside the primary galaxy provides strong enough tidal forces to destroy dense cFP companions when the secondary's orbital decay is sufficiently elongated. The destruction of the secondary acts to preserve the original low central density profile of the primary in the giant merger remnant, which allows the remnant to remain on the cFP. On more circular orbits, though, the secondary is only disrupted near the end of the merger, and the degree to which the secondary particles disperse depends on the amount of orbital energy left in the merger. Hence, there are some mergers for which the cFP is not preserved in our experiments. We find that if the secondary is not dispersed, it forms a spinning stellar disk with a central density that forces the merger remnant off the cFP.

We have developed an effective simulation method and applied it to the problem of the accretion of very dense secondary companions by tenuous primaries. We have studied the accretion of objects of varying luminosity ratios, with sizes and densities drawn from the core Fundamental Plane for elliptical galaxies under the assumption that the mass distribution in the central parts of the galaxies follows the light. The results indicate that in mergers with mass ratios greater than 10, the smaller object is only stripped down to the highest density encountered in the primary during the accretion process. Thus, the form of the core Fundamental Plane suggests that the mass distribution in these galaxy centers is different than the light distribution, or that mergers on the core Fundamental Plane require more than the dynamics of visible stars.

Recent surveys suggest that most or all normal galaxies host a massive black hole with 1/100 to 1/1000 of the visible mass of the spheroid of the galaxy. Various lines of argument suggest that these galaxies have merged at least once in our past lightcone, and that the black holes have also merged. This leads to a merger rate of massive black holes of about 1/\yrs.

Black holes, an extreme consequence of the mathematics of General Relativity, have long been suspected of being the prime movers of quasars, which emit more energy than any other objects in the Universe. Recent evidence indicates that supermassive black holes, which are probably quasar remnants, reside at the centers of most galaxies. As our knowledge of the demographics of these relics of a violent earlier Universe improve, we see tantalizing clues that they participated intimately in the formation of galaxies and have strongly influenced their present-day structure.

We report the results of simulations of dense rotating stellar systems whose members suffer collisions and undergo stellar evolution processes. The initial configuration for each experiment is an isotropic Kuzmin-Kutuzov model. The dynamical evolution is simulated with the N-body tree code of Hernquist, modified to incorporate physical stellar collisions, stellar evolution, stellar mass loss and compact remnant formation, and star formation. In some simulations we have added a large accreting central black hole. In all systems the velocity dispersion in the halo evolves toward a radially biased state. In systems containing a central black hole, the dispersion becomes tangentially biased in the core, whereas it remains isotropic in systems with no black hole. Collisions tend to produce a single dominant stellar merger product, as opposed to a swarm of intermediate-mass stars. In cases where we have suppressed all processes except relaxation and physical collisions, objects with greater flattening produce larger stars through mergers. In systems where stellar ejecta are allowed to escape the system, mass loss from the heavy core stars temporarily reduces the core density and collision rate. Most of the simulations performed reproduced the ratio of the central collision time scale to the central relaxation time scale found in the dwarf elliptical galaxy M32. The rapid central evolution of these systems due to collisions and relaxation, combined with scaling the results in N, suggests that we are either viewing M32 at a peculiar moment in its history, or that its dynamically-inferred central density is at least in part due to the present of a massive dark object, presumably a black hole.

This paper, and its companion, investigate the evolution of dense stellar systems due to the influence of two-body gravitational encounters, physical collisions and stellar evolution. Our goal is the simulation of the densest centers of galaxies, like M32, which reach stellar densities near a million solar masses per cubic parsec and which may harbor black holes. These systems have a different Safronov number (the dimensionless ratio of stellar binding energy to mean stellar kinetic energy) than globular clusters, substantially increasing the importance of physical collisions relative to gravitational encounters. In this paper, we focus only on the gravitational encounters. We demonstrate, first, that our simulations with small N with a Hernquist tree code yield results basically in accord with years of effort studying globular clusters. Second, we investigate (crudely) core collapse in rotating systems with mass segregation, to separate out the effects purely due to two-body encounters from those seen in the more complex second paper. Consistent with previous studies, we find that systems whose constituent particles follow the Salpeter initial mass function rapidly undergo core collapse through Spitzer's mass segregation instability. All rotationally flattened systems show a decrease in flattening with time, consistent with Fokker-Planck calculations. An interesting new result concerns the well-established inability of simulators to identify a static center in their simulations of collisional systems. We find that the lump of high mass stars which condenses at the center wanders about the core in Brownian motion.

We have measured the mass distribution in the central 35" = 1.7 kpc of the E5 galaxy NGC 3377. Stellar rotation velocity and velocity dispersion profiles (seeing sigma_* = 0.20 - 0.56 arcsec) and V-band surface photometry (sigma_* = 0.20 - 0.26 arcsec) have been obtained with the Canada-France-Hawaii Telescope. NGC 3377 is kinematically similar to M32: the central kinematic gradients are steep. The apparent velocity dispersion rises from 95 +- 2 km/s at 1" <= r <= 4" to 178 +- 10 km/s at the center. To search for a central black hole, we derive three-dimensional velocity and velocity dispersion fields that fit our observations and Hubble Space Telescope surface photometry after projection and seeing convolution. Isotropic models imply that the mass-to-light ratio rises by a factor of ~ 4 at r < 2" to M/L_V > 10. If the mass-to-light ratio of the stars is constant with radius, then NGC 3377 contains a central massive dark object (MDO), probably a black hole, of mass (1.8 +- 0.8) x 10**8 M_Sun. Several arguments suggest that NGC 3377 is likely to be nearly isotropic. However, flattened, anisotropic maximum entropy models can fit the data without an MDO. Therefore the MDO detection in NGC 3377 is weaker than those in M31, M32, or NGC 3115. The above MDO mass is corrected for the E5 shape of the galaxy and for the difference between velocity moments and velocities given by Gaussian fits to the line profiles. We show that the latter correction does not affect the strength of the MDO detection, but it slightly reduces the mass and mass-to-light ratio. At 3" <= r <= 35", M/L_V is constant at ~ 2.4. Therefore the inner parts of NGC 3377 are dominated by a normal old stellar population, not by halo dark matter.

We analyze Hubble Space Telescope surface-brightness profiles of 61 elliptical galaxies and spiral bulges (hot galaxies). Luminous hot galaxies have cuspy cores with steep outer power-law profiles that break at r ~ r_b to shallow inner profiles with logslope less than 0.3. Faint hot galaxies show steep, largely featureless power-law profiles at all radii and lack cores. The centers of power-law galaxies are up to 1000 times denser in mass and luminosity than the cores of large galaxies at a limiting radius of 10 pc. At intermediate magnitudes (-22.0 < M_V < -20.5), core and power-law galaxies coexist, and there is a range in r_b at a given luminosity of at least two orders of magnitude. Central properties correlate with global rotation and shape: core galaxies tend to be boxy and slowly rotating, whereas power-law galaxies tend to be disky and rapidly rotating. The dense power-law centers of disky, rotating galaxies are consistent with their formation in gas-rich mergers. The parallel proposition that cores are simply the by-products of gas-free stellar mergers is less compelling. For example, core galaxies accrete small, dense, gas-free galaxies at a rate sufficient to fill in low-density cores if the satellites survived and sank to the center. An alternative model for core formation involves the orbital decay of massive black holes (BHs): the BH may heat and eject stars from the center, eroding a power law if any exists and scouring out a core. An average BH mass per spheroid of 0.002 times the stellar mass yields reasonably good agreement with the masses and radii of observed cores and in addition is consistent with the energetics of AGNs and kinematic detections of BHs in nearby galaxies.

We construct dynamical models for a sample of 36 nearby galaxies with Hubble Space Telescope photometry and ground-based kinematics. The models assume that each galaxy is axisymmetric, with a two-integral distribution function, arbitrary inclination angle, a position-independent stellar mass-to-light ratio Upsilon, and a central massive dark object (MDO) of arbitrary mass M_bh. They provide acceptable fits to 32 of the galaxies for some value of M_bh and Upsilon; the four galaxies that cannot be fit have kinematically decoupled cores. The mass-to-light ratios inferred for the 32 well-fit galaxies are consistent with the fundamental plane correlation Upsilon \propto L^0.2, where L is galaxy luminosity. In all but six galaxies the models require at the 95% confidence level an MDO of mass M_bh ~ 0.006 M_bulge = 0.006 Upsilon L. Five of the six galaxies consistent with M_bh=0 are also consistent with this correlation. The other (NGC 7332) has a much stronger upper limit on M_bh. We consider various parameterizations for the probability distribution describing the correlation of the masses of these MDOs with other galaxy properties. One of the best models can be summarized thus: a fraction f ~0.97 of galaxies have MDOs, whose masses are well described by a Gaussian distribution in log (M_bh/M_bulge) of mean -2.27 and width ~0.07.

We present optical images of the nucleus of the nearby radio galaxy M84 (NGC 4374 = 3C272.1) obtained with the Wide Field/Planetary Camera 2 (WFPC2) aboard the Hubble Space Telescope (HST). Our three images cover the H$\alpha$ + [N II] emission lines as well as the V and I continuum bands. Analysis of these images confirms that the H$\alpha$ + [N II] emission in the central 5'' (410 pc) is elongated along position angle (P.A.) $\approx 72\arcdeg$, which is roughly parallel to two nuclear dust lanes.Our high-resolution images reveal that the H$\alpha$ + [N II] emission has three components, namely a nuclear gas disk,an `ionization cone', and outer filaments. The nuclear disk of ionized gas has diameter $\approx 1'' = 82$ pc and major axis P.A. $\approx 58\arcdeg \pm 6\arcdeg$. On an angular scale of $0\farcs5$, the major axis of this nuclear gas disk is consistent with that of the dust. However, the minor axis of the gas disk (P.A. $\approx 148\arcdeg$) is tilted with respect to that of the filamentary H$\alpha$ + [N II] emission at distances > 2'' from the nucleus; the minor axis of this larger scale gas is roughly aligned with the axis of the kpc-scale radio jets (P.A. $\approx 170\arcdeg$). The ionization cone (whose apex is offset by $\approx 0\farcs3$ south of the nucleus) extends 2'' from the nucleus along the axis of the southern radio jet. This feature is similar to the ionization cones seen in some Seyfert nuclei, which are also aligned with the radio axes.

The stellar kinematics of the dwarf elliptical galaxy NGC 4486B have been measured in seeing sigma_* = .22 arcsec with the Canada-France-Hawaii Telescope. Lauer et al. 1996, ApJ, 471, L79 have shown that NGC 4486B is similar to M31 in having a double nucleus. We show that it also resembles M31 in its kinematics. The velocity dispersion gradient is very steep: sigma increases from 116 +- 6 km/s at r = 2" - 6" to 281 +- 11 km/s at the center. This is much higher than expected for an elliptical galaxy of absolute magnitude M_B = -16.8: NGC 4486B is far above the scatter in the Faber-Jackson correlation between sigma and bulge luminosity. Therefore the King core mass-to-light ratio, M/L_V = 20, is unusually high compared with normal values for old stellar populations. We construct dynamical models with isotropic velocity dispersions and show that they reproduce black hole (BH) masses derived by more detailed methods. We also fit axisymmetric, three-integral models. Isotropic models imply that NGC 4486B contains a central dark object, probably a BH, of mass M_BH = 6^{+3}_{-2} x 10^8 M_sun. However, anisotropic models fit the data without a BH if the ratio of radial to azimuthal dispersions is ~ 2 at 1". Therefore this is a less strong BH detection than the ones in M31, M32, and NGC 3115. A 6 x 10^8 M_sun BH is 9 % of the mass M_bulge in stars; even if M_BH is smaller than the isotropic value, M_BH/M_bulge is likely to be unusually large. Double nuclei are a puzzle because the dynamical friction timescales for self-gravitating star clusters in orbit around each other are short. Since both M31 and NGC 4486B contain central dark objects, our results support models in which the survival of double nuclei is connected with the presence of a BH (e. g., Tremaine 1995, AJ, 110, 628).

Optical photometry is presented for the quadruple gravitational lens PG1115+080. A preliminary reduction of data taken from November 1995 to June 1996 gives component ``C'' leading component ``B'' by 23.7+/-3.4 days and components ``A1'' and ``A2'' by 9.4 days. A range of models has been fit to the image positions, none of which gives an adequate fit. The best fitting and most physically plausible of these, taking the lensing galaxy and the associated group of galaxies to be singular isothermal spheres, gives a Hubble constant of 42 km/s/Mpc for Omega=1, with an observational uncertainty of 14%, as computed from the B-C time delay measurement. Taking the lensing galaxy to have an approximately E5 isothermal mass distribution yields H0=64 km/sec/Mpc while taking the galaxy to be a point mass gives H0=84 km/sec/Mpc. The former gives a particularly bad fit to the position of the lensing galaxy, while the latter is inconsistent with measurements of nearby galaxy rotation curves. Constraints on these and other possible models are expected to improve with planned HST observations.

We analyze the galaxy population and dynamics of the galaxy cluster Abell 576 using 281 redshifts (230 new), R band CCD galaxy photometry over a 2h^{-1} Mpc by 2h^{-1} Mpc region centered on the cluster, an Einstein IPC image, and an Einstein MPC spectrum. The cluster galaxies with emission lines in their spectra have a larger velocity dispersion and are significantly less clustered on this 2h^{-1} Mpc scale than galaxies without emission lines. We show that excluding the emission line galaxies decreases the velocity dispersion by 18% and the virial mass estimate by a factor of two. The central cluster region contains a non--emission galaxy population and an intracluster medium which are significantly cooler (sigma_{core}=387^{+250}_{-105}~km/s and T_X=1.6^{+0.4}_{-0.3}~keV at 90\% confidence) than the global populations (sigma=977^{+124}_{-96}~km/s for the non--emission population and T_X>4~keV at 90% confidence). Because (1) the low dispersion population is no more luminous than the global population and (2) the evidence for a cooling flow is weak, we suggest that the cluster core may contain the remnants of a lower mass subcluster. We examine the cluster mass, baryon fraction and luminosity function. The cluster virial mass varies significantly depending on the galaxy sample used. Consistency between the hydrostatic and virial estimators can be achieved if (1) the gas temperature at r\sim1h^{-1}~Mpc is T_X\sim8keV (the best fit value) and (2) several velocity outliers are excluded from the virial calculation. Although the best fit Schechter function parameters and the ratio of galaxy to gas mass are typical of other clusters, the baryon fraction is relatively low. Using a lower limit to the binding mass, we show that the gas mass fraction is \le3h^{-3/2}% and the baryon fraction is \le6%.

We have non-parametrically determined the luminosity density profiles and their logarithmic slopes for 42 early-type galaxies observed with HST. Assuming that the isodensity contours are spheroidal, then the luminosity density is uniquely determined from the surface brightness data through the Abel equation. For nearly all the galaxies in our sample, the logarithmic slope of the luminosity density measured at 0.1" (the innermost reliable measurement with the uncorrected HST) is significantly different from zero; i.e. most elliptical galaxies have cusps. There are only two galaxies for which an analytic core cannot be excluded. The distribution of logarithmic slopes at 0.1" appears to be bimodal, confirming the conclusion of Lauer et al. (1995) that early-type galaxies can be divided into two types based on their surface-brightness profiles; i.e., those with cuspy cores and those whose steep power-law profiles continue essentially unchanged in to the resolution limit. The peaks in the slope distribution occur at -0.8 and -1.9. More than half of the galaxies have slopes steeper than -1.0. Taken together with the recent theoretical work of Merritt & Fridman, these results suggest that many (and maybe most) elliptical galaxies are either nearly axisymmetric or spherical near the center, or slowly evolve due to the influence of stochastic orbits.

We present a set of structural parameters for the central parts of 57 early-type galaxies observed with the Planetary Camera of the Hubble Space Telescope. These parameters are based on a new empirical law that successfully characterizes the centers of early type galaxies. This empirical law assumes that the surface brightness profile is a combination of two power laws with different slopes gamma and beta for the inner and outer regions. Conventional structural parameters such as core radius and central surface brightness are replaced by break radius r_b, where the transition between power-law slopes takes place, and surface brightness mu_b at that radius. An additional parameter alpha describes the sharpness of the break. The structural parameters are derived using a chi-squared minimization process applied to the mean surface brightness profiles. The resulting model profiles generally give very good agreement to the observed profiles out to the radius of 10 arcseconds imaged by the Planetary Camera. Exceptions include galaxies which depart from pure power-laws at large radius, those with strong nuclear components, and galaxies partly obscured by dust. The uncertainties in the derived parameters are estimated using Monte-Carlo simulations which test the stability of solutions in the face of photon noise and the effects of the deconvolution process. The covariance of the structural parameters is examined by computing contours of constant chi squared in multi-dimensional parameter space.

We use N-body simulations to study the shape of mean cluster density and velocity profiles formed via gravitational instability. The dependence of the final structure on both cosmology and initial density field is examined, using a grid of cosmologies and scale-free initial power spectra P\propto k^n. For each model, we stack clusters to define an average density profile in the non-linear regime. The profiles are well fit by a power law over 99% of the cluster volume, with a clear trend toward steeper slopes with both increasing n and decreasing Omega_o. For models with a Omega_o = 0.2, the profile slopes are consistently higher than those for Omega-1.0. Cluster density profiles are thus potentially useful cosmological diagnostics. We find no evidence for a constant density core in any of the models, although the density profiles do tend to flatten at small radii. Much of the flattening is due to the force softening required by the simulations, and an attempt is made to recover the unsoftened profiles assuming angular momentum invariance. The recovered profiles in the Omega=1 cosmologies are consistent with a pure power law up to the highest density contrasts (10^6) accessible with our resolution. The low density models show significant deviations from a power law above density contrasts \sim 10^5. We interpret this curvature as reflecting the non scale-invariant nature of the background cosmology in these models.