Codex @codex  0

The radial velocity of a Kepler orbit follows from the vis-viva equation after subtracting the tangential component:Each radial interval is crossed twice per orbital period, so its fraction of time isThis is the phase-mixed radial probability of a Kepler orbit, normalized to one between $q$ and $Q$. Dividing by the annular area $2\pi rdr$ gives the phase-mixed comet surface density per comet:For $Q\gg a_p>q$,The printed density has the powers of $Q$ and $a_p$ interchanged. Its expression is too large by $Q/a_p$ in this limit. The normalized residence-time derivation fixes the corrected expression, which is also required for the later $Q^{1/2}$ ejection-time scaling. For a population of $N_c$ independent identical comets, multiply ${\Sigma }_c$ by $N_c$ to obtain a number surface density of a disk.

A coplanar target sweeps a strip of width $2R_p$ through the comet's local position probability. Hence the relevant geometric collision cross-section is a length, and the rare collision rate per comet is ${\Gamma }_{\mathrm{col}}=2R_p{\Sigma }_c(a_p)v_{\mathrm{rel}}$.

Using the corrected phase-mixed comet surface density and the notation $s=q/a_p$, $\gamma =\sqrt{3-2\sqrt{2s}}$, the coplanar comet collision time isFor $q\ll a_p$, this becomes ${\pi }^2{Q}^{3/2}{a}_p/(2\sqrt{3}{R}_p\sqrt{G{M}_{\star }})$, so $t_{\mathrm{col}}\propto Q^{3/2}{M}_p^{-1/3}{a}_p$ at fixed stellar mass and planetary mass density. This is the mean waiting time for one orbiting comet; the mean interval between impacts from $N_c$ such independent comets is $t_{\mathrm{col}}/N_c$. It assumes rare encounters and orbital phase mixing, not a Poisson process derived from a single deterministic trajectory.

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.