Verifying black hole orbits with gravitational spectroscopy

Gravitational waves from test masses bound to geodesic orbits of rotating black holes are simulated, using Teukolsky's black hole perturbation formalism, for about ten thousand generic orbital configurations. Each binary radiates power exclusively in modes with frequencies that are integer-linear combinations of the orbit's three fundamental frequencies. General spectral properties are found with a survey of orbits about a black hole taken to be rotating at 80% of the maximal spin. The orbital eccentricity is varied from 0.1 to 0.9. Inclination ranges from 20\ifmmode^\circ\else\textdegree\fi{} to 160\ifmmode^\circ\else\textdegree\fi{} and comes to within 20\ifmmode^\circ\else\textdegree\fi{} of polar. Semilatus rectum is varied from 1.2 to 3 times the value at the innermost stable circular orbits. The following general spectral properties are found: (i) 99% of the radiated power is typically carried by a few hundred modes, and at most by about a thousand modes, (ii) the dominant frequencies can be grouped into a small number of families defined by fixing two of the three integer frequency multipliers, and (iii) the specifics of these trends can be qualitatively inferred from the geometry of the orbit under consideration. Detections using triperiodic analytic templates modeled on these general properties would constitute a verification of radiation from an adiabatic sequence of black hole orbits and would recover the evolution of the fundamental orbital frequencies. In an analogy with ordinary spectroscopy, this would compare to observing the Bohr model's atomic hydrogen spectrum without being able to rule out alternative atomic theories or nuclei. The suitability of such a detection technique is demonstrated using snapshots computed at 12-hour intervals throughout the last three years before merger of a kludged inspiral. The system chosen is typical of those thought to occur in galactic nuclei and to be observable with space-based gravitational wave detectors like LISA. Because of circularization, the number of excited modes decreases as the binary evolves. A hypothetical detection algorithm that tracks mode families dominating the first 12 hours of the inspiral would capture 98% of the total power over the remaining three years.