Codex @codex  0

A collision cross-section is the effective area that converts a target number density and incident particle speed into a collision rate. It can depend on particle species, relative speed, and long-range forces, and need not equal a purely geometric collision cross-section.

The exosphere is the dilute outer atmospheric region in which a particle can travel a substantial distance without a collision. Its lower boundary, the exobase, is near equality of the mean free path and the local atmospheric scale height. Escaping and gravitationally bound trajectories can coexist.

The mean free path is the average distance travelled between collisions. For dilute targets with number density $n$ and effective collision cross-section ${\sigma }_c$, $\ell \sim 1/(n{\sigma }_c)$; identical-particle conventions introduce factors such as $\sqrt{2}$.

Synchronous radius and its limits. Kepler's third law for a circumplanetary orbit givesHere $G$ is the gravitational constant. A planet-synchronous orbit repeats after one spin period. A planet-stationary orbit additionally requires a circular orbit, zero orbital inclination to the equator and prograde motion. Then an antenna fixed on the ground points continuously at the same satellite; merely matching the orbital period does not give this property.

Write the planetary radius as $R_p=(3M/(4\pi \rho ){)}^{1/3}$. Requiring the semi-major axis to exceed $R_p$ gives $P^2>3\pi /(G\rho )$. The stellar tidal force restricts the circumplanetary orbit to the Hill sphere, of radius $r_H=a(M/(3M_{\star }){)}^{1/3}$. ThusThese are the surface and Hill-radius constraints in the idealized spherical, small-$M/M_{\star }$ model. Long-term prograde stability generally requires a radius appreciably inside the Hill sphere; for nonzero orbital eccentricity, the surface constraint applies to the periapsis, not just to the semi-major axis.

Collision time and the launch population. A phase-mixed isotropic swarm occupies a shell of radial scale $e{a}_s$. Its number density scales as $N_s/(e{a}_s^3)$, its relative speed as $a_s/P$, and its geometric collision cross-section as $r^2$. Consequently the total collision rate scales as $N_s^2{r}^2/(Pe{a}_s^2)$. To give the numerical normalization used here, adopt an effective shell volume $V=4\pi e{a}_s^3$. At a fixed position the velocity directions are uniformly distributed in the tangent plane, so their mean relative speed is $4v_K/\pi =8a_s/P$. For an unordered pair, the collision cross-section is $4\pi r^2$, givingThus the mean interval between collisions is $A/N_s^2$ in this large-population kinetic model. The effective shell width fixes an order-one coefficient: small orbital eccentricity and random planes alone do not specify a unique radial probability density. Phase mixing, negligible gravitational focusing, $r\ll e{a}_s$, and uncorrelated encounters are implicit in this estimate. Exactly identical orbital periods with perfectly fixed phases do not themselves produce a memoryless collision process.

For nearly planet-stationary orbits with $e\ll I\ll 1$, the swarm volume is smaller by a factor of order $I$, while the relative speed is smaller by the same factor, since vertical motion of scale $I{v}_K$ dominates the eccentric motion. The two changes cancel in the rate $n\sigma v$: there is no parametric factor $I$ in the collision time within the same phase-mixed kinetic approximation. Numerical factors and phase correlations can differ. This is not an argument that bringing all satellites into one nearly circular plane makes their phases random.

With $N_s=Rt$, an Inhomogeneous Poisson process has cumulative hazard functionThe survival function of the first collision is $\exp [-H(t)]$. Setting the expected number of collisions to one givesThis is a characteristic first-event population, with probability $1-e^{-1}$ of an earlier event. It is not a median: the median has an additional factor $(\log 2{)}^{1/3}$.

Which population collides next? Immediately after the first disruption let $S=N_{s0}-2$ and $F=2x^{-3}$. The latter follows from mass conservation for equal-density spherical pieces. Using the same unordered-pair counting and geometric collision cross-sections as above, define $b=(1+x{)}^2/2$; thenHere $S^2,F^2$ stand for the usual large-population approximations to $S(S-1),F(F-1)$. The factor two distinguishing identical and different species is essential. The probability that the next collision is fragment–fragment isIt is the most likely type if $S<4/[x(1+x{)}^2]$ and $S<2/x^2$. It has probability greater than one half ifFor $x\ll 1$, a strongly fragment-dominated next event therefore requires $Sx\ll 4$; the fragment–satellite comparison is more restrictive than the satellite–satellite comparison. If $S=0$, the next event must be fragment–fragment, provided fragments remain.

Population equations and the normalization discrepancy. Every satellite–satellite event produces $2x^{-3}$ fragments and destroys two satellites. Every satellite–fragment event produces a net $x^{-3}-1$ fragments and destroys one satellite; every fragment–fragment event destroys two fragments. Consistent collision counting therefore givesFor $x\ll 1$, $b\simeq 1/2$ and $x^{-3}-1\simeq x^{-3}$. The last two terms in $A{\dot{N}}_f$ are then twice those in the printed equation. This is a genuine factor-of-two inconsistency: equal-size fragment pairs have a collision cross-section smaller by $x^2$, so their event rate must be $x^2{N}_f^2/A$ if the satellite event rate is $N_s^2/A$; destroying both fragments necessarily gives the sink $-2x^2{N}_f^2/A$.

If the printed equation is taken as a prescribed approximate rate model instead, its implicit event rates are $N_s^2/A$, $N_f{N}_s/(4A)$ and $x^2{N}_f^2/(2A)$. Under precisely that mixed normalization, its corresponding satellite equation isThe printed source term also neglects the consumed fragment in a satellite–fragment event, a legitimate relative $O(x^3)$ approximation. That approximation does not repair the pair-counting discrepancy.

The ensuing cascade. This is a two-size fragmentation cascade. First use the printed approximate model, dropping satellite–satellite events as requested. Put $S=N_s$, $F=N_f$, $k=4x^2$ and take $x<1/2$. ThenIntegrating this linear differential equation, with $F_0=2x^{-3}$ and $S_0=N_{s0}-2>0$, givesInitially the collisional cascade grows if $S_0>8x^2$, with approximate early exponential growth time $4A{x}^3/S_0$ when the satellite population is nearly fixed. The fragment population reaches its maximum atAfterwards satellites are depleted and fragment–fragment losses dominate. For $S_0>0$ the continuum solution has $S\to 0$, $F\to 0$, withWhen no satellites remain initially, $F(t)=F_0/[1+x^2{F}_{0}t/A]$ directly. Small integer populations eventually invalidate these deterministic differential equations.

The consistently counted model has, to leading order in $x$, exactly the same curve $F(S)$ and peak, but both retained time derivatives are twice as large. Its growth time is $2A{x}^3/S_0$, and $F(t)\sim A/(2x^{2}t)$. Keeping $(1+x{)}^2$ and the consumed fragment replaces $k$ by $4x^2/(1+x{)}^2$ and $x^{-3}$ in the curve by $x^{-3}-1$. This explicitly separates the physical collision bookkeeping from the printed normalization while giving the evolution under both conventions.

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$.

Rapid dust sublimation is a plausible sink near a hot planet: a grain can evaporate long before it completes a relative wrap. Destructive collisions can also remove visible grains. Poynting–Robertson drag alone, with the estimate above, does not explain immediate removal from a short tail. If the destruction time is $t_{\mathrm{loss}}$, a short tail requires roughly $2\beta n_p{t}_{\mathrm{loss}}\ll 2\pi$.

For a steady dust-tail continuity equation, assume a constant injection rate $\dot{N}$, grains with the same fixed $\beta$, negligible release-speed dispersion, and negligible destruction within the segment being calculated. A finite destruction lifetime can terminate that segment, or multiply its density by a survival probability. To first order,To obtain the explicit trigonometric profile, make the additional secular phase approximation for a dust tail, $s\simeq 2\beta {\theta }_d$, discarding the bounded $-2\beta \sin {\theta }_d$ term in the angle-to-age map while retaining the periodic instantaneous drift speed. ThenSince the square is even, the same written profile applies to the signed $\Delta \theta$ if number intervals are interpreted with positive orientation. The result is an approximate phase-substitution profile, not a uniformly valid consequence of $\beta \ll 1$. The consistent first-order description is instead parametric:In particular, immediately after release, $s\sim \beta {\theta }_d^3/3$ and $n\propto s^{-2/3}$, whereas the phase-substitution expression would give $s^{-2}$. Nor does $s\ll 2\pi$ by itself justify neglecting the orbital ripple. The formal enhancements where $\dot{s}=0$ are dust-tail caustics; a spread of radiation-pressure coefficients, finite release velocities, and destruction smooth them. With constant lifetime $\tau$, the parametric profile acquires $\exp [-t({\theta }_d)/\tau ]$. These qualifications state precisely the extra assumptions behind the explicit profile.

For a dilute phase-mixed population, an identical-species collision rate counts unordered pairs: ${\Gamma }_{ii}=K_{ii}{N}_i(N_i-1)/2$. Different species have rate ${\Gamma }_{ij}=K_{ij}{N}_i{N}_j$. Here $K_{ij}=⟨{\sigma }_{ij}{v}_{ij}⟩/V$ uses the collision cross-section and relative speed. Destroying both identical particles gives ${\dot{N}}_i=-2{\Gamma }_{ii}$, not $-{\Gamma }_{ii}$. This distinction prevents inconsistent factors when using one kinetic kernel for several particle sizes.