research32 min studyLesson 849 in Physics

Wave extraction and asymptotic observables

Track tetrad conventions, curvature-to-strain integration and finite-radius extrapolation errors.

On this page

What you’ll understand

  • State a null-tetrad and polarization convention for the outgoing curvature scalar.
  • Integrate a nonzero-frequency toy curvature signal with explicit constants.
  • Perform a stipulated inverse-radius extrapolation and identify its limits.

The metric evolved near a binary is not directly the strain recorded by a distant detector. Numerical wave extraction constructs curvature or perturbative radiation variables, identifies an asymptotic frame, and converts them into waveform observables. Each step has conventions and errors. A clean oscillation on a finite extraction sphere is useful evidence, but it does not automatically equal the waveform at future null infinity or settle the uncertainty in its integrated strain.

Declare the frame before naming a curvature component

Retain signature minus-plus-plus-plus. In a locally orthonormal asymptotic frame, choose outgoing and ingoing null legs from the time and radial directions with reciprocal square-root-two normalization. Choose the complex angular leg from the theta direction plus i times the phi direction. The null legs have mutual product minus one, and the complex angular leg has product one with its conjugate. Spin rotations and null boosts change tetrad components even when the underlying curvature is unchanged.

Use the Riemann convention with derivative order mu then nu in the upper-index curvature definition, and define the outgoing Newman-Penrose scalar with a minus Weyl contraction on two ingoing legs and two conjugate angular legs. Fix polarizations so the leading asymptotic relation is the second retarded-time derivative of h-plus minus i h-cross. A different curvature sign or angular-leg convention can flip or conjugate familiar formulas, so comparisons must translate conventions explicitly.

ℓ=et+er2,n=et−er2,m=eθ+ieφ2,Ψ4=−Cabcdnamˉbncmˉd\ell=\frac{e_t+e_r}{\sqrt2},\quad n=\frac{e_t-e_r}{\sqrt2},\quad m=\frac{e_\theta+ie_\varphi}{\sqrt2},\qquad \Psi_4=-C_{abcd}n^a\bar m^bn^c\bar m^d
The Weyl tensor uses R-rho-sigma-mu-nu equal to partial-mu Gamma-rho-nu-sigma minus partial-nu Gamma-rho-mu-sigma plus the corresponding connection products. These choices fix the outgoing leading strain relation below.↔ Scroll wide equations horizontally.

Integrate curvature with a controlled asymptotic interpretation

At future null infinity in an appropriate asymptotic frame, multiply strain and curvature scalar by areal radius before taking their leading limits. Let H denote the radius-weighted complex strain and Q the corresponding radius-weighted curvature. Q is the second derivative of H with respect to retarded time. H has length units, Q inverse-length units, and the ordinary strain is dimensionless. A nonzero-frequency sinusoid therefore acquires a negative inverse-frequency-squared factor when integrated twice.

Two integration constants remain, corresponding to a constant and a term linear in retarded time. Finite noisy time series introduce further difficulties: low-frequency contamination can create large secular drift even after those two constants are chosen. Frequency-domain prescriptions can reduce such artifacts but require a cutoff and assessment of physical low-frequency content. Removing every slow component without analysis can also remove genuine signal information, including relevant memory contributions.

H=lim⁡r→∞r(h+−ih×),Q=lim⁡r→∞rΨ4=d2Hdu2,Hω=−Qωω2(ω≠0)H=\lim_{r\to\infty}r(h_+-ih_\times),\qquad Q=\lim_{r\to\infty}r\Psi_4=\frac{d^2H}{du^2},\qquad H_\omega=-\frac{Q_\omega}{\omega^2}\quad(\omega\ne0)
The relation assumes the declared asymptotic polarization and retarded-time frame. It does not hold without finite-radius and gauge corrections on every coordinate sphere.↔ Scroll wide equations horizontally.

Keep extraction error separate from grid error

Finite-radius extraction carries near-zone, tetrad and coordinate effects. Comparing spheres should use a consistent retarded-time prescription and a radius appropriate to the assumed expansion, often an areal rather than a raw coordinate radius. Extrapolating in inverse radius is a model for asymptotic behavior; its order, sphere selection and gauge dependence require checks. More points in a fit do not automatically remove an incorrect asymptotic ansatz.

Characteristic extraction can instead propagate data on null hypersurfaces toward future null infinity, introducing its own worldtube, initial-data and numerical errors. Either approach still needs resolution studies, mode decomposition conventions and treatment of startup radiation. The two-radius arithmetic below tests one stipulated leading correction, not a proof that a real waveform has only that correction. A converged finite-radius answer can remain systematically different from the desired asymptotic observable.

Q(r,u)=Q∞(u)+a(u)r+O(r−2),Q∞≈2Q(2r,u)−Q(r,u)Q(r,u)=Q_\infty(u)+\frac{a(u)}r+O(r^{-2}),\qquad Q_\infty\approx2Q(2r,u)-Q(r,u)
The two-radius combination cancels a leading inverse-radius term at a matched retarded time. Its remaining error, frame dependence and validity must be assessed separately.↔ Scroll wide equations horizontally.

A periodic curvature-to-strain calculation

For an explicitly periodic toy asymptotic signal take Q(u)=(0.008/L)cos⁡(0.2u/L)Q(u)=(0.008/L)\cos(0.2u/L). Choose the constant and linear integration terms to vanish. Find H and the leading strain at areal radius 100L100L. State the role of the negative sign.

  1. The angular frequency is 0.2/L0.2/L, whose square is 0.04/L20.04/L^2. The amplitude of H is therefore (0.008/L)/(0.04/L2)=0.2L(0.008/L)/(0.04/L^2)=0.2L.

  2. Two derivatives of a cosine produce minus frequency squared, so H(u)=−0.2Lcos⁡(0.2u/L)H(u)=-0.2L\cos(0.2u/L). The negative sign is a phase difference of pi in the declared convention, not a negative oscillation magnitude.

  3. Dividing by radius gives leading complex strain h=−0.002cos⁡(0.2u/L)h=-0.002\cos(0.2u/L), real for this chosen polarization. The zero integration terms define this periodic toy solution; they are not automatically determined for an arbitrary finite numerical time series.

The result

The radius-weighted strain is minus 0.2 L times the cosine, and its leading strain amplitude at radius 100 L is 0.002.

Explore the model

Interactive lab

Interactive lab · waves

A phase convention can change a sum

Reset and set phase difference zero and amplitude ratio one; predict resultant amplitude two. Compare a fresh Reset run at phase 180 degrees and amplitude ratio one; predict cancellation. Both components have the same fixed frequency and propagation form. The lab does not extract curvature, integrate strain, vary frequency, or extrapolate a waveform to null infinity.

A phase convention can change a sumTwo waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). Current readings: Sum at x = 0: 0.866; Resultant amplitude: 1.732; Phase difference: 60°. All parameters are adjustable with the labeled controls below.−3.14−1.5701.573.14−3.4−1.701.73.4x · radiansAmplitude
SumFirst componentSecond component
Sum at x = 0
0.866
Resultant amplitude
1.732
Phase difference
60°
t = 0 s

Two waves share the same frequency and wave number. Playback advances t in sin(2x − t) + a sin(2x − t + phase). 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 forRelative phase conventions affect waveform comparisons; the simple interference model supplies no asymptotic-extraction calculation.

Check your understanding

Why can a well-converged finite-radius curvature waveform still differ from the desired asymptotic strain?

Try it yourself

At the same prescribed retarded time, a toy radius-weighted scalar is Q(50L)=0.012/LQ(50L)=0.012/L and Q(100L)=0.010/LQ(100L)=0.010/L. Assume exactly Q(r)=Q∞+a/rQ(r)=Q_\infty+a/r for this exercise. Find Q-infinity and a, then state why two radii cannot verify the assumption.

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

Give me a hint

Subtract the two equations to isolate a. The constant asymptotic term is twice the farther value minus the nearer value.

Compare with the solution

The combination gives Q∞=(2×0.010−0.012)/L=0.008/LQ_\infty=(2\times0.010-0.012)/L=0.008/L. Using the nearer value gives a=(0.012−0.008)(50)=0.2a=(0.012-0.008)(50)=0.2, dimensionless in these units. Two values determine the two stipulated coefficients exactly but leave no independent test of an omitted inverse-radius-squared term or a bad time/frame match. Additional radii and alternative extraction controls are needed for a physical uncertainty estimate.

Take it with you

  • Tetrad, curvature and polarization conventions belong in a waveform definition.
  • Curvature-to-strain integration needs constants and low-frequency error control.
  • Grid convergence and asymptotic extraction are different limits.

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.