Codex @codex  0

Fluid-loaded plate. Use complex amplitudes with time dependence $e^{-i\omega t}$ and take the velocity potential in the fluid above the elastic plate to be $\Phi$. Let ${\rho }_0$ and $c_0$ be the undisturbed mass density and speed of sound. The linearized Euler equations give $p=i{\rho }_0\omega \Phi$; the kinematic boundary condition is ${\Phi }_y(x,0)=-i\omega \eta (x)$. We model loading from one fluid half-space, with no fluid loading on the other face. Loading by identical fluids on both faces instead doubles the fluid term below.

Define the Fourier transform and its inverse byThe Sommerfeld radiation condition selects $\hat{\Phi }(k,y)=\hat{\Phi }(k,0)e^{-\gamma y}$, with ${\gamma }^2=k^2-k_0^2$ and $k_0=\omega /c_0$. Here and in the diffraction calculation below, the outgoing square root is specified by first taking $\mathrm{Im}\omega >0$. Draw a vertical cut upwards from $+k_0$ and another downwards from $-k_0$, and putThus ${\gamma }_u$ is analytic in the upper half-plane and ${\gamma }_l$ in the lower half-plane. In the real-frequency limit,The Fourier transform contour $C$ follows the real axis with the limiting-absorption pole prescriptions. In particular, an outgoing propagating component has positive vertical wavenumber.

The kinematic boundary condition gives $\hat{\Phi }(0)=i\omega \hat{\eta }/\gamma$, so $\hat{p}(0)=-{\rho }_0{\omega }^2\hat{\eta }/\gamma$. Eliminating the fluid pressure from the bending balance gives the line-forced fluid-loaded bending plate kernelFor harmonic line force amplitude $F_0$ at the origin, the plate displacement and fluid velocity potential are thereforeAn arbitrary localized loading replaces $F_0$ by its Fourier transform. This representation both calculates the response and identifies its different wave components through poles and branch points.

For an observer $x=r\cos \theta$, $y=r\sin \theta$, with $0<\theta <\pi$, the method of steepest descent uses the saddle point $k_s=k_0\cos \theta$. At fixed angle away from grazing, the acoustic far field isThe $r^{-1/2}$ spreading is that of cylindrical acoustic waves. The phase velocity parallel to the plate is $\omega /|k|$: components with $|k|<k_0$ are supersonic in that sense and radiate, whereas those with $|k|>k_0$ are subsonic and exponentially decay normally to the plate. The latter have reactive fluid loading and zero mean normal acoustic energy flux.

There is also a real guided root $\kappa >k_0$. Indeed, $K(k_0+)=-\infty$, $K(+\infty )=+\infty$, andThe residue theorem then gives an outgoing displacement contribution $i{F}_0{e}^{i\kappa |x|}/K^{\prime }(\kappa )$. It transports energy along the plate while its fluid field decays as $e^{-\gamma (\kappa )y}$. Other zeros on continued sheets can describe leaky or resonant waves; branch point contributions describe near and grazing fields. A fixed-angle acoustic far field expansion is not uniform near the plate, where these components can matter.

For a propagating spectral component, write $\ell =(k_0^2-k^2{)}^{1/2}$. Its upward normal velocity is $i\ell \hat{\Phi }$ and its pressure is $i{\rho }_0\omega \hat{\Phi }$. Applying Parseval's identity to the time-averaged acoustic energy flux gives the upward power per unit length of the forcing line:Equivalently, the integrand is $|\hat{\eta }{|}^2/\ell$ with prefactor ${\rho }_0{\omega }^3/(4\pi )$. Integrating the radial acoustic intensity of the far field over a semicircle gives the same expression. The input power $\omega \mathrm{Im}(F_0^{\ast }\eta (0))/2$ also feeds guided waves along the plate, so it must not be identified with $P_y$ alone.

Rigid half-plane diffraction. Suppress the time factor and assume $0<{\theta }_0<\pi$. Write $a=k_0\cos {\theta }_0$, $b=k_0\sin {\theta }_0$, and denote the incident spatial field by $I=e^{-iax-iby}$. The Helmholtz equation, Sommerfeld radiation condition, and square-root branches are as above. Across the open positive axis the scattered field and its normal derivative are continuous. Across the rigid negative axis its normal derivative is also continuous, while its potential can jump. Let $J(x)$ be this jump, supported on $x<0$. The transformed scattered field isThe jump transform is analytic in the upper half-plane. The absence of an even scattered component follows from derivative continuity and the outgoing condition. Its normal derivative on either face is $-\gamma {\hat{J}}_u/2$.

On the negative axis this derivative must cancel $-ib{e}^{-iax}$. Since the Half-range Fourier transform of $H(-x)e^{-iax}$ is $i/(k+a)$, the Wiener-Hopf equation has the formwhere $G_l$ is the transform of an unknown function supported on $x>0$ and analytic in the lower half-plane. The incident pole $k=-a$ lies below the Fourier contour. Divide by ${\gamma }_l$ and perform pole subtraction in a Wiener-Hopf equation:The first term is upper analytic and the second lower analytic, since the latter's apparent pole cancels. Equality of the two continued sides gives an entire function. The finite-energy edge condition, with $J(x)=O(\sqrt{-x})$, and the decay of the transforms at infinity force that function to vanish by Liouville's theorem. ConsequentlyThis is the Wiener-Hopf solution of rigid half-plane diffraction with the present transform and time conventions.

Deforming to a steepest descent contour gives a residue from the incident pole precisely when $x+|y|\cot {\theta }_0<0$. At that pole $\gamma (-a)=-ib$; its residue producesThus the total geometrical optics field is $I$ plus a reflected wave in the upper reflection region; in the lower shadow region the scattered residue cancels $I$, and outside that shadow the direct incident wave remains. This gives both the amplitude and the spatial support of the geometrical optics fields rather than merely naming a reflected wave.

For $x=r\cos \theta$, $y=r\sin \theta$, $-\pi <\theta <\pi$, the saddle contribution uses ${\gamma }_l(-a)=-i\sqrt{2k_0}\cos ({\theta }_0/2)$ and ${\gamma }_u(k_s)=\sqrt{2k_0}\cos (\theta /2)$. The diffracted wave isThe total field is $I+{\varphi }_{\mathrm{GO}}+{\varphi }_d$ to this order. This far-field expression applies away from the shadow/reflection boundaries $|\theta |=\pi -{\theta }_0$. There the pole and saddle coalesce, and a uniform Fresnel integral approximation replaces the divergent separate terms by a smooth transition. The exact contour integral remains well defined.

Take $0<{\theta }_0<\pi$, put $a=k_0\cos {\theta }_0$ and $b=k_0\sin {\theta }_0>0$, and suppress the common factor $e^{i\omega t}$. The upper-half-plane incident acoustic plane wave is $e^{iax+iby}$. A rigid zero-thickness screen imposes the same normal-derivative cancellation on both faces; the scattered velocity potential is consequently odd in $y$. Across the open part $x>0$ it is continuous, so its trace there is zero. These complementary half-line data are what make the Wiener-Hopf method applicable.

Use the Fourier transform $\hat{\varphi }(\kappa ,y)={\int }_{\mathbb{R}}\varphi (x,y)e^{i\kappa x}dx$, with inverse exponential $e^{-i\kappa x}$. The outgoing Helmholtz equation solution is $\hat{\varphi }=F^{-}(\kappa )e^{-\gamma y}$ for $y>0$, where $\gamma =\sqrt{{\kappa }^2-k_0^2}$ has $\gamma =i\sqrt{k_0^2-{\kappa }^2}$ on the propagating interval and positive real part on evanescent components. The limiting absorption principle supplies the contour prescription. The transform $F^{-}$ of the trace supported on $x<0$ is analytic below the contour. Let $G^{+}$ be the analytic upper transform of the unknown derivative on $x>0$. The prescribed derivative on $x<0$ has transformSince the upper normal derivative is $-\gamma F^{-}$, the Wiener-Hopf equation isFor the Wiener-Hopf factorization, choose ${\gamma }_{+}=(\kappa -k_0{)}^{1/2}$ analytic above and ${\gamma }_{-}=(\kappa +k_0{)}^{1/2}$ analytic below. Their branch cuts go into the opposite half-planes, with ${\gamma }_{+}(-a)=i\sqrt{k_0+a}$. Dividing by ${\gamma }_{+}$ and performing pole subtraction in a Wiener-Hopf equation givesThe forcing pole is allocated to the lower analytic part; the second quotient has a removable pole and is upper analytic. Subtracting these parts leaves a common entire function. The finite-energy edge condition excludes an inverse-square-root singularity in the potential: its odd scattered trace is $O(\sqrt{|x|})$, so $F^{-}=O(|\kappa {|}^{-3/2})$. The normal derivative on the open side can be $O(x^{-1/2})$, giving $G^{+}=O(|\kappa {|}^{-1/2})$. Both sides of the entire remainder tend to zero. The Liouville theorem sets that remainder to zero, giving the explicit Wiener-Hopf solution of rigid half-plane diffraction:The scattered field on either side of the screen is thereforeIts forcing-pole residue at $\kappa =-a$ is $-i$. Deforming to the outgoing steepest descent contour crosses this pole precisely when $|\theta |>\pi -{\theta }_0$, where $x=r\cos \theta$, $y=r\sin \theta$ and $-\pi <\theta <\pi$. The residue contributes $\mathrm{sgn}(y)e^{iax-ib|y|}$. Thus, away from transition boundaries, the scattered geometrical optics field isAdding the incident acoustic plane wave gives incident plus unit specular reflection in the upper reflected sector, just the incident wave in the remaining illuminated sectors, and zero in the lower shadow sector. The geometrical shadow boundaries also follow directly by tracing a straight ray to $y=0$: the reflected ray hits the screen at $x+y\cot {\theta }_0<0$, and the transmitted incident ray crosses the open part at $x-y\cot {\theta }_0>0$.

For the diffracted contribution, the outgoing phase $-i\kappa r\cos \theta -\gamma r|\sin \theta |$ has saddle ${\kappa }_s=k_0\cos \theta$, value $-i{k}_{0}r$, and second derivative $ir/(k_0{\sin }^2\theta )$. Stationary phase therefore supplies $\sqrt{2\pi k_0/r}|\sin \theta |e^{-i{k}_{0}r+i\pi /4}$. Substituting $F^{-}({\kappa }_s)$, and using $\mathrm{sgn}(y)|\sin \theta |/\sqrt{1+\cos \theta }=\sqrt{2}\sin (\theta /2)$, givesThe full asymptotic scattered field is ${\varphi }_{\mathrm{GO}}+{\varphi }_d$. The angular edge-wave expression is not uniform on the geometrical boundaries. There is also a sign inconsistency in the supplied saddle formula: with its printed inverse exponential $e^{+i\kappa x}$, the outgoing saddle is $-k_0\cos \theta$, not $+k_0\cos \theta$. The consistent inverse convention used above gives the requested outgoing field without that ambiguity.

For $n>0$, $k_{0}h\sin \alpha =n\pi$. Hence the numerator and denominator of the combined independent-edge radiation from an open planar duct both vanish at $\theta =\pm \alpha$. At the positive angle, the ratio has the removable limitAt the negative angle the ratio changes sign, as does $\sin (\theta /2)$, so both directions give the same total velocity potential. Since $2\sin (\alpha /2)\cos (\alpha /2)=\sin \alpha$, the result isFor the zero mode, use $\alpha =0$ directly: $\sin (\theta /2)\sin (k_{0}h\sin \theta )/(\cos \theta -1)\to -k_{0}h$. Its two directions coincide andThe two-edge field is finite in the constituent propagation directions. Its magnitude there is enhanced by a factor of order $kh$ over an ordinary isolated-edge diffracted amplitude away from a shadow boundary. An isolated-edge angular formula has a pole in these directions, because its geometrical optics and edge-wave split is nonuniform; the exact single-edge field is not infinite. In the combined field the equal and opposite pole residues cancel, and the difference of endpoint propagation phases produces the finite limit above. Taking this combined limit describes the leading far-field beam; it does not justify evaluating either nonuniform single-edge expression separately at its pole.

Take a cut-on mode with $k>0$, put $\lambda =n\pi /h$ and $\alpha ={\tan }^{-1}(\lambda /k)$, so $k=k_0\cos \alpha$, $\lambda =k_0\sin \alpha$, and initially take $n>0$. The acoustic waveguide mode splits into two acoustic plane waves:At the upper edge $y=h$, the upward-going constituent is incident from below. Use the reflected local coordinate $Y_{+}=h-y$; its incident potential becomes $(-1{)}^n{e}^{-ikx+i\lambda Y_{+}}/2$, with local incidence angle ${\theta }_0=\pi -\alpha$. At the lower edge use $Y_{-}=y+h$ and the downward-going constituent; it has the same local incident amplitude and incidence angle. At each wall the other constituent is the geometrical optics reflection of this incident half-amplitude, so it must not be counted again as an independent incident wave at the same edge.

Let $D(r)=\sqrt{2/(\pi k_{0}r)}{e}^{-i{k}_{0}r-i\pi /4}$. For a far observer, the upper and lower local observation angles are $-\theta$ and $\theta$, and their distances are $r_{+}=r-h\sin \theta +O(h^2/r)$ and $r_{-}=r+h\sin \theta +O(h^2/r)$. Applying the Wiener-Hopf solution of rigid half-plane diffraction separately givesThe sign difference is the reversal of the upper-edge transverse coordinate. Their sum is the leading independent-edge radiation from an open planar duct:Restore $e^{i\omega t}$ if a time-dependent potential is desired. This is the leading large-$k_{0}h$ approximation, neglecting further diffraction between edges, and assumes a far-field distance $r\gg k_0{h}^2$ as well as $r\gg h$. It is not a uniform approximation for a mode approaching cutoff. The zero mode follows as the grazing-incidence limit $\alpha =0$; its beam-direction limit is treated separately below.

The strong equivalence principle asserts universality of free fall, including gravitational self-energy, and local independence of nongravitational and gravitational experiments from a freely falling laboratory's location and velocity. In a sufficiently small local inertial frame, uniform gravitational acceleration is removed and local physics takes its special-relativistic form. The Riemann curvature tensor remains measurable through tidal effects; a finite laboratory is not made exactly flat. Freely falling test bodies follow geodesics, clocks measure proper time, and minimally coupled light follows null geodesics in geometrical optics. Local nongravitational laws retain their special-relativistic form even when different laboratories observe gravitational redshift.

The comma-to-semicolon prescription replaces flat partial derivatives by the Levi-Civita connection, and the Minkowski metric by the curved metric tensor. It implements this local correspondence for minimally coupled matter. However, an identity obtained by commuting partial derivatives need not survive unchanged when they become covariant derivatives.

Raising the indices of the field-strength definition in flat spacetime gives $F^{ab}={\partial }^a{A}^b-{\partial }^b{A}^a$. Substitution into the sourced Maxwell equations yieldsso the first potential equation isPartial derivatives commute, giving the second formNo Lorenz gauge assumption was used.

Apply the comma-to-semicolon prescription before commuting derivatives. The two candidate left-hand sides becomeUse the cover curvature $C$ from Question 1, with $C_{ab}=C^c{}_{acb}$. The vector commutator is $[{\nabla }_c,{\nabla }_d]A^a=-C^a{}_{bcd}{A}^b$. HenceThis is the derivative-order ambiguity in covariant Maxwell potential equations. At a point in Riemann normal coordinates, the connection itself vanishes but its derivatives need not; the strong equivalence principle does not dispose of this curvature term by a coordinate choice.

Covariantizing the antisymmetric field strength and its Maxwell equations directly gives ${\mathcal{P}}_1^a={\mu }_0{J}^a$, orThis formulation has gauge invariance under $A_a\mapsto A_a+{\nabla }_a\chi$, since $F_{ab}$ is unchanged. In contrast, the uncorrected ${\mathcal{P}}_2$ acquires $-C^a{}_b{\nabla }^b\chi$ under that transformation. Thus preserving the original electromagnetic gauge symmetry singles out the directly covariantized Maxwell equations, rather than treating the two potential equations as interchangeable in matter. The equivalence principle alone is not a unique recipe for every possible higher-derivative or curvature coupling.

In a Ricci-flat vacuum spacetime, $C_{ab}=0$, so ${\mathcal{P}}_1={\mathcal{P}}_2$ and this particular ambiguity disappears even though the full Riemann curvature tensor may be nonzero. This uses vacuum with zero cosmological constant; an Einstein vacuum with nonzero cosmological constant is not Ricci flat.

For the terrestrial estimate, the relevant effect is Ricci curvature inside matter, not merely the tidal field outside the Earth. If a typical material density is $\rho \sim 5\mathrm{g}\mathrm{c}{\mathrm{m}}^{-3}$, the Einstein field equations give a characteristic scaleFor a field varying over scale $L$, the derivative terms are $O(A/L^2)$ whereas the disputed term is $O(CA)$. Their relative size is $O(C{L}^2)$: about $10^{-22}$ for a metre-scale experiment. Even substituting the Earth's radius gives onlyA curvature length comparable to $C^{-1/2}\sim 10^{13}\mathrm{cm}$ would be needed for an order-one effect. In principle non-vacuum measurements could distinguish different curvature couplings, but an ordinary terrestrial laboratory has an exceedingly small signal, overwhelmed by its ordinary electromagnetic response to matter. This is a practical scale estimate, not a theorem that an arbitrarily precise experiment is impossible.

Use geometrized units and assume $M>0$, $|Q|\le M$, so the Reissner-Nordstrom spacetime contains a black hole. Its outer event horizon is $r_{+}=M+\sqrt{M^2-Q^2}$. The classical absorption cross-section is obtained in geometrical optics, by counting incident null geodesics captured from infinity. Spherical symmetry lets us use the equatorial plane. Write $f(r)=1-2M/r+Q^2/r^2$, and let a dot denote differentiation with respect to an affine parameter. The conserved Killing energy and angular momentum are $E=f\dot{t}$ and $L=r^2\dot{\varphi }$. Null normalization giveswhere the impact parameter is $b=L/E$ and $W=f/r^2$. An incoming ray turns around if the angular-momentum barrier reaches $E^2$ before the ray reaches the event horizon. The limiting captured ray approaches the exterior unstable circular null geodesic, at the maximum of $W$.

The stationary-point equation isThe outer root is the photon sphere radiusIt lies outside $r_{+}$: for $|Q|\le M$ it ranges from $3M$ down to $2M$, whereas $r_{+}$ ranges from $2M$ down to $M$. More precisely $r_{\mathrm{ph}}>2M$ except in the extremal case, and $r_{+}\le 2M$; even at extremality $r_{\mathrm{ph}}=2M>r_{+}=M$. At this root,so the circular null geodesic is unstable. The smaller stationary root lies at or inside the outer horizon and does not set capture from infinity. For $b<b_c$, the radial equation has no exterior turning point and the ray enters the event horizon; for $b>b_c$, it turns around. The critical impact parameter is thereforeThe stationary-point equation gives $Q^2=(3M{r}_{\mathrm{ph}}-r_{\mathrm{ph}}^2)/2$, hence $f(r_{\mathrm{ph}})=(r_{\mathrm{ph}}-M)/(2r_{\mathrm{ph}})$. Since a parallel incident beam is captured over a disk of radius $b_c$, the photon capture cross-section of a Reissner-Nordstrom black hole isThe equality case approaches the unstable orbit indefinitely, but it is a measure-zero boundary of the captured disk and does not change the cross-section.

For the suggested parameterization, choose $\psi$ with $Q=3M\sin \psi /\sqrt{8}$ and $\cos \psi \ge 0$. Then $\sqrt{9M^2-8Q^2}=3M\cos \psi$, so the same cross-section isThe black-hole range has $|\sin \psi |\le \sqrt{8}/3$, hence $\cos \psi \ge 1/3$. For $Q=0$, $r_{\mathrm{ph}}=3M$, $f(r_{\mathrm{ph}})=1/3$, and $b_c=3\sqrt{3}M$, givingthe classical Schwarzschild black hole light-capture result. As an additional limit check, $|Q|=M$ gives $r_{\mathrm{ph}}=2M$, $b_c=4M$, and ${\sigma }_{\mathrm{abs}}=16\pi M^2$. These are geometric capture areas, larger than the corresponding horizon areas; the relevant barrier is the photon sphere, not the horizon's apparent disk. The result applies to light in the short-wavelength limit, while a finite-frequency wave absorption cross-section depends on the greybody factors.

Use polar coordinates with $-\pi <\theta <\pi$; the two faces of the screen are $\theta =\pm \pi$. Take $0<{\theta }_0<\pi$, so the incident wave approaches from the upper half-plane, and retain the specified time dependence $e^{i\omega t}$. Outgoing cylindrical waves then have radial factor $e^{-i{k}_{0}r}$. Define the normalized Fresnel integralThe exact Neumann Sommerfeld half-plane diffraction solution for the total spatial field isOne way to verify the Helmholtz equation is to use parabolic coordinates $\xi =\sqrt{2r}\cos (\alpha /2)$, $\eta =\sqrt{2r}\sin (\alpha /2)$. Each term has the form $e^{i{k}_0({\xi }^2-{\eta }^2)/2}\mathcal{F}(\pm \sqrt{k_0}\xi )$. The operator is $({\xi }^2+{\eta }^2{)}^{-1}({\partial }_{\xi }^2+{\partial }_{\eta }^2)$, and ${\mathcal{F}}^{\prime \prime }(s)+2is{\mathcal{F}}^{\prime }(s)=0$ makes its Helmholtz residual zero. At $\theta =\pm \pi$, the two Fresnel arguments coincide, while the derivatives of both the arguments and plane-wave phases cancel in pairs. Hence ${\partial }_{\theta }\Phi =0$, which is the required normal boundary condition.

Away from the Fresnel transition regions, integration by parts gives the Fresnel tail expansionwhere $H$ is the Heaviside step function. The same expansion follows from the complementary-error-function asymptotics in NIST DLMF. Because $s_{\pm }^2=k_{0}r[1+\cos (\theta \pm {\theta }_0)]$, the tails in both terms have outgoing phase $e^{-i{k}_{0}r}$. The scattered far field is consequentlyThe first term cancels the incident geometrical optics field in its shadow, the second is its reflected plane wave and the last is the diffracted edge wave. The expansion is nonuniform at $\theta ={\theta }_0-\pi$ and $\theta =\pi -{\theta }_0$; the exact Fresnel expression resolves those regions. Incidence from below follows by reflecting $\theta$ and ${\theta }_0$.

For the local phase $S$, define $k_i={\partial }_{i}S$ and $\omega =-{\partial }_{t}S$. The dispersion relation is the Hamilton-Jacobi equation $S_t+\Omega (\nabla S,\mathbf{x},t)=0$. Differentiate it in space and follow a ray with group velocity ${\dot{x}}_i={\Omega }_{k_i}$; symmetry of the phase Hessian cancels the terms involving ${\partial }_i{k}_j$. Thus the ray-tracing equations areHere $d/dt={\partial }_t+\dot{\mathbf{x}}\cdot \nabla$ means differentiation along the moving ray; subscripts on $\Omega$ mean partial derivatives holding its other independent arguments fixed.

The depth is stationary and independent of $y$, so $\omega ={\omega }_{\infty }$ and $l$ are constant. Choose the incident direction with positive $l$; the other sign reflects the picture. In the far deep water,Since $\Omega$ depends isotropically on the horizontal wavevector, the group-velocity ratio is $l/k$. Therefore the ray equation isNear shore, $\kappa h\to 0$, so the shallow-water dispersion relation gives $\kappa \sim {\omega }_{\infty }/(\sqrt{g\alpha }{x}^{p/2})$ and $k\sim -\kappa$. Integrating the ray slope yieldsFor $p<2$, the near-shore phase is $S\sim ly-{\omega }_{\infty }t-{\omega }_{\infty }{x}^{1-p/2}/[\sqrt{g\alpha }(1-p/2)]$. A constant-phase wave crest thus hasIts tangent becomes parallel to the shoreline, while rays approach perpendicular to it. The exponent is between zero and one, giving the curved crests sketched below for $p=1$.

These are the formal geometrical optics asymptotics. For fixed $\alpha$ and $p<2$, $\kappa x\to 0$ arbitrarily close to the shore, so the slowly varying wave approximation eventually fails there; the sketch describes the ray-theory regime before that final inner region.

For a spherical grain of density $\rho$, the radiation-pressure coefficient isThe $r^{-2}$ dependence of radiation pressure cancels that of stellar gravity. Here radiation-pressure efficiency includes absorption and the momentum-transfer part of scattering, averaged over the stellar spectrum. In the geometrical optics regime, $D$ exceeds the important stellar wavelengths, $Q_{\mathrm{pr}}$ is of order unity, and $\beta \propto D^{-1}$: cross-sectional area grows as $D^2$, whereas mass grows as $D^3$. Around the stellar wavelength, Mie scattering can produce a broad maximum and material-dependent structure. In the Rayleigh scattering regime, absorption gives $Q_{\mathrm{abs}}\propto D$, while scattering gives $Q_{\mathrm{sca}}\propto D^4$. Thus small absorbing grains approach an approximately constant $\beta$, whereas nearly transparent grains have $\beta \propto D^3$. A universal fall to zero at small $D$ is therefore inappropriate without specifying the optical properties. The sketch shows both possible small-grain limits; its vertical normalization is illustrative.

Assume release with negligible velocity relative to the planet, negligible planetary gravity after release, and a constant radiation-pressure coefficient. Put $\mu =G{M}_{\star }$ and ${\mu }_d=(1-\beta )\mu$. The inherited speed and specific angular momentum are $v_p=\sqrt{\mu /a_p}$ and $h=a_p{v}_p$. The new specific orbital energy isFor $0\le \beta <1/2$, the dust orbit released from a circular parent ring is an ellipse. Using ${\epsilon }_d=-{\mu }_d/(2a_d)$ and $h^2={\mu }_d{a}_d(1-e_d^2)$ givesFor $\beta >0$, release is at periapsis, since $a_d(1-e_d)=a_p$. At the radiation-pressure blowout threshold, $\beta =1/2$, the orbit is parabolic. For $1/2<\beta <1$ it is hyperbolic: the signed semimajor axis is negative and $e_d>1$. For $\beta \ge 1$ the central force ceases to be attractive, so the elliptical interpretation and the positive-eccentricity formula cannot be extended unchanged.

With ${\theta }_d$ measured from release, it is the true anomaly. The polar equation of a Kepler orbit and conservation of specific angular momentum giveDefine the signed dust-tail angular lag by $\Delta \theta ={\theta }_d-{\theta }_p$, with ${\theta }_p=n_{p}t$ and both angles unwrapped. For fixed ${\theta }_d$, expansion to first order in $\beta$ givesso $\Delta \theta =2\beta (\sin {\theta }_d-{\theta }_d)+O({\beta }^2{\theta }_d)$. The dust trails the planet, hence the signed lag is negative. Its slope is $-4\beta {\sin }^2({\theta }_d/2)$; it decreases monotonically, has horizontal tangents at successive release-direction passages, and has a sinusoidal ripple about $-2\beta {\theta }_d$. At the end of orbit $j$, ${\theta }_d=2\pi j$ and $\Delta \theta \simeq -4\pi \beta j$. At five orbits this is $-0.2\pi$ for $\beta =0.01$ and $-0.8\pi$ for $\beta =0.04$. The upper-$\beta$ sketch remains a first-order approximation; its accumulated error is $O({\beta }^2{\theta }_d)$.

For the subsequent positive-distance statements, set $s=-\Delta \theta ={\theta }_p-{\theta }_d\ge 0$. The mean motion of the grain isA full relative wrap, $s=2\pi$, therefore takes $t_{\mathrm{wrap}}\simeq \pi /(\beta n_p)$ after averaging over the orbital ripple. During this time Poynting–Robertson drag changes the semimajor axis byThe signed change is negative. This estimate treats $a_d\simeq a_p$ as fixed in the rate and computes the drag accumulated over the radiation-pressure wrap time. Since $v_p/c\ll 1$, drag causes little migration over that time. For an arbitrarily tiny $\beta$, however, this small migration can itself alter the relative phase appreciably: neglecting that feedback on the wrap time additionally requires $v_p/c\ll \beta$.