Time-Dependent Models for Dark Matter at the Galactic Center

The prospects of indirect detection of dark matter at the galactic center depend sensitively on the mass profile within the inner parsec. We calculate the distribution of dark matter on sub-parsec scales by integrating the time-dependent Fokker-Planck equation, including the effects of self-annihilations, scattering of dark matter particles by stars, and capture in the supermassive black hole. We consider a variety of initial dark matter distributions, including models with very high densities ("spikes") near the black hole, and models with "adiabatic compression" of the baryons. The annihilation signal after 10 Gyr is found to be substantially reduced from its initial value, but in dark matter models with an initial spike, order-of-magnitude enhancements can persist compared with the rate in spike-free models, with important implications for indirect dark matter searches with GLAST and Air Cherenkov Telescopes like HESS and CANGAROO.Comment: Four page

Instability of the Gravitational N-Body Problem in the Large-N Limit

We use a systolic N-body algorithm to evaluate the linear stability of the gravitational N-body problem for N up to 1.3 x 10^5, two orders of magnitude greater than in previous experiments. For the first time, a clear ~ln N-dependence of the perturbation growth rate is seen. The e-folding time for N = 10^5 is roughly 1/20 of a crossing time.Comment: Accepted for publication in The Astrophysical Journa

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

Nonlinear Fluid Computations in a Distributed Environment

The performance of a loosely and tightly-coupled workstation cluster is compared against a conventional vector supercomputer for the solution the Reynolds- averaged Navier-Stokes equations. The application geometries include a transonic airfoil, a tiltrotor wing/fuselage, and a wing/body/empennage/nacelle transport. Decomposition is of the manager-worker type, with solution of one grid zone per worker process coupled using the PVM message passing library. Task allocation is determined by grid size and processor speed, subject to available memory penalties. Each fluid zone is computed using an implicit diagonal scheme in an overset mesh framework, while relative body motion is accomplished using an additional worker process to re-establish grid communication

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.Comment: 4 pages, 4 figs, revtex

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

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

Long-Term Evolution of Massive Black Hole Binaries. II. Binary Evolution in Low-Density Galaxies

We use direct-summation N-body integrations to follow the evolution of binary black holes at the centers of galaxy models with large, constant-density cores. Particle numbers as large as 400K are considered. The results are compared with the predictions of loss-cone theory, under the assumption that the supply of stars to the binary is limited by the rate at which they can be scattered into the binary's influence sphere by gravitational encounters. The agreement between theory and simulation is quite good; in particular, we are able to quantitatively explain the observed dependence of binary hardening rate on N. We do not verify the recent claim of Chatterjee, Hernquist & Loeb (2003) that the hardening rate of the binary stabilizes when N exceeds a particular value, or that Brownian wandering of the binary has a significant effect on its evolution. When scaled to real galaxies, our results suggest that massive black hole binaries in gas-poor nuclei would be unlikely to reach gravitational-wave coalescence in a Hubble time.Comment: 13 pages, 8 figure