skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Forced motion near black holes

Journal Article · · Physical Review. D, Particles Fields
 [1];  [2];  [3];  [4];  [3]
  1. Institute of Astronomy, Madingley Road, Cambridge CB3 0HA (United Kingdom)
  2. Center for Radiophysics and Space Research, Cornell University, Ithaca, New York 14853 (United States)
  3. Albert-Einstein-Institut, Max-Planck-Institut fuer Gravitationsphysik, D-14476 Golm (Germany)
  4. Theoretical Astrophysics, California Institute of Technology, Pasadena, California 91125 (United States)

We present two methods for integrating forced geodesic equations in the Kerr spacetime. The methods can accommodate arbitrary forces. As a test case, we compute inspirals caused by a simple drag force, mimicking motion in the presence of gas. We verify that both methods give the same results for this simple force. We find that drag generally causes eccentricity to increase throughout the inspiral. This is a relativistic effect qualitatively opposite to what is seen in gravitational-radiation-driven inspirals, and similar to what others have observed in hydrodynamic simulations of gaseous binaries. We provide an analytic explanation by deriving the leading order relativistic correction to the Newtonian dynamics. If observed, an increasing eccentricity would thus provide clear evidence that the inspiral was occurring in a nonvacuum environment. Our two methods are especially useful for evolving orbits in the adiabatic regime. Both use the method of osculating orbits, in which each point on the orbit is characterized by the parameters of the geodesic with the same instantaneous position and velocity. Both methods describe the orbit in terms of the geodesic energy, axial angular momentum, Carter constant, azimuthal phase, and two angular variables that increase monotonically and are relativistic generalizations of the eccentric anomaly. The two methods differ in their treatment of the orbital phases and the representation of the force. In the first method, the geodesic phase and phase constant are evolved together as a single orbital phase parameter, and the force is expressed in terms of its components on the Kinnersley orthonormal tetrad. In the second method, the phase constants of the geodesic motion are evolved separately and the force is expressed in terms of its Boyer-Lindquist components. This second approach is a direct generalization of earlier work by Pound and Poisson [A. Pound and E. Poisson, Phys. Rev. D 77, 044013 (2008).] for planar forces in a Schwarzschild background.

OSTI ID:
21513103
Journal Information:
Physical Review. D, Particles Fields, Vol. 83, Issue 4; Other Information: DOI: 10.1103/PhysRevD.83.044037; (c) 2011 American Institute of Physics; ISSN 0556-2821
Country of Publication:
United States
Language:
English