Integrable Hamiltonian system

Liouville integrability means having $n$ almost-everywhere independent first integrals in mutual Poisson involution.

On the regular set of two uncoupled oscillators with irrational frequency ratio, every continuous globally defined integral of motion is constant on each fixed-action flat torus, because every orbit is dense there. A differentiable integral is consequently a function of the two actions, so its differential of a smooth map lies in their span. Local angle combinations can still give extra local integrals; they fail to be single-valued on the torus.

For a compact connected regular level of $n$ independent Poisson-commuting functions, the periods of their joint flow form a full-rank discrete subgroup $\Lambda \subset {\mathbb{R}}^n$. The level is ${\mathbb{R}}^n/\Lambda$. A local smooth basis of these periods determines action integrals and the normalized periodic fields used to construct action-angle variables.

A Stäckel potential allows additive separation of the Hamilton-Jacobi equation in an orthogonal coordinate system. In axisymmetric prolate spheroidal coordinates, the relative-potential form is $\psi =[G(\lambda )-G(\mu )]/(\lambda -\mu )$. It gives energy, axial angular momentum and an independent third separation integral. Here $G$ is a one-variable generating function, not the gravitational constant. Generic bounded orbits are regular invariant-torus motions; resonance and coordinate degeneracy require separate treatment.

For fixed $E,l,I$, separated motion obeys $p_{\tau }^2=[E\tau +G(\tau )-I+(\alpha -\beta )l^2/(2(\tau +\alpha ))]/[2(\tau +\alpha )(\tau +\beta )]$. Allowed intervals are where this is nonnegative, with turning points and coordinate-fold endpoints treated separately. A generic nonresonant bound orbit with $l\ne 0$ circulates azimuthally and densely samples a projected volume bounded by confocal spheroids and hyperboloids. It is not uniformly distributed in volume. Resonant or zero-angular-momentum orbits can project to lower-dimensional regions; integrability alone does not imply that every orbit fills a three-dimensional shell.

For $H=\frac{1}{2}\sum p_j^2/h_j^2-\psi$, $l=p_{\varphi }$ and a Stäckel potential, multiplying the separated Hamilton-Jacobi equation by $\lambda -\mu$ gives $2(\tau +\alpha )(\tau +\beta )p_{\tau }^2-(\alpha -\beta )l^2/[2(\tau +\alpha )]-G(\tau )-E\tau =-I$ for both $\tau =\lambda ,\mu$. Eliminating $E$ gives $I=\frac{1}{2}[\mu v_{\lambda }^2+\lambda v_{\mu }^2+(\lambda +\mu +\alpha )v_{\varphi }^2]-[\mu G(\lambda )-\lambda G(\mu )]/(\lambda -\mu )$. Physical components are $v_j=h_j{\dot{q}}_j$ and $p_j=h_j{v}_j$. Separation proves conservation locally; direct Poisson brackets extend it across regular coordinate patches.

First integrals are in involution when every pair has zero Poisson bracket. An $n$-degree-of-freedom Hamiltonian system is Liouville integrable when it has $n$ almost-everywhere independent first integrals in involution.

An invariant torus is a subset homeomorphic to a torus and preserved by a flow. In action-angle variables, an integrable Hamiltonian flow moves linearly in the angular coordinates of each regular invariant torus.

A quasiperiodic orbit on a torus has angular coordinates $\theta (t)=\theta (0)+\omega t$ modulo $2\pi$. It is periodic exactly when the nonzero frequencies are commensurate; rationally independent frequencies produce a dense orbit.

For two uncoupled harmonic oscillators with nonzero amplitudes and irrational frequency ratio ${\omega }_1/{\omega }_2$, the orbit is bounded but not periodic. It winds quasiperiodically on a two-dimensional invariant torus; a common period would force the frequency ratio to be rational.

The direct product of integrable Hamiltonian systems is integrable. Pull each family's commuting integrals of motion back to the product phase space. Integrals from different factors have zero Poisson bracket because they use disjoint canonical coordinates, and their differentials are independent wherever the differentials within each factor are independent.

Forthe functions $I_j=(q_j^2+p_j^2)/2$ are independent almost everywhere, Poisson commute, and sum to $H$. A regular common level $I_j=c_j>0$ is the product of circles of radii $\sqrt{2c_j}$, hence an invariant $n$-torus. Allowing some $c_j=0$ gives lower-dimensional singular tori.

A compact connected regular common level of $n$ commuting integrals is an $n$-torus with local action-angle coordinates.

For a one-degree-of-freedom closed orbit, action-angle coordinates make the Hamiltonian a function $H(I)$ and give $\dot{I}=0$, $\dot{\varphi }=dH/dI$.

Near a regular Lagrangian torus, choose a one-form $\lambda$ with $d\lambda =-\omega$ and a basis of fibre cycles ${\gamma }_j$. The action $I_j=(2\pi {)}^{-1}{\int }_{{\gamma }_j}\lambda$ has differential of a smooth map $d{I}_j=(2\pi {)}^{-1}{\sum }_i{t}_i^{(j)}d{F}_i$, where $t^{(j)}$ is the corresponding joint-flow period. Nonsingularity of the period matrix makes these actions local coordinates transverse to the tori.

With angle period $2\pi$, the action is the enclosed phase-space area divided by $2\pi$:

When a Hamiltonian parameter varies slowly compared with the orbital period and no separatrix is crossed, the action changes only at higher adiabatic order.

For a one-dimensional particle of energy $E$ between walls separated by $L$, the unnormalized action isAdiabatic conservation gives $E(t)L(t{)}^2=\text{constant}$.

If an action has parameter dependence $I\propto \sqrt{m/|E|}$, adiabatic conservation implies $|E|\propto m$.

When a gravitational potential changes slowly relative to an orbital period, the orbital actions remain constant provided no separatrix is crossed. In a spherical potential, angular momentum is conserved exactly and the radial action is adiabatically conserved.

Adiabatic invariance of an orbital action preserves action values during slow regular changes, but frequencies $\partial H/\partial J_i$ can vary. Closure requires rational frequency ratios, which need not persist. A slowly added harmonic component in a Kepler orbit potential produces generic apsidal precession.

For radial motion between turning radii $r_{-}$ and $r_{+}$ in a central potential, the radial action per unit mass isIt vanishes on a circular orbit and is canonically conjugate to the radial angle.

For a bound Kepler orbit of semimajor axis $a$ around mass $M$,Adiabatic central-mass change therefore preserves $Ma$ and orbital eccentricity.

For $H=p^2/(2m)-1/|q|=-|E|$, collision continuation between the two turning points gives $I=\sqrt{2m/|E|}$.