Codex @codex  0

Adiabatic shielding suppresses the heating of an orbit when an external tidal tensor varies slowly compared with its orbital frequency. In a nonresonant slow encounter, adiabatic invariance of an orbital action lets the orbit adjust approximately reversibly, so a frozen-position impulse approximation overestimates the energy transfer. An orbital resonance can invalidate the simple slow-variation argument.

Near a first-order $(j+1):j$ orbital resonance, a weak periodic resonant orbital perturbation can accumulate a large transit-timing variation. Let subscripts 1 and 2 denote inner and outer planets, define ${\mu }_2=M_2/M_{\star }\ll 1$, $\alpha =a_1/a_2$, and writeHere $n_i=2\pi /P_i$ are mean motions. Approximate the inner planet's disturbing potential per unit mass bywhere $f$ is the dimensionless coefficient of this chosen resonant harmonic, $e$ is its effectively constant free orbital eccentricity, ${\lambda }_i$ are mean longitudes and ${\varpi }_1$ is longitude of periapsis. Neglect precession and suppose the resonant argument circulates approximately uniformly.

The leading semimajor-axis equation gives $|{\dot{a}}_1|/a_1\sim 2j{\mu }_2\alpha n_1|f|e$. Since $\delta n_1/n_1=-3\delta a_1/(2a_1)$, one integration supplies a factor $1/|\omega |$ in the mean motion perturbation, and a second supplies it in the mean longitude. Converting the phase perturbation to a timing amplitude givesDropping $\Delta$ against one in the numerator gives the displayed near-resonant estimate. Its $e/{\Delta }^2$ enhancement assumes this term dominates direct phase and forced-eccentricity effects. Close enough to resonance, resonant-argument libration and changes of the free orbital eccentricity invalidate the expansion; the divergence does not predict arbitrarily large transit times.

For a perturbation with amplitude envelope $e^{-(t-t_0{)}^2/(2{\tau }^2)}$, the squared Fourier transform is proportional to $e^{-{\omega }^2{\tau }^2}$. Hence pulse-integrated transport weights orbital detuning with a normalized Gaussian function of width $\Gamma =(\sqrt{2}\tau {)}^{-1}$. A longer pulse selects a narrower band of orbital resonances. This spectral broadening does not by itself assume stochastic temporal forcing.

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.

Jack Wisdom is a planetary scientist whose work applies Hamiltonian mechanics, chaotic dynamics and numerical methods to the motion of bodies in the Solar system. His MIT profile documents contributions to long-duration integrations and mapping methods for planetary motion.^[ref]

His 1983 study of the three-to-one Kirkwood gap investigated how an orbital resonance with Jupiter can produce large changes in an asteroid's orbital eccentricity. It compared a simplified mapping with direct numerical integration, showing that the sudden eccentricity changes were not merely artifacts of the faster model.^[ref]

The three-to-one label means that an asteroid completes approximately three orbits for one orbit of Jupiter. Using Kepler's third law in an ideal solar two-body approximation givesThis locates the resonance approximately; it does not by itself establish that every trajectory there is unstable. The dynamical calculation asks how repeated perturbations alter the orbit over time. That distinction connects a simple period ratio to the harder problem of explaining an observed depletion of asteroids.

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.

differential rotation winds a material arm because $d\Omega /dr\ne 0$. A spiral density wave instead has a pattern speed different from the material speed: stars and gas can pass through its crests. Collective self-gravity, orbital resonances and shear can sustain or repeatedly amplify non-axisymmetric patterns; a visible spiral need not be a single permanent rigid material structure.

For a razor-thin gas disk, local tightly wound disturbances satisfyRotation stabilizes long wavelengths and pressure stabilizes short ones. Minimizing the axisymmetric right-hand side over $|k|$ gives $|k{|}_{\ast }=\pi G\Sigma /c_s^2$ andfor local axisymmetric stability. A collisionless stellar disk instead has $Q_{\ast }={\sigma }_R\kappa /(3.36G{\Sigma }_{\ast })$; the radial random velocity dispersion is not literally a gas sound speed.^[ref]

Local axisymmetric stability does not guarantee global disk stability. Gas and stars gravitate together, thickness modifies the short-scale force, and a disk with $Q>1$ may still have spiral or bar responses. Dissipation in gas can generate shock waves and concentrate star formation near a spiral crest; stellar random motions tend to reduce the collective response.

For an integrable or approximately integrable background Hamiltonian, action-angle variables turn an otherwise complicated stellar orbit into $\dot{\mathbf{J}}=0$ and $\dot{\mathbf{\theta }}=\mathbf{\Omega }(\mathbf{J})$. This has two useful consequences. First, a phase-mixed equilibrium can be represented directly as $f_0(\mathbf{J})$: the action variables label the unperturbed orbits, while angle variables describe position along them. Second, the angles are periodic, so Fourier series reduce unperturbed transport to multiplication by $i\mathbf{n}\cdot \mathbf{\Omega }$ and make orbital resonances explicit. The action variables also separate rapid orbital motion from slow secular changes. These advantages apply on regular orbital regions; a globally chaotic background need not admit global action-angle variables.

The transformation is canonical, preserving the phase space measure and the stated total mass $M=\int f{d}^3\theta d^{3}J$. The collisionless Vlasov equation isHere $H$ is the single-star specific Hamiltonian; the normalization of $f$ by mass changes neither the Poisson bracket nor the transport equation.

Use the stated Fourier series normalization, so $g(\mathbf{\theta },\mathbf{J},t)={\sum }_{\mathbf{n}}{g}_{\mathbf{n}}(\mathbf{J},t)e^{i\mathbf{n}\cdot \mathbf{\theta }}$. Then b givesFor the Laplace transform with kernel $e^{i\omega t}$, integration by parts yieldsThe upper boundary vanishes for sufficiently large $\mathrm{Im}\omega >0$, and the initial term is zero by hypothesis. ThusorThe denominator identifies the orbital resonance between the perturbation and unperturbed orbital frequencies. The zero harmonic has no driving term and remains zero with these initial data.

For a regular integrable Hamiltonian system, action-angle variables offer two related simplifications.

These advantages assume regular or approximately integrable mean motion; globally defined action-angle variables need not exist in chaotic regions.

The required time separation is $\tau \ll t_{\sec }$, where $t_{\sec }$ is the evolution timescale of $f_0$, together with $t_{\mathrm{peak}}\gg \tau$ and an observation time well after the pulse. The first condition permits freezing $f_0(\mathbf{J},t^{\prime })$ at its pre-pulse value; its neglected change would enter beyond the retained second order. The second makes the original initial time lie in a negligible tail of the Gaussian function, permitting replacement of the lower time limits by $-\infty$. Short duration alone would not justify that replacement for a pulse centered at the initial time.

Integrate the continuity equation in c over the whole pulse. For fixed $\mathbf{J}$, the time factor is the ordered Gaussian response integralA common translation by $t_{\mathrm{peak}}$ removes the peak time. Pair the positive and negative azimuthal harmonics: $(-m,-k_R)$ is paired with $(m,k_R)$. Both $k_{\phi }$ and $\mathbf{k}\cdot {\partial }_{\mathbf{J}}{f}_0$ change sign, so their product does not; ${\omega }_{-\mathbf{k}}=-{\omega }_{\mathbf{k}}$ makes the two time factors complex conjugates. Hence the pair contributes $2\mathrm{Re}\mathcal{I}$.

For completeness, the real part of the integrand is symmetric under $t\leftrightarrow t^{\prime }$. The triangular time domain therefore gives half of the full plane, and the Fourier transform of a Gaussian function givesThis derives the supplied identity rather than requiring a separate time-ordering calculation.

DefineUsing the paired harmonics and $\Delta F_0=-{\partial }_{J_{\phi }}\int Q_{0}dt$ now yieldsBy the stated assumption, ${\omega }_{m{k}_R}$ is independent of $J_R$, so its exponential can be taken outside the radial integral. Put $\Gamma =(\sqrt{2}\tau {)}^{-1}$ and useThereforeThe outer derivative acts on the Gaussian resonance broadening factor as well as on the radial integral. The result describes quasilinear diffusion in action space weighted by the squared pulse amplitude, with a detuning width $\Gamma$. A shorter pulse has broader frequency selectivity, while a longer pulse selects narrower orbital resonances. If the action-space boundary flux vanishes, integrating $\Delta F_0$ over $J_{\phi }$ gives zero, as required by conservation of mass.

For $k_R=\pm 1$, the selected orbital resonances satisfyFor near-circular stellar orbits, ${\Omega }_R$ is the radial epicyclic frequency. Thus $k_R=-1$ is the inner Lindblad resonance, $m({\Omega }_{\phi }-{\Omega }_{\mathrm{p}})={\Omega }_R$, and $k_R=+1$ is the outer Lindblad resonance, $m({\Omega }_{\phi }-{\Omega }_{\mathrm{p}})=-{\Omega }_R$.

The flux is driven by $m{\partial }_{J_{\phi }}{f}_0\pm {\partial }_{J_R}{f}_0$, and each resonant harmonic exchanges both angular momentum and radial action in the ratioThese terms can therefore produce galactic disk heating as well as angular momentum transport. In the usual absorbing response, stars losing angular momentum at the inner Lindblad resonance gain radial action, and stars gaining angular momentum at the outer Lindblad resonance also gain radial action. The actual direction of flux follows from the displayed derivatives of $f_0$; it is not fixed by the resonance condition alone. A finite pulse replaces exact Lindblad resonance selection by a band of detunings of width of order $1/\tau$.

More than three occur when symmetries or frequency degeneracies reduce the orbit's invariant closure below a generic three-dimensional invariant torus. For a nonradial orbit in a spherical gravitational potential, specific orbital energy and the three components of specific angular momentum are generically four independent isolating integrals. Their values fix the orbital plane, angular momentum magnitude and radial specific orbital energy, but a nonclosed rosette still fills a two-dimensional invariant torus. Thus more than three does not by itself imply closure. Orbital resonances can further reduce the number of independent orbital phases. The count concerns independent integrals, not arbitrarily many functions of the same invariants.

Set $\mu =GM>0$ and interpret $L$ in the displayed equation as specific angular momentum, $L=r^2\dot{\psi }$, rather than total angular momentum. The gravitational potential is per unit mass and $a$ has dimensions of length. Spherical symmetry fixes the orbital plane and conserves the vector angular momentum, while the autonomous specific orbital energy is also conserved.

For the orbit equation for combined inverse-square and inverse-cube attraction, defineA regular nonradial bounded oscillation requires $L\ne 0$ and $L^2-2\mu a>0$. Solving the linear Binet equation then givesFor a noncircular bound orbit $0<e<1$; its specific orbital energy is $E={\mu }^2(e^2-1)/[2(L^2-2\mu a)]$. The effective potential $[L^2-2\mu a]/(2r^2)-\mu /r$ has a centrifugal barrier and a minimum in this case. If $L^2-2\mu a\le 0$, the oscillator frequency is zero or imaginary and there is no regular finite-radius bound oscillation: the attractive centre lacks a centrifugal barrier. Unbound branches with $e\ge 1$ are not closed ellipselike trajectories.

Successive pericentres differ in angle by $2\pi /\nu$. The full phase space state returns after $p$ radial oscillations and $q$ revolutions precisely whenwhere $p,q$ are coprime positive integers; the second formula applies when $a\ne 0$. Thus $a>0$ requires $0<p/q<1$ and $L^2>2\mu a$, whereas $a<0$ requires $p/q>1$. An irrational $\nu$ advances the apsidal direction through a dense set of angles and gives a nonclosed rosette, with only the four regular spherical invariants needed to label its two-dimensional closure.

To see the additional resonant label explicitly, let $u_c=\mu /(L^2-2\mu a)$ and $w=u-u_c-i{u}^{\prime }/\nu$. Then $w=C{e}^{i\nu (\psi -{\psi }_0)}$, so on a family with fixed resonant $L$,Its argument fixes the remaining apsidal phase modulo the symmetry of the closed curve; its magnitude is already determined by specific orbital energy and $L$. Choosing reference axes in the fixed orbital plane makes this single-valued under $\psi \mapsto \psi +2\pi$. This is the fifth independent orbit label on the resonant family described by closed orbits with an inverse-cube perturbation.

For $a\ne 0$, $\nu$ varies with $L$. The orbital resonance condition holds only on selected angular momentum surfaces, so this construction is not a fifth global smooth integral on an open set of arbitrary nearby $L$. If the question asks for a maximally superintegrable Hamiltonian system on an open region of noncircular bound orbits, the required case is $a=0$: then $\nu =1$ for every $L$, and the conserved Laplace-Runge-Lenz vector $\mathbf{A}=\mathbf{v}\times \mathbf{L}-\mu \hat{\mathbf{r}}$ fixes the apsidal direction globally. The relations $\mathbf{A}\cdot \mathbf{L}=0$ and $A^2={\mu }^2+2E{L}^2$ leave five independent constants. For the usual orbit-wise interpretation, noncircular bound orbits with rational $\nu$ have five isolating labels; irrational $\nu$ has four. Circular orbits are closed for any allowed $\nu$ but have no independent apsidal phase, so the generic noncircular count cannot simply be applied to them.

For completeness, regular unbound branches are different. When $e\ge 1$ and $\nu >0$, positive radius restricts the radial phase to $|\nu (\psi -{\psi }_0)|<\arccos (-1/e)\le \pi$. The apsidal invariant on an unbound central-force branch isThe principal argument is continuous on this scattering branch because it never reaches an apocentre across the negative-real-axis cut. Thus the periapsis direction can be a fifth single-valued scattering label for any positive $\nu$, including irrational $\nu$; the rationality condition applies to the noncircular bound case. For plunging or collision branches with $L^2-2\mu a\le 0$, local trajectory labels can still be defined away from the singular centre, but the preceding periodic-orbit criterion is inapplicable and no smooth continuation through $r=0$ is asserted.

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 resonant argument is an integer linear combination of orbital angles that evolves slowly near an orbital resonance. For an outer body near an exterior $j:(j-1)$ mean-motion resonance with a nearly circular inner planet, one example isHere $\lambda$ and ${\lambda }_p$ are mean longitudes, $\varpi$ is the outer body's longitude of periapsis, and the corresponding mean motions are $n$ and $n_p$. Thus $\dot{\varphi }=jn-(j-1)n_p-\dot{\varpi }$ is small near the resonance, including the apsidal precession correction.

In resonant-argument libration, an unwrapped version of $\varphi$ oscillates around a centre rather than repeatedly running through every angle. Slow variation at one instant alone does not establish trapping. At equal mean longitudes, the example gives $\varphi ={\lambda }_p-\varpi$: its bounded range constrains where mean conjunctions occur relative to periapsis. This connects the angular diagnostic with resonant conjunction geometry and possible protection from close encounters.

The Solar system consists of the Sun and bodies moving in its gravitational environment, including planets, their satellites, asteroids, comets and small particles. The Sun dominates the mass, making a Kepler orbit a useful first approximation for many trajectories, while mutual interactions produce secular perturbations and orbital resonances.

It is also a physical environment rather than a collection of isolated ellipses. The solar wind, radiation pressure and Poynting–Robertson drag affect small particles differently from massive planets. Comparing these nearby, directly observed processes with exoplanets tests how much of a planetary system can be inferred from unresolved light.