We introduce a new Markov Chain Monte Carlo (MCMC) algorithm with parallel tempering for fitting theoretical models of horizon-scale images of black holes to the interferometric data from the Event Horizon Telescope (EHT). The algorithm implements forms of the noise distribution in the data that are accurate for all signal-to-noise ratios. In addition to being trivially parallelizable, the algorithm is optimized for high performance, achieving 1 million MCMC chain steps in under 20 seconds on a single processor. We use synthetic data for the 2017 EHT coverage of M87 that are generated based on analytic as well as General Relativistic Magnetohydrodynamic (GRMHD) model images to explore several potential sources of biases in fitting models to sparse interferometric data. We demonstrate that a very small number of data points that lie near salient features of the interferometric data exert disproportionate influence on the inferred model parameters. We also show that the preferred orientations of the EHT baselines introduce significant biases in the inference of the orientation of the model images. Finally, we discuss strategies that help identify the presence and severity of such biases in realistic applications.

Interferometers, such as the Event Horizon Telescope (EHT), do not directly observe the images of sources but rather measure their Fourier components at discrete spatial frequencies up to a maximum value set by the longest baseline in the array. Construction of images from the Fourier components or analysis of them with high-resolution models requires careful treatment of fine source structure nominally beyond the array resolution. The primary EHT targets, Sgr A* and M87, are expected to have black-hole shadows with sharp edges and strongly filamentary emission from the surrounding plasma on scales much smaller than those probed by the currently largest baselines. We show that for aliasing not to affect images reconstructed with regularized maximum likelihood methods and model images that are directly compared to the data, the sampling of these images (i.e., their pixel spacing) needs to be significantly finer than the scale probed by the largest baseline in the array. Using GRMHD simulations of black-hole images, we estimate the maximum allowable pixel spacing to be approximately equal to (1/8)GM/c^2; for both of the primary EHT targets, this corresponds to an angular pixel size of <0.5 microarcseconds. With aliasing under control, we then advocate use of the second-order Butterworth filter with a cut-off scale equal to the maximum array baseline as optimal for visualizing the reconstructed images. In contrast to the traditional Gaussian filters, this Butterworth filter retains most of the power at the scales probed by the array while suppressing the fine image details for which no data exist.

The Kuiper Belt is a distant region of the Solar System. On 1 January 2019, the New Horizons spacecraft flew close to (486958) 2014 MU69, a Cold Classical Kuiper Belt Object, a class of objects that have never been heated by the Sun and are therefore well preserved since their formation. Here we describe initial results from these encounter observations. MU69 is a bi-lobed contact binary with a flattened shape, discrete geological units, and noticeable albedo heterogeneity. However, there is little surface color and compositional heterogeneity. No evidence for satellites, ring or dust structures, gas coma, or solar wind interactions was detected. By origin MU69 appears consistent with pebble cloud collapse followed by a low velocity merger of its two lobes.

Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

The Cold Classical Kuiper Belt, a class of small bodies in undisturbed orbits beyond Neptune, are primitive objects preserving information about Solar System formation. The New Horizons spacecraft flew past one of these objects, the 36 km long contact binary (486958) Arrokoth (2014 MU69), in January 2019. Images from the flyby show that Arrokoth has no detectable rings, and no satellites (larger than 180 meters diameter) within a radius of 8000 km, and has a lightly-cratered smooth surface with complex geological features, unlike those on previously visited Solar System bodies. The density of impact craters indicates the surface dates from the formation of the Solar System. The two lobes of the contact binary have closely aligned poles and equators, constraining their accretion mechanism.

The New Horizons spacecraft's encounter with the cold classical Kuiper belt object (486958) Arrokoth (formerly 2014 MU69) revealed a contact-binary planetesimal. We investigate how it formed, finding it is the product of a gentle, low-speed merger in the early Solar System. Its two lenticular lobes suggest low-velocity accumulation of numerous smaller planetesimals within a gravitationally collapsing, solid particle cloud. The geometric alignment of the lobes indicates the lobes were a co-orbiting binary that experienced angular momentum loss and subsequent merger, possibly due to dynamical friction and collisions within the cloud or later gas drag. Arrokoth's contact-binary shape was preserved by the benign dynamical and collisional environment of the cold classical Kuiper belt, and so informs the accretion processes that operated in the early Solar System.

On January 1st 2019, the New Horizons spacecraft flew by the classical Kuiper belt object (486958) Arrokoth (provisionally designated 2014 MU69), possibly the most primitive object ever explored by a spacecraft. The I/F of Arrokoth is analyzed and fit with a photometric function that is a linear combination of the Lommel-Seeliger (lunar) and Lambert photometric functions. Arrokoth has a geometric albedo of p_V = 0.21_(-0.04)^(+0.05) at a wavelength of 550 nm and ~0.24 at 610 nm. Arrokoth's geometric albedo is greater than the median but consistent with a distribution of cold classical Kuiper belt objects whose geometric albedos were determined by fitting a thermal model to radiometric observations. Thus, Arrokoth's geometric albedo adds to the orbital and spectral evidence that it is a cold classical Kuiper belt object. Maps of the normal reflectance and hemispherical albedo of Arrokoth are presented. The normal reflectance of Arrokoth's surface varies with location, ranging from ~0.10-0.40 at 610 nm with an approximately Gaussian distribution. Both Arrokoth's extrema dark and extrema bright surfaces are correlated to topographic depressions. Arrokoth has a bilobate shape and the two lobes have similar normal reflectance distributions: both are approximately Gaussian, peak at ~0.25 at 610 nm, and range from ~0.10-0.40, which is consistent with co-formation and co-evolution of the two lobes. The hemispherical albedo of Arrokoth varies substantially with both incidence angle and location, the average hemispherical albedo at 610 nm is 0.063 +/- 0.015. The Bond albedo of Arrokoth at 610 nm is 0.062 +/- 0.015.

The LOng Range Reconnaissance Imager (LORRI) is a panchromatic (360--910 nm), narrow-angle (field of view = 0.29 deg), high spatial resolution (pixel scale = 1.02 arcsec) visible light imager used on NASA's New Horizons (NH) mission for both science observations and optical navigation. Calibration observations began several months after the NH launch on 2006 January 19 and have been repeated annually throughout the course of the mission, which is ongoing. This paper describes the in-flight LORRI calibration measurements, and the results derived from our analysis of the calibration data. LORRI has been remarkably stable over time with no detectable changes (at the 1% level) in sensitivity or optical performance since launch. By employing 4 by 4 re-binning of the CCD pixels during read out, a special spacecraft tracking mode, exposure times of 30 sec, and co-addition of approximately 100 images, LORRI can detect unresolved targets down to V = 22 (SNR=5). LORRI images have an instantaneous dynamic range of 3500, which combined with exposure time control ranging from 0ms to 64,967 ms in 1ms steps supports high resolution, high sensitivity imaging of planetary targets spanning heliocentric distances from Jupiter to deep in the Kuiper belt, enabling a wide variety of scientific investigations. We describe here how to transform LORRI images from raw (engineering) units into scientific (calibrated) units for both resolved and unresolved targets. We also describe various instrumental artifacts that could affect the interpretation of LORRI images under some observing circumstances.

The flyby of Pluto and Charon by the New Horizons spacecraft provided high-resolution images of cratered surfaces embedded in the Kuiper belt, an extensive region of bodies orbiting beyond Neptune. Impact craters on Pluto and Charon were formed by collisions with other Kuiper belt objects (KBOs) with diameters from ~40 kilometers to ~300 meters, smaller than most KBOs observed directly by telescopes. We find a relative paucity of small craters less than approximately 13 kilometers in diameter, which cannot be explained solely by geological resurfacing. This implies a deficit of small KBOs (less than 1 to 2 kilometers in diameter). Some surfaces on Pluto and Charon are likely greater than 4 billion years old, thus their crater records provide information on the size-frequency distribution of KBOs in the early Solar System.

The deep nitrogen-covered Sputnik Planitia (SP; informal name) basin on Pluto is located very close to the longitude of Pluto's tidal axis[1] and may be an impact feature [2], by analogy with other large basins in the solar system[3,4]. Reorientation[5-7] due to tidal and rotational torques can explain SP's location, but requires it to be a positive gravity anomaly[7], despite its negative topography. Here we argue that if SP formed via impact and if Pluto possesses a subsurface ocean, a positive gravity anomaly would naturally result because of shell thinning and ocean uplift, followed by later modest N2 deposition. Without a subsurface ocean a positive gravity anomaly requires an implausibly thick N2 layer (greater than 40 km). A rigid, conductive ice shell is required to prolong such an ocean's lifetime to the present day[8] and maintain ocean uplift. Because N2 deposition is latitude-dependent[9], nitrogen loading and reorientation may have exhibited complex feedbacks[7].

The vast, deep, volatile-ice-filled basin informally named Sputnik Planum is central to Pluto's geological activity[1,2]. Composed of molecular nitrogen, methane, and carbon monoxide ices[3], but dominated by N2-ice, this ice layer is organized into cells or polygons, typically ~10-40 km across, that resemble the surface manifestation of solid state convection[1,2]. Here we report, based on available rheological measurements[4], that solid layers of N2 ice approximately greater than 1 km thick should convect for estimated present-day heat flow conditions on Pluto. More importantly, we show numerically that convective overturn in a several-km-thick layer of solid nitrogen can explain the great lateral width of the cells. The temperature dependence of N2-ice viscosity implies that the SP ice layer convects in the so-called sluggish lid regime[5], a unique convective mode heretofore not definitively observed in the Solar System. Average surface horizontal velocities of a few cm/yr imply surface transport or renewal times of ~500,000 years, well under the 10 Myr upper limit crater retention age for Sputnik Planum[2]. Similar convective surface renewal may also occur on other dwarf planets in the Kuiper belt, which may help explain the high albedos of some of them.

Pluto's atmospheric haze settles out rapidly compared with geological timescales. It needs to be accounted for as a surface material, distinct from Pluto's icy bedrock and from the volatile ices that migrate via sublimation and condensation on seasonal timescales. This paper explores how a steady supply of atmospheric haze might affect three distinct provinces on Pluto. We pose the question of why they each look so different from one another if the same haze material is settling out onto all of them. Cthulhu is a more ancient region with comparatively little present-day geological activity, where the haze appears to simply accumulate over time. Sputnik Planitia is a very active region where glacial convection, as well as sublimation and condensation rapidly refresh the surface, hiding recently deposited haze from view. Lowell Regio is a region of intermediate age featuring very distinct coloration from the rest of Pluto. Using a simple model haze particle as a colorant, we are not able to match the colors in both Lowell Regio and Cthulhu. To account for their distinct colors, we propose that after arrival at Pluto's surface, haze particles may be less inert than might be supposed from the low surface temperatures. They must either interact with local materials and environments to produce distinct products in different regions, or else the supply of haze must be non-uniform in time and/or location, such that different products are delivered to different places.

A unique feature of Pluto's large satellite Charon is its dark red northern polar cap. Similar colours on Pluto's surface have been attributed to organic macromolecules produced by energetic radiation processing of hydrocarbons. The polar location of this material on Charon implicates the temperature extremes that result from Charon's high obliquity and long seasons. The escape of Pluto's atmosphere provides a potential feed stock for production of complex chemistry. Gas from Pluto that is transiently cold-trapped and processed at Charon's winter pole was proposed as an explanation on the basis of an image of Charon's northern hemisphere, but not modelled quantitatively. Here we report images of the southern hemisphere illuminated by Pluto-shine and also images taken during the approach phase showing the northern polar cap over a range of longitudes. We model the surface thermal environment on Charon, the supply and temporary cold-trapping of material escaping from Pluto, and, while cold-trapped, its photolytic processing into more complex and less volatile molecules. The model results are consistent with the proposed mechanism producing the observed colour pattern on Charon.

The DESI Legacy Imaging Surveys are a combination of three public projects (the Dark Energy Camera Legacy Survey, the Beijing-Arizona Sky Survey, and the Mayall z-band Legacy Survey) that will jointly image approximately 14,000 deg^2 of the extragalactic sky visible from the northern hemisphere in three optical bands (g, r, and z) using telescopes at the Kitt Peak National Observatory and the Cerro Tololo Inter-American Observatory. The combined survey footprint is split into two contiguous areas by the Galactic plane. The optical imaging is conducted using a unique strategy of dynamically adjusting the exposure times and pointing selection during observing that results in a survey of nearly uniform depth. In addition to calibrated images, the project is delivering a catalog, constructed by using a probabilistic inference-based approach to estimate source shapes and brightnesses. The catalog includes photometry from the grz optical bands and from four mid-infrared bands (at 3.4, 4.6, 12 and 22 micorons) observed by the Wide-field Infrared Survey Explorer (WISE) satellite during its full operational lifetime. The project plans two public data releases each year. All the software used to generate the catalogs is also released with the data. This paper provides an overview of the Legacy Surveys project.

The New Horizons encounter with the cold classical Kuiper Belt object (KBO) 2014 MU69 (informally named 'Ultima Thule,' hereafter Ultima) on 1 January 2019 will be the first time a spacecraft has ever closely observed one of the free-orbiting small denizens of the Kuiper Belt. Related to but not thought to have formed in the same region of the Solar System as the comets that been explored so far, it will also be the largest, most distant, and most primitive body yet visited by spacecraft. In this letter we begin with a brief overview of cold classical KBOs, of which Ultima is a prime example. We give a short preview of our encounter plans. We note what is currently known about Ultima from earth-based observations. We then review our expectations and capabilities to evaluate Ultima's composition, surface geology, structure, near space environment, small moons, rings, and the search for activity.

We present the study of stellar populations in the central 5.5' (~1.2 kpc) of the M31 bulge by using the optical color magnitude diagram derived from HST ACS WFC/HRC observations. In order to enhance image quality and then obtain deeper photometry, we construct Nyquist-sampled images and use a deconvolution method to detect sources and measure their photometry. We demonstrate that our method performs better than DOLPHOT in the extremely crowded region. The resolved stars in the M31 bulge have been divided into nine annuli and the color magnitude diagram fitting is performed for each of them. We confirm that the majority of stars (> 70%) in the M31 bulge are indeed very old (>5 Gyr) and metal-rich ([Fe/H] > 0.3). At later times, the star formation rate decreased and then experienced a significant rise around 1 Gyr ago, which pervaded the entire M31 bulge. After that, stars formed at less than 500 Myr ago in the central 130". Through simulation, we find that these intermediate-age stars cannot be the artifacts introduced by the blending effect. Our results suggest that although the majority of the M31 bulge are very old, the secular evolutionary process still continuously builds up the M31 bulge slowly. We compare our star formation history with an older analysis derived from the spectral energy distribution fitting, which suggests that the latter one is still a reasonable tool for the study of stellar populations in remote galaxies.

We explore the use of principal component analysis (PCA) to characterize high-fidelity simulations and interferometric observations of the millimeter emission that originates near the horizons of accreting black holes. We show mathematically that the Fourier transforms of eigenimages derived from PCA applied to an ensemble of images in the spatial-domain are identical to the eigenvectors of PCA applied to the ensemble of the Fourier transforms of the images, which suggests that this approach may be applied to modeling the sparse interferometric Fourier-visibilities produced by an array such as the Event Horizon Telescope (EHT). We also show that the simulations in the spatial domain themselves can be compactly represented with a PCA-derived basis of eigenimages allowing for detailed comparisons between variable observations and time-dependent models, as well as for detection of outliers or rare events within a time series of images. Furthermore, we demonstrate that the spectrum of PCA eigenvalues is a diagnostic of the power spectrum of the structure and, hence, of the underlying physical processes in the simulated and observed images.

Ten weeks of daily imaging of the giant elliptical galaxy M87 with the Hubble Space Telescope (HST) has yielded 41 nova light curves of unprecedented quality for extragalactic cataclysmic variables. We have recently used these light curves to demonstrate that the observational scatter in the so-called Maximum-Magnitude Rate of Decline (MMRD) relation for classical novae is so large as to render the nova-MMRD useless as a standard candle. Here we demonstrate that a modified Buscombe - de Vaucouleurs hypothesis, namely that novae with decline times t2 > 10 days converge to nearly the same absolute magnitude about two weeks after maximum light in a giant elliptical galaxy, is supported by our M87 nova data. For 13 novae with daily-sampled light curves, well determined times of maximum light in both the F606W and F814W filters, and decline times $t2 > 10 days we find that M87 novae display M(606W,15) = -6.37 +/- 0.46 and M(814W,15) = -6.11 +/- 0.43. If very fast novae with decline times t2 < 10 days are excluded, the distances to novae in elliptical galaxies with stellar binary populations similar to those of M87 should be determinable with 1 sigma accuracies of +/-20% with the above calibrations.

A search for temporal changes on Pluto and Charon was motivated by (1) the discovery of young surfaces in the Pluto system that imply ongoing or recent geologic activity, (2) the detection of active plumes on Triton during the Voyager 2 flyby, and (3) the abundant and detailed information that observing geologic processes in action provides about the processes. A thorough search for temporal changes using New Horizons images was completed. Images that covered the same region were blinked and manually inspected for any differences in appearance. The search included full-disk images such that all illuminated regions of both bodies were investigated and higher resolution images such that parts of the encounter hemispheres were investigated at finer spatial scales. Changes of appearance between different images were observed but in all cases were attributed to variability of the imaging parameters (especially geometry) or artifacts. No differences of appearance that are strongly indicative of a temporal change were found on the surface or in the atmosphere of either Pluto or Charon. Limits on temporal changes as a function of spatial scale and temporal interval during the New Horizons encounter are determined. The longest time interval constraint is one Pluto/Charon rotation period (~6.4 Earth days). Contrast reversal and high-phase bright features that change in appearance with solar phase angle are identified. The change of appearance of these features is most likely due to the change in phase angle rather than a temporal change. Had active plumes analogous to the plumes discovered on Triton been present on the encounter hemispheres of either Pluto or Charon, they would have been detected. The absence of active plumes may be due to temporal variability (i.e., plumes do occur but none were active on the encounter hemispheres during the epoch of the New Horizons encounter ...

We searched for dust or debris rings in the Pluto-Charon system before, during, and after the New Horizons encounter. Methodologies included searching for back-scattered light during the approach to Pluto (phase $\sim15^\circ$), in situ detection of impacting particles, a search for stellar occultations near the time of closest approach, and by forward-scattered light during departure (phase $\sim165^\circ$). A search using HST prior to the encounter also contributed to the results. No rings, debris, or dust features were observed, but our detection limits provide an improved picture of the environment throughout the Pluto-Charon system. Searches for rings in back-scattered light covered 35,000-250,000 km from the system barycenter, a zone that starts interior to the orbit of Styx, and extends to four times the orbital radius of Hydra. We obtained our firmest limits using the NH LORRI camera in the inner half of this region. Our limits on the normal $I/F$ of an unseen ring depends on the radial scale of the rings: $2\times10^{-8}$ ($3\sigma$) for 1500 km wide rings, $1\times10^{-8}$ for 6000 km rings, and $7\times10^{-9}$ for 12,000 km rings. Beyond $\sim100,000$ km from Pluto, HST observations limit normal $I/F$ to $\sim8\times10^{-8}$. Searches for dust from forward-scattered light extended from the surface of Pluto to the Pluto-Charon Hill sphere ($r_{\rm Hill}=6.4\times10^6$ km). No evidence for rings or dust was detected to normal $I/F$ limits of $\sim8.9\times10^{-7}$ on $\sim10^4$ km scales. Four occulation observations also probed the space interior to Hydra, but again no dust or debris was detected. Elsewhere in the solar system, small moons commonly share their orbits with faint dust rings. Our results suggest that small grains are quickly lost from the system due to solar radiation pressure, whereas larger particles are unstable due to perturbations by the known moons.

We present VLA images and HST/STIS spectra of sources within the center of the brightest cluster galaxy (BCG) in Abell 2261. These observations were obtained to test the hypothesis that its extremely large, flat core reflects the ejection of its supermassive black hole. Spectra of three of the four most luminous "knots" embedded in the core were taken to test whether one may represent stars bound to a displaced massive black hole. The three knots have radial velocity offsets dV < ~150 km/s from the BCG. Knots 2 and 3 show kinematics, colors, and stellar masses consistent with infalling low-mass galaxies or larger stripped cluster members. Large errors in the stellar velocity dispersion of Knot 1, however, mean that we cannot rule out the hypothesis that it hosts a high-mass black hole. A2261-BCG has a compact, relic radio-source offset by 6.5 kpc (projected) from the optical core's center, but no active radio core that would pinpoint the galaxy's central black hole to a tight 10 GHz flux limit <3.6 uJy. Its spectrum and morphology are suggestive of an active galactic nucleus that switched off >48 Myr ago, with an equipartition condition magnetic field of 15 uG. These observations are still consistent with the hypothesis that the nuclear black hole has been ejected from its core, but the critical task of locating the supermassive black hole or demonstrating that A2261-BCG lacks one remains to be done.

Several trace constituents can reach saturation vapor pressure in Pluto's present day atmosphere. As such, we describe a search for discrete could features in Pluto's atmosphere using New Horizons data obtained on 14-15 July 2015, during the Pluto flyby closest approach. We report that Pluto's present day atmosphere is at least largely (>99% by surface area) free of discrete clouds. We also report a handful of features that may plausibly be clouds, all of which were detected obliquely and at high phase observing geometry. No cloud candidates were identified away from the terminators or in low phase (backscattering geometry) images.

The extensive grid of numerical simulations of nova eruptions of Yaron et al.(2005) first predicted that some classical novae might deviate significantly from the Maximum Magnitude - Rate of Decline (MMRD) relation, which purports to characterise novae as standard candles. Kasliwal et al. (2011) have announced the observational detection of an apparently new class of faint, fast classical novae in the Andromeda galaxy. These objects deviate strongly from the MMRD relationship, exactly as predicted by Yaron et al. (2005). Shara et al. (2016) recently reported the first detections of faint, fast novae in M87. These previously overlooked objects are as common in the giant elliptical galaxy M87 as they are in the giant spiral M31; they comprise about 40% of all classical nova eruptions and greatly increase the observational scatter in the MMRD relation. We use the extensive grid of nova simulations of Yaron et al. (2005) to identify the underlying causes of the existence of faint, fast novae. These are systems which have accreted, and can thus eject, only very low mass envelopes, of order 10^-7 - 10^-8 Msun, on massive white dwarfs. Such binaries include, but are not limited to, the recurrent novae. These same models predict the existence of ultrafast novae which display decline times t2 as short as five hours. We outline a strategy for their future detection.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

NASA's New Horizons spacecraft observed (15810) 1994 JR$_1$, a 3:2 resonant Kupier Belt Object (KBO), using the LOng Range Reconnaissance Imager (LORRI) on November 2, 2015 from a distance of 1.85 AU, and again on April 7, 2016 from a distance of 0.71 AU. These were the first close observations of any KBO other than Pluto. Combining ground-based and Hubble Space Telecope (HST) observations at small phase angles and the LORRI observations at higher phase angles, we produced the first disk-integrated solar phase curve of a typical KBO from $\alpha$=0.6-58$^\circ$. Observations at these geometries, attainable only from a spacecraft in the outer Solar System, constrain surface properties such as macroscopic roughness and the single particle phase function. 1994 JR$_1$ has a rough surface with a 37$\pm$5$^\circ$ mean topographic slope angle and has a relatively rapid rotation period of 5.47$\pm$0.33 hours. 1994 JR$_1$ is currently 2.7 AU from Pluto; our astrometric points enable high-precision orbit determination and integrations which show that it comes this close to Pluto every 2.4 million years (10$^4$ heliocentric orbits), causing Pluto to perturb 1994 JR$_1$. During the November spacecraft observation, the KBO was simultaneously observed using HST in two colors, confirming its very red spectral slope. These observations have laid the groundwork for numerous potential future distant KBO observations in the New Horizons-Kuiper Belt Extended Mission.

The Hubble Space Telescope has imaged the central part of M87 over a 10 week span, leading to the discovery of 32 classical novae and nine fainter, likely very slow and/or symbiotic novae. In this first in a series of papers we present the M87 nova finder charts, and the light and color curves of the novae. We demonstrate that the rise and decline times, and the colors of M87 novae are uncorrelated with each other and with position in the galaxy. The spatial distribution of the M87 novae follows the light of the galaxy, suggesting that novae accreted by M87 during cannibalistic episodes are well-mixed. Conservatively using only the 32 brightest classical novae we derive a nova rate for M87: $363_{-45}^{+33}$ novae/yr. We also derive the luminosity-specific classical nova rate for this galaxy, which is $7.88_{-2.6}^{+2.3} /yr/ 10^{10}L_\odot,_{K}$. Both rates are 3-4 times higher higher than those reported for M87 in the past, and similarly higher than those reported for all other galaxies. We suggest that most previous ground-based surveys for novae in external galaxies, including M87, miss most faint, fast novae, and almost all slow novae near the centers of galaxies.

We report on a variety of standard techniques used by New Horizons including a solar ultraviolet occultation, ultraviolet airglow observations, and high-phase look-back particulate search imaging to search for an atmosphere around Pluto's large moon Charon during its flyby in July 2015. Analyzing these datasets, no evidence for a present day atmosphere has been found for 14 potential atomic and molecular species, all of which are now constrained to have pressures below 0.3 nanobar, as we describe below, these are much more stringent upper limits than the previously available 15-110 nanobar constraints (e.g., Sicardy et al. 2006); for example, we find a 3$\sigma$ upper limit for an N$_2$ atmosphere on Charon is 4.2 picobars and a 3$\sigma$ upper limit for the brightness of any atmospheric haze on Charon of I/F=2.6x10$^{-5}$. A radio occultation search for an atmosphere around Charon was also conducted by New Horizons but will be published separately by other authors.

Approach images taken by the LORRI imaging system during the New Horizons spacecraft encounter have been used to determine the mean radii and shapes of Pluto and Charon. The primary observations are limb locations derived using three independent approaches. The resulting mean radii of Pluto and Charon are 1188.3 +/- 1.6 km and 606.0 +/- 1.0 km, respectively (2-sigma). The corresponding densities are 1854 +/- 11 kg m-3 and 1701 +/- 33 kg m-3 (2-sigma). The Charon radius value is consistent with previous Earth-based occultation estimates. The Pluto radius estimate is consistent with solar occultation measurements performed by the ALICE and Fine Sun Sensor instruments on New Horizons. Neither Pluto nor Charon show any evidence for tidal/rotational distortions; upper bounds on the oblateness are <0.6% and <0.5%, respectively.

NASA's New Horizons spacecraft has revealed the complex geology of Pluto and Charon. Pluto's encounter hemisphere shows ongoing surface geological activity centered on a vast basin containing a thick layer of volatile ices that appears to be involved in convection and advection, with a crater retention age no greater than $\approx$10 Ma. Surrounding terrains show active glacial flow, apparent transport and rotation of large buoyant water-ice crustal blocks, and pitting, likely by sublimation erosion and/or collapse. More enigmatic features include tall mounds with central depressions that are conceivably cryovolcanic, and ridges with complex bladed textures. Pluto also has ancient cratered terrains up to ~4 Ga old that are extensionally fractured and extensively mantled and perhaps eroded by glacial or other processes. Charon does not appear to be currently active, but experienced major extensional tectonism and resurfacing (probably cryovolcanic) nearly 4 billion years ago. Impact crater populations on Pluto and Charon are not consistent with the steepest proposed impactor size-frequency distributions proposed for the Kuiper belt.

The New Horizons mission has provided resolved measurements of Pluto's moons Styx, Nix, Kerberos, and Hydra. All four are small, with equivalent spherical diameters of $\approx$40 km for Nix and Hydra and ~10 km for Styx and Kerberos. They are also highly elongated, with maximum to minimum axis ratios of $\approx$2. All four moons have high albedos ( $\approx$50-90 %) suggestive of a water-ice surface composition. Crater densities on Nix and Hydra imply surface ages $\gtrsim$ 4 Ga. The small moons rotate much faster than synchronous, with rotational poles clustered nearly orthogonal to the common pole directions of Pluto and Charon. These results reinforce the hypothesis that the small moons formed in the aftermath of a collision that produced the Pluto-Charon binary.

We map the dust distribution in the central 180" (~680 pc) region of the M31 bulge, based on HST/WFC3 and ACS observations in ten bands from near-ultraviolet (2700 A) to near-infrared (1.5 micron). This large wavelength coverage gives us great leverage to detect not only dense dusty clumps, but also diffuse dusty molecular gas. We fit a pixel-by-pixel spectral energy distributions to construct a high-dynamic-range extinction map with unparalleled angular resolution (~0.5" , i.e., ~2 pc) and sensitivity (the extinction uncertainty, \delta A_V~0.05). In particular, the data allow to directly fit the fractions of starlight obscured by individual dusty clumps, and hence their radial distances in the bulge. Most of these clumps seem to be located in a thin plane, which is tilted with respect to the M31 disk and appears face-on. We convert the extinction map into a dust mass surface density map and compare it with that derived from the dust emission as observed by Herschel . The dust masses in these two maps are consistent with each other, except in the low-extinction regions, where the mass inferred from the extinction tends to be underestimated. Further, we use simulations to show that our method can be used to measure the masses of dusty clumps in Virgo cluster early-type galaxies to an accuracy within a factor of ~2.

The Pluto system was recently explored by NASA's New Horizons spacecraft, making closest approach on 14 July 2015. Pluto's surface displays diverse landforms, terrain ages, albedos, colors, and composition gradients. Evidence is found for a water-ice crust, geologically young surface units, surface ice convection, wind streaks, volatile transport, and glacial flow. Pluto's atmosphere is highly extended, with trace hydrocarbons, a global haze layer, and a surface pressure near 10 microbars. Pluto's diverse surface geology and long-term activity raise fundamental questions about how small planets remain active many billions of years after formation. Pluto's large moon Charon displays tectonics and evidence for a heterogeneous crustal composition, its north pole displays puzzling dark terrain. Small satellites Hydra and Nix have higher albedos than expected.

We map the distribution of dust in M31 at 25pc resolution, using stellar photometry from the Panchromatic Hubble Andromeda Treasury. We develop a new mapping technique that models the NIR color-magnitude diagram (CMD) of red giant branch (RGB) stars. The model CMDs combine an unreddened foreground of RGB stars with a reddened background population viewed through a log-normal column density distribution of dust. Fits to the model constrain the median extinction, the width of the extinction distribution, and the fraction of reddened stars. The resulting extinction map has >4 times better resolution than maps of dust emission, while providing a more direct measurement of the dust column. There is superb morphological agreement between the new map and maps of the extinction inferred from dust emission by Draine et al. 2014. However, the widely-used Draine & Li (2007) dust models overpredict the observed extinction by a factor of ~2.5, suggesting that M31's true dust mass is lower and that dust grains are significantly more emissive than assumed in Draine et al. (2014). The discrepancy we identify is consistent with similar findings in the Milky Way by the Planck Collaboration (2015), but has a more complex dependence on parameters from the Draine & Li (2007) dust models. We also show that the discrepancy with the Draine et al. (2014) map is lowest where the interstellar radiation field has a harder spectrum than average. We discuss possible improvements to the CMD dust mapping technique, and explore further applications.

We explore the assembly history of the M31 bulge within a projected major-axis radius of 180" (~680 pc) by studying its stellar populations in Hubble Space Telescope WFC3 and ACS observations. Colors formed by comparing near-ultraviolet vs. optical bands are found to become bluer with increasing major-axis radius, which is opposite to that predicted if the sole sources of near-ultraviolet light were old extreme horizontal branch stars with a negative radial gradient in metallicity. Spectral energy distribution fits require a metal-rich intermediate-age stellar population (300 Myr to 1 Gyr old, ~solar metallicity) in addition to the dominant old population. The radial gradients in age and metallicity of the old stellar population are consistent with those in previous works. For the intermediate-age population, we find an increase in age with radius and a mass fraction that increases up to 2% at 680 pc away from the center. We exclude contamination from the M31 disk and/or halo as the main origin for this population. Our results thus suggest that intermediate-age stars exist beyond the central 5" (19 pc) of M31 and contribute ~1% of the total stellar mass in the bulge. These stars could be related to the secular growth of the M31 bulge.

We have analyzed new HST/ACS and HST/WFC3 imaging in F475W and F814W of two previously-unobserved fields along the M31 minor axis to confirm our previous constraints on the shape of M31's inner stellar halo. Both of these new datasets reach a depth of at least F814W$<$27 and clearly detect the blue horizontal branch (BHB) of the field as a distinct feature of the color-magnitude diagram. We measure the density of BHB stars and the ratio of BHB to red giant branch stars in each field using identical techniques to our previous work. We find excellent agreement with our previous measurement of a power-law for the 2-D projected surface density with an index of 2.6$^{+0.3}_{-0.2}$ outside of 3 kpc, which flattens to $\alpha <$1.2 inside of 3 kpc. Our findings confirm our previous suggestion that the field BHB stars in M31 are part of the halo population. However, the total halo profile is now known to differ from this BHB profile, which suggests that we have isolated the metal-poor component. This component appears to have an unbroken power-law profile from 3-150 kpc but accounts for only about half of the total halo stellar mass. Discrepancies between the BHB density profile and other measurements of the inner halo are therefore likely due to the different profile of the metal-rich halo component, which is not only steeper than the profile of the met al-poor component, but also has a larger core radius. These profile differences also help to explain the large ratio of BHB/RGB stars in our observations.

We have observed 433 z<=0.08 brightest cluster galaxies (BCGs) in a full-sky survey of Abell clusters. The BCG Hubble diagram is consistent to within 2% of a Omega_m=0.3, Lambda=0.7 Hubble relation. The L_m-alpha relation for BCGs, which uses alpha, the log-slope of the BCG photometric curve of growth, to predict metric luminosity, L_m, has 0.27 mag residuals. We measure central stellar velocity dispersions, sigma, of the BCGs, finding the Faber-Jackson relation to flatten as the metric aperture grows to include an increasing fraction of the total BCG luminosity. A 3-parameter "metric plane" relation using alpha and sigma together gives the best prediction of L_m, with 0.21 mag residuals. The projected spatial offset, r_x, of BCGs from the X-ray-defined cluster center is a gamma=-2.33 power-law over 1<r_x<10^3 kpc. The median offset is ~10 kpc, but ~15% of the BCGs have r_x>100 kpc. The absolute cluster-dispersion normalized BCG peculiar velocity |Delta V_1|/sigma_c follows an exponential distribution with scale length 0.39+/-0.03. Both L_m and alpha increase with sigma_c. The alpha parameter is further moderated by both the spatial and velocity offset from the cluster center, with larger alpha correlated with the proximity of the BCG to the cluster mean velocity or potential center. At the same time, position in the cluster has little effect on L_m. The luminosity difference between the BCG and second-ranked galaxy, M2, increases as the peculiar velocity of the BCG within the cluster decreases. Further, when M2 is a close luminosity "rival" of the BCG, the galaxy that is closest to either the velocity or X-ray center of the cluster is most likely to have the larger alpha. We conclude that the inner portions of the BCGs are formed outside the cluster, but interactions in the heart of the galaxy cluster grow and extend the envelopes of the BCGs.

We have measured stellar photometry with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) and Advanced Camera for Surveys (ACS) in near ultraviolet (F275W, F336W), optical (F475W, F814W), and near infrared (F110W, F160W) bands for 117 million resolved stars in M31. As part of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, we measured photometry with simultaneous point spread function fitting across all bands and at all source positions after precise astrometric image alignment (<5-10 milliarcsecond accuracy). In the outer disk, the photometry reaches a completeness-limited depth of F475W~28, while in the crowded, high surface brightness bulge, the photometry reaches F475W~25. We find that simultaneous photometry and optimized measurement parameters significantly increase the detection limit of the lowest resolution filters (WFC3/IR) providing color-magnitude diagrams that are up to 2.5 magnitudes deeper when compared with color-magnitude diagrams from WFC3/IR photometry alone. We present extensive analysis of the data quality including comparisons of luminosity functions and repeat measurements, and we use artificial star tests to quantify photometric completeness, uncertainties and biases. We find that largest sources of systematic error in the photometry are due to spatial variations in the point spread function models and charge transfer efficiency corrections. This stellar catalog is the largest ever produced for equidistant sources, and is publicly available for download by the community.

We measure the extinction curve in the central 200 pc of M31 at mid-ultraviolet to near-infrared wavelengths (from 1928A to 1.5{\mu}m), using Swift/UVOT and HST WFC3/ACS observations in thirteen bands. Taking advantage of the high angular resolution of the HST WFC3 and ACS detectors, we develop a method to simultaneously determine the relative extinction and the fraction of obscured starlight for five dusty complexes located in the circumnuclear region. The extinction curves of these clumps (RV =2.4-2.5) are steeper than the average Galactic one (RV =3.1), but are similar to optical and near-infrared curves recently measured toward the Galactic Bulge (RV~2.5). This similarity suggests that steep extinction curves may be common in the inner bulge of galaxies. In the ultraviolet, the extinction curves of these clumps are also unusual. We find that one dusty clump (size <2 pc) exhibits a strong UV bump (extinction at 2175A), more than three standard deviation higher than that predicted by common models. Although the high stellar metallicity of the M31 bulge indicates that there are sufficient carbon and silicon to produce large dust grains, the grains may have been destroyed by supernova explosions or past activity of the central super-massive black hole, resulting in the observed steepened extinction curve.

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.

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

I have combined the Emsellem et al. ATLAS3D rotation measures of a large sample of early-type galaxies with HST-based classifications of their central structure to characterize the rotation velocities of galaxies with cores. "Core galaxies" rotate slowly, while "power-law galaxies" (galaxies that lack cores) rotate rapidly, confirming the analysis of Faber et al. Significantly, the amplitude of rotation sharply discriminates between the two types in the -19 > Mv > -22 domain over which the two types coexist. The slow rotation in the small set of core galaxies with Mv > -20, in particular, brings them into concordance with the more massive core galaxies. The ATLAS3D "fast-rotating" and "slow-rotating" early-type galaxies are essentially the same as power-law and core galaxies, respectively, or the Kormendy & Bender two families of elliptical galaxies based on rotation, isophote shape, and central structure. The ATLAS3D fast rotators do include roughly half of the core galaxies, but their rotation-amplitudes are always at the lower boundary of that subset. Essentially all core galaxies have ATLAS3D rotation-amplitudes lambda_(R_e/2) <= 0.25, while all galaxies with lambda_(R_e/2) > 0.25 and figure eccentricity > 0.2 lack cores. Both figure rotation and the central structure of early-type galaxies should be used together to separate systems that appear to have formed from "wet" versus "dry" mergers.

We attempt to constrain the shape of M31's inner stellar halo by tracing the surface density of blue horizontal branch (BHB) stars at galactocentric distances ranging from 2 kpc to 35 kpc. Our measurements make use of resolved stellar photometry from a section of the Panchromatic Hubble Andromeda Treasury (PHAT) survey, supplemented by several archival Hubble Space Telescope observations. We find that the ratio of BHB to red giant stars is relatively constant outside of 10 kpc, suggesting that the BHB is as reliable a tracer of the halo population as the red giant branch. In the inner halo, we do not expect BHB stars to be produced by the high metallicity bulge and disk, making BHB stars a good candidate to be a reliable tracer of the stellar halo to much smaller galactocentric distances. If we assume a power-law profile r^(-\alpha) for the 2-D projected surface density BHB distribution, we obtain a high-quality fit with a 2-D power-law index of \alpha=2.6^{+0.3}_{-0.2} outside of 3 kpc, which flattens to \alpha<1.2 inside of 3 kpc. This slope is consistent with previous measurements but is anchored to a radial baseline that extends much farther inward. Finally, assuming azimuthal symmetry and a constant mass-to-light ratio, the best-fitting profile yields a total halo stellar mass of 2.1^{+1.7}_{-0.4} x 10^9 M_sun. These properties are comparable with both simulations of stellar halo formation formed by satellite disruption alone, and with simulations that include some in situ formation of halo stars.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. Part of the original charge was to produce an interim design reference mission by mid-2011. That document was delivered to NASA and widely circulated within the astronomical community. In late 2011 the Astrophysics Division augmented its original charge, asking for two design reference missions. The first of these, DRM1, was to be a finalized version of the interim DRM, reducing overall mission costs where possible. The second of these, DRM2, was to identify and eliminate capabilities that overlapped with those of NASA's James Webb Space Telescope (henceforth JWST), ESA's Euclid mission, and the NSF's ground-based Large Synoptic Survey Telescope (henceforth LSST), and again to reduce overall mission cost, while staying faithful to NWNH. This report presents both DRM1 and DRM2.

Hubble Space Telescope images of the galaxy cluster Abell 2261, obtained as part of the Cluster Lensing And Supernova survey with Hubble, show that the brightest galaxy in the cluster, A2261-BCG, has the largest core yet detected in any galaxy. The cusp radius of A2261-BCG is 3.2 kpc, twice as big as the next largest core known, and ~3x bigger than those typically seen in the most luminous BCGs. The morphology of the core in A2261-BCG is also unusual, having a flat or even slightly-depressed interior surface brightness profile, rather than the typical shallow cusp. This implies that the galaxy has a core with constant or even centrally decreasing stellar density. Interpretation of the core as an end product of the "scouring" action of a binary supermassive black hole implies a total black hole mass ~1E+10 M_sun from the extrapolation of most relationships between core structure and black hole mass. The core falls 1-sigma above the cusp-radius versus galaxy luminosity relation. Its large size in real terms, and the extremely large black hole mass required to generate it, raise the possibility that the core has been enlarged by additional processes, such as the ejection of the black holes that originally generated the core. The flat central stellar density profile is consistent with this hypothesis. The core is also displaced by 0.7 kpc from the center of the surrounding envelope, consistent with a local dynamical perturbation of the core.

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

We present time-series photometry of two fields near M32 using archival observations from ACS/WFC onboard HST. One field is centered about 2 arcmin from M32 while the other is located 15 arcmin to the southeast of M31. We identify a total of 1139 RR Lyrae variables of which 821 are ab-type and 318 are c-type. In the field near M32, we find a radial gradient in the density of RR Lyraes relative to the center of M32. This gradient is consistent with the surface brightness profile of M32 suggesting that a significant number of the RR Lyraes in this region belong to M32. This provides further confirmation that M32 contains an ancient stellar population formed around the same time as the oldest population in M31 and the Milky Way. The RR Lyrae stars in M32 exhibit a mean metal abundance of [Fe/H] ~ -1.42 +/- 0.02, which is ~15 times lower than the metal abundance of the overall M32 stellar population. Moreover, the abundance of RR Lyrae stars normalized to the luminosity of M32 in the field analyzed further indicates that the ancient metal-poor population in M32 represents only a very minor component of this galaxy, consistent with the 1% to 4.5% in mass inferred from the CMD analysis of Monachesi et al. In the other field, we find unprecedented evidence for two populations of RR Lyraes in M31 as shown by two distinct sequences among the ab-type variables in the Bailey Diagram. When interpreted in terms of metal abundance, one population exhibits a peak at [Fe/H] ~ -1.3 and the other is at [Fe/H] ~ -1.9. One possible interpretation of this result is that the more metal-rich population represents the dominant M31 halo, while the metal-poorer group could be a disrupted dwarf satellite galaxy orbiting M31. If true, this represents a further indication that the formation of the M31 spheroid has been significantly influenced by the merger and accretion of dwarf galaxy satellites. [abridged]

As part of the Panchromatic Hubble Andromeda Treasury (PHAT) multi-cycle program, we observed a 12' \times 6.5' area of the bulge of M31 with the WFC3/UVIS filters F275W and F336W. From these data we have assembled a sample of \sim4000 UV-bright, old stars, vastly larger than previously available. We use updated Padova stellar evolutionary tracks to classify these hot stars into three classes: Post-AGB stars (P-AGB), Post-Early AGB (PE-AGB) stars and AGB-manqu\'e stars. P-AGB stars are the end result of the asymptotic giant branch (AGB) phase and are expected in a wide range of stellar populations, whereas PE-AGB and AGB-manqu\'e (together referred to as the hot post-horizontal branch; HP-HB) stars are the result of insufficient envelope masses to allow a full AGB phase, and are expected to be particularly prominent at high helium or {\alpha} abundances when the mass loss on the RGB is high. Our data support previous claims that most UV-bright sources in the bulge are likely hot (extreme) horizontal branch stars (EHB) and their progeny. We construct the first radial profiles of these stellar populations, and show that they are highly centrally concentrated, even more so than the integrated UV or optical light. However, we find that this UV-bright population does not dominate the total UV luminosity at any radius, as we are detecting only the progeny of the EHB stars that are the likely source of the UVX. We calculate that only a few percent of MS stars in the central bulge can have gone through the HP-HB phase and that this percentage decreases strongly with distance from the center. We also find that the surface density of hot UV-bright stars has the same radial variation as that of low-mass X-ray binaries. We discuss age, metallicity, and abundance variations as possible explanations for the observed radial variation in the UV-bright population.

The early Universe at redshift z\sim6-11 marks the reionization of the intergalactic medium, following the formation of the first generation of stars. However, those young galaxies at a cosmic age of \lesssim 500 million years (Myr, at z \gtrsim 10) remain largely unexplored as they are at or beyond the sensitivity limits of current large telescopes. Gravitational lensing by galaxy clusters enables the detection of high-redshift galaxies that are fainter than what otherwise could be found in the deepest images of the sky. We report the discovery of an object found in the multi-band observations of the cluster MACS1149+22 that has a high probability of being a gravitationally magnified object from the early universe. The object is firmly detected (12 sigma) in the two reddest bands of HST/WFC3, and not detected below 1.2 {\mu}m, matching the characteristics of z\sim9 objects. We derive a robust photometric redshift of z = 9.6 \pm 0.2, corresponding to a cosmic age of 490 \pm 15Myr (i.e., 3.6% of the age of the Universe). The large number of bands used to derive the redshift estimate make it one of the most accurate estimates ever obtained for such a distant object. The significant magnification by cluster lensing (a factor of \sim15) allows us to analyze the object's ultra-violet and optical luminosity in its rest-frame, thus enabling us to constrain on its stellar mass, star-formation rate and age. If the galaxy is indeed at such a large redshift, then its age is less than 200 Myr (at the 95% confidence level), implying a formation redshift of zf \lesssim 14. The object is the first z>9 candidate that is bright enough for detailed spectroscopic studies with JWST, demonstrating the unique potential of galaxy cluster fields for finding highly magnified, intrinsically faint galaxies at the highest redshifts.

The Panchromatic Hubble Andromeda Treasury (PHAT) is an on-going HST Multicycle Treasury program to image ~1/3 of M31's star forming disk in 6 filters, from the UV to the NIR. The full survey will resolve the galaxy into more than 100 million stars with projected radii from 0-20 kpc over a contiguous 0.5 square degree area in 828 orbits, producing imaging in the F275W and F336W filters with WFC3/UVIS, F475W and F814W with ACS/WFC, and F110W and F160W with WFC3/IR. The resulting wavelength coverage gives excellent constraints on stellar temperature, bolometric luminosity, and extinction for most spectral types. The photometry reaches SNR=4 at F275W=25.1, F336W=24.9, F475W=27.9, F814W=27.1, F110W=25.5, and F160W=24.6 for single pointings in the uncrowded outer disk; however, the optical and NIR data are crowding limited, and the deepest reliable magnitudes are up to 5 magnitudes brighter in the inner bulge. All pointings are dithered and produce Nyquist-sampled images in F475W, F814W, and F160W. We describe the observing strategy, photometry, astrometry, and data products, along with extensive tests of photometric stability, crowding errors, spatially-dependent photometric biases, and telescope pointing control. We report on initial fits to the structure of M31's disk, derived from the density of RGB stars, in a way that is independent of the assumed M/L and is robust to variations in dust extinction. These fits also show that the 10 kpc ring is not just a region of enhanced recent star formation, but is instead a dynamical structure containing a significant overdensity of stars with ages >1 Gyr. (Abridged)

We obtained U_330 and B band images of the M31 nucleus using the High Resolution Camera of the Advanced Camera for Surveys on board the Hubble Space Telescope (HST). The spatial resolution in the U_330-band, 0.03" FWHM, or 0.1 pc at M31, is sufficient to resolve the outskirts of the compact cluster (P3) of UV-bright stars surrounding the M31 black hole. The center of the cluster is marked by an extended source that is both brighter and redder than the other point sources within P3; it is likely to be a blend of several bright stars. We hypothesize that it marks the location of the M31 black hole. Both stellar photometry and a surface brightness fluctuation analysis, show that the P3 stellar population is consistent with early-type main sequence stars formed in a ~100 - ~200 Myr old starburst population. Evolutionary tracks of post early asymptotic giant-branch stars, associated with late-stage evolution of an old population, also traverse the U and U-B domain occupied by the P3 stars; but we argue that only a few stars could be accounted for that way. PEAGB evolution is very rapid, and there is no progenitor population of red giants associated with P3. The result that P3 comprises young stars is consistent with inferences from earlier HST observations of the integrated light of the cluster. Like the Milky Way, M31 harbors a black hole closely surrounded by apparently young stars.

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

We use deep HST ACS/HRC observations of a field within M32 (F1) and an M31 background field (F2) to determine the star formation history (SFH) of M32 from its resolved stellar population. We find that 2-5Gyr old stars contribute \som40%+/- 17% of M32's mass, while 55%+/-21% of M32's mass comes from stars older than 5 Gyr. The mass-weighted mean age and metallicity of M32 at F1 are <Age>=6.8+/-1.5 Gyr and <[M/H]>=-0.01+/-0.08 dex. The SFH additionally indicates the presence of young (<2 Gyr old), metal-poor ([M/H]\sim-0.7) stars, suggesting that blue straggler stars contribute ~2% of the mass at F1; the remaining \sim3% of the mass is in young metal-rich stars. Line-strength indices computed from the SFH imply a light-weighted mean age and metallicity of 4.9 Gyr and [M/H] = -0.12 dex, and single-stellar-population-equivalent parameters of 2.9+/-0.2 Gyr and [M/H]=0.02+/-0.01 dex at F1 (~2.7 re). This contradicts spectroscopic studies that show a steep age gradient from M32's center to 1re. The inferred SFH of the M31 background field F2 reveals that the majority of its stars are old, with \sim95% of its mass already acquired 5-14 Gyr ago. It is composed of two dominant populations; \sim30%+/-7.5% of its mass is in a 5-8 Gyr old population, and \sim65%+/-9% of the mass is in a 8-14 Gyr old population. The mass-weighted mean age and metallicity of F2 are <Age>=9.2+/-1.2 Gyr and <[M/H]>=-0.10+/-0.10 dex, respectively. Our results suggest that the inner disk and spheroid populations of M31 are indistinguishable from those of the outer disk and spheroid. Assuming the mean age of M31's disk at F2 (\sim1 disk scale length) to be 5-9 Gyr, our results agree with an inside-out disk formation scenario for M31's disk.

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.

In December 2010, NASA created a Science Definition Team (SDT) for WFIRST, the Wide Field Infra-Red Survey Telescope, recommended by the Astro 2010 Decadal Survey as the highest priority for a large space mission. The SDT was chartered to work with the WFIRST Project Office at GSFC and the Program Office at JPL to produce a Design Reference Mission (DRM) for WFIRST. This paper describes an Interim DRM. The DRM will be completed in 2012.

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

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

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

We present the deepest optical color-magnitude diagram (CMD) to date of the local elliptical galaxy M32. We have obtained F435W and F555W photometry based on HST ACS/HRC images for a region 110" from the center of M32 and a background field about 320" away from M32 center. Due to the high resolution of our Nyquist-sampled images, the small photometric errors, and the depth of our data we obtain the most detailed resolved photometric study of M32 yet. Deconvolution of HST images proves to be superior than other standard methods to derive stellar photometry on extremely crowded HST images. The location of the strong red clump in the CMD suggests a mean age between 8 and 10 Gyr for [Fe/H] = -0.2 in M32. We detect for the first time a red giant branch bump and an asymptotic giant branch bump in M32 which indicate that the mean age of M32's dominant population at ~2' from its center is between 5 and 10 Gyr. We see evidence of an intermediate-age population in M32 mainly due to the presence of bright asymptotic giant branch stars. Our detection of a blue component of stars (blue plume) may indicate for the first time the presence of a young stellar population, with ages of the order of 0.5 Gyr, in our M32 field. However, it is likely that the brighter stars of this blue plume belong to the disk of M31 rather than to M32. The fainter stars populating the blue plume indicate the presence of stars not younger than 1 Gyr and/or blue straggler stars in M32. M32's dominant population of 8--10 Gyr implies a formation redshift of 1 < z_f < 2, precisely when observations of the specific star formation rates and models of "downsizing" imply galaxies of M32's mass ought to be forming their stars. Our CMD therefore provides a "ground-truth" of downsizing scenarios at z=0. Our background field data represent the deepest optical observations yet of the inner disk and bulge of M31. [Abridged]

The Karhunen-Loeve (KL) transform can compactly represent the information contained in large, complex datasets, cleanly eliminating noise from the data and identifying elements of the dataset with extreme or inconsistent characteristics. We develop techniques to apply the KL transform to the 4000-5700A region of 9,800 QSO spectra with z < 0.619 from the SDSS archive. Up to 200 eigenspectra are needed to fully reconstruct the spectra in this sample to the limit of their signal/noise. We propose a simple formula for selecting the optimum number of eigenspectra to use to reconstruct any given spectrum, based on the signal/noise of the spectrum, but validated by formal cross-validation tests. We show that such reconstructions can boost the effective signal/noise of the observations by a factor of 6 as well as fill in gaps in the data. The improved signal/noise of the resulting set will allow for better measurement and analysis of these spectra. The distribution of the QSO spectra within the eigenspace identifies regions of enhanced density of interesting subclasses, such as Narrow Line Seyfert 1s (NLS1s). The weightings, as well as the inability of the eigenspectra to fit some of the objects, also identifies "outliers," which may be objects that are not valid members of the sample or objects with rare or unique properties. We identify 48 spectra from the sample that show no broad emission lines, 21 objects with unusual [O III] emission line properties, and 9 objects with peculiar H-beta emission line profiles. We also use this technique to identify a binary supermassive black hole candidate. We provide the eigenspectra and the reconstructed spectra of the QSO sample.

We observed two fields near M32 with the Advanced Camera for Surveys/High Resolution Channel (ACS/HRC) on board the Hubble Space Telescope (HST). The main field, F1, is 1.8 arcmin from the center of M32; the second field, F2, constrains the M31 background, and is 5.4 arcmin distant. Each field was observed for 16-orbits in each of the F435W (narrow B) and F555W (narrow V) filters. The duration of the observations allowed RR Lyrae stars and other short-period variables to be detected. A population of RR Lyrae stars determined to belong to M32 would prove the existence of an ancient population in that galaxy, a subject of some debate. We detected 17 RR Lyrae variables in F1 and 14 in F2. A 1-sigma upper limit of 6 RR Lyrae variables belonging to M32 is inferred from these two fields alone. Use of our two ACS/WFC parallel fields provides better constraints on the M31 background, however, and implies that $7_{-3}^{+4}$ (68 % confidence interval) RR Lyrae variables in F1 belong to M32. We have therefore found evidence for an ancient population in M32. It seems to be nearly indistinguishable from the ancient population of M31. The RR Lyrae stars in the F1 and F2 fields have indistinguishable mean V-band magnitudes, mean periods, distributions in the Bailey diagram and ratios of RRc to RR(tot) types. However, the color distributions in the two fields are different, with a population of red RRab variables in F1 not seen in F2. We suggest that these might be identified with the detected M32 RR Lyrae population, but the small number of stars rules out a definitive claim.

Observations suggest that effective radii of high-z massive spheroids are as much as a factor ~6 smaller than low-z galaxies of comparable mass. Given the apparent absence of low-z counterparts, this has often been interpreted as indicating that the high density, compact red galaxies must be 'puffed up' by some mechanism. We compare the ensemble of high-z observations with large samples of well-observed low-z ellipticals. At the same physical radii, the stellar surface mass densities of low and high-z systems are comparable. Moreover, the abundance of high surface density material at low redshift is comparable to or larger than that observed at z>1-2, consistent with the continuous buildup of spheroids over this time. The entire population of compact, high-z red galaxies may be the progenitors of the high-density cores of present-day ellipticals, with no need for a decrease in stellar density from z=2 to z=0. The primary difference between low and high-z systems is thus the observed low-density material at large radii in low-z spheroids (rather than the high-density material in high-z spheroids). Such low-density material may either (1) assemble at z<2 or (2) be present, but not yet detected, at z>2. Mock observations of low-z massive systems show that the high-z observations do not yet probe sufficiently low surface brightness material to detect the low surface density 'wings' (if present). Thus, if the high-z galaxies resemble the most massive systems today, their inferred effective radii could be under-estimated by factors ~2-4. This difference arises because massive systems at low redshift are not well-fit by single Sersic profiles. We discuss implications of our results for physical models of galaxy evolution.

We observed two fields near M32 with the ACS/HRC on board the Hubble Space Telescope, located at distances of about 1.8' and 5.4' (hereafter F1 and F2, respectively) from the center of M32. To obtain a very detailed and deep color-magnitude diagram (CMD) and to look for short period variability, we obtained time-series imaging of each field in 32-orbit-long exposures using the F435W (B) and F555W (V) filters, spanning a temporal range of 2 days per filter. We focus on our detection of variability on RR Lyrae variable stars, which represents the only way to obtain information about the presence of a very old population (larger than 10 Gyr) in M32 from optical data. Here we present results obtained from the detection of 31 RR Lyrae in these fields: 17 in F1 and 14 in F2.

We present HST WFPC2/PC images and KPNO 4-m longslit spectroscopy of the QSO SDSS J153636.22+044127.0, which we advanced as a candidate binary supermassive black hole. The images reveal a close companion coincident with the radio source identified by Wrobel & Laor (2009). It appears to be consistent with a M_g ~ -21.4 elliptical galaxy, if it is at the QSO redshift. The spectroscopy, however, shows no spatial offset of the red or blue Balmer line subcomponents. The companion is thus not the source of either the red or blue broad line systems; SDSS J153636.22+044127.0 cannot be explained as a chance superposition of objects, or as an ejected black hole. Over the Delta T=0.75 yr difference between the rest frame epochs of the present and SDSS spectroscopy, we find no velocity shift to within 40 km/s, nor any amplitude change in either broad line system. The lack of a shift can be admitted under the binary hypothesis if the implied radial velocity is a larger component of the full orbital velocity than was assumed in our earlier work. A strong test of the binary hypothesis requires yet longer temporal baselines. The lack of amplitude variations is unusual for the alternative explanation of this object as a "double-peaked" emitter; we further argue that SDSS J153636.22+044127.0 has unique spectral features that have no obvious analogue with other members of this class.

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 present Hubble Space Telescope observations taken with the Advanced Camera for Surveys Wide Field Channel of two fields near M32 - between four and six kpc from the center of M31. The data cover a time baseline sufficient for the identification and characterization of 681 RR Lyrae variables of which 555 are ab-type and 126 are c-type. The mean magnitude of these stars is <V>=25.29 +/- 0.05 where the uncertainty combines both the random and systematic errors. The location of the stars in the Bailey Diagram and the ratio of c-type RR Lyraes to all types are both closer to RR Lyraes in Oosterhoff type I globular clusters in the Milky Way as compared with Oosterhoff II clusters. The mean periods of the ab-type and c-type RR Lyraes are <P(ab)>=0.557 +/- 0.003 and <P(c)>=0.327 +/- 0.003, respectively, where the uncertainties in each case represent the standard error of the mean. When the periods and amplitudes of the ab-type RR Lyraes in our sample are interpreted in terms of metallicity, we find the metallicity distribution function to be indistinguishable from a Gaussian with a peak at <[Fe/H]>=-1.50 +/- 0.02, where the quoted uncertainty is the standard error of the mean. Using a relation between RR Lyrae luminosity and metallicity along with a reddening of E(B-V) = 0.08 +/- 0.03, we find a distance modulus of (m-M)o=24.46 +/- 0.11 for M31. We examine the radial metallicity gradient in the environs of M31 using published values for the bulge and halo of M31 as well as the abundances of its dwarf spheroidal companions and globular clusters. In this context, we conclude that the RR Lyraes in our two fields are more likely to be halo objects rather than associated with the bulge or disk of M31, in spite of the fact that they are located at 4-6 kpc in projected distance from the center.

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 investigate how 'extra' central light in the surface brightness profiles of cusp ellipticals relates to the profiles of ellipticals with cores. Cusp elliptical envelopes are formed by violent relaxation in mergers acting on stars in progenitor disks, while their centers are structured by dissipational starbursts. Core ellipticals are formed by subsequent merging of (now gas-poor) cusp ellipticals, with the fossil starburst components combining to preserve a compact component in the remnant (although the 'transition' is smoothed). Comparing hydrodynamical simulations and observed profiles, we show how to observationally isolate the relic starburst components in core ellipticals. We demonstrate that these survive re-mergers and reliably trace the dissipation in the initial gas-rich merger(s). The typical degree of dissipation is a strong function of stellar mass, tracing observed disk gas fractions. We find a correlation between dissipation and effective radius: systems with more dissipation are more compact. The survival of this component and scattering of stars into the envelope naturally explain high-Sersic index profiles characteristic of massive core ellipticals. This is also closely related to the kinematics and isophotal shapes: only systems with matched starburst components from their profile fits also reproduce the observed kinematics of boxy/core ellipticals. We show that it is critical to adopt physically motivated profiles when attempting to quantify how much mass has been 'scoured' or scattered out of the inner regions by binary black holes. Estimates of scoured mass ignoring multi-component structure can be strongly biased, potentially explaining observed systems with large inferred core masses in apparent conflict with core-scouring models.

We identify SDSS J153636.22+044127.0, a QSO discovered in the Sloan Digital Sky Survey, as a promising candidate for a binary black hole system. This QSO has two broad-line emission systems separated by 3500 km/sec. The redder system at z=0.3889 also has a typical set of narrow forbidden lines. The bluer system (z=0.3727) shows only broad Balmer lines and UV Fe II emission, making it highly unusual in its lack of narrow lines. A third system, which includes only unresolved absorption lines, is seen at a redshift, z=0.3878, intermediate between the two emission-line systems. While the observational signatures of binary nuclear black holes remain unclear, J1536+0441 is unique among all QSOs known in having two broad-line regions, indicative of two separate black holes presently accreting gas. The interpretation of this as a bound binary system of two black holes having masses of 10^8.9 and 10^7.3 solar masses, yields a separation of ~ 0.1 parsec and an orbital period of ~100 years. The separation implies that the two black holes are orbiting within a single narrow-line region, consistent with the characteristics of the spectrum. This object was identified as an extreme outlier of a Karhunen-Loeve Transform of 17,500 z < 0.7 QSO spectra from the SDSS. The probability of the spectrum resulting from a chance superposition of two QSOs with similar redshifts is estimated at 2X10^-7, leading to the expectation of 0.003 such objects in the sample studied; however, even in this case, the spectrum of the lower redshift QSO remains highly unusual.

We present an analysis of a 50 orbit HST ACS observation of the M87 globular cluster system. We use the extraordinary depth of this dataset to test whether the colors and magnitudes show evidence for a mass-metallicity relation in globular cluster populations. We find only a weak or absent relation between the colors and magnitudes of the metal poor subpopulation of globular clusters. The weakness or absence of a color-magnitude relation is established over a wide range in luminosity from $M_V=-11$ to $M_V=-6$, encompassing most of the M87 globular clusters. The constancy of the colors of the metal-poor subpopulation seen in our 50 orbit observation is in contrast to suggestions from single orbit ACS data that the metal-poor globular clusters in M87 and several other galaxies show a "blue tilt." The formal best fit for the mass-metallicity relation for the metal-poor subpopulation in our much deeper data is $Z\propto M^{0.08\pm0.05}$. Our analysis of these data also shows a possible small "red tilt" in the metal-rich globular cluster subpopulation. While either of these small tilts may be real, they may also illustrate the limit to which mass-metallicity relations can be determined, even in such extraordinarily deep data. We specifically test for a wide range of systematic effects and find that while small tilts cannot be confirmed or rejected, the data place a strong upper limit to any tilt of $|0.20|\pm0.05$. This upper limit is much smaller than some earlier claims from single orbit data, and strongly limits self-enrichment within globular clusters. This mass-metallicity relation for globular clusters is also shallower than the relation for galaxies, suggesting that the formation mechanisms for these two types of objects are different.

We study the origin and properties of 'extra' or 'excess' central light in the surface brightness profiles of cusp or power-law ellipticals. Dissipational mergers give rise to two-component profiles: an outer profile established by violent relaxation acting on stars present in the progenitors prior to the final merger, and an inner stellar population comprising the extra light, formed in a compact starburst. Combining a large set of hydrodynamical simulations with data that span a broad range of profiles and masses, we show that this picture is borne out -- cusp ellipticals are indeed 'extra light' ellipticals -- and examine how the properties of this component scale with global galaxy properties. We show how to robustly separate the 'extra' light, and demonstrate that observed cusps are reliable tracers of the degree of dissipation in the spheroid-forming merger. We show that the typical degree of dissipation is a strong function of stellar mass, tracing observed disk gas fractions at each mass. We demonstrate a correlation between extra light content and effective radius at fixed mass: systems with more dissipation are more compact. The outer shape of the light profile does not depend on mass, with a mean outer Sersic index ~2.5. We explore how this relates to shapes, kinematics, and stellar population gradients. Simulations with the gas content needed to match observed profiles also reproduce observed age, metallicity, and color gradients, and we show how these can be used as tracers of the degree of dissipation in spheroid formation.

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

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

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

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

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

Destiny is a simple, direct, low cost mission to determine the properties of dark energy by obtaining a cosmologically deep supernova (SN) type Ia Hubble diagram. Operated at L2, its science instrument is a 1.65m space telescope, featuring a grism-fed near-infrared (NIR) (0.85-1.7micron) survey camera/spectrometer with a 0.12 square degree field of view. During its two-year primary mission, Destiny will detect, observe, and characterize ~3000 SN Ia events over the redshift interval 0.4<z<1.7 within a 3 square degree survey area. In conjunction with ongoing ground-based SN Ia surveys for z<0.8, Destiny mission data will be used to construct a high-precision Hubble diagram and thereby constrain the dark energy equation of state from a time when it was strongly matter-dominated to the present when dark energy dominates. The grism-images simultaneously provide broad-band photometry, redshifts, and SN classification, as well as time-resolved diagnostic data for investigating additional SN luminosity diagnostics. Destiny will be used in its third year as a high resolution, wide-field imager to conduct a multicolor NIR weak lensing (WL) survey covering 1000 square degrees. The large-scale mass power spectrum derived from weak lensing distortions of field galaxies as a function of redshift will provide independent and complementary constraints on the dark energy equation of state. The combination of SN and WL is much more powerful than either technique on its own. Used together, these surveys will have more than an order of magnitude greater sensitivity than will be provided by ongoing ground-based projects. The dark energy parameters, w_0 and w_a, will be measured to a precision of 0.05 and 0.2 respectively.

We present the luminosity function to very faint magnitudes for the globular clusters in M87, based on a 30 orbit \textit{Hubble Space Telescope (HST)} WFPC2 imaging program. The very deep images and corresponding improved false source rejection allow us to probe the mass function further beyond the turnover than has been done before. We compare our luminosity function to those that have been observed in the past, and confirm the similarity of the turnover luminosity between M87 and the Milky Way. We also find with high statistical significance that the M87 luminosity function is broader than that of the Milky Way. We discuss how determining the mass function of the cluster system to low masses can constrain theoretical models of the dynamical evolution of globular cluster systems. Our mass function is consistent with the dependence of mass loss on the initial cluster mass given by classical evaporation, and somewhat inconsistent with newer proposals that have a shallower mass dependence. In addition, the rate of mass loss is consistent with standard evaporation models, and not with the much higher rates proposed by some recent studies of very young cluster systems. We also find that the mass-size relation has very little slope, indicating that there is almost no increase in the size of a cluster with increasing mass.

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

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

We present HST/WFPC2 observations of the five bluest E+A galaxies (z~0.1) in the Zabludoff et al. sample to study whether their detailed morphologies are consistent with late-to-early type evolution and to determine what drives that evolution. The morphologies of four galaxies are disturbed, indicating that a galaxy-galaxy merger is at least one mechanism that leads to the E+A phase. Two-dimensional image fitting shows that the E+As are generally bulge-dominated systems, even though at least two E+As may have underlying disks. In the Fundamental Plane, E+As stand apart from the E/S0s mainly due to their high effective surface brightness. Fading of the young stellar population and the corresponding increase in their effective radii will cause these galaxies to migrate toward the locus of E/S0s. E+As have profiles qualitatively like those of normal power-law early-type galaxies, but have higher surface brightnesses. This result provides the first direct evidence supporting the hypothesis that power-law ellipticals form via gas-rich mergers. In total, at least four E+As are morphologically consistent with early-type galaxies. We detect compact sources, possibly young star clusters, associated with the galaxies. These sources are much brighter (M_R ~ -13) than Galactic globular clusters, have luminosities consistent with the brightest clusters in nearby starburst galaxies, and have blue colors consistent with the ages estimated from the E+A galaxy spectra (several 10^8 yr). Further study of such young star cluster candidates might provide the elusive chronometer needed to break the age/burst-strength degeneracy for these post-merger galaxies.

Only one certain classical nova eruption has ever been detected inside a globular cluster - nova 1860 A.D. (T Sco) in M80. During a survey of M87 we have detected an erupting star coincident (to within 0.08 pixels) with a globular cluster of that giant elliptical galaxy. We are able to discount variables in the foreground or background of M87. The light curve and color of the erupting star match those expected for a nova at the distance of M87. The chance superposition of an M87 field nova on the globular cluster is very unlikely but cannot be completely ruled out.Our detection hints at a globular cluster nova frequency $f \sim .004$ novae/cluster/year, much higher than previous observations have suggested.

The position of the giant elliptical galaxy M87 at the center of the Virgo Cluster means that the inferred column density of dark matter associated with both the cluster halo and the galaxy halo is quite large. This system is thus an important laboratory for studying massive dark objects in elliptical galaxies and galaxy clusters by gravitational microlensing, strongly complementing the studies of spiral galaxy halos performed in the Local Group. We have performed a microlensing survey of M87 with the WFPC2 instrument on the Hubble Space Telescope. Over a period of thirty days, with images taken once daily, we discover seven variable sources. Four are variable stars of some sort, two are consistent with classical novae, and one exhibits an excellent microlensing lightcurve, though with a very blue color implying the somewhat disfavored possibility of a horizontal branch source being lensed. Based on sensitivity calculations from artificial stars and from artificial lightcurves, we estimate the expected microlensing rate. We find that the detection of one event is consistent with a dark halo with a 20% contribution of microlensing objects for both M87 and the Virgo Cluster, similar to the value found from observations in the Local Group. Further work is required to test the hypothesized microlensing component to the cluster.

We have measured the clustering of z<0.9 red galaxies and constrained models of the evolution of large-scale structure using the initial 1.2 sq. degree data release of the NOAO Deep Wide-Field Survey (NDWFS). The area and BwRI passbands of the NDWFS allow samples of >1000 galaxies to be selected as a function of spectral type, absolute magnitude, and photometric redshift. Spectral synthesis models can be used to predict the colors and luminosities of a galaxy population as a function of redshift. We have used PEGASE2 models, with exponentially declining star formation rates, to estimate the observed colors and luminosity evolution of galaxies and to connect, as an evolutionary sequence, related populations of galaxies at different redshifts. A red galaxy sample, with present-day rest-frame Vega colors of Bw-R>1.44, was chosen to allow comparisons with the 2dF Galaxy Redshift Survey and Sloan Digital Sky Survey. We find the spatial clustering of red galaxies to be a strong function of luminosity, with r0 increasing from 4.4+/-0.4 Mpc/h at M_R=-20 to 11.2+/-1.0 Mpc/h at M_R=-22. Clustering evolution measurements using samples where the rest-frame selection criteria vary with redshift, including all deep single-band magnitude limited samples, are biased due to the correlation of clustering with rest-frame color and luminosity. The clustering of M_R=-21, Bw-R>1.44 galaxies exhibits no significant evolution over the redshift range observed with r0= 6.3+/-0.5 Mpc/h in comoving coordinates. This is consistent with recent LCDM models where the bias of L* galaxies undergoes rapid evolution and r0 evolves very slowly at z<2.

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 present a new sample of 35 quasars selected from the FIRST radio survey and the Deeprange I-band survey (Postman et al. 1998, 2002). A comparison with the FIRST Bright Quasar survey samples reveals that this I-band selected sample is redder by 0.25-0.5 mag in B-R, and that the color difference is not explained by the higher mean redshift of this sample but must be intrinsic. Our small sample contains five quasars with unusually red colors, including three that appear very heavily reddened. Our data are fitted well with normal blue quasar spectra attenuated by more than 2.5 magnitudes of extinction in the I-band. The red quasars are only seen at low redshifts (z<1.3). Even with a magnitude limit I<20.5, our survey is deep enough to detect only the most luminous of these red quasars at z~1; similar objects at higher redshifts would fall below our in I-band limit. Indeed, the five most luminous objects (using dereddened magnitudes) with z<1.3 are all red. Our data strongly support the hypothesis that radio quasars are dominated by a previously undetected population of red, heavily obscured objects. Unless highly reddened quasars are preferentially also highly luminous, there must be an even larger, as yet undiscovered, population of red quasars at lower luminosity. We are likely to be finding only the most luminous tip of the red quasar iceberg. Five of our six z<1 quasars are associated with Deeprange cluster candidates of similar estimated redshifts. This association is very unlikely to be the result of chance. It has some surprising implications, including the possibility that up to half of the Deeprange clusters at z~1 have associated quasars.

We used the HST WFPC2 to obtain I-band images of the centers of 81 brightest cluster galaxies (BCGs), drawn from a volume-limited sample of nearby BCGs. The images show a rich variety of morphological features, including multiple or double nuclei, dust, stellar disks, point source nuclei, and central surface brightness depressions. High resolution surface brightness profiles could be inferred for 60 galaxies. Of those, 88% have well-resolved cores. Twelve percent of the BCG sample lacks a well-resolved core; all but one of these BCGs have ``power-law'' profiles. Some of these galaxies have higher luminosities than any power-law galaxy identified by Faber et al. (1997), and have physical upper limits on the break radius well below the values observed for core galaxies of the same luminosity. These results support the idea that the central structure of early-type galaxies is bimodal in its physical properties, but also suggest that there exist high luminosity galaxies with power-law profiles (or unusually small cores). The BCGs in the latter category tend to fall at the low end of the BCG luminosity function and tend to have low values of the quantity alpha (the logarithmic slope of the metric luminosity as a function of radius, at 10 kpc). Since theoretical calculations have shown that the luminosities and alpha values of BCGs grow with time as a result of accretion, this suggests a scenario in which elliptical galaxies evolve from power-law profiles to core profiles through accretion and merging. This is consistent with theoretical scenarios that invoke the formation of massive black hole binaries during merger events (Abridged).

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.

The Galactic Exoplanet Survey Telescope (GEST) will observe a 2 square degree field in the Galactic bulge to search for extra-solar planets using a gravitational lensing technique. This gravitational lensing technique is the only method employing currently available technology that can detect Earth-mass planets at high signal-to-noise, and can measure the frequency of terrestrial planets as a function of Galactic position. GEST's sensitivity extends down to the mass of Mars, and it can detect hundreds of terrestrial planets with semi-major axes ranging from 0.7 AU to infinity. GEST will be the first truly comprehensive survey of the Galaxy for planets like those in our own Solar System.

Application of deconvolution algorithms to astronomical images is often limited by variations in PSF structure over the domain of the images. One major difficulty is that Fourier methods can no longer be used for fast convolutions over the entire images. However, if the PSF is modeled as a sum of orthogonal functions that are individually constant in form over the images, but whose relative amplitudes encode the PSF spatial variability, then separation of variables again allows global image operations to be used. This approach is readily adapted to the Lucy-Richardson deconvolution algorithm. Use of the Karhunen-Loeve transform allows for a particularly compact orthogonal expansion of the PSF. These techniques are demonstrated on the deconvolution of Gemini/Hokupa'a adaptive optics images of the galactic center.

We have conducted an automated search for galaxy clusters within a contiguous 16 square degree I-band survey in the north Galactic hemisphere. A matched filter detection algorithm identifies 444 cluster candidates in the range 0.2 <= z <= 1.2. The full catalog is presented along with the results from a follow-up spectroscopic survey. The estimated redshift distribution of the cluster candidates is consistent with a constant comoving density over the range 0.2 <= z_est <= 0.8. A decline in the cluster space density by more than a factor of 3 over this redshift range is rejected at >99.9% confidence level. We find that the space density of Lambda_CL >= 40 clusters in our survey lies in the range (1.6 - 1.8) x 10^{-5} h_{75}^{3} Mpc^{-3}, ~1.5 times higher than the local distribution of comparably rich Abell RC >= 0 clusters. The Lambda_CL distribution is consistent with a power-law. The discrepancy between the space density of Abell clusters and the clusters in this survey declines quickly as Lambda_CL increases, suggesting that the difference at lower richness is due to significant incompleteness in the Abell catalog. A percolation analysis reveals that 10 - 20% of the spectroscopically confirmed distant clusters are linked into superclusters at overdensities between 10 < {\delta\rho \over \rho} < 50, similar to what is seen in the local cluster distribution. This suggests that there has been little evolution of the cluster-cluster correlation length for z <= 0.5.

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

Observations of nearby galaxies reveal a strong correlation between the mass of the central dark object M and the velocity dispersion sigma of the host galaxy, of the form log(M/M_sun) = a + b*log(sigma/sigma_0); however, published estimates of the slope b span a wide range (3.75 to 5.3). Merritt & Ferrarese have argued that low slopes (<4) arise because of neglect of random measurement errors in the dispersions and an incorrect choice for the dispersion of the Milky Way Galaxy. We show that these explanations account for at most a small part of the slope range. Instead, the range of slopes arises mostly because of systematic differences in the velocity dispersions used by different groups for the same galaxies. The origin of these differences remains unclear, but we suggest that one significant component of the difference results from Ferrarese & Merritt's extrapolation of central velocity dispersions to r_e/8 (r_e is the effective radius) using an empirical formula. Another component may arise from dispersion-dependent systematic errors in the measurements. A new determination of the slope using 31 galaxies yields b=4.02 +/- 0.32, a=8.13 +/- 0.06, for sigma_0=200 km/s. The M-sigma relation has an intrinsic dispersion in log M that is no larger than 0.3 dex. In an Appendix, we present a simple model for the velocity-dispersion profile of the Galactic bulge.

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

We measured infrared surface brightness fluctuation (SBF) distances to an isotropically-distributed sample of 16 distant galaxies with redshifts reaching 10,000 km/s using the near-IR camera and multi-object spectrometer (NICMOS) on the Hubble Space Telescope (HST). The excellent spatial resolution, very low background, and brightness of the IR fluctuations yielded the most distant SBF measurements to date. Twelve nearby galaxies were also observed and used to calibrate the F160W (1.6 micron) SBF distance scale. Of these, three have Cepheid variable star distances measured with HST and eleven have optical I-band SBF distance measurements. A distance modulus of 18.5 mag to the Large Magellanic Cloud was adopted for this calibration. We present the F160W SBF Hubble diagram and find a Hubble constant Ho=76 +/- 1.3 (1-sigma statistical) +/- 6 (systematic) km/s/Mpc. This result is insensitive to the velocity model used to correct for local bulk motions. Restricting the fit to the six most distant galaxies yields the smallest value of Ho=72 +/- 2.3 km/s/Mpc consistent with the data. This 6% decrease in the Hubble constant is consistent with the hypothesis that the Local Group inhabits an under-dense region of the universe, but is also consistent with the best-fit value of Ho=76 km/s/Mpc at the 1.5-sigma level.

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.

Constraints on the validity of the hierarchical gravitational instability theory and the evolution of biasing are presented based upon measurements of higher order clustering statistics in the Deeprange Survey, a catalog of $\sim710,000$ galaxies with $I_{AB} \le 24$ derived from a KPNO 4m CCD imaging survey of a contiguous $4^{\circ} \times 4^{\circ}$ region. We compute the 3-point and 4-point angular correlation functions using a direct estimation for the former and the counts-in-cells technique for both. The skewness $s_3$ decreases by a factor of $\simeq 3-4$ as galaxy magnitude increases over the range $17 \le I \le 22.5$ ($0.1 \lesssim z \lesssim 0.8$). This decrease is consistent with a small {\it increase} of the bias with increasing redshift, but not by more than a factor of 2 for the highest redshifts probed. Our results are strongly inconsistent, at about the $3.5-4 \sigma$ level, with typical cosmic string models in which the initial perturbations follow a non-Gaussian distribution - such models generally predict an opposite trend in the degree of bias as a function of redshift. We also find that the scaling relation between the 3-point and 4-point correlation functions remains approximately invariant over the above magnitude range. The simplest model that is consistent with these constraints is a universe in which an initially Gaussian perturbation spectrum evolves under the influence of gravity combined with a low level of bias between the matter and the galaxies that decreases slightly from $z \sim 0.8$ to the current epoch.

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

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

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

An undersampled point spread function may interact with the microstructure of a solid-state detector such that the total flux detected can depend sensitively on where the PSF center falls within a pixel. Such intra-pixel sensitivity variations will not be corrected by flat field calibration and may limit the accuracy of stellar photometry conducted with undersampled images, as are typical for Hubble Space Telescope observations. The total flux in a stellar image can vary by up to 0.03 mag in F555W WFC images depending on how it is sampled, for example. For NIC3, these variations are especially strong, up to 0.39 mag, strongly limiting its use for stellar photometry. Intra-pixel sensitivity variations can be corrected for, however, by constructing a well-sampled PSF from a dithered data set. The reconstructed PSF is the convolution of the optical PSF with the pixel response. It can be evaluated at any desired fractional pixel location to generate a table of photometric corrections as a function of relative PSF centroid. A caveat is that the centroid of an undersampled PSF can also be affected by the pixel response function, thus sophisticated centroiding methods, such as cross-correlating the observed PSF with its fully-sampled counterpart, are required to derive the proper photometric correction.

Undersampled images, such as those produced by the HST WFPC-2, misrepresent fine-scale structure intrinsic to the astronomical sources being imaged. Analyzing such images is difficult on scales close to their resolution limits and may produce erroneous results. A set of ``dithered'' images of an astronomical source generally contains more information about its structure than any single undersampled image, however, and may permit reconstruction of a ``superimage'' with Nyquist sampling. I present a tutorial on a method of image reconstruction that builds a superimage from a complex linear combination of the Fourier transforms of a set of undersampled dithered images. This method works by algebraically eliminating the high order satellites in the periodic transforms of the aliased images. The reconstructed image is an exact representation of the data-set with no loss of resolution at the Nyquist scale. The algorithm is directly derived from the theoretical properties of aliased images and involves no arbitrary parameters, requiring only that the dithers are purely translational and constant in pixel-space over the domain of the object of interest. I show examples of its application to WFC and PC images. I argue for its use when the best recovery of point sources or morphological information at the HST diffraction limit is of interest.

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.

WFPC-2 images are used to study the central structure of M31, M32, and M33. The dimmer peak, P2, of the M31 double nucleus is centered on the bulge to 0.1", implying that it is the dynamical center of M31. P2 contains a compact source discovered by King et al. (1995) at 1700 A. This source is resolved, with r_{1/2} approx0.2 pc. It dominates the nucleus at 3000 A, and is consistent with late B-early A stars. This probable cluster may consist of young stars and be an older version of the cluster of hot stars at the center of the Milky Way, or it may consist of heavier stars built up from collisions in a possible cold disk of stars orbiting P2. In M32, the central cusp rises into the HST limit with gamma approx0.5, and the central density rho_0>10^7M_sol pc^-3. The V-I and U-V color profiles are flat, and there is no sign of an inner disk, dust, or any other structure. This total lack of features seems at variance with a nominal stellar collision time of 2 X 10^10 yr, which implies that a significant fraction of the light in the central pixel should come from blue stragglers. InM33, the nucleus has an extremely steep gamma=1.49 power-law profile for 0.05"<r<0.2" that becomes shallower as the HST resolution limit is approached. The profile for r<0.04" has either a gamma approx 0.8 cusp or a small core with r_c ~<0.13 pc. The central density is rho_0 > 2 10^6M_sol pc^-3, and the implied relaxation time is only ~3 X 10^6 yr, indicating that the nucleus is highly relaxed. The accompanying short collision time of 7 X 10^9 yr predicts a central blue straggler component quantitatively consistent with the strong V-I and B-R color gradients seen with HST and from the ground.

We present constraints on the evolution of large-scale structure from a catalog of 710,000 galaxies with I_AB <= 24 derived from a KPNO 4m CCD imaging survey of a contiguous 4 deg x 4 deg region. The advantage of using large contiguous surveys for measuring clustering properties on even modest angular scales is substantial: the effects of cosmic scatter are strongly suppressed. We provide highly accurate measurements of the two-point angular correlation function, w(theta), as a function of magnitude on scales up to 1.5 degrees. The amplitude of w(theta) declines by a factor of ~10 over the range 16 <= I <= 20 but only by a factor of 2 - 3 over the range 20 < I <= 23. For a redshift dependence of the spatial correlation function, xi(r), parameterized as xi(r,z)=(r/r_o)^(-gamma)(1 + z)^(-[3+epsilon]), we find r_o=5.2 +/- 0.4 Mpc/h, and epsilon >= 0 for I <= 20. This is in good agreement with the results from local redshift surveys. At I > 20, our best fit values shift towards lower r_o and more negative epsilon. A strong covariance between r_o and epsilon prevent us from rejecting epsilon > 0 even at faint magnitudes but if epsilon > 1, we strongly reject r_o <= 4/h Mpc (co-moving). The above expression for xi(r,z) and our data give a correlation length of r_o(z=0.5) approx 3.0 +/- 0.4 Mpc/h, about a factor of 2 larger than the correlation length at z = 0.5 derived from the Canada--France Redshift Survey (CFRS). The small volume sampled by the CFRS and other deep redshift probes, however, make these spatial surveys strongly susceptible to cosmic scatter and will tend to bias their derived correlation lengths low. Our galaxy counts agree well with those from the HDF survey and, thus, argue against a significant inclusion of sub-galactic components in the latter census for I < 24.

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 used HST to obtain surface brightness fluctuation (SBF) observations of four nearby brightest cluster galaxies (BCG) to calibrate the BCG Hubble diagram of Lauer & Postman (1992). This BCG Hubble diagram contains 114 galaxies covering the full celestial sphere and is volume limited to 15,000 km/s, providing excellent sampling of the far field Hubble flow. The SBF zero point is based on the Cepheid calibration of the ground I_KC method (Tonry et al. 1997) as extended to the WFPC2 F814W filter by Ajhar et al. (1997). The BCG globular cluster luminosity functions give distances essentially identical to the SBF results. Using the velocities and SBF distances of the four BCG alone gives H_0 = 82 +/- 8 km/s/Mpc in the CMB frame, valid on ~4,500 km/s scales. Use of BCG as photometric redshift estimators allows the BCG Hubble diagram to be calibrated independently of recession velocities, yielding a far field H_0 = 89 +/- 10 km/s/Mpc with an effective depth of ~11,000 km/s. The error in this case is dominated by the photometric cosmic scatter in using BCG as distance estimators. The concordance of the present results with other recent H_0 determinations, and a review of theoretical treatments on perturbations in the near field Hubble flow, argue that going to the far field removes an important source of uncertainty, but that there is not a large systematic error to be corrected for to begin with. Further improvements in H_0 depend more on understanding nearby calibrators than on improved sampling of the distant flow.

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 present high-resolution, broad- and narrow-band, pre-refurbishment images of the central region of M51 taken with the Planetary Camera of the Hubble Space Telescope. The V-band images show a rather chaotic distribution of dust lanes, though some are oriented radially, roughly aligned with the major axis of the bar, and may be transporting gas to the AGN in the nucleus. The dust lane obscuring the nucleus of the galaxy, which was previously thought to be an edge-on accretion disk feeding the AGN, is not centered on the nucleus. It is unlikely that this is a stable configuration, suggesting that the material has only recently entered the nuclear region. The nucleus is contained within a cluster of stars having a total luminosity of order 5x10^7 L_0. Fitting a King model to the least obscured portions of the cluster yields a maximum core radius of 14 pc. The morphology apparent in the forbidden-line images of the extra-nuclear cloud is consistent with a narrow jet striking and scattering off the boundary of a relatively dense cocoon of gas in the disk of the galaxy. The emission-line regions are concentrated along the inner borders of dust filaments, supporting the view that the nuclear jet is ramming into and stirring up the ISM of the disk.

We present a V-I color-magnitude diagram for a region 1'-2' from the center of M32 based on Hubble Space Telescope WFPC2 images. The broad color-luminosity distribution of red giants shows that the stellar population comprises stars with a wide range in metallicity. This distribution cannot be explained by a spread in age. The blue side of the giant branch rises to M_I ~ -4.0 and can be fitted with isochrones having [Fe/H] ~ -1.5. The red side consists of a heavily populated and dominant sequence that tops out at M_I ~ -3.2, and extends beyond V-I=4. This sequence can be fitted with isochrones with -0.2 < [Fe/H] < +0.1, for ages running from 15 Gyr to 5 Gyr respectively. We do not find the optically bright asymptotic giant branch stars seen in previous ground-based work and argue that the majority of them were artifacts of crowding. Our results are consistent with the presence of the infrared-luminous giants found in ground-based studies, though their existence cannot be directly confirmed by our data. There is little evidence for an extended or even a red horizontal branch, but we find a strong clump on the giant branch itself. If the age spread is not extreme, the distribution of metallicities in M32 is considerably narrower than that of the closed-box model of chemical evolution, and also appears somewhat narrower than that of the solar neighborhood. Overall, the M32 HST color-magnitude diagram is consistent with the average luminosity-weighted age of 8.5 Gyr and [Fe/H] = -0.25 inferred from integrated spectral indices.

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

We present HST Planetary Camera V and I~band images of the central region of the peculiar giant elliptical galaxy NGC 1316. The inner profile is well fit by a nonisothermal core model with a core radius of 0.41" +/- 0.02" (34 pc). At an assumed distance of 16.9 Mpc, the deprojected luminosity density reaches $\sim 2.0 \times 10^3 L_{\sun}$ pc$^{-3}$. Outside the inner two or three arcseconds, a constant mass-to-light ratio of $\sim 2.2 \pm 0.2$ is found to fit the observed line width measurements. The line width measurements of the center indicate the existence of either a central dark object of mass $2 \times 10^9 M_{\sun}$, an increase in the stellar mass-to-light ratio by at least a factor of two for the inner few arcseconds, or perhaps increasing radial orbit anisotropy towards the center. The mass-to-light ratio run in the center of NGC 1316 resembles that of many other giant ellipticals, some of which are known from other evidence to harbor central massive dark objects (MDO's). We also examine twenty globular clusters associated with NGC 1316 and report their brightnesses, colors, and limits on tidal radii. The brightest cluster has a luminosity of $9.9 \times 10^6 L_{\sun}$ ($M_V = -12.7$), and the faintest detectable cluster has a luminosity of $2.4 \times 10^5 L_{\sun}$ ($M_V = -8.6$). The globular clusters are just barely resolved, but their core radii are too small to be measured. The tidal radii in this region appear to be $\le$ 35 pc. Although this galaxy seems to have undergone a substantial merger in the recent past, young globular clusters are not detected.

We model the surface brightness profiles of a sample of 119 Abell Brightest Cluster Galaxies (BCG), finding a generalised deVaucouleurs R^{1/n} law, where n is a free parameter, to be appropriate. Departures from the R^{1/4} law are shown to be a real feature of galaxy profiles, not due to observational errors or coupling of n with the other model parameters. BCG typically have values of n greater than 4. The shape parameter n is shown to correlate with effective half-light radius, such that the larger BCG have larger values of n. This continues a trend noticed amongst ordinary elliptical galaxies and dwarf ellipticals, such that the fainter galaxies have smaller values of n.

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

We analyze post-refurbishment Hubble Space Telescope images of four globular clusters in M31. The ability to resolve stars to below the horizontal branch permits us to use star counts to extend the surface brightness profiles determined using aperture photometry to almost 5 orders of magnitude below the central surface density. Three of the resulting cluster profiles are reasonably well-fit using single-mass King models, with core and tidal radii typical of those seen in Galactic globular clusters. We confirm an earlier report of the discovery of a cluster which has apparently undergone core collapse. Three of the four clusters show departures in their outskirts from King model behavior which, based on recent results for Galactic globulars, may indicate the presence of tidal tails.

Lauer \& Postman (LP) observe that all Abell clusters with redshifts less than 15,000\kms\ appear to be participating in a bulk flow of 689 km s$^{-1}$ with respect to the Cosmic Microwave Background. We find this result difficult to reconcile with all popular models for large-scale structure formation that assume Gaussian initial conditions. This conclusion is based on Monte-Carlo realizations of the LP data, drawn from large Particle-Mesh $N$-body simulations. We have taken special care to treat properly the longest-wavelength components of the power spectra. Bulk flows with amplitude as large as that reported by LP are not uncommon in the Monte-Carlo datasets. However, the $\chi^2$ of the observed bulk flow, taking into account the anisotropy of the error ellipsoid, is much more difficult to match in the simulations. The models examined are ruled out at confidence levels between 94\% and 98\%. Any model that has {\it intrinsic} flows of less than 480\kms\ on the scales probed by LP scales can be ruled out at a similar level.

We describe a one-parameter family of models of stable spherical stellar systems in which the phase-space distribution function depends only on energy. The models have similar density profiles in their outer parts ($\rho\propto r^{-4}$) and central power-law density cusps, $\rho\propto r^{3-\eta}$, $0<\eta\le 3$. The family contains the Jaffe (1983) and Hernquist (1990) models as special cases. We evaluate the surface brightness profile, the line-of-sight velocity dispersion profile, and the distribution function, and discuss analogs of King's core-fitting formula for determining mass-to-light ratio. We also generalize the models to a two-parameter family, in which the galaxy contains a central black hole; the second parameter is the mass of the black hole. Our models can be used to estimate the detectability of central black holes and the velocity-dispersion profiles of galaxies that contain central cusps, with or without a central black hole.