Codex @codex  0

Write $h=r^2\dot{\theta }$ for the conserved specific angular momentum. Zero specific orbital energy in the inverse-square Newtonian gravitational field givesDefine $r_0=h^2/(2GM)$, the periapsis radius for a nonradial parabolic trajectory. ThenOn the outgoing branch, integrate from the periapsis time $t_0$ and put $s=\sqrt{r-r_0}$:The incoming branch has the opposite sign. Squaring gives the radius-time relation for a parabolic Kepler orbit, a version of the Barker equation:Here $r_0$ and $t_0$ are fixed by angular momentum and the time origin at periapsis. The relation applies on both branches, with $r\ge r_0$; it is not necessary to solve for the polar angle.

For a fixed central gravitational source, the radial and transverse components of Newton's second law in polar coordinates areThe second equation gives $r^2\dot{\theta }=h$, the constant specific angular momentum. Take $h>0$ by choosing the orientation of $\theta$, and put $u=1/r$. Writing $u^{\prime }$ and $u^{\prime \prime }$ for derivatives with respect to $\theta$,Substituting into the radial equation and dividing by $-h^2{u}^2$ proves the Binet equationAt perihelion the velocity is tangential. The escape velocity is $v_0=\sqrt{2GM/r_0}$, so $h=r_0{v}_0$ and $h^2=2GM{r}_0$. Consequently $k=1/(2r_0)$. The initial conditions are $u(0)=1/r_0$ and $u^{\prime }(0)=0$. Solving the linear ordinary differential equation gives the parabolic Kepler orbitUsing angular momentum conservation again,On the outgoing branch $r=2r_0$ is first reached at $\theta =\pi /2$. Integrate from release at perihelion, using $w=\tan (\theta /2)$:More generally the same integral gives $t=(2r_0^2/h)(w+w^3/3)$, the Barker equation for this parabolic orbit.

Put $k=p+1$, $\mu =G{M}_{\star }$ and $n_{\mathrm{pl}}=\sqrt{\mu /a_{\mathrm{pl}}^3}$. To leading order the mean-motion resonance fixes $a=a_{\mathrm{pl}}{k}^{-2/3}$. The specific angular momentum isThe asteroid's instantaneous angular speed is $h/r^2$. Thus the high-eccentricity angular-speed crossover satisfies $h/r_x^2=n_{\mathrm{pl}}$, givingHere $q\ll r_x\ll a$, so a parabolic Kepler orbit gives the local motion accurately. Its polar equation is $r=2q/(1+\cos f)$. Write $f_x=\pi -{ϵ}_x$. Since $1+\cos f_x\simeq {ϵ}_x^2/2$,This is the outward crossing; the inward crossing has true anomaly $-f_x$ modulo $2\pi$.

For the pericentre-to-crossover flight time use $D=\tan (f/2)$. The parabolic Kepler orbit has $r=q(1+D^2)$ and $dt=r^{2}df/h$, whence direct integration gives the Barker equationAt the crossover $D_x\simeq \sqrt{r_x/q}\gg 1$. Consequently the planet's angular displacement isThe error tends to zero as $q/a_{\mathrm{pl}}\to 0$ at fixed $k$; it includes the finite binding energy of the elliptic Kepler orbit as well as the subleading term in the Barker equation.

At the $j$th pericentre passage the asteroid's unwrapped mean longitude is $\lambda =\varpi +2\pi j$. The resonant argument therefore givesThus $\varphi /k$ specifies the planet's direction relative to the asteroid's longitude of pericentre when the asteroid passes pericentre, with $k$ branches separated by $2\pi /k$. It is not the instantaneous true-longitude separation throughout the orbit. At high orbital eccentricity the asteroid quickly sweeps to a direction nearly opposite pericentre, while the planet scarcely moves. During the long outer excursion the asteroid's direction changes slowly and the planet advances substantially. This makes an astronomical conjunction during the outer excursion much more dangerous than an astronomical conjunction close to pericentre.

At the following apocentre, the planet has advanced $\pi /k$. The asteroid's direction in the rotating reference frame is consequentlyThese $k$ apocentre directions are equally spaced. Phase protection of interior integer resonances places the planet in the largest gap between them, rather than at an apocentre direction. For $k=2$, $\varphi =\pi$ puts an apocentre toward the planet, whereas $\varphi =0$ puts the apocentres at $\pm \pi /2$. For $k=3$, $\varphi =0$ puts an apocentre toward the planet, whereas $\varphi =\pi$ puts the apocentres at $\pi /3,\pi ,5\pi /3$. Hence the geometrically favored resonant-argument libration centres in this regime of high orbital eccentricity areThis is an encounter-avoidance argument for likely stability, not a calculation of a resonant Hamiltonian or a claim that every orbit near either centre is stable at arbitrary orbital eccentricity.

The following rotating-frame drawing uses $a_{\mathrm{pl}}=1$, $a=3^{-2/3}$, $q=0.002$ and $\varphi =\pi$. It follows three successive Kepler orbits, which fill one planet orbital period and close the pattern. The blue segment is just the first outward half-orbit. The marked crossover is a maximum of its rotating-frame polar angle, since $d\psi /dt=h/r^2-n_{\mathrm{pl}}$ changes from positive to negative there. The planet stays at $(1,0)$; the three apocentres avoid that direction.

To estimate angular-momentum kicks to nearly radial orbits, note that $q\simeq h^2/(2\mu )$, so an order-one fractional change in $q$ requires $|\Delta h|$ of order $h$. A characteristic distant encounter at radii of order $a_{\mathrm{pl}}$ has tangential acceleration of order $G{M}_{\mathrm{pl}}/a_{\mathrm{pl}}^2$ and duration of order $n_{\mathrm{pl}}^{-1}$. Its torque per unit asteroid mass therefore producesWith the most favorable coherent signs, the crude encounter count isOrder-one factors depend on $p$ and encounter geometry. Uncorrelated signs instead give a random walk count of order $(M_{\star }/M_{\mathrm{pl}}{)}^{2}q/a_{\mathrm{pl}}$ when this is large. A librating phase-protected orbit can suppress the kicks even further.

There is an important limit to calling this a minimum. The question does not specify an encounter impact parameter. For a weak close encounter with impact parameter $b$ and relative speed $u\sim \sqrt{G{M}_{\star }/a_{\mathrm{pl}}}$, the impulse estimate gives $|\delta h|\sim G{M}_{\mathrm{pl}}{a}_{\mathrm{pl}}/(bu)$ and hence $N\sim (M_{\star }/M_{\mathrm{pl}})(b/a_{\mathrm{pl}})\sqrt{q/a_{\mathrm{pl}}}$, capped below by one. In the 2:1 case the Kepler orbit can cross the planet's orbit, so one suitably close encounter can suffice: there is no impact-parameter-independent large lower bound. For $k\ge 3$, $Q\simeq 2a_{\mathrm{pl}}{k}^{-2/3}<a_{\mathrm{pl}}$ bounds the separation below by $a_{\mathrm{pl}}-Q$, and supplies a geometric factor at fixed $k$. The boxed count is the usual orbital-scale, maximally coherent estimate; the assumptions are essential.

Treat the just-bound orbit as a parabolic Kepler orbit with $q\simeq a_p$. Its specific angular momentum is $h=\sqrt{2G{M}_{\ast }q}$, so its instantaneous angular speed is $h/r^2$. Equating this to the planet's mean motion $n_p=\sqrt{G{M}_{\ast }/a_p^3}$ givesThe parabolic orbit equation $r=2q/(1+\cos f)$ then gives $f\simeq 47^{\circ }$, or $50^{\circ }$ to the nearest ten degrees.

If the bodies were in conjunction at periapsis, Barker equation givesThe planetesimal is therefore at its maximum lead, approximatelyahead of the planet as seen from the star.

On a nonradial zero-energy Kepler orbit, conserved specific angular momentum $h$ gives periapsis $r_p=h^2/(2GM)$ and ${\dot{r}}^2=2GM(r-r_p)/r^2$. Integrating with $s=\sqrt{r-r_p}$ yields $t-t_p=\pm \sqrt{2/(GM)}(r+2r_p)\sqrt{r-r_p}/3$. Squaring gives the displayed version of the Barker equation, valid on both incoming and outgoing branches. The time origin $t_p$ is passage through periapsis, not an arbitrary point on the orbit.

“Arrives at infinity with zero speed” is an asymptotic statement. For outward radial motion at zero specific orbital energy,Every finite distance takes finite time, but $t(r)\to \infty$ as $r\to \infty$. There is no instant of rest at an infinitely distant destination. Differentiating the displayed solution gives $\ddot{r}=-\mu /r^2$, independently checking Newton's second law. In the synthetic example, reaching $r=4\mathrm{m}$ takes exactly $(4-\sqrt{2})/3\mathrm{s}\approx 0.861929\mathrm{s}$. For nonzero specific angular momentum, the Barker equation gives the corresponding infinite-time limit along a parabolic trajectory. Positive-energy escape also takes infinite time to reach infinity; it retains nonzero asymptotic speed.