Evolution of the Dark Matter Distribution at the Galactic Center

Annihilation radiation from neutralino dark matter at the Galactic center (GC) would be greatly enhanced if the dark matter were strongly clustered around the supermassive black hole (SBH). The existence of a dark-matter "spike" is made plausible by the observed, steeply-rising stellar density near the GC SBH. Here the time-dependent equations describing gravitational interaction of the dark matter particles with the stars are solved. Scattering of dark matter particles by stars would substantially lower the dark matter density near the GC SBH over 10^10 yr, due both to kinetic heating, and to capture of dark matter particles by the SBH. This result suggests that enhancements in the dark matter density around a SBH would be modest whether or not the host galaxy had experienced the scouring effects of a binary SBH.Comment: 5 pages, 3 figures. Submitted to Physical Review Letter

Interaction of Recoiling Supermassive Black Holes with Stars in Galactic Nuclei

Supermassive black hole binaries (SMBHBs) are the products of frequent galaxy mergers. The coalescence of the SMBHBs is a distinct source of gravitational wave (GW) radiation. The detections of the strong GW radiation and their possible electromagnetic counterparts are essential. Numerical relativity suggests that the post-merger supermassive black hole (SMBH) gets a kick velocity up to 4000 km/s due to the anisotropic GW radiations. Here we investigate the dynamical co-evolution and interaction of the recoiling SMBHs and their galactic stellar environments with one million direct N-body simulations including the stellar tidal disruption by the recoiling SMBHs. Our results show that the accretion of disrupted stars does not significantly affect the SMBH dynamical evolution. We investigate the stellar tidal disruption rates as a function of the dynamical evolution of oscillating SMBHs in the galactic nuclei. Our simulations show that most of stellar tidal disruptions are contributed by the unbound stars and occur when the oscillating SMBHs pass through the galactic center. The averaged disruption rate is ~10^{-6} M_\odot yr^{-1}, which is about an order of magnitude lower than that by a stationary SMBH at similar galactic nuclei. Our results also show that a bound star cluster is around the oscillating SMBH of about ~ 0.7% the black hole mass. In addition, we discover a massive cloud of unbound stars following the oscillating SMBH. We also investigate the dependence of the results on the SMBH masses and density slopes of the galactic nuclei.Comment: 38 pages, 10 figues; accepted for publication in Ap

The Stellar Kinematic Signature of Massive Black Hole Binaries

The stalling radius of a merging massive binary black hole (BBH) is expected to be below 0".1 even in nearby galaxies (Yu 2002), and thus BBHs are not expected to be spatially resolved in the near future. However, as we show below, a BBH may be detectable through the significantly anisotropic stellar velocity distribution it produces on scales 5-10 times larger than the binary separation. We calculate the velocity distribution of stable orbits near a BBH by solving the restricted three body problem for a BBH embedded in a bulge potential. We present high resolution maps of the projected velocity distribution moments, based on snapshots of ~ 10^8 stable orbits. The kinematic signature of a BBH in the average velocity maps is a counter rotating torus of stars outside the BBH Hill spheres. The velocity dispersion maps reveal a dip in the inner region, and an excess of 20-40% further out, compared to a single BH of the same total mass. More pronounced signatures are seen in the third and fourth Gauss-Hermite velocity moments maps. The detection of these signatures may indicate the presence of a BBH currently, or at some earlier time, which depends on the rate of velocity phase space mixing following the BBH merger.Comment: Accepted to MNRA

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*.Comment: ApJ, accepte

Secular Stellar Dynamics near a Massive Black Hole

The angular momentum evolution of stars close to massive black holes (MBHs) is driven by secular torques. In contrast to two-body relaxation, where interactions between stars are incoherent, the resulting resonant relaxation (RR) process is characterized by coherence times of hundreds of orbital periods. In this paper, we show that all the statistical properties of RR can be reproduced in an autoregressive moving average (ARMA) model. We use the ARMA model, calibrated with extensive N-body simulations, to analyze the long-term evolution of stellar systems around MBHs with Monte Carlo simulations. We show that for a single-mass system in steady-state, a depression is carved out near an MBH as a result of tidal disruptions. Using Galactic center parameters, the extent of the depression is about 0.1 pc, of similar order to but less than the size of the observed "hole" in the distribution of bright late-type stars. We also find that the velocity vectors of stars around an MBH are locally not isotropic. In a second application, we evolve the highly eccentric orbits that result from the tidal disruption of binary stars, which are considered to be plausible precursors of the "S-stars" in the Galactic center. We find that RR predicts more highly eccentric (e > 0.9) S-star orbits than have been observed to date.Comment: 24 pages, 31 figures; final version as published in Ap

The business of farming : a youth publication

Cover title

Self consistent model for the evolution of eccentric massive black hole binaries in stellar environments: implications for gravitational wave observations

We construct evolutionary tracks for massive black hole binaries (MBHBs) embedded in a surrounding distribution of stars. The dynamics of the binary is evolved by taking into account the erosion of the central stellar cusp bound to the massive black holes, the scattering of unbound stars feeding the binary loss cone, and the emission of gravitational waves (GWs). Stellar dynamics is treated in a hybrid fashion by coupling the results of numerical 3-body scattering experiments of bound and unbound stars to an analytical framework for the evolution of the stellar density distribution and for the efficiency of the binary loss cone refilling. Our main focus is on the behaviour of the binary eccentricity, in the attempt of addressing its importance in the merger process and its possible impact for GW detection with the planned Laser Interferometer Space Antenna ({\it LISA}), and ongoing and forthcoming pulsar timing array (PTA) campaigns. We produce a family of evolutionary tracks extensively sampling the relevant parameters of the system which are the binary mass, mass ratio and initial eccentricity, the slope of the stellar density distribution, its normalization and the efficiency of loss cone refilling. We find that, in general, stellar dynamics causes a dramatic increase of the MBHB eccentricity, especially for initially already mildly eccentric and/or unequal mass binaries. When applied to standard MBHB population models, our results predict eccentricities in the ranges $10^{-3}-0.2$ and $0.03-0.3$ for sources detectable by {\it LISA} and PTA respectively. Such figures may have a significant impact on the signal modelling, on source detection, and on the development of parameter estimation algorithms.Comment: 15 pages, 9 figures, accepted for publication in the Astrophysical Journa

The Tidal Disruption of Giant Stars and Their Contribution to the Flaring Supermassive Black Hole Population

Sun-like stars are thought to be regularly disrupted by supermassive black holes (SMBHs) within galactic nuclei. Yet, as stars evolve off the main sequence their vulnerability to tidal disruption increases drastically as they develop a bifurcated structure consisting of a dense core and a tenuous envelope. Here we present the first hydrodynamic simulations of the tidal disruption of giant stars and show that the core has a substantial influence on the star's ability to survive the encounter. Stars with more massive cores retain large fractions of their envelope mass, even in deep encounters. Accretion flares resulting from the disruption of giant stars should last for tens to hundreds of years. Their characteristic signature in transient searches would not be the $t^{-5/3}$ decay typically associated with tidal disruption events, but a correlated rise over many orders of magnitude in brightness on months to years timescales. We calculate the relative disruption rates of stars of varying evolutionary stages in typical galactic centers, then use our results to produce Monte Carlo realizations of the expected flaring event populations. We find that the demographics of tidal disruption flares are strongly dependent on both stellar and black hole mass, especially near the limiting SMBH mass scale of $\sim 10^8 M_\odot$. At this black hole mass, we predict a sharp transition in the SMBH flaring diet beyond which all observable disruptions arise from evolved stars, accompanied by a dramatic cutoff in the overall tidal disruption flaring rate. Black holes less massive than this limiting mass scale will show observable flares from both main sequence and evolved stars, with giants contributing up to 10% of the event rate. The relative fractions of stars disrupted at different evolutionary states can constrain the properties and distributions of stars in galactic nuclei other than our own.Comment: 18 pages, 18 figures, submitted to Ap

Long Term Evolution of Massive Black Hole Binaries

The long-term evolution of massive black hole binaries at the centers of galaxies is studied in a variety of physical regimes, with the aim of resolving the ``final parsec problem,'' i.e., how black hole binaries manage to shrink to separations at which emission of gravity waves becomes efficient. A binary ejects stars by the gravitational slingshot and carves out a loss cone in the host galaxy. Continued decay of the binary requires a refilling of the loss cone. We show that the standard treatment of loss cone refilling, derived for collisionally relaxed systems like globular clusters, can substantially underestimate the refilling rates in galactic nuclei. We derive expressions for non-equilibrium loss-cone dynamics and calculate time scales for the decay of massive black hole binaries following galaxy mergers, obtaining significantly higher decay rates than heretofore. Even in the absence of two-body relaxation, decay of binaries can persist due to repeated ejection of stars returning to the nucleus on eccentric orbits. We show that this recycling of stars leads to a gradual, approximately logarithmic dependence of the binary binding energy on time. We derive an expression for the loss cone refilling induced by the Brownian motion of a black hole binary. We also show that numerical N-body experiments are not well suited to probe these mechanisms over long times due to spurious relaxation.Comment: Replaced to match the accepted version, ApJ, 596 (2003