Past exam of the mathematics course of the University of Cambridge / 2017 / iii / Paper 316 / 2

Treat $a,e$ as instantaneous osculating orbital elements of the fixed-central-parameter Kepler orbit. The perturbing work changes the specific orbital energy according toThe specific angular momentum and conic equation giveSubstituting proves the Gauss planetary equation for the semi-major axis:Here $f$ denotes true anomaly, not the fractional luminosity used in Question 1. The equation follows from work alone; a normal perturbing acceleration has no instantaneous contribution because the velocity lies in the orbital plane. For slowly averaged evolution, the perturbation must be weak over one orbital period.

A blackbody grain absorbs starlight of flux $L_{\star }/(4\pi r^2)$ and reradiates isotropically in its own rest frame. Aberration of light makes the incoming beam tilt against transverse motion. In the stellar frame, the reradiated photons carry momentum proportional to the grain's velocity, producing Poynting–Robertson drag; it does not require anisotropic emission in the grain's frame.

Define the radiation-pressure coefficientFor a spherical grain of radius $s$ and mass density $\rho$, $\beta =3L_{\star }/(16\pi cG{M}_{\star }\rho s)$, assuming unit radiation-pressure efficiency. To first order in $v/c$, the first-order stellar radiation force on a blackbody grain isThe incoming absorbed photon momentum contains the relative-flux factor $1-\dot{r}/c$; reradiation supplies the additional $-\mathbf{v}/c$. Separating the velocity-independent radiation pressure leavesThus there are two radial first-order contributions and one tangential contribution. The approximation assumes $v\ll c$, steady absorption/reradiation and constant grain properties; a non-blackbody requires its actual optical efficiencies.

Put $A=\beta \mu /c$. The Poynting–Robertson drag acceleration has $\bar{R}=-2A\dot{r}/r^2$ and $\bar{T}=-A{v}_{\theta }/r^2$. Its work is strictly dissipative:For an unperturbed elliptic orbit, use $dt=r^{2}df/h$, $\dot{r}=(\mu /h)e\sin f$ and $v_{\theta }=(\mu /h)(1+e\cos f)$. Averaging over its orbital period $P$ givesSince $\dot{a}=2a^2\dot{\epsilon }/\mu$, $P=2\pi a^{3/2}/\sqrt{\mu }$ and $h^3={\mu }^{3/2}{a}^{3/2}(1-e^2{)}^{3/2}$, the Poynting–Robertson decay of a circumstellar orbit isFor a circular Kepler orbit this becomes $-2A/a$. This is a secular first-order average, not an exact instantaneous decay law. If static radiation pressure is appreciable, use a bound conservative orbit with central parameter ${\mu }_{\mathrm{eff}}=(1-\beta )\mu$, requiring $0\le \beta <1$, and define its orbital elements accordingly; repeating the work calculation cancels ${\mu }_{\mathrm{eff}}$ and leaves the same coefficient $A=\beta G{M}_{\star }/c$.

Let ${\mu }_p=G{M}_p$, $A=\beta G{M}_{\star }/c$ and $\kappa =A/a_p^2$. A tightly bound circumplanetary orbit lies well inside the Hill sphere, so $a/a_p\ll 1$ and its orbital period is short compared with the planet's year. Treat the stellar flux and stellar direction $\hat{\mathbf{R}}$ as constant during one dust orbit. With stellar-frame velocity ${\mathbf{V}}_p+\mathbf{u}$, the velocity-dependent acceleration isThe terms independent of $\mathbf{u}$, including the leading static radiation pressure, do no net work on an unperturbed closed circular orbit. HenceFor a coplanar circular orbit, $u^2={\mu }_p/a$ and the component along the stellar direction has mean square $u^2/2$. Thus $\dot{a}=2a^2{\dot{\epsilon }}_p/{\mu }_p$ yields the Poynting–Robertson decay of a circumplanetary orbitThe coefficient three assumes coplanarity, which the PDF does not state. For orbital normal $\hat{\mathbf{k}}$, the general circular covariance is $⟨u_i{u}_j⟩=(u^2/2)({\delta }_{ij}-k_i{k}_j)$, givingAn orbital plane normal to the stellar radial direction has coefficient two during that orbit, a counterexample to an orientation-independent coefficient three. Averaging also over the planet's circular stellar orbit, with fixed dust-plane orbital inclination $I$ to it, gives coefficient $3-\frac{1}{2}{\sin }^{2}I$. The conservative forces must be weak enough for the assumed approximately circular, planet-bound orbit to persist; the small planet-to-star mass ratio alone does not ensure this for arbitrary grain $\beta$.

Keep $A=\beta G{M}_{\star }/c$ constant. A circular circumstellar orbit stays circular in the secular drag approximation, and $da/dt=-2A/a$ integrates to $a^2=a_p^2-4At$. The inspiral time under Poynting–Robertson drag to the stellar surface isThe final expression treats the star as a point, or assumes $R_{\star }\ll a_p$.

For the coplanar circumplanetary orbit used in the preceding result, $da/dt=-3Aa/a_p^2$, so $a=a_0\exp (-3At/a_p^2)$. The inspiral time of circumplanetary dust to the planet's surface, for $a_0>R_p$, isThus the timescales have the same dependence on stellar flux and $\beta$, but planetary arrival contains a logarithm of the initial planetocentric radius. It is not always shorter: for a point star, $t_{\mathrm{cp}}<t_{\mathrm{cs}}$ only if $a_0/R_p<e^{3/4}$. A constant tilted-orbit average replaces three by its appropriate orientation coefficient. Sublimation or other removal can terminate the evolution before either idealized arrival time.

Radiative drag is only one of several loss mechanisms. Around a star, radiation-pressure blowout can eject small fragments; its $\beta =1/2$ threshold applies specifically to zero-kick release from a circular parent orbit. Stellar-wind drag and gas drag can drive planetary migration, while sublimation destroys grains approaching high-temperature regions. Collisional cascades destroy or fragment grains and can feed the unbound size range. Planetary scattering can cause ejection, collision with a planet, or a stellar impact; resonant trapping of dust can instead delay planetary migration.

For circumplanetary orbits, collisions with the planet or its satellites, disruption in collisions, and escape under stellar tidal forces are additional losses. Orbits near or outside the Hill sphere need not remain planet-bound. Radiation pressure on circumplanetary dust can excite planetocentric orbital eccentricity or unbind very small grains; it need not act only through slow Poynting–Robertson drag. For charged grains, the Lorentz force in stellar or planetary magnetic fields can alter or destabilize an orbit. Shadowing of circumplanetary dust changes the radiation-force average and can reduce the quoted decay rate. Which mechanism dominates depends on grain size and composition, environment, orbit orientation and the available collision or gas density.