research33 min studyLesson 834 in Physics

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

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.

Iij=∫ρ(xixj−13δijr′2)d3x,hijTT(t,R)=2Gc4RΛij,kl(R^)I¨kl(t−R/c)I_{ij}=\int\rho\left(x_ix_j-\frac13\delta_{ij}r'^2\right)d^3x,\qquad h^{\mathrm{TT}}_{ij}(t,\mathbf R)=\frac{2G}{c^4R}\Lambda_{ij,kl}(\widehat{\mathbf R})\ddot I_{kl}(t-R/c)
I has units kg m², R is observer distance, r′ is source radius and Λ is the spatial transverse-traceless projector. h is dimensionless. The leading formula assumes retarded outgoing propagation, weak/slow source dynamics in the basic derivation, and a far radiation zone; it omits higher multipoles and propagation corrections.↔ Scroll wide equations horizontally.

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.

F=c332πG⟨h˙ijTTh˙ijTT⟩=c316πG⟨h˙+2+h˙×2⟩,P=G5c5⟨I...ijI...ij⟩F=\frac{c^3}{32\pi G}\left\langle\dot h^{\mathrm{TT}}_{ij}\dot h^{\mathrm{TT}}_{ij}\right\rangle=\frac{c^3}{16\pi G}\left\langle\dot h_+^2+\dot h_\times^2\right\rangle,\qquad P=\frac{G}{5c^5}\left\langle\dddot I_{ij}\dddot I_{ij}\right\rangle
F is energy flux in W/m² and P total radiated power in W. The strain derivative is in s⁻¹ and angular brackets indicate the specified averaging. These are leading quadrupole/high-frequency energy formulas, not instantaneous gauge-invariant local densities or a detector's projected strain squared.↔ Scroll wide equations horizontally.

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.

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.

  1. 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)².

  2. The derivative scale squared is 10⁻⁴²×(200π)²≈3.948×10⁻³⁷ s⁻². Multiplication by the given prefactor gives approximately 1.59×10⁻³ W/m².

  3. 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 result

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.

Explore the model

Interactive lab

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.

Energy flux is quadratic in amplitudeThis 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 forThe parabola illustrates one scaling law. It contains no source quadrupole, TT projection, cycle average, physical energy units, memory or radiation-reaction dynamics.

Check your understanding

Why does the polarization form of the flux have coefficient c³/(16πG) rather than c³/(32πG)?

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.

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.