Source author record

David Merritt

David Merritt 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

69works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

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

preprint2016arXiv

On the origins of enigmatic stellar populations in Local Group galactic nuclei

We consider the origins of enigmatic stellar populations in four Local Group galactic nuclei, specifically the Milky Way, M31, M32 and M33. These are centrally concentrated blue stars, found in three out of the four nuclear star clusters (NSCs) considered here. Their origins are unknown, but could include blue straggler (BS) stars, extended horizontal branch stars and young recently formed stars. Here, we calculate order-of-magnitude estimates for various collision rates, as a function of the host NSC environment and distance from the cluster centre. These rates are sufficiently high that BSs, formed via collisions between main sequence (MS) stars, could contribute non-negligibly ($\sim$ 1-10% in mass) to every surface brightness profile, with the exception of the Milky Way. Stellar evolution models show that the envelopes of red giant branch (RGB) stars must be nearly completely stripped to significantly affect their photometric appearance, which requires multiple collisions. Hence, the collision rates for individual RGB stars are only sufficiently high in the inner $\lesssim$ 0.1 pc of M31 and M32 for RGB destruction to occur. Collisions between white dwarfs and MS stars, which should ablate the stars, could offer a steady and significant supply of gas in every NSC in our sample. The gas could either fragment to form new stars, or accrete onto old MS stars already present. Thus, collisional processes could contribute significantly to the observed blue excesses in M31 and M33; future studies should be aimed at better constraining theoretical predictions to compliment existing and future observational data.

preprint2015arXiv

Gravitational Encounters and the Evolution of Galactic Nuclei. I. Method

An algorithm is described for evolving the phase-space density of stars or compact objects around a massive black hole at the center of a galaxy. The technique is based on numerical integration of the Fokker-Planck equation in energy-angular momentum space, f(E,L,t), and includes, for the first time, diffusion coefficients that describe the effects of both random and correlated encounters (resonant relaxation), as well as energy loss due to emission of gravitational waves. Destruction or loss of stars into the black hole are treated by means of a detailed boundary-layer analysis. Performance of the algorithm is illustrated by calculating two-dimensional, time-dependent and steady-state distribution functions and their corresponding loss rates.

preprint2015arXiv

Gravitational Encounters and the Evolution of Galactic Nuclei. II. Classical and Resonant Relaxation

Direct numerical integrations of the Fokker-Planck equation in energy-angular momentum space are carried out for stars orbiting a supermassive black hole (SBH) at the center of a galaxy. The algorithm, which was described in detail in an earlier paper, includes diffusion coefficients that describe the effects of both random ("classical") and correlated ("resonant") encounters. Steady-state solutions are similar to the Bahcall-Wolf solution but are modified at small radii due to the higher rate of diffusion in angular momentum, which results in a low-density core. The core radius is a few percent of the influence radius of the SBH. The corresponding phase-space density f(E,L) drops nearly to zero at low energies, implying almost no stars on tightly-bound orbits about the SBH. Steady-state rates of stellar disruption are presented, and a simple analytic expression is found that reproduces the numerical feeding rates with good accuracy. The distribution of periapsides of disrupted stars is also computed. Time-dependent solutions are also computed, starting from initial conditions similar to those produced by a binary SBH. In these models, feeding rates evolve on two timescales: rapid evolution during which the region evacuated by the massive binary is refilled by angular-momentum diffusion; and slower evolution as diffusion in energy causes the density profile at large radii to attain the Bahcall-Wolf form.

preprint2015arXiv

Gravitational Encounters and the Evolution of Galactic Nuclei. III. Anomalous Relaxation

This paper is the third in a series presenting the results of direct numerical integrations of the Fokker-Planck equation for stars orbiting a supermassive black hole (SBH) at the center of a galaxy. The algorithm of Paper II included diffusion coefficients that described the effects of random ("classical") and correlated ("resonant") relaxation. In this paper, the diffusion coefficients of Paper II have been generalized to account for the effects of "anomalous relaxation," the qualitatively different way in which eccentric orbits evolve in the regime of rapid relativistic precession. Two functional forms for the anomalous diffusion coefficients are investigated, based on power-law or exponential modifications of the resonant diffusion coefficients. The parameters defining the modified coefficients are first constrained by comparing the results of Fokker-Planck integrations with previously-published N-body integrations. Steady-state solutions are then obtained via the Fokker-Planck equation for models with properties similar to those of the Milky Way nucleus. Inclusion of anomalous relaxation leads to the formation of less prominent cores than in the case of resonant relaxation alone, due to the lengthening of diffusion timescales for eccentric orbits. Steady-state capture rates of stars by the SBH are found to always be less, by as much as an order of magnitude, than capture rates in the presence of resonant relaxation alone.

preprint2015arXiv

Gravitational Encounters and the Evolution of Galactic Nuclei. IV. Captures Mediated by Gravitational-Wave Energy Loss

Direct numerical integrations of the two-dimensional Fokker-Planck equation are carried out for compact objects orbiting a supermassive black hole (SBH) at the center of a galaxy. As in Papers I-III, the diffusion coefficients incorporate the effects of the lowest-order post-Newtonian corrections to the equations of motion. In addition, terms describing the loss of orbital energy and angular momentum due to the 5/2-order post-Newtonian terms are included. In the steady state, captures are found to occur in two regimes that are clearly differentiated in terms of energy, or semimajor axis; these two regimes are naturally characterized as "plunges" (low binding energy) and "EMRIs," or extreme-mass-ratio inspirals (high binding energy). The capture rate, and the distribution of orbital elements of the captured objects, are presented for two steady-state models based on the Milky Way: one with a relatively high density of remnants and one with a lower density. In both models, but particularly in the second, the steady-state energy distribution and the distribution of orbital elements of the captured objects are substantially different than if the Bahcall-Wolf energy distribution were assumed. The ability of classical relaxation to soften the blocking effects of the Schwarzschild barrier is quantified.These results, together with those of Papers I-III, suggest that a Fokker-Planck description can adequately represent the dynamics of collisional loss cones in the relativistic regime.

preprint2015arXiv

On the persistence of two small-scale problems in ΛCDM

We investigate the degree to which the inclusion of baryonic physics can overcome two long-standing problems of the standard cosmological model on galaxy scales: (i) the problem of satellite planes around Local Group galaxies, and (ii) the "too big to fail" problem. By comparing dissipational and dissipationless simulations, we find no indication that the addition of baryonic physics results in more flattened satellite distributions around Milky-Way-like systems. Recent claims to the contrary are shown to derive in part from a non-standard metric for the degree of flattening, which ignores the satellites' radial positions. If the full 3D positions of the satellite galaxies are considered, none of the simulations we analyse reproduce the observed flattening nor the observed degree of kinematic coherence of the Milky Way satellite system. Our results are consistent with the expectation that baryonic physics should have little or no influence on the structure of satellite systems on scales of hundreds of kiloparsecs. Claims that the "too big to fail" problem can be resolved by the addition of baryonic physics are also shown to be problematic.

preprint2015arXiv

The final-parsec problem in the collisionless limit

A binary supermassive black hole loses energy via ejection of stars in a galactic nucleus, until emission of gravitational waves becomes strong enough to induce rapid coalescence. Evolution via the gravitational slingshot requires that stars be continuously supplied to the binary, and it is known that in spherical galaxies the reservoir of such stars is quickly depleted, leading to stalling of the binary at parsec-scale separations. Recent N-body simulations of galaxy mergers and isolated nonspherical galaxies suggest that this stalling may not occur in less idealized systems. However, it remains unclear to what degree these conclusions are affected by collisional relaxation, which is much stronger in the numerical simulations than in real galaxies. In this study, we present a novel Monte Carlo method that can efficiently deal with both collisional and collisionless dynamics, and with galaxy models having arbitrary shapes. We show that without relaxation, the final-parsec problem may be overcome only in triaxial galaxies. Axisymmetry is not enough, but even a moderate departure from axisymmetry is sufficient to keep the binary shrinking. We find that the binary hardening rate is always substantially lower than the maximum possible, "full-loss-cone" rate, and that it decreases with time, but that stellar-dynamical interactions are nevertheless able to drive the binary to coalescence on a timescale <=1 Gyr in any triaxial galaxy.

preprint2014arXiv

An Eccentricity-Mass Relation for Galaxies from Tidally Disrupting Satellites

We infer the past orbit of the Sagittarius (Sgr) dwarf galaxy in the Milky Way halo by integrating backwards from its observed position and proper motions, including the effects of dynamical friction. Given measured proper motions, we show that there is a relation between the eccentricity ($e$) of Sgr's orbit and the mass of the Milky Way ($M_{T}$) in the limit of no dynamical friction. That relation can be fit by a power-law of the form: $e \approx 0.49 \left(M_{T}/10^{12} M_{\odot}\right)^{-0.88}$. At a fixed Milky Way mass, the dynamical friction term increases the mean eccentricity of the orbit and lowers the spread in eccentricities in proportion to the mass of the Sgr dwarf. We explore the implications of various observational constraints on Sgr's apocenter on the $e-M$ relation; Sgr masses outside the range $10^{9} M_{\odot} \la M_{\rm Sgr} \la 5 \times 10^{10} M_{\odot}$ are precluded, for Milky Way masses $\sim 1 - 2.5 \times 10^{12} M_{\odot}$. If Belokurov et al.'s (2014) observations represent the farthest point of Sgr's stream, then Milky Way masses in excess of $2 \times 10^{12} M_{\odot}$ are excluded for $M_{\rm Sgr} \la 10^{10} M_{\odot}$. Deeper observations of Sgr's tidal debris, from upcoming surveys such as GAIA, will allow better measurement of the Milky Way mass and of the Sgr dwarf.

preprint2014arXiv

Co-orbiting satellite galaxy structures are still in conflict with the distribution of primordial dwarf galaxies

Both major galaxies in the Local Group host planar distributions of co-orbiting satellite galaxies, the Vast Polar Structure (VPOS) of the Milky Way and the Great Plane of Andromeda (GPoA). The $Λ$CDM cosmological model did not predict these features. However, according to three recent studies the properties of the GPoA and the flattening of the VPOS are common features among sub-halo based $Λ$CDM satellite systems, and the GPoA can be naturally explained by satellites being acquired along cold gas streams. We point out some methodological issues in these studies: either the selection of model satellites is different from that of the observed ones, or an incomplete set of observational constraints has been considered, or the observed satellite distribution is inconsistent with basic assumptions. Once these issues have been addressed, the conclusions are different: features like the VPOS and GPoA are very rare (each with probability $\lesssim 10^{-3}$, and combined probability $< 10^{-5}$) if satellites are selected from a $Λ$CDM simulation combined with semi-analytic modelling, and accretion along cold streams is no natural explanation of the GPoA. The origin of planar dwarf galaxy structures remains unexplained in the standard paradigm of galaxy formation.

preprint2014arXiv

Relativistic dynamics of stars near a supermassive black hole

General relativistic precession limits the ability of gravitational encounters to increase the eccentricity $e$ of orbits near a supermassive black hole (SBH). This "Schwarzschild barrier" (SB) has been shown to play an important role in the orbital evolution of stars like the galactic center S-stars. However, the evolution of orbits below the SB, $e>e_\mathrm{SB}$, is not well understood; the main current limitation is the computational complexity of detailed simulations. Here we present an $N$-body algorithm that allows us to efficiently integrate orbits of test stars around a SBH including general relativistic corrections to the equations of motion and interactions with a large ($\gtrsim 10^3$) number of field stars. We apply our algorithm to the S-stars and extract diffusion coefficients describing the evolution in angular momentum $L$. We identify three angular momentum regimes, in which the diffusion coefficients depend in functionally different ways on $L$. Regimes of lowest and highest $L$ are well-described in terms of non-resonant relaxation (NRR) and resonant relaxation (RR), respectively. In addition, we find a new regime of "anomalous relaxation" (AR). We present analytic expressions, in terms of physical parameters, that describe the diffusion coefficients in all three regimes, and propose a new, empirical criterion for the location of the SB in terms of the $L$-dependence of the diffusion coefficients. Subsequently we apply our results to obtain the steady-state distribution of angular momentum for orbits near a SBH.

preprint2014arXiv

The final-parsec problem in non-spherical galaxies revisited

We consider the evolution of supermassive black hole binaries at the center of spherical, axisymmetric, and triaxial galaxies, using direct N-body integrations as well as analytic estimates. We find that the rates of binary hardening exhibit a significant N-dependence in all the models, at least for N in the investigated range of 10^5<=N<=10^6. Binary hardening rates are also substantially lower than would be expected if the binary loss cone remained full, as it would be if the orbits supplying stars to the binary were being efficiently replenished. The difference in binary hardening rates between the spherical and nonspherical models is less than a factor of two even in the simulations with the largest N. By studying the orbital populations of our models, we conclude that the rate of supply of stars to the binary via draining of centrophilic orbits is indeed expected to be much lower than the full-loss-cone rate, consistent with our simulations. We argue that the binary's evolution in the simulations is driven in roughly equal amounts by collisional and collisionless effects, even at the highest N-values currently accessible. While binary hardening rates would probably reach a limiting value for large N, our results suggest that we cannot approach that rate with currently available algorithms and computing hardware. The extrapolation of results from N-body simulations to real galaxies is therefore not straightforward, casting doubt on recent claims that triaxiality or axisymmetry alone are capable of solving the final-parsec problem in gas-free galaxies.

preprint2013arXiv

Loss Cone Dynamics

Supermassive black holes can capture or disrupt stars that come sufficiently close. This article reviews the dynamical processes by which stars or stellar remnants are placed onto loss-cone orbits and the implications for feeding rates. The capture rate is well defined for spherical galaxies with nuclear relaxation times that are shorter than the galaxy's age. However, even the dense nucleus of the Milky Way may be less than one relaxation time old, and this is certainly the case for more massive galaxies; the capture rate in such galaxies is an initial-value problem with poorly-known initial conditions and the rate can be much higher, or much lower, than the rate in a collisionally relaxed nucleus. In nonspherical (axisymmetric, triaxial) galaxies, torquing of orbits by the mean field can dominate perturbations due to random encounters, leading to much higher capture rates than in the spherical geometry, particularly in (massive) galaxies with long central relaxation times. Relativistic precession plays a crucial role in mediating the capture of compact remnants from regions very near to the black hole, by destroying the orbital correlations that would otherwise dominate the torques. The complex dynamics of relativistic loss cones are not yet well enough understood for accurate estimates of compact-object (EMRI) capture rates to be made.

preprint2013arXiv

Relativity and the evolution of the Galactic center S-star orbits

We consider the orbital evolution of the S-stars, the young main-sequence stars near the supermassive black hole (SBH) at the Galactic center (GC), and put constraints on competing models for their origin. Our analysis includes for the first time the joint effects of Newtonian and relativistic perturbations to the motion, including the dragging of inertial frames by a spinning SBH as well as torques due to finite-N asymmetries in the field-star distribution (resonant relaxation, RR). The evolution of the S-star orbits is strongly influenced by the Schwarzschild barrier (SB), the locus in the (E,L) plane where RR is ineffective at driving orbits to higher eccentricities. Formation models that invoke tidal disruption of binary stars by the SBH tend to place stars below (i.e., at higher eccentricities than) the SB; some stars remain below the barrier, but most stars are able to penetrate it, after which they are subject to RR and achieve a nearly thermal distribution of eccentricities. This process requires roughly 50 Myr in nuclear models with relaxed stellar cusps, or >~10 Myr, regardless of the initial distribution of eccentricities, in nuclear models that include a dense cluster of 10 M_Sun black holes. We find a probability of <~1% for any S-star to be tidally disrupted by the SBH over its lifetime.

preprint2013arXiv

The loss cone problem in axisymmetric nuclei

We consider the problem of consumption of stars by a supermassive black hole (SBH) at the center of an axisymmetric galaxy. Inside the SBH sphere of influence, motion of stars in the mean field is regular and can be described analytically in terms of three integrals of motion: the energy E, the z-component of angular momentum L_z, and the secular Hamiltonian H. There exist two classes of orbits, tubes and saucers; saucers occupy the low-angular-momentum parts of phase space and their fraction is proportional to the degree of flattening of the nucleus. Perturbations due to gravitational encounters lead to diffusion of stars in integral space, which can be described using the Fokker-Planck equation. We calculate the diffusion coefficients and solve this equation in the two-dimensional phase space (L_z, H), for various values of the capture radius and the degree of flattening. Capture rates are found to be modestly higher than in the spherical case, up to a factor of a few, and most captures take place from saucer orbits. We also carry out a set of collisional N-body simulations to confirm the predictions of the Fokker-Planck models. We discuss the implications of our results for rates of tidal disruption and capture in the Milky Way and external galaxies.

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.

preprint2012arXiv

Dissipationless Formation and Evolution of the Milky Way Nuclear Star Cluster

Abridged: In one widely discussed model for the formation of nuclear star clusters (NSCs), massive globular clusters spiral into the center of a galaxy and merge to form the nucleus. It is now known that at least some NSCs coexist with supermassive black holes (SBHs); this is the case, for instance, in the Milky Way (MW). In this paper, we investigate how the presence of a SMBH at the center of the MW impacts the merger hypothesis for the formation of its NSC. Starting from a model consisting of a low-density nuclear stellar disk and the SMBH, we use N-body simulations to follow the successive inspiral and merger of globular clusters. The clusters are started on circular orbits of radius 20 pc, and their initial masses and radii are set up in such a way as to be consistent with the galactic tidal field at that radius. The total accumulated mass by ~10 clusters is about 1.5x10^7 Solar masses. Each cluster is disrupted by the SMBH at a distance of roughly one parsec. The density profile that results after the final inspiral event is characterized by a core of roughly this radius, and an envelope with density that falls off as 1/r^2. These properties are similar to those of the MW NSC, with the exception of the core size, which in the MW is a little smaller. But by continuing the evolution of the model after the final inspiral event, we find that the core shrinks substantially via gravitational encounters in a time (when scaled to the MW) of 10 Gyr as the stellar distribution evolves toward a Bahcall-Wolf cusp. We also show that the luminosity function of the MW NSC is consistent with the hypothesis that a large fraction of the mass comes from (~10Gyr) old stars, brought in by globular clusters. We conclude that a model in which a large fraction of the mass of the MW NSC arose from infalling globular clusters is consistent with existing observational constraints.

preprint2012arXiv

Spin Evolution of Supermassive Black Holes and Galactic Nuclei

The spin angular momentum S of a supermassive black hole (SBH) precesses due to torques from orbiting stars, and the stellar orbits precess due to dragging of inertial frames by the spinning hole. We solve the coupled post-Newtonian equations describing the joint evolution of S and the stellar angular momenta Lj, j = 1...N in spherical, rotating nuclear star clusters. In the absence of gravitational interactions between the stars, two evolutionary modes are found: (1) nearly uniform precession of S about the total angular momentum vector of the system; (2) damped precession, leading, in less than one precessional period, to alignment of S with the angular momentum of the rotating cluster. Beyond a certain distance from the SBH, the time scale for angular momentum changes due to gravitational encounters between the stars is shorter than spin-orbit precession times. We present a model, based on the Ornstein-Uhlenbeck equation, for the stochastic evolution of star clusters due to gravitational encounters and use it to evaluate the evolution of S in nuclei where changes in the Lj are due to frame dragging close to the SBH and to encounters farther out. Long-term evolution in this case is well described as uniform precession of the SBH about the cluster's rotational axis, with an increasingly important stochastic contribution when SBH masses are small. Spin precessional periods are predicted to be strongly dependent on nuclear properties, but typical values are 10-100 Myr for low-mass SBHs in dense nuclei, 100 Myr - 10 Gyr for intermediate mass SBHs, and > 10 Gyr for the most massive SBHs. We compare the evolution of SBH spins in stellar nuclei to the case of torquing by an inclined, gaseous accretion disk.

preprint2011arXiv

Dynamical Friction around Supermassive Black Holes

The density of stars in galactic bulges is often observed to be flat or slowly rising inside the influence radius of the supermassive black hole (SMBH). Attributing the dynamical friction force to stars moving more slowly than the test body, as is commonly done, is likely to be a poor approximation in such a core since there are no stars moving more slowly than the local circular velocity. We have tested this prediction using large-scale N-body experiments. The rate of orbital decay never drops precisely to zero, because stars moving faster than the test body also contribute to the frictional force. When the contribution from the fast-moving stars is included in the expression for the dynamical friction force, and the changes induced by the massive body on the stellar distribution are taken into account, Chandrasekhar's theory is found to reproduce the rate of orbital decay remarkably well. However, this rate is still substantially smaller than the rate predicted by Chandrasekhar's formula in its most widely-used forms, implying longer time scales for inspiral. Motivated by recent observations that suggest a parsec-scale core around the Galactic center SMBH, we investigate the evolution of a population of stellar-mass black holes (BHs) as they spiral in to the center of the Galaxy. After ~10Gyr, we find that the density of BHs can remain substantially less than the density in stars at all radii; we conclude that it would be unjustified to assume that the spatial distribution of BHs at the Galactic center is well described by steady-state models. One consequence is that rates of capture of BHs by the SMBH at the Galactic center (EMRIs) may be much lower than in standard models. We finally study the orbital decay of satellite galaxies into the central region of giant ellipticals and discuss the formation of multiple nuclei and multiplet of black holes in such systems.

preprint2011arXiv

Efficient Merger of Binary Supermassive Black Holes in Merging Galaxies

In spherical galaxies, binary supermassive black holes (SMBHs) have difficulty reaching sub-parsec separations due to depletion of stars on orbits that intersect the massive binary - the final-parsec problem. Galaxies that form via major mergers are substantially nonspherical, and it has been argued that the centrophilic orbits in triaxial galaxies might provide stars to the massive binary at a high enough rate to avoid stalling. Here we test that idea by carrying out fully self-consistent merger simulations of galaxies containing central SMBHs. We find hardening rates of the massive binaries that are indeed much higher than in spherical models, and essentially independent of the number of particles used in the simulations. Binary eccentricities remain high throughout the simulations. Our results constitute a fully stellar-dynamical solution to the final-parsec problem and imply a potentially high rate of events for low-frequency gravitational wave detectors like LISA.

preprint2011arXiv

Long-term evolution of massive black hole binaries. IV. Mergers of galaxies with collisionally relaxed nuclei

We simulate mergers between galaxies containing collisionally-relaxed nuclei around massive black holes (MBHs). Our galaxies contain four mass groups, representative of old stellar populations; a primary goal is to understand the distribution of stellar-mass black holes (BHs) after the merger. Mergers are followed using direct-summation N-body simulations, assuming a mass ratio of 1:3 and two different orbits. Evolution of the binary MBH is followed until its separation has shrunk by a factor of 20 below the hard-binary separation. During the galaxy merger, large cores are carved out in the stellar distribution, with radii several times the influence radius of the massive binary. Much of the pre-existing mass segregation is erased during this phase. We follow the evolution of the merged galaxies for approximately three, central relaxation times after coalescence of the massive binary; both standard, and top-heavy, mass functions are considered. The cores that were formed in the stellar distribution persist, and the distribution of the stellar-mass black holes evolves against this essentially fixed background. Even after one central relaxation time, these models look very different from the relaxed, multi-mass models that are often assumed to describe the distribution of stars and stellar remnants near a massive BH. While the stellar BHs do form a cusp on roughly a relaxation time-scale, the BH density can be much smaller than in those models. We discuss the implications of our results for the EMRI problem and for the existence of Bahcall-Wolf cusps.

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.

preprint2011arXiv

The HST/ACS Coma Cluster Survey IV. Intergalactic Globular Clusters and the Massive Globular Cluster System at the Core of the Coma Galaxy Cluster

Intracluster stellar populations are a natural result of tidal interactions in galaxy clusters. Measuring these populations is difficult, but important for understanding the assembly of the most massive galaxies. The Coma cluster is one of the nearest truly massive galaxy clusters, and is host to a correspondingly large system of globular clusters (GCs). We use imaging from the HST/ACS Coma Cluster Survey to present the first definitive detection of a large population of intracluster GCs (IGCs) that fills the Coma cluster core and is not associated with individual galaxies. The GC surface density profile around the central massive elliptical galaxy, NGC 4874, is dominated at large radii by a population of IGCs that extend to the limit of our data (R<520 kpc). We estimate that there are 47000+/-1600 (random) +4000/-5000 (systematic) IGCs out to this radius, and that they make up ~70% of the central GC system, making this the largest GC system in the nearby Universe. Even including the GC systems of other cluster galaxies, IGCs still make up ~30-45% of the GCs in the cluster core. Observational limits from previous studies of the intracluster light (ICL) suggest that the IGC population has a high specific frequency. If the IGC population has a specific frequency similar to high-S_N dwarf galaxies, then the ICL has a total stellar mass of ~10^12 M_sun within the cluster core. The ICL makes up approximately half of the stellar luminosity and one-third of the stellar mass of the central (NGC4874+ICL) system. The color distribution of the IGC population is bimodal, with blue, metal-poor GCs outnumbering red, metal-rich GCs by a ratio of 4:1. The fraction of red IGCs (20%), and the red color of those GCs, implies that IGCs can originate from the halos of relatively massive, L* galaxies, and not solely from the disruption of dwarf galaxies. (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

Tidal breakup of binary stars at the Galactic Center. II. Hydrodynamic simulations

In Paper I, we followed the evolution of binary stars as they orbited near the supermassive black hole (SMBH) at the Galactic center, noting the cases in which the two stars would come close enough together to collide. In this paper we replace the point-mass stars by fluid realizations, and use a smoothed-particle hydrodynamics (SPH) code to follow the close interactions. We model the binary components as main-sequence stars with initial masses of 1, 3 and 6 Solar masses, and with chemical composition profiles taken from stellar evolution codes. Outcomes of the close interactions include mergers, collisions that leave both stars intact, and ejection of one star at high velocity accompanied by capture of the other star into a tight orbit around the SMBH. For the first time, we follow the evolution of the collision products for many ($\gtrsim 100$) orbits around the SMBH. Stars that are initially too small to be tidally disrupted by the SMBH can be puffed up by close encounters or collisions, with the result that tidal stripping occurs in subsequent periapse passages. In these cases, mass loss occurs episodically, sometimes for hundreds of orbits before the star is completely disrupted. Repeated tidal flares, of either increasing or decreasing intensity, are a predicted consequence. In collisions involving a low-mass and a high-mass star, the merger product acquires a high core hydrogen abundance from the smaller star, effectively resetting the nuclear evolution "clock" to a younger age. Elements like Li, Be and B that can exist only in the outermost envelope of a star are severely depleted due to envelope ejection during collisions and due to tidal forces from the SMBH. In the absence of collisions, tidal spin-up of stars is only important in a narrow range of periapse distances, $r_t/2\lesssim r_per \lesssim r_t$ with $r_t$ the tidal disruption radius.

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.

preprint2010arXiv

Orbits Around Black Holes in Triaxial Nuclei

We discuss the properties of orbits within the influence sphere of a supermassive black hole (BH), in the case that the surrounding star cluster is nonaxisymmetric. There are four major orbit families; one of these, the pyramid orbits, have the interesting property that they can approach arbitrarily closely to the BH. We derive the orbit-averaged equations of motion and show that in the limit of weak triaxiality, the pyramid orbits are integrable: the motion consists of a two-dimensional libration of the major axis of the orbit about the short axis of the triaxial figure, with eccentricity varying as a function of the two orientation angles, and reaching unity at the corners. Because pyramid orbits occupy the lowest angular momentum regions of phase space, they compete with collisional loss cone repopulation and with resonant relaxation in supplying matter to BHs. General relativistic advance of the periapse dominates the precession for sufficiently eccentric orbits, and we show that relativity imposes an upper limit to the eccentricity: roughly the value at which the relativistic precession time is equal to the time for torques to change the angular momentum. We argue that this upper limit to the eccentricity should apply also to evolution driven by resonant relaxation, with potentially important consequences for the rate of extreme-mass-ratio inspirals in low-luminosity galaxies. In giant galaxies, we show that capture of stars on pyramid orbits can dominate the feeding of BHs, at least until such a time as the pyramid orbits are depleted; however this time can be of order a Hubble time.

preprint2010arXiv

The Distribution of Stars and Stellar Remnants at the Galactic Center

Motivated by recent observations that suggest a low density of old stars around the Milky Way supermassive black hole, models for the nuclear star cluster are considered that have not yet reached a steady state under the influence of gravitational encounters. A core of initial radius 1-1.5 pc evolves to a size of approximately 0.5 pc after 10 Gyr, roughly the size of the observed core. The absence of a Bahcall-Wolf cusp is naturally explained in these models, without the need for fine-tuning or implausible initial conditions. In the absence of a cusp, the time for a 10-solar-mass black hole to spiral in to the Galactic center from an initial distance of 5 pc can be much greater than 10 Gyr. Assuming that the stellar black holes had the same phase-space distribution initially as the stars, their density after 5-10 Gyr is predicted to rise very steeply going into the stellar core, but could remain substantially below the densities inferred from steady-state models that include a steep density cusp in the stars. Possible mechanisms for the creation of the parsec-scale initial core include destruction of stars on centrophilic orbits in a pre-existing triaxial nucleus, inhibited star formation near the supermassive black hole, or ejection of stars by a massive binary. The implications of these models are discussed for the rates of gravitational-wave inspiral events, as well as other physical processes that depend on a high density of stars or stellar mass black holes near Sagittarius A*.

preprint2010arXiv

The Galactic Centre star S2 as a dynamical probe for intermediate-mass black holes

We study the short-term effects of an intermediate mass black hole (IBH) on the orbit of star S2 (S02), the shortest-period star known to orbit the supermassive black hole (MBH) in the centre of the Milky Way. Near-infrared imaging and spectroscopic observations allow an accurate determination of the orbit of the star. Given S2's short orbital period and large eccentricity, general relativity (GR) needs to be taken into account, and its effects are potentially measurable with current technology. We show that perturbations due to an IBH in orbit around the MBH can produce a shift in the apoapsis of S2 that is as large or even larger than the GR shift. An IBH will also induce changes in the plane of S2's orbit at a level as large as one degree per period. We apply observational orbital fitting techniques to simulations of the S-cluster in the presence of an IBH and find that an IBH more massive than about 1000 solar masses at the distance of the S-stars will be detectable at the next periapse passage of S2, which will occur in 2018.

preprint2010arXiv

The HST/ACS Coma Cluster Survey III. Structural Parameters of Galaxies using single-Sérsic Fits

We present a catalogue of structural parameters for 8814 galaxies in the 25 fields of the HST/ACS Coma Treasury Survey. Parameters from Sérsic fits to the two-dimensional surface brightness distributions are given for all galaxies from our published Coma photometric catalogue with mean effective surface brightness brighter than 26.0 mag/sq. arcsec and brighter than 24.5 mag (equivalent to absolute magnitude - 10.5), as given by the fits, all in F814W(AB). The sample comprises a mixture of Coma members and background objects; 424 galaxies have redshifts and of these 163 are confirmed members. The fits were carried out using both the Gim2D and Galfit codes. We provide the following parameters: Galaxy ID, RA, DEC, the total corrected automatic magnitude from the photometric catalogue, the total magnitude of the model (F814W_AB), the geometric mean effective radius Re, the mean surface brightness within the effective radius <μ>_e, the Sérsic index n, the ellipticity and the source position angle. The selection limits of the catalogue and the errors listed for the Sérsic parameters come from extensive simulations of the fitting process using synthetic galaxy models. The agreement between Gim2D and Galfit parameters is sensitive to details of the fitting procedure; for the settings employed here the agreement is excellent over the range of parameters covered in the catalogue. We define and present two goodness-of-fit indices which quantify the degree to which the image can be approximated by a Sérsic model with concentric, coaxial elliptical isophotes; such indices may be used to objectively select galaxies with more complex structures such as bulge-disk, bars or nuclear components. We make the catalog available in electronic format at Astro-WISE and MAST.

preprint2010arXiv

The HST/ACS Coma Cluster Survey. II. Data Description and Source Catalogs

The Coma cluster was the target of a HST-ACS Treasury program designed for deep imaging in the F475W and F814W passbands. Although our survey was interrupted by the ACS instrument failure in 2007, the partially completed survey still covers ~50% of the core high-density region in Coma. Observations were performed for 25 fields that extend over a wide range of cluster-centric radii (~1.75 Mpc) with a total coverage area of 274 arcmin^2. The majority of the fields are located near the core region of Coma (19/25 pointings) with six additional fields in the south-west region of the cluster. In this paper we present reprocessed images and SExtractor source catalogs for our survey fields, including a detailed description of the methodology used for object detection and photometry, the subtraction of bright galaxies to measure faint underlying objects, and the use of simulations to assess the photometric accuracy and completeness of our catalogs. We also use simulations to perform aperture corrections for the SExtractor Kron magnitudes based only on the measured source flux and half-light radius. We have performed photometry for ~73,000 unique objects; one-half of our detections are brighter than the 10-sigma point-source detection limit at F814W=25.8 mag (AB). The slight majority of objects (60%) are unresolved or only marginally resolved by ACS. We estimate that Coma members are 5-10% of all source detections, which consist of a large population of unresolved objects (primarily GCs but also UCDs) and a wide variety of extended galaxies from a cD galaxy to dwarf LSB galaxies. The red sequence of Coma member galaxies has a constant slope and dispersion across 9 magnitudes (-21<M_F814W<-13). The initial data release for the HST-ACS Coma Treasury program was made available to the public in 2008 August. The images and catalogs described in this study relate to our second data release.

preprint2010arXiv

Towards relativistic orbit fitting of Galactic center stars and pulsars

The S stars orbiting the Galactic center black hole reach speeds of up to a few percent the speed of light during pericenter passage. This makes, for example, S2 at pericenter much more relativistic than known binary pulsars, and opens up new possibilities for testing general relativity. This paper develops a technique for fitting nearly-Keplerian orbits with perturbations from Schwarzschild curvature, frame dragging, and spin-induced torque, to redshift measurements distributed along the orbit but concentrated around pericenter. Both orbital and light-path effects are taken into account. It turns out that absolute calibration of rest-frame frequency is not required. Hence, if pulsars on orbits similar to the S stars are discovered, the technique described here can be applied without change, allowing the much greater accuracies of pulsar timing to be taken advantage of. For example, pulse timing of 3 microsec over one hour amounts to an effective redshift precision of 30 cm/s, enough to measure frame dragging and the quadrupole moment from an S2-like orbit, provided problems like the Newtonian "foreground" due to other masses can be overcome. On the other hand, if stars with orbital periods of order a month are discovered, the same could be accomplished with stellar spectroscopy from the E-ELT at the level of 1 km/s.

preprint2009arXiv

A Counterpart to the Radial Orbit Instability in Triaxial Stellar Systems

Self-consistent solutions for triaxial mass models are highly non-unique. In general, some of these solutions might be dynamically unstable, making them inappropriate as descriptions of steady-state galaxies. Here we demonstrate for the first time the existence in triaxial galaxy models of an instability similar to the radial-orbit instability of spherical models. The instability manifests itself when the number of box orbits, with predominantly radially motions, is sufficiently large. N-body simulations verify that the evolution is due neither to chaotic orbits nor to departures of the model from self-consistency, but rather to a collective mode. The instability transforms the triaxial model into a more prolate, but still triaxial, configuration. Stable triaxial models are obtained when the mass contribution of radial orbits is reduced. The implications of our results for the shapes of dark-matter halos are discussed.

preprint2009arXiv

A VLA radio-survey of early-type galaxies in the Virgo cluster

We present the results of a 8.4 GHz Very Large Array radio survey of early-type galaxies extracted from the ACS Virgo Cluster Survey. The aim of this survey is to investigate the origin of radio emission in early-type galaxies and its link with the host properties in an unexplored territory toward the lowest levels of both radio and optical luminosities. Radio images, available for all 63 galaxies with BT < 14.4, show the presence of a compact radio source in 12 objects, with fluxes spanning from 0.13 to 2700 mJy. The remaining 51 galaxies, undetected at a flux limit of ~0.1 mJy, have radio luminosities L < 4 10E18 W/Hz . The fraction of radio-detected galaxies are a strong function of stellar mass, in agreement with previous results: none of the 30 galaxies with stellar mass M(star) < 1.7 10E10 M(sun) is detected, while 8 of the 11 most massive galaxies have radio cores. There appears to be no simple relation between the presence of a stellar nucleus and radio emission. A multiwavelength analysis of the active galactic nucleus (AGN) emission, combining radio and X-ray data, confirms the link between optical surface brightness profile and radio loudness in the sense that the bright core galaxies are associated with radio-loud AGNs, while non-core galaxies host radio-quiet AGNs. Not all radio-detected galaxies have a X-ray nuclear counter part (and vice-versa). A complete census of AGNs (and supermassive black holes, SMBHs) thus requires observations, at least, in both bands. Nonetheless, there are massive galaxies in the sample, expected to host a large SMBH (M(BH) ~ 10E8 M(sun)), whose nuclear emission eludes detection despite their proximity and the depth and the spatial resolution of the available observations. Most likely this is due to an extremely low level of accretion onto the central SMBH.

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

Intrinsic Shapes of Elliptical Galaxies

Tests for the intrinsic shape of the luminosity distribution in elliptical galaxies are discussed, with an emphasis on the uncertainties. Recent determinations of the ellipticity frequency function imply a paucity of nearly spherical galaxies, and may be inconsistent with the oblate hypothesis. Statistical tests based on the correlation of surface brightness, isophotal twisting, and minor axis rotation with ellipticity have so far not provided strong evidence in favor of the nearly oblate or nearly prolate hypothesis, but are at least qualitatively consistent with triaxiality. The possibility that the observed deviations of elliptical galaxy isophotes form ellipses are due to projection effects is evaluated. Dynamical instabilities may explain the absence of elliptical galaxies flatter than about E6, and my also play a role in the lack of nearly-spherical galaxies.

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.

preprint2009arXiv

Tidal break-up of binary stars at the Galactic center and its consequences

The tidal breakup of binary star systems by the supermassive black hole (SMBH) in the center of the galaxy has been suggested as the source of both the observed sample of hypervelocity stars (HVSs) in the halo of the Galaxy and the S-stars that remain in tight orbits around Sgr A*. Here, we use a post-Newtonian N-body code to study the dynamics of main-sequence binaries on highly elliptical bound orbits whose periapses lie close to the SMBH, determining the properties of ejected and bound stars as well as collision products. Unlike previous studies, we follow binaries that remain bound for several revolutions around the SMBH, finding that in the case of relatively large periapses and highly inclined binaries the Kozai resonance can lead to large periodic oscillations in the internal binary eccentricity and inclination. Collisions and mergers of the binary elements are found to increase significantly for multiple orbits around the SMBH, while HVSs are primarily produced during a binary's first passage. This process can lead to stellar coalescence and eventually serve as an important source of young stars at the galactic center.

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'' 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'', 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 "mass deficit," 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'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'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' 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'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' 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'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". 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 "final parsec problem"). 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'. 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's anisotropic model; Einasto's model (identical in functional form to Sersic'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's R^{1/n} law. Overall, the best-fitting model to the LCDM halos is Einasto'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'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', 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 "turns over" 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

Massive Black Hole Binary Evolution

Coalescence of binary supermassive black holes (SBHs) would constitute the strongest sources of gravitational waves to be observed by LISA. While the formation of binary SBHs during galaxy mergers is almost inevitable, coalescence requires that the separation between binary components first drop by a few orders of magnitude, due presumably to interaction of the binary with stars and gas in a galactic nucleus. This article reviews the observational evidence for binary SBHs and discusses how they would evolve. No completely convincing case of a bound, binary SBH has yet been found, although a handful of systems (e.g. interacting galaxies; remnants of galaxy mergers) are now believed to contain two SBHs at projected separations of roughly one kpc. N-body studies of binary evolution in gas-free galaxies have reached large enough particle numbers to reproduce the slow, "diffusive" refilling of the binary's loss cone that is believed to characterize binary evolution in real galactic nuclei. While some of the results of these simulations -- e.g. the binary hardening rate and eccentricity evolution -- are strongly $N$-dependent, others -- e.g. the "damage" inflicted by the binary on the nucleus -- are not. Luminous early-type galaxies often exhibit depleted cores with masses of 1-2 times the mass of their nuclear SBHs, consistent with the predictions of the binary model. Studies of the interaction of massive binaries with gas are still in their infancy, although much progress is expected in the near future. Binary coalescence has a large influence on the spins of SBHs, even for mass ratios as extreme as 10:1, and evidence of spin-flips may have been observed.

preprint2005arXiv

The ACS Virgo Cluster Survey IX: The Color Distributions of Globular Cluster Systems in Early-Type Galaxies

We present the color distributions of globular cluster (GC) systems for 100 Virgo cluster early-type galaxies observed in the ACS Virgo Cluster Survey. The color distributions of individual GC systems are consistent with continuous trends across galaxy luminosity, color, and stellar mass. On average, almost all galaxies possess a component of metal-poor GCs, with the average fraction of metal-rich GCs ranging from 15 to 60%. The colors of both subpopulations correlate with host galaxy luminosity and color, with the red GCs having a steeper slope. To convert color to metallicity, we also introduce a preliminary (g-z)-[Fe/H] relation calibrated to Galactic, M49 and M87 GCs. This relation is nonlinear with a steeper slope for [Fe/H] < -0.8. As a result, the metallicities of the metal-poor and metal-rich GCs vary similarly with respect to galaxy luminosity and stellar mass, with relations of [Fe/H]_MP ~ L^0.16 ~ M_star^0.17 and [Fe/H]_MR ~ L^0.26 ~ M_star^0.22, respectively. Although these relations are shallower than the mass-metallicity relation predicted by wind models and observed for dwarf galaxies, they are very similar to the mass-metallicity relation for star forming galaxies in the same mass range. The offset between the two GC populations varies slowly (~ M_star^0.05) and is approximately 1 dex across three orders of magnitude in mass, suggesting a nearly universal amount of enrichment between the formation of the two populations of GCs. We also find that although the metal-rich GCs show a larger dispersion in color, it is the *metal-poor GCs* that have an equal or larger dispersion in metallicity. Like the color-magnitude relation, these relations derived from globular clusters present stringent constraints on the formation and evolution of early-type galaxies. (Abridged)

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