In mathematics, the geodesic equations are second-order non-linear differential equations, and are commonly presented in the form of Euler–Lagrange equations of motion. However, they can also be presented as a set of coupled first-order equations, in the form of Hamilton's equations. This latter formulation is developed in this article.
It is frequently said that geodesics are "straight lines in curved space". By using the Hamilton–Jacobi approach to the geodesic equation, this statement can be given a very intuitive meaning: geodesics describe the motions of particles that are not experiencing any forces. In flat space, it is well known that a particle moving in a straight line will continue to move in a straight line if it experiences no external forces; this is Newton's first law. The Hamiltonian describing such motion is well known to be
H=p2/2m
\gamma:I\toM
that maps an interval I of the real number line to the manifold M, one writes the energy
E(\gamma)= | 1 |
2 |
\intIg(
|
where
\gamma(t) |
\gamma
t\inI
g( ⋅ , ⋅ )
d2xa | |
dt2 |
+
a | |
\Gamma | |
bc |
dxb | |
dt |
dxc | |
dt |
=0
where the xa(t) are the coordinates of the curve γ(t),
a | |
\Gamma | |
bc |
Geodesics can be understood to be the Hamiltonian flows of a special Hamiltonian vector field defined on the cotangent space of the manifold. The Hamiltonian is constructed from the metric on the manifold, and is thus a quadratic form consisting entirely of the kinetic term.
The geodesic equations are second-order differential equations; they can be re-expressed as first-order equations by introducing additional independent variables, as shown below. Note that a coordinate neighborhood U with coordinates xa induces a local trivialization of
*M| | |
T | |
U |
\simeqU x Rn
η\in
*M| | |
T | |
U |
η=padxa
(x,pa)\inU x Rn
H(x,p)= | 1 |
2 |
gab(x)papb.
Here, gab(x) is the inverse of the metric tensor: gab(x)gbc(x) =
a | |
\delta | |
c |
x |
a=
\partialH | |
\partialpa |
=gab(x)pb
and
p |
a=-
\partialH | |
\partialxa |
=-
1 | |
2 |
\partialgbc(x) | |
\partialxa |
pbpc.
The flow determined by these equations is called the cogeodesic flow; a simple substitution of one into the other obtains the Euler–Lagrange equations, which give the geodesic flow on the tangent bundle TM. The geodesic lines are the projections of integral curves of the geodesic flow onto the manifold M. This is a Hamiltonian flow, and the Hamiltonian is constant along the geodesics:
dH | |
dt |
=
\partialH | |
\partialxa |
x |
a+
\partialH | |
\partialpa |
p |
a=-
p |
a
x |
a+
x |
a
p |
a=0.
Thus, the geodesic flow splits the cotangent bundle into level sets of constant energy
ME=\{(x,p)\inT*M:H(x,p)=E\}
for each energy E ≥ 0, so that
*M=cup | |
T | |
E\ge0 |
ME