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.Comment: Accepted for publication in The Astrophysical Journa

Long-Term Evolution of Massive Black Hole Binaries. III. Binary Evolution in Collisional Nuclei

[Abridged] In galactic nuclei with sufficiently short relaxation times, binary supermassive black holes can evolve beyond their stalling radii via continued interaction with stars. We study this "collisional" evolutionary regime using both fully self-consistent N-body integrations and approximate Fokker-Planck models. The N-body integrations employ particle numbers up to 0.26M and a direct-summation potential solver; close interactions involving the binary are treated using a new implementation of the Mikkola-Aarseth chain regularization algorithm. Even at these large values of N, two-body scattering occurs at high enough rates in the simulations that they can not be simply scaled to the large-N regime of real galaxies. The Fokker-Planck model is used to bridge this gap; it includes, for the first time, binary-induced changes in the stellar density and potential. The Fokker-Planck model is shown to accurately reproduce the results of the N-body integrations, and is then extended to the much larger N regime of real galaxies. Analytic expressions are derived that accurately reproduce the time dependence of the binary semi-major axis as predicted by the Fokker-Planck model. Gravitational wave coalescence is shown to occur in <10 Gyr in nuclei with velocity dispersions below about 80 km/s. Formation of a core results from a competition between ejection of stars by the binary and re-supply of depleted orbits via two-body scattering. Mass deficits as large as ~4 times the binary mass are produced before coalescence. After the two black holes coalesce, a Bahcall-Wolf cusp appears around the single hole in one relaxation time, resulting in a nuclear density profile consisting of a flat core with an inner, compact cluster, similar to what is observed at the centers of low-luminosity spheroids.Comment: 21 page

Self-consistent models of cuspy triaxial galaxies with dark matter haloes

We have constructed realistic, self-consistent models of triaxial elliptical galaxies embedded in triaxial dark matter haloes. We examined three different models for the shape of the dark matter halo: (i) the same axis ratios as the luminous matter (0.7:0.86:1); (ii) a more prolate shape (0.5:0.66:1); (iii) a more oblate shape (0.7:0.93:1). The models were obtained by means of the standard orbital superposition technique introduced by Schwarzschild. Self-consistent solutions were found in each of the three cases. Chaotic orbits were found to be important in all of the models,and their presence was shown to imply a possible slow evolution of the shapes of the haloes. Our results demonstrate for the first time that triaxial dark matter haloes can co-exist with triaxial galaxies.Comment: Latex paper based on the AASTEX format, 20 pages, 11 figures, 2 tables. Paper submitted to Ap

Australian Glow-worms in Caves

Glow-worms are the larvae of a fly from the family Keroplatidae. Their closest relatives are the “fungus flies” that seek out mushrooms for their larvae to consume. Glow-worms have gone out on an evolutionary limb, albeit a successful one. They have lost their association with fungi and have instead become carnivorous. The unique feature of glow-worms is their ability to bioluminesce—to produce light

Spin Flips and Precession in Black-Hole-Binary Mergers

We use the `moving puncture' approach to perform fully non-linear evolutions of spinning quasi-circular black-hole binaries with individual spins not aligned with the orbital angular momentum. We evolve configurations with the individual spins (parallel and equal in magnitude) pointing in the orbital plane and 45-degrees above the orbital plane. We introduce a technique to measure the spin direction and track the precession of the spin during the merger, as well as measure the spin flip in the remnant horizon. The former configuration completes 1.75 orbits before merging, with the spin precessing by 98-degrees and the final remnant horizon spin flipped by ~72-degrees with respect to the component spins. The latter configuration completes 2.25 orbits, with the spins precessing by 151-degrees and the final remnant horizon spin flipped ~34-degrees with respect to the component spins. These simulations show for the first time how the spins are reoriented during the final stage of binary black hole mergers verifying the hypothesis of the spin-flip phenomenon. We also compute the track of the holes before merger and observe a precession of the orbital plane with frequency similar to the orbital frequency and amplitude increasing with time.Comment: Revtex4, 17 figures, 14 pages. Accepted for publication in PR

Is there a Supermassive Black Hole at the Center of the Milky Way?

This review outlines the observations that now provide an overwhelming scientific case that the center of our Milky Way Galaxy harbors a supermassive black hole. Observations at infrared wavelength trace stars that orbit about a common focal position and require a central mass (M) of 4 million solar masses within a radius of 100 Astronomical Units. Orbital speeds have been observed to exceed 5,000 km/s. At the focal position there is an extremely compact radio source (Sgr A*), whose apparent size is near the Schwarzschild radius (2GM/c^2). This radio source is motionless at the ~1 km/s level at the dynamical center of the Galaxy. The mass density required by these observations is now approaching the ultimate limit of a supermassive black hole within the last stable orbit for matter near the event horizon.Comment: Invited review submitted to International Journal of Modern Physics D; 23 pages; 10 figure

"Kludge" gravitational waveforms for a test-body orbiting a Kerr black hole

One of the most exciting potential sources of gravitational waves for low-frequency, space-based gravitational wave (GW) detectors such as the proposed Laser Interferometer Space Antenna (LISA) is the inspiral of compact objects into massive black holes in the centers of galaxies. The detection of waves from such "extreme mass ratio inspiral" systems (EMRIs) and extraction of information from those waves require template waveforms. The systems' extreme mass ratio means that their waveforms can be determined accurately using black hole perturbation theory. Such calculations are computationally very expensive. There is a pressing need for families of approximate waveforms that may be generated cheaply and quickly but which still capture the main features of true waveforms. In this paper, we introduce a family of such "kludge" waveforms and describe ways to generate them. We assess performance of the introduced approximations by comparing "kludge" waveforms to accurate waveforms obtained by solving the Teukolsky equation in the adiabatic limit (neglecting GW backreaction). We find that the kludge waveforms do extremely well at approximating the true gravitational waveform, having overlaps with the Teukolsky waveforms of 95% or higher over most of the parameter space for which comparisons can currently be made. Indeed, we find these kludges to be of such high quality (despite their ease of calculation) that it is possible they may play some role in the final search of LISA data for EMRIs.Comment: 29 pages, 11 figures, requires subeqnarray; v2 contains minor changes for consistency with published versio

Chaotic mixing in noisy Hamiltonian systems

This paper summarises an investigation of the effects of low amplitude noise and periodic driving on phase space transport in 3-D Hamiltonian systems, a problem directly applicable to systems like galaxies, where such perturbations reflect internal irregularities and.or a surrounding environment. A new diagnsotic tool is exploited to quantify how, over long times, different segments of the same chaotic orbit can exhibit very different amounts of chaos. First passage time experiments are used to study how small perturbations of an individual orbit can dramatically accelerate phase space transport, allowing `sticky' chaotic orbits trapped near regular islands to become unstuck on suprisingly short time scales. Small perturbations are also studied in the context of orbit ensembles with the aim of understanding how such irregularities can increase the efficacy of chaotic mixing. For both noise and periodic driving, the effect of the perturbation scales roughly in amplitude. For white noise, the details are unimportant: additive and multiplicative noise tend to have similar effects and the presence or absence of a friction related to the noise by a Fluctuation- Dissipation Theorem is largely irrelevant. Allowing for coloured noise can significantly decrease the efficacy of the perturbation, but only when the autocorrelation time, which vanishes for white noise, becomes so large that t here is little power at frequencies comparable to the natural frequencies of the unperturbed orbit. This suggests strongly that noise-induced extrinsic diffusion, like modulational diffusion associated with periodic driving, is a resonance phenomenon. Potential implications for galaxies are discussed.Comment: 15 pages including 18 figures, uses MNRAS LaTeX macro