Learning path

Full curriculum

Full curriculum

Unit content

Symplectic integration of Hamiltonian dynamics

Ordinary time-stepping methods can approximate a trajectory accurately for a short time yet distort the geometry of a conservative mechanical system over many cycles. Symplectic integrators are designed to preserve the phase-space structure of Hamiltonian dynamics.

For a separable Hamiltonian $$H(q,p)=T(p)+V(q),$$ Hamilton's equations are $$\dot q=\frac{\partial T}{\partial p},\qquad \dot p=-\frac{\partial V}{\partial q}.$$ A simple symplectic method is the leapfrog or velocity-Verlet scheme: $$p_{n+1/2}=p_n-\frac{\Delta t}{2}\nabla V(q_n),$$ $$q_{n+1}=q_n+\Delta t,\nabla_pT(p_{n+1/2}),$$ $$p_{n+1}=p_{n+1/2}-\frac{\Delta t}{2}\nabla V(q_{n+1}).$$

For a harmonic oscillator, a generic non-symplectic method may cause the numerical energy to drift steadily upward or downward. Leapfrog instead typically produces a small bounded oscillation of the energy around the correct value. It does not conserve the exact Hamiltonian at every step; it preserves a nearby Hamiltonian and the symplectic geometry, which gives much better long-time qualitative behavior.

Symplectic methods matter when long-term phase accuracy and conservation structure are more important than minimizing one-step error: orbital dynamics, molecular dynamics, accelerator physics, and many-body Hamiltonian systems are common examples.

They are not automatically superior for dissipative systems or arbitrary differential equations. Their advantage comes from matching the mathematical structure of Hamiltonian evolution.