Codex @codex  0

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

A directed stream of molecules used to prepare samples for spectroscopy.

A molecular beam also includes a stream of chemical atoms, then called an atomic beam. An oven supplies particles; aligned slits restrict their directions, and a vacuum limits subsequent collisions. Narrowing the directions does not by itself make every particle's speed equal. This gives a useful separation of tasks: prepare a path, select or measure speed, apply an interaction, and read the resulting particle counts or deflection.

The advantage is unusually tangible. A magnetic field can separate spin components in the Stern-Gerlach experiment; even gravity can sort a beam by flight time. The effusive beam speed distribution explains which particles leave the source, while gravity deflection of an atomic beam explains how that distribution becomes a pattern on a detector.

Momentum conservation states that total momentum stays constant when the net force exerted from outside the chosen collection is zero. Use a closed collection, meaning the same fixed-mass point masses are followed throughout, and an inertial frame. ThenTo see why, sum Newton's second law over the members. If ${\mathbf{F}}_{ij}$ is the force on member $i$ from member $j$, Newton's third law assumes ${\mathbf{F}}_{ij}=-{\mathbf{F}}_{ji}$. Each pair cancels in the total:The internal forces act on different bodies and can change their individual momenta; their contributions cancel only in the collection's sum. An isolated system has no external interaction, but balanced external forces also suffice for momentum conservation. A closed system merely excluding matter exchange need not have conserved momentum. ^[ref]

The usefulness is that the detailed internal forces during a collision need not be known. Both an elastic collision and a perfectly inelastic collision obey the same isolated momentum balance; kinetic energy conservation is an additional condition for the elastic collision. If the external impulse is small but nonzero, the total momentum change is correspondingly small. For a variable-mass system exchanging matter, the entering or leaving momentum must also be included.

In a translation-invariant field theory, integrating over the spacetime zero mode enforces momentum conservation at interactions and scattering amplitudes.

The momentum relaxation time is the exponential relaxation time of mean directed momentum after the driving electric field is removed. In the simplest Drude model, independent collisions occur with probability $dt/\tau$ during a short time interval $dt$ and reset the mean velocity of the collided population to zero. The remaining fraction preserves its previous mean momentum. To first order in $dt$, averaging gives $d⟨\mathbf{p}⟩/dt=q\mathbf{E}-⟨\mathbf{p}⟩/\tau$. This is a population equation: it does not assert that every particle stops at equally spaced instants.^[ref]

For real materials, a collision that changes a particle's direction only slightly need not erase its directed momentum, so the momentum relaxation time need not equal the time between all scattering events. This transport time is also distinct from the bulk charge relaxation time $ϵ/\sigma$: one describes loss of directed motion, the other redistribution of charge density.

In 2001, Andrew G. Truscott, Randall G. Hulet, and collaborators observed degeneracy pressure in trapped lithium. A Zeeman slower fed a dual magneto-optical trap; optical pumping prepared atoms for a cloverleaf apparatus for magnetic trapping of neutral atoms. Microwave evaporation removed energetic lithium-7 bosons; elastic collisions cooled lithium-6 fermions by sympathetic cooling.^[ref]

Sequential absorption imaging recorded both trapped clouds on a camera. Fits to axial profiles used the boson cloud as a thermometer. At common temperature, the nearly equal magnetic trapping potentials provided a comparison: the clouds had similar sizes near $T/T_F=1$, but the fermion cloud remained broader as cooling reached $T/T_F=0.25$. Figure 3 compared its radius with ideal Fermi gas and classical predictions.^[ref]

Varying the wait before imaging tested thermal equilibrium. Finite imaging resolution restricted the analysis to the axial direction; fits neglected interaction-induced density shifts. Cooling eventually suffered from the dwindling boson heat capacity. The result supports quantum occupation effects, without directly resolving individual occupied states or reaching zero temperature.^[ref]

These were neutral alkali atoms in a spatially varying trap. The Fermi energy of a model electron gas instead uses uniform electrons, two equally populated spin states, and the electron mass. The measured trapped-atom radii are not measurements of that electron model's $p_F$ or $E_F$; comparing their numerical Fermi energies requires the appropriate particle mass, internal-state populations, and confining potential.

Orbit crossing occurs when two geometrical orbital paths intersect. It permits collisions or close encounters when bodies also reach the intersection at compatible times; crossing alone does not imply a collision. For nearly circular coplanar nested orbits, a useful local test follows by differentiating their radial positions at fixed longitude.

A plasma contains mobile charged species whose collective electromagnetic response matters. On scales large compared with the Debye length, positive and negative charge densities nearly cancel: this is quasineutrality, not the absence of electric currents or electric fields. For electrons of temperature $T_e$ and number density $n_e$, the characteristic plasma screening and oscillation scales areA collective plasma description requires many particles within a screening volume; long-wavelength, sufficiently slow behavior can often be described by magnetohydrodynamics with an appropriate fluid closure.

An individual particle feels the Lorentz force $q(\mathbf{E}+\mathbf{v}\times \mathbf{B})$. In a uniform magnetic field, its transverse motion gyrates while its parallel motion remains free, giving a helical path when both components are present. The nonrelativistic cyclotron frequency and gyroradius areWhen the system scale is much larger than $r_L$ and the relevant time scale longer than ${\Omega }_c^{-1}$, the field strongly organizes transport across it, provided collisions do not interrupt the gyromotion. A perpendicular electric field also gives the common guiding-centre drift ${\mathbf{v}}_E=\mathbf{E}\times \mathbf{B}/B^2$.

For a conducting fluid, the macroscopic Lorentz force density has the useful decompositionThese terms are magnetic pressure and magnetic tension. Magnetic tension resists bending, while magnetic-pressure gradients resist nonuniform compression. The plasma beta $\beta =2{\mu }_{0}p/B^2$ compares thermal and magnetic stresses.

With uniform magnetic diffusivity ${\eta }_m$, the resistive induction equation isIf the magnetic Reynolds number $R_m=UL/{\eta }_m\gg 1$, ideal magnetohydrodynamics applies approximately and magnetic flux freezing links the field geometry to fluid motion. Finite diffusivity permits field slippage and magnetic reconnection. Field tension supports Alfvén waves with $v_A=B/\sqrt{{\mu }_0\rho }$; coupled magnetic and gas compression also supports magnetosonic waves. These effects allow a field to transmit forces and angular momentum between distant parts of an astrophysical plasma.

Use Gaussian units, plane waves proportional to $e^{i\mathbf{k}\cdot \mathbf{r}-i\omega t}$, and the signed cyclotron frequency ${\Omega }_s=q_s{B}_0/(m_{s}c)$. Write ${\Omega }_e=e{B}_0/(m_{e}c)>0$ for its electron magnitude, so the signed electron value is $-{\Omega }_e$. Define the species plasma frequency by ${\omega }_{ps}^2=4\pi n_{0s}{q}_s^2/m_s$. The linearization of the cold Lorentz force equations givesSolving this linear system and using the species current density ${\mathbf{j}}_s=n_{0s}{q}_s{\mathbf{u}}_s$ gives the conductivity tensorThe magnetized cold-plasma dielectric tensor is $\mathbf{ϵ}=\mathbf{I}+(4\pi i/\omega ){\sum }_s{\mathbf{\sigma }}_s$. ThusAt the stated electron frequency scale, the ion contribution is smaller by $O(Z{m}_e/m_i)$ for $q_i=Ze$ and quasineutrality. To this accuracy,These expressions require negligible collisions and the cold ordering $k{v}_{\mathrm{th},s}\ll |\omega |$; sufficiently close to a cyclotron resonance, thermal broadening or collisions must also be retained.

Write the ordinary plasma wave asThe eikonal function $S$ is dimensionless; it varies by order ${\epsilon }^{-1}$ across a background scale $L_{\mathrm{bg}}$, where $\epsilon \sim \lambda /L_{\mathrm{bg}}\ll 1$. The amplitudes and equilibrium plasma parameters vary on $L_{\mathrm{bg}}$, so $|\nabla {\mathbf{a}}_j|\ll k|{\mathbf{a}}_j|$ and $|\nabla k|\ll k^2$. The equilibrium is stationary, so $\omega$ is constant along propagation. Assume infinitesimal fields, negligible collisions, $k{v}_{\mathrm{th},s}\ll \omega$, negligible ion response, and a fixed ordinary plasma wave polarization with ${\mathbf{B}}_0\parallel \hat{\mathbf{z}}$ and spatial variation along $x$. The WKB approximation is used away from wave turning points and plasma wave resonances.

Replacing derivatives on the fast factor by $i\mathbf{k}$ and $-i\omega$ in the Maxwell equations reproduces the local magnetized cold-plasma dielectric tensor equation of b. In this one-dimensional ordinary plasma wave geometry, the Maxwell equations and the local cold parallel current density in fact give the scalar equation directly:For $\delta E_z=a(x)e^{iS(x)}$, substitution givesThe leading eikonal equation is $(S^{\prime }{)}^2=k^2(x)$, and the next order gives $2k{a}^{\prime }+k^{\prime }a=0$, hence the flux-normalized WKB amplitude $a\propto k^{-1/2}$. In particular, require $|k^{\prime }|/|k{|}^2\ll 1$ and $|k^{\prime \prime }|/|k{|}^3\ll 1$; these fail at the wave turning points.

The long-wavelength limit in d is resonant heating. Most electrons screen the slowly varying spatial field. The noise nevertheless excites collective Langmuir waves, whose sharp temporal spectrum lies near $\omega =\pm {\omega }_{pe}$. Electrons with $v\simeq \omega /k$ remain in phase with a wave and exchange energy with it. A stable tail absorbs net wave energy through Landau damping. The heating is initially deposited in narrow tail regions of velocity space, even though the volume-integrated rate equals that in c.

The short-wavelength limit in c is stochastic heating in the terminology intended here. With $ϵ\simeq 1$, the imposed electric field is essentially unscreened. Its temporally uncorrelated kicks cause a random walk of electron velocity, so $⟨(\Delta v{)}^2⟩=2D_0\Delta t$. Electrons throughout the thermal population contribute to the resulting kinetic energy increase; no narrow collective wave peak selects a small tail population.

Both mechanisms are driven by white noise, and the formal condition $\omega =kv$ enters the stationary quasilinear diffusion calculation in both cases. The physical distinction is the sharp collective resonance in d versus the broad direct random forcing in c, rather than the mere presence of randomness. In collisionless theory this heating is an increase of the averaged kinetic energy; subsequent collisions can establish thermal equilibrium.

For the strong-collision ion-neutral Alfvén modes, put $a=k_{\parallel }{v}_A$ and assume $0<|a|/\mu \ll 1$. Rewrite the cubic asThe two slow roots obey $\omega =O(a)$, so the first approximation is ${\omega }_0^2=a^2/2$. Set $\omega ={\omega }_0+\delta \omega$ and retain the next order, of size $a^3$:ConsequentlyThese are weakly damped Alfvén waves. Their phase speed is $v_A/\sqrt{2}=B_0/\sqrt{4\pi ({\rho }_i+{\rho }_n)}$, since collisions make the neutral mass participate in the oscillation. Their ambipolar damping rate is ${\Gamma }_A=a^2/(8\mu )$, with ${\Gamma }_A/|\mathrm{Re}{\omega }_{\pm }|\ll 1$.

For the remaining root, at $a=0$ the cubic is ${\omega }^2(\omega +2i\mu )$, so start from ${\omega }_{d0}=-2i\mu$. Denote the cubic by $F(\omega ,a)$. Its value and derivative at this root areThe first correction is therefore $\delta {\omega }_d=i{a}^2/(4\mu )$, givingThis branch is purely damped. Indeed, the real growth-rate polynomial for $s=-i\omega$ is $s^3+2\mu s^2+a^{2}s+\mu a^2$, and its simple real root near $-2\mu$ stays real for sufficiently small real $a^2$. It represents rapid decay of ion-neutral relative motion. In the $a\to 0$ limit, the other two frequencies tend to zero; all three branches are accounted for.

The neutral displacement equation gives, for each nonzero-frequency root,The vector displacements are parallel as complex polarization vectors, with this scalar amplitude and phase ratio.

For the rapidly damped root, $\mu -i{\omega }_d=-\mu +a^2/(4\mu )+O(a^4/{\mu }^3)$, and henceThe two fluids move almost equally and oppositely; collisions remove the counterflow at a rate close to $2\mu$.

For either weakly damped Alfvén wave, expand the ratio around unity:Equivalently, the exact slippage relation is ${\mathbf{\xi }}_i-{\mathbf{\xi }}_n=-(i\omega /\mu ){\mathbf{\xi }}_n$. Thus the relative slippage has size $|a|/(\sqrt{2}\mu )\ll 1$ and is primarily a small phase difference. The fluids almost share a common displacement, but magnetic tension acts directly on the ions and drag accelerates the neutrals. This imperfect locking causes ambipolar damping.

The sign of the damping is also checked by the perturbation energy balance. With periodic or vanishing-flux boundaries,It follows by dotting each linear momentum equation with its velocity and the induction equation with $\delta \mathbf{B}/(4\pi )$; magnetic-tension work cancels magnetic-energy change, pressure contributes only boundary flux, and the two drag terms combine into the negative square. Even though spatial diffusion was neglected, ion-neutral friction still dissipates the wave energy.

At each fixed semimajor axis, let $B=b_{3/2}^2(\alpha )/b_{3/2}^1(\alpha )$ and $Z_f=B{e}_{\mathrm{pl}}$. The planet's longitude of pericentre is the zero of longitude. The secular perturbation equation becomes $\dot{z}=iA(z-Z_f)$. Integrating it with the initial condition $z(0)=0$ givesThe forced eccentricity vector is the constant real $Z_f$; the free eccentricity vector is $-Z_f{e}^{iAt}$. Consequently the complex eccentricity follows a circle in the Argand plane centered at $Z_f$, with radius $Z_f$, starting at the origin and moving counterclockwise. Its initial velocity is downward, $\dot{z}(0)=-iA{Z}_f$. The secular period is $2\pi /A$, andThe longitude of pericentre is undefined at the origin, but the complex eccentricity passes smoothly through it. At $At=\pi$, $z=2Z_f$ and the orbit is maximally eccentric and aligned with the planet.

Using the small-$\alpha$ Laplace coefficient expansion givesAlso $n=n_{\mathrm{pl}}{\alpha }^{-3/2}$ to leading order in the planet-to-star mass ratio. The secular forcing by a distant outer planet is thereforeA relative $O({\alpha }^2)$ correction to $A$ can accumulate a phase error at very long times; the leading-frequency formula is not a uniform-in-time expansion.

If $a_2=a_1(1+\delta )$ with $\delta =\delta a/a_1\ll 1$, then $Z_{f,2}/Z_{f,1}=1+\delta +O({\delta }^2)$ and $A_2/A_1=1+3\delta /2+O({\delta }^2)$. The outer planetesimal therefore has a slightly larger forced eccentricity circle and a shorter secular period. Its accumulated secular phase lead is $\Delta (At)\simeq 3A_{1}t\delta /2$. Initially both circles start at the origin, but their eccentricity vectors progressively lose alignment. Linearizing the exponential in $\delta$ additionally requires $|A_{1}t\delta |\ll 1$.

The orbit crossing is caused by this accumulating differential secular phase, even though individual orbital eccentricities remain small. To first order in orbital eccentricity, the radius of an orbit at fixed inertial longitude $\theta$ isFor an infinitesimally neighboring orbit,This derivative is evaluated within the first-order orbital-shape model; derivatives of the neglected quadratic terms give relative $O(Z_f)$ corrections near crossing. The least separation over longitude is consequently $\delta a[1-|z+a{\partial }_{a}z|]$. The differential secular orbit-crossing criterion is $|z+a{\partial }_{a}z|=1$. Since $Z_f\propto a$ and $A\propto a^{3/2}$,The first term stays bounded by $4Z_f\ll 1$, whereas the second grows in proportion to time. To leading order in $Z_f$, crossing begins when $(3/2)Z_{f}At=1$. ThusMore precisely, in this first-order orbital-shape approximation the exact first crossing lies within a fractional $O(Z_f)$ window of this value, since the magnitude differs from $(3/2)Z_{f}At$ by at most $4Z_f$. Higher-order Laplace coefficients give additional corrections. The local derivation takes the neighboring-orbit limit before the long-time limit; for a finite separation it requires $\delta a/a_1\ll Z_f$ near crossing. A broad continuous disk contains such neighboring orbits. The crossings first appear in directions selected by the phase of $z+a{\partial }_{a}z$, not simultaneously at every longitude.

For the collision energy estimate adopt the phase-averaged mean orbital eccentricity $\overline{e}=4Z_f/\pi \simeq 5\alpha e_{\mathrm{pl}}/\pi$. The specified velocity estimate givesUsing the projectile kinetic energy per target mass convention for specific impact energy, equal-size and equal-density bodies have equal masses and $Q=\frac{1}{2}{v}_{\mathrm{col}}^2$. A catastrophic planetesimal collision requires $Q\gtrsim Q_D^{\ast }$, soThe numerical coefficient is only illustrative because the collision speed was specified only to order of magnitude. If the disruption threshold is instead defined using the centre of mass kinetic energy per combined mass, equal masses give $Q_R=v_{\mathrm{col}}^2/8$ and the coefficient becomes $8{\pi }^2/25$ for a threshold expressed in that convention. The robust result is a lower bound of order $Q_D^{\ast }{a}_{\mathrm{pl}}/(G{M}_{\star }{e}_{\mathrm{pl}}^2)$.

Within the inner-disk approximation, catastrophic collisions require sufficiently large $a_1/a_{\mathrm{pl}}$, not arbitrarily small radii. Although the Keplerian speed rises inward, the secularly induced orbital eccentricity falls faster: $v_{\mathrm{col}}^2\propto \alpha$. The lower bound must be much smaller than one to leave a domain compatible with $\alpha \ll 1$. A small planet orbital eccentricity, a large disruption threshold, or a distant planet can eliminate that domain. The planet mass controls the time to crossing but cancels from the eventual forced eccentricity amplitude and this collision-energy criterion. The threshold alone does not ensure that crossing has occurred within the disk's available lifetime.

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.

Use nonrelativistic mechanics and let $\mathbf{V}=(m_1{\mathbf{v}}_1+m_2{\mathbf{v}}_2)/(m_1+m_2)$ be the constant centre-of-mass velocity. In the centre-of-mass frame, conservation of momentum gives $m_1{\mathbf{u}}_1+m_2{\mathbf{u}}_2=0$ before and after the elastic collision. Thus ${\mathbf{u}}_2=-(m_1/m_2){\mathbf{u}}_1$, and the kinetic energy isThe decomposition $T_{\mathrm{lab}}=T_{\ast }+\frac{1}{2}(m_1+m_2)V^2$ and constant centre-of-mass velocity show that the elastic collision preserves $T_{\ast }$, so the positive constant multiplying $u_1^2$ forces that speed to remain unchanged. The same momentum relation fixes the second speed:Only their directions can change; the two velocities remain antiparallel in this frame. This derivation does not require a head-on collision.