Planetary system dynamics

Planetary system dynamics explains the changing orbits of planets and smaller bodies under gravity and dissipative forces. A Kepler orbit supplies a useful two-body reference; interactions with other bodies produce secular perturbations, orbital resonances and close encounters.

Conservative interactions exchange energy and angular momentum within the system, while tides, gas drag and radiation forces can dissipate or transfer them to other reservoirs. Comparing orbital times, migration times and collision times identifies which approximation applies. The resulting theory connects present configurations with possible formation histories without assuming that every observed orbit has stayed fixed since birth.

A temporary reversal in a planet's apparent direction across the background sky, arising from the observer's changing position. For an outer planet observed from the Earth, the Earth overtakes it along its smaller, faster orbit. Different apparent reversals for two planets can produce three apparent conjunctions around one Sun-centred encounter, even when both follow prograde circular orbits.

The systematic gravitational drag experienced by a massive body interacting with a population of lighter bodies. In a planetesimal swarm, it tends to damp a massive embryo's random motion while transferring random kinetic energy to the smaller bodies. Its detailed strength depends on the background velocity distribution and geometry.

The stellar density and Newtonian gravitational potential perturbation induced around a moving massive object. The stellar response includes both a forced component that follows the object and switch-on collective transients. In a homogeneous stable damped model, each late-time spatial harmonic is modified by the gravitational dielectric function at the object's forcing frequency. The wake's force back on the object produces dynamical friction.

Orbit crossing occurs when two geometrical orbital paths intersect. It permits collisions or close encounters when bodies also reach the intersection at compatible times; crossing alone does not imply a collision. For nearly circular coplanar nested orbits, a useful local test follows by differentiating their radial positions at fixed longitude.

For a cold coplanar disk with complex eccentricity $z(a,t)$, the first-order radial shape is $r=a[1-\mathrm{Re}(z{e}^{-i\theta })]$. Infinitesimally neighboring orbits first cross when $|z+a{\partial }_{a}z|=1$. Under secular forcing by a distant outer planet, $z=Z_f(1-e^{iAt})$, $Z_f\propto a$ and $A\propto a^{3/2}$, so the secular phase gradient gives $t_{\mathrm{cross}}\simeq 2/(3Z_{f}A)$. The finite-neighbor version requires the accumulated phase difference to remain small when replacing a difference by a derivative.

Isotropic stellar mass loss changes the central gravitational parameter $\mu (t)$ while preserving the central direction of the force, provided the escaping material causes no recoil or drag on the orbiting body. The relative equation is $\ddot{\mathbf{r}}=-\mu (t)\mathbf{r}/r^3$. Specific angular momentum remains exactly constant; specific orbital energy obeys $\dot{C}=-\dot{\mu }/r$. Orbital expansion is adiabatic only when the mass-loss law varies slowly on the orbital period.

Suppose an initially circular Kepler orbit of radius $R$ surrounds central mass $M$, and a fraction $0<f<1/2$ is removed quickly compared with the orbital period. Assume no recoil or wind drag, so the satellite's position and velocity are continuous through the change. Its specific angular momentum stays $\sqrt{GMR}$, while its new specific orbital energy is $-GM(1-2f)/(2R)$. The resulting elliptic Kepler orbit has semi-major axis $a=(1-f)R/(1-2f)$, pericentre $r_p=R$, apocentre $r_a=R/(1-2f)$ and orbital eccentricity $e=f/(1-f)$. The limiting loss fraction $f=1/2$ gives a parabolic Kepler orbit; larger losses unbind the initially circular satellite. Slow removal instead produces adiabatic orbital expansion under isotropic mass loss, which preserves circularity to leading adiabatic order.

Time averaging over a Kepler orbit gives $⟨\cos f⟩=-e$, as follows from $dt=(1-e\cos E)dE/n$ and $\cos f=(\cos E-e)/(1-e\cos E)$. Consequently slow isotropic stellar mass loss preserves the secular orbital eccentricity. Conserved specific angular momentum then gives $\mu a=\mathrm{constant}$, so pericentre, apocentre and semimajor axis expand inversely with the remaining gravitational parameter.

The signed instantaneous contribution of apsidal precession to azimuthal motion during isotropic stellar mass loss is $\dot{\varpi }/\dot{\theta }=-(\dot{\mu }/\mu n)(1-e^2{)}^{3/2}\sin f/[e(1+e\cos f{)}^2]$. Its maximum absolute magnitude is small when $n|\mu /\dot{\mu }|\gg F(e)$, where $F(e)\sim 1/e$ near a circular orbit and $F(e)\to 3\sqrt{3}/4$ as $e\to 1$. The low-orbital eccentricity coordinate singularity can be avoided by evolving the eccentricity vector.

For isotropic stellar mass loss the eccentricity vector obeys $\dot{\mathbf{e}}=-(\dot{\mu }/\mu )(\mathbf{e}+\hat{\mathbf{r}})$. With true anomaly $f$ and $e>0$, this gives $\dot{e}=-(\dot{\mu }/\mu )(e+\cos f)$ and $\dot{\varpi }=-(\dot{\mu }/\mu )\sin f/e$. The osculating pericentre satisfies $\dot{q}/q=-(\dot{\mu }/\mu )(1-\cos f)/(1+e)$, hence cannot decrease during mass loss and is instantaneously stationary at pericentre. The vector form is nonsingular at a circular orbit.

Anisotropic thermal radiation produces recoil that changes the semi-major axis of a rotating solid body. Finite thermal inertia provides the temperature lag. Large bodies have small acceleration per unit mass, while very small bodies become nearly isothermal. The Yarkovsky effect can supply bodies to dynamical escape regions.

Drift driven by the Yarkovsky effect can carry a body out of a belt or into an unstable orbital region. A useful two-limit removal law is $t_y=BD+C/D$, with slow escape at both large and small diameters.

Let $n_b(D)$ be the differential number of belt objects per size interval, and $t_y(D)$ their removal timescale through Yarkovsky removal. Suppose removed objects feed an external population with size-independent mean residence time $T$. A simple population balance isFor a maintained belt inventory, the steady solution is $n_{\mathrm{out}}=T{n}_b/t_y$. Starting from an empty external population adds the transient factor $1-e^{-t/T}$. The steady formula therefore requires a sufficiently long-lived supply; it is not an instantaneous identity for arbitrary initial data.

If $n_b\propto D^{-p}$ and $t_y\propto D^a$, the steady differential number spectrum is $n_{\mathrm{out}}\propto D^{-(p+a)}$. If the belt is itself in a removal-dominated steady state with fragment injection rate $S(D)$, its balance instead gives $n_b/t_y\simeq S$. The external spectrum then follows $TS(D)$. These two balances explain why a removal law can either change a supplied population's size slope or disappear from its final steady spectrum.

For $t_c=(A/M)D^{1/2}$, the ratio $t_y/t_c$ has its minimum at $D_{\ast }=\sqrt{3C/B}$. Removal beats collisions somewhere exactly when $M<M_y=(3^{3/4}/4)A{B}^{-3/4}{C}^{-1/4}$, provided the relevant diameters are present. This criterion uses the collision-only cascade as a reference.

The law $t_y=BD+C/D$ for positive $B,C$ has a minimum at $D_0=\sqrt{C/B}$ with value $2\sqrt{BC}$. This size differs from the minimum of the ratio of removal and collision times.

A passing star perturbs a planetary or cometary system through its differential Newtonian gravitational field. A distant encounter is predominantly tidal; a sufficiently close one can substantially change specific orbital energy and angular momentum, eject a comet, or shift its periapsis into or out of the planetary region.

Planetary migration changes a planet's orbital semimajor axis through torques from a disc or exchanges of energy and angular momentum with other bodies, often followed by tidal dissipation. Short-period giant planets, eccentric progenitors, resonant architectures, and stellar obliquity provide complementary constraints rather than one unique migration diagnostic.

Orbital migration of a low-mass embedded planet driven by its interaction with disk gas, without a deep gap. Spiral wave torques and corotation-region torques exchange orbital angular momentum. The sign depends on the disk gradients and thermodynamics; this regime is not restricted to inward migration.

The scale $M_{\mathrm{th}}=h^3{M}_{\star }$ makes the Bondi radius comparable with the disk scale height. It organizes when planet-induced disturbances become nonlinear; it is not identical to an exact universal gap-opening mass.

For a fixed-mass planet in a circular orbit of radius $a$, let $q=M_p/M_{\star }$ be the planet-to-star mass ratio, $\Sigma$ the local surface density, and $h=H/a$ the disk aspect ratio defined by disk scale height $H$. The local angular velocity is $\Omega =\sqrt{G{M}_{\star }/a^3}$. Suppose the torque on the planet iswhere the dimensionless coefficient $C$ can have either sign and depends on the disk model.

The planet's orbital angular momentum is $L=M_p{a}^2\Omega =M_p\sqrt{G{M}_{\star }a}$. Its derivative gives $\dot{L}=(M_{p}a\Omega /2)\dot{a}=\Gamma$, hencePositive $\Gamma$ in this convention produces outward planetary migration; negative $\Gamma$ produces inward planetary migration. The torque exerted on the disk has the opposite sign. The radial evolution timescale is half $L/|\Gamma |$, because $L\propto \sqrt{a}$. This expression assumes the low-mass Type I migration regime and the stated torque law, rather than a universal direction or rate.

A location where the net disk torque on a planet vanishes and, for a stable trap, restoring torque changes sign across the location. A torque zero alone does not guarantee stable trapping as the disk or planet evolves.

The orbital torque associated with disk material near the planet's corotation radius. Vortensity and entropy gradients matter because gas exchanged across the orbit carries different conserved properties. Horseshoe turns can make the response nonlinear even while the planet remains in the low-mass migration regime.

Repeated mixing or libration in an isolated corotation region can erase the gradient that drives its torque. Viscosity and thermal diffusion can help restore relevant gradients; excessive diffusion also changes the response. Saturation concerns the sustained torque, not cessation of all disk waves.

An entropy contrast carried by gas executing a horseshoe turn creates a density asymmetry around an embedded planet. The resulting torque can oppose the usual negative wave torque. Its persistence depends on thermal transport and the replenishment of the background entropy gradient.

A disk-gas trajectory that turns from one side of a planet's orbital radius to the other in the corotating frame. Its angular-momentum exchange contributes to the corotation torque; repeated libration without restoring gradients can saturate that contribution.

A rotating gravitational perturbation drives an epicyclic disk response where $m(\Omega -{\Omega }_p)=\pm \kappa$. The launched density waves carry angular momentum. For a planet, inner and outer waves exert torques of opposite signs, and their imbalance contributes to migration.

The $k_R=+1$ orbital resonance of an $m$-fold perturbation of a galactic disk. For near-circular stellar orbits, ${\Omega }_R$ is the radial epicyclic frequency. The harmonic changes radial action and azimuthal angular momentum in the ratio $\Delta J_R/\Delta J_{\phi }=1/m$, so an orbit gaining angular momentum gains random radial motion.

The $k_R=-1$ orbital resonance of an $m$-fold perturbation of a galactic disk. For near-circular stellar orbits, ${\Omega }_R$ is the radial epicyclic frequency. The harmonic changes radial action and azimuthal angular momentum in the ratio $\Delta J_R/\Delta J_{\phi }=-1/m$, so an orbit losing angular momentum gains random radial motion.

A planet transfers angular momentum to disk material through its spiral wakes, tending to push gas away from its orbital annulus. A persistent deep gap requires this clearing to compete with pressure and disk transport. A gap can reduce local supply without completely shutting off gas streams or fixing the final planet mass.

For a fixed mass on a Keplerian circular orbit, $L_s=M_s\sqrt{G{M}_{\ast }{r}_0}$ has derivative $M_s{r}_0\Omega /2$. The torque on the satellite is minus the torque $\Gamma$ on the disk. This gives the displayed drift law. In the density-slope satellite torque model, ${\dot{r}}_0=-(8\chi /9)\beta q({\Sigma }_0{r}_0^2/M_{\ast })r_0\Omega (H/r_0{)}^{-2}$.

The stellar obliquity is the angle between a star's spin and a planet's orbital angular momentum. The Rossiter-McLaughlin effect constrains the angle $\lambda$ between the sky projections of those two axes, rather than a component of $\psi$. If the orbital and stellar-spin inclinations to the line of sight are $i_{\mathrm{orb}}$ and $i_{\ast }$, then $\cos \psi =\cos i_{\ast }\cos i_{\mathrm{orb}}+\sin i_{\ast }\sin i_{\mathrm{orb}}\cos \lambda$; $\lambda$ alone generally does not determine $\psi$. Large misalignments constrain formation and planetary migration, but can also reflect a tilted birth disc.

A transiting body hides different velocity regions of a rotating stellar disc and causes an anomalous measured Doppler shift. Its time dependence constrains the sky-projected spin-orbit angle, which is not generally the full stellar obliquity.

Tidal locking settles a body's spin into a sustained relation with its orbital motion through tidal forces and tidal dissipation. The simplest case is synchronous rotation, a $1:1$ spin-orbit resonance. On an ideal circular orbit, with spin axis perpendicular to the orbital plane, the body keeps the same face toward its companion, as the Moon approximately does toward Earth. ^[ref].

This article retains the broader sense of tidal locking as capture into a spin-orbit resonance; the term is also commonly used specifically for synchronous rotation. Tidally evolved rotation need not be $1:1$: Mercury has a $3:2$ spin-orbit resonance with the Sun, turning three times in two orbital periods. The distinction matters when inferring illumination or climate from a spin period. ^[ref].

Tidal deformation is the change in a body's shape caused by differential gravity across it, that is, by a tidal force. The familiar two-bulge picture represents an elongated response along the companion direction, rather than a separate force pulling only the nearest surface. If spin and orbital motion differ, the forcing sweeps through the body's material. With tidal dissipation, the response is out of alignment with the instantaneous forcing, allowing a net torque. ^[ref].

A secular tidal torque is the torque remaining after averaging the changing tidal interaction over fast orbital and rotational motion. The misaligned tidal deformation lets the companion's gravity change the spin angular momentum; the orbit receives the opposite torque. Tidal dissipation converts mechanical energy into heat, while total angular momentum is conserved in an isolated system. Thus the spin-orbit coupling exchanges angular momentum rather than destroying it. ^[ref].

For intuition near synchronous rotation on a circular orbit, hold mean motion $n$ and moment of inertia $I$ fixed and adopt the illustrative linear dissipative lawHere $\omega$ is spin angular velocity, $\theta$ is a marked surface direction in an inertial frame, and $K$ is a response coefficient, not a universal material constant. A faster spin is slowed and a slower prograde spin is accelerated. The solution and the associated mechanical power loss areThe orbital contribution to this power is $-n⟨\tau ⟩$, and the spin contribution is $\omega ⟨\tau ⟩$; their sum becomes heat. Holding $n$ fixed is a local approximation to the coupled evolution. Actual secular tidal torque depends on the deformation response and orbital eccentricity, so this simple law does not establish which spin-orbit resonance an eccentric body captures into.

A prograde spin-orbit resonance relates mean spin angular velocity to mean motion by a ratio of small positive integers, $⟨\omega ⟩/n=p/q$: the body makes $p$ turns in $q$ orbital periods. Synchronous rotation is $1:1$; Mercury's $3:2$ spin-orbit resonance is nonsynchronous. The relevant combination of spin and orbital angles must stay bounded around a preferred orientation, not merely have a near-integer frequency ratio. For an elongated body on an eccentric orbit, the companion's gravitational torque on the permanent shape can sustain the resonance while tidal dissipation permits capture. ^[ref].

Synchronous rotation means one prograde turn about the body's own axis per orbital period, measured against fixed directions in an inertial frame. This $1:1$ spin-orbit resonance approximately describes the Moon's synchronous rotation toward Earth. ^[ref].

The striking geometric idea is that one turn in space can keep one face toward the companion: a fixed orientation in the frame turning with the orbit is a changing orientation in an inertial frame. This distinction helps us predict illumination: a synchronously rotating exoplanet orbiting its illuminating star can have persistent day and night hemispheres, whereas the Moon, synchronized to Earth, receives changing illumination from a different body, the Sun.

View a circular orbit along its perpendicular spin axis. Let $f$ be the orbital angle measured from a fixed direction, choosing an arbitrary zero for the circular limit of true anomaly, and let $\theta$ be the direction of an outward arrow painted on the body, measured from the same fixed direction. The direction from the body to its companion is $f+\pi$. Thus an inward-facing arrow satisfiesMore generally synchronous rotation keeps $\theta -f$ constant; its particular value selects which surface direction faces the companion. Over one orbital period $T$, both $f$ and the unwrapped $\theta$ increase by $2\pi$. This is the short geometric explanation for one turn per orbit. For a nonspinning body, $\dot{\theta }=0$ instead, and the companion direction sweeps through all surface longitudes. On an eccentric orbit, mean $1:1$ spin does not force $\theta -f$ to stay exactly constant: true anomaly advances nonuniformly, producing the longitudinal component of libration.

Facing a companion and facing the source of light are different relations. The Moon keeps its near side of the Moon approximately toward Earth, while the direction to the distant Sun is nearly fixed during the idealized orbit shown below. A fixed surface patch therefore turns toward, across and away from sunlight. The far side of the Moon has day and night too; its name does not mean permanent darkness. ^[ref].

If the Sun direction has fixed angle $\alpha$ and the marked surface patch has outward normal vector at angle $\theta$, its direct illumination is proportional to $\max (0,\cos (\theta -\alpha ))$. With synchronous rotation, $\theta =f+\pi$, so the patch goes through day and night even though it always faces Earth. In reality the Earth-Moon system advances around the Sun: the lunar spin and orbital cycle relative to distant stars is about $27.3$ Earth days, while the lunar phase cycle has a synodic period of about $29.5$ Earth days. These numbers describe the actual system, not the fixed-sunlight cartoon. ^[ref].

For a synchronously rotating exoplanet whose orbital companion is itself the illuminating star, the direction to the star turns with the orbit. On an ideal circular orbit, with perpendicular spin axis and no additional oscillations, this direction is fixed relative to the surface, giving persistent day and night hemispheres. Day-night heat redistribution then helps determine the nightside temperature; permanent lack of direct starlight does not imply zero temperature. Orbital eccentricity and departures from this ideal spin geometry modify that simple picture.

Tidal heating converts mechanical energy of tidal deformation into internal heat through tidal dissipation. Sustained orbital eccentricity or a spin axis tilted relative to the orbital plane's perpendicular can maintain changing tidal deformation and heat an exoplanet interior. Thus synchronous rotation alone does not eliminate tidal heating.

Tidal dissipation converts mechanical energy associated with a time-varying tidal deformation into heat. It can reduce orbital eccentricity, drive synchronous rotation, and contribute deep planetary heating. A circular orbit with exact synchronous rotation, perpendicular spin axis and no additional forcing has a static tidal deformation in the body's frame and need not sustain substantial tidal heating.

An astronomical conjunction occurs when two orbiting bodies have the same apparent longitude as viewed from their primary or from an observer.

For two uniform angular motions with $P_b<P_c$, an inner exoplanet transit at zero and next outer exoplanet transit at $\tau$, the next conjunction solves $(n_b-n_c)t+n_c\tau =2\pi$. Hence its longitude modulo $2\pi$ isFor an eccentric outer body this is only an approximation. Let $f_L,M_L$ be its line-of-sight anomalies and ${\Delta }_e(t)=f(t)-M(t)-(f_L-M_L)$. Exact true alignment instead requires $(n_b-n_c)t+n_c\tau -{\Delta }_e(t)=2\pi$. The correction follows directly from uniform mean anomaly and nonuniform true anomaly.

Keplerian shear is differential orbital motion caused by $n\propto a^{-3/2}$: neighboring material at smaller semi-major axis advances in longitude relative to material farther out.

Phase mixing disperses an initially localized ensemble because its members have slightly different orbital frequencies, causing their phases to spread with time.

A planetesimal is a solid body participating in the growth and collisional evolution of a planetary system. Its collisions can build larger bodies or fragment existing ones; the result depends on relative velocity, material strength and self-gravity. Gravitational encounters can increase random motion, while gas drag and dynamical friction can reduce it.

The category marks a physical role rather than one sharply defined diameter. A population can feed planetary embryos during planet formation or act as the parent bodies of a later debris disk. The same measured total mass can have very different collision behaviour when distributed among many small bodies or a few large ones.

For an impactor number density per size $n(D_{\mathrm{im}})$, relative speed $v_{\mathrm{rel}}$, and collision cross-section $\sigma (D,D_{\mathrm{im}})$, the differential collision rate of one target is $n(D_{\mathrm{im}})v_{\mathrm{rel}}\sigma (D,D_{\mathrm{im}})d{D}_{\mathrm{im}}$.

The specific impact energy is an impactor's kinetic energy divided by the target mass. For equal-density spherical bodies,

A binary planetesimal is a gravitationally bound pair of planetesimals. Its internal two-body problem and its centre of mass orbit respond differently to external forces, including gas drag.

For linear gas drag with aerodynamic stopping times $t_1,t_2$, masses $m_1,m_2$, total mass $m=m_1+m_2$, and relative velocity $\mathbf{v}$, defineIf $\mathbf{w}$ is the gas velocity relative to the centre of mass, the relative drag acceleration is $-\kappa \mathbf{v}+\Delta \mathbf{w}$. The first term dissipates internal specific orbital energy; the second is a differential headwind that vanishes for equal aerodynamic stopping times.

For a weakly perturbed circular binary planetesimal, the constant differential headwind does zero work averaged over the unperturbed orbit. Orbit-averaged drag work givesWith constant aerodynamic stopping times, $a(t)=a_0{e}^{-t/\tau }$. This averaging assumes drag is slow compared with the orbital period and does not appreciably distort the circular orbit within one period.

For a nearly circular stellar orbit of radius $a_b$ in gas moving at $(1-\eta )$ times the local circular speed, the averaged tangential gas drag on the centre of mass is $-\eta {\kappa }_b{v}_K$. Orbit-averaged drag work givesIf $m_2/m_1=ϵ\ll 1$ and aerodynamic stopping time scales as $m^k$, $0<k<1$, then $\tau /{\tau }_b\simeq \eta {ϵ}^k\ll 1$. Internal contraction can outpace stellar migration. Tightly gas-coupled drift requires solving radial and azimuthal motion together instead of this weak-drag approximation.

For a small body with orbital elements $a,e,I$ relative to a planet on a circular orbit of semi-major axis $a_p$, the Tisserand parameter is approximatelyIt is derived from the Jacobi constant and is approximately conserved across separated weak encounters in the circular restricted three-body problem.

For bound planet-crossing Kepler orbits, $T=2/(q+Q)+2\sqrt{2qQ/(q+Q)}\cos I$, with the planet radius taken as one. At fixed $q<1$, $T$ is maximized by $I=0$ and $Q=1$. For $2<T<3$ this gives $q\ge (1-\sqrt{3-T}{)}^2/[1+2\sqrt{3-T}-(3-T)]$. This form remains regular at $T=\sqrt{8}$. It is an invariant-based accessibility bound, not a guarantee that a scattering history attains it; for $0<T\le 2$ it gives no positive universal pericentre floor.

In units where a circular planet has orbital radius and speed one, the particle's relative speed $u$ outside a close encounter satisfies $u^2=3-T$, to leading order in the planet-to-star mass ratio. This follows from the vis-viva equation and $u^2=v^2+1-2h_z$ at the encounter radius. A short conservative encounter changes the direction of this relative velocity while preserving its magnitude.

The disturbing function is the perturbing part of the gravitational potential used to derive the evolution of orbital elements in celestial mechanics.

For a particle far inside a planet's stellar orbit, the planet-centred disturbing function is $G{M}_{\star }/|{\mathbf{r}}_{\star }-\mathbf{r}|-G{M}_{\star }\mathbf{r}\cdot {\mathbf{r}}_{\star }/r_{\star }^3$. Expanding the first term cancels its dipole against the indirect disturbing function. Dropping the spatially constant monopole leaves $G{M}_{\star }{r}^2(3{\cos }^2\psi -1)/(2r_{\star }^3)$, where $\psi$ is the angle between the two offset vectors. The omitted instantaneous terms are $O(G{M}_{\star }{r}^3/r_{\star }^4)$. The gradient of the retained expression is the leading differential stellar tidal acceleration.

For two circular Kepler orbits with mutual inclination $I_m$, independently average the stellar quadrupole disturbing function over both mean anomalies. A circular ring with unit normal $\mathbf{k}$ has $⟨{\hat{r}}_i{\hat{r}}_j⟩=({\delta }_{ij}-k_i{k}_j)/2$. Contracting the two ring tensors gives $⟨{\cos }^2\psi ⟩=[1+(\mathbf{k}\cdot {\mathbf{k}}_p{)}^2]/4$, proving the displayed exact angular factor. For small orbital inclinations relative to a common plane, $I_m^2=I^2+I_p^2-2I{I}_p\cos ({\Omega }_p-\Omega )+O(I^4)$. Thus the secular angular factor is $1/4-3(I^2+I_p^2)/8+3I{I}_p\cos ({\Omega }_p-\Omega )/4$. This potential generates the slow stellar contribution to satellite nodal precession.

The Lagrange planetary equations express changes of osculating orbital elements in terms of derivatives of a disturbing function. Their small-eccentricity secular limit becomes a linear equation for complex eccentricity.

If retained Lagrange planetary equations have $\dot{a}=-2(p+q)a{e}^q{C}_r\sin \varphi$ and $\dot{e}=-q{e}^{q-1}{C}_r\sin \varphi$, then $d(e^2)/d\log a=q/(p+q)$. Thus $e^2-[q/(p+q)]\log (a/a_{\mathrm{ref}})$ is constant. A resonantly oscillating orbital eccentricity–semi-major axis trajectory retraces one curve rather than enclosing area in that projection.

For disturbing function $\mathcal{R}=(G{M}_2/a_2)f(\alpha )e\cos \varphi$, $\alpha =a/a_2$, the leading planar Lagrange planetary equations give $\dot{\lambda }=n+n(M_2/M_{\star })eg(\alpha )\cos \varphi$. The derivative with respect to $a$ holds the other osculating orbital elements fixed. If one writes $\lambda =n(t)t+ϵ(t)$, its literal derivative includes $t\dot{n}$; the physical mean-longitude equation cannot be obtained by silently dropping that term. A consistent accumulated-phase parametrization uses $\lambda ={\int }_0^{t}n(s)ds+{ϵ}_{\mathrm{int}}$.

For perturbing acceleration $R\hat{\mathbf{r}}+T\hat{\mathbf{\theta }}$ on an osculating orbital element system with fixed central parameter $\mu$,To prove it, differentiate the specific orbital energy $-\mu /(2a)$, obtaining $\dot{a}=(2a^2/\mu )(R\dot{r}+Th/r)$. Substitute $\dot{r}=\mu e\sin f/h$, $h/r=\mu (1+e\cos f)/h$ and $h^2=\mu a(1-e^2)$. Here $f$ is the true anomaly; a normal perturbing acceleration does no instantaneous work.

The direct disturbing function is the gravitational potential of the perturbing body evaluated at the disturbed body. In heliocentric coordinates it is proportional to the inverse separation of the two bodies.

The indirect disturbing function accounts for the acceleration of a non-inertial origin, such as a star-centered coordinate frame accelerated by a perturbing planet.

The D'Alembert characteristic constrains each disturbing-function harmonic so that the integer coefficients of its angular arguments sum to zero, with corresponding parity and eccentricity-inclination order restrictions imposed by rotational symmetry.

A Laplace coefficient is the Fourier coefficientIt packages the radial dependence of a disturbing-function expansion.

A short-period orbital perturbation is driven by rapidly varying orbital-angle combinations and usually averages to zero over many conjunctions away from resonance.

A resonant orbital perturbation has a slowly varying angular argument because a nearly integer relation between orbital frequencies prevents successive conjunction effects from averaging away.

A secular perturbation survives averaging over fast orbital phases and drives slow precession and exchange of eccentricity or inclination.

For a quantity on a periodic Kepler orbit, its orbit average is $⟨x⟩=P^{-1}{\int }_0^{P}x(t)dt$. Equal increments of true anomaly are not equal increments of time. Eccentric anomaly gives the weight $dt=(1-e\cos E)dE/n$, which yields $⟨\cos f⟩=-e$. Slow orbital perturbations can be averaged using elements held fixed over a single orbital period.

The difference between the angular momentum of circular coplanar orbits at fixed semi-major axes and the system's actual angular momentum along the reference normal. With circular momenta $L_j=m_j\sqrt{G{M}_{\star }{a}_j}$ it is ${\sum }_j{L}_j(1-\sqrt{1-e_j^2}\cos I_j)$. For small planar orbital eccentricities its quadratic part is $\frac{1}{2}{\sum }_j{L}_j|z_j{|}^2$. Its conservation in $\dot{z}=iAz$ requires $L_j{A}_{jk}=L_k{A}_{kj}$.

The proper eccentricity is the magnitude of the homogeneous part of the complex eccentricity in linear secular perturbation theory. It remains constant in conservative theory, while its proper longitude of periapsis advances at the free precession frequency.

The proper longitude of periapsis is the phase of the free, or proper, component of a body's complex eccentricity. In a linear secular model with constant forced eccentricity vector $z_f$,The proper eccentricity $e_p$ is the radius of the corresponding circle in an Argand diagram, and ${\varpi }_p$ specifies position around that circle.

It generally differs from the actual longitude of periapsis $\varpi =\arg z$. Random proper phases therefore do not automatically mean uniformly random actual apsidal orientations around the star. This distinction matters when constructing a proper-eccentricity annulus or inferring the width of an eccentric debris ring. If the proper component is zero, its phase is undefined.

Laplace-Lagrange secular theory linearizes the orbit-averaged disturbing function to second order in eccentricities and inclinations. It represents planetary eccentricity evolution by a constant matrix.

For a test particle interior to a distant eccentric planet, write $\alpha =a/a_p\ll 1$. The leading Laplace coefficients are $b_{3/2}^1=3\alpha$ and $b_{3/2}^2=15{\alpha }^2/4$. Thus its forced eccentricity is $Z_f=5\alpha e_p/4$ and its secular frequency is $A=(3/4)n_p(m_p/M){\alpha }^{3/2}$. More distant inner particles have larger forced orbital eccentricities and faster secular precession, producing a growing gradient of complex eccentricity.

Secular forcing of a test particle is the slow change of its orbital eccentricity or orbital inclination caused by a perturber after averaging over fast orbital phases. In a simple planar Laplace-Lagrange secular theory model, a fixed external eccentricity vector supplies a constant forced eccentricity $z_f$ andHere $z$ is complex eccentricity and $g$ the free precession rate. The solution is a circle centred on $z_f$ in an Argand diagram.

An initially circular orbit has $z(0)=0$ and consequently reaches eccentricity $2|z_f|$ in this idealized model. A changing perturber, multiple forcing frequencies or a secular resonance changes the response. Secular forcing by a distant outer planet supplies one concrete approximation to the forcing and precession coefficients.

If both precession and forcing in $\dot{z}=ig(t)(z-z_f)$ share a common mass factor while $z_f$ stays fixed, then $z=z_f+(z(0)-z_f)\exp (i{\int }_0^{t}g(t\prime )dt\prime )$. The Argand diagram circle is unchanged; only its traversal speed varies. With $g(t)=g_{f}t/t_p$, the phase is $g_f{t}^2/(2t_p)$.

The eigenvalues of the Laplace-Lagrange secular matrix are apsidal precession frequencies, and its eigenvectors give the relative eccentricities and apsidal phases of the corresponding normal modes.

The complex eccentricity packages eccentricity and longitude of periapsis as $z=e{e}^{i\varpi }$, turning linear secular dynamics into a complex linear system.

Choose the fixed planetary pericentre as the real axis, with positive forced amplitude $Z_f$. Then the linear secular equation $\dot{z}=iA(z-Z_f)$ has solution $z=Z_f+(z(0)-Z_f)e^{iAt}$. An initially circular orbit follows $z=Z_f(1-e^{iAt})$, a circle in the complex plane centered at the forced eccentricity $Z_f$, with free eccentricity amplitude $Z_f$. Its maximum orbital eccentricity is $2Z_f$ and its phase-averaged mean orbital eccentricity is $4Z_f/\pi$. The traversal is counterclockwise for $A>0$ and clockwise for $A<0$; at $A=0$ the initial state is stationary. The phase averages assume a nonzero secular frequency.

A secular eigenmode is a normal mode in which all complex eccentricities precess at one eigenfrequency with fixed amplitude ratios and relative apsidal phases.

For a real two-planet Laplace-Lagrange secular matrix with negative off-diagonal entries, the upper secular eigenmode has anti-aligned apsides and the lower has aligned apsides. Their positive planet-1 to planet-2 amplitude ratios obey $r_p{r}_m=A_{12}/A_{21}$. Under the conventional angular momentum deficit normalization this is $L_2/L_1$.

Complex inclination packages small orbital inclination $I$ and longitude of ascending node $\Omega$. It differs from complex eccentricity. Linear Laplace-Lagrange secular theory gives $\dot{\mathbf{y}}=iB\mathbf{y}$, whose matrix exponential solution decomposes into nodal secular eigenmodes.

The particular, forced part of a complex inclination solution. For a fixed, circular, inclined inner planet and an exterior test ring, linear Laplace-Lagrange secular theory gives $\dot{w}=-i\nu (w-w_p)$, so $w_f=w_p$ and $w=w_p+(w(0)-w_p)e^{-i\nu t}$. The ring describes a circle about the planetary tilt in an Argand diagram; relative to the planet, its mutual inclination stays constant while its ascending node precesses.

For circular planetary orbits with angular momenta $L_j=M_j\sqrt{G{M}_{\star }{a}_j}$, the inclination matrix has positive off-diagonal entries, zero row sums, and $L_j{B}_{jk}=L_k{B}_{kj}$. Positive diagonal symmetrization of a matrix with $W=\mathrm{diag}(L_j)$ makes $W^{1/2}B{W}^{-1/2}$ a symmetric matrix. Moreover,All eigenvalues are real and nonpositive. The zero eigenvalue is a common tilt of every orbital plane; nonzero modes describe relative nodal precession.

If two nearby inner planets couple much more strongly to each other than to a distant third planet, their slow secular eigenmode has a common complex inclination. With $L_{12}=L_1+L_2$, the effective inner-to-outer coupling isAngular momentum balance makes the reverse coupling $c{L}_{12}/L_3$. The reduced matrix has eigenvalues $0$ and $-c(1+L_{12}/L_3)$. If $L_3\gg L_{12}$, the nonzero frequency is approximately $-c$, while the fast inner differential mode has frequency $-(B_{12}+B_{21})$.