Skip to content

Hamiltonian Mechanics: Advanced Topics

10.1 Canonical Transformations

A canonical transformation from (q,p)(q, p) to (Q,P)(Q, P) preserves the form of Hamilton’s equations. It is generated by a generating function FF.

Type 1 (F1(q,Q,t)F_1(q, Q, t)): p=F1/qp = \partial F_1/\partial q, P=F1/QP = -\partial F_1/\partial Q, K=H+F1/tK = H + \partial F_1/\partial t.

Type 2 (F2(q,P,t)F_2(q, P, t)): p=F2/qp = \partial F_2/\partial q, Q=F2/PQ = \partial F_2/\partial P, K=H+F2/tK = H + \partial F_2/\partial t.

Type 3 (F3(p,Q,t)F_3(p, Q, t)): q=F3/pq = -\partial F_3/\partial p, P=F3/QP = -\partial F_3/\partial Q, K=H+F3/tK = H + \partial F_3/\partial t.

Type 4 (F4(p,P,t)F_4(p, P, t)): q=F4/pq = -\partial F_4/\partial p, Q=F4/PQ = \partial F_4/\partial P, K=H+F4/tK = H + \partial F_4/\partial t.

The transformation is canonical if and only if:

ipidqiiPidQi=dF\sum_i p_i\,dq_i - \sum_i P_i\,dQ_i = dF

Poisson bracket invariance. Canonical transformations preserve the Poisson bracket: {F,G}q,p={F,G}Q,P\{F, G\}_{q,p} = \{F, G\}_{Q,P}.

Example. The point transformation Qi=qiQ_i = q_i (just relabelling coordinates) with Pi=piP_i = p_i is trivially canonical. A non-trivial example: the rotation Q=qcosθ+psinθQ = q\cos\theta + p\sin\theta, P=qsinθ+pcosθP = -q\sin\theta + p\cos\theta is canonical (it preserves dqdp=dQdPdq \wedge dp = dQ \wedge dP).

10.2 Action-Angle Variables

For a periodic system with frequency ω\omega, define the action:

Ji=pidqiJ_i = \oint p_i\,dq_i

The conjugate angle variable θi\theta_i evolves linearly: θi(t)=ωit+θi(0)\theta_i(t) = \omega_i t + \theta_i(0).

Properties:

  • JiJ_i is adiabatically invariant (conserved under slow parameter changes).
  • The motion is trivial in action-angle variables: Ji=constJ_i = \text{const}, θ˙i=ωi\dot{\theta}_i = \omega_i.
  • The transformation to action-angle variables is a canonical transformation.

Application: Kepler problem. In the Kepler problem, the action variables are: J1=LJ_1 = L (angular momentum), J2=L+LzJ_2 = L + L_z (related to eccentricity), J3=E/ωJ_3 = -E/\omega (related to energy). The angle variables give the orbital elements.

10.3 Hamilton—Jacobi Equation

The Hamilton—Jacobi equation is a reformulation of Hamiltonian mechanics as a PDE for Hamilton’s principal function S(q,α,t)S(q, \alpha, t):

H ⁣(q,Sq,t)+St=0H\!\left(q, \frac{\partial S}{\partial q}, t\right) + \frac{\partial S}{\partial t} = 0

If SS can be found by separation of variables, the transformation to new coordinates makes all momenta constant, effectively solving the problem.

Separation for the harmonic oscillator. For H=p2/(2m)+mω2q2/2H = p^2/(2m) + m\omega^2 q^2/2:

S(q,E,t)=2m(E12mω2q2)dqEtS(q, E, t) = \int \sqrt{2m(E - \frac{1}{2}m\omega^2 q^2)}\,dq - Et

The new momentum P1=α1=EP_1 = \alpha_1 = E (constant), and the new coordinate Q1=S/EQ_1 = \partial S/\partial E gives a linear function of time.

Maupertuis principle. For energy-conserving systems, the Hamilton—Jacobi equation simplifies to:

2m[EV(q)]=(Wq)22m\left[E - V(q)\right] = \left(\frac{\partial W}{\partial q}\right)^2

where WW is the abbreviated action. This is equivalent to Fermat’s principle in optics.

10.4 Adiabatic Invariants

An adiabatic invariant is a quantity that remains constant when a parameter of the system is changed slowly compared to the period of motion.

For a harmonic oscillator with slowly varying ω(t)\omega(t):

Eω=const(adiabatic invariant)\frac{E}{\omega} = \text{const} \quad \text{(adiabatic invariant)}

This has important applications:

  • Fermi acceleration: Cosmic rays gaining energy from moving magnetic fields
  • Magnetic mirrors: Charged particles trapped between magnetic mirrors conserve μ=mv2/(2B)\mu = mv_\perp^2/(2B) (the magnetic moment)
  • Planetary orbits: The semi-major axis is an adiabatic invariant under slow mass loss from the Sun
  • Pendulum with slowly changing length: If the length l(t)l(t) changes slowly, E/g/lE/\sqrt{g/l} is conserved.

Derivation (sketch). For a periodic Hamiltonian H(q,p,λ(t))H(q, p, \lambda(t)) with λ˙ω\dot{\lambda} \ll \omega, the action J=pdqJ = \oint p\,dq satisfies J˙=O(λ˙/ω)\dot{J} = O(\dot{\lambda}/\omega), so JJ is approximately constant over many periods.

10.5 Liouville’s Theorem and Phase Space

Liouville’s theorem: The phase space distribution function ρ(q,p,t)\rho(q, p, t) is constant along trajectories:

dρdt=ρt+{ρ,H}=0\frac{d\rho}{dt} = \frac{\partial\rho}{\partial t} + \{\rho, H\} = 0

This means phase space volume is conserved: a region of phase space evolves like an incompressible fluid.

Consequences:

  • The phase space density of an ensemble of systems cannot increase
  • This constrains the focusing of particle beams in accelerators
  • The theorem underpins the ergodic hypothesis in statistical mechanics
  • Poincare recurrence: bounded Hamiltonian systems return arbitrarily close to their initial state

Proof. The flow ϕt\phi_t generated by Hamilton’s equations satisfies det((q,p)/(q,p))=1\det(\partial(q', p')/\partial(q, p)) = 1 (canonical transformation), so the Jacobian is 1. By the change of variables formula, ρ\rho is conserved along trajectories. \square

10.6 Perturbation Theory and KAM Theorem

Perturbation theory addresses Hamiltonian systems close to integrable ones: H=H0+εH1H = H_0 + \varepsilon H_1.

Lindstedt—Poincare method. Expand both the solution and the frequency in powers of ε\varepsilon to avoid secular terms.

KAM theorem (Kolmogorov—Arnold—Moser). For sufficiently small ε\varepsilon, most invariant tori of the integrable system persist under perturbation, provided the frequency vector satisfies a Diophantine condition (sufficiently irrational).

The surviving tori have positive measure for small ε\varepsilon, but the remaining phase space contains chaotic regions. This explains the stability of the solar system over billions of years (integrable Kepler problem + small planetary perturbations).

Common Pitfalls

  1. Using the wrong generating function type. Each type of generating function expresses different variables as independent. Choosing the wrong type leads to implicit equations that may not be solvable.

  2. Forgetting the time dependence in generating functions. If FF depends explicitly on tt, the new Hamiltonian is K=H+F/tK = H + \partial F/\partial t, not K=HK = H.

  3. Confusing canonical with point transformations. Not all coordinate changes are canonical. A transformation is canonical only if it preserves the symplectic structure dpdqdp \wedge dq.

  4. Applying perturbation theory to non-integrable systems. The KAM theorem requires the unperturbed system to be integrable. For chaotic systems, perturbation theory fails.

Worked Example 10.1: Action-Angle Variables for the Harmonic Oscillator

For the 1D harmonic oscillator: H=p2/(2m)+12mω2q2H = p^2/(2m) + \frac{1}{2}m\omega^2 q^2.

The action variable:

J=pdq=2mEm2ω2q2dqJ = \oint p\,dq = \oint \sqrt{2mE - m^2\omega^2 q^2}\,dq

The contour is the ellipse p2/(2mE)+q2/(2E/mω2)=1p^2/(2mE) + q^2/(2E/m\omega^2) = 1 with semi-axes 2mE\sqrt{2mE} and 2E/(mω2)\sqrt{2E/(m\omega^2)}.

The area (and hence the action):

J=π×2mE×2Emω2=2πEωJ = \pi \times \sqrt{2mE} \times \sqrt{\frac{2E}{m\omega^2}} = \frac{2\pi E}{\omega}

So E=Jω/2E = J\omega/2 and the Hamiltonian in action-angle form is:

H(J)=JωH(J) = J\omega

The angle variable evolves as θ˙=H/J=ω\dot{\theta} = \partial H/\partial J = \omegaGiving θ(t)=ωt+θ0\theta(t) = \omega t + \theta_0.

The frequency is ω=H/J=const\omega = \partial H/\partial J = \text{const}Independent of JJ (harmonic oscillator has no frequency shift with amplitude --- a special property).