Precessing Waveforms

Misaligned spins torque the orbital plane and rotate the radiation pattern relative to an inertial observer. The resulting waveform has amplitude, phase, and polarization modulations on a precession timescale between the orbital and radiation-reaction timescales. A rotating-frame description separates this geometric mode mixing from intrinsic higher multipoles and from the radial oscillations of an eccentric orbit.

Helpful background. Einstein Equations supplies spin–curvature dynamics; Waves, Green Functions, and Geometric Optics supplies polarization and phase transport; Relativistic Two-Body Dynamics supplies the accepted spin and orbital precession; Gravitational Radiation and Waveforms supplies strain modes and frame conventions; and Computational Relativity supplies strong-field precessing waveforms.

Three timescales and the moving radiation axis

Let $\boldsymbol L$ be the orbital angular momentum and $\boldsymbol S_1,\boldsymbol S_2$ the spins, with $\boldsymbol\chi_A=\boldsymbol S_A/(Gm_A^2)$. On the conservative precession timescale, spin–orbit and spin–spin couplings generate accepted evolution equations of the form

$$\dot{\boldsymbol L} =\boldsymbol\Omega_L\times\boldsymbol L, \qquad \dot{\boldsymbol S}_A =\boldsymbol\Omega_A\times\boldsymbol S_A.$$

These equations describe changes of direction while preserving each vector’s magnitude. The inspiral also requires a dissipative angular-momentum-loss term in $\dot{\boldsymbol L}$; a cross product alone cannot shrink the orbit’s angular momentum.

During a post-Newtonian inspiral, a common hierarchy is

$$T_{\mathrm{orb}}\ll T_{\mathrm{prec}}\ll T_{\mathrm{rr}}.$$

Orbital motion creates the fast carrier, precession changes the radiation axis and viewing geometry, and radiation reaction drifts the masses’ separation and frequencies. This hierarchy can narrow or fail near merger and during transitional precession, so it is an approximation regime rather than a definition.

Apostolatos et al. 1994, pp. 6274–6297 derives the classic spin-modulated waveform under $\Omega_{\mathrm{prec}}\ll\Omega_{\mathrm{orb}}$, while Blanchet 2024, §§3.6.2–3.6.4, pp. 242–252 gives the modern PN spin and orbital-precession organization.

Inertial, coprecessing, and corotating frames

Use the source-frame distance-rescaled strain

$$H=r(h_+-ih_\times) =\sum_{\ell,m}H_{\ell m}\,{}_{-2}Y_{\ell m},$$

with unit-normalized harmonics and $e^{-i\omega t}$. Three frames serve different purposes:

• an inertial frame has fixed axes at infinity and is the frame in which a fixed observer sees modulations;
• a coprecessing frame rotates its $z$ axis with a preferred radiation or orbital axis but retains the orbital carrier;
• a corotating frame additionally removes the dominant orbital rotation, making the strongest modes vary more slowly.

These names are not interchangeable. A coprecessing axis also needs a rotation about itself. For $z$–$y$–$z$ Euler angles

$$R(t)=R_z(\alpha)R_y(\beta)R_z(\gamma),$$

the minimal-rotation condition fixes that residual freedom by

$$\dot\gamma=-\dot\alpha\cos\beta.$$

It prevents an arbitrary spin of the transverse basis from being mistaken for physical phase evolution.

Define $R$ to map coprecessing basis vectors into the inertial frame. The modes transform at fixed $\ell$ as

$$H^{\mathrm I}_{\ell m} =\sum_{m'=-\ell}^{\ell} D^\ell_{mm'}[R(t)]\, H^{\mathrm C}_{\ell m'}.$$

This equation fixes the active-rotation and index convention used here. A source dominated by $H^{\mathrm C}_{2,\pm2}$ can populate every inertial-frame $m$ through the time-dependent Wigner matrix. The total fixed-$\ell$ norm is unchanged because $D^\ell$ is unitary. Boyle et al. distinguish the source, coprecessing, corotating, and inertial descriptions and test their mode symmetries in Boyle et al. 2014, §§II–VI.

Unitarity preserves the instantaneous strain norm. To compute news or luminosity from moving-frame coefficients, also differentiate the frame:

$$\dot{\boldsymbol H}^{\mathrm I} =D^\ell\left(\dot{\boldsymbol H}^{\mathrm C} +\mathcal A^\ell\boldsymbol H^{\mathrm C}\right), \qquad \mathcal A^\ell=(D^\ell)^\dagger\dot D^\ell.$$

The matrix $\mathcal A^\ell$ is anti-Hermitian, as follows by differentiating $(D^\ell)^\dagger D^\ell=I$. Thus the physical flux uses the norm of $\dot{\boldsymbol H}^{\mathrm C}+\mathcal A^\ell\boldsymbol H^{\mathrm C}$, not the ordinary derivative of coprecessing coefficients alone. A corotating frame can make a carrier nearly constant without eliminating the radiation it carries. See Modes, Harmonics, and Conventions for the accompanying field and frame dictionary.

Modulations generated by rotation

For a fixed observer, precession changes the instantaneous inclination and polarization dyad. The same underlying coprecessing carrier therefore acquires:

• amplitude modulation as the radiation lobes sweep toward and away from the line of sight;
• phase modulation from both the changing projection and the minimal-rotation angle;
• polarization modulation as plus and cross axes rotate;
• inertial-frame sidebands separated by combinations of orbital and precession frequencies.

These are geometrical effects of a time-dependent rotation, but the rotation itself is driven by physical spin couplings. “Twisting up” a nonprecessing waveform can capture the leading geometry only when the coprecessing-frame modes and precession dynamics are accurate in the stated parameter regime. The intrinsic spin-dependent asymmetries and orbital-timescale structure found in coprecessing and co-orbital modes are not removable by a Wigner rotation Boyle et al. 2014, §§IV–VI.

In the aligned-spin limit, every spin is parallel or antiparallel to $\boldsymbol L$. Then $\boldsymbol\Omega_L\times\boldsymbol L=0$, $\beta$ is constant and may be chosen zero, and $D^\ell$ becomes diagonal up to a fixed phase convention. Inertial and coprecessing modes coincide after that constant alignment, recovering a nonprecessing waveform.

Precession, eccentricity, and higher modes

The three effects can overlap in a time-frequency representation, but they have different organizing principles:

Effect	Intrinsic clock	Modal signature	Limiting check
Spin precession	$T_{\mathrm{prec}}$	Leading Wigner mixing plus intrinsic spin-dependent coprecessing changes	Time-dependent rotation vanishes for aligned spins
Eccentricity	Radial period $T_r$	Periastron bursts and $m\Omega_\phi+n\Omega_r$ comb	Radial sidebands vanish as $e\to0$
Higher intrinsic modes	Orbital phase and strong-field scale	Nonzero coprecessing $H_{\ell m}$ beyond $(2,\pm2)$	Obeys source symmetries even in a fixed frame

Precession does not create new radiative degrees of freedom, and higher modes do not require a moving orbital plane. Eccentricity can modulate amplitude on an orbital-scale radial period even when spins are zero. A robust diagnosis compares several $\ell$ sectors in a declared frame rather than assigning every sideband to one mechanism.

All frequencies above are source-frame quantities. Cosmological redshift multiplies orbital, precession, and radiation-reaction times by the same $1+z$, so their dimensionless hierarchy is unchanged. Detailed binary equations of motion, production numerical algorithms, and observational parameter estimation are not treated here.

Common pitfalls

Calling any preferred-axis frame coprecessing. Different axis-selection rules can disagree, particularly through merger. The axis and residual rotation condition must be specified.

Interpreting mode mixing as extra radiation. A unitary rotation redistributes fixed-$\ell$ coefficients without changing their total norm. Physical differences must persist after frames are aligned.

Assuming clean timescale separation everywhere. Simple-precession intuition can fail near transitional configurations and strong-field merger. A model must expose where its multiscale approximation stops.

Exercise

Prove that the fixed-$\ell$ modal norm is invariant under the inertial-to-coprecessing rotation. Explain why this is a useful check and what happens in the aligned-spin limit.

Solution

Write the modes as vectors $\boldsymbol H^{\mathrm I}=D^\ell\boldsymbol H^{\mathrm C}$. Wigner matrices are unitary, so $(D^\ell)^\dagger D^\ell=I$. Therefore $\sum_m\lvert H^{\mathrm I}_{\ell m}\rvert^2=(\boldsymbol H^{\mathrm C})^\dagger(D^\ell)^\dagger D^\ell\boldsymbol H^{\mathrm C}=\sum_m\lvert H^{\mathrm C}_{\ell m}\rvert^2$. Failure of this equality signals an index, conjugation, or normalization error rather than physical mode amplification. For aligned spins, the tilt may be chosen $\beta=0$ and the rotation reduces to a constant or diagonal phase, so no time-dependent mixing or precessional modulation remains.

References

• Apostolatos, Theocharis A., Curt Cutler, Gerald J. Sussman, and Kip S. Thorne. “Spin-Induced Orbital Precession and Its Modulation of the Gravitational Waveforms from Merging Binaries.” Physical Review D 49 (1994) 6274–6297. DOI.
• Blanchet, Luc. “Post-Newtonian Theory for Gravitational Waves.” Living Reviews in Relativity 27 (2024) 4. DOI. Open PDF.
• Boyle, Michael, et al. “Gravitational-Wave Modes from Precessing Black-Hole Binaries.” Physical Review D 90 (2014) 124010. DOI. Open PDF.