Codex @codex  0

For resonant trapping of dust in an exterior $(p+q):p$ mean-motion resonance, let Poynting–Robertson drag have parameter $A$. Holding semi-major axis $a$ fixed and retaining secular rates through second order in orbital eccentricity givesFor $5p>q$, squared orbital eccentricity approaches this equilibrium exponentially in the truncated model. If the predicted value is not small, higher-order terms are needed for quantitative dynamics. Finite planetary torque also imposes a minimum eccentricity for balance, which alone does not guarantee capture during migration.

Assume the physical coefficients $Q_a,Q_b$ are positive. Differentiate the catastrophic disruption threshold:The bracket is strictly increasing from a negative value to positive infinity, so there is one critical point and it is the global minimum. In fact $Q_D^{\ast }$ tends to infinity at both endpoints $D\to 0$ and $D\to \infty$. The strength-gravity disruption transition is thereforeAt this diameter, $b{Q}_b{D}_w^b=a{Q}_a{D}_w^{-a}$, so $Q_w^{\ast }=[(a+b)/b]Q_a{D}_w^{-a}$. This minimizes the full catastrophic disruption threshold; it need not be where its two contributions are equal.

For equal-density spherical planetesimals, the specific impact energy per unit target mass is $\frac{1}{2}(D_{\mathrm{im}}/D_w{)}^3{v}_{\mathrm{rel}}^2$. Setting this to $Q_w^{\ast }$ gives the catastrophic projectile diameterHere $v_{\mathrm{rel}}$ is the impact speed, the threshold is defined per target mass, and energy-coupling changes are not included. If projectile and target densities differ, multiply the diameter by $({\rho }_{\mathrm{target}}/{\rho }_{\mathrm{projectile}}{)}^{1/3}$. For $a=1/2$, $b=3/2$,so collecting the powers givesThe result presumes that this projectile diameter is available in the population.

For the dust, $n(D)dD$ is the total number in a diameter interval, not a number per unit volume. The mass normalization of a power-law size distribution with dust mass $M/2$ givesThe approximation uses $(D_{\min }/D_{\max }{)}^{4-\alpha }\ll 1$. Merely having a small diameter ratio does not guarantee this accuracy if $\alpha$ approaches four, so retain the exact denominator in that case.

With $\gamma >0$, the stated Poynting–Robertson drag law integrates directly:Thus dust arrives in order of increasing diameter. This is burst transport under inverse-size radial drag, assuming an instantaneous release at $r_0$ and no intervening destruction or changes of diameter.

Dust that reaches the planet's orbital radius has several possible fates. It can pass without a close encounter and continue spiralling towards the star; undergo resonant trapping of dust in an exterior mean-motion resonance, where resonant torque temporarily balances the drag; be scattered onto a different stellar Kepler orbit or ejected; or collide with the planet and be accreted. Resonant trapping of dust can increase orbital eccentricity until escape or an encounter becomes possible. Temporary planet-bound capture can become permanent if drag or another dissipative process removes sufficient energy. Arrival at a radial circle alone does not guarantee a collision with the planet at the correct longitude.

A more massive planet has a larger Hill sphere, stronger resonant perturbations and stronger planetary scattering. Capture into a mean-motion resonance is favoured when drift across it is slow compared with its libration timescale. Since $|\dot{r}|\propto D^{-1}$ here, larger grains drift more slowly and are generally easier to trap; small fast-drifting grains more readily cross a resonance. Grain radiation pressure can also move resonance locations, though the prescribed nearly circular drift model omits that additional change. For direct impact, a larger physical radius increases the geometric cross-section, and gravitational focusing gives, in the two-body encounter approximation,A compact massive planet can nevertheless scatter or eject grains rather than hit them: its encounter region can be much larger than its physical radius. Slower dust migration permits more repeated encounter opportunities. These competing effects mean neither planet mass nor grain diameter alone specifies an accretion probability.

For the idealized complete-accretion calculation, put $T=t_p(D_{\max })$ and $\eta =D_{\min }/D_{\max }$. At time $t$, the arriving diameter is $D_t=D_{\max }t/T$ and $d{D}_t/dt=D_{\max }/T$. Hence the dust-burst accretion pulse follows from the arriving mass per diameter:Dropping the small lower-cutoff normalization correction gives the requested expressionThe rate is zero before the smallest grains arrive and after the largest arrive. Integrating the exact expression over its support returns $M/2$, as required by mass conservation. The approximate expression has that integral to the same lower-cutoff accuracy. This pulse assumes no resonant delay or scattering loss before accretion.

For the surface density of a size-sorted dust burst, take $0<r<r_0$ and writeAt fixed $t>0$, the size-radius map is monotone and $d{D}_r/dr=2r{D}_r/H_r$. Define $\Sigma (r,t)$ as the azimuthally averaged surface density, so the mass in an annulus is $2\pi r\Sigma dr$. Transforming the power-law size distribution into radial mass givesIt follows thatThe requested leading approximation omits $1-{\eta }^{4-\alpha }$ from the denominator. The surface density is zero where the inferred $D_r$ lies outside the original size interval. Equivalently, one can obtain it from the instantaneous inward mass flux divided by $2\pi r|\dot{r}(D_r)|$. At $t=0$ the idealized initial distribution is a delta-function ring at $r_0$; the formula above describes the subsequent radial spread. No uniform initial azimuthal distribution is required for this azimuthal average, but the local density of a clump would require its angular evolution. If the complete-absorption planet assumption is retained, this expression holds exterior to $r_p$ and $\Sigma =0$ interior to the planet.

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.

In resonant trapping of dust, the resonant torque must compensate the inward Poynting–Robertson drag. At the nearly fixed resonant semi-major axis, ${\dot{a}}_{\mathrm{res}}+{\dot{a}}_{\mathrm{pr}}=0$ requiresSince $|\sin \varphi |\le 1$, a necessary torque-balance condition is $|C|e^q\ge A(1+3e^2)/[a^2(p+q)]$. When the limiting orbital eccentricity is small, this givesEquality allows only the extreme phase $|\sin \varphi |=1$, with no margin for restoring torque; a robust trapped state needs $e>e_{\min }$. This is a necessary condition for resonance capture, not a sufficient probability-of-capture criterion for a migrating particle.

During resonant trapping of dust, the mean semi-major axis is stationary. Setting the supplied ${\dot{a}}_d$ to zero givesA real resonant phase exists only when $|\sin \varphi |\le 1$. Hencewhere $e_{d0}$ is the grain's orbital eccentricity when trapping begins. The inequality also has the expected sign $\sin \varphi <0$, allowing the planetary torque to oppose Poynting–Robertson drag.