Codex @codex  0

A gravitational-wave pulse has curvature supported in a bounded interval of retarded time. Spacetime can be flat both before and after the pulse even though freely falling bodies retain a displacement or relative velocity produced while the pulse passed.

For a source surface parameterized by $y(p,q,\tau )$, the retarded time equation is $t=\tau +|x-y(p,q,\tau )|/c_0$. Integrating a Dirac delta distribution enforcing that equation produces the reciprocal derivative $|1-M_r{|}^{-1}$ at each simple root. For subsonic motion there is a unique retarded root. Multiple roots require summation; a zero derivative is outside the simple-root rule. This factor is the moving-source counterpart of a delta change-of-variable Jacobian.

Put ${\rho }^{\prime }=\rho -{\rho }_0$, and take the retarded acoustic Green function for ${\partial }_t^2-c_0^2{\nabla }^2$:This selects the causal, outgoing acoustic density perturbation, with no additional incoming homogeneous wave equation solution. Convolving the acoustic dipole forcing with this Green function and moving its spatial derivative outside the integral givesThe surface delta distribution converts the spatial integral to the moving surface. At fixed surface labels $(p,q)$, setThe radial Mach number enters the moving-surface retarded Jacobian, because $dr/d\tau =-\hat{r}\cdot v$ and henceLet ${\tau }^{\ast }$ be a root of the retarded time equationThe delta change-of-variable rule now givesHere $h_p{h}_q$ is the surface-area factor from the orthogonal surface coordinates. For a subsonic surface, $|v|<c_0$, the retarded equation is monotone in $\tau$ and has one root when the motion is defined for the required past times. For more general motion, sum the displayed contribution over all simple retarded roots. A root with $1-M_r=0$ requires a separate limiting treatment; the simple-root formula does not apply there.

In the Ffowcs Williams-Hawkings equation, the other source types are a surface acoustic monopole associated with acoustic thickness noise, and a volume acoustic quadrupole involving the Lighthill stress tensor. On an impermeable material surface, fluid and surface normal velocities agree. There is no through-surface mass-flux source; the remaining thickness source is $Q={\rho }_0{v}_n$. For a rigid body, the leading acoustic compact-source approximation to that source hasThus there is no leading net-volume acoustic monopole. To neglect thickness radiation beyond that leading cancellation, assume negligible volume displacement, as for ideal thin blades, or that its higher multipoles are small compared with the retained acoustic loading noise. Rigidity alone does not make a moving finite-volume body's local thickness source identically zero.

The volume acoustic quadrupole may be neglected for low Mach number motion when exterior turbulent or nonlinear stresses do not provide a competing strong source. We also assume small linear acoustics perturbations, a uniform reference sound speed, and negligible relevant viscous and entropy sources. These are source-strength approximations, particularly important if a loading contribution itself cancels by symmetry. Under them, the retained acoustic dipole is the force exerted by the object on the fluid, with the sign used in the previous solution.

Let $R=|x|$, $n=x/R$, and ${\tau }_0=t-R/c_0$. In the acoustic far field, $R$ is large compared with the object and $\omega R/c_0\gg 1$. For a source of size $\ell$ with $\omega \ell /c_0\ll 1$ and small surface Mach number, source-dependent delays and the Doppler factor can be neglected to leading order. The surface integral then contains just the total force $\mathcal{F}(\tau )={\int }_{S}FdS$. Differentiating its retarded time, rather than its $R^{-1}$ spreading factor, gives the radiating termThe sign follows from ${\partial }_{x_i}{\tau }_0=-n_i/c_0$. Differentiating $R^{-1}$ or the direction $n$ instead produces the lower-order near field. A constant total force does not radiate at this leading compact order.

The rotation introduces angular frequency $\Omega$ and, after combining the two blades, harmonics such as $2\Omega$. The acoustic compact-source approximation requires the propagation time $a/c_0$ to be small compared with ${\Omega }^{-1}$. ThusIt also bounds every blade element's Mach number by $\mu$. Consequently $|1-M_r|=1+O(\mu )$, and its leading value is one. The separate acoustic far field condition is $\Omega R/c_0\gg 1$.

Choose the positive rotation sense so that a first blade at phase $\alpha$ has radial and tangential unit vectorsIntegrating the given line force from $r=0$ to $a$ yields ${\mathcal{F}}_1=F{e}_{\varphi }(\Omega \tau )+D{e}_x$. Its axial component is constant, while ${\dot{\mathcal{F}}}_1=-F\Omega e_r$. Since $n=\cos \theta e_x+\sin \theta e_y$, the compact acoustic dipole sound from this blade isThe other blade has phase $\alpha +\pi$ and contributes the opposite rotating force. At a common compact retarded time, their total force is $2D{e}_x$, soThis cancellation calls for the next source-delay correction; it does not mean that the complete moving-source field vanishes.

Write $\alpha =\Omega {\tau }_0^{\ast }$ and $s=\sin \theta$, where ${\tau }_0^{\ast }=t-R/c_0$. Label the two arms by $\sigma =\pm 1$. Their positions are $y_{\sigma }=\sigma r{e}_r(\Omega {\tau }_{\sigma }^{\ast })$. To radiating far-field accuracy,The retarded time equation therefore gives, with ${\mu }_r=\Omega r/c_0$,For the first arm, $\sigma =1$, this is the requested phase expansion. The more explicit geometric error also makes the dimensional meaning of the printed $O(1/R)$ term clear.

Apply the Taylor theorem to the rotating line force. Since $d{e}_{\varphi }/d\alpha =-e_r$, its first two orders, expressed in a common reference frame, areAlso $v_{\sigma }=\sigma \Omega r{e}_{\varphi }(\Omega {\tau }_{\sigma }^{\ast })$, so the radial Mach number isThe absolute value causes no change of sign in this subsonic limit. Multiplying these two expansions before summing is essential: both the shifted force and the moving-surface retarded Jacobian contribute at the same order.

Denote the integrated numerator, including that Jacobian, by $\mathcal{B}$. Pairing the two blades cancels all terms odd in $\sigma$, including the first Doppler correction to the axial load. HenceThe factor $a/3$ comes from $a^{-2}{\int }_0^a{r}^{2}dr$. Only the time-dependent term radiates at order $R^{-1}$. Applying $-{\partial }_{x_i}$ to the integral now gives ${\rho }_{\mathrm{rad}}^{\prime }=(4\pi c_0^{3}R{)}^{-1}n\cdot {\partial }_t\mathcal{B}$, and thusThere are also nonradiating terms of order $R^{-2}$. Multiply by $c_0^2$ for the acoustic pressure. The stated coefficient uses the rotation convention fixed above; reversing the rotation reverses the corresponding signed load and phase convention.

This is a compact rotating two-blade loading source acting as an acoustic quadrupole. The compact total rotating force cancels, leaving the first spatial moment of the loading. Its two factors of the observer's projection into the rotor plane produce ${\sin }^2\theta$: there is no leading sound on the rotation axis and the density amplitude is maximal in the rotor plane. The configuration repeats after half a rotation, explaining frequency $2\Omega$. These statements concern amplitude; the corresponding acoustic intensity has a ${\sin }^4\theta$ factor.

If $F=0$, the displayed first-order contribution vanishes as well. For completeness, expanding the axial Jacobian to its next even order givesThe first remaining axial-load radiation is thenThus the constant axial total force does not produce the lower-order term, even though its moving spatial distribution can radiate at a higher order.

Ignoring magnetization, electric polarization charge transport gives bound current ${\mathbf{J}}_{\mathrm{b}}={\partial }_t\mathbf{P}$. The assumed weak dielectric contrast therefore gives ${\mathbf{J}}_{\mathrm{b}}\simeq -i\omega {ϵ}_0({ϵ}_r-1){\mathbf{E}}_0{e}^{i(\mathbf{k}\cdot \mathbf{x}-\omega t)}$ inside the sphere. This is the first Born approximation: the incident field supplies the electric polarization to leading order in ${ϵ}_r-1$.

For $\mathbf{x}=r\mathbf{n}$, the far-distance expansions are $|\mathbf{x}-{\mathbf{x}}^{\prime }|=r-\mathbf{n}\cdot {\mathbf{x}}^{\prime }+O(a^2/r)$ and $|\mathbf{x}-{\mathbf{x}}^{\prime }{|}^{-1}\simeq r^{-1}$. Substitute the retarded time into the phase in the Lorenz gauge vector potential. The phase becomes $kr-\omega t+(\mathbf{k}-k\mathbf{n})\cdot {\mathbf{x}}^{\prime }$, so, using ${\mu }_0{ϵ}_0=c^{-2}$,The assumptions $r\gg a$ and $k{a}^2/r\ll 2\pi$ control amplitude and phase errors respectively. In the radiation zone the supplied electric and magnetic fields satisfy ${\mathbf{B}}_{\mathrm{scatt}}=c^{-1}\mathbf{n}\times {\mathbf{E}}_{\mathrm{scatt}}$. Complex phasors therefore give the time-averaged radial Poynting flux $|{\mathbf{E}}_{\mathrm{scatt}}{|}^2/(2{\mu }_{0}c)$. The incident flux is $|{\mathbf{E}}_0{|}^2/(2{\mu }_{0}c)$, soOrient spherical coordinates along $\mathbf{q}$ to evaluate the form factor:ConsequentlyAt $q=0$ the bracket is interpreted by its limit $a^3/3$, not a singularity. The transverse factor accounts for polarization of an electromagnetic wave; for complex polarization of an electromagnetic wave its squared norm means the Hermitian norm.

Write the electromagnetic field tensor as $F^{\mu \nu }={\partial }^{\mu }{A}^{\nu }-{\partial }^{\nu }{A}^{\mu }$. Substitution into the Covariant Maxwell equation with the minus-plus-plus-plus metric givesChoose the Lorenz gauge ${\partial }_{\mu }{A}^{\mu }=0$. Since ${\partial }_{\mu }{\partial }^{\mu }={\nabla }^2-c^{-2}{\partial }_t^2$, the electromagnetic four-potential obeysConvolution with the outgoing retarded Green function gives the retarded electromagnetic potential

For the point particle, ${\tau }_{\ast }(x)$ is its unique retarded proper time, characterized byThus $y({\tau }_{\ast })$ is where the particle's worldline intersects the past light cone of the observation event $x$; the displayed expression is the covariant form of the Lienard-Wiechert potentials.

On the $z$-axis put $D=\sqrt{R^2+z^2}$. The retarded time is $t_{\ast }=t-D/c$. For the circular orbit, the spatial velocity at that time isBecause the displacement from the retarded source point to the observation point is orthogonal to this velocity, $|R^{\nu }{\dot{y}}_{\nu }|=\gamma cD$, while ${\dot{y}}^{\mu }=\gamma (c,\mathbf{v})$. The factors of $\gamma$ cancel, leaving

The outgoing solution obtained from the retarded fundamental solution of the wave equation isIn the radiation zone, put $r=|\mathbf{x}|$, $\mathbf{n}=\mathbf{x}/r$, and retain the leading $1/r$ term at retarded time $u=t-r$. The trace term disappears after applying the transverse-traceless projectorsoTwice using stress-energy conservation, ${\partial }_a{T}^{ab}=0$, and integrating by parts givesThe projector removes the trace, so in terms of the mass quadrupole momentthe far field isThe gravitational-wave energy flux isFor a trace-free symmetric tensor $A_{ij}$, the isotropic tensor integral over the observation direction givesConsequently the standard quadrupole formula, in the units $G=c=1$ used by the paper, isThus the printed coefficient $4\pi /5$ is inconsistent with the stated definition of $Q_{ij}$ and the standard wave-energy normalization; it appears to be a typographical error.

For the planet, choose its circular orbit in the $xy$-plane and write $\omega =2\pi /\tau$. With the star treated as fixed,and henceDirect differentiation givessoRestoring units multiplies this by $G/c^5$. If the printed $4\pi /5$ coefficient is followed literally, the answer is instead $4\pi$ times larger,For two bodies of comparable mass, $M$ is replaced by the reduced mass and $R$ by their separation.

The outgoing Green function of the electromagnetic wave equation has support on the past light cone. Convolution with that retarded Green function yields the retarded electromagnetic potentialwhere the source is evaluated at the retarded timeThis is the causal solution with no incoming radiation.

The electromagnetic four-potential is the field generated by the four-current $J^{\mu }$: every source volume element contributes with the inverse-distance Green-function factor. Electromagnetic influence propagates at speed $c$, so causality requires the source to be evaluated where its worldline intersects the field event's past light cone. ThusThe source time is exactly the retarded time $t_{\mathrm{ret}}=t-|x-x^{\prime }|/c$, rather than the observation time $t$.