The galaxy catalogs generated from low-resolution emission line surveys often contain both foreground and background interlopers due to line misidentification, which can bias the cosmological parameter estimation. In this paper, we present a method for correcting the interloper bias by using the joint-analysis of auto- and cross-power spectra of the main and the interloper samples. In particular, we can measure the interloper fractions from the cross-correlation between the interlopers and survey galaxies, because the true cross-correlation must be negligibly small. The estimated interloper fractions, in turn, remove the interloper bias in the cosmological parameter estimation. For example, in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) low-redshift ($z<0.5$) [O II] $\lambda3727${\AA} emitters contaminate high-redshift ($1.9<z<3.5$) Lyman-$\alpha$ line emitters. We demonstrate that the joint-analysis method yields a high signal-to-noise ratio measurement of the interloper fractions while only marginally increasing the uncertainties in the cosmological parameters relative to the case without interlopers. We also show the same is true for the high-latitude spectroscopic survey of Wide-Field Infrared Survey Telescope (WFIRST) mission where contamination occurs between the Balmer-$\alpha$ line emitters at lower redshifts ($1.1<z<1.9$) and Oxygen ([O III] $\lambda5007${\AA}) line emitters at higher redshifts ($1.7<z<2.8$).

Supermassive black holes are located at the center of most, if not all, massive galaxies. They follow close correlations with global properties of their host galaxies (scaling relations), and are thought to play a crucial role in galaxy evolution. Yet, we lack a complete understanding of fundamental aspects of their growth across cosmic time. In particular, we still do not understand: (1) whether black holes or their host galaxies grow faster and (2) what is the maximum mass that black holes can reach. The high angular resolution capability and sensitivity of 30-m class telescopes will revolutionize our understanding of the extreme end of the black hole and galaxy mass scale. With such facilities, we will be able to dynamically measure masses of the largest black holes and characterize galaxy properties out to redshift $z \sim 1.5$. Together with the evolution of black hole-galaxy scaling relations since $z \sim 1.5$, the maximum mass black hole will shed light on the main channels of black hole growth.

We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist.

By utilizing large-scale graph analytic tools implemented in the modern Big Data platform, Apache Spark, we investigate the topological structure of gravitational clustering in five different universes produced by cosmological $N$-body simulations with varying parameters: (1) a WMAP 5-year compatible $\Lambda$CDM cosmology, (2) two different dark energy equation of state variants, and (3) two different cosmic matter density variants. For the Big Data calculations, we use a custom build of stand-alone Spark/Hadoop cluster at Korea Institute for Advanced Study (KIAS) and Dataproc Compute Engine in Google Cloud Platform (GCP) with the sample size ranging from 7 millions to 200 millions. We find that among the many possible graph-topological measures, three simple ones: (1) the average of number of neighbors (the so-called average vertex degree) $\alpha$, (2) closed-to-connected triple fraction (the so-called transitivity) $\tau_\Delta$, and (3) the cumulative number density $n_{s\ge5}$ of subcomponents with connected component size $s \ge 5$, can effectively discriminate among the five model universes. Since these graph-topological measures are in direct relation with the usual $n$-points correlation functions of the cosmic density field, graph-topological statistics powered by Big Data computational infrastructure opens a new, intuitive, and computationally efficient window into the dark Universe.

We investigate the spatial distribution of Lyman alpha emitting galaxies (LAEs) at $z \approx 2.67$, selected from the NOAO Deep Wide-Field Survey (NDWFS), using two-point statistics and topological diagnostics adopted from network science. We measure the clustering length, $r_0 \approx 4 h^{-1}$ Mpc, and the bias, $b_{LAE} = 2.2^{+0.2}_{-0.1}$. Fitting the clustering with halo occupation distribution (HOD) models results in two disparate possibilities: (1) where the fraction of central galaxies is $<$1% in halos of mass $>10^{12}$$M_\odot$; and (2) where the fraction is $\approx$20%. We refer to these two scenarios as the `Dusty Core Scenario' for Model#1 since most of central galaxies in massive halos are dead in Ly$\alpha$ emission, and the `Pristine Core Scenario' for Model#2 since the central galaxies are bright in Ly$\alpha$ emission. Traditional two-point statistics cannot distinguish between these disparate models given the current data sets. To overcome this degeneracy, we generate mock catalogs for each HOD model using a high resolution $N$-body simulation and adopt a network statistics approach, which provides excellent topological diagnostics for galaxy point distributions. We find three topological anomalies from the spatial distribution of observed LAEs, which are not reproduced by the HOD mocks. We find that Model#2 matches better all network statistics than Model#1, suggesting that the central galaxies in $> 10^{12} h^{-1} M_\odot$ halos at $z \approx 2.67$ need to be less dusty to be bright as LAEs, potentially implying some replenishing channels of pristine gas such as the cold mode accretion.

Present-day clusters are massive halos containing mostly quiescent galaxies, while distant protoclusters are extended structures containing numerous star-forming galaxies. We investigate the implications of this fundamental change in a cosmological context using a set of N-body simulations and semi-analytic models. We find that the fraction of the cosmic volume occupied by all (proto)clusters increases by nearly three orders of magnitude from z=0 to z=7. We show that (proto)cluster galaxies are an important, and even dominant population at high redshift, as their expected contribution to the cosmic star-formation rate density rises (from 1% at z=0) to 20% at z=2 and 50% at z=10. Protoclusters thus provide a significant fraction of the cosmic ionizing photons, and may have been crucial in driving the timing and topology of cosmic reionization. Internally, the average history of cluster formation can be described by three distinct phases: at z~10-5, galaxy growth in protoclusters proceeded in an inside-out manner, with centrally dominant halos that are among the most active regions in the Universe; at z~5-1.5, rapid star formation occurred within the entire 10-20 Mpc structures, forming most of their present-day stellar mass; at z<~1.5, violent gravitational collapse drove these stellar contents into single cluster halos, largely erasing the details of cluster galaxy formation due to relaxation and virialization. Our results motivate observations of distant protoclusters in order to understand the rapid, extended stellar growth during Cosmic Noon, and their connection to reionization during Cosmic Dawn.

Dedicated photometric and spectroscopic surveys have provided unambiguous evidence for a strong stellar mass-size evolution of galaxies within the last 10 Gyr. The likely progenitors of today's most massive galaxies are remarkably small, disky, passive and have already assembled much of their stellar mass at redshift z=2. An in-depth analysis of these objects, however, is currently not feasible due to the lack of high-quality, spatially-resolved photometric and spectroscopic data. In this paper, we present a sample of nearby compact elliptical galaxies (CEGs), which bear resemblance to the massive and quiescent galaxy population at earlier times. Hubble Space Telescope (HST) and wide-field integral field unit (IFU) data have been obtained, and are used to constrain orbit-based dynamical models and stellar population synthesis (SPS) fits, to unravel their structural and dynamical properties. We first show that our galaxies are outliers in the present-day stellar mass-size relation. They are, however, consistent with the mass-size relation of compact, massive and quiescent galaxies at redshift z=2. The compact sizes of our nearby galaxies imply high central stellar mass surface densities, which are also in agreement with the massive galaxy population at higher redshift, hinting at strong dissipational processes during their formation. Corroborating evidence for a largely passive evolution within the last 10 Gyr is provided by their orbital distribution as well as their stellar populations, which are difficult to reconcile with a very active (major) merging history. This all supports that we can use nearby CEGs as local analogues of the high-redshift, massive and quiescent galaxy population, thus providing additional constraints for models of galaxy formation and evolution.

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.

We present a Bayesian approach to the redshift classification of emission-line galaxies when only a single emission line is detected spectroscopically. We consider the case of surveys for high-redshift Lyman-alpha-emitting galaxies (LAEs), which have traditionally been classified via an inferred rest-frame equivalent width (EW) greater than 20 angstrom. Our Bayesian method relies on known prior probabilities in measured emission-line luminosity functions and equivalent width distributions for the galaxy populations, and returns the probability that an object in question is an LAE given the characteristics observed. This approach will be directly relevant for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which seeks to classify ~10^6 emission-line galaxies into LAEs and low-redshift [O II] emitters. For a simulated HETDEX catalog with realistic measurement noise, our Bayesian method recovers 86% of LAEs missed by the traditional EW > 20 angstrom cutoff over 2 < z < 3, outperforming the EW cut in both contamination and incompleteness. This is due to the method's ability to trade off between the two types of binary classification error by adjusting the stringency of the probability requirement for classifying an observed object as an LAE. In our simulations of HETDEX, this method reduces the uncertainty in cosmological distance measurements by 14% with respect to the EW cut, equivalent to recovering 29% more cosmological information. Rather than using binary object labels, this method enables the use of classification probabilities in large-scale structure analyses. It can be applied to narrowband emission-line surveys as well as upcoming large spectroscopic surveys including Euclid and WFIRST.

We present post-cryogenic Spitzer imaging at 3.6 and 4.5 micron with the Infrared Array Camera (IRAC) of the Spitzer/HETDEX Exploratory Large-Area (SHELA) survey. SHELA covers $\sim$deg$^2$ of the Sloan Digital Sky Survey "Stripe 82" region, and falls within the footprints of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) and the Dark Energy Survey. The HETDEX blind R $\sim$ 800 spectroscopy will produce $\sim$ 200,000 redshifts from the Lyman-$\alpha$ emission for galaxies in the range 1.9 < z < 3.5, and an additional $\sim$200,000 redshifts from the [OII] emission for galaxies at z < 0.5. When combined with deep ugriz images from the Dark Energy Camera, K-band images from NEWFIRM, and other ancillary data, the IRAC photometry from Spitzer will enable a broad range of scientific studies of the relationship between structure formation, galaxy stellar mass, halo mass, AGN, and environment over a co-moving volume of $\sim$0.5 Gpc$^3$ at 1.9 < z < 3.5. Here, we discuss the properties of the SHELA IRAC dataset, including the data acquisition, reduction, validation, and source catalogs. Our tests show the images and catalogs are 80% (50%) complete to limiting magnitudes of 22.0 (22.6) AB mag in the detection image, which is constructed from the weighted sum of the IRAC 3.6 and 4.5 micron images. The catalogs reach limiting sensitivities of 1.1 $\mu$Jy at both 3.6 and 4.5 micron (1$\sigma$, for R=2 arcsec circular apertures). As a demonstration of science, we present IRAC number counts, examples of highly temporally variable sources, and galaxy surface density profiles of rich galaxy clusters. In the spirit of Spitzer Exploratory programs we provide all images and catalogs as part of the publication.

(abridged) The large-scale distribution of galaxies is generally analyzed using the two-point correlation function. However, this statistic does not capture the topology of the distribution, and it is necessary to resort to higher order correlations to break degeneracies. We demonstrate that an alternate approach using network analysis can discriminate between topologically different distributions that have similar two-point correlations. We investigate two galaxy point distributions, one produced by a cosmological simulation and the other by a L\'evy walk. For the cosmological simulation, we adopt the redshift $z = 0.58$ slice from Illustris (Vogelsberger et al. 2014A) and select galaxies with stellar masses greater than $10^8$$M_\odot$. The two point correlation function of these simulated galaxies follows a single power-law, $\xi(r) \sim r^{-1.5}$. Then, we generate L\'evy walks matching the correlation function and abundance with the simulated galaxies. We find that, while the two simulated galaxy point distributions have the same abundance and two point correlation function, their spatial distributions are very different; most prominently, \emph{filamentary structures}, absent in L\'evy fractals. To quantify these missing topologies, we adopt network analysis tools and measure diameter, giant component, and transitivity from networks built by a conventional friends-of-friends recipe with various linking lengths. Unlike the abundance and two point correlation function, these network quantities reveal a clear separation between the two simulated distributions; therefore, the galaxy distribution simulated by Illustris is not a L\'evy fractal quantitatively. We find that the described network quantities offer an efficient tool for discriminating topologies and for comparing observed and theoretical distributions.

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.

We investigate the compact, early-type galaxy NGC 1281 with integral field unit observations to map the stellar LOSVD out to 5 effective radii and construct orbit-based dynamical models to constrain its dark and luminous matter content. Under the assumption of mass-follows-light, the H-band stellar mass-to-light ratio (M/L) is {\Upsilon} = 2.7(+-0.1) {\Upsilon}_{sun}, higher than expected from our stellar population synthesis fits with either a canonical Kroupa ({\Upsilon} = 1.3 {\Upsilon}_{sun}) or Salpeter ({\Upsilon} = 1.7 {\Upsilon}_{sun}) stellar initial mass function. Such models also cannot reproduce the details of the LOSVD. Models with a dark halo recover the kinematics well and indicate that NGC 1281 is dark matter dominated, making up ~ 90 per cent of the total enclosed mass within the kinematic bounds. Parameterised as a spherical NFW profile, the dark halo mass is 11.5 < log(M_{DM}/M_{sun}) < 11.8 and the stellar M/L is 0.6 < {\Upsilon} < 1.1. However, this stellar M/L is lower than predicted by its old stellar population. Moreover, the halo mass within the kinematic extent is ten times larger than expected based on {\Lambda}CDM predictions, and an extrapolation yields cluster sized dark halo masses. Adopting {\Upsilon} = 1.7 {\Upsilon}_{sun} yields more moderate dark halo virial masses, but these models fit the kinematics worse. A non-NFW model might solve the discrepancy between the unphysical consequences of the best-fitting dynamical models and models based on more reasonable assumptions for the dark halo and stellar mass-to-light ratio, which are disfavoured according to our parameter estimation.

Recently, two independent groups found very different results when measuring the central velocity dispersion of the galactic globular cluster NGC 6388 with different methods. While L\"utzgendorf et al. (2011) found a rising profile and a high central velocity dispersion (23.3 km/s), measurements obtained by Lanzoni et al. (2013) showed a value 40% lower. The value of the central velocity dispersion has a serious impact on the mass and possible presence of an intermediate-mass black hole at the center of NGC 6388. We use a photometric catalog of NGC 6388 to create a simulated SINFONI and ARGUS dataset. The construction of the IFU data cube is done with different observing conditions reproducing the conditions reported for the original observations as closely as possible. In addition, we produce an N-body realization of a 10^6 M_SUN stellar cluster with the same photometric properties as NGC 6388 to account for unresolved stars. We find that the individual radial velocities, i.e. the measurements from the simulated SINFONI data, are systematically biased towards lower velocity dispersions. The reason is that due to the wings in the point spread function the velocities get biased towards the mean cluster velocity. This study shows that even with AO supported observations, individual radial velocities in crowded fields are likely to be biased. The ARGUS observations do not show this kind of bias but were found to have larger uncertainties than previously obtained. We find a bias towards higher velocity dispersions in the ARGUS pointing when fixing the extreme velocities of the three brightest stars but find those variations are within the determined uncertainties. We rerun Jeans models and fit the kinematic profile with the new uncertainties. This yields a BH mass of M_BH = (2.8 +- 0.4) x 10^4 M_SUN and M/L ratio M/L = (1.6 +- 0.1) M_SUN/L_SUN, consistent with our previous results.

We present a dynamical analysis to infer the structural parameters and properties of the two nearby, compact, high velocity dispersion galaxies MRK1216 & NGC1277. Combining deep HST imaging, wide-field IFU stellar kinematics, and complementary long-slit spectroscopic data out to 3 R_e, we construct orbit-based models to constrain their black hole masses, dark matter content and stellar mass-to-light ratios. We obtain a black hole mass of log(Mbh/Msun) = 10.1(+0.1/-0.2) for NGC1277 and an upper limit of log(Mbh/Msun) = 10.0 for MRK1216, within 99.7 per cent confidence. The stellar mass-to-light ratios span a range of Upsilon_V = 6.5(+1.5/-1.5) in NGC1277 and Upsilon_H = 1.8(+0.5/-0.8) in MRK1216 and are in good agreement with SSP models of a single power-law Salpeter IMF. Even though our models do not place strong constraints on the dark halo parameters, they suggest that dark matter is a necessary ingredient in MRK1216, with a dark matter contribution of 22(+30/-20) per cent to the total mass budget within 1 R_e. NGC1277, on the other hand, can be reproduced without the need for a dark halo, and a maximal dark matter fraction of 13 per cent within the same radial extent. In addition, we investigate the orbital structures of both galaxies, which are rotationally supported and consistent with photometric multi-S\'ersic decompositions, indicating that these compact objects do not host classical, non-rotating bulges formed during recent (z <= 2) dissipative events or through violent relaxation. Finally, both MRK 1216 and NGC 1277 are anisotropic, with a global anisotropy parameter delta of 0.33 and 0.58, respectively. While MRK 1216 follows the trend of fast-rotating, oblate galaxies with a flattened velocity dispersion tensor in the meridional plane of the order of beta_z = delta, NGC 1277 is highly tangentially anisotropic and seems to belong kinematically to a distinct class of objects.

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.

Galaxy proto-clusters at z >~ 2 provide a direct probe of the rapid mass assembly and galaxy growth of present day massive clusters. Because of the need of precise galaxy redshifts for density mapping and the prevalence of star formation before quenching, nearly all the proto-clusters known to date were confirmed by spectroscopy of galaxies with strong emission lines. Therefore, large emission-line galaxy surveys provide an efficient way to identify proto-clusters directly. Here we report the discovery of a large-scale structure at z = 2.44 in the HETDEX Pilot Survey. On a scale of a few tens of Mpc comoving, this structure shows a complex overdensity of Lya emitters (LAE), which coincides with broad-band selected galaxies in the COSMOS/UltraVISTA photometric and zCOSMOS spectroscopic catalogs, as well as overdensities of intergalactic gas revealed in the Lya absorption maps of Lee et al. (2014). We construct mock LAE catalogs to predict the cosmic evolution of this structure. We find that such an overdensity should have already broken away from the Hubble flow, and part of the structure will collapse to form a galaxy cluster with 10^14.5 +- 0.4 M_sun by z = 0. The structure contains a higher median stellar mass of broad-band selected galaxies, a boost of extended Lya nebulae, and a marginal excess of active galactic nuclei relative to the field, supporting a scenario of accelerated galaxy evolution in cluster progenitors. Based on the correlation between galaxy overdensity and the z = 0 descendant halo mass calibrated in the simulation, we predict that several hundred 1.9 < z < 3.5 proto-clusters with z = 0 mass of > 10^14.5 M_sun will be discovered in the 8.5 Gpc^3 of space surveyed by the Hobby Eberly Telescope Dark Energy Experiment.

We present new wide-field photometry and spectroscopy of the globular clusters (GCs) around NGC 4649 (M60), the third brightest galaxy in the Virgo cluster. Imaging of NGC 4649 was assembled from a recently-obtained HST/ACS mosaic, and new Subaru/Suprime-Cam and archival CFHT/MegaCam data. About 1200 sources were followed up spectroscopically using combined observations from three multi-object spectrographs: Keck/DEIMOS, Gemini/GMOS and MMT/Hectospec. We confirm 431 unique GCs belonging to NGC 4649, a factor of 3.5 larger than previous datasets and with a factor of 3 improvement in velocity precision. We confirm significant GC colour bimodality and find that the red GCs are more centrally concentrated, while the blue GCs are more spatially extended. We infer negative GC colour gradients in the innermost 20 kpc and flat gradients out to large radii. Rotation is detected along the galaxy major axis for all tracers: blue GCs, red GCs, galaxy stars and planetary nebulae. We compare the observed properties of NGC 4649 with galaxy formation models. We find that formation via a major merger between two gas-poor galaxies, followed by satellite accretion, can consistently reproduce the observations of NGC 4649 at different radii. We find no strong evidence to support an interaction between NGC 4649 and the neighbouring spiral galaxy NGC 4647. We identify interesting GC kinematic features in our data, such as counter-rotating subgroups and bumpy kinematic profiles, which encode more clues about the formation history of NGC 4649.

We have conducted an optical long-slit spectroscopic survey of 1022 galaxies using the 10m Hobby-Eberly Telescope (HET) at McDonald Observatory. The main goal of the HET Massive Galaxy Survey (HETMGS) is to find nearby galaxies that are suitable for black hole mass measurements. In order to measure accurately the black hole mass, one should kinematically resolve the region where the black hole dominates the gravitational potential. For most galaxies, this region is much less than an arcsecond. Thus, black hole masses are best measured in nearby galaxies with telescopes that obtain high-spatial resolution. The HETMGS focuses on those galaxies predicted to have the largest sphere-of-influence, based on published stellar velocity dispersions or the galaxy fundamental plane. To ensure coverage over galaxy types, the survey targets those galaxies across a face-on projection of the fundamental plane. We present the sample selection and resulting data products from the long-slit observations, including central stellar kinematics and emission line ratios. The full dataset, including spectra and resolved kinematics, is available online. Additionally, we show that the current crop of black hole masses are highly biased towards dense galaxies and that especially large disks and low dispersion galaxies are under-represented. This survey provides the necessary groundwork for future systematic black hole mass measurement campaigns.

The study of intermediate-mass black holes (IMBHs) is a young and promising field of research. Formed by runaway collisions of massive stars in young and dense stellar clusters, intermediate-mass black holes could still be present in the centers of globular clusters, today. Our group investigated the presence of intermediate-mass black holes for a sample of 10 Galactic globular clusters. We measured the inner kinematic profiles with integral-field spectroscopy and determined masses or upper limits of central black holes in each cluster. In combination with literature data we further studied the positions of our results on known black-hole scaling relations (such as M_bh - sigma) and found a similar but flatter correlation for IMBHs. Applying cluster evolution codes, the change in the slope could be explained with the stellar mass loss occurring in clusters in a tidal field over its life time. Furthermore, we present results from several numerical simulations on the topic of IMBHs and integral field units (IFUs). We ran N-body simulations of globular clusters containing IMBHs in a tidal field and studied their effects on mass-loss rates and remnant fractions and showed that an IMBH in the center prevents core collapse and ejects massive objects more rapidly. These simulations were further used to simulate IFU data cubes. For the specific case of NGC 6388 we simulated two different IFU techniques and found that velocity dispersion measurements from individual velocities are strongly biased towards lower values due to blends of neighbouring stars and background light. In addition, we use the Astrophysical Multipurpose Software Environment (AMUSE) to combine gravitational physics, stellar evolution and hydrodynamics to simulate the accretion of stellar winds onto a black hole.

The Hobby-Eberly Dark Energy Experiment pilot survey identified 284 [O II] 3727 emitting galaxies in a 169 square-arcminute field of sky in the redshift range 0 < z < 0.57. This line flux limited sample provides a bridge between studies in the local universe and higher-redshift [O II] surveys. We present an analysis of the star formation rates (SFRs) of these galaxies as a function of stellar mass as determined via spectral energy distribution fitting. The [O II] emitters fall on the "main sequence" of star-forming galaxies with SFR decreasing at lower masses and redshifts. However, the slope of our relation is flatter than that found for most other samples, a result of the metallicity dependence of the [O II] star formation rate indicator. The mass specific SFR is higher for lower mass objects, supporting the idea that massive galaxies formed more quickly and efficiently than their lower mass counterparts. This is confirmed by the fact that the equivalent widths of the [O II] emission lines trend smaller with larger stellar mass. Examination of the morphologies of the [O II] emitters reveals that their star formation is not a result of mergers, and the galaxies' half-light radii do not indicate evolution of physical sizes.

We present the results from a VLT/SINFONI and Keck/NIRSPEC near-infrared spectroscopic survey of 16 Lyman-alpha emitters (LAEs) at $z$ = 2.1 - 2.5 in the COSMOS and GOODS-N fields discovered from the HETDEX Pilot Survey. We detect rest-frame optical nebular lines (H$\alpha$ and/or [OIII]$\lambda$5007) for 10 of the LAEs and measure physical properties, including the star formation rate (SFR), gas-phase metallicity, gas-mass fraction, and Ly$\alpha$ velocity offset. We find that LAEs may lie below the mass-metallicity relation for continuum-selected star-forming galaxies at the same redshift. The LAEs all show velocity shifts of Ly$\alpha$ relative to the systemic redshift ranging between +85 and +296 km s$^{-1}$ with a mean of +180 km s$^{-1}$. This value is smaller than measured for continuum-selected star-forming galaxies at similar redshifts. The Ly$\alpha$ velocity offsets show a moderate correlation with the measured star formation rate (2.5$\sigma$), but no significant correlations are seen with the SFR surface density, specific SFR, stellar mass, or dynamical mass ($\lesssim$ 1.5$\sigma$). Exploring the role of dust, kinematics of the interstellar medium (ISM), and geometry on the escape of Ly$\alpha$ photons, we find no signature of selective quenching of resonantly scattered Ly$\alpha$ photons. However, we also find no evidence that a clumpy ISM is enhancing the Ly$\alpha$ equivalent width. Our results suggest that the low metallicity in LAEs may be responsible for yielding an environment with a low neutral hydrogen column density as well as less dust, easing the escape of Ly$\alpha$ photons over that in continuum-selected star-forming galaxies.

We present new constraints on the density profiles of dark matter (DM) halos in seven nearby dwarf galaxies from measurements of their integrated stellar light and gas kinematics. The gas kinematics of low mass galaxies frequently suggest that they contain constant density DM cores, while N-body simulations instead predict a cuspy profile. We present a data set of high resolution integral field spectroscopy on seven galaxies and measure the stellar and gas kinematics simultaneously. Using Jeans modeling on our full sample, we examine whether gas kinematics in general produce shallower density profiles than are derived from the stars. Although 2/7 galaxies show some localized differences in their rotation curves between the two tracers, estimates of the central logarithmic slope of the DM density profile, gamma, are generally robust. The mean and standard deviation of the logarithmic slope for the population are gamma=0.67+/-0.10 when measured in the stars and gamma=0.58+/-0.24 when measured in the gas. We also find that the halos are not under concentrated at the radii of half their maximum velocities. Finally, we search for correlations of the DM density profile with stellar velocity anisotropy and other baryonic properties. Two popular mechanisms to explain cored DM halos are an exotic DM component or feedback models that strongly couple the energy of supernovae into repeatedly driving out gas and dynamically heating the DM halos. We investigate correlations that may eventually be used to test models. We do not find a secondary parameter that strongly correlates with the central DM density slope, but we do find some weak correlations. Determining the importance of these correlations will require further model developments and larger observational samples. (Abridged)

We have obtained radial velocity measurements for 51 new globular clusters around the Sombrero galaxy. These measurements were obtained using spectroscopic observations from the AAOmega spectrograph on the Anglo-Australian Telescope and the Hydra spectrograph at WIYN. Combined with our own past measurements and velocity measurements obtained from the literature we have constructed a large database of radial velocities that contains a total of 360 confirmed globular clusters. Previous studies' analyses of the kinematics and mass profile of the Sombrero globular cluster system have been constrained to the inner ~9' (~24 kpc or ~5 effective radii), but our new measurements have increased the radial coverage of the data, allowing us to determine the kinematic properties of M104 out to ~15' (~41 kpc or ~9 effective radii). We use our set of radial velocities to study the GC system kinematics and to determine the mass profile and V-band mass-to-light profile of the galaxy. We find that the V-band mass-to-light ratio increases from 4.5 at the center to a value of 20.9 at 41 kpc (~9 effective radii or 15'), which implies that the dark matter halo extends to the edge of our available data set. We compare our mass profile at 20 kpc (~4 effective radii or ~7.4') to the mass computed from x-ray data and find good agreement. We also use our data to look for rotation in the globular cluster system as a whole, as well as in the red and blue subpopulations. We find no evidence for significant rotation in any of these samples.

We use broadband photometry extending from the rest-frame UV to the near-IR to fit the individual spectral energy distributions (SEDs) of 63 bright (L(Ly-alpha) > 10^43 ergs/s) Ly-alpha emitting galaxies (LAEs) in the redshift range 1.9 < z < 3.6. We find that these LAEs are quite heterogeneous, with stellar masses that span over three orders of magnitude, from 7.5 < log M < 10.5. Moreover, although most LAEs have small amounts of extinction, some high-mass objects have stellar reddenings as large as E(B-V) ~0.4. Interestingly, in dusty objects the optical depths for Ly-alpha and the UV continuum are always similar, indicating that Ly-alpha photons are not undergoing many scatters before escaping their galaxy. In contrast, the ratio of optical depths in low-reddening systems can vary widely, illustrating the diverse nature of the systems. Finally, we show that in the star formation rate (SFR)-log mass diagram, our LAEs fall above the "main-sequence" defined by z ~ 3 continuum selected star-forming galaxies. In this respect, they are similar to sub-mm-selected galaxies, although most LAEs have much lower mass.

We have measured the line-of-sight velocity distribution from integrated stellar light at two points in the outer halo of M87 (NGC 4486), the second-rank galaxy in the Virgo Cluster. The data were taken at R = 480" ($\sim 41.5$ kpc) and R = 526" ($\sim 45.5$ kpc) along the SE major axis. The second moment for a non-parametric estimate of the full velocity distribution is $420 \pm 23$ km/s and $577 \pm 35$ km/s respectively. There is intriguing evidence in the velocity profiles for two kinematically distinct stellar components at the position of our pointing. Under this assumption we employ a two-Gaussian decomposition and find the primary Gaussian having rest velocities equal to M87 (consistent with zero rotation) and second moments of $383 \pm 32$ km/s and $446 \pm 43$ km/s respectively. The asymmetry seen in the velocity profiles suggests that the stellar halo of M87 is not in a relaxed state and confuses a clean dynamical interpretation. That said, either measurement (full or two component model) shows a rising velocity dispersion at large radii, consistent with previous integrated light measurements, yet significantly higher than globular cluster measurements at comparable radial positions. These integrated light measurements at large radii, and the stark contrast they make to the measurements of other kinematic tracers, highlight the rich kinematic complexity of environments like the center of the Virgo Cluster and the need for caution when interpreting kinematic measurements from various dynamical tracers.

We present the first set of maps and band-merged catalog from the Herschel Stripe 82 Survey (HerS). Observations at 250, 350, and 500 micron were taken with the Spectral and Photometric Imaging Receiver (SPIRE) instrument aboard the Herschel Space Observatory. HerS covers 79 deg$^2$ along the SDSS Stripe 82 to a depth of 13.0, 12.9, and 14.8 mJy beam$^{-1}$ (including confusion) at 250, 350, and 500 micron, respectively. HerS was designed to measure correlations with external tracers of the dark matter density field --- either point-like (i.e., galaxies selected from radio to X-ray) or extended (i.e., clusters and gravitational lensing) --- in order to measure the bias and redshift distribution of intensities of infrared-emitting dusty star-forming galaxies and AGN. By locating HeRS in Stripe 82, we maximize the overlap with available and upcoming cosmological surveys. The band-merged catalog contains 3.3x10$^4$ sources detected at a significance of >3 $\sigma$ (including confusion noise). The maps and catalog are available at http://www.astro.caltech.edu/hers/

For a complete picture of galaxy cluster formation, it is important that we start probing the early epoch of z~2-7 during which clusters and their galaxies first began to form. Because the study of these so-called "proto-clusters" is currently limited by small number statistics, widely varying selection techniques and assumptions, we have performed a systematic study of cluster formation utilizing cosmological simulations. We use the Millennium Simulations to track the evolution of dark matter and galaxies in ~3,000 clusters from the earliest times to z=0. We define an effective radius R_e for proto-clusters and characterize their growth in size and mass. We show that the progenitor regions of galaxy clusters (M>10^14 M_sun/h) can already be identified at least up to z~5, provided that the galaxy overdensities, delta_gal, are measured on a sufficiently large scale (R_e~5-10 cMpc). We present the overdensities in matter, DM halos, and galaxies as functions of present-day cluster mass, redshift, bias, and window size that can be used to interpret the structures found in real surveys. We derive the probability that a structure having a delta_gal, defined by a set of observational selection criteria, is indeed a proto-cluster, and show how their z=0 masses can already be estimated long before virialization. Galaxy overdensity profiles as a function of radius are presented. We further show how their projected surface overdensities decrease as the uncertainties in redshift measurements increase. We provide a table of proto-cluster candidates selected from the literature, and discuss their properties in the light of our simulations predictions. This work provides the general framework that will allow us to extend the study of cluster formation out to much higher redshifts using the large number of proto-clusters that are expected to be discovered in, e.g., the upcoming HETDEX and HSC surveys.

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.

To demonstrate the feasibility of studying the epoch of massive galaxy cluster formation in a more systematic manner using current and future galaxy surveys, we report the discovery of a large sample of proto-cluster candidates in the 1.62 deg^2 COSMOS/UltraVISTA field traced by optical/IR selected galaxies using photometric redshifts. By comparing properly smoothed 3D galaxy density maps of the observations and a set of matched simulations incorporating the dominant observational effects (galaxy selection and photometric redshift uncertainties), we first confirm that the observed ~15 comoving Mpc scale galaxy clustering is consistent with LCDM models. Using further the relation between high-z overdensity and the present day cluster mass calibrated in these matched simulations, we found 36 candidate structures at 1.6<z<3.1, showing overdensities consistent with the progenitors of M_z=0 ~10^15 M_sun clusters. Taking into account the significant upward scattering of lower mass structures, the probabilities for the candidates to have at least M_z=0 ~10^14 M_sun are ~70%. For each structure, about 15%-40% of photometric galaxy candidates are expected to be true proto-cluster members that will merge into a cluster-scale halo by z=0. With solely photometric redshifts, we successfully rediscover two spectroscopically confirmed structures in this field, suggesting that our algorithm is robust. This work generates a large sample of uniformly-selected proto-cluster candidates, providing rich targets for spectroscopic follow-up and subsequent studies of cluster formation. Meanwhile, it demonstrates the potential for probing early cluster formation with upcoming redshift surveys such as the Hobby-Eberly Telescope Dark Energy Experiment and the Subaru Prime Focus Spectrograph survey.

Survey observations of the three-dimensional locations of galaxies are a powerful approach to measure the distribution of matter in the universe, which can be used to learn about the nature of dark energy, physics of inflation, neutrino masses, etc. A competitive survey, however, requires a large volume (e.g., Vsurvey is roughly 10 Gpc3) to be covered, and thus tends to be expensive. A "sparse sampling" method offers a more affordable solution to this problem: within a survey footprint covering a given survey volume, Vsurvey, we observe only a fraction of the volume. The distribution of observed regions should be chosen such that their separation is smaller than the length scale corresponding to the wavenumber of interest. Then one can recover the power spectrum of galaxies with precision expected for a survey covering a volume of Vsurvey (rather than the volume of the sum of observed regions) with the number density of galaxies given by the total number of observed galaxies divided by Vsurvey (rather than the number density of galaxies within an observed region). We find that regularly-spaced sampling yields an unbiased power spectrum with no window function effect, and deviations from regularly-spaced sampling, which are unavoidable in realistic surveys, introduce calculable window function effects and increase the uncertainties of the recovered power spectrum. While we discuss the sparse sampling method within the context of the forthcoming Hobby-Eberly Telescope Dark Energy Experiment, the method is general and can be applied to other galaxy surveys.

Using the VIMOS Integral Field Unit (IFU) spectrograph on the Very Large Telescope (VLT), we have spatially mapped the kinematic properties of 10 nearby Brightest Cluster Galaxies (BCGs) and 4 BCG companion galaxies located within a redshift of $z=0.1$. In the hierarchical formation model, these massive galaxies $(10^{10.5} M_{\odot} < M_{dyn} < 10^{11.9} M_{\odot})$ are expected to undergo more mergers than lower mass galaxies, and simulations show that dry minor mergers can remove angular momentum. We test whether BCGs have low angular momenta by using the $\lambda_{Re}$ parameter developed by the SAURON and ATLAS\textsuperscript{3D} teams and combine our kinematics with Sloan Digital Sky Survey (SDSS) photometry to analyze the BCGs' merger status. We find that 30% (3/10) of the BCGs and 100% of the companion galaxies (4/4) are fast rotators as defined by the ATLAS\textsuperscript{3D} criteria. Our fastest rotating BCG has a $\lambda_{Re}=0.35\pm0.05$. We increase the number of BCGs analyzed from 1 in the combined SAURON and ATLAS\textsuperscript{3D} surveys to 11 BCGs total and find that above $M_{dyn}\sim11.5 M_{\odot}$, virtually all galaxies regardless of environment are slow rotators. To search for signs of recent merging, we analyze the photometry of each system and use the $G-M_{20}$ selection criteria to identify mergers. We find that $40\pm20$% of our BCGs are currently undergoing or have recently undergone a merger (within 0.2 Gyrs). Surprisingly, we find no correlation between galaxies with high angular momentum and morphological signatures of merging.

The discovery of cosmic acceleration has inspired ambitious experimental and observational efforts to understand its origin. Many of these take the form of large astronomical surveys, sometimes using new, special-purpose instrumentation, and in some cases entirely new facilities. This Report summarizes some of the major ongoing and planned dark energy experiments, focusing on those in which the U.S. community has a leading or significant supporting role. It provides background information for other Reports from the "Dark Energy and CMB" working group of the APS Division of Particles and Fields long-term planning exercise. We provide 1-2 page summaries of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) and its SDSS-IV successor eBOSS, the Dark Energy Survey (DES), the Subaru Hyper-Suprime Camera (HSC) and Prime Focus Spectrograph (PFS), the Hobby-Eberley Telescope Dark Energy Experiment (HEDTEX), the Dark Energy Spectroscopic Instrument (DESI), the Large Synpotic Survey Telescope (LSST), and the Euclid and AFTA-WFIRST (Wide Field Infrared Survey Telescope) space missions. Over the next one to two decades, these experiments seek to improve the precision of measurements of the cosmic expansion history and cosmic structure growth by factors of 10-100, while simultaneously increasing their redshift range and tightening control and cross-checks of sytematic effects. These measurements will dramatically improve the empirical constraints on theoretical explanations of cosmic acceleration and could lead to breakthroughs in our understanding of the nature of gravity and spacetime or the forms of energy present in the universe.

Using a newly-developed modeling technique, we present orbit-based dynamical models of the Carina, Draco, Fornax, Sculptor, and Sextans dwarf spheroidal (dSph) galaxies. These models calculate the dark matter profiles non-parametrically without requiring any assumptions to be made about their profile shapes. By lifting this restriction we discover a host of dark matter profiles in the dSphs that are different from the typical profiles suggested by both theorists and observers. However, when we scale these profiles appropriately and plot them on a common axis they appear to follow an approximate r^-1 power law with considerable scatter.

We study the radial dependence in stellar populations of 33 nearby early-type galaxies with central stellar velocity dispersions sigma* > 150 km/s. We measure stellar population properties in composite spectra, and use ratios of these composites to highlight the largest spectral changes as a function of radius. Based on stellar population modeling, the typical star at 2 R_e is old (~10 Gyr), relatively metal poor ([Fe/H] -0.5), and alpha-enhanced ([Mg/Fe]~0.3). The stars were made rapidly at z~1.5-2 in shallow potential wells. Declining radial gradients in [C/Fe], which follow [Fe/H], also arise from rapid star formation timescales due to declining carbon yields from low-metallicity massive stars. In contrast, [N/Fe] remains high at large radius. Stars at large radius have different abundance ratio patterns from stars in the center of any present-day galaxy, but are similar to Milky Way thick disk stars. Our observations are thus consistent with a picture in which the stellar outskirts are built up through minor mergers with disky galaxies whose star formation is truncated early (z~1.5-2).

We present new results on the spectrally resolved Lyman-alpha (LyA) emission of three LyA emitting field galaxies at z~2.4 with high LyA equivalent width (>100 Angstroms) and LyA luminosity (~10^43 erg/s). At 120 km/s (FWHM) spectral resolution, the prominent double-peaked LyA profile straddles the systemic velocity, where the velocity zero-point is determined from spectroscopy of the galaxies' rest-frame optical nebular emission lines. The average velocity offset from systemic of the stronger redshifted emission component for our sample is 176 km/s while the average total separation between the redshifted and main blueshifted emission components is 380 km/s. These measurements are a factor of ~2 smaller than for UV continuum-selected galaxies that show LyA in emission with lower LyA equivalent width. We compare our LyA spectra to the predicted line profiles of a spherical "expanding shell" LyA radiative transfer grid that models large-scale galaxy outflows. Specifically blueward of the systemic velocity where two galaxies show a weak, highly blueshifted (by ~1000 km/s) tertiary emission peak, the model line profiles are a relatively poor representation of the observed spectra. Since the neutral gas column density has a dominant influence over the shape of the LyA line profile, we caution against equating the observed LyA velocity offset with a physical outflow velocity, especially at lower spectral resolution where the unresolved LyA velocity offset is a convoluted function of several degenerate parameters. Referring to rest-frame ultraviolet and optical Hubble Space Telescope imaging, we find that galaxy-galaxy interactions may play an important role in inducing a starburst that results in copious LyA emission, as well as perturbing the gas distribution and velocity field which have strong influence over the LyA emission line profile.

For galaxies hosting supermassive black holes (SMBHs), it has been observed that the mass of the central black hole (M_BH) tightly correlates with the effective or central velocity dispersion (sigma) of the host galaxy. The origin of this M_BH - sigma scaling relation is assumed to lie in the merging history of the galaxies but many open questions about its origin and the behavior in different mass ranges still need to be addressed. The goal of this work is to study the black-hole scaling relations for low black-hole masses, where the regime of intermediate-mass black holes (IMBHs) in globular clusters (GCs) is entered. We collect all existing reports of dynamical black-hole measurements in globular clusters, providing black-hole masses or upper limits for 14 candidates. We plot the black-hole masses versus different cluster parameters including total mass, velocity dispersion, concentration and half-mass radius. We search for trends and test the correlations in order to quantify their significance using a set of different statistical approaches. For correlations showing a large significance we perform a linear fit, accounting for uncertainties and upper limits. We find a clear correlation between the mass of the IMBH and the velocity dispersion of the globular cluster. As expected, the total mass of the globular cluster then also correlates with the mass of the IMBH. While the slope of the M_BH - sigma correlation differs strongly from the one observed for SMBHs, the other scaling relations M_BH - M_TOT, and M_BH - L are similar to the correlations in galaxies. Significant correlations of black-hole mass with other cluster properties were not found in the present sample.

We present an analysis of the luminosities and equivalent widths of the 284 z < 0.56 [O II]-emitting galaxies found in the 169 square arcmin pilot survey for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). By combining emission-line fluxes obtained from the Mitchell spectrograph on the McDonald 2.7-m telescope with deep broadband photometry from archival data, we derive each galaxy's de-reddened [O II] 3727 luminosity and calculate its total star formation rate. We show that over the last ~5 Gyr of cosmic time there has been substantial evolution in the [O II] emission-line luminosity function, with L* decreasing by ~0.6 +/-0.2 dex in the observed function, and by ~0.9 +/-0.2 dex in the de-reddened relation. Accompanying this decline is a significant shift in the distribution of [O II] equivalent widths, with the fraction of high equivalent-width emitters declining dramatically with time. Overall, the data imply that the relative intensity of star formation within galaxies has decreased over the past ~5 Gyr, and that the star formation rate density of the universe has declined by a factor of ~2.5 between z ~ 0.5 and z ~ 0. These observations represent the first [O II]-based star formation rate density measurements in this redshift range, and foreshadow the advancements which will be generated by the main HETDEX survey.

Intermediate-mass black holes (IMBHs, 10^2-10^5 M_sun) fill the gap between stellar-mass black holes and supermassive black holes (SMBHs). Simulations have shown that IMBHs may form in dense star clusters, and therefore may still be present in these smaller stellar systems. We investigate the Galactic globular cluster NGC 5286 for indications of a central IMBH using spectroscopic data from VLT/FLAMES, velocity measurements from the Rutgers Fabry Perot at CTIO, and photometric data from HST. We run analytic spherical and axisymmetric Jeans models with different central black-hole masses, anisotropy, mass-to-light ratio, and inclination. Further, we compare the data to a grid of N-body simulations without tidal field. Additionally, we use one N-body simulation to check the results of the spherical Jeans models for the total cluster mass. Both the Jeans models and the N-body simulations favor the presence of a central black hole in NGC 5286 and our detection is at the 1- to 1.5-sigma level. From the spherical Jeans models we obtain a best fit with black-hole mass M_BH=(1.5+-1.0)x10^3 M_sun. The error is the 68% confidence limit from Monte Carlo simulations. Axisymmetric models give a consistent result. The best fitting N-body model is found with a black hole of 0.9% of the total cluster mass (4.38+-0.18)x10^5 M_sun, which results in an IMBH mass of M_BH=(3.9+-2.0)x10^3 M_sun. Jeans models give lower values for the total cluster mass. Our test of the Jeans models with N-body simulation data shows that this discrepancy has two reasons: The influence of a radially varying M/L profile, and underestimation of the velocity dispersion as the measurements are limited to bright stars. We conclude that detection of IMBHs in Galactic globular clusters remains a challenging task unless their mass fractions are above those found for SMBHs in nearby galaxies. [abridged]

We have tested the effect of spatial gradients in stellar mass-to-light ratio (Y) on measurements of black hole masses (MBH) derived from stellar orbit superposition models. Such models construct a static gravitational potential for a galaxy and its central black hole, but typically assume spatially uniform Y. We have modeled three giant elliptical galaxies with gradients alpha = d(log Y)/d(log r) from -0.2 to +0.1. Color and line strength gradients suggest mildly negative alpha in these galaxies. Introducing a negative (positive) gradient in Y increases (decreases) the enclosed stellar mass near the center of the galaxy and leads to systematically smaller (larger) MBH measurements. For models with alpha = -0.2, the best-fit values of MBH are 28%, 27%, and 17% lower than the constant-Y case, in NGC 3842, NGC 6086, and NGC 7768, respectively. For alpha = +0.1, MBH are 14%, 22%, and 17% higher than the constant-Y case for the three respective galaxies. For NGC 3842 and NGC 6086, this bias is comparable to the statistical errors from individual modeling trials. At larger radii, negative (positive) gradients in Y cause the total stellar mass to decrease (increase) and the dark matter fraction within one effective radius to increase (decrease).

We present the survey design, data reduction, and spectral fitting pipeline for the VIRUS-P Exploration of Nearby Galaxies (VENGA). VENGA is an integral field spectroscopic survey, which maps the disks of 30 nearby spiral galaxies. Targets span a wide range in Hubble type, star formation activity, morphology, and inclination. The VENGA data-cubes have 5.6'' FWHM spatial resolution, ~5A FWHM spectral resolution, sample the 3600A-6800A range, and cover large areas typically sampling galaxies out to ~0.7 R_25. These data-cubes can be used to produce 2D maps of the star formation rate, dust extinction, electron density, stellar population parameters, the kinematics and chemical abundances of both stars and ionized gas, and other physical quantities derived from the fitting of the stellar spectrum and the measurement of nebular emission lines. To exemplify our methods and the quality of the data, we present the VENGA data-cube on the face-on Sc galaxy NGC 628 (a.k.a. M 74). The VENGA observations of NGC 628 are described, as well as the construction of the data-cube, our spectral fitting method, and the fitting of the stellar and ionized gas velocity fields. We also propose a new method to measure the inclination of nearly face-on systems based on the matching of the stellar and gas rotation curves using asymmetric drift corrections. VENGA will measure relevant physical parameters across different environments within these galaxies, allowing a series of studies on star formation, structure assembly, stellar populations, chemical evolution, galactic feedback, nuclear activity, and the properties of the interstellar medium in massive disk galaxies.

For a sample of nine Galactic globular clusters we measured the inner kinematic profiles with integral-field spectroscopy that we combined with existing outer kinematic measurements and HST luminosity profiles. With this information we are able to detect the crucial rise in the velocity-dispersion profile which indicates the presence of a central black hole. In addition, N-body simulations compared to our data will give us a deeper insight in the properties of clusters with black holes and stronger selection criteria for further studies. For the first time, we obtain a homogeneous sample of globular cluster integral- field spectroscopy which allows a direct comparison between clusters with and without an intermediate-mass black hole.

Recent hydrostatic X-ray studies of the hot interstellar medium (ISM) in early-type galaxies underestimate the gravitating mass as compared to stellar dynamics, implying modest, but significant deviations from exact hydrostatic equilibrium. We present a method for combining X-ray measurements and stellar dynamical constraints in the context of Bayesian statistics that allows the radial distribution of the implied nonthermal pressure or bulk motions in the hot ISM to be constrained. We demonstrate the accuracy of the method with hydrodynamical simulations tailored to produce a realistic galaxy model. Applying the method to the nearby elliptical galaxy NGC4649, we find a significant but subdominant nonthermal pressure fraction (0.27+/-0.06) in the central (<5 kpc) part of the galaxy, similar to the level of deviations from hydrostatic equilibrium expected in galaxy clusters. Plausible sources of systematic error, if important, may reduce this fraction. This would imply >360 km/s random turbulence or a magnetic field B=(39+/-6)(n_e/0.1 cm^{-3})^{0.59+/-0.09} muG, whereas gas rotation alone is unlikely to explain the detailed nonthermal profile. Future observations with Astro-H will allow turbulence or gas rotation at this level to be detected.

We measure the radial profile of the 12CO(1-0) to H_2 conversion factor (Xco) in NGC 628. The H\alpha emission from the VENGA integral field spectroscopy is used to map the star formation rate surface density (\Sigma_{SFR}). We estimate the molecular gas surface density (\Sigma_{H2}) from \Sigma_{SFR} by inverting the molecular star formation law (SFL), and compare it to the CO intensity to measure Xco. We study the impact of systematic uncertainties by changing the slope of the SFL, using different SFR tracers (H\alpha vs. far-UV plus 24\mu m), and CO maps from different telescopes (single-dish and interferometers). The observed Xco profile is robust against these systematics, drops by a factor of 2 from R~7 kpc to the center of the galaxy, and is well fit by a gradient \Delta log(Xco)=0.06\pm0.02 dex kpc^-1. We study how changes in Xco follow changes in metallicity, gas density, and ionization parameter. Theoretical models show that the gradient in Xco can be explained by a combination of decreasing metallicity, and decreasing \Sigma_{H2} with radius. Photoelectric heating from the local UV radiation field appears to contribute to the decrease of Xco in higher density regions. Our results show that galactic environment plays an important role at setting the physical conditions in star forming regions, in particular the chemistry of carbon in molecular complexes, and the radiative transfer of CO emission. We caution against adopting a single Xco value when large changes in gas surface density or metallicity are present.

The formation of supermassive black holes at high redshift still remains a puzzle to astronomers. Their growth becomes reasonable only when starting from a massive seed black hole with mass of the order of 10^2 - 10^5 M_SUN. Intermediate-mass black holes (IMBHs) are therefore an important field of research. Especially the possibility of finding them in the centers of globular clusters has recently drawn attention. The search for IMBHs in the centers of globular clusters could therefore shed light on the process of black-hole formation and cluster evolution. We are investigating six galactic globular clusters for the presence of an IMBH at their centers. Based on their kinematic and photometric properties, we selected the globular clusters NGC 1851, NGC 1904 (M79), NGC 5694, NGC 5824, NGC 6093 (M80) and NGC 6266 (M62). We use integral field spectroscopy in order to obtain the central velocity-dispersion profile of each cluster. We compute the cluster photometric center and the surface brightness profile using HST data. After combining these datasets we compare them to analytic Jeans models. We use varying M/L_V profiles for clusters with enough data points in order to reproduce their kinematic profiles in an optimal way. Finally, we vary the mass of the central black hole and test whether the cluster is better fitted with or without an IMBH. We present the statistical significance, including upper limits, of the black-hole mass for each cluster. NGC 1904 and NGC 6266 provide the highest significance for a black hole. Jeans models in combination with a M/L_V profile obtained from N-body simulations (in the case of NGC 6266) predict a central black hole of M_BH = (3 +- 1) x 10^3 M_SUN for NGC 1904 and M_BH = (2 +- 1) x 10^3 M_SUN for NGC 6266. Furthermore, we discuss the possible influence of dark remnants and mass segregation at the center of the cluster on the detection of an IMBH.

All massive galaxies likely have supermassive black holes at their centers, and the masses of the black holes are known to correlate with properties of the host galaxy bulge component. Several explanations have been proposed for the existence of these locally-established empirical relationships; they include the non-causal, statistical process of galaxy-galaxy merging, direct feedback between the black hole and its host galaxy, or galaxy-galaxy merging and the subsequent violent relaxation and dissipation. The empirical scaling relations are thus important for distinguishing between various theoretical models of galaxy evolution, and they further form the basis for all black hole mass measurements at large distances. In particular, observations have shown that the mass of the black hole is typically 0.1% of the stellar bulge mass of the galaxy. The small galaxy NGC4486B currently has the largest published fraction of its mass in a black hole at 11%. Here we report observations of the stellar kinematics of NGC 1277, which is a compact, disky galaxy with a mass of 1.2 x 10^11 Msun. From the data, we determine that the mass of the central black hole is 1.7 x 10^10 Msun, or 59% its bulge mass. Five other compact galaxies have properties similar to NGC 1277 and therefore may also contain over-sized black holes. It is not yet known if these galaxies represent a tail of a distribution, or if disk-dominated galaxies fail to follow the normal black hole mass scaling relations.

We introduce a novel implementation of orbit-based (or Schwarzschild) modeling that allows dark matter density profiles to be calculated non-parametrically in nearby galaxies. Our models require no assumptions to be made about velocity anisotropy or the dark matter profile. The technique can be applied to any dispersion-supported stellar system, and we demonstrate its use by studying the Local Group dwarf spheroidal (dSph) galaxy Draco. We use existing kinematic data at larger radii and also present 12 new radial velocities within the central 13 pc obtained with the VIRUS-W integral field spectrograph on the 2.7m telescope at McDonald Observatory. Our non-parametric Schwarzschild models find strong evidence that the dark matter profile in Draco is cuspy for 20 < r < 700 pc. The profile for r > 20 pc is well-fit by a power law with slope \alpha=-1.0 +/- 0.2, consistent with predictions from Cold Dark Matter (CDM) simulations. Our models confirm that, despite its low baryon content relative to other dSphs, Draco lives in a massive halo.

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 report the detection of five high-velocity stars in the core of the globular cluster NGC 2808. The stars lie on the the red giant branch and show total velocities between 40 and 45 km/s. For a core velocity dispersion sigma_c = 13.4 km/s, this corresponds to up to 3.4 sigma_c. These velocities are close to the estimated escape velocity (~ 50 km/s) and suggest an ejection from the core. Two of these stars have been confirmed in our recent integral field spectroscopy data and we will discuss them in more detail here. These two red giants are located at a projected distance of ~ 0.3 pc from the center. According to their positions on the color magnitude diagram, both stars are cluster members. We investigate several possible origins for the high velocities of the stars and conceivable ejection mechanisms. Since the velocities are close to the escape velocity, it is not obvious whether the stars are bound or unbound to the cluster. We therefore consider both cases in our analysis. We perform numerical simulations of three-body dynamical encounters between binaries and single stars and compare the resulting velocity distributions of escapers with the velocities of our stars. We compare the predictions for a single dynamical encounter with a compact object with those of a sequence of two-body encounters due to relaxation. If the stars are unbound, the encounter must have taken place recently, when the stars were already in the giant phase. After including binary fractions and black-hole retention fractions, projection effects, and detection probabilities from Monte-Carlo simulations, we estimate the expected numbers of detections for all the different scenarios. Based on these numbers, we conclude that the most likely scenario is that the stars are bound and were accelerated by a single encounter between a binary of main-sequence stars and a ~ 10 M_sun black hole.

Globular clusters are an excellent laboratory for stellar population and dynamical research. Recent studies have shown that these stellar systems are not as simple as previously assumed. With multiple stellar populations as well as outer rotation and mass segregation they turn out to exhibit high complexity. This includes intermediate-mass black holes which are proposed to sit at the centers of some massive globular clusters. Today's high angular resolution ground based spectrographs allow velocity-dispersion measurements at a spatial resolution comparable to the radius of influence for plausible IMBH masses, and to detect changes in the inner velocity-dispersion profile. Together with high quality photometric data from HST, it is possible to constrain black-hole masses by their kinematic signatures. We determine the central velocity-dispersion profile of the globular cluster NGC 2808 using VLT/FLAMES spectroscopy. In combination with HST/ACS data our goal is to probe whether this massive cluster hosts an intermediate-mass black hole at its center and constrain the cluster mass to light ratio as well as its total mass. We derive a velocity-dispersion profile from integral field spectroscopy in the center and Fabry Perot data for larger radii. High resolution HST data are used to obtain the surface brightness profile. Together, these data sets are compared to dynamical models with varying parameters such as mass to light ratio profiles and black-hole masses. Using analytical Jeans models in combination with variable M/L profiles from N-body simulations we find that the best fit model is a no black hole solution. After applying various Monte Carlo simulations to estimate the uncertainties, we derive an upper limit of the back hole mass of M_BH < 1 x 10^4 M_SUN (with 95 % confidence limits) and a global mass-to-light ratio of M/L_V = (2.1 +- 0.2) M_SUN/L_SUN.

We use the Mitchell Spectrograph (formerly VIRUS-P) on the McDonald Observatory 2.7m Harlan J. Smith Telescope to search for the chemical signatures of massive elliptical galaxy assembly. The Mitchell Spectrograph is an integral-field spectrograph with a uniquely wide field of view (107x107 sq arcsec), allowing us to achieve remarkably high signal-to-noise ratios of ~20-70 per pixel in radial bins of 2-2.5 times the effective radii of the eight galaxies in our sample. Focusing on a sample of massive elliptical galaxies with stellar velocity dispersions sigma* > 150 km/s, we study the radial dependence in the equivalent widths (EWs) of key metal absorption lines. By twice the effective radius, the Mgb EWs have dropped by ~50%, and only a weak correlation between sigma* and Mgb EW remains. The Mgb EWs at large radii are comparable to those seen in the centers of elliptical galaxies that are approximately an order of magnitude less massive. We find that the well-known metallicity gradients often observed within an effective radius continue smoothly to 2.5R_e, while the abundance ratio gradients remain flat. Much like the halo of the Milky Way, the stellar halos of our galaxies have low metallicities and high alpha-abundance ratios, as expected for very old stars formed in small stellar systems. Our observations support a picture in which the outer parts of massive elliptical galaxies are built by the accretion of much smaller systems whose star formation history was truncated at early times.

This paper presents further results from our spectroscopic study of the globular cluster (GC) system of the group elliptical NGC 3923. From observations made with the GMOS instrument on the Gemini South telescope, an additional 50 GC and Ultra Compact Dwarf (UCD) candidates have been spectroscopically confirmed as members of the NGC 3923 system. When the recessional velocities of these GCs are combined with the 29 GC velocities reported previously, a total sample of 79 GC/UCD velocities is produced. This sample extends to over 6 arcmin (>6 Re \sim30 kpc) from the centre of NGC 3923, and is used to study the dynamics of the GC system and the dark matter content of NGC 3923. It is found that the GC system of NGC 3923 displays no appreciable rotation, and that the projected velocity dispersion is constant with radius within the uncertainties. The velocity dispersion profiles of the integrated light and GC system of NGC 3923 are indistinguishable over the region in which they overlap. We find some evidence that the diffuse light and GCs of NGC 3923 have radially biased orbits within \sim130". The application of axisymmetric orbit-based models to the GC and integrated light velocity dispersion profiles demonstrates that a significant increase in the mass-to-light ratio (from M/Lv = 8 to 26) at large galactocentric radii is required to explain these observations. We therefore confirm the presence of a dark matter halo in NGC 3923. We find that dark matter comprises 17.5% of the mass within 1 Re, 41.2% within 2 Re, and 75.6% within the radius of our last kinematic tracer at 6.9 Re. The total dynamical mass within this radius is found to be 1.5 x 10^12 solar masses. In common with other studies of large ellipticals, we find that our derived dynamical mass profile is consistently higher than that derived by X-ray observations, by a factor of around 2.

This work aims to study the distribution of luminous and dark matter in Coma early-type galaxies. Dynamical masses obtained under the assumption that mass follows light do not match with the masses of strong gravitational lens systems of similar velocity dispersions. Instead, dynamical fits with dark matter halos are in good agreement with lensing results. We derive mass-to-light ratios of the stellar populations from Lick absorption line indices, reproducing well the observed galaxy colours. Even in dynamical models with dark matter halos the amount of mass that follows the light increases more rapidly with galaxy velocity dispersion than expected for a constant stellar initial mass function (IMF). While galaxies around sigma ~ 200 km/s are consistent with a Kroupa IMF, the same IMF underpredicts luminous dynamical masses of galaxies with sigma ~ 300 km/s by a factor of two and more. A systematic variation of the stellar IMF with galaxy velocity dispersion could explain this trend with a Salpeter IMF for the most massive galaxies. If the IMF is instead constant, then some of the dark matter in high velocity dispersion galaxies must follow a spatial distribution very similar to that of the light. A combination of both, a varying IMF and a component of dark matter that follows the light is possible as well. For a subsample of galaxies with old stellar populations we show that the tilt in the fundamental plane can be explained by systematic variations of the total (stellar + dark) mass inside the effective radius. We tested commonly used mass estimator formulae, finding them accurate at the 20-30% level.

We construct axisymmetric Schwarzschild models to measure the mass profile of the local group dwarf galaxy Fornax. These models require no assumptions to be made about the orbital anisotropy of the stars, as is the case for commonly used Jeans models. We test a variety of parameterizations of dark matter density profiles and find cored models with uniform density rho_c = (1.6 +/- 0.1) x 10^-2 M_sun pc^-3 fit significantly better than the cuspy halos predicted by cold dark matter simulations. We also construct models with an intermediate-mass black hole, but are unable to make a detection. We place a 1-sigma upper limit on the mass of a potential intermediate-mass black hole at M_BH < 3.2 x 10^4 M_sun.

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.

Supermassive black holes (SMBHs) are fundamental keys to understand the formation and evolution of their host galaxies. However, the formation and growth of SMBHs are not yet well understood. One of the proposed formation scenarios is the growth of SMBHs from seed intermediate-mass black holes (IMBHs, 10^2 to 10^5 M_{\odot}) formed in star clusters. In this context, and also with respect to the low mass end of the M-sigma relation for galaxies, globular clusters are in a mass range that make them ideal systems to look for IMBHs. Among Galactic star clusters, the massive cluster $\omega$ Centauri is a special target due to its central high velocity dispersion and also its multiple stellar populations. We study the central structure and dynamics of the star cluster $\omega$ Centauri to examine whether an IMBH is necessary to explain the observed velocity dispersion and surface brightness profiles. We perform direct N-body simulations to follow the dynamical evolution of $\omega$ Centauri. The simulations are compared to the most recent data-sets in order to explain the present-day conditions of the cluster and to constrain the initial conditions leading to the observed profiles. We find that starting from isotropic spherical multi-mass King models and within our canonical assumptions, a model with a central IMBH mass of 2% of the cluster stellar mass, i.e. a 5x10^4 M_{\odot} IMBH, provides a satisfactory fit to both the observed shallow cusp in surface brightness and the continuous rise towards the center of the radial velocity dispersion profile. In our isotropic spherical models, the predicted proper motion dispersion for the best-fit model is the same as the radial velocity dispersion one. (abridged)

We study the mass distribution in the late-type dwarf galaxy NGC 2976 through stellar kinematics obtained with the VIRUS-P integral-field spectrograph and anisotropic Jeans models as a test of cosmological simulations and baryonic processes that putatively alter small-scale structure. Previous measurements of the H-alpha emission-line kinematics have determined that the dark matter halo of NGC 2976 is most consistent with a cored density profile. We find that the stellar kinematics are best fit with a cuspy halo. Cored dark matter halo fits are only consistent with the stellar kinematics if the stellar mass-to-light ratio is significantly larger than that derived from stellar population synthesis, while the best-fitting cuspy model has no such conflict. The inferred mass distribution from a harmonic decomposition of the gaseous kinematics is inconsistent with that of the stellar kinematics. This difference is likely due to the gas disk not meeting the assumptions that underlie the analysis such as no pressure support, a constant kinematic axis, and planar orbits. By relaxing some of these assumptions, in particular the form of the kinematic axis with radius, the gas-derived solution can be made consistent with the stellar kinematic models. A strong kinematic twist in the gas of NGC 2976's center suggests caution, and we advance the mass model based on the stellar kinematics as more reliable. The analysis of this first galaxy shows promising evidence that dark matter halos in late-type dwarfs may in fact be more consistent with cuspy dark matter distributions than earlier work has claimed.

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.

Intermediate-mass black holes (IMBHs) are of interest in a wide range of astrophysical fields. In particular, the possibility of finding them at the centers of globular clusters has recently drawn attention. IMBHs became detectable since the quality of observational data sets, particularly those obtained with HST and with high resolution ground based spectrographs, advanced to the point where it is possible to measure velocity dispersions at a spatial resolution comparable to the size of the gravitational sphere of influence for plausible IMBH masses. We present results from ground based VLT/FLAMES spectroscopy in combination with HST data for the globular cluster NGC 6388. The aim of this work is to probe whether this massive cluster hosts an intermediate-mass black hole at its center and to compare the results with the expected value predicted by the $M_{\bullet} - \sigma$ scaling relation. The spectroscopic data, containing integral field unit measurements, provide kinematic signatures in the center of the cluster while the photometric data give information of the stellar density. Together, these data sets are compared to dynamical models and present evidence of an additional compact dark mass at the center: a black hole. Using analytical Jeans models in combination with various Monte Carlo simulations to estimate the errors, we derive (with 68% confidence limits) a best fit black-hole mass of $ (17 \pm 9) \times 10^3 M_{\odot}$ and a global mass-to-light ratio of $M/L_V = (1.6 \pm 0.3) \ M_{\odot}/L_{\odot}$.

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.

This paper presents deep high quality photometry of globular cluster (GC) systems belonging to five early-type galaxies covering a range of mass and environment. Photometric data were obtained with the Gemini North and Gemini South telescopes in the filter passbands g', r', and i'. The combination of these filters with good seeing conditions allows an excellent separation between GC candidates and unresolved field objects. Bimodal GC colour distributions are found in all five galaxies. Most of the GC systems appear bimodal even in the (g' -r') vs (r' -i') plane. A population of resolved/marginally resolved GC and Ultra Compact Dwarf candidates was found in all the galaxies. A search for the so-called "blue tilt" in the colour-magnitude diagrams reveals that NGC 4649 clearly shows that phenomenon although no conclusive evidence was found for the other galaxies in the sample. This "blue tilt" translates into a mass-metallicity relation given by Z \propto M^0.28\pm0.03 . This dependence was found using a new empirical (g' -i') vs [Z/H] relation which relies on an homogeneous sample of GC colours and metallicities. This paper also explores the radial trends in both colour and surface density for the blue (metal-poor) and red (metal-rich) GC subpopulations. As usual, the red GCs show a steeper radial distribution than the blue ones. Evidence of galactocentric colour gradients is found in some of the GC systems, being more significant for the two S0 galaxies in the sample. Red GC subpopulations show similar colours and gradients to the galaxy halo stars in their inner region. A GC mean colour-galaxy luminosity relation, consistent with [Z/H] \propto L_B ^0.26\pm0.08, is present for the red GCs. An estimate of the total GC populations and specific frequency SN values is presented for NGC 3115, NGC 3379, NGC 3923 and NGC 4649.

We study the escape of Ly-alpha photons from Ly-alpha emitting galaxies (LAEs) and the overall galaxy population using a sample of 99 LAEs at 1.9<z<3.8 detected through integral-field spectroscopy of blank fields by the HETDEX Pilot Survey. For 89 LAEs showing counterparts in deep broad-band images we measure the rest-frame UV luminosity and the UV slope, which we use to estimate E(B-V) under the assumption of a constant intrinsic UV slope for LAEs. These two quantities are used to measure the dust-corrected star formation rate (SFR). A comparison between the observed Ly-alpha luminosity and that predicted by the dust-corrected SFR yields the Ly-alpha escape fraction. We also measure the Ly-alpha luminosity function. Integration of the luminosity function provides a measurement of the Ly-alpha luminosity density across our redshift range. We combine our data with that from other surveys at 0.3<z<7.7 to trace the evolution of the Ly-alpha luminosity density. We then compare it to that expected from the star-formation history of the universe in order to characterize the evolution of the Ly-alpha escape fraction of the overall galaxy population [abriged]

Stellar dynamical modeling is a powerful method to determine the mass of black holes in quiescent galaxies. However, in previous work the presence of a dark matter halo has been ignored in the modeling. Gebhardt & Thomas (2009) showed that accounting for a dark matter halo increased the black hole mass of the massive galaxy M87 by a factor of two. We used a sample of 12 galaxies to investigate the effect of accounting for a dark matter halo in the dynamical modeling in more detail, and also updated the masses using improved modeling. The sample of galaxies possesses Hubble Space Telescope and ground based observations of stellar kinematics. Their black hole masses have been presented before, but without including a dark matter halo in the models. Without a dark halo, we find a mean increase in the estimated mass of 1.5 for the whole sample compared to previous results. We attribute this change to using a more complete orbit library. When we include a dark matter halo, along with the updated models, we find an additional increase in black hole mass by a factor of 1.2 in the mean, much less than for M87. We attribute the smaller discrepancy in black hole mass to using data that better resolves the black hole's sphere of influence. We redetermined the M-sigma and M-L relationships using our updated black hole masses and found a slight increase in both normalization and intrinsic scatter.

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 2-D stellar kinematics of M87 out to R = 238" taken with the integral field spectrograph VIRUS-P. We run a large set of axisymmetric, orbit-based dynamical models and find clear evidence for a massive dark matter halo. While a logarithmic parameterization for the dark matter halo is preferred, we do not constrain the dark matter scale radius for an NFW profile and therefore cannot rule it out. Our best-fit logarithmic models return an enclosed dark matter fraction of 17.2 +/- 5.0 % within one effective radius (R_e ~ 100"), rising to 49.4 (+7.2,-8.8) % within 2 R_e. Existing SAURON data (R < 13"), and globular cluster kinematic data covering 145" < R < 540" complete the kinematic coverage to R = 47 kpc. At this radial distance the logarithmic dark halo comprises 85.3 (+2.5,-2.4) % of the total enclosed mass of 5.7^(+1.3)_(-0.9) X 10^(12) M_sun making M87 one of the most massive galaxies in the local universe. Our best-fit logarithmic dynamical models return a stellar mass-to-light ratio of 9.1^(+0.2)_(-0.2) (V-band), a dark halo circular velocity of 800^(+75)_(-25) kms, and a dark halo scale radius of 36^(+7)_(-3) kpc. The stellar M/L, assuming an NFW dark halo, is well constrained to 8.20^(+0.05)_(-0.10) (V-band). The stars in M87 are found to be radially anisotropic out to R ~ 0.5 R_e then isotropic or slightly tangentially anisotropic to our last stellar data point at R = 2.4 R_e where the anisotropy of the stars and globular clusters are in excellent agreement. The globular clusters then become radially anisotropic in the last two modeling bins at R = 3.4 R_e and R = 4.8 R_e. As one of the most massive galaxies in the local universe, constraints on both the mass distribution of M87 and anisotropy of its kinematic components strongly informs our theories of early-type galaxy formation and evolution in dense environments.

We present Keck/NIRSPEC spectroscopic observations of three Lyman alpha emitting galaxies (LAEs) at z ~ 2.3 discovered with the HETDEX pilot survey. We detect Halpha, [OIII], and Hbeta emission from two galaxies at z = 2.29 and 2.49, designated HPS194 and HPS256, respectively, representing the first detection of multiple rest-frame optical emission lines in galaxies at high-redshift selected on the basis of their Lyman alpha emission. The redshifts of the Lyman alpha emission from these galaxies are offset redward of the systemic redshifts by Delta_v = 162 +/- 37 (photometric) +/- 42 (systematic) km/s for HPS194, and Delta_v = 36 +/- 35 +/- 18 km/s for HPS256. An interpretation for HPS194 is that a large-scale outflow may be occurring in its interstellar medium. The emission line ratios imply that neither LAE hosts an active galactic nucleus. Using the upper limits on the [NII] emission we place meaningful constraints on the gas-phase metallicities in these two LAEs of Z < 0.17 and < 0.28 Zsol (1 sigma). Measuring the stellar masses of these objects via spectral energy distribution (SED) fitting (~ 10^10 and 6 x 10^8 Msol, respectively), we study the nature of LAEs in a mass-metallicity plane. At least one of these two LAEs appears to be more metal poor than continuum-selected star-forming galaxies at the same redshift and stellar mass, implying that objects exhibiting Lyman alpha emission may be systematically less chemically enriched than the general galaxy population. We use the SEDs of these two galaxies to show that neglecting the emission lines when fitting stellar population models to the observed photometry can result in overestimates of the population age by orders of magnitude, and the stellar mass by a factor of ~ 2. This effect is particularly important at z > 7, where similarly strong emission lines may masquerade in the photometry as a 4000 A break (abridged).

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.

We present a catalog of emission-line galaxies selected solely by their emission-line fluxes using a wide-field integral field spectrograph. This work is partially motivated as a pilot survey for the upcoming Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We describe the observations, reductions, detections, redshift classifications, line fluxes, and counterpart information for 397 emission-line galaxies detected over 169 sq.arcmin with a 3500-5800 Ang. bandpass under 5 Ang. full-width-half-maximum (FWHM) spectral resolution. The survey's best sensitivity for unresolved objects under photometric conditions is between 4-20 E-17 erg/s/sq.cm depending on the wavelength, and Ly-alpha luminosities between 3-6 E42 erg/s are detectable. This survey method complements narrowband and color-selection techniques in the search for high redshift galaxies with its different selection properties and large volume probed. The four survey fields within the COSMOS, GOODS-N, MUNICS, and XMM-LSS areas are rich with existing, complementary data. We find 104 galaxies via their high redshift Ly-alpha emission at 1.9<z<3.8, and the majority of the remainder objects are low redshift [OII]3727 emitters at z<0.56. The classification between low and high redshift objects depends on rest frame equivalent width, as well as other indicators, where available. Based on matches to X-ray catalogs, the active galactic nuclei (AGN) fraction amongst the Ly-alpha emitters (LAEs) is 6%. We also analyze the survey's completeness and contamination properties through simulations. We find five high-z, highly-significant, resolved objects with full-width-half-maximum sizes >44 sq.arcsec which appear to be extended Ly-alpha nebulae. We also find three high-z objects with rest frame Ly-alpha equivalent widths above the level believed to be achievable with normal star formation, EW(rest)>240 Ang.

The Galactic globular cluster omega Centauri is a prime candidate for hosting an intermediate mass black hole. Recent measurements lead to contradictory conclusions on this issue. We use VLT-FLAMES to obtain new integrated spectra for the central region of omega Centauri. We combine these data with existing measurements of the radial velocity dispersion profile taking into account a new derived center from kinematics and two different centers from the literature. The data support previous measurements performed for a smaller field of view and show a discrepancy with the results from a large proper motion data set. We see a rise in the radial velocity dispersion in the central region to 22.8+-1.2 km/s, which provides a strong sign for a central black hole. Isotropic dynamical models for omega Centauri imply black hole masses ranging from 3.0 to 5.2x10^4 solar masses depending on the center. The best-fitted mass is 4.7+-1.0x10^4 solar masses.

We present line-of-sight stellar velocity distributions of elliptical galaxy NGC 821 obtained to approximately 100 arcsec (over 2 effective radii) with long-slit spectroscopy from the Hobby-Eberly Telescope. Our measured stellar line-of-sight velocity distributions are larger than the planetary nebulae measurements at similar radii. We fit axisymmetric orbit-superposition models with a range of dark halo density profiles, including two-dimensional kinematics at smaller radii from SAURON data. Within our assumptions, the best-fitted model gives a total enclosed mass of 2.0x10^11 M_sun within 100 arcsec, with an accuracy of 2%; this mass is equally divided between halo and stars. At 1R_e the best-fitted dark matter halo accounts for 13% of the total mass in the galaxy. This dark halo is inconsistent with previous claims of little to no dark matter halo in this galaxy from planetary nebula measurements. We find that a power-law dark halo with a slope 0.1 is the best-fitted model; both the no dark halo and NFW models are worse fits at a greater than 99% confidence level. NGC 821 does not appear to have the expected dark halo density profile. The internal moments of the stellar velocity distribution show that the model with no dark halo is radially anisotropic at small radii and tangentially isotropic at large radii, while the best-fitted halo models are slightly radially anisotropic at all radii. We test the potential effects of model smoothing and find that there are no effects on our results within the errors. Finally, we run models using the planetary nebula kinematics and assuming our best-fitted halos and find that the planetary nebulae require radial orbits throughout the galaxy.

VENGA is a large-scale extragalactic IFU survey, which maps the bulges, bars and large parts of the outer disks of 32 nearby normal spiral galaxies. The targets are chosen to span a wide range in Hubble types, star formation activities, morphologies, and inclinations, at the same time of having vast available multi-wavelength coverage from the far-UV to the mid-IR, and available CO and 21cm mapping. The VENGA dataset will provide 2D maps of the SFR, stellar and gas kinematics, chemical abundances, ISM density and ionization states, dust extinction and stellar populations for these 32 galaxies. The uniqueness of the VIRUS-P large field of view permits these large-scale mappings to be performed. VENGA will allow us to correlate all these important quantities throughout the different environments present in galactic disks, allowing the conduction of a large number of studies in star formation, structure assembly, galactic feedback and ISM in galaxies.

We analyze six X-ray bright elliptical galaxies, observed with Chandra and XMM-Newton, and approximate their gravitational potentials by isothermal spheres phi(r)=v_c^2 ln(r) over a range of radii from ~0.5 to ~25 kpc. We then compare the circular speed v_c derived from X-ray data with the estimators available from optical data. In particular we discuss two simple and robust procedures for evaluating the circular speed of the galaxy using the observed optical surface brightness and the line-of-sight velocity dispersion profiles. The best fitting relation between the circular speeds derived from optical observations of stars and X-ray observations of hot gas is v_{c,opt}~ \eta * v_{c,X}, where \eta=1.10-1.15 (depending on the method), suggesting, albeit with large statistical and systematic uncertainties, that non-thermal pressure on average contributes ~20-30% of the gas thermal pressure.

We apply the axisymmetric orbit superposition modeling to estimate the mass of the supermassive black hole and dark matter halo profile of NGC 4649. We have included data sets from the Hubble Space Telescope, stellar, and globular cluster observations. Our modeling gives the black hole mass = 4.5 \pm 1.0 10^9 \Msun and M/L = 8.7 \pm 1.0 (or 8.0 \pm 0.9 after foreground Galactic extinction is corrected). We confirm the presence of a dark matter halo, but the stellar mass dominates inside the effective radius. The parameters of the dark halo are less constrained due to the sparse globular cluster data at large radii. We find that in NGC 4649 the dynamical mass profile from our modeling is consistently larger than that derived from the X-ray data over most of the radial range by roughly 60% to 80%. It implies that either some forms of non-thermal pressure need to be included, the assumed hydrostatic equilibrium may not be a good approximation in the X-ray modelings of NGC 4649, or our assumptions used in the dynamical models are biased. Our new black hole mass is about two times larger than the previous published value; the earlier model did not adequately sample the orbits required to match the large tangential anisotropy in the galaxy center. If we assume that there is no dark matter, the results on the black hole mass and M/L do not change significantly, which we attribute to the inclusion of HST spectra, the sparse globular cluster kinematics, and a diffuse dark matter halo. Without the HST data, the significance of the black hole detection is greatly reduced.

We investigate the relation between the star formation rate surface density (Sigma_SFR) and the mass surface density of gas (Sigma_gas) in NGC 5194. VIRUS-P integral field spectroscopy is used to measure H-alpha, H-beta [NII]6548,6584, and [SII]6717,6731 emission line fluxes for 735 regions ~170 pc in diameter, and derive extinction corrected Sigma_SFRs. HI 21cm and CO maps are used to measure the HI and H_2 gas surface density for each region. We present a new method for fitting the Star Formation Law (SFL), which includes the intrinsic scatter in the relation as a free parameter, allows the inclusion of non-detections, and is free of the systematics involved in performing linear correlations over incomplete data in logarithmic space. We use the [SII]/H-alpha ratio to separate the H-alpha flux contribution from the diffuse ionized gas (DIG). After removing the DIG, we measure a slope N=0.82+/-0.05, and an intrinsic scatter epsilon=0.43+/-0.02 dex for the molecular gas SFL. We also measure a typical depletion timescale tau~2 Gyr, in good agreement with Bigiel et al. (2008). The HI density shows no correlation with the SFR, and the total gas SFL closely follows the molecular gas SFL. We assess the validity of corrections applied in narrow-band H-alpha measurements to overcome issues related to continuum subtraction, underlying photospheric absorption, and contamination by the [NII] doublet. The disagreement with the super-linear molecular SFL measured by Kennicutt et al. (2007) is due to differences in the fitting method. Our results support a low and close to constant star formation efficiency (SFE = 1/tau) in the molecular ISM. The data shows excellent agreement with the model of Krumholz et al (2009). The large intrinsic scatter may imply the existence of other important parameters setting the SFR.

We present new mass measurements for the supermassive black holes (SMBHs) in the centres of three early-type galaxies. The gas pressure in the surrounding, hot interstellar medium (ISM) is measured through spatially resolved spectroscopy with the Chandra X-ray observatory, allowing the SMBH mass (Mbh) to be inferred directly under the hydrostatic approximation. This technique does not require calibration against other SMBH measurement methods and its accuracy depends only on the ISM being close to hydrostatic, which is supported by the smooth X-ray isophotes of the galaxies. Combined with results from our recent study of the elliptical galaxy NGC4649, this brings to four the number of galaxies with SMBHs measured in this way. Of these, three already have mass determinations from the kinematics of either the stars or a central gas disc, and hence join only a handful of galaxies with Mbh measured by more than one technique. We find good agreement between the different methods, providing support for the assumptions implicit in both the hydrostatic and the dynamical models. The stellar mass-to-light ratios for each galaxy inferred by our technique are in agreement with the predictions of stellar population synthesis models assuming a Kroupa initial mass function (IMF). This concurrence implies that no more than ~10-20% of the ISM pressure is nonthermal, unless there is a conspiracy between the shape of the IMF and nonthermal pressure. Finally, we compute Bondi accretion rates, finding that the two galaxies with the highest rates exhibit little evidence of X-ray cavities, suggesting that the correlation with the AGN jet power takes time to be established.

We model the dynamical structure of M87 (NGC4486) using high spatial resolution long-slit observations of stellar light in the central regions, two-dimensional stellar light kinematics out to half of the effective radius, and globular cluster velocities out to 8 effective radii. We simultaneously fit for four parameters, black hole mass, dark halo core radius, dark halo circular velocity, and stellar mass-to-light ratio. We find a black hole mass of 6.4(+-0.5)x10^9 Msun(the uncertainty is 68% confidence marginalized over the other parameters). The stellar M/L_V=6.3+-0.8. The best-fitted dark halo core radius is 14+-2 kpc, assuming a cored logarithmic potential. The best-fitted dark halo circular velocity is 715+-15 km/s. Our black hole mass is over a factor of 2 larger than previous stellar dynamical measures, and our derived stellar M/L ratio is 2 times lower than previous dynamical measures. When we do not include a dark halo, we measure a black hole mass and stellar M/L ratio that is consistent with previous measures, implying that the major difference is in the model assumptions. The stellar M/L ratio from our models is very similar to that derived from stellar population models of M87. The reason for the difference in the black hole mass is because we allow the M/L ratio to change with radius. The dark halo is degenerate with the stellar M/L ratio, which is subsequently degenerate with the black hole mass. We argue that dynamical models of galaxies that do not include the contribution from a dark halo may produce a biased result for the black hole mass. This bias is especially large for a galaxy with a shallow light profile such as M87, and may not be as severe in galaxies with steeper light profiles unless they have a large stellar population change with radius.

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.

We use oblate axisymmetric dynamical models including dark halos to determine the orbital structure of intermediate mass to massive Coma early-type galaxies. We find a large variety of orbital compositions. Averaged over all sample galaxies the unordered stellar kinetic energy in the azimuthal and the radial direction are of the same order, but they can differ by up to 40 percent in individual systems. In contrast, both for rotating and non-rotating galaxies the vertical kinetic energy is on average smaller than in the other two directions. This implies that even most of the rotating ellipticals are flattened by an anisotropy in the stellar velocity dispersions. Using three-integral axisymmetric toy models we show that flattening by stellar anisotropy maximises the entropy for a given density distribution. Collisionless disk merger remnants are radially anisotropic. The apparent lack of strong radial anisotropy in observed early-type galaxies implies that they may not have formed from mergers of disks unless the influence of dissipational processes was significant.

Axisymmetric, orbit-based dynamical models are used to derive dark matter scaling relations for Coma early-type galaxies. From faint to bright galaxies halo core-radii and asymptotic circular velocities increase. Compared to spirals of the same brightness, the majority of Coma early-types -- those with old stellar populations -- have similar halo core-radii but more than 2 times larger asymptotic halo velocities. The average dark matter density inside 2 reff decreases with increasing luminosity and is 6.8 times larger than in disk galaxies of the same B-band luminosity. Compared at the same stellar mass, dark matter densities in ellipticals are 13.5 times higher than in spirals. Different baryon concentrations in ellipticals and spirals cannot explain the higher dark matter density in ellipticals. Instead, the assembly redshift (1+z) of Coma early-type halos is likely about two times larger than of comparably bright spirals. Assuming that local spirals typically assemble at a redshift of one, the majority of bright Coma early-type galaxy halos must have formed around z = 2-3. For about half of our Coma galaxies the assembly redshifts match with constraints derived from stellar populations. We find dark matter densities and estimated assembly redshifts of our observed Coma galaxies in reasonable agreement with recent semi-analytic galaxy formation models.

We study a sample of 43 early-type galaxies, selected from the SDSS because they appeared to have velocity dispersion > 350 km/s. High-resolution photometry in the SDSS i passband using HRC-ACS on board the HST shows that just less than half of the sample is made up of superpositions of two or three galaxies, so the reported velocity dispersion is incorrect. The other half of the sample is made up of single objects with genuinely large velocity dispersions. None of these objects has sigma larger than 426 +- 30 km/s. These objects define rather different relations than the bulk of the early-type galaxy population: for their luminosities, they are the smallest, most massive and densest galaxies in the Universe. Although the slopes of the scaling relations they define are rather different from those of the bulk of the population, they lie approximately parallel to those of the bulk "at fixed sigma". These objects appear to be of two distinct types: the less luminous (M_r>-23) objects are rather flattened and extremely dense for their luminosities -- their properties suggest some amount of rotational support and merger histories with abnormally large amounts of gaseous dissipation. The more luminous objects (M_r<-23) tend to be round and to lie in or at the centers of clusters. Their properties are consistent with the hypothesis that they are BCGs. Models in which BCGs form from predominantly radial mergers having little angular momentum predict that they should be prolate. If viewed along the major axis, such objects would appear to have abnormally large sigma for their sizes, and to be abnormally round for their luminosities. This is true of the objects in our sample once we account for the fact that the most luminous galaxies (M_r<-23.5), and BCGs, become slightly less round with increasing luminosity.

The radio galaxy Fornax A (NGC 1316) is a prominent merger remnant in the outskirts of the Fornax cluster. Its giant radio lobes suggest the presence of a powerful AGN and thus a central supermassive black hole (SMBH). We present high-resolution adaptive optics assisted integral-field data of Fornax A, taken with SINFONI at the Very Large Telescope in the K band. We use axisymmetric orbit models to determine the mass of the SMBH in the centre of Fornax A. The three-dimensional nature of our data provides the possibility to directly test the consistency of the data with axisymmetry by modelling each of the four quadrants separately. According to our dynamical models, consistent SMBH masses and dynamical Ks band mass-to-light ratios are obtained for all quadrants, with <M_BH>=1.3x10^8 M_\odot (rms(M_BH)=0.4x10^8 Msun) and <M/L>=0.68 (rms(M/L)=0.03), confirming the assumption of axisymmetry. For the folded and averaged data we find M_BH=(1.5+0.75-0.8)x10^8 Msun and M/L=(0.65+0.075-0.05) (3-sigma errors). Thus the black-hole mass of Fornax A is consistent within the error with the Tremaine (2002) M-sigma relation, but is a factor ~4 smaller than expected from its bulge mass and the Marconi&Hunt (2003) relation.

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.

We present Hobby-Eberly Telescope (HET) observations for galaxies at redshift z < 0.3 from the Sloan Digital Sky Survey (SDSS) showing large velocity dispersions while appearing to be single galaxies in HST images. The high signal-to-noise HET spectra provide more definitive velocity dispersions. The maximum velocity dispersion we find is 444 km/s. Emission-line widths in QSOs indicate that black holes can exist with masses exceeding 5 billion solar masses, implying velocity dispersions greater than 500 km/s by the local black hole mass - velocity dispersion relationship. This suggests either that QSO black hole masses are overestimated or that the black hole - bulge relationship changes at high black hole mass. The latter option is consistent with evidence that the increase in velocity dispersion with luminosity levels off for the brightest elliptical galaxies.

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) will outfit the 10 m HET with a new wide field and an array of 150 integral-field spectrographs to survey a 420 sq. deg. area in the north Galactic cap. Each fiber-coupled unit spectrograph will cover 350-550 nm, simultaneously. This instrument, called VIRUS, will produce ~34,000 spectra per exposure, and will open up the emission-line universe to large surveys for the first time. The survey will detect 0.8 million Lyman-alpha emitting (LAE) galaxies with 1.9<z<3.5 and more than a million [OII] emitting galaxies with z<0.5. The 3-D map of LAE galaxies in 9 cubic Gpc volume will be used to measure the expansion history at this early epoch using baryonic acoustic oscillations and the shape of the power spectrum. The aim of HETDEX is to provide a direct detection of dark energy at z~3. The measurement will constrain the evolution of dark energy and will also provide 0.1%-level accuracy on the curvature of the Universe, ten times better than current. The prototype of the VIRUS unit spectrograph (VIRUS-P) is a powerful instrument in its own right. Used on the McDonald 2.7 m, it covers the largest area of any integral field spectrograph, and reaches wavelengths down to 340 nm. VIRUS-P is being used for a pilot survey to better measure the properties of LAE galaxies in support of HETDEX. We report initial results from this survey.

We present a Chandra study of the hot ISM in the giant elliptical galaxy NGC4649. In common with other group-centred ellipticals, its temperature profile rises with radius in the outer parts of the galaxy, from ~0.7keV at 2kpc to ~0.9keV by 20kpc. However, within the central ~2kpc the trend reverses and the temperature peaks at ~1.1keV within the innermost 200pc. Under the assumption of hydrostatic equilibrium, we demonstrate that the central temperature spike arises due to the gravitational influence of a quiescent central super-massive black hole. We constrain the black hole mass (MBH) to $(3.35^{+0.67}_{-0.95})\times 10^9$Msun (90% confidence), in good agreement with stellar kinematics measurements. This is the first direct measurement of MBH based on studies of hydrostatic X-ray emitting gas, which are sensitive to the most massive black holes, and is a crucial validation of both mass-determination techniques. This agreement clearly demonstrates the gas must be close to hydrostatic, even in the very centre of the galaxy, which is consistent with the lack of morphological disturbances in the X-ray image. NGC4649 is now one of only a handful of galaxies for which MBH has been measured by more than one method. At larger radii, we were able to decompose the gravitating mass profile into stellar and dark matter (DM) components. Unless one accounts for the DM, a standard Virial analysis of the stars dramatically over-estimates the stellar mass of the galaxy. We find the measured J-band stellar mass-to-light ratio, 1.37+/-0.10 Msun/Lsun, is in good agreement with simple stellar population model calculations for this object.

The globular cluster omega Centauri is one of the largest and most massive members of the galactic system. However, its classification as a globular cluster has been challenged making it a candidate for being the stripped core of an accreted dwarf galaxy; this together with the fact that it has one of the largest velocity dispersions for star clusters in our galaxy makes it an interesting candidate for harboring an intermediate mass black hole. We measure the surface brightness profile from integrated light on an HST}/ACS image of the center, and find a central power-law cusp of logarithmic slope -0.08. We also analyze Gemini GMOS-IFU kinematic data for a 5x5 arcsec field centered on the nucleus of the cluster, as well as for a field 14 arcsecaway. We detect a clear rise in the velocity dispersion from 18.6 km/s at 14 arcsec to 23 km/s in the center. A rise in the velocity dispersion could be due to a central black hole, a central concentration of stellar remnants, or a central orbital structure that is radially biased. We discuss each of these possibilities. An isotropic, spherical dynamical model implies a black hole mass of 4.0^{+0.75}_{-1.0} times 10^4 M_sun, and excludes the no black hole case at greater than 99% significance. We have also run flattened, orbit-based models and find similar results. While our preferred model is the existence of a central black hole, detailed numerical simulations are required to confidently rule out the other possibilities.

We present a technique to extract ultra-deep diffuse-light spectra from the standard multi-object spectroscopic observations used to investigate extragalactic globular cluster (GC) systems. This technique allows a clean extraction of the spectrum of the host galaxy diffuse light from the same slitlets as the GC targets. We show the utility of the method for investigating the kinematics and stellar populations of galaxies at radii much greater than usually probed in longslit studies, at no additional expense in terms of telescope time. To demonstrate this technique we present Gemini/GMOS spectroscopy of 29 GCs associated with the elliptical galaxy NGC 3923. We compare the measured stellar population parameters of the GC system with those of the spheroid of NGC 3923 at the same projected radii, and find the GCs to have old ages (> 10 Gyr), [alpha/Fe]~0.3 and a range of metallicities running from [Z/H] = -1.8 to +0.35. The diffuse light of the galaxy is found to have ages, metallicities and [alpha/Fe] abundance ratios indistinguishable from those of the red GCs.

The observational evidence for central black holes in globular clusters has been argued extensively, and their existence has important consequences for both the formation and evolution of the cluster. Most of the evidence comes from dynamical arguments, but the interpretation is difficult, given the short relaxation times and old ages of the clusters. One of the most robust signatures for the existence of a black hole is radio and/or X-ray emission. We observed three globular clusters, NGC6093 (M80), NGC6266 (M62), and NGC7078 (M15), with the VLA in the A and C configuration with a 3-sigma noise of 36, 36 and 25 microJy, respectively. We find no statistically-significant evidence for radio emission from the central region for any of the three clusters. NGC6266 shows a 2-sigma detection. It is difficult to infer a mass from these upper limits due to uncertainty about the central gas density, accretion rate, and accretion model.

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.

Dynamical models for 17 Coma early-type galaxies are presented. The galaxy sample consists of flattened, rotating as well as non-rotating early-types including cD and S0 galaxies with luminosities between M=-18.79 and M=-22.56. Kinematical long-slit observations cover at least the major and minor axis and extend to 1-4 effective radii. Axisymmetric Schwarzschild models are used to derive stellar mass-to-light ratios and dark halo parameters. In every galaxy models with a dark matter halo match the data better than models without. The statistical significance is over 95 percent for 8 galaxies, around 90 percent for 5 galaxies and for four galaxies it is not significant. For the highly significant cases systematic deviations between observed and modelled kinematics are clearly seen; for the remaining galaxies differences are more statistical in nature. Best-fit models contain 10-50 percent dark matter inside the half-light radius. The central dark matter density is at least one order of magnitude lower than the luminous mass density. The central phase-space density of dark matter is often orders of magnitude lower than in the luminous component, especially when the halo core radius is large. The orbital system of the stars along the major-axis is slightly dominated by radial motions. Some galaxies show tangential anisotropy along the minor-axis, which is correlated with the minor-axis Gauss-Hermite coefficient H4. Changing the balance between data-fit and regularisation constraints does not change the reconstructed mass structure significantly. Model anisotropies tend to strengthen if the weight on regularisation is reduced, but the general property of a galaxy to be radially or tangentially anisotropic, respectively, does not change. (abridged)

We use Hubble Space Telescope archival images to measure central surface brightness profiles of globular clusters around satellite galaxies of the Milky Way. We report results for 21 clusters around the LMC, 5 around the SMC, and 4 around the Fornax dwarf galaxy. The profiles are obtained using a recently developed technique based on measuring integrated light, which is tested on an extensive simulated dataset. Our results show that for 70% of the sample, the central photometric points of our profiles are brighter than previous measurements using star counts with deviations as large as 2 mag/arcsec^2. About 40% of the objects have central profiles deviating from a flat central core, with central logarithmic slopes continuously distributed between -0.2 and -1.2. These results are compared with those found for a sample of Galactic clusters using the same method. We confirm the known correlation in which younger clusters tend to have smaller core radii, and we find that they also have brighter central surface brightness values. This seems to indicate that globular clusters might be born relatively concentrated, and that a profile with extended flat cores might not be the ideal choice for initial profiles in theoretical models.

The near-infrared integral field spectrograph SINFONI at the ESO VLT opens a new window for the study of central supermassive black holes. With a near-IR spatial resolution similar to HST optical and the ability to penetrate dust it provides the possibility to explore the low-mass end of the M-sigma relation (sigma<120km/s) where so far very few black hole masses were measured with stellar dynamics. With SINFONI we observed the central region of the low-luminosity elliptical galaxy NGC4486a at a spatial resolution of ~0.1arcsec in the K band. The stellar kinematics was measured with a maximum penalised likelihood method considering the region around the CO absorption band heads. We determined a black hole mass of M_BH=1.25^{+0.75}_{-0.79} x 10^7 M_sun (90% C.L.) using the Schwarzschild orbit superposition method including the full 2-dimensional spatial information. This mass agrees with the predictions of the M-sigma relation, strengthening its validity at the lower sigma end.

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 investigate the relationship between black hole mass and host galaxy velocity dispersion for QSOs in Data Release 3 of the Sloan Digital Sky Survey. We derive black hole mass from the broad Hbeta line width and continuum luminosity, and the bulge stellar velocity dispersion from the [OIII] narrow line width. At higher redshifts, we use MgII and [OII] in place of Hbeta and [OIII]. For redshifts z < 0.5, our results agree with the black hole mass - bulge velocity dispersion relationship for nearby galaxies. For 0.5 < z < 1.2, this relationship appears to show evolution with redshift in the sense that the bulges are too small for their black holes. However, we find that part of this apparent trend can be attributed to observational biases, including a Malmquist bias involving the QSO luminosity. Accounting for these biases, we find ~0.2 dex evolution in the black hole mass-bulge velocity dispersion relationship between now and redshift z ~ 1.

This article was co-authored by all invited speakers at the Joint Discussion on `Neutron Stars and Black Holes in Star Clusters,' which took place during the IAU General Assembly in Prague, Czech Republic, on August 17 and 18, 2006. Each section presents a short summary of recent developments in a key area of research, incorporating the main ideas expressed during the corresponding panel discussion at the meeting. Our meeting, which had close to 300 registered participants, was broadly aimed at the large community of astronomers around the world working on the formation and evolution of compact objects and interacting binary systems in dense star clusters, such as globular clusters and galactic nuclei. Great advances have occurred in this field during the past few years, including the introduction of fundamentally new theoretical paradigms for the formation and evolution of compact objects in binaries as well as countless new discoveries by astronomers that have challenged many accepted models. Some of the highlights include: a nearly complete census of all the millisecond pulsars in 47 Tuc; first detections of many new radio pulsars in other clusters, particularly Terzan 5; detailed studies of X-ray binary populations and their luminosity functions in many galaxies and extragalactic globular clusters; increasing evidence for intermediate-mass black holes in clusters and greatly improved theoretical understanding of their possible formation processes.

Accretion disks around supermassive black holes are widely believed to be the dominant source of the optical-ultraviolet continuum in many classes of active galactic nuclei (AGN). We study here the relationship between the continuum colors of AGN and the characteristic accretion disk temperature (T_max). Based on NLTE models of accrection disks in AGN computed as described by Hubeny et al. (2000), we find that continuum intensity ratios for several pairs of wavelengths between 1350 and 5100 A should show a trend of bluer colors for higher T_max, notwithstanding random disk inclinations. We compare this theoretical expectation with observed colors of QSOs in the Sloan Digital Sky Survey,deriving black hole mass and thence T_max from the width of the Mg II broad emission line. The observed colors generally do not show the expected trend and in some cases show a reverse trend of redder colors with increasing T_max. The cause of this discrepancy does not appear to be dust reddening or galaxy contamination but may relate to the accretion rate, as the offset objects are accreting above ~30 % of the Eddington limit. The derived disk temperature depends primarily on line width, with little or no dependence on luminosity.

Spectra of galaxies contain an enormous amount of information about the relative mixture of ages and metallicities of constituent stars. We present a comprehensive study designed to extract the maximum information from spectra of data quality typical in large galaxy surveys. We test fitting techniques using the Bruzual-Charlot 2003 high resolution simple stellar population synthesis models to simultaneously estimate the ages and metallicities of 101 globular clusters in M31 and the Magellanic Clouds. The clusters cover a wide range of ages and metallicities, 4 Myr to 20 Gyr and -1.6 < [Fe/H] < +0.3, estimated by other methods in the literature. We compare results from model fits to both the spectra and photometry and find that fits to continuum-normalized (CN) spectra over the entire range available, typically 350-1000 nm for this sample, provides the best results. For clusters older than 1 Gyr we agree with literature ages to 0.16 dex (35%) and [Fe/H] to 0.12 dex. For younger clusters we agree with literature ages to 0.3 dex (63%), but cannot constrain the metallicity. It is particularly important to use the entire CN spectrum to avoid problems with model continua for young objects and to break age-metallicity degeneracies of broadband photometry. Our required S/N = 15-30 for 20% age uncertainties and 30-55 for 10% uncertainties over spectral resolutions of 5-25 Angstroms. This technique should work well for the age-metallicity parameter space expected for early-type galaxies at z~1, although individual galaxy spectral S/N may require the coaddition of a few like objects.

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.

The phases as well as the amplitudes of baryonic acoustic oscillations can be extracted out of a galaxy power spectrum and used as a standard ruler in a cosmological test. A non-oscillating phenomenological fitting function is used to extract the oscillatory part of the galaxy power spectrum and to disentangle phase information from amplitude information. The method (FITEX) is tested with simulated data of the Hubble Volume Simulation to include redshift space effects, non-linear structure growth and biasing. A cosmological test is introduced, which compares the extracted oscillations against a theoretical model template to derive constraints on the $w$ parameter. The phenomenological fitting function is found to model the various distortions of the galaxy power spectrum to the sub-percent level. The various distortions only boost the amplitude of the oscillations. The theoretical template is accurate enough to test for small deviations in the phases of the oscillations, resulting from different $w$ values. A cosmological test, using the baryonic acoustic oscillations as a standard ruler, is able to constrain $w$ in a robust way.

From observations with the GMOS multi-slit spectrograph on the Gemini North telescope, we have obtained spectra for 39 globular cluster candidates in the Virgo giant elliptical galaxy NGC 4649 (M60), of which 38 are confirmed globular clusters. The clusters extend out to a radius of 260" (3.5 effective radii). We find no rotation of the globular cluster system, with an upper limit of v/sigma < 0.6 at a confidence level of 95%. The globular cluster velocity dispersion is constant with radius, within the uncertainties. We fit isotropic models to the globular cluster and stellar kinematics; these models yield a M/L_V around 16 at 200" radius (16 kpc), an increase of a factor of two from the central M/L. We also use the mass profile as derived from X-rays to determine the orbital structure. Using axisymmetric orbit-based models and the X-ray mass profile, we find the orbital distribution is close to isotropic within 100", and becomes tangentially biased beyond. Furthermore, when using the X-ray profile, we find a better fit to the kinematics compared to using a constant M/L model. Thus, both isotropic and axisymmetric orbit-based models give support for the presence of a dark matter halo in NGC 4649.

We have used HST imaging of the central regions (R<100 arcsec, about 5 core radii) of the globular cluster 47 Tucanae to derive proper motions and U- and V-band magnitudes for 14,366 cluster members. We also present a catalogue of astrometry and F475W photometry for nearly 130,000 stars in a rather larger central area. These data are made available in their entirety, in the form of downloadable electronic tables. We use them first to obtain a new estimate for the position of the cluster center and to define the stellar density profile into essentially zero radius. We then search in particular for any very fast-moving stars, such as might be expected to result from very close stellar encounters. Likely fewer than 0.1% (and no more than about 0.3%) of stars have total speeds above the nominal central escape velocity in 47 Tuc, and at lower speeds the velocity distribution is described very well by a regular King model. Considerations of only the proper-motion velocity dispersion then lead to a number of results: (1) Blue stragglers in the core of 47 Tuc have a velocity dispersion lower than that of the cluster giants by a factor of sqrt{2}. (2) The velocity distribution in the cluster center is essentially isotropic, as expected. (3) Using a sample of radial velocities for stars in the core, we estimate the distance to 47 Tuc: D = 4.0 +/- 0.35 kpc. And (4) we infer a 1-sigma upper limit of M<1000-1500 solar masses for any central, intermediate-mass black hole. We can neither confirm nor refute the hypothesis that 47 Tuc might lie on an extension of the M-sigma relation observed for supermassive black holes in galaxy bulges. [Abridged]

We combine Hubble Space Telescope spectroscopy and ground-based integral-field data from the SAURON and OASIS instruments to study the central black hole in the nearby elliptical galaxy NGC 3379. From these data, we obtain kinematics of both the stars and the nuclear gaseous component. Axisymmetric three-integral models of the stellar kinematics find a black hole of mass 1.4 (+2.6 / -1.0) x 10^8 M_sun (3 sigma errors). These models also probe the velocity distribution in the immediate vicinity of the black hole and reveal a nearly isotropic velocity distribution throughout the galaxy and down to the black hole sphere of influence R_BH. The morphology of the nuclear gas disc suggests that it is not in the equatorial plane; however the core of NGC 3379 is nearly spherical. Inclined thin-disc models of the gas find a nominal black hole of mass 2.0 (+/- 0.1) x 10^8 M_sun (3 sigma errors), but the model is a poor fit to the kinematics. The data are better fit by introducing a twist in the gas kinematics (with the black hole mass assumed to be 2.0 x 10^8 M_sun), although the constraints on the nature and shape of this perturbation are insufficient for more detailed modelling. Given the apparent regularity of the gas disc's appearance, the presence of such strong non-circular motion indicates that caution must be used when measuring black hole masses with gas dynamical methods alone.

We identify active galactic nuclei (AGN) in the Groth-Westphal Survey Strip (GSS) using the independent and complementary selection techniques of optical spectroscopy and photometric variability. We discuss the X-ray properties of these AGN using Chandra/XMM data for this region. From a sample of 576 galaxies with high quality spectra we identify 31 galaxies with AGN signatures. Seven of these have broad emission lines (Type 1 AGNs). We also identify 26 galaxies displaying nuclear variability in HST WFPC2 images of the GSS separated by ~7 years. The primary overlap of the two selected AGN samples is the set of broad-line AGNs, of which 80% appear as variable. Only a few narrow-line AGNs approach the variability threshold. The broad-line AGNs have an average redshift of z~1.1 while the other spectroscopic AGNs have redshifts closer to the mean of the general galaxy population (z~0.7). Eighty percent of the identified broad-line AGNs are detected in X-rays and these are among the most luminous X-ray sources in the GSS. Only one narrow-line AGN is X-ray detected. Of the variable nuclei galaxies within the X-ray survey, 27% are X-ray detected. We find that 1.9+/-0.6% of GSS galaxies to V=24 are broad-line AGNs, 1.4+/-0.5% are narrow-line AGNs, and 4.5+/-1.4% contain variable nuclei. The fraction of spectroscopically identified BLAGNs and NLAGNs at z~1 reveals a marginally significant increase of 1.3+/-0.9% when compared to the local population.

Hubble Space Telescope allows us to study the central surface brightness profiles for globular clusters at unprecedented detail. We have mined the HST archives to obtain 38 WFPC2 images of galactic globular clusters with adequate exposure times and filters, which we use to measure their central structure. We outline a reliable method to obtain surface brightness profiles from integrated light that we test on an extensive set of simulated images. Most clusters have central surface brightness about 0.5 mag brighter than previous measurements made from ground-based data, with the largest differences around 2 magnitudes. Including the uncertainties in the slope estimates, the surface brightness slope distribution is consistent with half of the sample having flat cores and the remaining half showing a gradual decline from 0 to -0.8 (dlog(Sigma)/dlogr). We deproject the surface brightness profiles in a non-parametric way to obtain luminosity density profiles. The distribution of luminosity density logarithmic slopes show similar features with half of the sample between -0.4 and -1.8. These results are in contrast to our theoretical bias that the central regions of globular clusters are either isothermal (i.e. flat central profiles) or very steep (i.e. luminosity density slope ~-1.6) for core-collapse clusters. With only 50% of our sample having central profiles consistent with isothermal cores, King models appear to poorly represent most globular clusters in their cores.

NGC 4649 (M60) is one of a handful of giant Virgo ellipticals. We have obtained Gemini/GMOS spectra for 38 GCs associated with this galaxy. Applying the multi-index chi^2 minimisation technique of Proctor & Sansom (2002) with the single stellar population models of Thomas, Maraston & Korn (2004) we derive ages, metallicities and alpha-element abundance ratios. We find several young (2--3 Gyr old) super-solar metallicity GCs, while the majority are old (>10 Gyrs), spanning a range of metallicities from solar to [Z/H]=-2. At least two of these young GCs are at large projected radii of 17-20 kpc. The galaxy itself shows no obvious signs of a recent starburst, interaction or merger. A trend of decreasing alpha-element ratio with increasing metallicity is found.

We construct orbit-based axisymmetric dynamical models for the globular cluster M15 which fit groundbased line-of-sight velocities and Hubble Space Telescope line-of-sight velocities and proper motions. This allows us to constrain the variation of the mass-to-light ratio M/L as a function of radius in the cluster, and to measure the distance and inclination of the cluster. We obtain a best-fitting inclination of 60+/-15 degrees, a dynamical distance of 10.3+/-0.4 kpc and an M/L profile with a central peak. The inferred mass in the central 0.05 parsec is 3400 Msun, implying a central density of at least 7.4x10^6 Msun pc^-3. We cannot distinguish the nature of the central mass concentration. It could be an IMBH or it could be large number of compact objects, or it could be a combination. The central 4 arcsec of M15 appears to contain a rapidly spinning core, and we speculate on its origin.

We have obtained Gemini/GMOS spectra for 22 GCs associated with NGC 3379. We derive ages, metallicities and alpha-element abundance ratios from simple stellar population models using the multi-index chi^2 minimisation method of Proctor & Sansom (2002). All of these GCs are found to be consistent with old ages, i.e. >10 Gyr, with a wide range of metallicities. A trend of decreasing alpha-element abundance ratio with increasing metallicity is indicated. The projected velocity dispersion of the GC system is consistent with being constant with radius. Non-parametric, isotropic models require a significant increase in the mass-to-light ratio at large radii. This result is in contrast to that of Romanowsky et al. (2003) who find a decrease in the velocity dispersion profile as determined from planetary nebulae. Our constant dispersion requires a normal sized dark halo, although without anisotropic models we cannot rigorously determine the dark halo mass. A two-sided chi^2 test over all radii, gives a 2 sigma difference between the mass profile derived from our GCs compared to the PN-derived mass model of Romanowsky et al. (2003). However, if we restrict our analysis to radii beyond one effective radius and test if the GC velocity dispersion is consistently higher, we determine a >3 sigma difference between the mass models, and hence favor the conclusion that NGC 3379 does indeed have dark matter at large radii in its halo. (abridged)

We present dynamical models for the massive globular cluster G1. The goal is to measure or place a significant upper limit on the mass of any central black hole. Whether or not globular clusters contain central massive black holes has important consequences for a variety of studies. We use new kinematic data obtained with Keck and new photometry from the Hubble Space Telescope. The Keck spectra allow us to obtain kinematics out to large radii that are required to pin down the mass-to-light ratio of the dynamical model and the orbital structure. The Hubble Space Telescope observations give us a factor of two better spatial resolution for the surface brightness profile. By fitting non-parametric, spherical, isotropic models we find a best-fit black hole mass of 1.7(+-0.3)e4 Msun. Fully general axisymmetric orbit-based models give similar results, with a black hole mass of 1.8(+-0.5)e4 Msun. The no-black hole model has Delta_chi^2=5 (marginalized over mass-to-light ratio), implying less than 3% significance. We have taken into account any change in the mass-to-light ratio in the center due to stellar remnants. These results are consistent with our previous estimate in Gebhardt, Rich & Ho (2002), and inconsistent with the analysis of Baumgardt et al. (2003) who claim that G1 does not show evidence for a black hole. These new results make G1 the best example of a cluster that contains an intermediate-mass black hole.

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 present the analysis of Lick absorption-line indices from three separate long-slit spectroscopic observations of the nearby isolated elliptical galaxy NGC 821. The three data sets present a consistent picture of the stellar population within one effective radius, in which strong gradients are evident in both luminosity-weighted age and metallicity. The central population exhibits a young age of ~4 Gyr and a metallicity ~3 times solar. At one effective radius the age has risen to ~12 Gyr and the metallicity fallen to less than 1/3 solar. The low metallicity population around one effective radius appears to have an exclusively red horizontal branch, with no significant contribution from the blue horizontal branch evident in some globular clusters of the same age and metallicity. Despite the strong central age gradient, we demonstrate that only a small fraction (<10%) of the galaxy's stellar mass can have been created in recent star formation events. We consider possible star formation histories for NGC 821 and find that the most likely cause of the young central population was a minor merger or tidal interaction that caused NGC 821 to consume its own gas in a centrally concentrated burst of star formation 1 to 4 Gyr ago.

The DEEP2 and COMBO-17 surveys are used to study the evolution of the luminosity function of red and blue galaxies to $z \sim 1$. Schechter function fits show that, since $z = 1$, $M^*_B$ dims by $\sim$ 1.3 mag per unit redshift for both color classes, $\phi^*$ of blue galaxies shows little change, while $\phi^*$ for red galaxies has formally nearly quadrupled. At face value, the number density of blue galaxies has remained roughly constant since $ z = 1$, whereas that of red galaxies has been rising. Luminosity densities support both conclusions, but we note that most red-galaxy evolution occurs between our data and local surveys and in our highest redshift bin, where the data are weakest. We discuss the implications of having most red galaxies emerge after $z = 1$ from precursors among the blue population, taking into account the properties of local and distant E/S0s. We suggest a ``mixed'' scenario in which some blue galaxies have their star-formation quenched in gas-rich mergers, migrate to the red sequence with a variety of masses, and merge further on the red sequence in one or more purely stellar mergers. E/S0s of a given mass today will have formed via different routes, in a manner that may help to explain the fundamental plane and other local scaling laws.

This is the second in a series of papers dedicated to unveil the mass structure and orbital content of a sample of flattened early-type galaxies in the Coma cluster. The ability of our orbit libraries to reconstruct internal stellar motions and the mass composition of a typical elliptical in the sample is investigated by means of Monte-Carlo simulations of isotropic rotator models. The simulations allow a determination of the optimal amount of regularization needed in the orbit superpositions. It is shown that under realistic observational conditions and with the appropriate regularization internal velocity moments can be reconstructed to an accuracy of about 15 per cent; the same accuracy can be achieved for the circular velocity and dark matter fraction. In contrast, the flattening of the halo remains unconstrained. Regularized orbit superpositions are applied to a first galaxy in our sample, NGC 4807, for which stellar kinematical observations extend to 3 Reff. The galaxy seems dark matter dominated outside 2 Reff. Logarithmic dark matter potentials are consistent with the data, as well as NFW-profiles, mimicking logarithmic potentials over the observationally sampled radial range. In both cases, the derived stellar mass-to-light ratio agrees well with independently obtained mass-to-light ratios from stellar population analysis. Kinematically, NGC 4807 is characterized by mild radial anisotropy outside r>0.5 Reff, becoming isotropic towards the center. Our orbit models hint at either a distinct stellar component or weak triaxiality in the outer parts of the galaxy.

The Deep Extragalactic Exploratory Probe (DEEP) is a multi-phase research program dedicated to the study of the formation and evolution of galaxies and of large scale structure in the distant Universe. This paper describes the first five-year phase, denoted DEEP1. A series of ten DEEP1 papers will discuss a range of scientific topics (e.g., the study of photometric and spectral properties of a general distant galaxy survey, the evolution observed in galaxy populations of varied morphologies). The observational basis for these studies is the Groth Survey Strip field, a 127 square arcminute region which has been observed with the Hubble Space Telescope in both broad I-band and V-band optical filters and with the Low Resolution Imaging Spectrograph on the Keck Telescopes. Catalogs of photometric and structural parameters have been constructed for 11,547 galaxies and stars at magnitudes brighter than 29, and spectroscopy has been conducted for a magnitude-color weighted subsample of 818 objects. We evaluate three independent techniques for constructing an imaging catalog for the field from the HST data, and discuss the depth and sampling of the resultant catalogs. The selection of the spectroscopic subsample is discussed, and we describe the multifaceted approach taken to prioritizing objects of interest for a variety of scientific subprograms. A series of Monte Carlo simulations then demonstrates that the spectroscopic subsample can be adequately modeled as a simple function of magnitude and color cuts in the imaging catalog.

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.

We present a sample of over 50 luminous field bulges (including ellipticals) found in the Groth Strip Survey (GSS), with 0.73< z < 1.04 and with bulge magnitudes I <= 23. The exponential disk light is removed via decomposition of HST images using GIM2D. We find that 85% of these bulges are nearly as red as local E/S0's and have a shallow slope and a small color dispersion in the color-luminosity relation, suggesting roughly coeval formation. The surface brightnesses of these bulges are about 1 mag higher than local bulges. These results are explained adopting a "drizzling" scenario where a metal-rich early formation is later polluted by small amounts of additional star formation. Almost all disks have the same or bluer colors than their accompanying bulges, regardless of the bulge-disk ratio and bulge luminosity, as expected from semi-analytic hierarchical galaxy formation models. We present evidence that the few blue bulge candidates are not likely to be genuine blue ellipticals or bulges. Our deeper, more extensive, and less disk-contaminated observations challenge prior claims that 30% to 50% of field bulges or ellipticals are in a blue, star-forming phase at z < 1. We conclude that field bulges and ellipticals at z ~ 1, like luminous early- type cluster galaxies at the same redshift, are already dominated by metal-rich, old stellar populations that have been fading from a formation epoch earlier than z ~ 1.5. (abridged)

The Deep Extragalactic Evolutionary Probe (DEEP) is a series of spectroscopic surveys of faint galaxies, targeted at the properties and clustering of galaxies at redshifts z ~ 1. We present the redshift catalog of the DEEP 1 GSS pilot phase of this project, a Keck/LRIS survey in the HST/WFPC2 Groth Survey Strip. The redshift catalog and data, including reduced spectra, are publicly available through a Web-accessible database. The catalog contains 658 secure galaxy redshifts with a median z=0.65, and shows large-scale structure walls to z = 1. We find a bimodal distribution in the galaxy color-magnitude diagram which persists to z = 1. A similar color division has been seen locally by the SDSS and to z ~ 1 by COMBO-17. For red galaxies, we find a reddening of only 0.11 mag from z ~ 0.8 to now, about half the color evolution measured by COMBO-17. We measure structural properties of the galaxies from the HST imaging, and find that the color division corresponds generally to a structural division. Most red galaxies, ~ 75%, are centrally concentrated, with a red bulge or spheroid, while blue galaxies usually have exponential profiles. However, there are two subclasses of red galaxies that are not bulge-dominated: edge-on disks and a second category which we term diffuse red galaxies (DIFRGs). The distant edge-on disks are similar in appearance and frequency to those at low redshift, but analogs of DIFRGs are rare among local red galaxies. DIFRGs have significant emission lines, indicating that they are reddened mainly by dust rather than age. The DIFRGs in our sample are all at z>0.64, suggesting that DIFRGs are more prevalent at high redshifts; they may be related to the dusty or irregular extremely red objects (EROs) beyond z>1.2 that have been found in deep K-selected surveys. (abridged)

We present Sloan g and i imaging from the GMOS instrument on the Gemini North telescope for the globular cluster (GC) system around the Virgo galaxy NGC 4649 (M60). Our three pointings, taken in good seeing conditions, cover an area of about 90 sq. arcmins. We detect 2,151 unresolved sources. Applying colour and magnitude selection criteria to this source list gives 995 candidate GCs that is greater than 90% complete to a magnitude of i = 23.6, with little contamination from background galaxies. We find fewer than half a dozen potential Ultra Compact Dwarf galaxies around NGC 4649. Foreground extinction from the nearby spiral NGC 4647 is limited to be A_V < 0.1. We confirm the bimodality in the GC colour distribution found by earlier work using HST/WFPC2 imaging. As is commonly seen in other galaxies, the red GCs are concentrated towards the centre of the galaxy, having a steeper number density profile than the blue GC subpopulation. The varying ratio of red-to-blue GCs with radius can largely explain the overall GC system colour gradient. The underlying galaxy starlight has a similar density profile slope and colour to the red GCs. This suggests a direct connection between the galaxy field stars and the red GC subpopulation. We estimate a total GC population of 3700 +/- 900, with the uncertainty dominated by the extrapolation to larger radii than observed. This total number corresponds to a specific frequency S_N = 4.1 +/- 1.0. Future work will present properties derived from GMOS spectra of the NGC 4649 GCs.

Black hole masses in active galactic nuclei (AGN) are difficult to measure using conventional dynamical methods, but can be determined using the technique of reverberation mapping. However, it is important to verify that the results of these different methods are equivalent. This can be done indirectly, using scaling relations between the black hole and the host galaxy spheroid. For this purpose, we have obtained new measurements of the bulge stellar velocity dispersion, sigma, in Seyfert 1 galaxies. These are used in conjunction with the M_bh -- sigma relation to validate nuclear black hole masses, M_bh, in active galaxies determined through reverberation mapping. We find that Seyfert galaxies follow the same M_bh -- sigma relation as non-active galaxies, indicating that reverberation mapping measurements of M_bh are consistent with those obtained using other methods. We also reconsider the relationship between bulge absolute magnitude, M_bulge, and black hole mass. We find that Seyfert galaxies are offset from non-active galaxies, but that the deviation can be entirely understood as a difference in bulge luminosity, not black hole mass; Seyfert hosts are brighter than normal galaxies for a given value of their velocity dispersion, perhaps as a result of younger stellar populations.

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.

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.

Radial velocities of 468 globular clusters around NGC 1399, the central galaxy in the Fornax cluster, have been obtained with FORS2 and the Mask Exchange Unit (MXU) at the ESO Very Large Telescope. This is the largest sample of globular cluster velocities around any galaxy obtained so far. The mean velocity uncertainty is 50 km/sec. This data sample is accurate and large enough to be used in studies of the mass distribution of NGC 1399 and the properties of its globular cluster system. Here we describe the observations, the reduction procedure, and discuss the uncertainties of the resulting velocities. The complete sample of cluster velocities which is used in a dynamical study of NGC 1399 is tabulated. A subsample is compared with previously published values.

We study the kinematics and dynamics of the globular cluster system of NGC 1399, the central galaxy of the Fornax cluster. The observational data consists of medium resolution spectra, obtained at the Very Large Telescope. Our sample comprises 468 radial velocities in the magnitude range 20 < m_R < 23. This is the largest sample of globular cluster velocities around any galaxy obtained so far. The radial range is 2 arcmin < r < 9 arcmin, corresponding to 11 kpc to 50 kpc of galactocentric distance. There is the possibility that unbound clusters and/or objects in the foreground contaminate the NGC 1399 cluster sample. Under strong error selection, practically no objects are found with velocities lower than 800 km/s or higher than 2000 km/s. Since the extreme velocities influence the velocity dispersion considerably, uncertainty regarding the exact value of the dispersion remains. Within the above velocity limits, we derive a projected velocity dispersion for the total sample of 274+-9 km/s which within the uncertainties remains constant over the entire radial range. Without any velocity restriction, it increases to 325 km/s. Blue and red clusters show different dispersions corresponding to their different surface density profiles. Spherical models point to a circular velocity of 415+-30$ km/s, assuming isotropy for the red clusters. This value is constant out to 40 kpc. The inferred dark halo potential can be well represented by a logarithmic potential. Also a halo of the NFW type fits well to the observations. Some mass profiles derived from X-ray analyses do not agree with a constant circular velocity within our radial range, irrespective of its exact value.

We present photometry from Hubble Space Telescope (HST)/Wide Field Planetary Camera 2 parallel imagery of a remote M31 field at a projected distance of about 34 kpc from the nucleus near the SW major axis. This field is near the globular cluster G1, and near one of the candidate tidal plumes identified by Ferguson et al. (2002). The F606W (V) and F814W (I) images were obtained in parallel with Space Telescope Imaging Spectrograph spectroscopy of G1 (GO-9099) and total 7.11 hours of integration time -- the deepest HST field in the outer disk of M31 to date, reaching to V ~ 28. The color-magnitude diagram of the field shows a clearly-defined red clump at V = 25.25 and a red giant branch consistent with [Fe/H] ~ -0.7. The lack of a blue horizontal branch contrasts with other M31 halo fields, the Andromeda dwarf spheroidals, and with the nearby globular cluster G1. Comparing the observed luminosity function to the Padova models, we find that at least some of the stellar population must be younger than 6 - 8 Gyr. The outermost detected neutral hydrogen gas disk of M31 lies only 2 kpc in projection from our field. The finding that some giants in the field have radial velocities close to that of the neutral hydrogen gas (Reitzel, Guhathakurta, & Rich 2003) leads us to conclude that our field samples the old, low-surface-brightness disk rather than the true Population II spheroid.

We present recent results from our long-term Gemini/GMOS study of globular clusters in early-type galaxies. To date, we have obtained photometry and spectroscopy for clusters in NGCs 3379, 4649, 524, 7332, and IC 1459. We find a clear bimodality in the NGC 4649 cluster colour distribution, with the fraction of blue/red clusters increasing with galactocentric radius. We derive ages and metallicities for 22 clusters in NGC 3379, finding that most of the clusters appear old (10-15 Gyr); however, there is a group of 4 metal-rich, younger clusters with ages of 2-6 Gyr. The NGC 3379 cluster velocity dispersion decreases with radius, as does the inferred (local) mass-to-light ratio: there is no evidence for a dark matter halo in NGC 3379 based on our globular cluster data.

Fundamental Plane studies provide an excellent means of understanding the evolutionary history of early-type galaxies. Using the Low Resolution Imaging Spectrograph on the Keck telescope, we obtained internal stellar kinematic information for 36 field galaxies in the Groth Strip--21 early-type and 15 disk galaxies. Their redshifts range from 0.3--1.0, with a median redshift 0.8. The slope of the relation shows no difference compared with the local slope. However, there is significant evolution in the zero-point offset; an offset due to evolution in magnitude requires a 2.4 magnitude luminosity brightening at z=1. We see little differences of the offset with bulge fraction, which is a good surrogate for galaxy type. Correcting for the luminosity evolution reduces the orthogonal scatter in the Fundamental Plane to 8%, consistent with the local scatter. This scatter is measured for our sample, and does not include results from other studies which may have different selection effects. The difference in the degree of evolution between our field sample and published cluster galaxies suggests a more recent formation epoch--around z=1.5 for field galaxies compared to z>2.0 for cluster galaxies. The magnitude difference implies that the field early-type galaxies are about 2 Gyr younger than the cluster ellipticals using standard single-burst models. However, the same models imply a significant change in the rest-frame U-B color from then to present, which is not seen in our sample. Continuous low-level star formation, however, would serve to explain the constant colors over this large magnitude change. A consistent model has 7% of the stellar mass created after the initial burst, using an exponentially decaying star formation rate with an e-folding time of 5 Gyr.

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.

We review the current state of dynamical modeling for galaxies in terms of being able to measure both the central black hole mass and stellar orbital structure. Both of these must be known adequately to measure either property. The current set of dynamical models do provide accurate estimates of the black hole mass and the stellar orbital distribution. Generally, these models are able to measure the black hole mass to about 20%--30% accuracy given present observations, and the stellar orbital structure to about 20% accuracy in the radial to tangential dispersions. The stellar orbital structure of the stars near the galaxy center show strong tangential velocity anisotropy for most galaxies studied. Theoretical models that best match this trend are black hole binary/merger models. There is also a strong correlation between black hole mass and the contribution of radial motion at large radii. This correlation may be an important aspect of galaxy evolution.

We have obtained spectra and measured the stellar kinematics in a sample of 25 nearby early-type galaxies (with velocity dispersions from less than 100 km/s to over 300 km/s) using the near-infrared CO absorption bandhead at 2.29 microns. Our median uncertainty for the dispersions is ~10%. We examine the effects of dust on existing optical kinematic measurements. We find that the near-infrared velocity dispersions are in general smaller than optical velocity dispersions, with differences as large as 30%. The median difference is 11% smaller, and the effect is of greater magnitude for higher dispersion galaxies. The lenticular galaxies (18 out of 25) appear to be causing the shift to lower dispersions while the classical ellipticals (7 out of 25) are consistent between the two wavelength regimes. If uniformly distributed dust causes these differences, we would expect to find a correlation between the relative amount of dust in a galaxy and the fractional change in dispersion, but we do not find such a correlation. We do see correlations both between velocity dispersion and CO bandhead equivalent width, and velocity dispersion and Mg2 index. The differences in dispersion are not well explained by current models of dust absorption. The lack of correlation between the relative amount of dust and shift in dispersion possibly suggets that dust does not have a similar distribution from galaxy to galaxy. The CO equivalent widths of these galaxies are quite high (>10 angstroms for almost all), requiring the light at these wavelengths to be dominated by very cool stars.

We analyze HST/STIS spectra (see Paper I) of the central region of the dense globular cluster M15. We infer the velocities of 64 individual stars, two-thirds of which have their velocity measured for the first time. This triples the number of stars with measured velocities in the central 1 arcsec of M15 and doubles the number in the central 2 arcsec. Combined with existing ground-based data we obtain the radial profiles of the projected kinematical quantities. The RMS velocity sigma_RMS rises to 14 km/s in the central few arcsec, somewhat higher than the values of 10-12 km/s inferred previously from ground-based data. To interpret the results we construct dynamical models based on the Jeans equation, which imply that M15 must have a central concentration of non-luminous material. If this is due to a single black hole, then its mass is M_BH = (3.9 +/- 2.2) x 10^3 solar masses. This is consistent with the relation between M_BH and sigma_RMS that has been established for galaxies. Also, the existence of intermediate-mass black holes in globular clusters is consistent with several scenarios for globular cluster evolution proposed in the literature. Therefore, these results may have important implications for our understanding of the evolution of globular clusters, the growth of black holes, the connection between globular cluster and galaxy formation, and the nature of the recently discovered `ultra-luminous' X-ray sources in nearby galaxies. Instead of a single black hole, M15 could have a central concentration of dark remnants (e.g., neutron stars) due to mass segregation. However, the best-fitting Fokker-Planck models that have previously been constructed for M15 do not predict a central mass concentration that is sufficient to explain the observed kinematics.[ABRIDGED]

It has been reported that there is an error in the figure in Dull et al. (1997, D97) that shows the radial M/L profile in Fokker-Planck models of M15. We discuss how this modifies the interpretation of our kinematical data. These imply the existence of a dark and compact mass component near the center of M15, either a single black hole (BH) or a collection of dark remnants that have sunk to the cluster center due to mass segregation. We previously showed that the latter interpretation is in conflict with the D97 M/L profile, which supported the BH interpretation. We repeat our analysis here with the corrected D97 profile. Models without a BH are now found to be statistically acceptable (within 1-sigma), although inclusion of a BH still provides a marginally better fit. It does not necessarily follow that dark remnants are now the preferred interpretation of the data. The D97 models, as well as N-body models by Baumgardt et al. (2002), assume that all neutron stars are retained during cluster evolution. This conflicts with predictions of the neutron-star retention rate (typically below 10 per cent) based on pulsar kick velocities. The presence of a single BH therefore continues to be a viable interpretation of the data. The best fit BH mass with the corrected D97 M/L profile is M_BH = 1.7^{+2.7}_{-1.7} x 10^3 solar masses, and with a constant M/L it is M_BH = 3.2^{+2.2}_{-2.2} x 10^3 solar masses. A model that includes both neutron star escape and mass segregation would probably yield a value between these numbers. This agrees with the correlation between velocity dispersion and BH mass inferred for galaxies. However, with the presently available models and data it is neither uniquely implied nor ruled out that M15 has an intermediate-mass BH.

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 use QSO emission-line widths to examine the black hole mass - sigma relationship as a function of redshift and to extend the relationship to larger masses. Supermassive black holes in galactic nuclei are closely related to the bulge of the host galaxy. The mass of the black hole, M_BH, increases with the bulge luminosity and with the velocity dispersion of the bulge stars, sigma_*. An important clue to the origin of this correlation would be an observational determination of the evolution, if any, in the M_BH - sigma_* relationship as a function of cosmic time. The high luminosity of QSOs affords the potential for studies at large redshifts. We derive black hole masses from the continuum luminosity and the width of the broad Hbeta line and sigma_* from the width of the narrow [O III] lines. We find that radio quiet QSOs conform to the established M_BH - sigma_* relationship up to values M_BH of 10 billion solar masses, with no discernible change in the relationship out to redshifts of $z \approx 3$. These results are consistent with the idea that the growth of supermassive black holes and massive bulges occurred simultaneously.

In this series of two papers we describe a project with STIS on HST to measure the line-of-sight velocities of stars in the central few arcsec of the dense globular cluster M15. This first paper focuses on the observations and reduction of the data. We `scanned' the central region of M15 spectroscopically by consecutively placing an 0.1 arcsec slit at 18 adjacent positions. Careful and complicated data reduction and analysis were required. The spectral pixel size exceeds the velocity dispersion of M15, which puts the project at the limit of what is feasible with STIS. We applied corrections for the following effects: drifts in the STIS wavelength scale during an HST orbit; the orbital velocity component of HST along the line-of-sight to the cluster, and its variations during the HST orbit; and the apparent wavelength shift that is perceived for a star that is not centered in the slit. To make these corrections we present a detailed study of the HST/STIS point-spread function and line-spread function. We also present a catalog with positions and U, B and V magnitudes of 31,983 stars in M15, created from existing HST/WFPC2 data. We extract 19,200 one-dimensional spectra from the STIS data, each for a different position in M15, and each with a velocity scale accurate to better than 2.5 km/s. We develop an algorithm that co-adds the spectra for individual apertures and use it to extract spectra of individual stars with minimum blending and maximum S/N. In Paper II we use these spectra to extract the line-of-sight velocities of individual stars. The velocities constrain the central structure, dynamics and mass distribution of the cluster, and provide evidence for the existence of an intermediate-mass black hole in the center of the cluster.[ABRIDGED]

We present the detection of a 2.0(+1.4,-0.8)x10^4 solar mass black hole (BH) in the stellar cluster G1 (Mayall II), based on data taken with the Space Telescope Imaging Spectrograph onboard the Hubble Space Telescope. G1 is one of the most massive stellar clusters in M31. The central velocity dispersion (25 kms) and the measured BH mass of G1 places it on a linear extrapolation of the correlation between BH mass and bulge velocity dispersion established for nearby galaxies. The detection of a BH in this low-mass stellar system suggests that (1) the most likely candidates for seed massive BHs come from stellar clusters, (2) there is a direct link between massive stellar clusters and normal galaxies, and (3) the formation process of both bulges and massive clusters is similar due to their concordance in the M_BH/sigma relation. Globular clusters in our Galaxy should be searched for central BHs.

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.

We performed spectroscopy of globular clusters associated with NGC 1399 and measured radial velocities of more than 400 clusters, the largest sample ever obtained for dynamical studies. In this progress report, we present the sample and the first preliminary results. Red and blue clusters have slightly different velocity dispersions in accordance with their different density profiles. Their velocity dispersions remain constant with radial distance, which differs from earlier work.

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.

Many ellipticals contain nuclear disks of dust and gas. Some ellipticals contain nuclear disks of stars that are distinct from the rest of the galaxy. We show that the dwarf E2 galaxy NGC 4486A contains both -- it is a "Rosetta stone" object that tells us how nuclear disks evolve. Its properties suggest that, as accreted gas dissipates and settles toward the center, it forms stars and builds a stellar disk. Secular growth may explain not only the most distinct nuclear disks such as those in NGC 4486A but also some of the disky distortions that are commonly seen in elliptical galaxies. That is, density distributions may grow secularly cuspier. This would result in chaotic mixing of stellar orbits in phase space and would tend to make an elliptical galaxy evolve toward a more nearly axisymmetric shape.

Several recent studies find that 10 - 50% of morphologically selected field early-type galaxies at redshifts z <= 1 have blue colors indicative of recent star formation. Such ``blue spheroids''might be massive early-type galaxies with active star formation, perhaps induced by recent merger events. Alternatively, they could be starbursting, low-mass spheroids. To distinguish between these two choices, we have selected 10 ``Blue Spheroid Candidates'' (hereafter, BSCs) from a quantitatively selected E/S0 sample and studied their properties, including kinematics from Keck spectra obtained as part of the DEEP Groth Strip Survey (GSS). Most BSCs (70%) turn out to belong to two broad categories, while the remaining objects are likely to be misclassified objects. Type-1 BSCs have underlying red stellar components with bluer inner components. Type-2 BSCs do not show an obvious sign of the underlying red stellar component, and their overall colors are quite blue ((U-B)_rest < 0). Both Type-1 and Type-2 BSCs have internal velocity dispersions measured from emission lines sigma < 80$ km/sec and estimated dynamical masses of only a few times 10^10 M_sun or less. For Type-1 BSCs, we estimate sigma of the red component using the fundamental plane relation of distant field absorption-line galaxies and find that these sigma estimates are similar to the sigma measured from emission lines. Overall, we conclude that our Type-1 and Type-2 BSCs are more likely to be star-forming low mass spheroids than star-forming, massive, early-type galaxies.

We review the motivation and search for supermassive black holes (BHs) in galaxies. Energetic nuclear activity provides indirect but compelling evidence for BH engines. Ground-based dynamical searches for central dark objects are reviewed in Kormendy & Richstone (1995, ARA&A, 33, 581). Here we provide an update of results from the HST. As of 2001 March, dynamical BH detections are available for at least 37 galaxies. The demographics of these objects lead to the following conclusions: (1) BH mass correlates with the luminosity of the bulge component of the host galaxy. The median BH mass fraction is 0.13% of the mass of the bulge. (2) BH mass correlates with the mean velocity dispersion of the bulge inside its effective radius. For the best mass determinations, the scatter is consistent with the measurement errors. (3) BH mass correlates with the luminosity of the high-density central component in disk galaxies independent of whether this is a real bulge (a mini-elliptical, believed to form via a merger-induced dissipative collapse and starburst) or a ``pseudobulge'' (believed to form by inward transport of disk material). (4) BH mass does not correlate with the luminosity of galaxy disks. If pure disks contain BHs, then their masses are much smaller than 0.13% of the mass of the disk. We conclude that present observations show no dependence of BH mass on the details of whether BH feeding happens rapidly during a collapse or slowly via secular evolution of the disk. The above results increasingly support the hypothesis that the major events that form a bulge or elliptical galaxy and the main growth phases of its BH -- when it shone like a quasar -- were the same events.

We use long-slit spectroscopic observations of the sample of E+A galaxies described by Zabludoff et al. to constrain the nature of the progenitors and remnants of the E+A phase of galaxy evolution. We measure spatially-resolved kinematic properties of the young (<~1 Gyr) and old (> few Gyr) stellar populations. The young stellar populations are more centrally concentrated than the older populations, but they are not confined to the galaxy core (radius <~ 1 kpc). The kinematics of the old stellar population place 16 of 20 of our E+As on a trend parallel to the Faber-Jackson relation that is offset by $\sim$ 0.6 mag in R. Eighteen of 20 E+As have v/sigma < 1. As the young stars in these systems evolve, the luminosity offset will disappear and the remnants will be pressure-supported systems that lie on the Faber-Jackson relation. Although Zabludoff et al. spectroscopically selected the most extreme E+A galaxies in the local volume, the sample is kinematically diverse: velocity dispersions range from <~ 30 km/s to ~ 200 km/s over a luminosity range of M_R = -19 to -22 + 5 log h. Combining these results with an estimate of the number of galaxies that experience an E+A phase, we conclude that the E+A phase of galaxy evolution is important in the development of a large fraction of spheroid-dominated galaxies over a wide range of luminosities and masses. Our kinematic observations, together with evidence that E+As have recently evolved from a vigorous star forming phase to a quiescent phase (e.g. Couch and Sharples, Caldwell et al.) and that many have tidal features consistent with disk-like progenitors (Zabludoff et al.), indicate that these galaxies are undergoing a transformation from star-forming, disk-dominated galaxies into quiescent, spheroid-dominated galaxies.

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 present the luminosity function and color-redshift relation of a magnitude-limited sample of 145 mostly red field E/S0 galaxies at z < 1 from the DEEP Groth Strip Survey (GSS). Most of the E/S0s (86%) form a red envelope in the redshift-color diagram, consistent with predictions of spectral synthesis models in which the dominant stellar population is formed at redshifts z > 1.5. Constructing a luminosity function of the full sample of 145 E/S0s, we find that there is about 1.1--1.9 magnitude brightening in rest-frame B band luminosity back to z = 0.8 from z=0, consistent with other studies. Together with the red colors, this brightening favors models in which the bulk of stars in red field E/S0s formed before z_{for} > 1.5 and have been evolving rather quiescently with few large starbursts since then. Evolution in the number density of field E/S0 galaxies is harder to measure, and uncertainties in the raw counts and their ratio to local samples might amount to as much as a factor of two. Within that uncertainty, the number density of red E/S0s to z = 0.8 seems relatively static, being comparable to or perhaps moderately less than that of local E/S0s depending on the assumed cosmology. A doubling of E/S0 number density since z = 1 can be ruled out with high confidence (97%) if \Omega_{m}=1. Taken together, our results are consistent with the hypothesis that the majority of luminous field E/S0s were already in place by z = 1, that the bulk of their stars were already fairly old, and that their number density has not changed by large amounts since then.

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 have used an Imaging Fabry-Perot Spectrophotometer with the Adaptive Optics Bonnette on the Canada-France-Hawaii Telescope to measure stellar radial velocities in the globular cluster M15. An average seeing of 0.15" full-width at half maximum, with the best-seeing image having 0.09", allowed us to measure accurately the velocities for five stars within 1" of the center of M15. Our estimate of the second moment of the velocity distribution inside a radius of 2" is 11.5 km/s, the same value we find out to a radius of about 6". However, the projected net rotation does increase dramatically at small radii, as our previous observations led us to suspect. The rotation amplitude inside a radius of 3.4" is v = 10.4 +- 2.7 km/s and the dispersion after removing the rotation is sigma = 10.3 +- 1.4 km/s, so v/sigma = 1 in this region. In addition, the position angle (PA) of the projected rotation axis differs by 100 degrees from that of the net cluster rotation at larger radii. Current theoretical models do not predict either this large an increase in the rotation amplitude or such a change in the PA. However, a central mass concentration, such as a black hole, could possibly sustain such a configuration. The rotation increase is consistent with the existence of a central dark mass concentration equal to 2500 M_solar.

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.

The [O II] lambda3727 emission lines and absorption features from stellar Balmer and Ca H&K lines are the most accessible kinematic diagnotics in galaxies at z~1. We investigate the kinematics of 22 local late-type galaxies using these spectral features and we compare the results to 21-cm neutral hydrogen spectra in order to assess the utility of each diagnostic for measuring galaxy masses. In order to simulate data at high redshift where only 1-dimensional velocity profiles are normally available, we study spatially integrated as well as spatially-resolved spectra. Although the studied galaxies span a wide range of morphological types, inclinations, and star formation rates, we find that the gaseous and stellar kinematic tracers yield comparable kinematic line widths and systemic velocities. The [O II] and H I line widths correlate most strongly, showing an intrinsic dispersion of ~20 km/s, or ~10% for a typical galaxy with a kinematic width of 200 km/s. In a few extreme cases, the [O II] line widths underestimate the neutral hydrogen width by 50%. Reliable velocity widths can also be obtained from the stellar Balmer and Ca H&K absorption lines, even for some of the very late-type galaxies which have strong emission lines. The intrinsic dispersion is <10% between the stellar absorption and H I line widths. We provide a prescription for using these strong stellar absorption and [O II] emission features to measure the kinematics, and thus masses, of galaxies in the distant universe.

We have compiled data for the globular cluster systems of 50 galaxies from the HST WFPC2 archive, of which 43 are type S0 or earlier. In this paper, we present the data set and derive the V-I color distributions. We derive the first four moments of the color distributions, as well as a measure for their non--unimodality. The number of globular clusters in each galaxy ranges from 18 (in NGC 2778) to 781 (NGC 5846). For those systems having more than 100 clusters, seven of sixteen (44%) show significant bimodality. Overall, roughly half of all the systems in our sample show hints of a bimodal color distribution. In general, the distributions of the faint galaxies are consistent with unimodality, whereas those of the brighter galaxies are not. In comparing the moments of the V-I distributions with various galaxy properties for the early-type galaxies, we find the following difference in the correlations between the field and cluster galaxy populations: the peak V-I color of the globular cluster distribution correlates well with the central velocity dispersion---and hence the Mg2 index and total luminosity---for galaxies in cluster environments; there exists no such correlation for field galaxies. This difference between cluster and field galaxies possibly reflects different formation scenarios for their globular cluster systems. Among the explanations for such a correlation, we consider either a larger age spread in the field populations or the possibility that cluster galaxies are always affected by significant accretion whereas some field galaxies could host pure ``in situ'' formed populations.

As part of the Deep Extragalactic Evolutionary Probe (DEEP) survey, a sample of 190 field galaxies (I_{814} <= 23.5) in the ``Groth Survey Strip'' has been used to analyze the magnitude-size relation over the range 0.1 < z < 1.1. The survey is statistically complete to this magnitude limit. All galaxies have photometric structural parameters, including bulge fractions (B/T), from Hubble Space Telescope images, and spectroscopic redshifts from the Keck Telescope. The analysis includes a determination of the survey selection function in the magnitude-size plane as a function of redshift, which mainly drops faint galaxies at large distances. Our results suggest that selection effects play a very important role. A first analysis treats disk-dominated galaxies with B/T < 0.5. If selection effects are ignored, the mean disk surface brightness (averaged over all galaxies) increases by ~1.3 mag from z = 0.1 to 0.9. However, most of this change is plausibly due to comparing low luminosity galaxies in nearby redshift bins to high luminosity galaxies in distant bins. If this effect is allowed for, no discernible evolution remains in the disk surface brightness of bright (M_B < -19) disk-dominated galaxies. A second analysis treats all galaxies by substituting half-light radius for disk scale length, with similar conclusions. Indeed, at all redshifts, the bulk of galaxies is consistent with the magnitude-size envelope of local galaxies, i.e., with little or no evolution in surface brightness. In the two highest redshift bins (z > 0.7), a handful of luminous, high surface brightness galaxies appears that occupies a region of the magnitude-size plane rarely populated by local galaxies. Their wide range of colors and bulge fractions points to a variety of possible origins.

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 revisit the kinematical data for 204 globular clusters in the halo of M87. Beyond 3 r_eff along the major axis of the galaxy light, these globular clusters exhibit substantial rotation (~ 300 +/- 70 km/s) that translates into an equally substantial spin (lambda ~ 0.18). The present appearance of M87 is most likely the product of a single major merger, since this event is best able to account for so sizable a spin. A rotation this large makes improbable any significant accretion of material after this merger, since that would have diluted the rotation signature. We see weak evidence for a difference between the kinematics of the metal-poor and metal-rich population, in the sense that the metal-poor globular clusters appear to dominate the rotation. If, as we suspect, the last major merger event of M87 was mainly dissipationless and did not trigger the formation of a large number of globular clusters, the kinematic difference between the two could reflect their orbital properties in the progenitor galaxies; these differences would be compatible with these progenitors having formed in dissipational mergers. However, to put strong kinematic constraints on the origin of the globular clusters themselves is difficult, given the complex history of the galaxy and its last dominant merger event.

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.

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).

We have used an Imaging Fabry-Perot Spectrophotometer with the Sub-arcsecond Imaging Spectrograph on the CFHT to measure velocities for 1534 stars in the globular cluster M15. Combined with previous velocity samples, the total number of stars with measured velocities in M15 is 1597. The velocity dispersion profile for M15 remains flat at a value of 11 km/s from a radius of 0.4' into our innermost reliable point at 0.02' (0.06 pc). This profile and the previously-published surface brightness profile can be equally well represented either by a stellar population whose M/L varies with radius from 1.7 in solar units at large radii to 3 in the central region, or by a population with a constant M/L of 1.7 and a central black hole of 1000 M_sol. A non-parametric mass model that assumes no black hole, no rotation, and isotropy constrains the mass density of M15 to better than 30% at a radius of 0.07 parsecs. Using the assumption of local thermodynamic equilibrium, we estimate the present-day mass function and infer a significant number of 0.6-0.7 M_sol objects in the central few parsecs, 85% of which may be in the form of stellar remnants. Not only do we detect rotation; we find that the position angle of the projected rotation axis in the central 10" is 100 degrees different from that of the whole sample. We also detect an increase in the amplitude of the rotation at small radii. Although this increase needs to be confirmed with better-seeing data, it may be the result of a central mass concentration.

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 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.

In spite of the promise of new techniques for constraining the mass distribution in galaxy clusters, much remains to be learned from galaxy orbital velocities. This article reviews the theory of potential estimation in hot dynamical systems like galaxy clusters. An analysis is presented of the Coma cluster, based on a sample of about 1500 galaxies with probable membership, of which 450 have measured velocities. We find no significant evidence for a core in either the galaxy number densities or the mass distribution.

B,V CCD frames have been used to derive surface brightness profiles for NGC 3201 out to ~18 arcmin. A total of 857 radial velocities with median precision ~1 km/s for 399 member giants have been used to trace the velocity dispersion profile out to 32' (the approximate tidal radius from fits of single-mass, isotropic King-Michie models to the cluster surface brightness profiles). The median difference in radial velocity for stars on either side of an imaginary axis moved through the cluster in 1 degree steps shows a significant maximum amplitude of 1.22+/-0.25 km/s. We discuss possible explanations of this result, including: (1) cluster rotation; (2) preferential stripping of stars on prograde orbits near the limiting radius; (3) the projection of the cluster space velocity onto the plane of the sky and (4) a slight drift in the velocity zero point. It is difficult to identify the primary cause of the observed velocity field structure unambiguously, and we suspect that all of the above processes may play a role. The B,V surface brightness profiles and radial velocities have been modeled with single- & multi-mass King-Michie models and nonparametric techniques. The density and M/L profiles show good agreement over 1.5<R<10 pc, and both approaches suggest a steady rise in M/L with distance from the cluster center. Due to the low cluster luminosity, we are unable to place useful constraints on the anisotropy of the velocity dispersion profile, though the global mass-to-light ratio is well-constrained by the models as ~2.0 +/-0.2 for the multi-mass and nonparametric models, compared to ~ 1.65 +/-0.15 for models having equal-mass stars. Our best-fit, multi-mass models have mass function slopes of x~0.75 +/-0.25, consistent with findings that mass function depends on the position relative to the potential of the Galaxy.

We use radial velocities of member stars and cluster surface brightness profiles to non-parametrically determine the mass density profiles and isotropic phase-space distribution functions $f(E)$ for the globular clusters M15 (NGC7078), 47~Tuc (NGC104), NGC~362, and NGC~3201. Assuming isotropy and using the velocity dispersion and surface brightness profiles, the Jeans equation uniquely determines the mass density profile. In the two centrally-concentrated clusters, M15 and 47~Tuc, we find that the mass-to-light ratios (M/L's) reach minima around 1\arcmin, and increase by more than a factor of four towards the cluster centers. For the two less centrally concentrated clusters, the M/L decreases monotonically all the way into the center. All four clusters exhibit an increase in the M/L's in their outer parts. If the variations in the M/L's are due to equipartition of energy between different mass stars, then we attribute the central increases to massive remnants and the outer increases to low-mass stars (m$<0.3$M$_\odot$). By applying the crude approximation of local thermodynamic equilibrium, we derive the present-day mass function for each cluster. In the central 2--3 parsecs, 0.7--1.5 M$_\odot$ objects provide the bulk of the cluster mass. The paper is in a uuencoded compressed tar file (260kb).

We report the first use of the Rutgers Imaging Fabry-Perot Spectrophotometer to study the dynamics of the cores of globular clusters. We have obtained velocities for cluster stars by tuning the Fabry-Perot to take a series of narrow-band images at different wavelengths across one of the Na D (5890 AA) absorption lines. Measuring the flux in every frame yields a short portion of the spectrum for each star simultaneously. This proves to be a very efficient method for obtaining accurate stellar velocities; in crowded regions we are able to measure hundreds of velocities in 3-4 hours of observing time. We have measured velocities with uncertainties of less than 5 km/s for 216 stars within 1.5' of the center of the globular cluster M15 (NGC 7078). The paper is a uuencoded compressed postscript file.