Codex @codex  0

Validity and the complete angular equation. The Lagrange planetary equations here describe a planar, circular-planet, small-orbital eccentricity, weakly perturbed exterior mean-motion resonance. Require $M_{\mathrm{pl}}/M_{\star }\ll 1$, $e\ll 1$, one isolated slowly varying resonant argument, perturbation times long compared with an orbital period, and no close planetary encounters. Other resonant harmonics and short-period terms must be negligible; coefficients may be frozen only across a narrow range of semi-major axis. Take $p>0$, $q=1$ or $2$ in a reduced integer ratio. Exactly $e=0$ makes the longitude of periapsis undefined, so the displayed angular variables must then be replaced by nonsingular eccentricity vector components.

Neglecting the perturbation to the accumulated mean longitude gives $\dot{\lambda }\simeq n(a)$; this is an accumulated-phase approximation, not differentiation of a fictitious expression $n(t)t$ while discarding $t\dot{n}$. With $n^2{a}^3=G{M}_{\star }$,Differentiating, using $\dot{n}=-3n\dot{a}/(2a)$ and the supplied orbital element equations, gives the full result within this constant-coefficient model:An equivalent formula replaces the last two terms in brackets by $q^2{C}_r{e}^{q-2}[(p+q)n-p{n}_{\mathrm{pl}}-2q{C}_s]-2q^3{C}_r^2{e}^{2q-4}\cos \varphi$. No derivatives of $C_r,C_s$ are included because they were specified as constants.

A full fixed point requires $\sin \varphi =0$ and $\dot{\varphi }=0$. Hence ${\varphi }_{\ast }=0$ or $\pi$, with its corresponding detuned semi-major axis fixed byAt either point $\dot{a}=\dot{e}=0$. The longitude of periapsis still precesses, so the fixed point refers to the reduced resonant dynamics.

Linear stability analysis and encounter geometry. Linearizing the full equation about a fixed point givesFor $q=2$, $\pi$ is a center equilibrium and zero is a saddle equilibrium. For $q=1$, $\pi$ is always a center equilibrium, while zero is a saddle equilibrium ifAt sufficiently small orbital eccentricity the full truncated equations instead admit a center at zero as well; equality is a degenerate case requiring higher-order analysis. Thus an unconditional instability claim at zero does not follow from the full equation. In the usual fixed-eccentricity weak-resonance regime the inequality holds and the stable libration center is $\pi$ for either order.

The physical explanation is resonance protection. At equal mean longitudes, the conjunction direction relative to periapsis obeys $q(\Lambda -\varpi )=\varphi$ modulo $2\pi$. For $q=1$, $\varphi =0$ places conjunction near periapsis, where the exterior particle comes closest to the planet; $\varphi =\pi$ places it near apoapsis. For $q=2$, zero includes both apsidal conjunctions, one of them near periapsis, whereas $\pi$ places the two conjunction branches near quadrature. The protected arrangements give restoring kicks in the positive-$C_r$ leading-harmonic model. This geometric explanation is conditional on its non-crossing, small-orbital eccentricity approximation; it does not override the low-$e$ term retained in the full stability calculation.

Small resonant-argument librations. Define the nominal resonant semi-major axis and the displaced center byThis follows by expanding $n=n_r[1-3(a-a_r)/(2a_r)]$ and taking $\cos (\pi +\delta )=-1+O({\delta }^2)$. It gives the requested initial semi-major axis to first order in the precession-induced detuning; $e_0=e(0)$ and $\delta =\Delta \varphi$.

For the pendulum approximation of a mean-motion resonance, also require fractional changes in $e$ small across a libration, and precession terms small compared with $n$. In particular $C_r/(n_r{e}_0^3)\ll 1$ for $q=1$ makes the extra full-equation curvature negligible; for $q=2$, small fractional eccentricity changes require $C_r/(n_r{e}_0^2)\ll 1$. Freeze $e=e_0$ in the leading restoring coefficient and setThen $u=\varphi -\pi$ obeys simple harmonic motion. Using $u(0)=\delta$, $\dot{u}(0)=0$ and integrating the remaining Lagrange planetary equations givesThe omitted apsidal modulation comes from higher perturbative orders and $u^2$. All expressions are leading resonant approximations, valid while the omitted fractional changes remain small. The small-libration phase portrait iswith period $2\pi /{\omega }_0$. The semi-major axis half-width is $2a_r{\omega }_0|\delta |/(3p{n}_{\mathrm{pl}})$ and the orbital eccentricity half-width is $q{C}_r{e}_0^{q-1}|\delta |/{\omega }_0$.

Dividing the original $\dot{e}$ and $\dot{a}$ equations gives the exact invariant of their retained terms,Thus the $e$–$a$ plot is a segment of this increasing curve traversed back and forth, not a closed ellipse. Locally its slope is $q/[2(p+q)e_0{a}_c]$.

Finite amplitude and the separatrix. Keep the same weak-resonance approximation but allow $\delta$ to be finite. The leading angular equation is $\ddot{\varphi }={\omega }_0^2\sin \varphi$, or $\ddot{u}=-{\omega }_0^2\sin u$. Multiplication by $\dot{u}$ givesThis resonant pendulum energy is constant in the pendulum approximation of a mean-motion resonance. It is not an exact integral of the earlier full equation when $e$ and its restoring coefficient vary: the $\dot{\varphi }\sin \varphi$ term alone prevents that conclusion in general.

For $|\delta |<\pi$ the phase portrait has closed libration curves around $\pi$, withwhere $K$ is the complete elliptic integral of the first kind. The largest speed is $2{\omega }_0\sin (|\delta |/2)$ and the leading semi-major axis half-width is $4a_r{\omega }_0\sin (|\delta |/2)/(3p{n}_{\mathrm{pl}})$. As $|\delta |\to \pi$, the energy approaches $2{\omega }_0^2$, the separatrix through the saddles at zero and $2\pi$; the particle spends increasingly long intervals near those saddles. With $ϵ=\pi -|\delta |$,At exactly $\delta =\pi$, $\dot{\varphi }=0$, the particle stays at the unstable fixed point in the ideal model. Nonstationary separatrix trajectories approach the saddle only in infinite time; above separatrix energy, the resonant argument circulates and resonance protection is lost. These distinctions matter when describing the limiting motion.

The transverse velocity is $v_{\theta }=h/r$. Subtracting its square from the vis-viva equation, and putting $u=a/r$, givesConsequentlyThe positive branch is outward motion and the negative branch is inward motion. The PDF's unsigned expression gives the radial speed, or the outward branch, rather than the signed radial velocity everywhere. For example, on the inward part of any noncircular elliptic orbit, $\dot{r}<0$.

The same result follows from $\dot{r}=(\mu e/h)\sin f$, with $f$ the true anomaly. The radicand is nonnegative precisely on $1/(1+e)\le u\le 1/(1-e)$; it vanishes at periapsis and apoapsis. For a circular Kepler orbit, $e=0$, $u=1$ and $\dot{r}=0$.

Use the separation vector $\mathbf{r}$ from the star to the comet, and let $m_{\mathrm{red}}=m{M}_{\star }/(m+M_{\star })$. The two-body problem separates into uniform centre of mass motion andTaking a cross product with $\mathbf{r}$ proves conservation of angular momentum; taking the scalar product with $\dot{\mathbf{r}}$ proves conservation of energy. The relative constants areThe total relative angular momentum and energy are $m_{\mathrm{red}}\mathbf{h}$ and $m_{\mathrm{red}}\epsilon$, respectively. The conserved eccentricity vector is $\mathbf{e}=\dot{\mathbf{r}}\times \mathbf{h}/\mu -\hat{\mathbf{r}}$. Its scalar product with $\hat{\mathbf{r}}$ gives $r=h^2/[\mu (1+e\cos f)]$, with $f$ the true anomaly. The periapsis and apoapsis distances are therefore $h^2/[\mu (1+e)]$ and $h^2/[\mu (1-e)]$. Their sum is twice the semi-major axis, soAt either apsis $v=h/r$; substituting an apsidal radius gives $\epsilon =-\mu /(2a)$. Thus the specific orbital energy implies the vis-viva equationThese are exact relative-coordinate formulas for an elliptic orbit, $0\le e<1$. In the barycentric frame the comet's speed is $M_{\star }/(M_{\star }+m)$ times the relative speed; the printed identification with the comet's stellar orbital speed uses $M_{\star }\gg m$.

In the weak, non-crossing, leading-eccentricity model the exterior resonant coefficient has positive sign: write its disturbing function as $R_{\mathrm{res}}=B_q{e}_c^q\cos {\varphi }_{bc}$, with $B_q>0$ for the orders considered. The physical feedback can be seen for first order. Before an outer conjunction near apoapsis, the faster inner planet pulls the outer one backwards; after passing it pulls forwards. On the outward side, the nearer, stronger earlier pull wins, reducing the outer semi-major axis and increasing its mean motion. Thus ${\varphi }_{bc}<\pi$ is driven upward. On the inward side, the later forward pull wins and drives ${\varphi }_{bc}>\pi$ downward. The paired impulses cancel at the symmetric centre.

For higher order, sum the feedback over all $q$ conjunction phases. The first surviving angular Fourier harmonic is $\sin {\varphi }_{bc}$; the lower Fourier harmonics cancel over the equally spaced directions. A local Hamiltonian reduction with canonical momentum $J$ hasThe negative curvature of the Kepler orbit Hamiltonian gives the displayed sign. Linear stability analysis at $\pi$ gives a restoring acceleration. The symmetric libration centres of an eccentricity resonance are thereforeThe corresponding conjunction phases relative to periapsis are $\{\pi \}$, $\{\pi /2,3\pi /2\}$ and $\{\pi /3,\pi ,5\pi /3\}$. They suppress repeatedly close passages and provide resonance protection.

These are expected centres of the specified leading-harmonic model, not centres fixed solely by the integer order for arbitrary orbital eccentricity and masses. Asymmetric resonant-argument libration can occur when higher Fourier harmonics matter. For example, if $R=A\cos \varphi +B\cos 2\varphi$, $A>0$, $B>A/4$, the symmetric $\pi$ equilibrium loses stability and stable centres satisfy $\cos \varphi =-A/(4B)$. Thus further dynamical information is needed beyond the PDF's unrestricted eccentric orbit to assert unique centres. Exterior asymmetric branches are documented by Winter and Murray, “Resonance and chaos. II. Exterior resonances and asymmetric libration”.

Now the eccentric planet $c$ is inside the circular perturber $d$. The interior eccentricity-type mean-motion resonance hasIn the mean-conjunction approximation, ${\varphi }_{cd}=n(\Lambda -{\varpi }_c)$, so ${\Lambda }_j={\varpi }_c+({\varphi }_{cd}+2\pi j)/n$. Successive conjunctions advance by $2\pi m/n$; coprime $m,n$ visit $n$ distinct phases.

The near-conjunction tangential pulls are reversed for an eccentric body inside its perturber. Just before conjunction the inner body is pulled forwards and afterwards backwards. Near periapsis, an inward-moving inner body experiences the stronger earlier forward pull; its semi-major axis grows, its mean motion falls, and ${\varphi }_{cd}$ increases because its coefficient of ${\lambda }_c$ is negative. On the outward side the feedback reverses, restoring the first-order centre zero. Higher-order feedback is strongest on the side closest to the external perturber, near apoapsis; the leading coefficient alternates in sign with order.

Writing $R_{\mathrm{res}}=B_n{e}_c^n\cos {\varphi }_{cd}$ gives $B_n<0$ for odd $n$ and $B_n>0$ for even $n$ in the leading non-crossing model. The same negative-curvature Hamiltonian argument therefore givesConjunction phases relative to periapsis are $\{0\}$, $\{\pi /2,3\pi /2\}$ and $\{0,2\pi /3,4\pi /3\}$, respectively, as shown in the lower row of the figure. As before, these are small-eccentricity symmetric predictions; additional resonant Fourier harmonics and coupled secular dynamics can change the stable branches.

For actual alignment of the eccentric body at true longitude $\Lambda$, the exact phase relation is instead ${\varphi }_{cd}=(m+n)x-my$, where $x=\Lambda -{\varpi }_c$ and $y=M_c(x)$. The simple equally spaced angles replace $y$ by $x$ at leading order.

Put ${\mu }_{\star }=G{M}_{\star }$. Conservation of specific orbital energy and specific angular momentum givesThe last comet expression is the vis-viva equation at the planet's orbit. The planet has only tangential velocity, while the comet's tangential component is $v_t=h/a_p$. Consequently its prograde relative velocity is exactly, before planetary deflection,In the near-parabolic encounter limit $Q\gg a_p$, this reduces toThe displayed approximation also requires $Q\gg a_p$. The assumptions $Q\gg q$ and $q\ll a_p$ alone do not guarantee it: an orbit with $Q=a_p$ encounters the planet near apoapsis and has $v_{\mathrm{rel}}\simeq v_p$, not $\sqrt{3}{v}_p$.

If all coplanar near-parabolic planet-crossing comets are admitted, so that $0\le q\le a_p$, the tangential component ranges from $-\sqrt{2}{v}_p$ to $+\sqrt{2}{v}_p$ when retrograde orbits are included. ThusFor prograde orbits alone the upper limit is $\sqrt{3}{v}_p$; for the stated deep-crossing limit $q\ll a_p$, both prograde and retrograde encounters approach $\sqrt{3}{v}_p$. These are the speeds outside the planet's gravitational well. The surface impact speed is $v_{\mathrm{imp}}=\sqrt{v_{\mathrm{rel}}^2+v_{\mathrm{esc},\mathrm{p}}^2}$, where $v_{\mathrm{esc},\mathrm{p}}$ is the planetary escape velocity; under negligible gravitational focusing the two speeds agree.

Assume negligible fragment release speed relative to a circular parent planetesimal. Its initial speed is $v_b=\sqrt{G{M}_{\star }/r_b}$, but radiation pressure reduces the grain's gravitational parameter to $(1-\beta )G{M}_{\star }$. Its specific orbital energy and specific angular momentum areFor $0\le \beta <1/2$, the resulting dust orbit released from a circular parent ring hasRelease is at periapsis: the inherited tangential speed exceeds the new circular speed. Low-$\beta$ grains remain near the dust birth ring; higher bound $\beta$ gives increasingly eccentric orbits extending well outside it. Uniform release longitude produces an approximately axisymmetric dust halo, with long residence times near apoapsis. At $\beta =1/2$ the conservative release orbit is parabolic; above this radiation-pressure blowout threshold it is initially unbound. For $\beta \ge 1$, zero or repulsive effective gravity also prevents a bound ellipse.

The left panel shows one release longitude; rotating it around the ring fills the halo. At longer times Poynting–Robertson drag also carries bound grains inward, so an initially outward-spanning release distribution does not imply an empty inner region in a continuously supplied, collision-free population.

In the proper-eccentricity annulus, averaging over random proper phases cancels the proper vector in the first-order angular variation. The mean line density is proportional to $1-e_f\cos \theta$, while the full ring width is independent of angle. Thus at corresponding radial locationsThere is a forced-apocentre enhancement in number surface density of a disk, together with the proper-excursion edge profile.

For a narrow aligned eccentric ring with locally constant $e_f$, its width has the same first-order factor $1-e_f\cos \theta$ as its line density. Their ratio is constant: the constant-eccentricity aligned ring has no first-order surface-density contrast. It still contains more particles per unit longitude at apoapsis, spread across a wider area.

The more general surface density contrast in an eccentric ring depends on the eccentricity gradient. With a smooth number distribution per unit semi-major axis and aligned apsides,Hence decreasing $e_f(a)$ can also make the aligned ring apocentre-dense; for $e_f\propto a^{-1}$ it has the same first-order contrast as the proper-phase case. A positive gradient gives a pericentre enhancement. The problem leaves this profile unspecified, so there is no universal observational distinction based solely on the two apse densities.

Resolved ring width and a measured forced-eccentricity profile can distinguish the constant-width and similar-ellipse constructions. Interpreting observed brightness additionally requires the grain temperature, illumination, emissivity, and viewing geometry; brightness is not itself particle number density. Inferences from one contrast alone can therefore be degenerate.

If the proper eccentricity is negligible and the semi-major axis ranges over $a\pm \Delta a$, particles lie on nearly aligned forced Kepler orbits. For a locally constant forced vector, these are geometrically similar ellipses, represented to first order by circles of radius $a^{\prime }$ centered at $-a^{\prime }{\mathbf{e}}_f$.

The boundaries of the aligned eccentric ring are thereforeIt is narrower at forced periapsis and wider at forced apoapsis, unlike the constant-width proper-eccentricity annulus. Its centers shift slightly between the two edges, as shown in the right panel of the figure. The condition $e_p/e_f\ll \Delta a/a$ makes the proper radial excursion negligible compared with the variation in forced center location.

The forced vector need not actually be constant across the interval. More generallyFor aligned apsides this becomes $2\Delta a[1-(e_f+a{e}_f^{\prime })\cos \theta ]$, provided the orbits remain nested. The radial gradient of forced eccentricity matters for the density comparison even in a narrow ring.

To first order in orbital eccentricity, the semi-minor axis is $a+O(e^2)$, while the ellipse's center is displaced from the star by $-a\mathbf{e}$. Thus each orbit is a circle of radius $a$ about $-az$, in complex spatial coordinates.

Write $z=z_f+e_p{e}^{i{\varpi }_p}$ at one secular epoch. The centers of the individual circles lie on a circle of radius $a{e}_p$ about $C_f=-a{z}_f$. Their union is the proper-eccentricity annulusIts full radial width is $2a{e}_p$ and its center is displaced by $a{e}_f$ toward the forced apoapsis direction. Relative to the star, if the forced longitude of periapsis is zero,The middle panel of the preceding figure shows this construction. For fully sampled orbital phases and proper phases, the annulus is occupied throughout, but its density need not be uniform; turning points in the proper radial excursion give integrable edge enhancements. The construction requires $e_f,e_p\ll 1$ and is an instantaneous secular geometry.

At the periapsis $r_1$ and apoapsis $r_2$, the given Kepler orbit satisfiesHenceTherefore

First follow the expression printed in the paper. Put $e=1-\delta$ and compare its numerator with $(1-\beta {)}^2$. To first order in $\delta$,Thus the grain is unbound, $e_d\ge 1$, whenSince the cometary apoapsis is $Q=a(1+e)=2a+O(\delta a)$,For completeness, the energy and angular-momentum-consistent formula derived in part (ii) instead gives $r_{\mathrm{blow}}=2\beta a\sim Q\beta$. The distinction is not a matter of approximation: it exposes the typographical inconsistency in the question.

The sudden potential change leaves the star's position and velocity unchanged. Its new energy is thereforeThus $E_1>0$ precisely when $r<a$. Parametrize the orbit by eccentric anomaly $\eta$:A phase-mixed orbit is uniform in the mean anomaly ${\theta }_r$. For $e>0$, the condition $r<a$ is $-\pi /2<\eta <\pi /2$ modulo one period. The corresponding mean-anomaly interval has lengthHenceThe energy increase $G(M_0-M_1)/r$ is largest near periapsis, so stars are more readily unbound there than near apoapsis. A perfectly circular orbit is the measure-zero marginal case $E_1=0$ at every phase.

For the 8:5 resonance, $p=5$ and $q=3$. Five planetesimal periods equal eight planetary periods, so the curve closes afterStarting at conjunction at apoapsis gives ${\varphi }_1=q\pi =3\pi \equiv \pi (\mathrm{mod}2\pi )$. In the frame rotating with the planet, let $\theta$ be the planetesimal's longitude relative to the fixed planet. At periapsis, $5\theta ={\varphi }_1$ modulo $2\pi$, so the five possible periapsis directions areAt apoapsis, ${\varphi }_1=5\theta +3\pi$, givingA sketch should therefore show a five-lobed rotating-frame rosette centered on the star, with the fixed planet on the $\theta =0$ ray, five outer turning points separated by $72^{\circ }$, and five inner turning points halfway between their rays. Successive orbits visit these points in resonance order before the fifth orbit closes the pattern.

Let $\psi ={\lambda }_d-{\lambda }_p$ be azimuth in the frame rotating with the planet, with the planet fixed at $\psi =0$. Since $n_d/n_p=4/5$, one synodic circuit contains four radial oscillations. To first order in $e_d$, choosing conjunction at apoapsis givesThe requested sketch is therefore a four-lobed closed curve centred on the star. For $e_d=0.2$, its radii range from $0.8a_d$ to $1.2a_d$; an outer lobe points along the star-planet line, and the planet lies at $(a_p,0)$. At that radial alignment $\varphi =\pi$, because ${\lambda }_d={\lambda }_p$ and ${\lambda }_d={\varpi }_d+\pi$ at apoapsis.

At exact conjunction the planet's force is radial and has no instantaneous torque, but the approach and departure do not cancel when conjunction is displaced from apoapsis. Write the conjunction longitude as ${\lambda }_c={\varpi }_d+\pi +\delta$. At conjunction,The resonant part of the given evolution law has ${\dot{a}}_d=-8e_d{a}_{d}D\sin \varphi$. Thus a conjunction after apoapsis, $\delta >0$, gives $\sin \varphi <0$ and transfers angular momentum to the dust, increasing $a_d$. A conjunction before apoapsis removes angular momentum, while one exactly at apoapsis has zero secular transfer by symmetry.

A higher orbital line density at apoapsis does not always give a higher surface density of a disk: the local ring width also changes. In a narrow aligned eccentric ring, neglecting disc twist and keeping first-order orbital eccentricity, ${\Sigma }_{\mathrm{peri}}\propto 1+ade/da$ and ${\Sigma }_{\mathrm{apo}}\propto 1-ade/da$. Constant eccentricity therefore gives constant surface density to this order.