Planetary scattering

Planetary scattering is the alteration of a smaller body's orbital energy and angular momentum by close encounters with a planet. Repeated encounters can eject, accrete, or detach the body.

For a nearly radial Kepler orbit, $q\simeq h^2/(2GM)$, so modest changes in specific angular momentum can strongly change pericentre. A planetary torque acting for a reference orbital time gives $\delta h\sim (m_p/M)\sqrt{GM{a}_p}$. Coherent kicks change $q$ by order itself after $N\sim (M/m_p)\sqrt{q/a_p}$ orbital-scale encounters, whereas uncorrelated kicks require an order-$N^2$ random-walk count. Both estimates are conditional on encounter geometry: a close encounter with a smaller impact parameter can give a much larger kick.

Model successive uncorrelated planetary kicks as a random walk in inverse semi-major axis $x=1/a$, proportional to negative specific orbital energy. If the root mean square kick per periapsis passage is $s_x$ and the initial value is $x_0$, the characteristic number of passages to reach the escape scale is $N\sim (x_0/s_x{)}^2$. For $s_x=10M_p/(a_p{M}_{\star })$ and passage interval $P(a_0)$,The quantity $s_x$ has inverse-length units, not the units of a variance-per-time diffusion coefficient. With $⟨(\Delta x{)}^2⟩=2D_{x}t$, that coefficient is $D_x=s_x^2/(2P)$. This estimate assumes many independent encounters; changing orbital periods and resonant correlations can invalidate the clock. It is not an exact mean first-passage time.

At fixed stellar mass $M_{\star }$, planetary mass density ${\rho }_p$, age $t_{\star }$ and initial ratio $\kappa =a_0/a_p$, the order-unity Safronov number comparison and comet energy diffusion give boundariesTheir logarithmic slopes in planetary mass versus radius are $-3/2$ and $3/4$. Above both, strong kicks and a short estimated diffusion time favour ejection; below the escape boundary, repeated weak encounters favour collision statistically. Below the age boundary, the stated diffusion mechanism has not completed on its characteristic clock. Collision lifetimes require additional cross-section and encounter-rate information, so this map alone cannot establish retention or accretion within the age.

For a planet on a circular stellar orbit, a common convention is $\Theta =\frac{1}{2}(v_{\mathrm{esc}}/v_k{)}^2$. The escape-to-orbital-speed ratio measures the size of a strong gravitational assist relative to the stellar binding energy of a nearby comet. Large values tend to favour ejection over collision; small values favour repeated weak encounters with more collision opportunities. This is a statistical comparison, dependent on the initial orbit and encounter geometry, rather than a universal outcome threshold.

A gravitational assist ejects a weakly bound comet when the increase in its specific orbital energy exceeds its binding energy. For near-parabolic coplanar encounters, the weak-deflection ejection impact parameter scales as $(M_p/M_{\star })Q$. Combining this with the phase-mixed comet surface density gives a one-encounter ejection time proportional to $Q^{1/2}$, within the local two-body encounter approximation.

For a particle initially having semimajor axis one and encountering a circular planet at radius one, the incoming stellar-frame speed is one. A conservative encounter at planet-relative speed $u$ changes the stellar-frame velocity by at most $2u$. Since stellar escape requires speed $\sqrt{2}$, a necessary bound is $u>(\sqrt{2}-1)/2$. Retaining the outgoing relative-velocity geometry sharpens the necessary bound to $u>\sqrt{2}-1$. Neither inequality by itself guarantees a realizable deflection at a specified impact parameter.

The geometric collision cross-section in a coplanar calculation is a length $2R_p$, rather than an area. For one comet with local surface probability ${\Sigma }_c$ and relative velocity $v_{\mathrm{rel}}$, the rare-event rate is $2R_p{\Sigma }_c{v}_{\mathrm{rel}}$. Without gravitational focusing, the resulting mean collision time scales as $Q^{3/2}{M}_p^{-1/3}{a}_p$ for deep planet-crossing orbits at fixed stellar mass and planetary mass density.

For a population with randomized vertical phases, orbital inclination $I$, and $R_p\ll a_{p}I$, the effective volume probability near the planetary plane is proportional to ${\Sigma }_c/(a_{p}I)$. The collision time then scales as $I{Q}^{3/2}{M}_p^{-2/3}{a}_p^2$, at fixed stellar mass and planetary mass density. An individual inclined Kepler orbit whose orbital nodes do not intersect the planet's orbit need not collide at all.