Effective-One-Body Dynamics

Effective-one-body (EOB) theory reorganizes the relativistic two-body problem as the motion of an effective particle in a mass-ratio-deformed black-hole geometry, then couples that motion to resummed radiation and a merger–ringdown completion. Its value is not that two gravitating bodies literally form one exact metric. Rather, EOB supplies a Hamiltonian language in which post-Newtonian (PN), post-Minkowskian (PM), self-force, black-hole-perturbation, and numerical-relativity information constrain complementary functions while exact limiting behavior remains visible.

Helpful background. The Einstein Equation identifies the nonlinear dynamics being approximated; Weak Metric Perturbations and Gauge Freedom supplies weak-field matching; Presymplectic Potentials and Currents supplies canonical transformations and Hamiltonian flow; and Black-Hole Perturbation Theory supplies exact test-body dynamics, fluxes, and ringdown.

Enter the effective one-body description

For masses $m_1$ and $m_2$, define

$$M=m_1+m_2,\qquad \mu=\frac{m_1m_2}{M},\qquad \nu=\frac{\mu}{M},$$

where $0\lt \nu\leq1/4$. The chapter uses $c=1$, retains $G$, and follows the site’s mostly-plus metric convention. Real total energy includes both rest masses; real binding energy has $M$ subtracted; effective energy includes the effective particle’s rest mass $\mu$. Keeping those three quantities distinct prevents many apparent sign and normalization failures.

The central conservative composition is

$$H_{\rm real} =M\sqrt{1+2\nu \left(\frac{H_{\rm eff}}{\mu}-1\right)}.$$

In one standard nonspinning EOB gauge,

$$\left(\frac{H_{\rm eff}}{\mu}\right)^2 =A(u;\nu) \left[ 1+j^2u^2+A\bar D\,p_r^2+Q(u,p_r;\nu) \right],$$

with

$$u=\frac{GM}{r}, \qquad p_r=\frac{P_r}{\mu}, \qquad j=\frac{P_\phi}{GM\mu}.$$

The functions $A$ and $\bar D$ define a deformed effective metric; $Q$ contains nongeodesic radial-momentum terms. At $\nu=0$, they must recover Schwarzschild,

$$A=1-2u,\qquad \bar D=1,\qquad Q=0.$$

At finite $\nu$, these are Hamiltonian potentials rather than metric functions of an exact binary spacetime. Canonical transformations can move terms among them while preserving observables. The original energy-level matching and effective-metric construction are developed in Buonanno and Damour 1999, §§ II–VI.

The six-page route

Read the chapter in order. Each page adds one structural layer and supplies an independent limit.

1. The EOB Energy Map derives the nonlinear real-to-effective energy dictionary from invariant scattering kinematics and PN action matching. Its Newtonian, test-mass, and weak-field circular checks fix the rest-energy conventions.
2. EOB Potentials and Effective Metrics defines $A$, $\bar D$, and $Q$, explains their canonical gauge dependence, and shows how circular, eccentric, and scattering invariants probe different combinations.
3. Spinning EOB Hamiltonians introduces the background and effective-particle spin combinations, gyrogravitomagnetic factors, deformed-Kerr structure, spin gauge, and conservative precession checks.
4. EOB Radiation Reaction and Wave Generation adds a balance-consistent force, factorized waveform modes, test-mass flux constraints, nonquasicircular corrections, and merger–ringdown attachment.
5. Resummation and Strong-Field Completion compares Taylor, Padé, logarithmic, factorized, centrifugal-radius, and PM-informed choices. It is a dated frontier assessment because no finite series uniquely fixes its strong-field continuation.
6. Calibration, Extrapolation, and Uncertainty treats numerical calibration as an inverse problem, documents results for a specified SEOBNRv5 version, separates in-sample and held-out tests, and gives tests that could change the conclusions about waveform systematics.

The first three pages build the conservative Hamiltonian. The fourth changes the problem by imposing retarded boundary conditions and physical loss. The final two distinguish analytically constrained structure from empirically selected completion.

Conservative dynamics and invariant outputs

For nonspinning circular motion, $p_r=0$, so

$$\widehat H_{\rm eff}^2 =A(u)\left(1+j^2u^2\right).$$

Extremizing at fixed $j$ gives

$$j^2(u) =-\frac{A'(u)} {u[2A(u)+uA'(u)]}.$$

In the Schwarzschild limit this becomes

$$j^2=\frac{1}{u(1-3u)}, \qquad \widehat H_{\rm eff} =\frac{1-2u}{\sqrt{1-3u}}.$$

The energy map converts the effective result to the real binding energy, and differentiating the real Hamiltonian yields the orbital frequency. Eliminating the coordinate $u$ produces the invariant relation $E(\Omega)$. This sequence illustrates the right comparison logic: potentials and radii depend on canonical representation; energy versus frequency does not.

Generic conservative motion provides additional invariants. Infinitesimal radial oscillations determine periastron advance and constrain $\bar D$; the quartic-in-$p_r$ potential $Q$ requires finite-eccentricity or scattering information. Hyperbolic motion determines the scattering angle $\chi(E,J)$, which constrains the mass shell at unrestricted velocity. A candidate potential that fits circular energies can still fail either test. Damour’s PM-to-EOB correspondence shows how invariant scattering data fix an effective description without equating gauge-dependent trajectories Damour 2016, §§ II–VI.

Spin adds a deformed Kerr background and an effective particle spin. At leading weak-field order,

$$H_{\rm SO} =\frac{G}{r^3}\mathbf L\cdot \left(2\mathbf S+\frac32\mathbf S_*\right),$$

where

$$\mathbf S=\mathbf S_1+\mathbf S_2, \qquad \mathbf S_*= \frac{m_2}{m_1}\mathbf S_1 +\frac{m_1}{m_2}\mathbf S_2.$$

A valid spinning model must recover zero-spin EOB, Kerr geodesics for a nonspinning test body, the spinning-test-body Hamiltonian at the retained order, PN gyro factors, and conservation of spin magnitudes under conservative flow.

Dissipative closure

Using physical canonical momenta $P_r,P_\phi$ and the binding Hamiltonian $H_{\rm EOB}=H_{\rm real}-M$, radiation reaction changes Hamilton’s momentum equations to

$$\dot P_r=-\frac{\partial H_{\rm EOB}}{\partial r} +\mathcal F_r, \qquad \dot P_\phi=\mathcal F_\phi.$$

For quasicircular motion in a standard reaction gauge,

$$\mathcal F_r\simeq0, \qquad \mathcal F_\phi=-\frac{\Phi_E}{\Omega},$$

where $\Phi_E$ is computed from the same complex waveform modes that define the observable strain. More generally, near-zone Schott energy and angular momentum enter balance, and the simple circular relation must not be used for eccentric or scattering motion without modification.

The inspiral–plunge modes are factorized schematically as

$$\pi_{\ell m}=(\ell+m)\bmod 2,$$

where $\pi_{\ell m}=0$ labels even-parity mass-type modes and $\pi_{\ell m}=1$ labels odd-parity current-type modes. Then

$$H_{\ell m} =H_{\ell m}^{\rm N} \widehat S_{\rm eff}^{(\pi_{\ell m})} T_{\ell m}e^{i\delta_{\ell m}} \left(\rho_{\ell m}\right)^\ell N_{\ell m}.$$

This separates Newtonian amplitude, strong-field source, leading tails, residual phase, residual amplitude, and nonquasicircular corrections. Re-expansion must recover every known PN coefficient; the test-mass limit must agree with black-hole perturbation theory under the same horizon and infinity boundary conditions. The factorized construction and its test-mass comparisons are given in Damour, Iyer, and Nagar 2009, § II and §§ III–IV.

Near merger, the inspiral–plunge mode is attached to a remnant-based ringdown representation. That step uses additional input: matching time, nonquasicircular parameters, remnant mass and spin, quasinormal frequencies, and often numerical calibration. It is not a consequence of the conservative energy map or flux balance alone.

Four information channels

EOB is best understood by asking what each source can actually constrain.

• PN theory gives a weak-field and slow-motion expansion at general mass ratio. It fixes Taylor coefficients, spin couplings, logarithms, and flux terms but not a unique strong-field continuation.
• PM theory and scattering retain unrestricted asymptotic velocity while expanding in $G$. They constrain energy-dependent mass shells and scattering invariants, but a resummation is required before using a finite PM result near capture.
• Self-force and black-hole perturbation theory give strong-field functions and fluxes as expansions around $\nu=0$. Circular data constrain chiefly the $A$ potential; eccentric data are needed to separate radial potentials; finite-$\nu$ completion remains necessary.
• Numerical relativity supplies nonperturbative comparable-mass waveforms and dynamics over a finite, nonuniform parameter set. It can calibrate a small number of flexibility parameters and test completed models, but numerical errors and sparse extrapolation regions remain.

These channels overlap through invariants. Agreement among PN-expanded circular energy, self-force redshift information, PM scattering angle, and numerical binding energy is more diagnostic than one final waveform mismatch because compensating sector errors are harder to hide.

Structural constraints versus empirical choices

Some EOB statements are exact within the declared representation: definitions of $M,\mu,\nu$; the chosen energy map; canonical equations; balance equations including Schott terms; and re-expansion through known analytic orders. Other statements are model choices: Padé degree, logarithmic factorization, centrifugal-radius deformation, effective higher-PN coefficients, nonquasicircular ansatz, and merger attachment.

A strong-field choice earns credibility by passing four layers:

1. algebraic recovery of every known PN, PM, spin, and self-force coefficient;
2. exact Newtonian, Schwarzschild or Kerr test-body, and exchange-symmetry limits;
3. regularity—no unintended poles, complex energies, lost turning points, or discontinuous parameter limits;
4. agreement with held-out invariant and waveform data in a declared domain.

Calibration is then reported separately from validation. A waveform used to fit an EOB coefficient measures in-sample performance, even if its mismatch is extremely small. A withheld configuration, another numerical code, a precession sector not used in fitting, or an independent scattering invariant provides stronger evidence. The published SEOBNRv5HM construction illustrates how high-order analytic, self-force, perturbative, and numerical information are combined, while also documenting the nonuniform calibration domain Pompili et al. 2023, §§ II–V.

Domain and breakdown

The baseline chapter concerns asymptotically flat general relativity and compact binaries described by point-particle Hamiltonians, with black-hole spin-induced multipoles when spin is included. Tides, material modes, horizon absorption, eccentricity, precession, and nonlocal conservative tails require additional terms and, often, additional state variables. A quasicircular aligned-spin validation cannot be silently extended to their simultaneous action.

Analytic truncation becomes stressed at large $u$, large radial momentum, near-capture scattering, and near-extremal spin. Numerical calibration becomes unsupported in sparse high-mass-ratio or high-spin corners and outside the waveform lengths, modes, and orbital families represented by the data. Resummation can extend usefulness but cannot manufacture missing coefficients. The 2025 LISA waveform review emphasizes accuracy, breadth, coupled physical effects, robustness, and computational cost as distinct remaining requirements LISA Consortium Waveform Working Group 2025, §§ 4.5.1–4.5.6.

Review the chapter

Integrated test-mass and Newtonian check

Assume the effective Hamiltonian has the Newtonian expansion

$$\frac{H_{\rm eff}}{\mu} =1+\frac{\mathbf P^2}{2\mu^2} -\frac{GM}{r}+O(v^4).$$

Use the energy map to recover the real Newtonian binding Hamiltonian. Then combine the circular Newtonian energy with the leading quadrupole flux to find $\dot x$, where $x=(GM\Omega)^{2/3}$.

Solution

Set

$$z=\frac{H_{\rm eff}}{\mu}-1=O(v^2).$$

The reduced real binding energy is

$$\frac{H_{\rm real}-M}{\mu} =\frac{\sqrt{1+2\nu z}-1}{\nu} =z-\frac{\nu}{2}z^2+\cdots.$$

The nonlinear term is $O(v^4)$, so at Newtonian order

$$H_{\rm real}-M =\frac{\mathbf P^2}{2\mu} -\frac{GM\mu}{r}+O(v^4).$$

For a circular orbit $E_{\rm N}=-\mu x/2$, while

$$\Phi_E^{\rm N} =\frac{32}{5G}\nu^2x^5.$$

Balance gives

$$\frac{dE_{\rm N}}{dx}\dot x=-\Phi_E^{\rm N},$$

and therefore

$$\dot x =\frac{64}{5GM}\nu x^5\gt 0.$$

The first result checks the rest-energy subtraction and energy map. The second checks the flux normalization and dissipative sign. In the test-mass limit, the binding energy per unit $\mu$ remains finite while the inspiral rate scales as $\nu$, so radiation reaction becomes slow relative to orbital motion.

Diagnose an apparent disagreement

Two EOB models have different functions $A(u)$ but the same circular relation $E(\Omega)$ through 4PN. Does this prove one is wrong?

Solution

No. The coordinate $u$, the canonical gauge, the energy map convention, and momentum-dependent terms must be compared together. A canonical transformation can change $A(u)$ while leaving $E(\Omega)$ unchanged through the retained order. The models become physically distinguishable only through an invariant beyond the shared information—such as higher-order circular energy, periastron advance, scattering angle, or a held-out waveform comparison—with all conventions aligned.

References

• Buonanno, Alessandra, and Thibault Damour. “Effective One-Body Approach to General Relativistic Two-Body Dynamics.” Physical Review D 59, 084006 (1999). DOI. Open PDF.
• Damour, Thibault, Bala R. Iyer, and Alessandro Nagar. “Improved Resummation of Post-Newtonian Multipolar Waveforms from Circularized Compact Binaries.” Physical Review D 79, 064004 (2009). DOI. Open PDF.
• Damour, Thibault. “Gravitational Scattering, Post-Minkowskian Approximation and Effective One-Body Theory.” Physical Review D 94, 104015 (2016). DOI. Open PDF.
• LISA Consortium Waveform Working Group. “Waveform Modelling for the Laser Interferometer Space Antenna.” Living Reviews in Relativity 28, 9 (2025). DOI.
• Pompili, Lorenzo, et al. “Laying the Foundation of the Effective-One-Body Waveform Models SEOBNRv5: Improved Accuracy and Efficiency for Spinning Nonprecessing Binary Black Holes.” Physical Review D 108, 124035 (2023). DOI. Open PDF.