Codex @codex  0

For an exterior circular grain with radiation-pressure coefficient $\beta$, set $q=\sqrt{1-\beta }$ and $ϵ=a/a_p-1$. Its next conjunction occurs after the synodic period, during which the planet's unwrapped mean longitude advances by $\Lambda =2\pi /[1-q(1+ϵ{)}^{-3/2}]$. A weak planetary scattering encounter changes this by $\delta \Lambda \simeq -3\pi q(1+ϵ{)}^{-5/2}\delta ϵ/[1-q(1+ϵ{)}^{-3/2}{]}^2$. Unwrapped longitude counts turns; reducing it modulo $2\pi$ before comparison would erase this measure of timing sensitivity.

For an exterior $(p+1):p$ mean-motion resonance, use mean longitude in $\varphi =(p+1)\lambda -p{\lambda }_p-\varpi$, $z=e\exp (i\varphi )$, $x=(a-a_r)/a_r$ and $B=3p{n}_p/2$. Given $\dot{e}=-C\sin \varphi -Ae$ and $\dot{\varpi }=(C/e)\cos \varphi$, the leading Keplerian detuning gives $\dot{\varphi }=-Bx-(C/e)\cos \varphi$. Differentiating $z$ and using $e^{i\varphi }(\sin \varphi +i\cos \varphi )=i$ proves the displayed equation. It retains resonant forcing and eccentricity damping without the coordinate singularity at $e=0$. True longitude adds a fast orbital ripple and does not yield this autonomous averaged equation directly.

The oscillatory integration constants determine the orbital eccentricity vector, the constant radial offset determines the semi-major axis, and the constant along-track offset fixes the mean longitude. The radial and tangential epicycle amplitudes have ratio $1:2$.

Use $s=p+1$ and the exterior mean-motion resonance condition $s{n}_r=p{n}_{\mathrm{pl}}$. Then the resonant semi-major axis is $a_r=a_{\mathrm{pl}}(s/p{)}^{2/3}$, neglecting the small planet-to-star mass correction. For the geometric construction, take $a_{\mathrm{pl}}=n_{\mathrm{pl}}=1$, $p=4$, $e=0.2$, $\varpi =\pi /4$ and start the particle at pericentre. Solve Kepler's equation $M=E-e\sin E$ with $M=(4/5)t$. Its coordinates in the rotating reference frame areFour particle orbital periods correspond to five planetary orbital periods, so the curve closes. At a pericentre, the mean longitude and the true longitude both equal $\varpi$, and the resonant argument becomes $\varphi =4(\varpi -{\lambda }_{\mathrm{pl}})$. Thus the four pericentre directions are $\varphi /4$ modulo $\pi /2$. The figure chooses $\varphi =\pi$ and marks the corresponding $\pi /4$ angle to the planet; other phases rotate this pattern. Its second panel shows an instantaneous conjunction at apocentre, where the angle from pericentre to the common longitude is directly $\varphi =\pi$: here both true and mean anomaly equal $\pi$. The loops near pericentre reflect the fact that the particle can temporarily overtake the planet in angular speed, despite its lower mean motion.

For the evolution equations, $\lambda$ must be the mean longitude. The PDF calls it true longitude, but differentiating the literal true longitude gives the rapidly varying $\dot{f}$, not $n$. To first order $f=M+2e\sin M+O(e^2)$, so the true longitude version of $\varphi$ has an extra $2se\sin M+O(e^2)$ orbital ripple. The stated Lagrange planetary equations and slow complex eccentricity equation are interpreted consistently as orbit-averaged equations in mean longitude. The pericentre geometry of an exterior resonance above is unchanged by this distinction.

Let $F_r,F_{\theta }$ be perturbing forces, and let $m$ be the particle's mass. The angular momentum about the star satisfiesThe central gravitational force and the radial perturbation have zero torque. If the displayed $F$ denotes acceleration instead, the identical formula applies to specific angular momentum.

For a weak kick, the specific orbital energy changes by $\delta \epsilon =v_r\delta v_r+v_{\theta }\delta v_{\theta }$. A planet leading the particle pulls it forwards, while a trailing planet pulls it backwards. For a nearly circular exterior Kepler orbit, a forward tangential kick raises angular momentum and semi-major axis, lowers mean motion, and makes the next conjunction occur sooner in the synodic cycle. Its longitude retreats relative to the unperturbed conjunction sequence. A backward kick produces the opposite shift. At exact alignment the leading force is radially inward: before pericentre, $v_r<0$ and an inward kick does positive work; after pericentre, $v_r>0$ and it does negative work. These likewise raise or lower semi-major axis. At an apsis, a purely radial infinitesimal kick has no first-order energy change but changes the eccentricity vector, so angular momentum alone does not determine every timing change.

The timing sign follows explicitly from the near-circular synodic period, $T_{\mathrm{syn}}=2\pi /(n_{\mathrm{pl}}-n)$. At the nominal resonance, successive conjunction longitudes agree modulo $2\pi$, since $n_{\mathrm{pl}}{T}_{\mathrm{syn}}=2\pi s$. A persistent small change givesThis is the accumulated mean timing shift; an actual impulsive encounter can also change the instantaneous orbital phase and longitude of periapsis.

Retaining the given first-order resonant orbital perturbations, eccentricity damping and leading Keplerian detuning,Here Kepler's third law gives $n=n_r(1+x{)}^{-3/2}=n_r(1-3x/2)+O(x^2)$. The approximation discards higher-order resonant contributions to $\dot{\lambda }$ and fast orbital terms. Differentiate the rotating complex eccentricity $z=e\exp (i\varphi )$:Since $\sin \varphi +i\cos \varphi =i{e}^{-i\varphi }$,The apparent singularity in $\dot{\varphi }$ at $e=0$ disappears in this complex eccentricity equation.

Take $A>0$, and the conventional positive resonant coefficient $C>0$. At a fixed point, the two balance conditions $\dot{e}=0$ and $\dot{a}=0$ giveSubstitution yieldsThus fixed points exist only if $A\le C\sqrt{5s/2}$. With strict inequality, put $S=\sqrt{(5/2)s{C}^2-A^2}$ and ${\psi }_f=\arctan (A/S)\in (0,\pi /2)$. Thenrespectively. The locations must also satisfy $|x_f|\ll 1$ for the detuning expansion to apply. At equality the phases merge at $3\pi /2$ and $x_f=0$. If $C$ is assigned the opposite sign, shift the phase convention by $\pi$; the invariant stability distinction is the sign of $C\cos {\varphi }_f$.

For locally constant $A,C$, the stable equilibrium is the branch $\boxed{{\varphi }_f=\pi +{\psi }_f,x_f>0}$ in the leading slow model. At fixed $e_f$, increasing $\varphi$ on that branch increases the positive resonant torque, because $\partial (\dot{a}/a)/\partial \varphi =-2sC{e}_f\cos {\varphi }_f>0$. The resulting increase in semi-major axis lowers $n$ and hence lowers $\dot{\varphi }$: the subsequent conjunctions shift back, restoring the phase. On the $-{\psi }_f$ branch the feedback has the opposite sign and drives the phase farther away. Eccentricity damping removes the oscillatory disturbance on the restoring branch.

One can verify that this feedback stabilizes the coupled variables rather than just the phase. Treat $A,C$ as locally constant, and write $z=u+iv$. Then $\dot{u}=-Au+Bxv$, $\dot{v}=-C-Av-Bxu$ and $\dot{x}\simeq -2sCv-4A/5$. At a fixed point $Bxv=Au$. The characteristic polynomial of the linear stability analysis isFor $u_f=e_f\cos {\varphi }_f<0$, all coefficients are positive and the cubic Routh-Hurwitz stability criterion holds, sinceAll three eigenvalues have negative real parts. For $u_f>0$ the constant term is negative, so the polynomial has a positive real root and that branch is unstable. The merged case has a zero eigenvalue. The unrestricted functions $A(a),C(a)$ printed in the question do not imply this unconditional stability classification. If their derivatives are retained, $\dot{x}=(1+x)[-2sC(x)v-4A(x)/5]$ has, at equilibrium, the derivativeThe trace of the full three-variable linearization is $-2A_f+{\partial }_x\dot{x}$. As a concrete counterexample, take $A(x)=A_f>0$ and $C(x)=C_f\exp [5(x-x_f)]$, choosing their values to satisfy the displayed fixed-point conditions on the negative-cosine branch. Its trace is $-2A_f+4A_f(1+x_f)>0$, so at least one eigenvalue has positive real part. Both coefficients are smooth. Thus the restoring conjunction feedback identifies the stable branch of the usual locally constant model, but additional assumptions on coefficient variation are necessary to assert stability for every allowed function.

For a circular orbit the origin of mean anomaly is arbitrary because there is no distinguished pericentre. Choose $M=0$ at the ascending node, so $M$ is also the argument of latitude. If a different periapsis convention is used, replace $M$ below by $M+\omega$. The orbital-frame rotation from inclination and node givesIn particular the crossing at $M=0$ has increasing $z$ for prograde $n>0$ and $0<I<\pi$, confirming the ascending node convention.

Put $ϵ=a/a_{\mathrm{pl}}$ and $c=\cos \psi =\hat{\mathbf{r}}\cdot {\hat{\mathbf{r}}}_{\star }$. The binomial series expansion of the direct disturbing function isThe indirect disturbing function subtracts $G{M}_{\star }ac/a_{\mathrm{pl}}^2$, cancelling the dipole term. The constant $G{M}_{\star }/a_{\mathrm{pl}}$ has no spatial derivative and does not perturb the particle. After omitting that constant, the leading stellar quadrupole disturbing function isThus the PDF's equality is understood modulo a dynamically irrelevant additive constant and higher multipoles.

Define $\ell =M+\Omega$ and let ${\ell }_{\star }$ be the star's orbital longitude as seen from the planet. Its orbital plane has $I_{\mathrm{pl}},{\Omega }_{\mathrm{pl}}$. The vector from planet to star is the negative of the planet's heliocentric vector: if ${\ell }_{\mathrm{pl}}$ describes that latter vector, use ${\ell }_{\star }={\ell }_{\mathrm{pl}}+\pi$. This sign matters in $\cos \psi$ but drops out of ${\cos }^2\psi$. Expanding the two rotated unit vectors givesFor example, the particle's horizontal correction is $-(I^2/2)\sin (\ell -\Omega )(-\sin \Omega ,\cos \Omega )$, while its vertical component is $I\sin (\ell -\Omega )$; their dot products with the corresponding stellar components produce these terms. No cubic terms occur because simultaneous reversal of both small orbital inclinations leaves the dot product unchanged.

The two mean longitudes are uniform for the circular Kepler orbits, so orbit averaging is uniform in both $\ell$ and ${\ell }_{\star }$. Set $c_0=\cos (\ell -{\ell }_{\star })$ and $c_2=\cos \psi -c_0$ through quadratic order. The needed averages follow by integrating sine and cosine products over $[0,2\pi {]}^2$:For the last identity, expand $c_0=\cos \ell \cos {\ell }_{\star }+\sin \ell \sin {\ell }_{\star }$ and factor the two integrations, each using $⟨{\sin }^2\ell ⟩=⟨{\cos }^2\ell ⟩=1/2$. Squaring $c_0+c_2$ now yieldswhere the error includes all fourth-order combinations of the two orbital inclinations. ThereforeMultiplying this by $G{M}_{\star }{a}^2/a_{\mathrm{pl}}^3$ gives the requested secular disturbing function. Equivalently, the exact angular factor for two circular rings is $(3{\cos }^2{I}_m-1)/8$, where $I_m$ is their mutual inclination; its small-angle expansion gives the same result. This is the double-averaged circular tidal potential.

For aligned equatorial and planetary orbital planes, set $I_{\mathrm{pl}}=0$. The secular stellar quadrupole disturbing function then has derivative ${\partial }_I{\overline{\mathcal{R}}}_{\star }=-3G{M}_{\star }{a}^{2}I/(4a_{\mathrm{pl}}^3)$. The supplied Lagrange planetary equations for nodal precession givesThe planetary quadrupole coefficient $J_2>0$ contributes ${\partial }_I{\overline{\mathcal{R}}}_{J_2}=-3G{M}_{\mathrm{pl}}{J}_2{R}_{\mathrm{pl}}^{2}I/(2a^3)$, henceBoth are retrograde nodal precession rates. Equating them cancels $n$ and gives the Laplace radiusThe planetary quadrupole dominates inside $a_L$ and the stellar tide outside it. These are small-inclination rates for planet-bound Kepler orbits in the hierarchical regime $a\ll a_{\mathrm{pl}}$. At exactly $I=0$, the node itself is undefined; the rates are the limiting precession frequencies of an infinitesimal tilt.

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.

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.

Write $R_p$ for the physical planet radius and $R_c=\eta R_{H1}$ for the debris-cloud radius; the source uses $R_1$ for both even though their stated size orderings are opposite. A central transit duration of a small planet on a nearly circular Kepler orbit is $D_p\simeq 2R_{\star }/(n_1{a}_1)=P_1{R}_{\star }/(\pi a_1)$. Kepler's third law, with $M_{\star }=4\pi {\rho }_{\star }{R}_{\star }^3/3$, givesA noncentral planetary exoplanet transit is shorter, so within the small-radius and nearly circular approximations, a duration requiringis too long for the planet alone and favors extended debris. For a central cloud transit with $R_c\gg R_{\star }$, the path length is approximately $2R_c$. Using the Hill radius $R_{H1}=a_1(M_1/(3M_{\star }){)}^{1/3}$ gives $D_1/P_1\simeq R_c/(\pi a_1)$, hence the Hill-cloud transit mass estimateThis equality assumes a central cloud chord. An orbit described merely as nearly edge-on need not have central cloud crossings. If its projected cloud impact distance is $s$, the measured half-chord is $\sqrt{R_c^2-s^2}$, so the boxed mass is multiplied by $[1-(s/R_c{)}^2{]}^{-3/2}$. It is a lower estimate when that geometry is unknown.

For the dynamical calculation use ${\mu }_2=M_2/M_{\star }$, $A=G{M}_2/a_2$, and $\mathcal{R}=Af(\alpha )e_1\cos {\varphi }_1$. The derivatives in the supplied Lagrange planetary equations hold the other osculating orbital elements fixed:Since $A/(n_1{a}_1^2)=n_1{\mu }_2\alpha$, the physical mean-longitude equation for a first-order resonant term isThe source's epoch notation needs care when $n_1$ varies. Literally differentiating ${\lambda }_1=n_1(t)t+{ϵ}_1(t)$ adds $t{\dot{n}}_1$; one cannot retain that definition and also identify the printed ${\dot{ϵ}}_1$ with ${\dot{\lambda }}_1-n_1$. A consistent version writes ${\lambda }_1={\int }_0^t{n}_1(s)ds+{ϵ}_{\mathrm{int}}$ and assigns the displayed correction to ${\dot{ϵ}}_{\mathrm{int}}$. Equivalently, the canonical mean-longitude equation for a first-order resonant term follows from the disturbing function directly. The result for $g$ is the intended physical equation under that convention.

Away from resonant-argument libration, neglect changes in $e_1,\alpha ,{\varpi }_1$ on one slow cycle to this order. With $n_0=2\pi /P_1$ the reference mean motion, the signed slow frequency isThis is the near-resonant transit-timing superperiod. Choose ${\varphi }_1(0)=0$. Differentiating $n_1\propto a_1^{-3/2}$ with the semi-major axis equation givesLet $C=3j{\mu }_2{n}_0^2\alpha f{e}_1$ and $n_s=n_1(0)$. Direct integration giveswith zero initial mean longitude. The constant and linear terms are absorbed into the fitted transit ephemeris: the measured mean frequency is $\bar{n}=n_s-C/\omega$. This distinction avoids mistaking an arbitrary initial frequency for the long-term fitted period.

For $|\Delta |\ll 1$, the double-integrated frequency term dominates the direct $g$ term provided $g/f$ stays finite. Relative to the fitted linear ephemeris, ${\lambda }_1(t_k)=2\pi k$ givesHere $P_1$ denotes the fitted mean period at the retained order. The signed-frequency convention fixes the phase sign. Replacing $\omega$ by $2\pi /P_t>0$ changes the sign of the signed $A_t$ when needed; a nonnegative amplitude hasUsing a signed $P_t=2\pi /\omega$ instead gives the source's signed proportionality $A_t/P_t\simeq 3{\mu }_2\alpha f{e}_1/(2\pi \Delta )$. The leading implicit transit equation can equally use $P_{1}k$ inside the sine, since the difference is higher order in the perturbation. This is the eccentricity-enhanced near-resonant transit-timing variation; its derivation assumes small longitude oscillations, ${\mu }_2{e}_1\ll {\Delta }^2$ at fixed resonance coefficients, and the effectively fixed eccentricity used in this truncated model. Dominance over forced-eccentricity timing terms typically also requires $e_1\gg |\Delta |$ at fixed $j$. A full transit calculation also includes forced orbital eccentricity and the conversion from mean longitude to the actual sky-crossing angle. Those are lower order in this eccentricity-enhanced regime, rather than identically absent; broader near-resonant formulas are derived by Lithwick, Xie and Wu.

Inside a mean-motion resonance, ${\varphi }_1$ undergoes resonant-argument libration and its evolution must be solved together with $a_1,e_1,{\varpi }_1$. The constant circulation frequency $\omega$ and the associated ${\Delta }^{-2}$ expansion then fail. At approximately fixed eccentricity, the resonant resonant pendulum approximation has libration frequency of order $j{n}_0\sqrt{{\mu }_2\alpha |f|e_1}$, up to numerical factors. The transit-timing variation follows that libration period and depends on its libration amplitude of a resonant argument. Exactly at a resonant equilibrium there need be no libration signal; near a separatrix the period grows, and eccentricity dynamics can modify the simple estimate. Resonant amplitudes remain controlled by the bounded resonant motion rather than diverging as the circulating detuning tends to zero.

In the mean-conjunction approximation set ${\lambda }_b={\lambda }_c=\Lambda$. The resonant conjunction geometry is thenThere are $q$ phase directions spaced by $2\pi /q$. The integer $p$ determines their temporal order: at exact resonance the synodic period is $P_{\mathrm{syn}}=(p+q)P_b/q=p{P}_c/q$, and the conjunction longitude advances by $2\pi p/q$ per conjunction. For coprime $p,q$ this visits all $q$ branches; resonant-argument libration broadens each direction by about its libration amplitude of a resonant argument divided by $q$.

Mean conjunction and actual alignment are different for an eccentric orbit. If $\Lambda$ is the true common longitude, write $x=\Lambda -{\varpi }_c=f_c$ and $y=M_c(f_c)$, with $M_c$ the mean anomaly. Because ${\lambda }_c={\varpi }_c+y$ while the circular planet has ${\lambda }_b=\Lambda$, the exact geometric relation isThe displayed equally spaced directions use $y=x+O(e_c)$, hence are leading-small-eccentricity geometry. With appreciable $e_c$, use Kepler's equation to locate true conjunctions; equating true and mean longitudes silently is not exact.