research33 min studyLesson 836 in Physics

Post-Newtonian expansion and radiation reaction

Organize circular-binary corrections with an invariant frequency parameter and distinguish conservative dynamics, averaged flux and dissipative reaction.

On this page

What you’ll understand

  • Identify the velocity ordering behind a post-Newtonian expansion.
  • Evaluate a first-order conservative energy correction.
  • Check adiabaticity and explain why truncation and instantaneous balance require care.

Post-Newtonian theory expands general relativity in slow motion and weak orbital fields, while allowing compact constituents to be strongly self-gravitating through appropriate matching. It supplies corrections to conservative binding energy, radiated flux and accumulated waveform phase. These are related calculations, but they should not be merged into one undifferentiated correction factor. A small velocity parameter organizes terms; its size alone does not provide a rigorous error bound or justify extrapolation through merger.

A frequency parameter labels circular dynamics

For a nonspinning nearly circular binary define x from total mass and orbital angular frequency. At Newtonian order it equals orbital speed squared over c squared and GM over separation c squared. Defining it through frequency is preferable to treating a coordinate-dependent separation as an invariant observable. One post-Newtonian order represents an additional factor of x relative to the leading conservative result. Half-integer orders can occur in radiation-related quantities.

The conservative binding energy has a known expansion in x and symmetric mass ratio. Its first correction modifies the relation between orbital energy and frequency, rather than adding a measured wave amplitude directly. Coordinate gauges and field redefinitions can change intermediate accelerations or potentials while preserving appropriately defined circular energy-frequency relations. Comparisons between calculations must therefore use matched invariant quantities and mass conventions.

x=(GMΩc3)2/3,E(x)=−12μc2x[1+(−34−ν12)x+O(x2)]x=\left(\frac{GM\Omega}{c^3}\right)^{2/3},\qquad E(x)=-\frac12\mu c^2x\left[1+\left(-\frac34-\frac{\nu}{12}\right)x+O(x^2)\right]
x and ν=μ/M are dimensionless, E is binding energy excluding rest mass, and Ω is orbital angular frequency in rad/s. The displayed 1PN conservative expression assumes nonspinning circular compact-binary dynamics in the PN regime; the omitted terms are not assigned an automatic numerical uncertainty.↔ Scroll wide equations horizontally.

Dissipation appears on a longer timescale

The leading gravitational-wave luminosity scales as the fifth power of x. Combining it with the derivative of binding energy gives slow orbital-frequency evolution. The leading dissipative acceleration is called a relative two-and-a-half-post-Newtonian effect compared with Newtonian acceleration. That label describes its velocity ordering, not a statement that waveform flux begins only at a high relative correction to its own leading radiation formula.

For adiabatic circular inspiral, orbital evolution is slow compared with one orbital period. A convenient local diagnostic is the magnitude of frequency derivative divided by frequency squared. Its leading expression grows strongly with x, so increasing orbital speed erodes the timescale hierarchy. Near plunge, orbit averaging and an expansion around a sequence of stable circular orbits become inadequate even if a truncated expression remains easy to evaluate numerically.

FN=325c5Gν2x5,dEdt≃−F,ϵad=Ω˙Ω2=965νx5/2at leading order\mathcal F_N=\frac{32}{5}\frac{c^5}{G}\nu^2x^5,\qquad \frac{dE}{dt}\simeq-\mathcal F,\qquad \epsilon_{\mathrm{ad}}=\frac{\dot\Omega}{\Omega^2}=\frac{96}{5}\nu x^{5/2}\quad\text{at leading order}
Flux F has power units and E energy units. The balance equation is the stated secular circular approximation. εad is dimensionless and should be small for slow orbital evolution; its smallness is distinct from accuracy of every PN coefficient or the absence of other loss channels.↔ Scroll wide equations horizontally.

Backreaction and truncation have controlled meanings

Instantaneous source-mechanical energy loss need not equal the far-zone flux at every time. Near-zone field energy can exchange reversibly with matter, producing additional total-derivative terms often called Schott terms. Orbit-averaged balance eliminates or accounts for such terms in the appropriate setting. Using a secular flux to invent an arbitrary local damping force would fail to specify the gauge, reaction dynamics and angular-momentum balance needed for consistency.

PN corrections include hereditary effects depending on past evolution, such as tails from wave propagation on the source's background curvature. Coefficients and logarithms can make successive terms larger or smaller than a naive power estimate suggests. Waveform phasing accumulates small local errors across many cycles, so convergence should be assessed with invariant diagnostics, comparisons between approximations and other controlled calculations. Numerical relativity and resummation complement the expansion rather than turning every truncated formula into an exact strong-field solution.

A first post-Newtonian binding correction

For a nonspinning equal-mass circular binary take ν=1/4 and x=0.040. Use the displayed 1PN energy to calculate E/(μc²), its Newtonian value, and the fractional first correction to the bracket.

  1. The first coefficient is −3/4−(1/4)/12=−37/48. Multiplying by x=0.040 gives −0.0308333, so the bracket is approximately 0.9691667.

  2. The Newtonian energy in these units is −x/2=−0.0200000. Multiplying by the corrected bracket gives E/(μc²)≈−0.0193833. The binding energy remains negative.

  3. The first bracket correction is −3.0833 percent. Its negative sign makes the displayed binding-energy magnitude smaller at fixed x. It is not a 3.0833 percent certified error bar on the true energy or on an accumulated waveform phase.

The result

The 1PN expression gives E/(μc²)≈−0.0193833 versus −0.0200000 at Newtonian order, with a −3.0833 percent bracket correction at the chosen frequency parameter.

Explore the model

Interactive lab

Interactive lab · function

Compare one omitted-order power

Use the parabola only for a unit-coefficient x² term. Reset and set both controls to Coefficient a=1 and Position x=0.05; predict y=0.0025. For a fresh comparison, Reset and set both controls to a=1 and x=0.10; predict y=0.01. Doubling x quadruples this nominal next-order power.

Compare one omitted-order powerThis lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Current readings: y at chosen x: 0.608; Tangent slope · 2ax: 1.35; Signed area · ax³/3: 0.182. All parameters are adjustable with the labeled controls below.−2.4−1.201.22.4−9−4.504.59xy(x, ax²)
y = ax²TangentSigned integral
y at chosen x
0.608
Tangent slope · 2ax
1.35
Signed area · ax³/3
0.182
t = 0 s

This lab uses only y = ax². Playback traces the curve within the fixed integration interval from 0 to x; it does not estimate numerical error or change the chosen point. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forPower counting organizes an expansion but does not determine coefficients or truncation error. The quadratic lab contains no binary, PN Hamiltonian, tail effect or radiation-reaction force.

Check your understanding

What does a small leading adiabaticity parameter establish?

Try it yourself

Use εad=(96/5)νx^(5/2) with ν=1/4. Evaluate it at x=0.010 and x=0.040. By what factor does the diagnostic increase, and do the results alone prove a waveform's accuracy?

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

The prefactor is 4.8. Increasing x by four multiplies x^(5/2) by 32.

Compare with the solution

At 0.010 the diagnostic is 4.8×10⁻⁵. At 0.040 it is 4.8×0.00032=0.001536, an increase by 32. Both are small compared with one, indicating slow leading orbital evolution at those points. They do not prove waveform accuracy: omitted PN terms, long accumulated phase, spins, eccentricity and model systematics require separate assessment. A timescale condition and a quantitative error estimate answer different questions.

Take it with you

  • The invariant circular frequency parameter x organizes PN corrections without treating coordinate radius as an observable.
  • Conservative energy, radiative flux and reaction dynamics have different order labels and roles.
  • Adiabaticity, averaged balance and truncation accuracy are separate requirements.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.