Parameterizing black hole orbits for adiabatic inspiral

Adiabatic binary inspiral in the small mass ratio limit treats the small body as moving along a geodesic of a large Kerr black hole, with the geodesic slowly evolving due to radiative backreaction. Up to initial conditions, geodesics are typically parameterized in two ways: using the integrals of motion energy E, axial angular momentum Lz, and Carter constant Q; or, using orbit geometry parameters semilatus rectum p, eccentricity e, and (cosine of) inclination xI≡cosI. The community has long known how to compute orbit integrals as functions of the orbit geometry parameters, i.e., as functions expressing E(p,e,xI), and likewise for Lz and Q. Mappings in the other direction—functions p(E,Lz,Q), and likewise for e and xI—have not yet been developed in general. In this note, we develop generic mappings from (E,Lz,Q) to (p,e,xI). The mappings are particularly simple for equatorial orbits (Q=0, xI=±1), and can be evaluated efficiently for generic cases. These results make it possible to more accurately compute adiabatic inspirals by eliminating the need to use a Jacobian which becomes singular as inspiral approaches the last stable orbit.