Source author record

Tod R. Lauer

Tod R. Lauer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

51works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

51 published item(s)

preprint2022arXiv

Anomalous Flux in the Cosmic Optical Background Detected With New Horizons Observations

We used New Horizons LORRI images to measure the optical-band ($0.4\lesssimλ\lesssim0.9{\rmμm}$) sky brightness within a high galactic-latitude field selected to have reduced diffuse scattered light from the Milky Way galaxy (DGL), as inferred from the IRIS all-sky $100~μ$m map. We also selected the field to significantly reduce the scattered light from bright stars (SSL) outside the LORRI field. Suppression of DGL and SSL reduced the large uncertainties in the background flux levels present in our earlier New Horizons COB results. The raw total sky level, measured when New Horizons was 51.3 AU from the Sun, is $24.22\pm0.80{\rm ~nW ~m^{-2} ~sr^{-1}}.$ Isolating the COB contribution to the raw total required subtracting scattered light from bright stars and galaxies, faint stars below the photometric detection-limit within the field, and the hydrogen plus ionized-helium two-photon continua. This yielded a highly significant detection of the COB at ${\rm 16.37\pm 1.47 ~nW ~m^{-2} ~sr^{-1}}$ at the LORRI pivot wavelength of 0.608 $μ$m. This result is in strong tension with the hypothesis that the COB only comprises the integrated light of external galaxies (IGL) presently known from deep HST counts. Subtraction of the estimated IGL flux from the total COB level leaves a flux component of unknown origin at ${\rm 8.06\pm1.92 ~nW ~m^{-2} ~sr^{-1}}.$ Its amplitude is equal to the IGL.

preprint2022arXiv

Large-scale cryovolcanic resurfacing on Pluto

The New Horizons spacecraft returned images and compositional data showing that terrains on Pluto span a variety of ages, ranging from relatively ancient, heavily cratered areas to very young surfaces with few-to-no impact craters. One of the regions with very few impact craters is dominated by enormous rises with hummocky flanks. Similar features do not exist anywhere else in the imaged solar system. Here we analyze the geomorphology and composition of the features and conclude this region was resurfaced by cryovolcanic processes, of a type and scale so far unique to Pluto. Creation of this terrain requires multiple eruption sites and a large volume of material (>104 km^3) to form what we propose are multiple, several-km-high domes, some of which merge to form more complex planforms. The existence of these massive features suggests Pluto's interior structure and evolution allows for either enhanced retention of heat or more heat overall than was anticipated before New Horizons, which permitted mobilization of water-ice-rich materials late in Pluto's history.

preprint2021arXiv

A Red-Noise Eigenbasis for the Reconstruction of Blobby Images

We demonstrate the use of an eigenbasis that is derived from principal component analysis (PCA) applied on an ensemble of random-noise images that have a "red" power spectrum; i.e., a spectrum that decreases smoothly from large to small spatial scales. The pattern of the resulting eigenbasis allows for the reconstruction of images with a broad range of image morphologies. In particular, we show that this general eigen basis can be used to efficiently reconstruct images that resemble possible astronomical sources for interferometric observations; even though the images in the original ensemble used to generate the PCA basis are significantly different from the astronomical images. We further show that the efficiency and fidelity of the image reconstructions depends only weakly on the particular parameters of the red-noise power spectrum used to generate the ensemble of images.

preprint2020arXiv

Chandra Observations of Abell 2261 Brightest Cluster Galaxy, a Candidate Host to a Recoiling Black Hole

We use Chandra X-ray observations to look for evidence of a recoiling black hole from the brightest cluster galaxy in Abell 2261 (A2261-BCG). A2261-BCG is a strong candidate for a recoiling black hole because of its large, flat stellar core, revealed by Hubble Space Telescope imaging observations. We took 100-ksec observations with Chandra and combined it with 35 ksec of archival observations to look for low-level accretion onto a black hole of expected mass $M\sim10^{10}\ M_{\scriptscriptstyle \odot}$ that could possibly be located in one of four off-center stellar knots near the galaxy's center or else in the optical center of the galaxy or in the location of radio emission. We found no X-ray emission arising from a point source in excess of the cluster gas and can place limits on the accretion of any black hole in the central region to a 2-7 keV flux below $4.3 \times 10^{-16}\ \mathrm{erg\ s^{-1}\ cm^{-2}}$, corresponding to a bolometric Eddington fraction of about $10^{-6}$. Thus there is either no $10^{10}\ M_{\scriptscriptstyle \odot}$ black hole in the core of A2261-BCG, or it is accreting at a low level. We also discuss the morphology of the X-ray emitting gas in the cluster and how its asymmetry is consistent with a large dynamic event.

preprint2020arXiv

Discretization and Filtering Effects on Black Hole Images Obtained with the Event Horizon Telescope

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.

preprint2020arXiv

Markov Chains for Horizons (MARCH). I. Identifying Biases in Fitting Theoretical Models to Event Horizon Telescope Observations

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.

preprint2020arXiv

New Horizons Observations of the Cosmic Optical Background

We used existing data from the New Horizons LORRI camera to measure the optical-band ($0.4\lesssimλ\lesssim0.9{\rmμm}$) sky brightness within seven high galactic latitude fields. The average raw level measured while New Horizons was 42 to 45 AU from the Sun is $33.2\pm0.5{\rm ~nW ~m^{-2} ~sr^{-1}}.$ This is $\sim10\times$ darker than the darkest sky accessible to the {\it Hubble Space Telescope}, highlighting the utility of New Horizons for detecting the cosmic optical background (COB). Isolating the COB contribution to the raw total requires subtracting scattered light from bright stars and galaxies, faint stars below the photometric detection-limit within the fields, and diffuse Milky Way light scattered by infrared cirrus. We remove newly identified residual zodiacal light from the IRIS $100μ$m all sky maps to generate two different estimates for the diffuse galactic light (DGL). Using these yields a highly significant detection of the COB in the range ${\rm 15.9\pm 4.2\ (1.8~stat., 3.7~sys.) ~nW ~m^{-2} ~sr^{-1}}$ to ${\rm 18.7\pm 3.8\ (1.8~stat., 3.3 ~sys.)~ nW ~m^{-2} ~sr^{-1}}$ at the LORRI pivot wavelength of 0.608 $μ$m. Subtraction of the integrated light of galaxies (IGL) fainter than the photometric detection-limit from the total COB level leaves a diffuse flux component of unknown origin in the range ${\rm 8.8\pm4.9\ (1.8 ~stat., 4.5 ~sys.) ~nW ~m^{-2} ~sr^{-1}}$ to ${\rm 11.9\pm4.6\ (1.8 ~stat., 4.2 ~sys.) ~nW ~m^{-2} ~sr^{-1}}$. Explaining it with undetected galaxies requires the galaxy-count faint-end slope to steepen markedly at $V>24$ or that existing surveys are missing half the galaxies with $V< 30.$

preprint2020arXiv

Photometry of Kuiper Belt Object (486958) Arrokoth from New Horizons LORRI

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.

preprint2016arXiv

A Hubble Space Telescope Survey for Novae in M87. I. Light and Color Curves, Spatial Distributions and the Nova Rate

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.

preprint2016arXiv

High-Resolution Mapping of Dust via Extinction in the M31 Bulge

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

preprint2016arXiv

The DESI Experiment Part I: Science,Targeting, and Survey Design

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

preprint2016arXiv

The DESI Experiment Part II: Instrument Design

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= λ/Δλ$ 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.

preprint2016arXiv

The First High-Phase Observations of a KBO: New Horizons Imaging of (15810) 1994 JR1 from the Kuiper Belt

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 $α$=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.

preprint2016arXiv

The Geology of Pluto and Charon Through the Eyes of New Horizons

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.

preprint2015arXiv

Photometric Evidence of an Intermediate-age Stellar Population in the Inner Bulge of M31

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.

preprint2015arXiv

The Panchromatic Hubble Andromeda Treasury VIII: A Wide-Area, High-Resolution Map of Dust Extinction in M31

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.

preprint2015arXiv

Tracing the Metal-Poor M31 Stellar Halo with Blue Horizontal Branch Stars

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

preprint2014arXiv

Brightest Cluster Galaxies at the Present Epoch

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.

preprint2014arXiv

The Panchromatic Hubble Andromeda Treasury X. Ultraviolet to Infrared Photometry of 117 Million Equidistant Stars

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.

preprint2014arXiv

The Panchromatic Hubble Andromeda Treasury. VII. The Steep Mid-Ultraviolet to Near-Infrared Extinction Curve in the Central 200 pc of the M31 Bulge

We measure the extinction curve in the central 200 pc of M31 at mid-ultraviolet to near-infrared wavelengths (from 1928A to 1.5μ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.

preprint2013arXiv

The Black Hole Mass and the Stellar Ring in NGC 3706

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 ($Υ_V$), and dark matter halo profile. The best-fit model parameters with 1$σ$ uncertainties are $M = (6.0^{+0.7}_{-0.9}) \times 10^8\ M_{\scriptscriptstyle \odot}$ and $Υ_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 $Δχ^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|/σ$, the ratio of mean velocity to velocity dispersion. Based on the observed low $|V|/σ$, 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.

preprint2013arXiv

The Effect of Spatial Gradients in Stellar Mass-to-Light Ratio on Black Hole Mass Measurements

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

preprint2012arXiv

A Brightest Cluster Galaxy with an Extremely Large Flat Core

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.

preprint2012arXiv

A highly magnified candidate for a young galaxy seen when the Universe was 500 Myrs old

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

preprint2012arXiv

An Ancient Metal-Poor Population in M32, and Halo Satellite Accretion in M31, Identified by RR Lyrae Stars

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]

preprint2012arXiv

Cores and the Kinematics of Early-Type Galaxies

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.

preprint2012arXiv

Dynamical Measurements of Black Hole Masses in Four Brightest Cluster Galaxies at 100 Mpc

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.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury I: Bright UV Stars in the Bulge of M31

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é 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é (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 α 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.

preprint2012arXiv

The Panchromatic Hubble Andromeda Treasury II. Tracing the Inner M31 Halo with Blue Horizontal Branch Stars

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^(-α) for the 2-D projected surface density BHB distribution, we obtain a high-quality fit with a 2-D power-law index of α=2.6^{+0.3}_{-0.2} outside of 3 kpc, which flattens to α<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.

preprint2011arXiv

Is there a black hole in NGC 4382?

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σ) or M_BH < 1.4 \times 10^8 M_sun (3σ) is consistent with the current M-σ relation, which predicts M_BH = 8.8 \times 10^7 M_sun at σ_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.

preprint2011arXiv

Orbit-Based Dynamical Models of the Sombrero Galaxy (NGC 4594)

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 σ_e puts NGC 4594 squarely on the M-σrelation. We also determine that NGC 4594 lies directly on the M-L relation.

preprint2011arXiv

The Black-Hole Mass in M87 from Gemini/NIFS Adaptive Optics Observations

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.

preprint2011arXiv

The Cluster of Blue Stars Surrounding the M31 Nuclear Black Hole

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.

preprint2011arXiv

The Star Formation History of M32

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.

preprint2011arXiv

Two ten-billion-solar-mass black holes at the centres of giant elliptical galaxies

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.

preprint2010arXiv

Exploring the Spectral Space of Low Redshift QSOs

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.

preprint2010arXiv

The Black Hole Mass in Brightest Cluster Galaxy NGC 6086

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.

preprint2010arXiv

The Deepest HST Color-Magnitude Diagram of M32: Evidence for Intermediate-Age Populations

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]

preprint2009arXiv

A Quintet of Black Hole Mass Determinations

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.

preprint2009arXiv

RR Lyrae variables in M32 and the disk of M31

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.

preprint2009arXiv

RR Lyrae Variables in Two Fields in the Spheroid of M31

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.

preprint2009arXiv

The M-sigma and M-L Relations in Galactic Bulges and Determinations of their Intrinsic Scatter

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.

preprint2008arXiv

A Stellar Dynamical Measurement of the Black Hole Mass in the Maser Galaxy NGC 4258

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.

preprint2008arXiv

Color Bimodality in M87 Globular Clusters

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.

preprint2007arXiv

The Masses of Nuclear Black Holes in Luminous Elliptical Galaxies and Implications for the Space Density of the Most Massive Black Holes

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.

preprint2006arXiv

Luminosity Function of Faint Globular Clusters in M87

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.

preprint2006arXiv

The Centers of Early-Type Galaxies with HST. VI. Bimodal Central Surface Brightness Profiles

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.

preprint2000arXiv

A Relationship Between Nuclear Black Hole Mass and Galaxy Velocity Dispersion

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.

preprint1996arXiv

Hubble Space Telescope Observations of M32: The Color-Magnitude Diagram

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.

preprint1996arXiv

Hubble Space Telescope Planetary Camera Images of NGC 1316

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.

preprint1996arXiv

The Nuclear Region of M51 Imaged with the HST Planetary Camera

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.