Geodesic flow as a Hamiltonian flow
Anchor (Master): Jacobi 1837 (originator); Abraham-Marsden *Foundations of Mechanics* §3.7; Klingenberg *Riemannian Geometry* Ch. 3; Cannas da Silva *Lectures on Symplectic Geometry* §1 + §2
Intuition [Beginner]
A geodesic on a curved surface is the path you take when you walk straight without ever turning the wheel. On a sphere it bends back into a great circle; on a saddle it spreads apart from a neighbour; on the flat plane it is an ordinary straight line. The geodesics of a curved space record the shape of the space through the way "going straight" deforms.
There is a second story for the same paths: a particle gliding with no friction and no force, only the constraint of staying on the surface, follows a geodesic. So geodesics are also the free-motion trajectories of a frictionless point mass on the manifold. Two pictures, one set of paths.
The Hamiltonian framing puts both pictures into the same phase-space language used for every other mechanical system. Position is a point on the manifold; momentum is a covector — a way to pair tangent directions with numbers. The energy is half the squared length of the momentum, measured by the metric. Hamilton's equations turn that energy function into an evolution rule, and the resulting motion projects back down to a geodesic on the manifold. The metric becomes an energy function, the symplectic form becomes the rule for converting energy gradients into motion, and geodesic flow falls out as one specific Hamiltonian system.
Visual [Beginner]
A schematic showing a curved surface with a geodesic curve drawn on it, paired with a phase-space trajectory in the cotangent bundle above. An arrow labelled "project" connects a covector path back down to the geodesic on the surface.
The picture marks the two layers — the surface where geodesics live and the phase space where Hamilton's equations run — and the projection that ties them together.
Worked example [Beginner]
Take the flat plane with the ordinary distance. The metric components are simply and the inverse is the same. The energy function on phase space is
the standard kinetic energy with mass . Pick a starting point at the origin with momentum pointing in the first coordinate direction: position and momentum .
Hamilton's equations say position changes at the rate of momentum and momentum changes at the rate of minus the position-derivative of the energy. Here the energy does not depend on position, so momentum stays put: for every time . Position then advances at the constant rate , giving .
Projecting the phase-space curve back to the plane discards the momentum coordinate and keeps the trajectory . That is a straight line at unit speed — exactly the geodesic of the flat metric starting at the origin and pointing east. The Hamiltonian flow on phase space and the geodesic on the manifold tell the same story.
Check your understanding [Beginner]
Formal definition [Intermediate+]
Let be a smooth Riemannian manifold of dimension . The cotangent bundle carries the canonical Liouville one-form and the canonical symplectic form 05.02.05. In coordinates adapted to a chart on , and .
The metric on pulls back to a fibrewise inner product on via the inverse metric . The kinetic-energy Hamiltonian is the function
This is half the squared norm of the covector measured by the metric pulled to the cotangent space. Sign convention: , matching the conventions of 05.02.01 and 05.02.05.
The Hamiltonian vector field on is the unique smooth vector field satisfying . The cogeodesic flow is the flow of on . The geodesic flow in the classical sense is the corresponding flow on obtained by transporting across the Legendre isomorphism , , attached to the kinetic Lagrangian .
Key theorem with proof [Intermediate+]
Theorem (geodesic flow = cogeodesic flow). *Let be a Riemannian manifold and the kinetic-energy Hamiltonian on . The integral curves of on project under to geodesics of . Conversely, every geodesic of is the projection of a unique integral curve of specified by its initial position and initial momentum.*
Proof. Hamilton's equations associated with in canonical coordinates read
The first equation is the Legendre relation inverted, identifying momentum as the metric-flat of velocity. Differentiating it in ,
Substitute in the first term and in the second, then convert all 's to 's using :
Apply the inverse-metric derivative identity , valid because implies . After substitution and contraction the right-hand side collapses to
which is the geodesic equation in the Levi-Civita connection. The projected curve therefore satisfies the geodesic equation, and conversely every geodesic admits a horizontal lift to by the same Legendre relation, unique once initial momentum is fixed.
Bridge. The identification here builds toward 05.10.01 (contact manifold), where the level set — the unit cotangent bundle — appears again as the canonical contact manifold whose Reeb vector field is exactly the geodesic spray. Putting these together with the foliation picture of 05.02.04, the foundational reason the geodesic flow is integrable on highly symmetric manifolds (round spheres, flat tori, ellipsoids) is exactly that the Hamiltonian framing exposes Killing vector fields as Poisson-commuting integrals; the central insight is that metric symmetry is the same data as Hamiltonian symmetry once one passes to . Geodesic flow therefore appears again in 05.09.01 (KAM theorem) as the unperturbed integrable system for KAM-style perturbation analysis, and the bridge between the Riemannian-geometry and Hamiltonian-mechanics layers is that single energy function .
Exercises [Intermediate+]
Lean formalization [Intermediate+]
[object Promise]The proof in Mathlib would proceed by establishing the Legendre transform as a smooth bundle isomorphism , transporting the Hamiltonian vector field across it, and identifying the resulting vector field on as the geodesic spray. Each step is independently a candidate Mathlib contribution, listed in the lean_mathlib_gap field above.
Advanced results [Master]
The Hamiltonian framing of geodesic flow on is the natural setting for several deeper structural results.
Restriction to the unit cotangent bundle. The level set is a smooth -dimensional submanifold. The Liouville one-form restricts to a contact form : is a volume form on because is non-degenerate and the level set is transverse to . The Reeb vector field of — defined by and — coincides on with the restricted Hamiltonian vector field . So the geodesic flow on is, intrinsically, a Reeb flow of the canonical contact form on the unit cotangent bundle [Cannas da Silva §1].
Conservation laws and Killing vector fields. Energy is conserved because the Hamiltonian is autonomous: . More generally, every Killing vector field on — every infinitesimal isometry — induces a conserved quantity for the geodesic flow. The cotangent lift on is the Hamiltonian vector field of the function , and on the nose because preserves . On with the round metric, the rotation group supplies independent Killing fields, far more than the needed for Liouville-Arnold integrability; the geodesic flow on the round sphere is therefore superintegrable, and every geodesic is a great circle. On a flat torus the constant vector fields are Killing and supply commuting integrals, so the geodesic flow is Liouville integrable with the action-angle coordinates 05.02.04.
Anosov geodesic flow on negative curvature. On a closed Riemannian manifold with strictly negative sectional curvature — most cleanly, a compact hyperbolic surface — the geodesic flow on is Anosov: the tangent bundle of the flow direction admits an invariant splitting into uniformly contracting and uniformly expanding subspaces. The flow is exponentially mixing with respect to the Liouville measure on , and orbits are distributed equidistributionally. This is the bridge between symplectic / Riemannian geometry and ergodic theory: the same Hamiltonian flow that integrates explicitly on becomes the canonical example of a chaotic dynamical system on [Klingenberg Ch. 3].
Liouville-Arnold integrable cases. Beyond constant curvature, the Jacobi-integrable geodesic flow on the triaxial ellipsoid in is the classical example of an integrable geodesic flow with no underlying isometry-group transitivity. Jacobi 1838 Vorlesungen über Dynamik showed that elliptic-coordinate separation gives three independent Poisson-commuting integrals on (ellipsoid) — energy plus two extra — and the geodesics close up or wrap densely on Liouville tori according to rationality of the frequency ratios. Generalisations: the Neumann system, Stäckel-separable metrics, and the Manakov integrable cases on .
Length spectrum and Laplace spectrum. The set of lengths of closed geodesics on — the length spectrum — is the trace of the geodesic flow at periodic orbits. The Selberg trace formula on a compact hyperbolic surface gives an exact identity between the length spectrum and the Laplace spectrum of . The marked length spectrum (lengths together with free-homotopy classes) determines the metric on negatively curved closed surfaces (Otal, Croke, 1990). The inverse spectral question popularised by Mark Kac ("Can one hear the shape of a drum?") is the analogous question for the Laplace spectrum and is delicate even on flat tori (Milnor 1964 produced a 16-dimensional counterexample).
Symplectic capacities of unit codisk bundles. The unit codisk bundle is a symplectic manifold with boundary . Its symplectic capacities — Gromov width, Hofer-Zehnder capacity, Ekeland-Hofer capacities — are dynamical invariants of defined through the periodic orbits of the geodesic flow on the boundary. Floer-theoretic computations relate them to the volume and the systole. This connects geodesic flow to 05.07.02 symplectic capacity, 05.07.01 non-squeezing, and the symplectic homology of cotangent bundles (Viterbo 1996, Salamon-Weber 2006: , the homology of the free loop space).
Maupertuis-Jacobi reformulation. For a mechanical Hamiltonian with potential , on each energy level where the trajectories project to geodesics of the Jacobi metric , parametrised by -arc length. This converts a mechanical system with potential into a pure geodesic-flow problem at the cost of a metric depending on . The reformulation is due to Maupertuis 1744 in its variational form ("principle of least action") and Jacobi 1837 in its differential-geometric form; it is the bridge from Newtonian mechanics with conservative forces to Riemannian-geometry techniques.
Synthesis. The data of a Riemannian metric and the data of a kinetic-energy Hamiltonian on the cotangent bundle are equivalent: each determines the other via the Legendre transform , and the dynamics agree under projection. Read in one direction, the construction puts geodesic flow into the symplectic toolkit, where Liouville volume preservation, action-angle coordinates, KAM perturbation theory, and Floer homology become available techniques. Read in the opposite direction, every kinetic Hamiltonian on a cotangent bundle is the geodesic flow of some metric — its Hamiltonian comes from a quadratic form on covectors, and any positive-definite quadratic form defines a Riemannian metric. The bridge between the analytic equation and the geometric conclusion that integral curves project to geodesics is the foundational reason the symplectic, contact, and Riemannian pictures of free-particle motion coincide. Putting these together, one sees that the moduli space of "kinetic Hamiltonians on " is the moduli space of Riemannian metrics on , and the central insight is that metric symmetry, Hamiltonian symmetry, and contact symmetry are the same data viewed from three angles.
Full proof set [Master]
Theorem (geodesic flow = cogeodesic flow). Restated: the integral curves of for on project to geodesics of on , and this correspondence is bijective on initial conditions specified by point and tangent vector via the Legendre transform .
Proof. Given in the Intermediate Key theorem section above. The substitution chain , , combined with the inverse-metric derivative identity, produces , the Levi-Civita geodesic equation. The bijection on initial conditions is the assertion that is a diffeomorphism, which holds because the kinetic Lagrangian is hyper-regular (the Hessian is positive-definite).
Theorem (Reeb = geodesic spray on unit cotangent bundle). *Let be the restriction of the canonical Liouville form to the unit cotangent bundle. Then is a contact form, and its Reeb vector field is the restriction of .*
Proof. Non-degeneracy of on together with transversality of to — which holds because but off the zero section — implies is a volume form on . So is contact. For the Reeb identification, on one has , and . Restrict to and observe , so . The pair uniquely determines the Reeb vector field, and satisfies both.
Theorem (Killing → conserved quantity). Let be a Killing vector field on , with cotangent lift function . Then , and is conserved along the cogeodesic flow.
Proof. In coordinates . Compute the Poisson bracket:
Symmetrise the first term in : the coefficient of is . The Killing equation (equivalently ) translates in the inverse metric to , which in coordinates is exactly . So the bracket vanishes. By Hamilton's equations .
Theorem (Maupertuis-Jacobi). On the energy level of the mechanical Hamiltonian , the trajectories of coincide as unparametrised curves with the trajectories of the cogeodesic flow of the Jacobi metric , restricted to .
Proof. Solution given in Exercise 7. The crucial identity is that on one has , and , so on this set, identifying it with . The vector fields and are parallel on the common level set, so trajectories coincide as point sets. Reparametrisation by -arc length identifies -geodesics with the original mechanical trajectories.
Connections [Master]
Cotangent bundle
05.02.05. The Hamiltonian framing of geodesic flow is built directly on the canonical symplectic structure of . The Liouville one-form , the symplectic form , and the projection all play load-bearing roles in the construction, and the kinetic-energy Hamiltonian is the simplest non-zero function on that requires the full Riemannian data.Hamiltonian vector field
05.02.01. Geodesic flow is the canonical example of a Hamiltonian vector field on a cotangent bundle. The general theory of — energy conservation, Liouville volume preservation, Cartan-formula identities — specialises to give the standard properties of geodesic flow without separate calculation.Action-angle coordinates
05.02.04. On a manifold whose isometry group acts transitively (round spheres, flat tori, ellipsoids treated by Jacobi), the geodesic flow is Liouville integrable, and action-angle coordinates linearise the flow on each Liouville torus. This is the direct route from Riemannian symmetry to Hamiltonian integrability.Contact manifold
05.10.01. The unit cotangent bundle is the canonical contact manifold attached to , and the geodesic flow is its Reeb flow. The contact-geometry / Riemannian-geometry bridge runs through this identification.KAM theorem
05.09.01. A small Hamiltonian perturbation of an integrable geodesic flow — for example, a small symmetric perturbation of the round metric on — produces a Hamiltonian system to which KAM applies: most invariant Liouville tori survive, and the perturbed geodesic flow remains quasi-periodic on a positive-measure set. Geodesic flow on highly symmetric metrics is the canonical "unperturbed system" for KAM-style analysis.Symplectic reduction
05.04.02. When a Lie group acts on by isometries, the cotangent lift acts on by symplectomorphisms preserving . The Marsden-Weinstein reduced spaces are then phase spaces for the reduced geodesic flow on the quotient orbifold. Examples: geodesic flow on as reduction of the free-particle flow on .Symplectic capacity
05.07.02and non-squeezing05.07.01. Symplectic capacities of unit codisk bundles of are quantitative invariants computed from periodic orbits of the geodesic flow. Floer-theoretic computations of capacities of recover information about the Riemannian metric — volume, systole, length spectrum.
Historical & philosophical context [Master]
Carl Gustav Jacob Jacobi gave the first identification of geodesic flow with Hamiltonian flow in 1837, in Note sur l'intégration des équations différentielles de la dynamique (Comptes Rendus Acad. Sci. Paris) [Jacobi 1837]. The 1837 note announced the principle: trajectories of a conservative mechanical system on a configuration manifold reduce, after Maupertuis reparametrisation, to geodesics of a metric depending on the energy level. The full proofs appeared in Jacobi's Vorlesungen über Dynamik (1842-43 lectures, published posthumously 1866), where elliptic-coordinate separation of variables produced the integration of the geodesic flow on the triaxial ellipsoid — the first integrable geodesic flow established by techniques other than direct symmetry reduction.
The modern coordinate-free framing — Hamiltonian as a function on , Liouville form , symplectic form , kinetic energy — was crystallised in the symplectic-geometry programme of Souriau, Smale, Abraham, and Marsden in the 1960s and 1970s. Abraham-Marsden Foundations of Mechanics (1978) [Abraham-Marsden] §3.7 gives the canonical exposition of the cogeodesic flow on and its Legendre relation to the geodesic spray on , in the language used by every subsequent textbook.
Klingenberg's Riemannian Geometry (1982) [Klingenberg] develops the geodesic spray and the symplectic structure of in the Riemannian-geometry tradition, providing the translation between the symplectic-mechanics and the Riemannian-geometry literatures. Cannas da Silva's Lectures on Symplectic Geometry [Cannas da Silva §1] places the cogeodesic flow as the headline example of a cotangent-bundle Hamiltonian system in the modern symplectic-geometry pedagogical canon.
The connection between geodesic flow and ergodic theory was developed by Hadamard 1898 (negative curvature on surfaces of finite type produces dense geodesics), Hopf 1939 (mixing of the geodesic flow on closed surfaces of negative curvature), and Anosov 1967 (uniform hyperbolicity in arbitrary dimension), establishing geodesic flow as the canonical example of a chaotic dynamical system. The Selberg trace formula (1956) and the Otal-Croke length-spectrum rigidity theorems (1990) place the closed-geodesic spectrum at the centre of inverse spectral geometry.