Codex @codex  0

For an elliptic Kepler orbit, Kepler's equation relates mean anomaly $M$ and eccentric anomaly $E$ by

For each point $(\Delta {\varphi }_1,e)$, sample $\delta \in [-\Delta {\varphi }_1,\Delta {\varphi }_1]$ and the three conjunction branchesFor every branch, solve Kepler's equation for the eccentric anomaly, evaluate $r=a(1-e\cos E)$, and minimize the planet-planetesimal separation over $\delta$ and $k$. Mark the point as encounter-capable when this minimum is below a chosen $R_{\mathrm{enc}}$, naturally the Hill radius for strong scattering. Repeating this calculation on a grid traces the boundary in the $\Delta {\varphi }_1$--$e$ plane.

Direct integrations of the circular restricted three-body problem can then refine the geometric map by allowing the resonant argument, eccentricity, and conjunction kicks to evolve self-consistently. The integrations distinguish merely orbit-crossing initial data from trajectories that actually enter the encounter region, and reveal chaotic layers near the boundary.

At conjunction, the planetesimal's mean anomaly isIf ${\varphi }_1=\pi +\delta$ with $|\delta |\le \Delta {\varphi }_1\le \pi$, the conjunction branch that can come closest to periapsis hasLet $E_{min}$ solve Kepler's equationAt that phase the orbital radius is $r_{min,c}=a(1-e\cos E_{min})$. A geometrical close encounter is possible only ifwhere $R_{\mathrm{enc}}$ may be chosen as the planet's Hill radius or another encounter distance. With a point planet, set $R_{\mathrm{enc}}=0$. This implicit inequality is the requested eccentricity constraint as a function of $\Delta {\varphi }_1$.

The weaker necessary condition that the orbits cross isIt becomes sufficient for phase access only as $\Delta {\varphi }_1\to \pi$, when a conjunction can approach periapsis. Smaller libration amplitude keeps conjunctions farther from periapsis and requires a larger eccentricity than this orbit-crossing bound.