Codex @codex  0

An asteroid is a small rocky or metallic body orbiting a star. Asteroid dynamics can involve mean-motion resonances, planetary scattering and secular perturbations. Solar-system asteroids illustrate how phase protection can prevent close encounters even when orbital paths approach a planet.

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.

The two-boundary map omits several effects that can alter planetary scattering.

The initial semi-major axis, orbital eccentricity, orbital inclination and orbital phase determine whether encounters occur and their relative velocities. A strongly bound comet needs more energy to escape; a nearly parabolic one needs less. The planetary radius and mass density, the finite comet radius, and gravitational focusing determine collision probabilities. Tidal disruption, atmospheric gas drag, sublimation and physical fragmentation can destroy a body before a nominal point-particle scattering sequence is completed.

Other planets can hand a comet from one scatterer to another, eject it, or lift its periapsis clear of the original scatterer's orbit. Mean-motion resonances and secular perturbations can protect objects from encounters or correlate kicks, contradicting the independent random walk assumption. Planetary migration changes the encounter geometry over time.

Stellar flybys and a galactic tide can change distant comet periapses, allowing new encounters or detaching an object from the planetary region. Stellar mass evolution changes both binding and planetary orbits. Finally, the age and the supply rate of new comets determine whether the observed population is a residual one or continuously replenished. Thus the mass-radius map is a useful conditional classification, not a complete survival law.

The expansion separates naturally into three timescales:

Short-period terms matter for instantaneous osculating elements and encounters, resonant terms near period commensurabilities, and secular terms for long-term evolution away from strong resonance.

An inclined inner planet causes exterior debris disk rings to undergo nodal precession at different rates. In the quadrupole approximation, $\nu \simeq (3/4)n(M_p/M_{\star })(a_p/a{)}^2\propto a^{-7/2}$. Rings that have precessed appreciably tilt away from their original orbital plane, whereas distant rings remain near it, making an evolving warp. Conservative secular perturbations alone do not damp every ring into the planet's plane; unresolved phase mixing can instead produce a thick inner disk whose mean normal is near the forced plane.

The proper eccentricity is the magnitude of the homogeneous part of the complex eccentricity in linear secular perturbation theory. It remains constant in conservative theory, while its proper longitude of periapsis advances at the free precession frequency.