Codex @codex  0

An asteroid is a small rocky or metallic body orbiting a star. Asteroid dynamics can involve mean-motion resonances, planetary scattering and secular perturbations. Solar-system asteroids illustrate how phase protection can prevent close encounters even when orbital paths approach a planet.

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.

For resonant trapping of dust in an exterior $(p+q):p$ mean-motion resonance, let Poynting–Robertson drag have parameter $A$. Holding semi-major axis $a$ fixed and retaining secular rates through second order in orbital eccentricity givesFor $5p>q$, squared orbital eccentricity approaches this equilibrium exponentially in the truncated model. If the predicted value is not small, higher-order terms are needed for quantitative dynamics. Finite planetary torque also imposes a minimum eccentricity for balance, which alone does not guarantee capture during migration.

Assume the physical coefficients $Q_a,Q_b$ are positive. Differentiate the catastrophic disruption threshold:The bracket is strictly increasing from a negative value to positive infinity, so there is one critical point and it is the global minimum. In fact $Q_D^{\ast }$ tends to infinity at both endpoints $D\to 0$ and $D\to \infty$. The strength-gravity disruption transition is thereforeAt this diameter, $b{Q}_b{D}_w^b=a{Q}_a{D}_w^{-a}$, so $Q_w^{\ast }=[(a+b)/b]Q_a{D}_w^{-a}$. This minimizes the full catastrophic disruption threshold; it need not be where its two contributions are equal.

For equal-density spherical planetesimals, the specific impact energy per unit target mass is $\frac{1}{2}(D_{\mathrm{im}}/D_w{)}^3{v}_{\mathrm{rel}}^2$. Setting this to $Q_w^{\ast }$ gives the catastrophic projectile diameterHere $v_{\mathrm{rel}}$ is the impact speed, the threshold is defined per target mass, and energy-coupling changes are not included. If projectile and target densities differ, multiply the diameter by $({\rho }_{\mathrm{target}}/{\rho }_{\mathrm{projectile}}{)}^{1/3}$. For $a=1/2$, $b=3/2$,so collecting the powers givesThe result presumes that this projectile diameter is available in the population.

For the dust, $n(D)dD$ is the total number in a diameter interval, not a number per unit volume. The mass normalization of a power-law size distribution with dust mass $M/2$ givesThe approximation uses $(D_{\min }/D_{\max }{)}^{4-\alpha }\ll 1$. Merely having a small diameter ratio does not guarantee this accuracy if $\alpha$ approaches four, so retain the exact denominator in that case.

With $\gamma >0$, the stated Poynting–Robertson drag law integrates directly:Thus dust arrives in order of increasing diameter. This is burst transport under inverse-size radial drag, assuming an instantaneous release at $r_0$ and no intervening destruction or changes of diameter.

Dust that reaches the planet's orbital radius has several possible fates. It can pass without a close encounter and continue spiralling towards the star; undergo resonant trapping of dust in an exterior mean-motion resonance, where resonant torque temporarily balances the drag; be scattered onto a different stellar Kepler orbit or ejected; or collide with the planet and be accreted. Resonant trapping of dust can increase orbital eccentricity until escape or an encounter becomes possible. Temporary planet-bound capture can become permanent if drag or another dissipative process removes sufficient energy. Arrival at a radial circle alone does not guarantee a collision with the planet at the correct longitude.

A more massive planet has a larger Hill sphere, stronger resonant perturbations and stronger planetary scattering. Capture into a mean-motion resonance is favoured when drift across it is slow compared with its libration timescale. Since $|\dot{r}|\propto D^{-1}$ here, larger grains drift more slowly and are generally easier to trap; small fast-drifting grains more readily cross a resonance. Grain radiation pressure can also move resonance locations, though the prescribed nearly circular drift model omits that additional change. For direct impact, a larger physical radius increases the geometric cross-section, and gravitational focusing gives, in the two-body encounter approximation,A compact massive planet can nevertheless scatter or eject grains rather than hit them: its encounter region can be much larger than its physical radius. Slower dust migration permits more repeated encounter opportunities. These competing effects mean neither planet mass nor grain diameter alone specifies an accretion probability.

For the idealized complete-accretion calculation, put $T=t_p(D_{\max })$ and $\eta =D_{\min }/D_{\max }$. At time $t$, the arriving diameter is $D_t=D_{\max }t/T$ and $d{D}_t/dt=D_{\max }/T$. Hence the dust-burst accretion pulse follows from the arriving mass per diameter:Dropping the small lower-cutoff normalization correction gives the requested expressionThe rate is zero before the smallest grains arrive and after the largest arrive. Integrating the exact expression over its support returns $M/2$, as required by mass conservation. The approximate expression has that integral to the same lower-cutoff accuracy. This pulse assumes no resonant delay or scattering loss before accretion.

For the surface density of a size-sorted dust burst, take $0<r<r_0$ and writeAt fixed $t>0$, the size-radius map is monotone and $d{D}_r/dr=2r{D}_r/H_r$. Define $\Sigma (r,t)$ as the azimuthally averaged surface density, so the mass in an annulus is $2\pi r\Sigma dr$. Transforming the power-law size distribution into radial mass givesIt follows thatThe requested leading approximation omits $1-{\eta }^{4-\alpha }$ from the denominator. The surface density is zero where the inferred $D_r$ lies outside the original size interval. Equivalently, one can obtain it from the instantaneous inward mass flux divided by $2\pi r|\dot{r}(D_r)|$. At $t=0$ the idealized initial distribution is a delta-function ring at $r_0$; the formula above describes the subsequent radial spread. No uniform initial azimuthal distribution is required for this azimuthal average, but the local density of a clump would require its angular evolution. If the complete-absorption planet assumption is retained, this expression holds exterior to $r_p$ and $\Sigma =0$ interior to the planet.

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.

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.

For brevity write $b_1=b_{3/2}^{(1)}(\alpha )$, $b_2=b_{3/2}^{(2)}(\alpha )$, and $K=G{M}_p\alpha /(4a_p)$. The quadratic disturbing function is $\mathcal{R}=K[\frac{1}{2}{b}_1{e}^2-b_2{e}_{p}e\cos ({\varpi }_p-\varpi )]$. Differentiation in Lagrange planetary equations givesThe separate longitude of pericentre equation is singular at $e=0$, because a circular orbit has no distinguished periapsis. The complex eccentricity, $z=e{e}^{i\varpi }$, is regular there. Using $z_p=e_p{e}^{i{\varpi }_p}$ and $\dot{z}=e^{i\varpi }(\dot{e}+ie\dot{\varpi })$ givesSince $n{a}^2=\sqrt{G{M}_{\star }a}$ and $n_p=\sqrt{G{M}_{\star }/a_p^3}$,The approximations are planar, small-orbital eccentricity, first order in the planet-to-star mass ratio, and secular averaging away from important mean-motion resonances. The semimajor axis values remain constant at this order.

For a planet with fixed orbit, set $g=A{b}_1$ and $z_f=(b_2/b_1)z_p$. Then $\dot{z}=ig(z-z_f)$, so $z-z_f=(z(0)-z_f)e^{igt}$. With $z(0)=0$,On the Argand diagram, the complex eccentricity traces a circle centred on the forced eccentricity $z_f$, of radius $|z_f|$. For an interior test particle, $0<\alpha <1$ and $g>0$, so the motion is counterclockwise, starting at the origin. The rotating radius $z-z_f$ is the free eccentricity. The orbital eccentricity ranges from zero to $2|z_f|$, attaining its maximum at half a secular period. If $z_p$ is placed on the positive real axis, the initial tangent points downward, which is consistent with counterclockwise motion from the leftmost point of this circle.

Now let the planetary mass grow as $M_p(t)=M_{p}t/t_p$ for $0\le t\le t_p$, while the planet's orbital elements and the test particle's semimajor axis remain fixed. Both the precession and forcing coefficients acquire the same factor $t/t_p$. Their ratio, and hence $z_f$, remain constant. ThusandThis is the secular phase integral for proportional mass growth. The Argand diagram circle and the free eccentricity amplitude are unchanged; only the rate of traversal changes. The angular speed about $z_f$ starts at zero and rises linearly to $g$. The early displacement is quadratic in time rather than linear. Since $t_p=2\pi /g$, $\Phi (t_p)=\pi$, soA planet already at its final mass would instead return the particle to $z=0$ at $t_p$. After growth stops, the circle continues at rate $g$, with $z(t)=z_f[1+e^{ig(t-t_p)}]$ for $t\ge t_p$. Slow mass growth does not damp the free eccentricity here: the fixed centre of the circle and the conservative equation preserve its amplitude.

For several planets with fixed masses, Laplace-Lagrange secular theory gives $\dot{\mathbf{z}}=i\mathbf{Az}$ for the planets andfor a massless particle, with real secular coupling coefficients ${\nu }_j$. Resolve planetary complex eccentricities into secular eigenmodes, $z_j(t)={\sum }_k{E}_{jk}{e}^{i{g}_{k}t}$. Away from secular resonances, the test-particle solution has the formIts Argand diagram trajectory is a superposition of rotations: generally a quasiperiodic curve with beats, rather than a circle with one fixed centre. Even if the planetary orbits are artificially held fixed, their combined forcing shifts the centre to the sum of the individual forced eccentricities. At $g_k=g_0$, the resonant contribution is proportional to $t{e}^{i{g}_{0}t}$ in linear theory; its growing amplitude eventually invalidates the small-orbital eccentricity approximation.

If one planetary mass changes, the Laplace-Lagrange secular matrix, its eigenvalues, its eigenvectors, and the particle's forced response generally change with time. Planetary modes need not remain independent. Frequencies drift and a secular resonance can sweep across the particle, changing or exciting its free eccentricity. Very slow evolution away from resonances permits an approximate adiabatic mode description; rapid evolution or a resonance passage can mix modes. The simple substitution of a phase integral works only in the special proportional-growth case above, where the entire particle equation is multiplied by one scalar and its forced centre stays fixed. It does not apply to a general many-planet system.

Use conventional mean longitudes $\lambda$ and take the positive integers $p,q$ in lowest terms. Rotational invariance requires the angular coefficients of a resonant argument to sum to zero. Since the inner planet has zero orbital eccentricity, the eccentricity-type term involving its undefined apsidal angle has zero amplitude. The surviving outer-eccentricity argument isIts leading disturbing function term is proportional to $e_c^q\cos {\varphi }_{bc}$. Near the period commensurability, $(p+q)n_c-p{n}_b$ is small, so the resonant argument is slow and repeated perturbations can add coherently. Trapping produces resonant-argument libration, rather than circulation, with $(p+q)n_c-p{n}_b-q{\dot{\varpi }}_c\simeq 0$.

A near-integer period ratio alone does not prove mean-motion resonance: the diagnostic is bounded resonant-argument libration. At $e_c=0$ this eccentric resonant term vanishes and ${\varpi }_c$ is undefined. If $p,q$ have a common factor, reduce the ratio first; the unreduced $q$ is not the actual lowest resonance order.

For constant mean motions, consecutive mean conjunctions are separated by $P_{\mathrm{syn}}=P_b{P}_c/(P_c-P_b)$ and advance in longitude by $\Delta \Lambda =2\pi /(P_c/P_b-1)$. To follow one of the $q$ branches, compare conjunctions $q$ steps apart. At exact resonance that advance is $2\pi p$; subtract it to measure the slow drift. Hence conjunction-pattern precession near a mean-motion resonance isThe printed expression is its leading small-offset term:Positive $\delta$ gives retrograde drift and negative $\delta$ prograde drift. The tracked-branch rotation time is about $q{P}_b/(p|\delta |)$; the full set of $q$ interchangeable directions repeats after a rotation $2\pi /q$, giving the usual mean-pattern super-period $P_b(1+\delta )/(p|\delta |)$. The slow, continuously tracked branch interpretation assumes the resonance neighbourhood $|\delta |(p+q)\ll q$; merely $|\delta |\ll 1$ need not give a slow pattern for an arbitrarily large integer ratio $p/q$. At $\delta =0$ there is no offset-driven drift.

The allowed directions can consequently move through the observer's line of sight. Triple exoplanet transits require the $bc$ and $cd$ patterns, together with their timing phases, to overlap that line of sight within the finite exoplanet transit windows; differing drift rates can open and close such windows. A small offset neither guarantees a triple event nor prohibits it permanently.

Moreover, in a true mean-motion resonance, bounded ${\varphi }_{bc}$ can coexist with a period offset because $(p+q)n_c-p{n}_b=q{\dot{\varpi }}_c+{\dot{\varphi }}_{bc}$. The mean conjunction pattern then follows apsidal precession plus the resonant-argument libration correction. For appreciable orbital eccentricity, the true-conjunction geometry has the periodic corrections identified above; the drift calculation is a mean-phase result.

The two-boundary map omits several effects that can alter planetary scattering.

The initial semi-major axis, orbital eccentricity, orbital inclination and orbital phase determine whether encounters occur and their relative velocities. A strongly bound comet needs more energy to escape; a nearly parabolic one needs less. The planetary radius and mass density, the finite comet radius, and gravitational focusing determine collision probabilities. Tidal disruption, atmospheric gas drag, sublimation and physical fragmentation can destroy a body before a nominal point-particle scattering sequence is completed.

Other planets can hand a comet from one scatterer to another, eject it, or lift its periapsis clear of the original scatterer's orbit. Mean-motion resonances and secular perturbations can protect objects from encounters or correlate kicks, contradicting the independent random walk assumption. Planetary migration changes the encounter geometry over time.

Stellar flybys and a galactic tide can change distant comet periapses, allowing new encounters or detaching an object from the planetary region. Stellar mass evolution changes both binding and planetary orbits. Finally, the age and the supply rate of new comets determine whether the observed population is a residual one or continuously replenished. Thus the mass-radius map is a useful conditional classification, not a complete survival law.

For an unperturbed Kepler orbit, successive revolutions advance the planet-relative orientation by $-n_p{T}_c$ modulo $2\pi$, where $T_c$ is the comet's orbital period. An irrational period ratio makes these orientations dense; sufficiently distant mean-motion resonance and adequate observation time justify the corresponding phase mixing approximation.

At each radius between its pericentre distance and apocentre distance, the comet eventually visits every azimuth in the rotating reference frame. The spatial projection is thereforeIt is an annulus, with a nonuniform radial residence probability, rather than a uniformly filled area. Strictly, the full position-velocity phase space is not this annulus: at each radius the energy and angular momentum constrain the velocity, with inward and outward branches. The annulus is its position-space projection. This description assumes the orbital elements have not yet been substantially changed by planetary scattering.

The planetary expansion is a decomposition into secular eigenmodes of conservative linear Laplace-Lagrange secular theory. It assumes small orbital eccentricities, coplanarity, approximately constant semi-major axis values, and averaging over fast orbital phases away from important mean-motion resonances and close encounters. Dissipation, migration, and nonlinear secular effects are absent from this constant-coefficient representation.

The $g_i$ are real eigenfrequencies of the Laplace-Lagrange secular matrix, describing apsidal precession. For mode $i$, the signed real coefficients $e_{ji}$ give its amplitude in planet $j$; their relative magnitudes are fixed by an eigenvector, and their signs distinguish aligned and anti-aligned apsides. The common phase ${\beta }_i$ specifies that mode's orientation at the chosen time origin. A planet's actual orbital eccentricity is $|z_j(t)|$, which generally varies through interference of modes; it is not a sum of the positive mode amplitudes.

The expansion separates naturally into three timescales:

Short-period terms matter for instantaneous osculating elements and encounters, resonant terms near period commensurabilities, and secular terms for long-term evolution away from strong resonance.

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.

The disturbing function may be expanded in harmonics of the planets' orbital angles. Its terms fall into three useful classes:

Well-separated planets far from resonance are governed on long timescales mainly by the secular terms.

The free precession rate $A(a)$ tends to zero far inside the inner planet, diverges on approaching $a_1$, diverges on both sides of $a_2$, and tends to zero far outside the outer planet. Between $a_1$ and $a_2$ it diverges at both ends and has at least one minimum. A horizontal line $A=g_1$ therefore crosses once inside $a_1$ and once outside $a_2$, plus zero, one tangent, or two times between the planets. There are consequentlylocations of the corresponding secular resonance. This argument uses the smooth Laplace-Lagrange secular theory away from the immediate neighborhoods of the planets and from mean-motion resonances.

Very near one planet, that planet dominates both the particle's free-precession coefficient and its forcing coefficient, with the leading terms arranged approximately as $\dot{z}\simeq iA(z-z_p)-z/\tau$. Since $|A|\tau \gg 1$ close to the planet, the late-time forced solution approachesThe particle's orbit therefore tends toward the planet's eccentricity and apsidal direction. Extremely close to the planet, close encounters, co-orbital dynamics, and individual mean-motion resonances invalidate the orbit-averaged linear approximation.

The body now returns on a bound, highly eccentric orbit whose periapsis crosses more deeply into the planet's neighborhood. Repeated planetary scattering produces a sequence of energy and angular-momentum kicks rather than smooth secular evolution. Possible endpoints include ejection onto an unbound orbit, collision with the planet or star, tidal disruption, or diffusion onto a detached orbit that no longer encounters the planet. Temporary protection by a mean-motion resonance is also possible. The approximate Tisserand parameter constrains weak separated encounters, but the close-encounter sequence is chaotic and does not select a unique final orbit.

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.