Osculating orbits in Schwarzschild spacetime, with an application to extreme mass-ratio inspirals

We present a method to integrate the equations of motion that govern bound, accelerated orbits in Schwarzschild spacetime. At each instant the true worldline is assumed to lie tangent to a reference geodesic, called an osculating orbit, such that the worldline evolves smoothly from one such geodesic to the next. Because a geodesic is uniquely identified by a set of constant orbital elements, the transition between osculating orbits corresponds to an evolution of the elements. In this paper we derive the evolution equations for a convenient set of orbital elements, assuming that the force acts only within the orbital plane; this is the only restriction that we impose on the formalism, and we do not assume that the force must be small. As an application of our method, we analyze the relative motion of two massive bodies, assuming that one body is much smaller than the other. Using the hybrid Schwarzschild/post-Newtonian equations of motion formulated by Kidder, Will, and Wiseman, we treat the unperturbed motion as geodesic in a Schwarzschild spacetime whose mass parameter is equal to the system's total mass. The force then consists of terms that depend on the system's reduced mass. We highlight the importance of conservative terms in this force, which cause significant long-term changes in the time-dependence and phase of the relative orbit. From our results we infer some general limitations of the radiative approximation to the gravitational self-force, which uses only the dissipative terms in the force.Comment: 18 pages, 6 figures, final version to be published in Physical Review

Intrinsic and extrinsic geometries of a tidally deformed black hole

A description of the event horizon of a perturbed Schwarzschild black hole is provided in terms of the intrinsic and extrinsic geometries of the null hypersurface. This description relies on a Gauss-Codazzi theory of null hypersurfaces embedded in spacetime, which extends the standard theory of spacelike and timelike hypersurfaces involving the first and second fundamental forms. We show that the intrinsic geometry of the event horizon is invariant under a reparameterization of the null generators, and that the extrinsic geometry depends on the parameterization. Stated differently, we show that while the extrinsic geometry depends on the choice of gauge, the intrinsic geometry is gauge invariant. We apply the formalism to solutions to the vacuum field equations that describe a tidally deformed black hole. In a first instance we consider a slowly-varying, quadrupolar tidal field imposed on the black hole, and in a second instance we examine the tide raised during a close parabolic encounter between the black hole and a small orbiting body.Comment: 27 pages, 4 figure

Gravitational waves from binary systems in circular orbits: Convergence of a dressed multipole truncation

The gravitational radiation originating from a compact binary system in circular orbit is usually expressed as an infinite sum over radiative multipole moments. In a slow-motion approximation, each multipole moment is then expressed as a post-Newtonian expansion in powers of v/c, the ratio of the orbital velocity to the speed of light. The bare multipole truncation of the radiation consists in keeping only the leading-order term in the post-Newtonian expansion of each moment, but summing over all the multipole moments. In the case of binary systems with small mass ratios, the bare multipole series was shown in a previous paper to converge for all values v/c < 2/e, where e is the base of natural logarithms. In this paper, we extend the analysis to a dressed multipole truncation of the radiation, in which the leading-order moments are corrected with terms of relative order (v/c)^2 and (v/c)^3. We find that the dressed multipole series converges also for all values v/c < 2/e, and that it coincides (within 1%) with the numerically ``exact'' results for v/c < 0.2.Comment: 9 pages, ReVTeX, 1 postscript figur

Can the post-Newtonian gravitational waveform of an inspiraling binary be improved by solving the energy balance equation numerically?

The detection of gravitational waves from inspiraling compact binaries using matched filtering depends crucially on the availability of accurate template waveforms. We determine whether the accuracy of the templates' phasing can be improved by solving the post-Newtonian energy balance equation numerically, rather than (as is normally done) analytically within the post-Newtonian perturbative expansion. By specializing to the limit of a small mass ratio, we find evidence that there is no gain in accuracy.Comment: 13 pages, RevTeX, 5 figures included via eps

Nonsingular Black Hole Evaporation and ``Stable'' Remnants

We examine the evaporation of two--dimensional black holes, the classical space--times of which are extended geometries, like for example the two--dimensional section of the extremal Reissner--Nordstrom black hole. We find that the evaporation in two particular models proceeds to a stable end--point. This should represent the generic behavior of a certain class of two--dimensional dilaton--gravity models. There are two distinct regimes depending on whether the back--reaction is weak or strong in a certain sense. When the back--reaction is weak, evaporation proceeds via an adiabatic evolution, whereas for strong back--reaction, the decay proceeds in a somewhat surprising manner. Although information loss is inevitable in these models at the semi--classical level, it is rather benign, in that the information is stored in another asymptotic region.Comment: 23 pages, 6 figures, harvmac and epsf, RU-93-12, PUPT-1399, NSF-ITP-93-5

Gravitational radiation from infall into a black hole: Regularization of the Teukolsky equation

The Teukolsky equation has long been known to lead to divergent integrals when it is used to calculate the gravitational radiation emitted when a test mass falls into a black hole from infinity. Two methods have been used in the past to remove those divergent integrals. In the first, integrations by parts are carried out, and the infinite boundary terms are simply discarded. In the second, the Teukolsky equation is transformed into another equation which does not lead to divergent integrals. The purpose of this paper is to show that there is nothing intrinsically wrong with the Teukolsky equation when dealing with non-compact source terms, and that the divergent integrals result simply from an incorrect choice of Green's function. In this paper, regularization of the Teukolsky equation is carried out in an entirely natural way which does not involve modifying the equation.Comment: ReVTeX, 23 page

Killing vectors and anisotropy

We consider an action that can generate fluids with three unequal stresses for metrics with a spacelike Killing vector. The parameters in the action are directly related to the stress anisotropies. The field equations following from the action are applied to an anisotropic cosmological expansion and an extension of the Gott-Hiscock cosmic string

Improved filters for gravitational waves from inspiralling compact binaries

The order of the post-Newtonian expansion needed, to extract in a reliable and accurate manner the fully general relativistic gravitational wave signal from inspiralling compact binaries, is explored. A class of approximate wave forms, called P-approximants, is constructed based on the following two inputs: (a) The introduction of two new energy-type and flux-type functions e(v) and f(v), respectively, (b) the systematic use of Pade approximation for constructing successive approximants of e(v) and f(v). The new P-approximants are not only more effectual (larger overlaps) and more faithful (smaller biases) than the standard Taylor approximants, but also converge faster and monotonically. The presently available O(v/c)^5-accurate post-Newtonian results can be used to construct P-approximate wave forms that provide overlaps with the exact wave form larger than 96.5% implying that more than 90% of potential events can be detected with the aid of P-approximants as opposed to a mere 10-15 % that would be detectable using standard post-Newtonian approximants.Comment: Latex ([prd,aps,eqsecnum,epsf]{revtex}), 40 pages including 12 encapsulated figures. (The paper, together with all the figures and tables is available from ftp://carina.astro.cf.ac.uk/pub/incoming/sathya/dis97.uu

A generalized Damour-Navier-Stokes equation applied to trapping horizons

An identity is derived from Einstein equation for any hypersurface H which can be foliated by spacelike two-dimensional surfaces. In the case where the hypersurface is null, this identity coincides with the two-dimensional Navier-Stokes-like equation obtained by Damour in the membrane approach to a black hole event horizon. In the case where H is spacelike or null and the 2-surfaces are marginally trapped, this identity applies to Hayward's trapping horizons and to the related dynamical horizons recently introduced by Ashtekar and Krishnan. The identity involves a normal fundamental form (normal connection 1-form) of the 2-surface, which can be viewed as a generalization to non-null hypersurfaces of the Hajicek 1-form used by Damour. This 1-form is also used to define the angular momentum of the horizon. The generalized Damour-Navier-Stokes equation leads then to a simple evolution equation for the angular momentum.Comment: Added subsection IV.D; corrected an error in Appendix A; added some references; accepted for publication in Phys. Rev. D (16 pages, 4 EPS figures

Mass Inflation in the Loop Black Hole

In classical general relativity the Cauchy horizon within a two-horizon black hole is unstable via a phenomenon known as mass inflation, in which the mass parameter (and the spacetime curvature) of the black hole diverges at the Cauchy horizon. Here we study this effect for loop black holes -- quantum gravitationally corrected black holes from loop quantum gravity -- whose construction alleviates the $r=0$ singularity present in their classical counterparts. We use a simplified model of mass inflation, which makes use of the generalized DTR relation, to conclude that the Cauchy horizon of loop black holes indeed results in a curvature singularity similar to that found in classical black holes. The DTR relation is of particular utility in the loop black hole because it does not directly rely upon Einstein's field equations. We elucidate some of the interesting and counterintuitive properties of the loop black hole, and corroborate our results using an alternate model of mass inflation due to Ori.Comment: Latex 20 pages, 7 figure