Quadrupole radiation and energy flux
Connect a slowly moving source quadrupole to far-zone TT strain and correctly normalized cycle-averaged energy transport.
On this page
The question for today
What you’ll understand
- State the assumptions and projection in the leading quadrupole waveform.
- Calculate a monochromatic flux using the polarization normalization.
- Distinguish averaged wave energy, source backreaction and memory information.
The leading gravitational radiation from an isolated slowly moving source is controlled by the changing mass quadrupole. A conserved monopole and center-of-mass dipole cannot supply the same radiative mechanism. The resulting far-zone strain and energy flux are related but different quantities: strain falls inversely with distance, while energy flux is quadratic in time derivatives of strain. Projection, index sums and averaging conventions determine the numerical factors, so a compact formula must be accompanied by its approximation regime.
Retarded source structure produces a transverse waveform
For a leading slow-motion treatment, take an isolated system with weak orbital gravity, nonrelativistic motion and size small compared with the emitted wavelength. Construct its trace-free mass quadrupole by integrating density times the coordinate product with its spatial trace removed. At an observer far beyond both source size and wavelength, the leading radiative field depends on the quadrupole's second time derivative at retarded time. The observer's propagation direction selects its TT projection.
The projection is essential: an arbitrary quadrupole component is not independently visible along every line of sight. For self-gravitating motion, the conservation argument must also include gravitational field stress; conserving matter stress alone is insufficient. Conservation laws remove independent mass-monopole and center-of-mass dipole radiation for the isolated system. Compact binaries with strongly gravitating constituents can still use a leading orbital quadrupole description through a matched post-Newtonian treatment of their exterior dynamics. That extension is more subtle than interpreting a Newtonian mass-density integral literally inside a black-hole horizon.
Energy transport is quadratic and averaged
At first order a wave solves the linear field equation, but its carried energy appears at second order in amplitude. In a short-wavelength expansion on a slowly varying background, averaging the quadratic derivatives defines an effective wave stress tensor. Average over many oscillations while remaining short compared with source evolution and background-curvature timescales. This procedure supplies a useful leading gauge-independent energy description, rather than a unique local gravitational energy density at every spacetime point.
For a locally outward plane wave, the tensor derivative contraction equals twice the sum of the plus and cross derivative squares with our normalization. The luminosity obtained by integrating over directions can also be written in terms of the third quadrupole derivative. Its repeated spatial indices include both off-diagonal entries. The outgoing luminosity is positive, while the corresponding orbit-averaged source energy change has the opposite sign when radiation is the relevant loss channel.
Memory and propagation contain additional information
Averaged energy balance produces slow backreaction on a source and underlies the inspiral chirp. It should not be mistaken for a force obtained by inserting a local wave-energy density into a Newtonian equation without a consistent approximation. Bound orbital energy, near-zone field energy and radiated flux have distinct roles; instantaneous balance can include additional reversible field-energy terms even when the secular loss is correctly predicted.
Gravitational-wave memory is a lasting change in an appropriate relative-displacement observable after a burst. Ordinary memory can arise from changes in unbound matter motion; nonlinear memory is sourced by radiated energy itself. A single monochromatic oscillatory quadrupole example does not determine that offset. Its calculation needs the relevant source history and asymptotic or propagation assumptions. Cosmological redshift and distance likewise belong to propagation modeling before an observer-frame strain is contracted with a detector response.
Work it through
Flux from one sinusoidal polarization
A local outgoing plane wave has h+=A cos(2πft), h×=0, A=10⁻²¹ and f=100 Hz. Use c³/(32πG)=4.02×10³³ in SI units to compute its cycle-averaged flux. Then add h×=A sin(2πft) with the same amplitude and frequency.
For the single plus component the derivative-square average is A²(2πf)²/2. The tensor contraction supplies a factor two, so the flux reduces to [c³/(32πG)]A²(2πf)².
The derivative scale squared is 10⁻⁴²×(200π)²≈3.948×10⁻³⁷ s⁻². Multiplication by the given prefactor gives approximately 1.59×10⁻³ W/m².
The added quadrature cross component has the same derivative-square average. The two polarization contributions add in the energy formula, giving twice the flux, about 3.17×10⁻³ W/m². A detector may project those components differently even though both carry energy.
The single-polarization averaged flux is 1.59 mW/m²; the equal-amplitude circularly polarized wave carries 3.17 mW/m². Using a single detector output in place of both physical polarizations would not generally recover this flux.
The living notebook
Explore the model
Interactive lab · function
Energy flux is quadratic in amplitude
At fixed frequency let the lab x represent A/A0 and y a normalized flux. Reset and set both controls to Coefficient a=1 and Position x=1; predict y=1. For a fresh comparison, Reset and set both controls to a=1 and x=2; predict y=4. The comparison isolates amplitude-squared scaling.
- 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 forThe parabola illustrates one scaling law. It contains no source quadrupole, TT projection, cycle average, physical energy units, memory or radiation-reaction dynamics.
Pause & reason
Check your understanding
Your turn
Try it yourself
For a Newtonian circular binary with fixed component masses, the leading radiated power scales as a⁻⁵, strain at fixed observer orientation and distance as a⁻¹, and gravitational-wave frequency as a⁻³ᐟ². If the separation doubles, give all three ratios.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Apply the separation ratio independently to each power. The strain and radiated power do not share the same separation exponent.
Compare with the solution
The new power is 2⁻⁵=1/32 of the old power. The strain scale is 1/2, while the frequency is 2⁻³ᐟ²=1/(2√2)≈0.3536 of its previous value. The flux at the fixed observer consequently also falls by 1/32: its strain-squared factor supplies 1/4 and frequency-squared factor supplies 1/8. The check connects both descriptions without dropping the frequency change.
Keep the idea
Take it with you
- The leading waveform requires a trace-free source quadrupole, TT projection and a far-zone slow-motion regime.
- Wave energy is a quadratic averaged quantity with explicit polarization normalization.
- Backreaction, memory and cosmological propagation require information beyond a monochromatic strain snapshot.
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.