Dynamics and Interactions of Binaries and Neutron Stars in Globular Clusters

We model the dynamics of test binaries in isotropic, multi-mass models of galactic globular clusters. The evolution of binary orbits through the cluster potentials is modeled, including second order diffusion terms, and probabilities for close encounters with field stars are calculated. We carry out Monte Carlo simulations of the effects of the binary--single star encounters on the binary population and distribution in the cluster, and estimate the collision rate for different stellar populations in globular clusters with different structural parameters. Assuming a Salpeter IMF, for low concentration clusters the core encounter rate is dominated by turnoff mass main--sequence stars and medium mass white dwarfs. For high concentration, high density clusters the encounter probabilities are increasingly dominated by neutron stars and heavy white dwarfs. Hence we predict a smaller ratio of blue stragglers and cataclysmic variables to pulsars in high concentration clusters. The total number of millisecond pulsars, and the ratio of single to binary pulsars, is broadly consistent with the observed population, suggesting the binary--single star encounters contribute significantly to the pulsar formation rate in globular clusters, for the whole range of globular cluster types. The number of millisecond pulsars and the ratio of pulsars in different globular clusters is best explained by a total binary fraction comparable to that of the galaxy, and a modest number of primordial neutron stars in the globular clusters.Comment: 59 pages, uuencoded compressed postscript, including 18 figures. Astrophysical Journal Supplements, in pres

Survival of Terrestrial Planets in the Presence of Giant Planet Migration

The presence of ``Hot Jupiters'', Jovian mass planets with very short orbital periods orbiting nearby main sequence stars, has been proposed to be primarily due to the orbital migration of planets formed in orbits initially much further from the parent star. The migration of giant planets would have profound effects on the evolution of inner terrestrial planets in these systems, and previous analyses have assumed that no terrestrial planets survive after migration has occurred. We present numerical simulations showing that a significant fraction of terrestrial planets could survive the migration process, eventually returning to circular orbits relatively close to their original positions. A fraction of the final orbits are in the Habitable Zone, suggesting that planetary systems with close-in giant planets are viable targets for searches for Earth-like habitable planets around other stars.Comment: 5 pages, 3 figures, emulateapj. ApJL in press, referee comments changes and edited for lengt

Transition from adiabatic inspiral to plunge into a spinning black hole

A test particle of mass mu on a bound geodesic of a Kerr black hole of mass M >> mu will slowly inspiral as gravitational radiation extracts energy and angular momentum from its orbit. This inspiral can be considered adiabatic when the orbital period is much shorter than the timescale on which energy is radiated, and quasi-circular when the radial velocity is much less than the azimuthal velocity. Although the inspiral always remains adiabatic provided mu << M, the quasi-circular approximation breaks down as the particle approaches the innermost stable circular orbit (ISCO). In this paper, we relax the quasi-circular approximation and solve the radial equation of motion explicitly near the ISCO. We use the requirement that the test particle's 4-velocity remain properly normalized to calculate a new contribution to the difference between its energy and angular momentum. This difference determines how a black hole's spin changes following a test-particle merger, and can be extrapolated to help predict the mass and spin of the final black hole produced in finite-mass-ratio black-hole mergers. Our new contribution is particularly important for nearly maximally spinning black holes, as it can affect whether a merger produces a naked singularity.Comment: 9 pages, 6 figures, final version published in PRD with minor change

Gravitational waveforms from a point particle orbiting a Schwarzschild black hole

We numerically solve the inhomogeneous Zerilli-Moncrief and Regge-Wheeler equations in the time domain. We obtain the gravitational waveforms produced by a point-particle of mass $\mu$ traveling around a Schwarzschild black hole of mass M on arbitrary bound and unbound orbits. Fluxes of energy and angular momentum at infinity and the event horizon are also calculated. Results for circular orbits, selected cases of eccentric orbits, and parabolic orbits are presented. The numerical results from the time-domain code indicate that, for all three types of orbital motion, black hole absorption contributes less than 1% of the total flux, so long as the orbital radius r_p(t) satisfies r_p(t)> 5M at all times.Comment: revtex4, 24 pages, 23 figures, 3 tables, submitted to PR

Models of Cuspy Triaxial Galaxies

We construct numerical models of mildly triaxial elliptical galaxies with central density cusps. Using a technique we call ``adiabatic squeezing'', we begin with a spherical gamma=1 Hernquist model and apply a drag to the velocities of the particles along each principle axis. The final models are stable in isolation, preserving their density structure and figure shape over many dynamical timescales. The density profile and axial ratios compare well to the observed properties of elliptical galaxies. The orbital structure of these models show a mixture of tubes, boxes, and boxlets, as expected for triaxial systems, with very few chaotic orbits. These N-body realizations of cuspy triaxial galaxies provide a basis for the study of the dynamical evolution of elliptical galaxies.Comment: 14 pages, 7 figures. Accepted by Ap

Maximum black-hole spin from quasi-circular binary mergers

Black holes of mass M must have a spin angular momentum S below the Kerr limit chi = S/M^2 < 1, but whether astrophysical black holes can attain this limiting spin depends on their accretion history. Gas accretion from a thin disk limits the black-hole spin to chi_gas < 0.9980 +- 0.0002, as electromagnetic radiation from this disk with retrograde angular momentum is preferentially absorbed by the black hole. Extrapolation of numerical-relativity simulations of equal-mass binary black-hole mergers to maximum initial spins suggests these mergers yield a maximum spin chi_eq < 0.95. Here we show that for smaller mass ratios q = m/M << 1, the superradiant extraction of angular momentum from the larger black hole imposes a fundamental limit chi_lim < 0.9979 +- 0.0001 on the final black-hole spin even in the test-particle limit q -> 0 of binary black-hole mergers. The nearly equal values of chi_gas and chi_lim imply that measurement of supermassive black-hole spins cannot distinguish a black hole built by gas accretion from one assembled by the gravitational inspiral of a disk of compact stellar remnants. We also show how superradiant scattering alters the mass and spin predicted by models derived from extrapolating test-particle mergers to finite mass ratios.Comment: final version accepted in PRD, new Fig.4 and discussio

Cool Customers in the Stellar Graveyard I: Limits to Extrasolar Planets Around the White Dwarf G29-38

We present high contrast images of the hydrogen white dwarf G 29-38 taken in the near infrared with the Hubble Space Telescope and the Gemini North Telescope as part of a high contrast imaging search for substellar objects in orbit around nearby white dwarfs. We review the current limits on planetary companions for G29-38, the only nearby white dwarf with an infrared excess due to a dust disk. We add our recent observations to these limits to produce extremely tight constraints on the types of possible companions that could be present. No objects $>$ 6 M$_{Jup}$ are detected in our data at projected separations $>$ 12 AU, and no objects $>$ 16 M$_{Jup}$ are detected for separations from 3 to 12 AU, assuming a total system age of 1 Gyr. Limits for companions at separations $<$ 3 AU come from a combination of 2MASS photometry and previous studies of G29-38's pulsations. Our imaging with Gemini cannot confirm a tentative claim for the presence of a low mass brown dwarf. These observations demonstrate that a careful combination of several techniques can probe nearby white dwarfs for large planets and low mass brown dwarfs.Comment: 20 pages, 4 figures, Accepted to Ap

Improved approximate inspirals of test-bodies into Kerr black holes

We present an improved version of the approximate scheme for generating inspirals of test-bodies into a Kerr black hole recently developed by Glampedakis, Hughes and Kennefick. Their original "hybrid" scheme was based on combining exact relativistic expressions for the evolution of the orbital elements (the semi-latus rectum p and eccentricity e) with approximate, weak-field, formula for the energy and angular momentum fluxes, amended by the assumption of constant inclination angle, iota, during the inspiral. Despite the fact that the resulting inspirals were overall well-behaved, certain pathologies remained for orbits in the strong field regime and for orbits which are nearly circular and/or nearly polar. In this paper we eliminate these problems by incorporating an array of improvements in the approximate fluxes. Firstly, we add certain corrections which ensure the correct behaviour of the fluxes in the limit of vanishing eccentricity and/or 90 degrees inclination. Secondly, we use higher order post-Newtonian formulae, adapted for generic orbits. Thirdly, we drop the assumption of constant inclination. Instead, we first evolve the Carter constant by means of an approximate post-Newtonian expression and subsequently extract the evolution of iota. Finally, we improve the evolution of circular orbits by using fits to the angular momentum and inclination evolution determined by Teukolsky based calculations. As an application of the improved scheme we provide a sample of generic Kerr inspirals and for the specific case of nearly circular orbits we locate the critical radius where orbits begin to decircularise under radiation reaction. These easy-to-generate inspirals should become a useful tool for exploring LISA data analysis issues and may ultimately play a role in source detection.Comment: 25 pages, 14 figures, some typos corrected, short section on conservative corrections added, minor changes for consistency with published versio

The Triple Pulsar System PSR B1620-26 in M4

The millisecond pulsar PSR B1620-26, in the globular cluster M4, has a white dwarf companion in a half-year orbit. Anomalously large variations in the pulsar's apparent spin-down rate have suggested the presence of a second companion in a much wider orbit. Using timing observations made on more than seven hundred days spanning eleven years, we confirm this anomalous timing behavior. We explicitly demonstrate, for the first time, that a timing model consisting of the sum of two non-interacting Keplerian orbits can account for the observed signal. Both circular and elliptical orbits are allowed, although highly eccentric orbits require improbable orbital geometries. The motion of the pulsar in the inner orbit is very nearly a Keplerian ellipse, but the tidal effects of the outer companion cause variations in the orbital elements. We have measured the change in the projected semi-major axis of the orbit, which is dominated by precession-driven changes in the orbital inclination. This measurement, along with limits on the rate of change of other orbital elements, can be used to significantly restrict the properties of the outer orbit. We find that the second companion most likely has a mass m~0.01 Msun --- it is almost certainly below the hydrogen burning limit (m<0.036 Msun, 95% confidence) --- and has a current distance from the binary of ~35 AU and orbital period of order one hundred years. Circular (and near-circular) orbits are allowed only if the pulsar magnetic field is ~3x10^9 G, an order of magnitude higher than a typical millisecond pulsar field strength. In this case, the companion has mass m~1.2x10^-3 Msun and orbital period ~62 years.Comment: 12 pages, 6 figures, 3 tables. Very minor clarifications and rewording. Accepted for publication in the Astrophys.