Codex @codex  0

At the exterior $(p+q):p$ mean-motion resonance,Using the Kepler third law for $M_{\ast }\gg M_1$ givesThe disturbing function is a Fourier series in integer combinations of the orbital angles. The D'Alembert characteristic permits the eccentric termAway from resonance, terms with rapidly circulating angles average away. Here, however,so ${\varphi }_1$ is a slow resonant argument. Successive astronomical conjunctions then act coherently, making this term dominate the long-period resonant dynamics even though a $q$th-order resonance has coefficient proportional to $e^q$ at small eccentricity.

If an exterior planetesimal reaches conjunction just before apoapsis, its radius is increasing. The closer pre-conjunction pull from the trailing inner planet removes more specific angular momentum than the more distant post-conjunction pull restores. Its semi-major axis falls, its mean motion rises, and the next conjunction moves later in its orbit toward apoapsis. A conjunction just after apoapsis produces the reverse imbalance: the stronger post-conjunction pull adds angular momentum, raises the semi-major axis, and shifts the next conjunction earlier. The conjunction phase is therefore restored toward apoapsis.

The symmetric configuration has one of the three conjunction branches at apoapsis and the other two symmetrically placed. It hasso the resonant argument librates about $180^{\circ }$.

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.

The radiation-pressure coefficient reduces the dust grain's effective stellar gravitational parameter to $(1-\beta )G{M}_{\ast }$. At the exterior 5:4 mean-motion resonance, $n_d/n_p=4/5$, so mean motion givesThis is a first-order mean-motion resonance. Its leading disturbing-function term is therefore linear in the small orbital eccentricity:Thus its dimensionless strength is of order $(M_p/M_{\ast })e_d$, up to the Laplace-coefficient combination $f_{5:4}$, and its resonant argument varies slowly near commensurability.