Codex @codex  0

In a planetary gravitational assist, the incoming and outgoing relative speeds have the same magnitude. If $\mathbf{u}$ is the planet-frame velocity, the heliocentric specific orbital energy change is consequentlynot an arbitrary addition of $|\Delta \mathbf{u}|$ to the heliocentric speed. The initial binding energy per unit mass is $G{M}_{\star }/(Q+q)\simeq G{M}_{\star }/Q$.

For a weak deflection, $b\gg b_{90}=G{M}_p/v_{\infty }^2$, the hyperbolic scattering formula gives $|\Delta \mathbf{u}|\simeq 2G{M}_p/(b{v}_{\infty })$. Its direction is perpendicular to the incoming relative velocity. At the planet, in the near-parabolic prograde limit,On the favorable side of the planet the energy gain is approximately $2G{M}_p{v}_p\kappa /(b{v}_{\infty })$. Equating this to the binding energy gives the ejection impact parameterOnly one sign of impact parameter gains energy for each radial branch. Noncolliding ejections have $R_p<|b|<b_{\mathrm{ej}}$, giving width $(b_{\mathrm{ej}}-R_p{)}_{+}$, where $x_{+}=\max (x,0)$. For $b_{\mathrm{ej}}\gg R_p$, the planetary ejection of a comet rate is approximately ${\Gamma }_{\mathrm{ej}}={\Sigma }_c{v}_{\infty }{b}_{\mathrm{ej}}$. With ${\mu }_p=M_p/M_{\star }$ this givesThe scaling applies when $b_{90}\ll b_{\mathrm{ej}}\ll r_H$ so that the local two-body encounter and weak-deflection approximations are consistent. The favorable-side coefficient assumes a coplanar orbit; averaging over other encounter geometries changes it. Cumulative diffusion by many subthreshold encounters is not included in this one-encounter estimate.