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
The question for today
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.
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.
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.
Work it through
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.
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.
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.
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 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.
The living notebook
Explore the model
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.
- y at chosen x
- 0.608
- Tangent slope · 2ax
- 1.35
- Signed area · ax³/3
- 0.182
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.
Pause & reason
Check your understanding
Your turn
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.
Keep the idea
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.