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

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.

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

Blue straggler stars in dwarf spheroidal galaxies

Blue straggler star (BSS) candidates have been observed in all old dwarf spheroidal galaxies (dSphs), however whether or not they are authentic BSSs or young stars has been a point of debate. To both address this issue and obtain a better understanding of the formation of BSSs in different environments, we have analysed a sample of BSS candidates in two nearby Galactic dSphs, Draco and Ursa Minor. We have determined their radial and luminosity distributions from wide field multicolour imaging data extending beyond the tidal radii of both galaxies. BSS candidates are uniformly distributed through the host galaxy, whereas a young population is expected to show a more clumpy distribution. Furthermore, the observed radial distribution of BSSs, normalized to both red giant branch (RGB) and horizontal branch (HB) stars, is almost flat, with a slight decrease towards the centre. Such a distribution is at odds with the predictions for a young stellar population, which should be more concentrated. Instead, it is consistent with model predictions for BSS formation by mass transfer in binaries (MT-BSSs). Such results, although not decisive, suggest that these candidates are indeed BSSs and that MT-BSSs form in the same way in Draco and Ursa Minor as in globular clusters. This favours the conclusion that Draco and Ursa Minor are truly ‘fossil' galaxies, where star formation ceased completely more than 8 billion years ag

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

Blue straggler stars in dwarf spheroidal galaxies

Blue straggler star (BSS) candidates have been observed in all old dwarf spheroidal galaxies (dSphs), however whether or not they are authentic BSSs or young stars has been a point of debate. To both address this issue and obtain a better understanding of the formation of BSSs in different environments we have analysed a sample of BSS candidates in two nearby Galactic dSphs, Draco and Ursa Minor. We have determined their radial and luminosity distributions from wide field multicolour imaging data extending beyond the tidal radii of both galaxies. BSS candidates are uniformly distributed through the host galaxy, whereas a young population is expected to show a more clumpy distribution. Furthermore, the observed radial distribution of BSSs, normalized to both red giant branch (RGB) and horizontal branch (HB) stars, is almost flat, with a slight decrease towards the centre. Such a distribution is at odds with the predictions for a young stellar population, which should be more concentrated. Instead, it is consistent with model predictions for BSS formation by mass transfer in binaries (MT-BSSs). Such results, although not decisive, suggest that these candidates are indeed BSSs and that MT-BSSs form in the same way in Draco and Ursa Minor as in globular clusters. This favours the conclusion that Draco and Ursa Minor are truly 'fossil' galaxies, where star formation ceased completely more than 8 billion years ago.Comment: 16 pages, 10 figures, MNRAS, in pres

Blue straggler stars in dwarf spheroidal galaxies - II. Sculptor and Fornax

The existence of blue straggler stars (BSSs) in dwarf spheroidal galaxies (dSphs) is still an open question. In fact, many BSS candidates have been observed in the Local Group dSphs, but it is unclear whether they are real BSSs or young stars. Shedding light on the nature of these BSS candidates is crucial in order to understand the star formation history of dSphs. In this paper, we consider BSS candidates in Sculptor and Fornax. In Fornax, there are strong hints that the BSS population is contaminated by young stars, whereas in Sculptor there is no clear evidence of recent star formation. We derive the radial and luminosity distribution of BSS candidates from wide field imaging data extending beyond the nominal tidal radius of these galaxies. The observations are compared with the radial distribution of BSSs expected from dynamical simulations. In Sculptor, the radial distribution of BSS candidates is consistent with that of red horizontal branch (RHB) stars and is in agreement with theoretical expectations for BSSs generated via mass transfer in binaries. On the contrary, in Fornax, the radial distribution of BSS candidates is more concentrated than that of all the considered stellar populations. This result supports the hypothesis that most of BSS candidates in Fornax are young stars, and this is consistent with previous studie