Source author record

Monica Valluri

Monica Valluri 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

39works
9topics
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

39 published item(s)

preprint2022arXiv

A Spectroscopic Road Map for Cosmic Frontier: DESI, DESI-II, Stage-5

In this white paper, we present an experimental road map for spectroscopic experiments beyond DESI. DESI will be a transformative cosmological survey in the 2020s, mapping 40 million galaxies and quasars and capturing a significant fraction of the available linear modes up to z=1.2. DESI-II will pilot observations of galaxies both at much higher densities and extending to higher redshifts. A Stage-5 experiment would build out those high-density and high-redshift observations, mapping hundreds of millions of stars and galaxies in three dimensions, to address the problems of inflation, dark energy, light relativistic species, and dark matter. These spectroscopic data will also complement the next generation of weak lensing, line intensity mapping and CMB experiments and allow them to reach their full potential.

preprint2022arXiv

Gas inflows in the polar ring of NGC 4111: the birth of an AGN

We have used Hubble Space Telescope (HST) images, SAURON Integral Field Spectroscopy (IFS) and adaptative optics assisted Gemini NIFS near-infrared K-band IFS to map the stellar and gas distribution, excitation and kinematics of the inner few kpc of the nearby edge-on S0 galaxy NGC 4111. The HST images map its $\approx$ 450 pc diameter dusty polar ring, with an estimated gas mass $\ge10^7$ M$_\odot$. The NIFS datacube maps the inner 110 pc radius at $\approx$ 7 pc spatial resolution revealing a $\approx$ 220 pc diameter polar ring in hot ($2267\pm166$ K) molecular H$_2$ 1-0 S(1) gas embedded in the polar ring. The stellar velocity field shows disk-dominated kinematics along the galaxy plane both in the SAURON large-scale and in the NIFS nuclear-scale data. The large-scale [O III] $\lambda5007$ Åvelocity field shows a superposition of two disk kinematics: one similar to that of the stars and another along the polar ring, showing non-circular motions that seem to connect with the velocity field of the nuclear H$_2$ ring, whose kinematics indicate accelerated inflow to the nucleus. The estimated mass inflow rate is enough not only to feed an Active Galactic Nucleus (AGN) but also to trigger circumnuclear star formation in the near future. We propose a scenario in which gas from the polar ring, which probably originated from the capture of a dwarf galaxy, is moving inwards and triggering an AGN, as supported by the local X-ray emission, which seems to be the source of the H$_2$ 1-0 S(1) excitation. The fact that we see neither near-UV nor Br$γ$ emission suggests that the nascent AGN is still deeply buried under the optically thick dust of the polar ring.

preprint2022arXiv

New constraints on the dark matter density profiles of dwarf galaxies from proper motions of globular cluster streams

The central density profiles in dwarf galaxy halos depend strongly on the nature of dark matter. Recently, in Malhan et al. (2021), we employed N-body simulations to show that the cuspy cold dark matter (CDM) subhalos predicted by cosmological simulations can be differentiated from cored subhalos using the properties of accreted globular cluster (GC) streams since these GCs experience tidal stripping within their parent halos prior to accretion onto the Milky Way. We previously found that clusters that are accreted within cuspy subhalos produce streams with larger physical widths and higher dispersions in line-of-sight velocity and angular momentum than streams that are accreted within cored subhalos. Here, we use the same suite of simulations to demonstrate that the dispersion in the tangential velocities of streams ($σ_{v_\mathrm{Tan}}$) is also sensitive to the central DM density profiles of their parent dwarfs and GCs that were accreted from: cuspy subhalos produce streams with larger $σ_{v_\mathrm{Tan}}$ than those accreted inside cored subhalos. Using Gaia EDR3 observations of multiple GC streams we compare their $σ_{v_\mathrm{Tan}}$ values with simulations. The measured $σ_{v_\mathrm{Tan}}$ values are consistent with both an ``in situ'' origin and with accretion inside cored subhalos of $M\sim 10^{8-9}M_{\odot}$ (or very low-mass cuspy subhalos of mass $\sim 10^8M_{\odot}$). Despite the large current uncertainties in $σ_{v_\mathrm{Tan}}$, we find a low probability that any of the progenitor GCs were accreted from cuspy subhalos of $M_{\rm subhalo}\buildrel > \over \sim$ $10^9 M_{\odot}$. The uncertainties on Gaia tangential velocity measurements are expected to decrease in future and will allow for stronger constraints on subhalo DM density profiles.

preprint2022arXiv

Non-Parametric Spherical Jeans Mass Estimation with B-splines

Spherical Jeans modeling is widely used to estimate mass profiles of systems from star clusters to galactic stellar haloes to clusters of galaxies. It derives the cumulative mass profile, M(<r), from kinematics of tracers of the potential under the assumptions of spherical symmetry and dynamical equilibrium. We consider the application of Jeans modeling to mapping the dark matter distribution in the outer reaches of the Milky Way using field halo stars. We present a novel non-parametric routine for solving the spherical Jeans equation by fitting B-splines to the velocity and density profiles of halo stars. While most implementations assume parametric forms for these profiles, B-splines provide non-parametric fitting curves with analytical derivatives. Our routine recovers the mass profiles of equilibrium systems with flattened haloes or a stellar disc and bulge excellently (<~ 10% error at most radii). Tests with non-equilibrium, Milky Way-like galaxies from the Latte suite of FIRE-2 simulations perform quite well (<~ 15% error for r <~ 100 kpc). We also create observationally motivated datasets for the Latte suite by imposing selection functions and errors on phase space coordinates characteristic of Gaia and the DESI Milky Way Survey. The resulting imprecise and incomplete data require us to introduce an MCMC-based subroutine to obtain deconvolved density and velocity dispersion profiles from the tracer population. With these observational effects taken into account, the accuracy of the Jeans mass estimate remains at the level 20% or better.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Prospects for obtaining Dark Matter Constraints with DESI

Despite efforts over several decades, direct-detection experiments have not yet led to the discovery of the dark matter (DM) particle. This has led to increasing interest in alternatives to the Lambda CDM (LCDM) paradigm and alternative DM scenarios (including fuzzy DM, warm DM, self-interacting DM, etc.). In many of these scenarios, DM particles cannot be detected directly and constraints on their properties can ONLY be arrived at using astrophysical observations. The Dark Energy Spectroscopic Instrument (DESI) is currently one of the most powerful instruments for wide-field surveys. The synergy of DESI with ESA's Gaia satellite and future observing facilities will yield datasets of unprecedented size and coverage that will enable constraints on DM over a wide range of physical and mass scales and across redshifts. DESI will obtain spectra of the Lyman-alpha forest out to z~5 by detecting about 1 million QSO spectra that will put constraints on clustering of the low-density intergalactic gas and DM halos at high redshift. DESI will obtain radial velocities of 10 million stars in the Milky Way (MW) and Local Group satellites enabling us to constrain their global DM distributions, as well as the DM distribution on smaller scales. The paradigm of cosmological structure formation has been extensively tested with simulations. However, the majority of simulations to date have focused on collisionless CDM. Simulations with alternatives to CDM have recently been gaining ground but are still in their infancy. While there are numerous publicly available large-box and zoom-in simulations in the LCDM framework, there are no comparable publicly available WDM, SIDM, FDM simulations. DOE support for a public simulation suite will enable a more cohesive community effort to compare observations from DESI (and other surveys) with numerical predictions and will greatly impact DM science.

preprint2022arXiv

The MegaMapper: A Stage-5 Spectroscopic Instrument Concept for the Study of Inflation and Dark Energy

In this white paper, we present the MegaMapper concept. The MegaMapper is a proposed ground-based experiment to measure Inflation parameters and Dark Energy from galaxy redshifts at $2<z<5$. In order to achieve path-breaking results with a mid-scale investment, the MegaMapper combines existing technologies for critical path elements and pushes innovative development in other design areas. To this aim, we envision a 6.5-m Magellan-like telescope, with a newly designed wide field, coupled with DESI spectrographs, and small-pitch robots to achieve multiplexing of at least 26,000. This will match the expected achievable target density in the redshift range of interest and provide a 10x capability over the existing state-of the art, without a 10x increase in project budget.

preprint2020arXiv

Robotic Reverberation Mapping of the Southern Seyfert NGC 3783

We present spectroscopic and photometric monitoring of NGC 3783 conducted throughout the first half of 2020. Time delays between the continuum variations and the response of the broad optical emission lines were clearly detected, and we report reverberation measurements for H$β$, HeII $λ4686$, H$γ$, and H$δ$. From the time delay in the broad H$β$ emission line and the line width in the variable portion of the spectrum, we derive a black hole mass of $M_{\rm BH} = 2.34^{+0.43}_{-0.43} \times 10^7$ M$_{\odot}$. This is slightly smaller than, but consistent with, previous determinations. However, our significantly improved time sampling ($T_{\rm med}=1.7$ days compared to $T_{\rm med}=4.0$ days) has reduced the uncertainties on both the time delay and the derived mass by $\sim 50$%. We also detect clear velocity-resolved time delays across the broad H$β$ profile, with shorter lags in the line wings and a longer lag in the line core. Future modeling of the full velocity-resolved time delay response will further improve the reverberation-based mass for NGC 3783, adding it to the small but growing sample of AGNs for which we have direct, primary black hole mass measurements. Upcoming MUSE observations at VLT will also allow NGC 3783 to join the smaller sample of black holes where reverberation masses and masses from stellar dynamical modeling may be directly compared.

preprint2019arXiv

A Cepheid-Based Distance to the Seyfert Galaxy NGC 6814

We present a Cepheid-based distance to the nearby Seyfert galaxy NGC\,6814 from {\it Hubble Space Telescope} observations. We obtained F555W and F814W imaging over the course of 12 visits with logarithmic time spacing in 2013 August$-$October. We detected and made photometric measurements for 16,469 unique sources across all images in both filters, from which we identify 90 excellent Cepheid candidates spanning a range of periods of $13-84$ days. We find evidence for incompleteness in the detection of candidates at periods <21 days. Based on the analysis of Cepheid candidates above the incompleteness limit, we determine a distance modulus for NGC\,6814 relative to the LMC of $μ_{\rm rel\,LMC}=13.200^{+0.031}_{-0.031}$ mag. Adopting the recent constraint of the distance modulus to the LMC determined by Pietrzynski et al., we find $m-M=31.677^{+0.041}_{-0.041}$ which gives a distance of $21.65 \pm 0.41$ Mpc to NGC 6814.

preprint2019arXiv

A new implementation of the Schwarzschild method for constructing observationally-driven dynamical models of galaxies of all morphological types

We present Forstand, a new code for constructing dynamical models of galaxies with the Schwarzschild orbit-superposition method. These models are constrained by line-of-sight kinematic observations and applicable to galaxies of all morphological types, including disks and triaxial rotating bars. Our implementation has several novel and improved features, is computationally efficient, and made publicly available. Using mock datasets taken from N-body simulations, we demonstrate that the pattern speed of a bar can be recovered with an accuracy of 10-20%, regardless of orientation, if the 3D shape of the galaxy is known or inferred correctly.

preprint2019arXiv

Schwarzschild modeling of barred galaxies

We review the Schwarzschild orbit-superposition approach and present a new implementation of this method, which can deal with a large class of systems, including rotating barred disk galaxies. We discuss two conceptuals problems in this field: the intrinsic degeneracy of determining the potential from line-of-sight kinematics, and the non-uniqueness of deprojection and related biases in potential inference, especially acute for triaxial bars. When applied to mock datasets with known 3d shape, our method correctly recovers the pattern speed and other potential parameters. However, more work is needed to systematically address these two problems for real observational datasets.

preprint2019arXiv

The shape of the dark matter halo revealed from a hypervelocity star

A recently discovered young, high-velocity giant star J01020100-7122208 is a good candidate of hypervelocity star ejected from the Galactic center, although it has a bound orbit. If we assume that this star was ejected from the Galactic center, it can be used to constrain the Galactic potential, because the deviation of its orbit from a purely radial orbit informs us of the torque that this star has received after its ejection. Based on this assumption, we estimate the flattening of the dark matter halo of the Milky Way by using the Gaia DR2 data and the circular velocity data from Eilers et al. (2019). Our Bayesian analysis shows that the orbit of J01020100-7122208 favors a prolate dark matter halo within $\sim$ 10 kpc from the Galactic center. The posterior distribution of the density flattening $q$ shows a broad distribution at $q \gtrsim1$ and peaks at $q \simeq 1.5$. Also, 98.5\% of the posterior distribution is located at $q>1$, highly disfavoring an oblate halo.

preprint2016arXiv

Chaotic Dispersal of Tidal Debris

Several long, dynamically cold stellar streams have been observed around the Milky Way Galaxy, presumably formed from the tidal disruption of globular clusters. In integrable potentials---where all orbits are regular---tidal debris phase-mixes close to the orbit of the progenitor system. However, the Milky Way's dark matter halo is expected not to be fully integrable; an appreciable fraction of orbits will be chaotic. This paper examines the influence of chaos on the phase-space morphology of cold tidal streams. Streams even in weakly chaotic regions look very different from those in regular regions. We find that streams can be sensitive to chaos on a much shorter time-scale than any standard prediction (from the Lyapunov or frequency-diffusion times). For example, on a weakly chaotic orbit with a chaotic timescale predicted to be $>$1000 orbital periods ($>$1000 Gyr), the resulting stellar stream is, after just a few 10's of orbits, substantially more diffuse than any formed on a nearby but regular orbit. We find that the enhanced diffusion of the stream stars can be understood by looking at the variance in orbital frequencies of orbit ensembles centered around the parent (progenitor) orbit. Our results suggest that long, cold streams around our Galaxy must exist only on regular (or very nearly regular) orbits; they potentially provide a map of the regular regions of the Milky Way potential. This suggests a promising new direction for the use of tidal streams to constrain the distribution of dark matter around our Galaxy.

preprint2016arXiv

Dark Stars: A Review

Dark Stars are stellar objects made (almost entirely) of hydrogen and helium, but powered by the heat from Dark Matter annihilation, rather than by fusion. They are in hydrostatic and thermal equilibrium, but with an unusual power source. Weakly Interacting Massive Particles (WIMPs), among the best candidates for dark matter, can be their own antimatter and can annihilate inside the star, thereby providing a heat source. Although dark matter constitutes only $\lesssim 0.1\%$ of the stellar mass, this amount is sufficient to power the star for millions to billions of years. Thus, the first phase of stellar evolution in the history of the Universe may have been dark stars. We review how dark stars come into existence, how they grow as long as dark matter fuel persists, and their stellar structure and evolution. The studies were done in two different ways, first assuming polytropic interiors and more recently using the MESA stellar evolution code; the basic results are the same. Dark stars are giant, puffy ($\sim$ 10 AU) and cool (surface temperatures $\sim$10,000 K) objects. We follow the evolution of dark stars from their inception at $\sim 1 M_\odot$ as they accrete mass from their surroundings to become supermassive stars, some even reaching masses $> 10^6 M_\odot$ and luminosities $>10^{10} L_\odot$, making them detectable with the upcoming James Webb Space Telescope. Once the dark matter runs out and the dark star dies, it may collapse to a black hole; thus dark stars may provide seeds for the supermassive black holes observed throughout the Universe and at early times. Other sites for dark star formation may exist in the Universe today in regions of high dark matter density such as the centers of galaxies. The current review briefly discusses dark stars existing today, but focuses on the early generation of dark stars.

preprint2016arXiv

Detailed Chemical Abundances in the r-Process-Rich Ultra-Faint Dwarf Galaxy Reticulum 2

The ultra-faint dwarf galaxy Reticulum 2 (Ret 2) was recently discovered in images obtained by the Dark Energy Survey. We have observed the four brightest red giants in Ret 2 at high spectral resolution using the Michigan/Magellan Fiber System. We present detailed abundances for as many as 20 elements per star, including 12 elements heavier than the Fe group. We confirm previous detection of high levels of r-process material in Ret 2 (mean [Eu/Fe]=+1.69+/-0.05) found in three of these stars (mean [Fe/H]=-2.88+/-0.10). The abundances closely match the r-process pattern found in the well-studied metal-poor halo star CS22892-052. Such r-process-enhanced stars have not been found in any other ultra-faint dwarf galaxy, though their existence has been predicted by at least one model. The fourth star in Ret 2 ([Fe/H]=-3.42+/-0.20) contains only trace amounts of Sr ([Sr/Fe]=-1.73+/-0.43) and no detectable heavier elements. One r-process enhanced star is also enhanced in C (natal [C/Fe]=+1.1). This is only the third such star known, which suggests that the nucleosynthesis sites leading to C and r-process enhancements are decoupled. The r-process-deficient star is enhanced in Mg ([Mg/Fe]=+0.81+/-0.14), and the other three stars show normal levels of alpha-enhancement (mean [Mg/Fe]=+0.34+/-0.03). The abundances of other alpha and Fe-group elements closely resemble those in ultra-faint dwarf galaxies and metal-poor halo stars, suggesting that the nucleosynthesis that led to the large r-process enhancements either produced no light elements or produced light-element abundance signatures indistinguishable from normal supernovae.

preprint2016arXiv

Imprints of radial migration on the Milky Way's metallicity distribution functions

Recent analysis of the SDSS-III/APOGEE Data Release 12 stellar catalogue has revealed that the Milky Way's metallicity distribution function (MDF) changes shape as a function of radius, transitioning from being negatively skewed at small Galactocentric radii to positively skewed at large Galactocentric radii. Using a high resolution, N-body+SPH simulation, we show that the changing skewness arises from radial migration - metal-rich stars form in the inner disk and subsequently migrate to the metal-poorer outer disk. These migrated stars represent a large fraction (> 50%) of the stars in the outer disk; they populate the high metallicity tail of the MDFs and are, in general, more metal-rich than the surrounding outer disk gas. The simulation also reproduces another surprising APOGEE result: the spatially invariant high-[alpha/Fe] MDFs. This arises in the simulation from the migration of a population formed within a narrow range of radii (3.2+/-1.2 kpc) and time (8.8+/-0.6 Gyr ago), rather than from spatially extended star formation in a homogeneous medium at early times. These results point toward the crucial role radial migration has played in shaping our Milky Way.

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 impact of baryons on the direct detection of dark matter

The spatial and velocity distributions of dark matter particles in the Milky Way Halo affect the signals expected to be observed in searches for dark matter. Results from direct detection experiments are often analyzed assuming a simple isothermal distribution of dark matter, the Standard Halo Model (SHM). Yet there has been skepticism regarding the validity of this simple model due to the complicated gravitational collapse and merger history of actual galaxies. In this paper we compare the SHM to the results of cosmological hydrodynamical simulations of galaxy formation to investigate whether or not the SHM is a good representation of the true WIMP distribution in the analysis of direct detection data. We examine two Milky Way-like galaxies from the MaGICC cosmological simulations (a) with dark matter only and (b) with baryonic physics included. The inclusion of baryons drives the shape of the DM halo to become more spherical and makes the velocity distribution of dark matter particles less anisotropic especially at large heliocentric velocities, thereby making the SHM a better fit. We also note that we do not find a significant disk-like rotating dark matter component in either of the two galaxy halos with baryons that we examine, suggesting that dark disks are not a generic prediction of cosmological hydrodynamical simulations. We conclude that in the Solar neighborhood, the SHM is in fact a good approximation to the true dark matter distribution in these cosmological simulations (with baryons) which are reasonable representations of the Milky Way, and hence can also be used for the purpose of dark matter direct detection calculations.

preprint2015arXiv

A unified framework for the orbital structure of bars and triaxial ellipsoids

We examine a large random sample of orbits in self-consistent simulations of N-body bars. Orbits in the bars are classified both visually and with a new automated orbit classification method based on frequency analysis. The well known prograde x1 orbit family originates from the same parent orbit as the box orbits in stationary and rotating triaxial ellipsoids. However only a small fraction of bar orbits ~4% have predominately prograde motion like their periodic parent orbit. Most bar orbits arising from the x1 orbit have little net angular momentum in the bar frame making them equivalent to box orbits in rotating triaxial potentials. A small fraction of bar orbits (~7%) are long axis tubes that behave exactly like those in triaxial ellipsoids:they are tipped about the intermediate-axis due to the Coriolis force, with the sense of tipping determined by the sign of their angular momentum about the long axis. No orbits parented by prograde periodic x2 orbits are found in the pure bar model, but a tiny population (~2%) of short axis tube orbits parented by retrograde x4 orbits are found. When a central point mass representing a supermassive black hole (SMBH) is grown adiabatically at the center of the bar, those orbits that lie in the immediate vicinity of the SMBH are transformed into precessing Keplerian orbits (PKOs) which belong the same major families (short axis tubes, long axis tubes and boxes) occupying the bar at larger radii. During the growth of a SMBH the inflow of mass and outward transport of angular momentum transforms some x1 and long axis tube orbits into prograde short axis tubes. This study has important implications for future attempts to constrain the masses of SMBHs in barred galaxies using orbit based methods like the Schwarzschild orbit superposition scheme and for understanding the observed features in barred galaxies.

preprint2015arXiv

Stellar orbital properties as diagnostics of the origin of the stellar halo

We examine metallicities, ages and orbital properties of halo stars in a Milky-Way like disk galaxy formed in the cosmological hydrodynamical MaGICC simulations. Halo stars were either accreted from satellites or they formed in situ in the disk or bulge of the galaxy and were then kicked up into the halo ("in situ/ kicked-up" stars). Regardless of where they formed both types show surprisingly similar orbital properties: the majority of stars of both types are on short-axis tubes with the same sense of rotation as the disk -- implying that a large fraction of satellites are accreted onto the halo with the same sense of angular momentum as the disk.

preprint2014arXiv

AMUSE-Field II. Nucleation of early-type galaxies in the field vs. cluster environment

The optical light profiles of nearby early type galaxies are known to exhibit a smooth transition from nuclear light deficits to nuclear light excesses with decreasing galaxy mass, with as much as 80 per cent of the galaxies with stellar masses below 10^10 Msun hosting a massive nuclear star cluster. At the same time, while all massive galaxies are thought to harbor nuclear super-massive black holes (SMBHs), observational evidence for SMBHs is slim at the low end of the mass function. Here, we explore the environmental dependence of the nucleation fraction by comparing two homogeneous samples of nearby field vs. cluster early type galaxies with uniform Hubble Space Telescope (HST) coverage. Existing Chandra X-ray Telescope data for both samples yield complementary information on low-level accretion onto nuclear SMBHs. Specifically, we report on dual-band (F475W & F850LP) Advanced Camera for Surveys (ACS) imaging data for 28 out of the 103 field early type galaxies that compose the AMUSE-Field Chandra survey, and compare our results against the companion HST and Chandra surveys for a sample of 100 Virgo cluster early types (ACS Virgo Cluster and AMUSE-Virgo surveys, respectively). We model the two-dimensional light profiles of the field targets to identify and characterize NSCs, and find a field nucleation fraction of 26% (+17%, -11%; at the 1-sigma level), consistent with the measured Virgo nucleation fraction of 30% (+17%, -12%), across a comparable mass distribution. Coupled with the Chandra result that SMBH activity is higher for the field, our findings indicate that, since the last epoch of star formation, the funneling of gas to the nuclear regions has been inhibited more effectively for Virgo galaxies, arguably via ram pressure stripping.

preprint2014arXiv

Systematic Problems With Using Dark Matter Simulations to Model Stellar Halos

The limits of available computing power have forced models for the structure of stellar halos to adopt one or both of the following simplifying assumptions: (1) stellar mass can be "painted" onto dark matter particles in progenitor satellites; (2) pure dark matter simulations that do not form a luminous galaxy can be used. We estimate the magnitude of the systematic errors introduced by these assumptions using a controlled set of stellar halo models where we independently vary whether we look at star particles or painted dark matter particles, and whether we use a simulation in which a baryonic disk galaxy forms or a matching pure dark matter simulation that does not form a baryonic disk. We find that the "painting" simplification reduces the halo concentration and internal structure, predominantly because painted dark matter particles have different kinematics than star particles even when both are buried deep in the potential well of the satellite. The simplification of using pure dark matter simulations reduces the concentration further, but increases the internal structure, and results in a more prolate stellar halo. These differences can be a factor of 1.5-7 in concentration (as measured by the half-mass radius) and 2-7 in internal density structure. Given this level of systematic uncertainty, one should be wary of overinterpreting differences between observations and the current generation of stellar halo models based on dark matter only simulations when such differences are less than an order of magnitude.

preprint2014arXiv

The Black Hole Mass of NGC 4151. II. Stellar Dynamical Measurement from Near-Infrared Integral Field Spectroscopy

We present a revised measurement of the mass of the central black hole (Mbh) in the Seyfert 1 galaxy NGC 4151. The new stellar dynamical mass measurement is derived by applying an axisymmetric orbit-superposition code to near-infrared integral field data obtained using adaptive optics with the Gemini NIFS spectrograph. When our models attempt to fit both the NIFS kinematics and additional low spatial resolution kinematics, our results depend sensitively on how chi-squared is computed--probably a consequence of complex bar kinematics that manifest immediately outside the nuclear region. The most robust results are obtained when only the high spatial resolution kinematic constraints in the nuclear region are included in the fit. Our best estimates for the BH mass and H-band mass-to-light ratio are Mbh~(3.76+/-1.15)E7 Msun (1-sigma error) and M/L(H-band)~0.34+/-0.03 Msun/Lsun (3-sigma error), respectively (the quoted errors reflect the model uncertainties). Our BH mass measurement is consistent with estimates from both reverberation mapping (3.57[+0.45/-0.37]E7 Msun) and gas kinematics (3.0[+0.75/-2.2]E7 Msun; 1-sigma errors), and our best-fit mass-to-light ratio is consistent with the photometric estimate of M/L(H-band)=0.4+/-0.2 Msun/Lsun. The NIFS kinematics give a central bulge velocity dispersion sigma_c=116+/-3 km/s, bringing this object slightly closer to the M-sigma relation for quiescent galaxies. Although NGC 4151 is one of only a few Seyfert 1 galaxies in which it is possible to obtain a direct dynamical BH mass measurement--and thus, an independent calibration of the reverberation mapping mass scale--the complex bar kinematics makes it less than ideally suited for this purpose.

preprint2013arXiv

Halo orbits in cosmological disk galaxies: tracers of formation history

We analyze the orbits of stars and dark matter particles in the halo of a disk galaxy formed in a cosmological hydrodynamical simulation. The halo is oblate within the inner ~20 kpc and triaxial beyond this radius. About 43% of orbits are short axis tubes - the rest belong to orbit families that characterize triaxial potentials (boxes, long-axis tubes and chaotic orbits), but their shapes are close to axisymmetric. We find no evidence that the self-consistent distribution function of the nearly oblate inner halo is comprised primarily of axisymmetric short-axis tube orbits. Orbits of all families, and both types of particles are highly eccentric with mean eccentricity >0.6. We find that randomly selected samples of halo stars show no substructure in "integrals of motion" space. However individual accretion events can be clearly identified in plots of metallicity versus formation time. Dynamically young tidal debris is found primarily on a single type of orbit. However, stars associated with older satellites become chaotically mixed during the formation process (possibly due to scattering by the central bulge and disk, and baryonic processes), and appear on all four types of orbits. We find that the tidal debris in cosmological hydrodynamical simulations experiences significantly more chaotic evolution than in collisionless simulations, making it much harder to identify individual progenitors using phase space coordinates alone. However by combining information on stellar ages and chemical abundances with the orbital properties of halo stars in the underlying self-consistent potential, the identification of progenitors is likely to be possible.

preprint2013arXiv

On the Offset of Barred Galaxies From the Black Hole $M_{BH}$-$σ$ Relationship

We use collisionless $N$-body simulations to determine how the growth of a supermassive black hole (SMBH) influences the nuclear kinematics in both barred and unbarred galaxies. In the presence of a bar, the increase in the velocity dispersion $σ$ (within the effective radius) due to the growth of an SMBH is on average $\lesssim 10%$, whereas the increase is only $\lesssim 4%$ in an unbarred galaxy. In a barred galaxy, the increase results from a combination of three separate factors (a) orientation and inclination effects; (b) angular momentum transport by the bar that results in an increase in the central mass density; (c) an increase in the vertical and radial velocity anisotropy of stars in the vicinity of the SMBH. In contrast the growth of the SMBH in an unbarred galaxy causes the velocity distribution in the inner part of the nucleus to become less radially anisotropic. The increase in $σ$ following the growth of the SMBH is insensitive to a variation of a factor of 10 in the final mass of the SMBH, showing that it is the growth process rather than the actual SMBH mass that alters bar evolution in a way that increases $σ$. We argue that using an axisymmetric stellar dynamical modeling code to measure SMBH masses in barred galaxies could result in a slight overestimate of the derived $M_{BH}$, especially if a constant M/L ratio is assumed. We conclude that the growth of a black hole in the presence of a bar could result in an increase in $σ$ which is roughly of $4-8%$ larger than the increase that occurs in an axisymmetric system. While the increase in $σ$ due to SMBH growth in a barred galaxy might partially account for the claimed offset of barred galaxies and pseudo bulges from the $M_{BH}-σ$ relation obtained for elliptical galaxies and classical bulges in unbarred galaxies, it is inadequate to account for all of the offset.

preprint2013arXiv

The effect of bars on the M-sigma relation: offset, scatter and residuals correlations

We analyse a set of collisionless disc galaxy simulations to study the consequences of bar formation and evolution on the M-sigma relation of supermassive black holes. The redistribution of angular momentum driven by bars leads to a mass increase within the central region, raising the velocity dispersion of the bulge, sigma, on average by ~12% and as much as ~20%. If a disc galaxy with a SMBH satisfying the M-sigma relation forms a bar, and the SMBH does not grow in the process, then the increase in sigma moves the galaxy off the M-sigma relation. We explore various effects that can affect this result including contamination from the disc and anisotropy. The displacement from the M-sigma relation for individual model barred galaxies correlates with both M(B)/M(B+D) and beta_phi(B+D) measured within the effective radius of the bulge. Overall, this process leads to an M-sigma for barred galaxies offset from that of unbarred galaxies, as well as an increase in its scatter. We assemble samples of unbarred and barred galaxies with classical bulges and find tentative hints of an offset between the two consistent with the predicted. Including all barred galaxies, rather than just those with a classical bulge, leads to a significantly larger offset.

preprint2013arXiv

What's Up in the Milky Way? The Orientation of the Disc Relative to the Triaxial Halo

Models of the Sagittarius Stream have consistently found that the Milky Way disc is oriented such that its short axis is along the intermediate axis of the triaxial dark matter halo. We attempt to build models of disc galaxies in such an `intermediate-axis orientation'. We do this with three models. In the first two cases we simply rigidly grow a disc in a triaxial halo such that the disc ends up perpendicular to the global intermediate axis. We also attempt to coax a disc to form in an intermediate-axis orientation by producing a gas+dark matter triaxial system with gas angular momentum about the intermediate axis. In all cases we fail to produce systems which remain with stellar angular momentum aligned with the halo's intermediate axis, even when the disc's potential flattens the inner halo such that the disc is everywhere perpendicular to the halo's local minor axis. For one of these unstable simulations we show that the potential is even rounder than the models of the Milky Way potential in the region probed by the Sagittarius Stream. We conclude that the Milky Way's disc is very unlikely to be in an intermediate axis orientation. However we find that a disc can persist off one of the principal planes of the potential. We propose that the disc of the Milky Way must be tilted relative to the principal axes of the dark matter halo. Direct confirmation of this prediction would constitute a critical test of Modified Newtonian Dynamics.

preprint2011arXiv

Observing Dark Stars with JWST

We study the capability of the James Webb Space Telescope (JWST) to detect Supermassive Dark Stars (SMDS). If the first stars are powered by dark matter heating in triaxial dark matter haloes, they may grow to be very large and very bright, visible in deep imaging with JWST and even Hubble Space Telescope (HST). We use HST surveys to place bounds on the numbers of SMDSs that may be detected in future JWST imaging surveys. We showed that SMDS in the mass range $10^6-10^7 M_\odot$ are bright enough to be detected in all the wavelength bands of the NIRCam on JWST . If SMDSs exist at z ~10, 12, and 14, they will be detectable as J-band, H-band, or K-band dropouts, respectively. With a total survey area of 150 arcmin^2 (assuming a multi-year deep parallel survey with JWST), we find that typically the number of $10^6 M_\odot$ SMDSs found as H or K-band dropouts is ~10^5\fsmds, where the fraction of early DM haloes hosting DS is likely to be small, \fsmds<<1. If the SDMS survive down to z=10 where HST bounds apply, then the observable number of SMDSs as H or K-band dropouts with JWST is ~1-30. While individual SMDS are bright enough to be detected by JWST, standard PopIII stars are not, and would only be detected in first galaxies with total stellar masses of ~$10^6-10^8 M_\odot$. Differentiating first galaxies at z>10 from SMDSs would be possible with spectroscopy: the SMDS (which are too cool produce significant nebular emission) will have only absorption lines while the galaxies are likely to produce emission lines as well. Of particular interest would be the 1640 HeII emission line as well as Hα lines which would be signatures of early galaxies rather than SMDSs. The detection of SMDSs would not only provide alternative evidence for WIMPs but would also provide possible seeds for the formation of supermassive black holes that power QSOs at z~6.

preprint2011arXiv

Probing the shape and history of the Milky Way halo with orbital spectral analysis

The phase space coordinates of individual halo stars obtained by Galactic surveys enable the computation of their full 3-dimensional orbits. Spectral analysis of halo orbits can be used to construct "frequency maps" which provide a compact representation of the 6-dimensional phase space distribution function. Frequency maps identify important major orbit families, and the orbital abundances reflect the shape and orientation of the dark matter halo relative to the disk. We apply spectral analysis to halo orbits in a series of controlled simulations of disk galaxies. Although the shape of the simulated halo varies with radius, frequency maps of local samples of halo orbits confined to the inner halo contain most of the information about the global shape of the halo and its major orbit families. Quiescent or adiabatic disk formation results in significant trapping of halo orbits in resonant orbit families (i.e. orbits with commensurable frequencies). If a good estimate of the Galactic potential in the inner halo (within ~50kpc) is available, the appearance of strong, stable resonances in frequency maps of halo orbits will allow us to determine the degree of resonant trapping induced by the disk potential. The locations and strengths of these resonant families are determined both by the global shape of the halo and its distribution function. Identification of such resonances in the Milky Way's stellar halo would therefore provide evidence of an extended period of adiabatic disk growth. If the Galactic potential is not known exactly, a measure of the diffusion rate of large sample of 10^4 halo orbits can help distinguish between the true potential and an incorrect potential. The orbital spectral analysis methods described in this paper provide a strong complementarity to existing methods for constraining the potential of the Milky Way halo and its stellar distribution function (ABRIDGED).

preprint2011arXiv

The Orbital Structure of Triaxial Galaxies with Figure Rotation

We survey the properties of all orbit families in the rotating frame of a family of realistic triaxial potentials with central supermassive black holes (SMBHs). In such galaxies, most regular box orbits (vital for maintaining triaxiality) are associated with resonances which occupy two-dimensional surfaces in configuration space. For slow figure rotation all orbit families are largely stable. At intermediate pattern speeds a significant fraction of the resonant box orbits as well as inner long-axis tubes are destabilized by the "envelope doubling" that arises from the Coriolis forces and are driven into the destabilizing center. Thus, for pattern rotation periods .2 Gyr < Tp < 5 Gyr, the two orbit families that are most important for maintaining triaxiality are highly chaotic. As pattern speed increases there is also a sharp decrease in the overall fraction of prograde short-axis tubes and a corresponding increase in the retrograde variety. At the highest pattern speeds (close to that of triaxial bars), box-like orbits undergo a sudden transition to a new family of stable retrograde loop-like orbits, which resemble orbits in three-dimensional bars, and circulate about the short axis. Our analysis implies that triaxial systems (with central cusps and SMBHs) can either have high pattern speeds like fast bars or low patten speeds like triaxial elliptical galaxies or dark matter halos found in N-body simulations. Intermediate pattern speeds produce a high level of stochasticity in both the box and inner long-axis tube orbit families implying that stable triaxial systems are unlikely to have such pattern speeds.

preprint2011arXiv

Using frequency maps to constrain the distribution function of the Milky Way stellar halo

Resolved surveys of the Milky Way's stellar halo can obtain all 6 phase space coordinates of tens of thousands of individual stars, making it possible to compute their 3-dimensional orbits. Spectral analysis of large numbers of halo orbits can be used to construct frequency maps which are a compact, yet informative representation of their phase space distribution function (DF). Such maps can be used to infer the major types of orbit families that constitute the DF of stellar halo and their relative abundances. The structure of the frequency maps, especially the resonant orbits, reflects the formation history and shape of the dark matter potential and its orientation relative to the disk. The application of frequency analysis to cosmological hydrodynamic simulations of disk galaxies shows that the orbital families occupied by halo stars and dark matter particles are very similar, implying that stellar halo orbits can be used to constrain the DF of the dark matter halo, possibly impacting future direct dark matter detection experiments. An application of these methods to a sample of \sim 16,000 Milky Way halo and thick disk stars from the SDSS-SEGUE survey yields a frequency map with strong evidence for resonant trapping of halo stars by the Milky Way disk, in a manner predicted by controlled simulations in which the disk grows adiabatically. The application of frequency analysis methods to current and future phase space data for Milky Way halo stars will provide new insights into the formation history of the dierent components of the Galaxy and the DF of the halo.

preprint2010arXiv

New tools for probing the phase space structure of dark matter halos

We summarize recent developments in the use of spectral methods for analyzing large numbers of orbits in N-body simulations to obtain insights into the global phase space structure of dark matter halos. The fundamental frequencies of oscillation of orbits can be used to understand the physical mechanism by which the shapes of dark matter halos evolve in response to the growth of central baryonic components. Halos change shape primarily because individual orbits change their shapes adiabatically in response to the growth of a baryonic component, with those at small radii become preferentially rounder. Chaotic scattering of orbits occurs only when the central point mass is very compact and is equally effective for centrophobic long-axis tube orbits as it is for centrophilic box orbits.

preprint2010arXiv

Supermassive Dark Stars: Detectable in JWST

The first phase of stellar evolution in the history of the Universe may be Dark Stars, powered by dark matter heating rather than by nuclear fusion. Weakly Interacting Massive Particles, which may be their own antipartners, collect inside the first stars and annihilate to produce a heat source that can power the stars for millions to billions of years. In this paper we show that these objects can grow to be supermassive dark stars (SMDS) with masses $\gtrsim (10^5-10^7) \msun$. The growth continues as long as dark matter heating persists, since dark stars are large and cool (surface temperature $\lesssim 5\times 10^4$K) and do not emit enough ionizing photons to prevent further accretion of baryons onto the star. The dark matter may be provided by two mechanisms: (1) gravitational attraction of dark matter particles on a variety of orbits not previously considered, and (2) capture of WIMPs due to elastic scattering. Once the dark matter fuel is exhausted, the SMDS becomes a heavy main sequence star; these stars eventually collapse to form massive black holes that may provide seeds for supermassive black holes in the Universe. SMDS are very bright, with luminosities exceeding $(10^9-10^{11}) L_\odot$. We demonstrate that for several reasonable parameters, these objects will be detectable with JWST. Such an observational discovery would confirm the existence of a new phase of stellar evolution powered by dark matter.

preprint2009arXiv

Evolution of Dark Matter Phase-Space Density Distributions in Equal-Mass Halo Mergers

We use dissipationless N-body simulations to investigate the evolution of the true coarse-grained phase-space density distribution f(x,v) in equal-mass mergers between dark matter (DM) halos. The halo models are constructed with various asymptotic power-law indices ranging from steep cusps to core-like profiles and we employ the phase-space density estimator ``Enbid'' developed by Sharma & Steinmetz to compute f(x,v). The adopted force resolution allows robust phase-space density profile estimates in the inner ~1% of the virial radii of the simulated systems. We confirm that mergers result in a decrease of the coarse-grained phase-space density in accordance with expectations from Mixing Theorems for collisionless systems. We demonstrate that binary mergers between identical DM halos produce remnants that retain excellent memories of the inner slopes and overall shapes of the phase-space density distribution of their progenitors. The robustness of the phase-space density profiles holds for a range of orbital energies, and a variety of encounter configurations including sequences of several consecutive merger events, designed to mimic hierarchical merging, and collisions occurring at different cosmological epochs. If the progenitor halos are constructed with appreciably different asymptotic power-law indices, we find that the inner slope and overall shape of the phase-space density distribution of the remnant are substantially closer to that of the initial system with the steepest central density cusp. These results explicitly demonstrate that mixing is incomplete in equal-mass mergers between DM halos, as it does not erase memory of the progenitor properties. Our results also confirm the recent analytical predictions of Dehnen (2005) regarding the preservation of merging self-gravitating central density cusps.

preprint2009arXiv

Evolution of the Dark Matter Phase-Space Density Distributions of LCDM Halos

We study the evolution of phase-space density during the hierarchical structure formation of LCDM halos. We compute both a spherically-averaged surrogate for phase-space density (Q) and the coarse-grained distribution function f(x,v) for dark matter particles that lie within~2 virial radii of four Milky-Way-sized dark matter halos. The estimated f(x,v) spans over four decades at any radius. Dark matter particles that end up within two virial radii of a Milky-Way-sized DM halo at $z=0$ have an approximately Gaussian distribution in log(f) at early redshifts, but the distribution becomes increasingly skewed at lower redshifts. The value corresponding to the peak of the Gaussian decreases as the evolution progresses and is well described by a power-law in (1+z). The highest values of f are found at the centers of dark matter halos and subhalos, where f can be an order of magnitude higher than in the center of the main halo. The power-law Q(r) profile likely reflects the distribution of entropy (K = sigma^2/rho^{2/3} \propto r^{1.2}), which dark matter acquires as it is accreted onto a growing halo. The estimated f(x, v), on the other hand, exhibits a more complicated behavior. Although the median coarse-grained phase-space density profile F(r) can be approximated by a power-law in the inner regions of halos and at larger radii the profile flattens significantly. This is because phase-space density averaged on small scales is sensitive to the high-f material associated with surviving subhalos, as well as relatively unmixed material (probably in streams) resulting from disrupted subhalos, which contribute a sizable fraction of matter at large radii. (ABRIDGED)

preprint2008arXiv

Signatures of LCDM substructure in tidal debris

In the past decade, surveys of the stellar component of the Galaxy have revealed a number of streams from tidally disrupted dwarf galaxies and globular clusters. Simulations of hierarchical structure formation in LCDM cosmologies predict that the dark matter halo of a galaxy like the Milky Way contains hundreds of subhalos with masses of ~10^8 solar masses and greater, and it has been suggested that the existence of coherent tidal streams is incompatible with the expected abundance of substructure. We investigate the effects of dark matter substructure on tidal streams by simulating the disruption of a self-gravitating satellite on a wide range of orbits in different host models both with and without substructure. We find that the halo shape and the specific orbital path more strongly determine the overall degree of disruption of the satellite than does the presence or absence of substructure, i.e., the changes in the large-scale properties of the tidal debris due to substructure are small compared to variations in the debris from different orbits in a smooth potential. Substructure typically leads to an increase in the degree of clumpiness of the tidal debris in sky projection, and in some cases a more compact distribution in line-of-sight velocity. Substructure also leads to differences in the location of sections of debris compared to the results of the smooth halo model, which may have important implications for the interpretation of observed tidal streams. A unique signature of the presence of substructure in the halo which may be detectable by upcoming surveys is identified. We conclude, however, that predicted levels of substructure are consistent with a detection of a coherent tidal stream from a dwarf galaxy.

preprint2006arXiv

On relaxation processes in collisionless mergers

We analyze N-body simulations of halo mergers to investigate the mechanisms responsible for driving mixing in phase-space and the evolution to dynamical equilibrium. We focus on mixing in energy and angular momentum and show that mixing occurs in step-like fashion following pericenter passages of the halos. This makes mixing during a merger unlike other well known mixing processes such as phase mixing and chaotic mixing whose rates scale with local dynamical time. We conclude that the mixing process that drives the system to equilibrium is primarily a response to energy and angular momentum redistribution that occurs due to impulsive tidal shocking and dynamical friction rather than a result of chaotic mixing in a continuously changing potential. We also analyze the merger remnants to determine the degree of mixing at various radii by monitoring changes in radius, energy and angular momentum of particles. We confirm previous findings that show that the majority of particles retain strong memory of their original kinetic energies and angular momenta but do experience changes in their potential energies owing to the tidal shocks they experience during pericenter passages. Finally, we show that a significant fraction of mass (~ 40%) in the merger remnant lies outside its formal virial radius and that this matter is ejected roughly uniformly from all radii outside the inner regions. This highlights the fact that mass, in its standard virial definition, is not additive in mergers. We discuss the implications of these results for our understanding of relaxation in collisionless dynamical systems.

preprint2005arXiv

The Low End of the Supermassive Black Hole Mass Function: Constraining the Mass of a Nuclear Black Hole in NGC 205 via Stellar Kinematics

Hubble Space Telescope (HST) images and spectra of the nucleated dwarf elliptical galaxy NGC 205 are combined with 3-integral axisymmetric dynamical models to constrain the mass (M_BH) of a putative nuclear black hole. This is only the second attempt, after M33, to use resolved stellar kinematics to search for a nuclear black hole with mass below 10^6 solar masses. We are unable to identify a best-fit value of M_BH in NGC 205; however, the data impose a upper limit of 2.2x10^4 M_sun (1sigma confidence) and and upper limit of 3.8x10^4 M_sun (3sigma confidence). This upper limit is consistent with the extrapolation of the M_BH-sigma relation to the M_BH < 10^6 M_sunregime. If we assume that NGC 205 and M33 both contain nuclear black holes, the upper limits on M_BH in the two galaxies imply a slope of ~5.5 or greater for the M_BH-sigma relation. We use our 3-integral models to evaluate the relaxation time (T_r) and stellar collision time (T_coll) in NGC 205; T_r~10^8 yr or less in the nucleus and T_coll~10^11 yr. The low value of T_r is consistent with core collapse having already occurred, but we are unable to draw conclusions from nuclear morphology about the presence or absence of a massive black hole.

preprint2003arXiv

Difficulties with Recovering The Masses of Supermassive Black Holes from Stellar Kinematical Data

We investigate the ability of three-integral, axisymmetric, orbit-based modeling algorithms to recover the parameters defining the gravitational potential (M/L ratio and black hole mass Mh) in spheroidal stellar systems using stellar kinematical data. We show that the potential estimation problem is generically under-determined when applied to long-slit kinematical data of the kind used in most black hole mass determinations to date. A range of parameters (M/L, Mh) can provide equally good fits to the data, making it impossible to assign best-fit values. We illustrate the indeterminacy using a variety of data sets derived from realistic models as well as published observations of the galaxy M32. In the case of M32, our reanalysis demonstrates that data published prior to 2000 are equally consistent with Mh in the range 1.5x10^6-5x10^6 solar masses, with no preferred value in that range. While the HST/STIS data for this galaxy may overcome the degeneracy in Mh, HST data for most galaxies do not resolve the black hole's sphere of influence and in these galaxies the degree of degeneracy allowed by the data may be substantial. We investigate the effect on the degeneracy of enforcing smoothness (regularization) constraints. However we find no indication that the true potential can be recovered simply by enforcing smoothness. For a given smoothing level, all solutions in the minimum-chisquare valley exhibit similar levels of noise. These experiments affirm that the indeterminacy is real and not an artifact associated with non-smooth solutions. (Abridged)