Researcher profile

David Merritt

David Merritt contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
33works
0followers
6topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

33 published item(s)

preprint2022arXiv

MOND and Methodology

In Logik der Forschung and later works, Karl Popper proposed a set of methodological rules for scientists. Among these were requirements that theories evolve in the direction of increasing content, and that new theories should only be accepted if some of their novel predictions are experimentally confirmed. There are currently two, viable theories of cosmology: the standard cosmological model, and a theory due to Mordehai Milgrom called MOND. Both theories can point to successes and failures, but only MOND has repeatedly made novel predictions that were subsequently found to be correct. Standard-model cosmologists, by contrast, have almost always responded to new observations in a post-hoc manner, adjusting or augmenting their theory as needed to obtain correspondence with the facts. I argue that these methodological differences render a comparison of the two theories in terms of their "truthlikeness" or "verisimilitude" essentially impossible since the two groups of scientists achieve correspondence with the facts in often very different ways, and I suggest that a better guide to the theories' progress toward truth might be the methodologies themselves.

preprint2013arXiv

The S-Star Cluster at the Center of the Milky Way: On the nature of diffuse NIR emission in the inner tenth of a parsec

Sagittarius A*, the super-massive black hole at the center of the Milky Way, is surrounded by a small cluster of high velocity stars, known as the S-stars. We aim to constrain the amount and nature of stellar and dark mass associated with the cluster in the immediate vicinity of Sagittarius A*. We use near-infrared imaging to determine the $K_\mathrm{s}$-band luminosity function of the S-star cluster members, and the distribution of the diffuse background emission and the stellar number density counts around the central black hole. This allows us to determine the stellar light and mass contribution expected from the faint members of the cluster. We then use post-Newtonian N-body techniques to investigate the effect of stellar perturbations on the motion of S2, as a means of detecting the number and masses of the perturbers. We find that the stellar mass derived from the $K_\mathrm{s}$-band luminosity extrapolation is much smaller than the amount of mass that might be present considering the uncertainties in the orbital motion of the star S2. Also the amount of light from the fainter S-cluster members is below the amount of residual light at the position of the S-star cluster after removing the bright cluster members. If the distribution of stars and stellar remnants is strongly enough peaked near Sagittarius A*, observed changes in the orbital elements of S2 can be used to constrain both their masses and numbers. Based on simulations of the cluster of high velocity stars we find that at a wavelength of 2.2 $μ$m close to the confusion level for 8 m class telescopes blend stars will occur (preferentially near the position of Sagittarius A*) that last for typically 3 years before they dissolve due to proper motions.

preprint2011arXiv

Nuclear Star-Forming Ring of the Milky Way: Simulations

We present hydrodynamic simulations of gas clouds in the central kpc region of the Milky Way that is modeled with a three-dimensional bar potential. Our simulations consider realistic gas cooling and heating, star formation, and supernova feedback. A ring of dense gas clouds forms as a result of X1-X2 orbit transfer, and our potential model results in a ring radius of ~200 pc, which coincides with the extraordinary reservoir of dense molecular clouds in the inner bulge, the Central Molecular Zone (CMZ). The gas clouds accumulated in the CMZ can reach high enough densities to form stars, and with an appropriate choice of simulation parameters, we successfully reproduce the observed gas mass and the star formation rate (SFR) in the CMZ, ~2x10^7 Msun and ~0.1 Msun/yr. Star formation in our simulations takes place mostly in the outermost X2 orbits, and the SFR per unit surface area outside the CMZ is much lower. These facts suggest that the inner Galactic bulge may harbor a mild version of the nuclear star-forming rings seen in some external disk galaxies. Furthermore, from the relatively small size of the Milky Way's nuclear bulge, which is thought to be a result of sustained star formation in the CMZ, we infer that the Galactic inner bulge probably had a shallower density profile or stronger bar elongation in the past.

preprint2011arXiv

Stellar Dynamics of Extreme-Mass-Ratio Inspirals

Inspiral of compact stellar remnants into massive black holes (MBHs) is accompanied by the emission of gravitational waves at frequencies that are potentially detectable by space-based interferometers. Event rates computed from statistical (Fokker-Planck, Monte-Carlo) approaches span a wide range due to uncertaintities about the rate coefficients. Here we present results from direct integration of the post-Newtonian N-body equations of motion descrbing dense clusters of compact stars around Schwarzschild MBHs. These simulations embody an essentially exact (at the post-Newtonian level) treatment of the interplay between stellar dynamical relaxation, relativistic precession, and gravitational-wave energy loss. The rate of capture of stars by the MBH is found to be greatly reduced by relativistic precession, which limits the ability of torques from the stellar potential to change orbital angular momenta. Penetration of this "Schwarzschild barrier" does occasionally occur, resulting in capture of stars onto orbits that gradually inspiral due to gravitational wave emission; we discuss two mechanisms for barrier penetration and find evidence for both in the simulations. We derive an approximate formula for the capture rate, which predicts that captures would be strongly disfavored from orbits with semi-major axes below a certain value; this prediction, as well as the predicted rate, are verified in the N-body integrations. We discuss the implications of our results for the detection of extreme-mass-ratio inspirals from galactic nuclei with a range of physical properties.

preprint2010arXiv

Dynamical Models of the Galactic Center

The distribution of late-type (old) stars in the inner parsec of the Milky Way is very different than expected for a relaxed population around a supermassive black hole. Instead of a density cusp, there is a 0.5 pc core. This article discusses what sorts of dynamical models might explain this "conundrum of old age." A straightforward interpretation is that the nucleus is unrelaxed, and that the distribution of the old giants reflects the distribution of fainter stars and stellar remnants generally in the core. On the other hand, a density cusp could be present in the unobserved populations, and the deficit of bright giants could be a result of interactions with these objects. At the present time, no model is clearly preferred.

preprint2009arXiv

Evolution of Nuclear Star Clusters

Two-body relaxation times of nuclear star clusters are short enough that gravitational encounters should substantially affect their structure in 10 Gyr or less. In nuclear star clusters without massive black holes, dynamical evolution is a competition between core collapse, which causes densities to increase, and heat input from the surrounding galaxy, which causes densities to decrease. The maximum extent of a nucleus that can resist expansion is derived numerically for a wide range of initial conditions; observed nuclei are shown to be compact enough to resist expansion, although there may have been an earlier generation of low-density nuclei that were dissolved. An evolutionary model for NGC 205 is presented which suggests that the nucleus of this galaxy has already undergone core collapse. Adding a massive black hole to a nucleus inhibits core collapse, and nuclear star clusters with black holes always expand, due primarily to heat input from the galaxy and secondarily to heating from stellar disruptions. The expansion rate is smaller for larger black holes due to the smaller temperature difference between galaxy and nucleus when the black hole is large. The rate of stellar tidal disruptions and its variation with time are computed for a variety of initial models. The disruption rate generally decreases with time due to the evolving nuclear density, particularly in the faintest galaxies, assuming that scaling relations derived for luminous galaxies can be extended to low luminosities.

preprint2009arXiv

Hypercompact Stellar Systems Around Recoiling Supermassive Black Holes

A supermassive black hole ejected from the center of a galaxy by gravitational wave recoil carries a retinue of bound stars - a "hypercompact stellar system" (HCSS). The numbers and properties of HCSSs contain information about the merger histories of galaxies, the late evolution of binary black holes, and the distribution of gravitational-wave kicks. We relate the structural properties of HCSSs to the properties of their host galaxies, in two regimes: collisional, i.e. short nuclear relaxation times; and collisionless, i.e. long nuclear relaxtion times. HCSSs are expected to be similar in size and luminosity to globular clustersbut in extreme cases their stellar mass can approach that of UCDs. However they differ from all other classes of compact stellar system in having very high internal velocities. We show that the kick velocity is encoded in the velocity dispersion of the bound stars. Given a large enough sample of HCSSs, the distribution of gravitational-wave kicks can therefore be empirically determined. We combine a hierarchical merger algorithm with stellar population models to compute the rate of production of HCSSs over time and the probability of observing HCSSs in the local universe as a function of their apparent magnitude, color, size and velocity dispersion, under two assumptions about the star formation history prior to the kick. We predict that roughly 100 should be detectable within 2 Mpc of the center of the Virgo cluster and that many of these should be bright enough that their high internal velocity dispersions could be measured with reasonable exposure times.

preprint2009arXiv

Perturbations of Intermediate-mass Black Holes on Stellar Orbits in the Galactic Center

We study the short- and long-term effects of an intermediate mass black hole (IMBH) on the orbits of stars bound to the supermassive black hole (SMBH) at the center of the Milky Way. A regularized N-body code including post-Newtonian terms is used to carry out direct integrations of 19 stars in the S-star cluster for 10 Myr. The mass of the IMBH is assigned one of four values from 400 Msun to 4000 Msun, and its initial semi-major axis with respect to the SMBH is varied from 0.3-30 mpc, bracketing the radii at which inspiral of the IMBH is expected to stall. We consider two values for the eccentricity of the IMBH/SMBH binary, e=(0,0.7), and 12 values for the orientation of the binary's plane. Changes at the level of 1% in the orbital elements of the S-stars could occur in just a few years if the IMBH is sufficiently massive. On time scales of 1 Myr or longer, the IMBH efficiently randomizes the eccentricities and orbital inclinations of the S-stars. Kozai oscillations are observed when the IMBH lies well outside the orbits of the stars. Perturbations from the IMBH can eject stars from the cluster, producing hypervelocity stars, and can also scatter stars into the SMBH; stars with high initial eccentricities are most likely to be affected in both cases. The distribution of S-star orbital elements is significantly altered from its currently-observed form by IMBHs with masses greater than 1000 Msun if the IMBH/SMBH semi-major axis lies in the range 3-10 mpc. We use these results to further constrain the allowed parameters of an IMBH/SMBH binary at the Galactic center.

preprint2009arXiv

Testing Properties of the Galactic Center Black Hole Using Stellar Orbits

The spin and quadrupole moment of the supermassive black hole at the Galactic center can in principle be measured via astrometric monitoring of stars orbiting at milliparsec (mpc) distances, allowing tests of general relativistic "no-hair" theorems (Will 2008). One complicating factor is the presence of perturbations from other stars, which may induce orbital precession of the same order of magnitude as that due to general relativistic effects. The expected number of stars in this region is small enough that full N-body simulations can be carried out. We present the results of a comprehensive set of such simulations, which include a post-Newtonian treatment of spin-orbit effects. A number of possible models for the distribution of stars and stellar remnants are considered. We find that stellar perturbations are likely to obscure the signal due to frame-dragging for stars beyond ~0.5 mpc from the black hole, while measurement of the quadrupole moment is likely to require observation of stars inside ~0.2 mpc. A high fraction of stellar remnants, e.g. 10-Solar-mass black holes, in this region would make tests of GR problematic at all radii. We discuss the possibility of separating the effects of stellar perturbations from those due to GR.

preprint2008arXiv

An Instability in Triaxial Stellar Systems

The radial-orbit instability is a collective phenomenon that has heretofore only been observed in spherical systems. We find that this instability occurs also in triaxial systems, as we checked by performing extensive N-body simulations whose initial conditions were obtained by sampling a self-consistent triaxial model of a cuspy galaxy composed of luminous and dark matter. N-body simulations show a time evolution of the galaxy that is not due to the development of chaotic motions but, rather, to the collective instability induced by an excess of box-like orbits. The instability quickly transforms such models into a more prolate configuration, with 0.64<b/a<0.77 and 0.6<c/a<0.7 for the dark halo and 0.64<b/a<0.77 and 0.59<c/a<0.67 for the luminous matter. Stable triaxial, cuspy galaxies with dark matter halos are obtained when the contribution of radially-biased orbits to the solution is reduced. These results constitute the first evidence of the radial-orbit instability in triaxial galaxy models.

preprint2008arXiv

Binary Black Hole Merger in Galactic Nuclei: Post-Newtonian Simulations

This paper studies the formation and evolution of binary supermassive black holes (SMBHs) in rotating galactic nuclei, focusing on the role of stellar dynamics. We present the first N-body simulations that follow the evolution of the SMBHs from kiloparsec separations all the way to their final relativistic coalescence, and that can robustly be scaled to real galaxies. The N-body code includes post-Newtonian (PN) corrections to the binary equations of motion up to order 2.5; we show that the evolution of the massive binary is only correctly reproduced if the conservative 1PN and 2PN terms are included. The orbital eccentricities of the massive binaries in our simulations are often found to remain large until shortly before coalescence. This directly affects not only their orbital evolution rates, but has important consequences as well for the gravitational waveforms emitted during the relativistic inspiral. We estimate gravitational wave amplitudes when the frequencies fall inside the band of the (planned) Laser Interferometer Space Antennae (LISA). We find significant contributions -- well above the LISA sensitivity curve -- from the higher-order harmonics.

preprint2008arXiv

Merger of Massive Black Holes using N-Body Simulations with Post-Newtonian Corrections

We present preliminary results from self-consistent, high resolution direct {\it N}-body simulations of massive black hole binaries in mergers of galactic nuclei. The dynamics of the black hole binary includes the full Post-Newtonian corrections (up to 2.5PN) to its equations of motion. We show that massive black holes starting at separations of 100 pc can evolve down to gravitational-wave-induced coalescence in less than a Hubble time. The binaries, in our models, often form with very high eccentricity and, as a result, reach separations of 50 Schwarzschild radius with eccentricities which are clearly distinct from zero -- even though gravitational wave emission damps the eccentricity during the inspiral. These deviations from exact circular orbits, at such small separations, may have important consequences for LISA data analysis.

preprint2008arXiv

Tidal Disruption Flares from Recoiling Supermassive Black Holes

Supermassive black holes ejected from galaxy nuclei by gravitational wave recoil will carry a retinue of bound stars, even in the absence of an accretion disk. We discuss the observable signatures related to these stars, with an emphasis on electromagnetic flares from stars that are tidally disrupted by the black hole. We calculate disruption rates for the bound, and the unbound, stars. The rates are smaller than, but comparable to, rates for non-recoiling black holes. A key observational consequence is the existence of powerful off-nuclear and intergalactic X-ray flares. We also discuss other observable signatures associated with the bound stars, including episodic X-ray emission from accretion due to stellar mass loss; intergalactic supernovae; and feedback trails.

preprint2007arXiv

Consequences of Triaxiality for Gravitational Wave Recoil of black holes

Coalescing binary black holes experience a ``kick&#39;&#39; due to anisotropic emission of gravitational waves with an amplitude as great as 200$ km/s. We examine the orbital evolution of black holes that have been kicked from the centers of triaxial galaxies. Time scales for orbital decay are generally longer in triaxial galaxies than in equivalent spherical galaxies, since a kicked black hole does not return directly through the dense center where the dynamical friction force is highest. We evaluate this effect by constructing self-consistent triaxial models and integrating the trajectories of massive particles after they are ejected from the center; the dynamical friction force is computed directly from the velocity dispersion tensor of the self-consistent model. We find return times that are several times longer than in a spherical galaxy with the same radial density profile, particularly in galaxy models with dense centers, implying a substantially grea

preprint2007arXiv

Impact of astrophysical processes on the gamma-ray background from dark matter annihilations

We study the impact of astrophysical processes on the gamma-ray background produced by the annihilation of dark matter particles in cosmological halos, with particular attention to the consequences of the formation of supermassive black holes. In scenarios where these objects form adiabatically from the accretion of matter on small seeds, dark matter is first compressed into very dense ``spikes&#39;&#39;, then its density progressively decreases due to annihilations and scattering off of stellar cusps. With respect to previous analyses, based on non-evolving halos, the predicted annihilation signal is higher and significantly distorted at low energies, reflecting the large contribution to the total flux from unevolved spikes at high redshifts. The peculiar spectral feature arising from the specific redshift distribution of the signal, would discriminate the proposed scenario from more conventional astrophysical explanations. We discuss how this affects the prospects for detection and demonstrate that the gamma-ray background from DM annihilations might be detectable even in absence of a signal from the Galactic center.

preprint2007arXiv

Mass Deficits, Stalling Radii, and the Merger Histories of Elliptical Galaxies

A binary supermassive black hole leaves an imprint on a galactic nucleus in the form of a &#34;mass deficit,&#34; a decrease in the mass of the nucleus due to ejection of stars by the binary. The magnitude of the mass deficit is in principle related to the galaxy&#39;s merger history, but the relation has never been quantified. Here, high-accuracy N-body simulations are used to calibrate this relation. Mass deficits are shown to be approximately 0.5M_{12}, with M_{12} the total mass of the binary; the coefficient in this relation depends only weakly on the binary mass ratio or on the galaxy&#39;s pre-existing density profile. Hence, after N mergers, the mass deficit is ~0.5 N M_h with M_h the final (current) black hole mass. When compared with observed mass deficits, this result implies between 1 and 3 mergers for most galaxies, in accord with hierarchical galaxy formation models. Implications for binary stalling radii, the origin of hyper-velocity stars, and the distribution of dark matter at the centers of galaxies are discussed.

preprint2007arXiv

Maximum gravitational recoil

Recent calculations of gravitational radiation recoil generated during black-hole binary mergers have reopened the possibility that a merged binary can be ejected even from the nucleus of a massive host galaxy. Here we report the first systematic study of gravitational recoil of equal-mass binaries with equal, but anti-aligned, spins parallel to the orbital plane. Such an orientation of the spins is expected to maximize the recoil. We find that recoil velocity (which is perpendicular to the orbital plane) varies sinusoidally with the angle that the initial spin directions make with the initial linear momenta of each hole and scales up to a maximum of ~4000 km/s for maximally-rotating holes. Our results show that the amplitude of the recoil velocity can depend sensitively on spin orientations of the black holes prior to merger.

preprint2007arXiv

The ACS Virgo Cluster Survey. XIII. SBF Distance Catalog and the Three-Dimensional Structure of the Virgo Cluster

The ACS Virgo Cluster Survey consists of HST ACS imaging for 100 early-type galaxies in the Virgo Cluster, observed in the F475W and F850LP filters. We derive distances for 84 of these galaxies using the method of surface brightness fluctuations (SBFs), present the SBF distance catalog, and use this database to examine the three-dimensional distribution of early-type galaxies in the Virgo Cluster. The SBF distance moduli have a mean (random) measurement error of 0.07 mag (0.5 Mpc), or roughly 3 times better than previous SBF measurements for Virgo Cluster galaxies. Five galaxies lie at a distance of ~23 Mpc and are members of the W&#39; cloud. The remaining 79 galaxies have a narrow distribution around our adopted mean distance of 16.5+/-0.1 (random mean error) +/-1.1 Mpc (systematic). The rms distance scatter of this sample is 0.6+/-0.1 Mpc, with little dependence on morphological type or luminosity class (i.e., 0.7+/-0.1 and 0.5+/-0.1 Mpc for the giants and dwarfs, respectively). The back-to-front depth of the cluster measured from our sample of galaxies is 2.4+/-0.4 Mpc (i.e., +/-2sigma of the intrinsic distance distribution). The M87 (cluster A) and M49 (cluster B) subclusters are found to lie at distances of 16.7+/-0.2 and 16.4+/-0.2 Mpc, respectively. There may be a third subcluster associated with M86. A weak correlation between velocity and line-of-sight distance may be a faint echo of the cluster velocity distribution not having yet completely virialized. In three dimensions, Virgo&#39;s early-type galaxies appear to define a slightly triaxial distribution, with axis ratios of (1:0.7:0.5). The principal axis of the best-fit ellipsoid is inclined ~20-40 deg. from the line of sight, while the galaxies belonging to the W&#39; cloud lie on an axis inclined by ~10-15 deg.

preprint2006arXiv

A Fundamental Relation Between Compact Stellar Nuclei, Supermassive Black Holes, and Their Host Galaxies

Imaging surveys with the Hubble Space Telescope (HST) have shown that roughly 50-80% of low- and intermediate-luminosity galaxies contain a compact stellar nucleus at their center, regardless of host galaxy morphological type. We combine HST imaging for early-type galaxies from the ACS Virgo Cluster Survey with ground-based long-slit spectra from KPNO to show that the masses of compact stellar nuclei in Virgo Cluster galaxies obey a tight correlation with the masses of the host galaxies. The same correlation is obeyed by the supermassive black holes (SBHs) found in predominantly massive galaxies. The compact stellar nuclei in the Local Group galaxies M33 and NGC205 are also found to fall along this same scaling relation. These results indicate that a generic by-product of galaxy formation is the creation of a Central Massive Object (CMO) -- either a SBH or a compact stellar nucleus -- that contains a mean fraction, ~ 0.2%, of the total galactic mass. In galaxies with masses greater than a few times 10^10 solar masses, SBHs appear to be the dominant mode of CMO formation.

preprint2006arXiv

Black Holes and Nuclear Dynamics

Supermassive black holes inhabit galactic nuclei, and their presence influences in crucial ways the evolution of the stellar distribution. The low-density cores observed in bright galaxies are probably a result of black hole infall, while steep density cusps like those at the Galactic center are a result of energy exchange between stars moving in the gravitational field of the single black hole. Loss-cone dynamics are substantially more complex in galactic nuclei than in collisionally-relaxed systems like globular clusters due to the wider variety of possible geometries and orbital populations. The rate of star-black hole interactions has begun to be constrained through observations of energetic events associated with stellar tidal disruptions.

preprint2006arXiv

Dynamical Cusp Regeneration

After being destroyed by a binary supermassive black hole, a stellar density cusp can regrow at the center of a galaxy via energy exchange between stars moving in the gravitational field of the single, coalesced hole. We illustrate this process via high-accuracy N-body simulations. Regeneration requires roughly one relaxation time and the new cusp extends to a distance of roughly one-fifth the black hole&#39;s influence radius, with density rho ~ r^{-7/4}; the mass in the cusp is of order 10% the mass of the black hole. Growth of the cusp is preceded by a stage in which the stellar velocity dispersion evolves toward isotropy and away from the tangentially-anisotropic state induced by the binary. We show that density profiles similar to those observed at the center of the Milky Way and M32 can regenerate themselves in several Gyr following infall of a second black hole; the presence of density cusps at the centers of these galaxies can therefore not be used to infer that no merger has occurred. We argue that Bahcall-Wolf cusps are ubiquitous in stellar spheroids fainter than M_V ~ -18.5 that contain supermassive black holes, but the cusps have not been detected outside of the Local Group since their angular sizes are less than 0.1&#34;. We show that the presence of a cusp implies a lower limit of \~10^{-4} per year on the rate of stellar tidal disruptions, and discuss the consequences of the cusps for gravitational lensing and the distribution of dark matter on sub-parsec scales.

preprint2006arXiv

Efficient Merger of Binary Supermassive Black Holes in Non-Axisymmetric Galaxies

Binary supermassive black holes form naturally in galaxy mergers, but their long-term evolution is uncertain. In spherical galaxies, N-body simulations show that binary evolution stalls at separations much too large for significant emission of gravitational waves (the &#34;final parsec problem&#34;). Here, we follow the long-term evolution of a massive binary in more realistic, triaxial and rotating galaxy models. We find that the binary does not stall. The binary hardening rates that we observe are sufficient to allow complete coalescence of binary SBHs in 10 Gyr or less, even in the absence of collisional loss-cone refilling or gas-dynamical torques, thus providing a potential solution to the final parsec problem.

preprint2006arXiv

Empirical models for Dark Matter Halos. I. Nonparametric Construction of Density Profiles and Comparison with Parametric Models

We use techniques from nonparametric function estimation theory to extract the density profiles, and their derivatives, from a set of N-body dark matter halos. We consider halos generated from LCDM simulations of gravitational clustering, as well as isolated, spherical collapses. The logarithmic density slopes gamma = d(log rho)/d(log r) of the LCDM halos are found to vary as power-laws in radius, reaching values of gamma ~ -1 at the innermost resolved radii (~0.01 r_virial). This behavior is significantly different from that of broken power-law models like the NFW profile, but similar to that of models like de Vaucouleurs&#39;. Accordingly, we compare the N-body density profiles with various parametric models to find which provide the best fit. We consider an NFW-like model with arbitrary inner slope; Dehnen & McLaughlin&#39;s anisotropic model; Einasto&#39;s model (identical in functional form to Sersic&#39;s model but fit to the space density); and the density model of Prugniel & Simien that was designed to match the deprojected form of Sersic&#39;s R^{1/n} law. Overall, the best-fitting model to the LCDM halos is Einasto&#39;s, although the Prugniel-Simien and Dehnen-McLaughlin models also perform well. With regard to the spherical collapse halos, both the Prugniel-Simien and Einasto models describe the density profiles well, with an rms scatter some four times smaller than that obtained with either the NFW-like model or the 3-parameter Dehnen-McLaughlin model. Finally, we confirm recent claims of a systematic variation in profile shape with halo mass.

preprint2006arXiv

Empirical Models for Dark Matter Halos. II. Inner profile slopes, dynamical profiles, and rho/sigma^3

We have recently shown that both the Prugniel-Simien model and Sersic&#39;s function (hereafter referred to as the Einasto model when applied to internal density profiles) describe simulated dark matter halos better than an NFW-like model with equal number of parameters. Here we provide analytical expressions for the logarithmic slopes of these models, and compare them with data from real galaxies. Depending on the Einasto parameters of the dark matter halo, one can expect an extrapolated, inner (0.01--1 kpc), logarithmic profile slope ranging from -0.2 to -1.5, with a typical value at 0.1 kpc around -0.7. Application of this (better fitting) model therefore alleviates some of the past disagreement with observations on this issue. We additionally provide useful expressions for the concentration and assorted scale radii: r_s, r_{-2}, r_e, R_e, r_virial, and r_max -- the radius where the circular velocity profile has its maximum value. We also present the circular velocity profiles and the radial behavior of rho(r)/sigma(r)^3 for both the Einasto and Prugniel-Simien models. We find this representation of the phase-space density profile to be well approximated by a power-law with slope slightly shallower than -2 near r=r_{-2}.

preprint2006arXiv

Empirical Models for Dark Matter Halos. III. The Kormendy relation and the log(rho_e)-log(R_e) relation

We have recently shown that the 3-parameter density-profile model from Prugniel & Simien provides a better fit to simulated, galaxy- and cluster-sized, dark matter halos than an NFW-like model with arbitrary inner profile slope gamma (Paper I). By construction, the parameters of the Prugniel-Simien model equate to those of the Sersic R^{1/n} function fitted to the projected distribution. Using the Prugniel-Simien model, we are therefore able to show that the location of simulated (10^{12} M_sun) galaxy-sized dark matter halos in the <mu>_e-log(R_e) diagram coincides with that of brightest cluster galaxies, i.e., the dark matter halos appear consistent with the Kormendy relation defined by luminous elliptical galaxies. These objects are also seen to define the new, and equally strong, relation log(rho_e) = 0.5 - 2.5log(R_e), in which rho_e is the internal density at r=R_e. Simulated (10^{14.5} M_sun) cluster-sized dark matter halos and the gas component of real galaxy clusters follow the relation log(rho_e) = 2.5[1 - log(R_e)]. Given the shapes of the various density profiles, we are able to conclude that while dwarf elliptical galaxies and galaxy clusters can have dark matter halos with effective radii of comparable size to the effective radii of their baryonic component, luminous elliptical galaxies can not. For increasingly large elliptical galaxies, with increasingly large profile shapes `n&#39;, to be dark matter dominated at large radii requires dark matter halos with increasingly large effective radii compared to the effective radii of their stellar component.

preprint2006arXiv

Model Atmospheres for Irradiated Giant Stars: Implications for the Galactic Center

Irradiation of a stellar atmosphere by an external source (e.g. an AGN) changes its structure and therefore its spectrum. Using a state-of-the-art stellar atmosphere code, we calculate the infrared spectra of such irradiated and transformed stars. We show that the original spectrum of the star, which is dominated by molecular bands, changes dramatically when irradiated even by a low-luminosity AGN ($L_{\rm X} = 10^{33}$ erg s$^{-1}$), becoming dominated by atomic lines in absorption. We study the changes in the spectrum of low-mass carbon- and oxygen-rich giant stars as they are irradiated by a modest AGN, similar to the one at the Galactic center (GC). The resulting spectra are similar to those of the faintest S-cluster stars observed in the GC. The spectrum of a star irradiated by a much brighter AGN, like that powered by a tidally disrupted star, is very different from that of any star currently observed near the GC. For the first time we have discovered that the structure of the atmosphere of an irradiated giant changes dramatically and induces a double inversion layer. We show that irradiation at the current level can explain the observed trend of CO band intensities decreasing as a function of increasing proximity to Sg $A^{*}$. This may indicate that (contrary to previous claims) there is no paucity of old giants in the GC, which coexist simultaneously with young massive stars.

preprint2006arXiv

Trends in the Globular Cluster Luminosity Function of Early-Type Galaxies

We present results from a study of the globular cluster luminosity function (GCLF) in a sample of 89 early-type galaxies observed as part of the ACS Virgo Cluster Survey. Using a Gaussian parametrization of the GCLF, we find a highly significant correlation between the GCLF dispersion, sigma, and the galaxy luminosity, M_B, in the sense that the GC systems in fainter galaxies have narrower luminosity functions. The GCLF dispersions in the Milky Way and M31 are fully consistent with this trend, implying that the correlation between sigma and galaxy luminosity is more fundamental than older suggestions that GCLF shape is a function of galaxy Hubble type. We show that the sigma - M_B relation results from a bonafide narrowing of the distribution of (logarithmic) cluster masses in fainter galaxies. We further show that this behavior is mirrored by a steepening of the GC mass function for relatively high masses, M >~ 3 x 10^5 M_sun, a mass regime in which the shape of the GCLF is not strongly affected by dynamical evolution over a Hubble time. We argue that this trend arises from variations in initial conditions and requires explanation by theories of cluster formation. Finally, we confirm that in bright galaxies, the GCLF &#34;turns over&#34; at the canonical mass scale of M_TO ~ 2 x 10^5 M_sun. However, we find that M_TO scatters to lower values (~1-2 x 10^5 M_sun) in galaxies fainter than M_B >~ -18.5, an important consideration if the GCLF is to be used as a distance indicator for dwarf ellipticals.

preprint2005arXiv

The ACS Virgo Cluster Survey. VII. Resolving the Connection Between Globular Clusters and Ultra-Compact Dwarf Galaxies

We investigate the connection between globular clusters and ultra-compact dwarf galaxies (UCDs) by examining the properties of several compact objects associated with M87, all of which were previously classified as globular clusters. Combining imaging from the Hubble Space Telescope with ground-based Keck spectroscopy, we find two objects to have half-light radii, velocity dispersions and mass-to-light ratios that are consistent with the predictions of population synthesis models for old, metal-rich, luminous globular clusters. Three other objects are much larger, with half-light radii of approximately 20pc, and have V-band mass-to-light ratios in the range 6-9. These objects, which we consider to be UCDs, resemble the nuclei of nucleated dwarf elliptical galaxies in Virgo, having similar mass-to-light ratios, luminosities and colors. These UCDs are found to obey the extrapolated scaling relations of galaxies more closely than those of Galactic globular clusters. There appears to be a transition between the two types of stellar systems at a mass of about two million solar masses. If the UCDs are gravitationally bound, then we suggest that the presence of dark matter is the fundamental property distinguishing globular clusters from UCDs. More than half of the UCD candidates uncovered in the ACS Virgo Cluster Survey are associated with a single galaxy -- M87 -- which suggests that proximity to the Virgo center may be of critical importance for the formation of these objects. These results show that distinguishing bonafide UCDs from bright globular clusters requires a careful analysis of their detailed structural and dynamical properties, particularly their mass-to-light ratios. (ABRIDGED)

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.

preprint2004arXiv

Consequences of gravitational radiation recoil

Coalescing binary black holes experience an impulsive kick due to anisotropic emission of gravitational waves. We discuss the dynamical consequences of the recoil accompanying massive black hole mergers. Recoil velocities are sufficient to eject most coalescing black holes from dwarf galaxies and globular clusters, which may explain the apparent absence of massive black holes in these systems. Ejection from giant elliptical galaxies would be rare, but coalescing black holes are displaced from the center and fall back on a time scale of order the half-mass crossing time. Displacement of the black holes transfers energy to the stars in the nucleus and can convert a steep density cusp into a core. Radiation recoil calls into question models that grow supermassive black holes from hierarchical mergers of stellar-mass precursors.

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&#39;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)

preprint2002arXiv

Tracing Black Hole Mergers Through Radio Lobe Morphology

Binary supermassive black holes are produced by galactic mergers as the black holes from the two galaxies fall to the center of the merged system and form a bound pair. The two black holes will eventually coalesce in an enormous burst of gravitational radiation. Here we show that the orientation of a black hole&#39;s spin axis would change dramatically even in a minor merger, leading to a sudden flip in the direction of any associated jet. We identify the winged or X-type radio sources with galaxies in which this has occurred. The implied coalescence rate is similar to the overall galaxy merger rate, suggesting that the prospects are good for observing gravitational waves from coalescing supermassive black holes.