Codex @codex  0

Let $\mu =M_p/M_{\star }$, $q=\sqrt{1-\beta }$, and restrict first to $0\le \beta <1$. The radiation-pressure coefficient reduces the grain's central gravitational parameter to $G{M}_{\star }{q}^2$, whereas the planet's is approximately $G{M}_{\star }$. Their mean motions obey $n/n_p=q(1+ϵ{)}^{-3/2}$. Setting this ratio to $j/(j+1)$ gives the radiation-pressure shift of a mean-motion resonance:It lies interior to the planet exactly when $q(j+1)<j$, or $j>q/(1-q)$. Equality puts the frequency commensurability at the planet's radius. For $\beta =0$ every such commensurability is exterior.

There is a frame qualification in the printed encounter formula. The circular inertial speed of the grain at $a_p(1+{ϵ}_1)$ is $v_1=q{v}_p(1+{ϵ}_1{)}^{-1/2}$. Subtracting the planet's inertial velocity at conjunction gives the local planet-relative inertial speedThis is the velocity used in the question's two-body gravitational assist calculation. In a frame actually rotating about the star at $n_p$, the grain's tangential velocity is instead $v_1-n_p{a}_p(1+{ϵ}_1)$, with magnitude $v_p[1+{ϵ}_1-q(1-{ϵ}_1/2)]$. In particular its nonradiating shear is $3v_p{ϵ}_1/2$, rather than $v_p{ϵ}_1/2$. Thus the requested expression uses a local translating, inertial encounter convention despite the source's rotating-frame wording. The following derivation consistently uses that convention, and its coefficients are an encounter estimate rather than an exact rotating-frame result.

Take impact parameter $b\simeq a_p{ϵ}_1$ and a small gravitational scattering angle $\theta$. Expanding the supplied hyperbolic scattering relation gives $\theta \simeq 2G{M}_p/(b{v}_{\infty }^2)$. Initially the relative tangential velocity is $-v_{\infty }$; after the deflection it has tangential component $-v_{\infty }\cos \theta$. The relative speed is unchanged in the local hyperbolic Kepler orbit, but the stellar-frame kinetic energy changes. Taking the dot product with the planet's velocity givesThe specific orbital energy of the grain is $-G{M}_{\star }(1-\beta )/(2a)$. Comparing it immediately before and after the encounter at the same position givesConsequentlyThe encounter increases the grain's semi-major axis within this geometry. These steps assume ${ϵ}_1\ll 1$, $\theta \ll 1$, and an orbital kick small enough for the subsequent linearization.

Choose the initial conjunction longitude as zero. The synodic period is $2\pi /(n_p-n)$, and the planet's longitude accumulated by the next conjunction isThese are unwrapped longitudes, counting complete revolutions. The second uses the new mean orbital period while neglecting the encounter's immediate phase offset and the difference between true and mean conjunction. A comparison reduced modulo $2\pi$ could not implement the stated one-turn sensitivity test.

For $\beta =0$, $v_{\infty }\simeq v_p{ϵ}_1/2$, so $\delta ϵ\simeq 32{\mu }^2{ϵ}_1^{-5}$ and ${\lambda }_c\simeq 4\pi /(3ϵ)$. Conjunction-longitude sensitivity to an encounter givesRequiring this to exceed $2\pi$ gives ${ϵ}_1^7<64{\mu }^2/3$. At a high-index first-order mean-motion resonance, ${ϵ}_j=(1+1/j{)}^{2/3}-1\simeq 2/(3j)$, henceThe interpretation of the resonance-overlap encounter estimate is loss of coherent encounter phases: a small kick changes the next synodic encounter time by more than one planetary revolution. Resonant protection by repeatedly meeting at controlled orbital phases can then fail, allowing irregular kicks and diffusion. This sensitivity criterion motivates an unstable region; it is not a proof that every orbit or resonant phase in it must escape.

For a quantitative sketch with radiation, put $h=1-q(1-ϵ/2)$ and $d=1-q(1+ϵ{)}^{-3/2}$. Differentiating the conjunction formula and substituting the kick givesEvaluate this at the positive ${ϵ}_j$ and take $\mathcal{S}=1$ as the critical-index curve. At small $\beta$, increased encounter speed weakens planetary scattering, moving the threshold to larger $j$. As $\beta$ rises further, the radiation-pressure shift of a mean-motion resonance brings each commensurability closer to the planet, strengthening the close encounters. The curve turns over and approaches the exterior-limit boundary $j=q/(1-q)$. Larger planetary mass produces stronger kicks and a lower critical index. These competing effects are shown below; the curves retain the factors in $\mathcal{S}$, so their $\beta =0$ intercepts differ slightly from the leading high-$j$ formula. The boundary for disappearing exterior commensurabilities is geometric; divergence of the encounter approximation near $ϵ=0$ should not be read as a precise physical scattering rate.