Hamiltonian Vector Fields, Poisson Brackets & Liouville's Theorem

The dictionary: a symplectic form turns functions into vector fields

From Guide 1 you carry one structure: a symplectic manifold (M, omega), an even-dimensional smooth manifold with a closed (d omega = 0) and nondegenerate 2-form omega. Nondegeneracy is the workhorse here, so let us squeeze everything out of it. At each point p, omega is an antisymmetric, invertible bilinear pairing on the tangent space T_p M. Invertible means it gives a perfect way to convert vectors into covectors: feed a vector X into one slot of omega and you get a covector omega(X, -) sitting in the cotangent space. Because omega is nondegenerate this map TM -> T*M, X -> omega(X, -), is a bundle isomorphism — it has an inverse. That single isomorphism is the engine of all of Hamiltonian mechanics.

Now run the engine backwards. Take any smooth function H: M -> R — physically an energy, an observable, a 'Hamiltonian'. Its differential dH is a 1-form, i.e. a section of T*M. Since omega lets us convert covectors back into vectors, dH names a unique vector field. We call it the Hamiltonian vector field of H, written X_H, defined by the single equation iota_{X_H} omega = dH, that is omega(X_H, -) = dH(-). Read it as a dictionary entry: the symplectic form is the dictionary, a function H is a word, and X_H is its translation into the language of flows. A Riemannian metric would do the analogous trick and hand you the gradient field grad H; the symplectic version is its skew twin, and the skewness is exactly what makes the dynamics conservative rather than dissipative.

Standard R^{2n} with omega_0 = sum dp_i ^ dq_i  (Darboux coordinates, Guide 1)

  X_H defined by:   omega_0(X_H, -) = dH
  works out to:     X_H = sum ( dH/dp_i  d/dq_i  -  dH/dq_i  d/dp_i )

  flow of X_H = Hamilton's equations:
        dq_i/dt =  dH/dp_i
        dp_i/dt = -dH/dq_i

In Darboux coordinates the abstract definition iota_{X_H} omega = dH unfolds into the textbook Hamilton equations — the sign-swapping antisymmetry of omega is exactly the minus sign that conserves energy.

Hamiltonian flows preserve both energy and the symplectic form

A vector field is only as good as the flow it generates, so let phi_t be the flow of X_H — push every point along X_H for time t. Two conservation laws drop out almost for free, and both are one-line consequences of antisymmetry. First, energy is conserved: the rate of change of H along its own flow is X_H(H) = dH(X_H) = omega(X_H, X_H) = 0, because omega is antisymmetric and any vector paired with itself gives zero. So H is constant along every trajectory — the flow stays glued to a level set {H = c}. This is why you draw phase portraits as curves living on energy shells, never crossing them.

Second, and more structurally, the flow preserves omega itself: each phi_t is a symplectomorphism, Guide 1's structure-preserving maps. The cleanest proof is the Cartan magic formula, L_X = iota_X d + d iota_X, applied to the Lie derivative of omega along X_H. Watch it collapse: L_{X_H} omega = iota_{X_H} (d omega) + d(iota_{X_H} omega) = iota_{X_H} 0 + d(dH) = 0 + 0 = 0. The first term vanishes because omega is closed (d omega = 0) — this is the hidden job of closedness — and the second vanishes because d^2 = 0. So the Lie derivative is zero, which means omega is dragged along unchanged by the flow. Hamiltonian dynamics is precisely the dynamics that respects the symplectic form.

The Poisson bracket: observables become a Lie algebra

We now have a map from functions to vector fields, H -> X_H. The next move is to bring it back down to functions and discover hidden algebra. Given two observables f and g, define their Poisson bracket by {f, g} = omega(X_f, X_g) — feed the two Hamiltonian vector fields into the symplectic form. Equivalently {f, g} = X_f(g) = dg(X_f), so it measures the rate of change of g along the Hamiltonian flow of f. In Darboux coordinates this is the familiar {f, g} = sum (df/dq_i dg/dp_i - df/dp_i dg/dq_i). The bracket is antisymmetric, {f, g} = -{g, f} (from omega's antisymmetry), and bilinear over R.

The deep fact is that the bracket satisfies the Jacobi identity, {f, {g, h}} + {g, {h, f}} + {h, {f, g}} = 0 — and this is precisely where closedness of omega is spent again, d omega = 0 being equivalent to Jacobi. So the space of smooth functions C^infinity(M) becomes an infinite-dimensional Lie algebra under the Poisson bracket, with the same abstract shape as the Lie algebras of Lie groups you met earlier in this volume. Better still, the assignment f -> X_f is a Lie algebra homomorphism: X_{{f,g}} = -[X_f, X_g], turning the Poisson bracket of functions into the Lie bracket of vector fields (the sign depends on your convention). Two layers of bracket, one dictionary connecting them.

This algebra is not decoration — it is the language of conserved quantities. Since {f, g} = X_f(g), a function g is conserved along the flow of f exactly when {f, g} = 0. In particular {H, g} = 0 says g is a constant of motion for the energy H. And {H, H} = 0 re-derives energy conservation in one symbol. Noether's theorem, the moment maps of Guide 4, and the whole bookkeeping of symmetry all live here: symmetries are functions that Poisson-commute with H, and the bracket is how you detect, combine, and propagate them. When you reach the moment map you will see the Poisson bracket reappear as the structure that makes a group action 'Hamiltonian'.

Liouville's theorem: phase-space volume cannot be squeezed

Now collect the dividend. On a 2n-dimensional symplectic manifold the top wedge power of omega, namely omega^n = omega ^ omega ^ ... ^ omega (n copies), is nowhere zero — nondegeneracy guarantees it — so it is a volume form, the Liouville volume (up to a constant 1/n!). We already proved that the Hamiltonian flow preserves omega; wedging omega with itself, it therefore preserves omega^n too. That is Liouville's theorem: the flow of any Hamiltonian vector field preserves the phase-space volume. A blob of initial conditions can stretch, fold, and filament into wild shapes, but its total 2n-dimensional volume is exactly the same at every later time.

Picture it in the smallest case, R^2 with coordinates (q, p) and omega_0 = dp ^ dq, where the Liouville form is just ordinary area. A Hamiltonian flow there is an area-preserving flow of the plane: the pendulum's phase portrait swirls, but any little parallelogram of initial states keeps its area forever. The harmonic oscillator H = (p^2 + q^2)/2 makes X_H rotate the plane rigidly — circles of constant energy, area trivially preserved. A more chaotic H shears the blob into a thin tendril wrapping the energy shell, yet the area is pinned. This single picture is the geometric heart of statistical mechanics: it is why a uniform probability density on phase space stays uniform, the foundation under the microcanonical ensemble.

Read Liouville's theorem with honest fine print. (1) It preserves the 2n-dimensional Liouville volume, NOT lengths, NOT areas of intermediate dimension, and NOT shape — blobs genuinely distort and that distortion is the whole content of chaos; conservation of volume coexists with violent mixing. (2) It needs the flow to be Hamiltonian; a dissipative system (friction) is governed by a non-Hamiltonian field whose flow shrinks volume onto attractors — that is the whole point of an attractor, and it is why real engines lose energy. (3) Liouville is about the FULL phase-space volume; the much harder statement that you cannot squeeze a symplectic ball through a thin neck even while keeping volume constant is Gromov's non-squeezing theorem, a genuinely 2-dimensional symplectic rigidity we only glimpse in Guide 5 — it shows omega remembers more than volume alone. (4) Sign and ordering conventions for {f, g} and for iota_{X_H} omega = +dH vs -dH differ across books (Arnold, da Silva, Marsden); fix yours before quoting a formula, as the same warning from Guide 1 still applies.

Putting the dictionary to work

Everything in this guide is one dictionary read in different directions, so it pays to rehearse the round trip explicitly. A function H gives a covector field dH; nondegeneracy of omega converts it into the Hamiltonian vector field X_H; the flow of X_H gives the dynamics and, by Darboux's theorem from Guide 1, that flow looks locally like Hamilton's equations in canonical coordinates. Pairing two such fields back through omega returns to functions as the Poisson bracket; closedness of omega makes that bracket a Lie algebra and the flows symplectomorphisms; and wedging omega up to top degree delivers Liouville's volume conservation. Hold the whole chain in view and the subject stops being a list of named theorems and becomes one idea seen from several sides.

Pick a Hamiltonian: choose a smooth function H on (M, omega) — an energy or observable — and compute its differential dH.
Solve for the field: find the unique X_H with iota_{X_H} omega = dH; in Darboux coordinates this just reads off Hamilton's equations dq/dt = dH/dp, dp/dt = -dH/dq.
Hunt conserved quantities: any g with {H, g} = 0 is a constant of motion; in particular {H, H} = 0 re-proves energy conservation, and Poisson-commuting functions stack into symmetries.
Track volume: form the Liouville form omega^n; since the flow preserves omega it preserves omega^n, so phase-space volume is invariant even as shapes distort — Liouville's theorem in action.

Where does this lead? The conserved-quantity story matures into Guide 4's moment map and Marsden-Weinstein reduction, where a whole symmetry group's worth of Poisson-commuting functions is packaged into a single map and used to shrink the phase space. The level sets {H = c} and the graphs of these flows are themselves geometric objects — and the most important ones are Lagrangian, which is exactly Guide 3's subject. Keep the picture of a function generating a volume-preserving flow firmly in mind; it is the steady drumbeat under everything that follows in this rung.