Riemannian Metrics and Geodesics

⚡ Key Takeaways

A Riemannian metric is a smoothly varying inner product on the tangent spaces of a manifold, and geodesics are its locally length-minimizing curves — the manifold's intrinsic notion of "straight lines." Key result: geodesics satisfy a second-order ODE (the geodesic equation) determined by the Christoffel symbols of the metric, and they are exactly the critical points of the energy/length functional. Applies wherever curved spaces appear: general relativity, geometric optics, optimization on manifolds, shape analysis, and pure differential geometry. Introduced by Bernhard Riemann (1854) and developed by Christoffel, Ricci, Levi-Civita, and Cartan.

💡 The Core Insight

Give each tangent space an inner product, and all of metric geometry follows: lengths via integration of speed, distance via infimum over paths, "straight lines" via the variational principle, and curvature via how nearby geodesics fail to stay parallel. The metric is the only input; the Levi-Civita connection, geodesics, and curvature are uniquely forced by compatibility and symmetry.

Formal Statement

DefinitionRiemannian Metric

Let MM be a smooth nn-dimensional manifold. A Riemannian metric gg is a smooth (0,2)(0,2)-tensor field such that each gp:TpM×TpM→Rg_p : T_pM \times T_pM \to \mathbb{R} is symmetric and positive-definite. In local coordinates (x1,…,xn)(x^1,\dots,x^n) with coordinate basis ∂i=∂/∂xi\partial_i = \partial/\partial x^i, the metric is encoded by the component matrix gijg_{ij}:

1
g=gij dxi⊗dxj,gij(p)=gp(∂i,∂j)=gji(p).\quad g = g_{ij}\, dx^i \otimes dx^j, \qquad g_{ij}(p) = g_p(\partial_i,\partial_j) = g_{ji}(p).

Derivation hint: Evaluate the bilinear form on coordinate basis vectors; symmetry of gg gives gij=gjig_{ij}=g_{ji}, positive-definiteness gives gijvivj>0g_{ij}v^iv^j>0 for v≠0v\neq 0. The pair (M,g)(M,g) is a Riemannian manifold. The squared norm and inner product of tangent vectors u,vu,v at pp are

2
⟨u,v⟩g=gij uivj,∣v∣g=gij vivj.\quad \langle u,v\rangle_g = g_{ij}\,u^i v^j, \qquad |v|_g = \sqrt{g_{ij}\,v^i v^j}.

Derivation hint: Expand u=ui∂iu=u^i\partial_i, v=vj∂jv=v^j\partial_j, apply bilinearity of gg, and use g(∂i,∂j)=gijg(\partial_i,\partial_j)=g_{ij}.

DefinitionLength, Energy, and Riemannian Distance

For a piecewise-smooth curve γ:[a,b]→M\gamma:[a,b]\to M, define its length and energy:

3
L(γ)=∫ab∣γ˙(t)∣g dt=∫abgij(γ) γ˙iγ˙j dt,E(γ)=12∫abgij(γ) γ˙iγ˙j dt.\quad L(\gamma)=\int_a^b |\dot\gamma(t)|_g\,dt=\int_a^b \sqrt{g_{ij}(\gamma)\,\dot\gamma^i\dot\gamma^j}\,dt, \qquad E(\gamma)=\tfrac12\int_a^b g_{ij}(\gamma)\,\dot\gamma^i\dot\gamma^j\,dt.

Derivation hint: LL integrates instantaneous speed ∣γ˙∣|\dot\gamma|; EE uses the squared integrand, which is smoother and reparametrization-sensitive. The Riemannian distance between p,q∈Mp,q\in M is the infimum of lengths over connecting curves:

4
dg(p,q)=inf⁡{L(γ):γ piecewise-smooth, γ(a)=p, γ(b)=q}.\quad d_g(p,q)=\inf\{L(\gamma):\gamma \text{ piecewise-smooth},\ \gamma(a)=p,\ \gamma(b)=q\}.

Derivation hint: Length is reparametrization-invariant, so the infimum is well-defined; (M,dg)(M,d_g) is a metric space whose topology equals the manifold topology.

DefinitionLevi-Civita Connection and Christoffel Symbols

There is a unique affine connection ∇\nabla that is torsion-free (∇XY−∇YX=[X,Y]\nabla_X Y - \nabla_Y X = [X,Y]) and metric-compatible (∇g=0\nabla g = 0): the Levi-Civita connection. Its Christoffel symbols in coordinates are

5
Γijk=12 gkℓ(∂igjℓ+∂jgiℓ−∂ℓgij),\quad \Gamma^k_{ij}=\tfrac12\, g^{k\ell}\big(\partial_i g_{j\ell}+\partial_j g_{i\ell}-\partial_\ell g_{ij}\big),

Derivation hint: Impose ∇g=0\nabla g=0 on cyclic permutations of indices (i,j,ℓ)(i,j,\ell), subtract one equation from the sum of the other two, and use Γijk=Γjik\Gamma^k_{ij}=\Gamma^k_{ji} (torsion-free) to solve; gkℓg^{k\ell} is the inverse metric. where gkℓg^{k\ell} denotes the inverse matrix, gkℓgℓm=δmkg^{k\ell}g_{\ell m}=\delta^k_m.

TheoremGeodesic Equation

A curve γ\gamma is a geodesic (zero covariant acceleration, ∇γ˙γ˙=0\nabla_{\dot\gamma}\dot\gamma=0) iff its components satisfy

6
γ¨ k+Γijk(γ) γ˙iγ˙j=0,k=1,…,n.\quad \ddot\gamma^{\,k} + \Gamma^k_{ij}(\gamma)\,\dot\gamma^i\dot\gamma^j = 0, \qquad k=1,\dots,n.

Derivation hint: Write ∇γ˙γ˙=(γ¨k+Γijkγ˙iγ˙j)∂k\nabla_{\dot\gamma}\dot\gamma=(\ddot\gamma^k+\Gamma^k_{ij}\dot\gamma^i\dot\gamma^j)\partial_k using ∇∂i∂j=Γijk∂k\nabla_{\partial_i}\partial_j=\Gamma^k_{ij}\partial_k and the chain rule for γ˙=γ˙i∂i\dot\gamma=\dot\gamma^i\partial_i. Geodesics are precisely the critical points of the energy functional EE, and (after arclength reparametrization) the locally length-minimizing curves. They have constant speed: ∣γ˙∣g|\dot\gamma|_g is conserved along γ\gamma.

Why It Holds / Derivation

ProofGeodesics from the Energy Functional

We show that critical points of EE satisfy the geodesic equation, justifying the second-order ODE intrinsically.

Step 1 — Variation setup. Let γs(t)=γ(t,s)\gamma_s(t)=\gamma(t,s) be a smooth one-parameter family with fixed endpoints, γ0=γ\gamma_0=\gamma, variation field V=∂sγ∣s=0V=\partial_s\gamma|_{s=0} vanishing at a,ba,b. The energy is

7
E(γs)=12∫abgij(γs) ∂tγsi ∂tγsj dt.\quad E(\gamma_s)=\tfrac12\int_a^b g_{ij}(\gamma_s)\,\partial_t\gamma_s^i\,\partial_t\gamma_s^j\,dt.

Derivation hint: Substitute the family into the energy integral [3]; differentiation in ss will produce the first variation.

Step 2 — Differentiate in ss. Using ∂sgij=∂ℓgij Vℓ\partial_s g_{ij}=\partial_\ell g_{ij}\,V^\ell and ∂s∂tγi=∂tVi\partial_s\partial_t\gamma^i=\partial_t V^i:

8
dds∣0E=∫ab(12 ∂ℓgij Vℓ γ˙iγ˙j+gij γ˙i V˙j) dt.\quad \frac{d}{ds}\Big|_{0}E=\int_a^b \Big(\tfrac12\,\partial_\ell g_{ij}\,V^\ell\,\dot\gamma^i\dot\gamma^j + g_{ij}\,\dot\gamma^i\,\dot V^j\Big)\,dt.

Derivation hint: Apply the product rule under the integral; the first term comes from varying the metric components, the second from varying velocity.

Step 3 — Integrate by parts. On the second term, ∫gijγ˙iV˙j=−∫ddt(gijγ˙i) Vj\int g_{ij}\dot\gamma^i\dot V^j = -\int \frac{d}{dt}(g_{ij}\dot\gamma^i)\,V^j (boundary terms vanish since V(a)=V(b)=0V(a)=V(b)=0). Relabel indices and collect coefficients of VkV^k:

9
dds∣0E=−∫ab(gkjγ¨ j+(∂igkj−12∂kgij)γ˙iγ˙j)Vk dt.\quad \frac{d}{ds}\Big|_{0}E=-\int_a^b \Big( g_{kj}\ddot\gamma^{\,j} + \big(\partial_i g_{kj}-\tfrac12\partial_k g_{ij}\big)\dot\gamma^i\dot\gamma^j\Big)V^k\,dt.

Derivation hint: Expand ddt(gijγ˙i)=∂ℓgijγ˙ℓγ˙i+gijγ¨i\frac{d}{dt}(g_{ij}\dot\gamma^i)=\partial_\ell g_{ij}\dot\gamma^\ell\dot\gamma^i+g_{ij}\ddot\gamma^i, then symmetrize ∂ℓgij\partial_\ell g_{ij} in (ℓ,i)(\ell,i).

Step 4 — Fundamental lemma. Since VkV^k is arbitrary, the bracket must vanish. Multiply by gmkg^{mk} and symmetrize to recover the Christoffel symbols:

10
γ¨ m+gmk 12(∂igkj+∂jgki−∂kgij)γ˙iγ˙j=γ¨ m+Γijmγ˙iγ˙j=0.\quad \ddot\gamma^{\,m} + g^{mk}\,\tfrac12\big(\partial_i g_{kj}+\partial_j g_{ki}-\partial_k g_{ij}\big)\dot\gamma^i\dot\gamma^j = \ddot\gamma^{\,m}+\Gamma^m_{ij}\dot\gamma^i\dot\gamma^j = 0.

Derivation hint: Use γ˙iγ˙j\dot\gamma^i\dot\gamma^j symmetry to replace ∂igkj\partial_i g_{kj} by 12(∂igkj+∂jgki)\tfrac12(\partial_i g_{kj}+\partial_j g_{ki}), then apply definition [5]. □\square

Constant speed. Differentiate ∣γ˙∣2=gijγ˙iγ˙j|\dot\gamma|^2=g_{ij}\dot\gamma^i\dot\gamma^j and substitute [6]; metric compatibility forces ddt∣γ˙∣2=0\frac{d}{dt}|\dot\gamma|^2=0.

∎
mathematician · Bernhard Riemann | 1826–1866

In his 1854 Göttingen habilitation lecture Über die Hypothesen, welche der Geometrie zu Grunde liegen ("On the Hypotheses which lie at the Foundations of Geometry"), delivered before an audience including Gauss, Riemann introduced the concept of an nn-dimensional manifold equipped with a positive-definite quadratic form measuring infinitesimal length — what we now call a Riemannian metric. He proposed that the geometry of physical space is an empirical question determined by its curvature, anticipating general relativity by six decades. Riemann gave the curvature in terms of the metric and treated geodesics as the analog of straight lines, though the explicit machinery (Christoffel symbols, covariant differentiation, the tensor calculus) was supplied later by Christoffel (1869), Ricci-Curbastro, and Levi-Civita (1900).

Building Intuition

🧠 Intuition

Imagine an ant living on a curved surface that cannot see the ambient space — it has access only to a ruler (the metric) measuring distances within the surface. To walk "as straight as possible," the ant minimizes turning: at each instant it never veers left or right intrinsically. That is exactly ∇γ˙γ˙=0\nabla_{\dot\gamma}\dot\gamma = 0 — the acceleration has no component tangent to the surface; any acceleration is purely the unavoidable bending needed to stay on the surface. The Christoffel symbols are the "correction terms" that tell the ant how the coordinate grid itself curves, so that subtracting them isolates true geometric acceleration. On a sphere, great circles are geodesics; on a flat plane, ordinary straight lines. Curvature reveals itself when you release two ants walking parallel from a baseline: on a sphere they converge (positive curvature), on a saddle they diverge (negative curvature). Geodesics minimize length only locally — over long distances a great circle path is beaten by going the other way around.

Figure
Geodesics on a sphere: two great-circle arcs from p; the shorter one is the minimizer. Parallel geodesics converge under positive curvature.
p q short arc long arc converge Geodesics = locally shortest; globally one of two arcs wins.

Worked Examples

ExampleGeodesics on the 2-Sphere of Radius R

Given. The round sphere SR2S^2_R in spherical coordinates (θ,ϕ)(\theta,\phi), colatitude θ∈(0,π)\theta\in(0,\pi), longitude ϕ\phi. The induced metric is

11
g=R2 dθ2+R2sin⁡2θ dϕ2,so gθθ=R2, gϕϕ=R2sin⁡2θ.\quad g = R^2\,d\theta^2 + R^2\sin^2\theta\,d\phi^2, \quad\text{so } g_{\theta\theta}=R^2,\ g_{\phi\phi}=R^2\sin^2\theta.

Derivation hint: Embed SR2⊂R3S^2_R\subset\mathbb{R}^3 via (Rsin⁡θcos⁡ϕ,Rsin⁡θsin⁡ϕ,Rcos⁡θ)(R\sin\theta\cos\phi, R\sin\theta\sin\phi, R\cos\theta) and pull back the Euclidean metric.

Step 1 — Christoffel symbols. With gϕϕ=R2sin⁡2θg_{\phi\phi}=R^2\sin^2\theta, the nonzero symbols are

12
Γϕϕθ=−sin⁡θcos⁡θ,Γθϕϕ=Γϕθϕ=cot⁡θ.\quad \Gamma^\theta_{\phi\phi}=-\sin\theta\cos\theta, \qquad \Gamma^\phi_{\theta\phi}=\Gamma^\phi_{\phi\theta}=\cot\theta.

Derivation hint: Apply [5]; only ∂θgϕϕ=2R2sin⁡θcos⁡θ\partial_\theta g_{\phi\phi}=2R^2\sin\theta\cos\theta is nonzero, giving these two via gθθ=1/R2g^{\theta\theta}=1/R^2, gϕϕ=1/(R2sin⁡2θ)g^{\phi\phi}=1/(R^2\sin^2\theta).

Step 2 — Geodesic equations.

13
θ¨−sin⁡θcos⁡θ ϕ˙2=0,ϕ¨+2cot⁡θ θ˙ϕ˙=0.\quad \ddot\theta-\sin\theta\cos\theta\,\dot\phi^2=0, \qquad \ddot\phi+2\cot\theta\,\dot\theta\dot\phi=0.

Derivation hint: Substitute [12] into [6] for k=θk=\theta and k=ϕk=\phi.

Step 3 — Solve along the equator. Try θ(t)≡π/2\theta(t)\equiv\pi/2: then θ˙=0\dot\theta=0 and the first equation reads −sin⁡(π/2)cos⁡(π/2)ϕ˙2=0-\sin(\pi/2)\cos(\pi/2)\dot\phi^2=0, automatically satisfied; the second gives ϕ¨=0\ddot\phi=0, so ϕ=ωt\phi=\omega t. This is a great circle traversed at constant angular speed. By rotational symmetry, every great circle is a geodesic.

Step 4 — Distance from p to q. For p,qp,q subtending central angle Δσ\Delta\sigma,

 dg(p,q)=R Δσ \boxed{\,d_g(p,q)=R\,\Delta\sigma\,}

Numbers. For Earth (R=6371 kmR=6371\text{ km}), the great-circle distance between two points 45∘=π/4 rad45^\circ=\pi/4\text{ rad} apart is d=6371×π/4≈5004 kmd = 6371\times \pi/4 \approx \boxed{5004\ \text{km}}.

ExampleGeodesics in the Hyperbolic Upper Half-Plane

Given. The Poincaré upper half-plane H={(x,y):y>0}\mathbb{H}=\{(x,y):y>0\} with metric

14
g=dx2+dy2y2,gxx=gyy=1y2.\quad g=\frac{dx^2+dy^2}{y^2}, \qquad g_{xx}=g_{yy}=\frac{1}{y^2}.

Derivation hint: This is the unique (up to scale) complete metric of constant curvature −1-1 on the half-plane; conformal to the Euclidean metric with factor 1/y21/y^2.

Step 1 — Christoffel symbols. With gxx=gyy=y−2g_{xx}=g_{yy}=y^{-2}, ∂ygxx=∂ygyy=−2y−3\partial_y g_{xx}=\partial_y g_{yy}=-2y^{-3}:

15
Γxyx=Γyxx=−1y,Γxxy=1y,Γyyy=−1y.\quad \Gamma^x_{xy}=\Gamma^x_{yx}=-\tfrac1y,\quad \Gamma^y_{xx}=\tfrac1y,\quad \Gamma^y_{yy}=-\tfrac1y.

Derivation hint: Apply [5] with gxx=gyy=y2g^{xx}=g^{yy}=y^2; only yy-derivatives are nonzero.

Step 2 — Geodesic equations.

16
x¨−2yx˙y˙=0,y¨+1y(x˙2−y˙2)=0.\quad \ddot x-\frac{2}{y}\dot x\dot y=0, \qquad \ddot y+\frac{1}{y}(\dot x^2-\dot y^2)=0.

Derivation hint: Substitute [15] into [6].

Step 3 — Identify the curves. Vertical lines x=constx=\text{const} solve the system (x˙=0⇒x¨=0\dot x=0\Rightarrow\ddot x=0, and y¨=y˙2/y\ddot y=\dot y^2/y). The other geodesics are semicircles centered on the xx-axis: (x−c)2+y2=r2(x-c)^2+y^2=r^2. One checks the first integral x˙/y2=\dot x/y^2= const (a conserved momentum), which combined with constant speed yields these circles.

Step 4 — Hyperbolic distance. For two points on a common vertical line x=cx=c, heights y1<y2y_1<y_2,

 dg=∫y1y2dyy=ln⁡y2y1 \boxed{\,d_g=\int_{y_1}^{y_2}\frac{dy}{y}=\ln\frac{y_2}{y_1}\,}

Numbers. From (0,1)(0,1) to (0,e2)(0,e^2): dg=ln⁡(e2/1)=2 (hyperbolic units)d_g=\ln(e^2/1)=\boxed{2\ \text{(hyperbolic units)}}. Note the Euclidean separation is e2−1≈6.39e^2-1\approx 6.39, but hyperbolic distance compresses heights logarithmically.

Common Misconceptions

⚠️ Common mistake

"Geodesics are always the shortest paths." False globally. Geodesics are only locally length-minimizing. Past a conjugate point or cut point, a geodesic ceases to minimize: on the sphere, a great-circle arc longer than πR\pi R is no longer the shortest path between its endpoints. Geodesics are critical points of energy, not necessarily minima.

"The metric determines distances, so it must be the Euclidean dot product in some coordinates." False in general. A metric can be brought to the identity at a single point (orthonormal frame) and even to "flat to second order" (normal coordinates), but it can be flattened in a whole neighborhood iff the Riemann curvature tensor vanishes there.

"Energy and length minimizers are the same." They give the same images but length is reparametrization-invariant while energy is not — energy minimizers are automatically constant-speed (affinely parametrized) geodesics, whereas length minimizers can be parametrized arbitrarily. This is why one usually varies EE, not LL.

Counterexample

Incompleteness wrecks existence of minimizers. Take the punctured plane M=R2∖{0}M=\mathbb{R}^2\setminus\{0\} with the flat Euclidean metric. The points p=(−1,0)p=(-1,0) and q=(1,0)q=(1,0) have dg(p,q)=2d_g(p,q)=2 (the infimum), but the straight segment realizing it passes through the missing origin. No geodesic of length 2 exists in MM; minimizing sequences leave every compact set. This shows the Hopf–Rinow hypothesis (geodesic or metric completeness) is essential: on a complete Riemannian manifold any two points are joined by a minimizing geodesic, but completeness cannot be dropped.

Connections

connection

Riemannian metrics sit at the crossroads of several theories. The metric uniquely determines the Levi-Civita Connection, from which the Riemann Curvature Tensor is built; curvature in turn governs geodesic behavior via the Jacobi Equation and conjugate points. The variational viewpoint links geodesics to Calculus of Variations and the Euler–Lagrange formalism, and the energy functional connects to Morse Theory (geodesics as critical points of EE on the path space, with index = number of conjugate points). The distance function dgd_g makes (M,dg)(M,d_g) a Metric Space, bridging to metric geometry and Gromov–Hausdorff Convergence. In physics, the Lorentzian analog underlies General Relativity, where geodesics are the worldlines of free-falling particles. Optimization on the manifold of positive-definite matrices, on Lie groups, and on Stiefel/Grassmann manifolds all use Riemannian Optimization built on this same machinery.

Real-World Applications

🔬 Where this appears

Historical Development

⏳ History
- 1827 — Gauss, Disquisitiones generales circa superficies curvas. Introduces the first fundamental form (metric of a surface) and proves the Theorema Egregium: Gaussian curvature is intrinsic, computable from the metric alone.
- 1854 — Riemann, habilitation lecture at Göttingen. Defines nn-dimensional manifolds with an infinitesimal length form and the notion of curvature for higher dimensions.
- 1869 — Christoffel introduces the symbols Γijk\Gamma^k_{ij} and the covariant transformation rules central to tensor calculus.
- 1888–1901 — Ricci-Curbastro & Levi-Civita develop the absolute differential calculus (tensor analysis), systematizing covariant differentiation and the curvature tensor.
- 1900 — Bianchi publishes the differential identities (Bianchi identities) for the curvature tensor.
- 1915 — Einstein uses Riemannian (Lorentzian) geometry as the language of general relativity, with geodesics modeling free fall and the field equations relating curvature to matter.
- 1917 — Levi-Civita defines parallel transport, giving geometric meaning to the connection and clarifying the role of curvature as holonomy.
- 1926 — Élie Cartan reformulates the subject with moving frames and differential forms; introduces symmetric spaces and the modern global perspective.
- 1931 — Hopf & Rinow prove the equivalence of geodesic completeness, metric completeness, and the geodesic-connectedness of complete manifolds.

Subtleties & Limits

open-problem

Even with a complete classical theory, deep questions persist. The closed-geodesic problem: does every closed Riemannian manifold admit infinitely many distinct prime closed geodesics? This is known for surfaces (Bangert–Franks, 1992, via dynamics) and many cases, but remains open in general for higher dimensions, e.g. it is not fully settled for all metrics on SnS^n, n≥3n\ge 3. Spectral rigidity / "Can one hear the shape of a drum?" asks how much of the metric is determined by the Laplace spectrum; isospectral non-isometric manifolds exist (Milnor 1964, Gordon–Webb–Wolpert 1992), so the metric is not fully recoverable from eigenvalues. The boundary rigidity and geodesic X-ray transform problems — whether boundary distance data determine the interior metric — are central in inverse problems and only partially resolved (Pestov–Uhlmann for simple surfaces). In optimal transport and metric measure spaces, extending Ricci curvature bounds to non-smooth spaces (the Lott–Sturm–Villani synthetic theory) actively generalizes geodesic geometry beyond manifolds.

See Also

Levi-Civita Connection · Riemann Curvature Tensor · Geodesic Completeness · Exponential Map · Sectional Curvature · Lorentzian Geometry · Information Geometry · Symmetric Spaces