Quadrupole formula

For a slowly moving isolated source, the leading gravitational-wave luminosity iswhere $Q_{ij}$ is the trace-free mass quadrupole moment.

In units $c=1$, a slowly moving isolated source has trace-free mass quadrupole moment $Q^{ij}=\int T^{00}(x^i{x}^j-{\delta }^{ij}|\mathbf{x}{|}^2/3)d^{3}x$. The retarded linearized gravity solution and the stress-energy quadrupole identity give the displayed physical spatial strain, with $\mathbf{n}=\mathbf{x}/r$ and the transverse-traceless projector $\Lambda$. If the spatial metric is written $g_{ij}=-{\delta }_{ij}-2E_{ij}$ for metric signature $(+---)$, its trace-free Lorenz-gauge coefficient before the transverse projection is $E_{ij}=G{\ddot{Q}}_{ij}(t-r)/r$. Thus $H^{\mathrm{TT}}=2E^{\mathrm{TT}}$. Far distance alone does not permit dropping source-dependent retardation: the quadrupole approximation also requires the source size to be small compared with the radiation wavelength.

For a weakly self-gravitating bound source of size $R$, mass $M$ and nonspherical dynamical mass quadrupole moment $Q=O(M{R}^2)$, the virial theorem gives $v^2=O(GM/R)=O(ϵ)$. Its time scale is $R/v$, so $\ddot{Q}=O(M{v}^2)$. The far-zone quadrupole strain then has a spatial half-strain coefficient $E=O(GM{v}^2/r)=O({ϵ}^{2}R/r)$ in units $c=1$. The physical transverse-traceless strain is $2E^{\mathrm{TT}}$. Spherical or stationary sources can have zero radiation, so this is a characteristic scale for a changing nonspherical quadrupole, not a positive lower bound. The radiation zone also requires $r\gg R/v$.

Two equal masses $m$ orbiting their centre of mass at radius $R$ have separation $2R$. Newton's law of universal gravitation gives the displayed orbital frequency. Their mass quadrupole moment varies at twice that frequency, producing gravitational radiation with two gravitational wave polarizations. For inclination $\iota$, its leading plus and cross amplitudes are proportional to $1+{\cos }^2\iota$ and $2\cos \iota$, respectively. Face-on radiation is circularly polarized; edge-on radiation is linearly polarized.

The dimensionless ratio of gravitational radius to orbital scale measures how relativistic a bound orbit is. In an equal-mass circular binary with individual orbital radius $R$, the leading quadrupole formula gives $P=(2/5)(c^5/G)[GM/(R{c}^2){]}^5$. Compact objects permit smaller noncontact orbits than extended stars and hence larger orbital compactness; near merger the Newtonian approximation must be replaced.

For two masses $M$ at opposite points of a circular orbit, each at radius $R$ from the centre of mass, ${\Omega }^2=GM/(4R^3)$. The trace-free mass quadrupole moment oscillates at twice the orbital frequency, and its third-derivative squared contraction is $128M^2{R}^4{\Omega }^6$. The quadrupole formula gives the displayed phase-independent luminosity. Using separation $2R$ rather than individual radius $R$ avoids a factor-of-32 error.

For two equal point masses at $\pm z$ falling from rest at $\pm z_0$ under Newtonian gravity, in geometrized units one has $\ddot{z}=-m/(4z^2)$ and $\dot{z}=-\sqrt{m(1/z-1/z_0)/2}$. The only nonzero component of the second mass moment tensor is $I_{zz}=2m{z}^2$, and its third derivative is $-m^2\dot{z}/z^2$. The trace-free mass quadrupole moment and quadrupole formula give the displayed power. From infinity, integrating down to $z=z_c$ gives $E=2\sqrt{2}{m}^{9/2}/(105z_c^{7/2})$. Relative to the total initial mass $2m$ the fraction is $\sqrt{2}(m/z_c{)}^{7/2}/105$. The weak-field slow-motion approximation must be distinguished from its extrapolation to a compact endpoint.

For a point mass moving as $z=L\cos \omega t$ along one fixed axis, the trace-free mass quadrupole moment is $m{L}^2{\cos }^2\omega t\mathrm{diag}(-1/3,-1/3,2/3)$. Its radiating component has twice the source angular frequency. The quadrupole formula gives $⟨P⟩=16G{m}^2{L}^4{\omega }^6/(15c^5)$ for this contribution in the leading slow-motion weak-field approximation. Any physical driver or support has its own quadrupole contribution; this expression is for the prescribed mass alone.

For a slowly evolving circular binary star with separation $a$, reduced mass $\mu$ and orbital angular frequency $\Omega$, the leading gravitational wave losses satisfy $\dot{E}=\Omega \dot{J}$. One direct check uses the standard quadrupole angular-momentum loss lawwhere $Q$ is the mass quadrupole moment and ${ϵ}_{ijk}$ is the Levi-Civita symbol. For $\mathbf{r}=a(\cos \Omega t,\sin \Omega t,0)$, put $A=\mu a^2$; then $Q_{xx}=A(1/6+\cos (2\Omega t)/2)$, $Q_{yy}=A(1/6-\cos (2\Omega t)/2)$, $Q_{xy}=A\sin (2\Omega t)/2$, and $Q_{zz}=-A/3$. Substitution gives ${\dot{J}}_z=-32G{A}^2{\Omega }^5/(5c^5)$ and the quadrupole formula gives $\dot{E}=-32G{A}^2{\Omega }^6/(5c^5)$. Their ratio proves the balance relation. The same relation follows from the rotating wave's angular harmonic number two and frequency $2\Omega$. The additional angular-momentum flux law, not the energy flux alone, is needed to justify preservation of secular circularity.

For a Newtonian Kepler orbit, let $M=M_1+M_2$ be the total mass and $\mu =M_1{M}_2/M$ the reduced mass; put $k=GM$, $\mathbf{n}=\mathbf{r}/r$, $u=\dot{r}$ and ${\mathbf{v}}_{\perp }=\dot{\mathbf{r}}-u\mathbf{n}$, with $v^2=\dot{\mathbf{r}}\cdot \dot{\mathbf{r}}$. In the convention $q_{ij}=\mu (3r_i{r}_j-r^2{\delta }_{ij})/2$, differentiating the acceleration $\ddot{\mathbf{r}}=-k\mathbf{r}/r^3$ givesThe two terms have zero cross tensor contraction, with squared norms $6u^2$ and $72v_{\perp }^2$. Therefore ${\stackrel{...}{q}}_{ij}{\stackrel{...}{q}}_{ij}=72{\mu }^2{k}^2(v^2-11u^2/12)/r^4$. Substitution into the quadrupole formula in this convention, $L=4G{\stackrel{...}{q}}_{ij}{\stackrel{...}{q}}_{ij}/(45c^5)$, proves the luminosity. It is nonnegative since $v^2-11u^2/12=v_{\perp }^2+u^2/12$. The formula uses a weak gravitational field, slow motion and leading radiation order; orbital averaging gives a secular luminosity.

For a rigid body rotating around its third principal axis, the anisotropic part of its second mass moment tensor oscillates at angular frequency $2\Omega$. The quadrupole formula gives the displayed luminosity. Equal transverse principal moments remove this rotating quadrupole and its leading gravitational wave emission.