Riemannian Metrics, Geodesics, and the Exponential Map

⚡ Key Takeaways

A Riemannian metric $g$ equips a smooth manifold with a smoothly varying inner product on each tangent space, turning it into a measuring instrument for lengths, angles, areas, and volumes. Geodesics are the locally length-minimizing curves — the curved-space analogue of straight lines — characterized as zero-acceleration paths of the Levi-Civita connection. The exponential map $\exp_p$ packages all geodesics through a point $p$ into a single diffeomorphism (locally) from the tangent space onto the manifold, providing the canonical "straightening" coordinates. Foundational to differential geometry (Riemann 1854, Levi-Civita 1917); ubiquitous in general relativity, optimization on manifolds, and shape analysis.

💡 The Core Insight

A Riemannian metric is infinitesimal data — an inner product at each point — yet it determines global geometric structure. The bridge from infinitesimal to global is integration: lengths integrate $\sqrt{g(\dot\gamma,\dot\gamma)}$ , the connection differentiates the metric, geodesics integrate the connection's ODE, and the exponential map integrates geodesics. Each step trades local linear-algebraic data for global differential-geometric structure. The exponential map is the keystone: it lets you transport the flat linear structure of $T_pM$ onto the curved manifold, and the failure of this transport (Jacobi fields, conjugate points, cut locus) is precisely what curvature measures.

DefinitionRiemannian Metric

Let $M$ be a smooth $n$ -manifold. A Riemannian metric $g$ is a smooth $(0,2)$ -tensor field such that for every $p \in M$ , $g_p : T_pM \times T_pM \to \mathbb{R}$ is symmetric and positive-definite. In local coordinates $(x^1,\dots,x^n)$ with coordinate basis $\partial_i = \partial/\partial x^i$ ,

1

\quad g = g_{ij}\, dx^i \otimes dx^j, \qquad g_{ij}(p) = g_p(\partial_i, \partial_j), \qquad g_{ij} = g_{ji}, \quad (g_{ij}) > 0.

Derivation hint: Expand a general $(0,2)$ -tensor in the coordinate cobasis $dx^i$ ; symmetry forces $g_{ij}=g_{ji}$ , positive-definiteness forces the matrix $(g_{ij})$ to have all positive eigenvalues at each point. The pair $(M,g)$ is a Riemannian manifold. Symbols: $g_{ij}$ is the metric matrix (smooth functions of $p$ ); $g^{ij}$ denotes its matrix inverse, satisfying $g^{ik}g_{kj} = \delta^i_j$ (Einstein summation throughout).

DefinitionLength, Energy, and Riemannian Distance

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

2

\quad L(\gamma) = \int_a^b \sqrt{g_{\gamma(t)}(\dot\gamma(t),\dot\gamma(t))}\;dt, \qquad E(\gamma) = \tfrac12\int_a^b g_{\gamma(t)}(\dot\gamma(t),\dot\gamma(t))\;dt.

Derivation hint: $L$ is the integral of the speed $\|\dot\gamma\|_g = \sqrt{g(\dot\gamma,\dot\gamma)}$ ; $E$ replaces the square root by the square, removing reparametrization invariance — by Cauchy–Schwarz $L^2 \le 2(b-a)E$ with equality iff constant speed. The Riemannian distance is

3

\quad d(p,q) = \inf\{\, L(\gamma) : \gamma \text{ piecewise-smooth}, \ \gamma(a)=p,\ \gamma(b)=q \,\}.

Derivation hint: Take the infimum of length [2] over all admissible curves joining $p$ to $q$ ; smoothness of $g$ guarantees the integrand is positive on nonconstant curves, so $d(p,q)=0 \iff p=q$ . This $d$ is a genuine metric inducing the manifold topology.

DefinitionLevi-Civita Connection

There exists a unique affine connection $\nabla$ on $(M,g)$ that is torsion-free ( $\nabla_X Y - \nabla_Y X = [X,Y]$ ) and metric-compatible ( $X\,g(Y,Z) = g(\nabla_X Y, Z) + g(Y,\nabla_X Z)$ ). Its Christoffel symbols are

4

\quad \Gamma^k_{ij} = \tfrac12\, g^{k\ell}\bigl(\partial_i g_{j\ell} + \partial_j g_{i\ell} - \partial_\ell g_{ij}\bigr).

Derivation hint: Write metric compatibility for the three cyclic permutations $X g(Y,Z)$ , $Y g(Z,X)$ , $Z g(X,Y)$ ; add the first two, subtract the third, and use torsion-freeness $[\partial_i,\partial_j]=0$ to solve algebraically for $g(\nabla_{\partial_i}\partial_j,\partial_\ell)$ — this is the <strong>Koszul formula</strong>. Here $\Gamma^k_{ij} = \Gamma^k_{ji}$ (symmetry from torsion-freeness), and $\nabla_{\partial_i}\partial_j = \Gamma^k_{ij}\partial_k$ .

DefinitionGeodesic

A curve $\gamma$ is a geodesic if its velocity is parallel along itself, $\nabla_{\dot\gamma}\dot\gamma = 0$ . In coordinates this is the geodesic equation:

5

\quad \ddot\gamma^k(t) + \Gamma^k_{ij}\bigl(\gamma(t)\bigr)\,\dot\gamma^i(t)\,\dot\gamma^j(t) = 0, \qquad k=1,\dots,n.

Derivation hint: Expand $\nabla_{\dot\gamma}\dot\gamma$ with $\dot\gamma = \dot\gamma^i\partial_i$ using the Leibniz rule and $\nabla_{\partial_i}\partial_j=\Gamma^k_{ij}\partial_k$ ; the chain rule on $\dot\gamma^k(t)$ produces $\ddot\gamma^k$ , the connection produces the quadratic $\Gamma$ -term. Geodesics have constant speed: $\tfrac{d}{dt}g(\dot\gamma,\dot\gamma) = 2g(\nabla_{\dot\gamma}\dot\gamma,\dot\gamma)=0$ .

TheoremExistence, Uniqueness, and the Exponential Map

For each $p\in M$ and $v\in T_pM$ there is a unique maximal geodesic $\gamma_v$ with $\gamma_v(0)=p$ , $\dot\gamma_v(0)=v$ , depending smoothly on $(p,v)$ . Geodesics rescale: $\gamma_{\lambda v}(t)=\gamma_v(\lambda t)$ . Define the exponential map on the open star-shaped domain $\mathcal{E}_p = \{v\in T_pM : \gamma_v \text{ defined on } [0,1]\}$ by

6

\quad \exp_p(v) = \gamma_v(1), \qquad \exp_p(tv) = \gamma_v(t).

Derivation hint: Equation [5] is a second-order ODE; recast it as a first-order system on $TM$ and apply Picard–Lindelöf for existence/uniqueness and smooth dependence on initial data. The homogeneity $\gamma_{\lambda v}(t)=\gamma_v(\lambda t)$ follows because both sides solve [5] with initial velocity $\lambda v$ . Then $d(\exp_p)_0 = \mathrm{id}_{T_pM}$ (under the canonical identification $T_0(T_pM)\cong T_pM$ ), so by the inverse function theorem $\exp_p$ is a diffeomorphism from a neighborhood of $0\in T_pM$ onto a neighborhood of $p$ .

TheoremGauss Lemma & Local Minimization

Let $v\in T_pM$ with $\exp_p$ defined at $v$ , and let $w\in T_pM\cong T_v(T_pM)$ . Then

7

\quad g_{\exp_p(v)}\bigl(d(\exp_p)_v(v),\, d(\exp_p)_v(w)\bigr) = g_p(v,w).

Derivation hint: Differentiate the energy of a one-parameter family of geodesics (a geodesic variation) and apply the first variation formula; the radial geodesics meet the geodesic spheres $\exp_p(\{|v|=r\})$ orthogonally precisely because the cross term vanishes.

Consequence: within a geodesic ball $\exp_p(B_r(0))$ , the radial geodesics are the unique shortest paths from $p$ , so $d(p,\exp_p(v)) = |v|_g$ for small $v$ .

TheoremHopf–Rinow

For a connected Riemannian manifold $(M,g)$ the following are equivalent: (i) $(M,d)$ is a complete metric space; (ii) $\exp_p$ is defined on all of $T_pM$ for some (hence every) $p$ (geodesic completeness); (iii) closed bounded subsets are compact. Any of these implies any two points are joined by a minimizing geodesic. Derivation hint: From geodesic completeness, follow the minimizing radial geodesic of length $d(p,q)$ and show via the Gauss lemma and a continuity argument that the set of $r$ for which the geodesic minimizes up to distance $r$ is open, closed, and nonempty, hence all of $[0,d(p,q)]$ .

ProofGeodesics as Critical Points of Energy

We show that solutions of [5] are exactly the critical points of the energy functional $E$ — establishing the variational characterization.

Step 1 — Set up a variation. Let $\gamma:[a,b]\to M$ with fixed endpoints and let $\gamma_s(t)$ , $s\in(-\varepsilon,\varepsilon)$ , be a smooth variation with $\gamma_0=\gamma$ and $\gamma_s(a)=\gamma(a)$ , $\gamma_s(b)=\gamma(b)$ . Write the variation field $V(t) = \partial_s\gamma_s(t)\big|_{s=0}$ , which vanishes at endpoints.

Step 2 — Differentiate energy. With $T=\partial_t\gamma_s$ and $V=\partial_s\gamma_s$ ,

\frac{d}{ds}E(\gamma_s) = \int_a^b g(\nabla_s T,\, T)\,dt.

Because the variation is smooth and $\nabla$ is torsion-free with $[\,\partial_s,\partial_t\,]=0$ pulled back, the symmetry $\nabla_s T = \nabla_t V$ holds.

Step 3 — Integrate by parts. Using metric compatibility, $g(\nabla_t V, T) = \partial_t\, g(V,T) - g(V,\nabla_t T)$ . Hence

\frac{d}{ds}\Big|_0 E(\gamma_s) = \big[g(V,T)\big]_a^b - \int_a^b g\bigl(V,\, \nabla_t \dot\gamma\bigr)\,dt.

Step 4 — Apply boundary conditions. Since $V(a)=V(b)=0$ , the boundary term drops:

\frac{d}{ds}\Big|_0 E(\gamma_s) = -\int_a^b g\bigl(V,\, \nabla_{\dot\gamma}\dot\gamma\bigr)\,dt.

Step 5 — Fundamental lemma. This vanishes for all admissible $V$ iff $\nabla_{\dot\gamma}\dot\gamma=0$ , which is [5]. $\blacksquare$ Length critical points coincide with energy critical points after constant-speed reparametrization (Cauchy–Schwarz $L^2\le 2(b-a)E$ , equality iff $\|\dot\gamma\|$ constant).

∎

mathematician · Bernhard Riemann | 1826–1866

In his 1854 Göttingen Habilitationsvortrag, "Über die Hypothesen, welche der Geometrie zu Grunde liegen" (delivered before Gauss, who was reportedly astonished), Riemann introduced the notion of an $n$ -dimensional manifold equipped with a quadratic differential form measuring infinitesimal length — the birth of Riemannian geometry. He defined curvature in arbitrary dimensions, foreseeing that geometry need not be a global a-priori but could vary point to point, governed by a metric. He died of tuberculosis at 39, having also revolutionized complex analysis, number theory, and the theory of integration. His vision of intrinsic, variable-curvature geometry became, sixty years later, the mathematical skeleton of Einstein's general relativity.

🧠 Intuition

Think of the metric $g_{ij}$ as a position-dependent ruler-and-protractor. At each point it tells you how to convert coordinate displacements into true lengths and how to read off angles. The Christoffel symbols [4] are the correction terms needed because your coordinate grid is, in general, stretched and curved — they tell a moving frame how to stay "parallel" as it slides. A geodesic is what a particle does when no force acts on it intrinsically: it keeps its velocity vector "as constant as the curved space allows" — $\nabla_{\dot\gamma}\dot\gamma=0$ says the acceleration tangent to the manifold is zero. On a sphere these are great circles; in the plane, straight lines. The exponential map is "fire an arrow and see where it lands." Stand at $p$ , choose a direction-and-speed vector $v$ in the flat tangent space, shoot a geodesic for unit time, and $\exp_p(v)$ is the landing point. Because $d(\exp_p)_0=\mathrm{id}$ , near $p$ the arrow's landing spot is faithfully proportional to the arrow's length and direction — normal coordinates make the manifold look Euclidean to first order, with all $\Gamma^k_{ij}(p)=0$ and $g_{ij}(p)=\delta_{ij}$ . Curvature is exactly the obstruction to extending this Euclidean illusion to second order.

Figure

Geodesics radiating from

p

via the exponential map; a geodesic sphere meets radial geodesics orthogonally (Gauss Lemma).

ExampleGeodesics and Distance on the Round Sphere

S^2(R)

Given: The sphere of radius $R$ with the induced (round) metric. In spherical coordinates $(\theta,\phi)$ , $\theta\in(0,\pi)$ colatitude, $\phi\in[0,2\pi)$ longitude:

g = R^2\,d\theta^2 + R^2\sin^2\theta\,d\phi^2, \qquad g_{\theta\theta}=R^2,\ \ g_{\phi\phi}=R^2\sin^2\theta.

Step 1 — Christoffel symbols. From [4], the nonzero ones are

\Gamma^\theta_{\phi\phi} = -\sin\theta\cos\theta, \qquad \Gamma^\phi_{\theta\phi}=\Gamma^\phi_{\phi\theta} = \cot\theta.

Step 2 — Geodesic equations [5]:

\ddot\theta - \sin\theta\cos\theta\,\dot\phi^2 = 0, \qquad \ddot\phi + 2\cot\theta\,\dot\theta\dot\phi = 0.

Step 3 — Solve. The equator $\theta=\pi/2$ , $\phi=t$ satisfies both ( $\dot\theta=0$ , $\cos\theta=0$ ). By the rotational symmetry of the round metric, every great circle is a geodesic; these are the only geodesics.

Step 4 — Distance. Take two points separated by central angle $\Delta\sigma$ . The minimizing geodesic is the shorter great-circle arc:

d(p,q) = R\,\Delta\sigma, \qquad \Delta\sigma\in[0,\pi].

Numbers: On Earth ( $R\approx 6371$ km), London $(51.5^\circ N, 0^\circ)$ to New York $(40.7^\circ N, 74^\circ W)$ gives $\cos\Delta\sigma = \sin\phi_1\sin\phi_2 + \cos\phi_1\cos\phi_2\cos\Delta\lambda \approx 0.563$ , so $\Delta\sigma\approx 0.974$ rad.

\boxed{d \approx 6371\times 0.974 \approx 6206\ \text{km}}

(matching the real great-circle flight distance to within rounding).

ExampleExponential Map and Normal Coordinates on the Hyperbolic Plane

\mathbb{H}^2

Given: The upper half-plane model $\mathbb{H}^2 = \{(x,y): y>0\}$ with

g = \frac{dx^2 + dy^2}{y^2}, \qquad g_{xx}=g_{yy} = y^{-2}.

This has constant curvature $K=-1$ and is geodesically complete.

Step 1 — Geodesics. They are vertical lines and semicircles centered on the $x$ -axis (the boundary at infinity). Along the vertical geodesic $x=x_0$ , parametrize by arclength: $y(s) = y_0 e^{s}$ since $ds = dy/y$ .

Step 2 — Distance along the vertical. From $p=(x_0,y_0)$ to $q=(x_0,y_1)$ :

d(p,q) = \left|\int_{y_0}^{y_1}\frac{dy}{y}\right| = \left|\ln\frac{y_1}{y_0}\right|.

Step 3 — Exponential map at $p=(0,1)$ . A unit tangent vector pointing straight up, $v = (0,a)$ in coordinates with $|v|_g = a/1 = a$ (since $g=I$ at $(0,1)$ ), generates the geodesic $\exp_p(t v) = (0, e^{ta})$ . So

\exp_{(0,1)}(0,a) = (0,e^{a}).

Step 4 — Check the Gauss-lemma length. The geodesic distance from $p$ to $\exp_p(v)$ equals $|v|_g$ :

d\bigl((0,1),(0,e^{a})\bigr) = |\ln e^{a}| = a = |v|_g.\ \checkmark

Numbers: Take $a = \ln 5 \approx 1.609$ . Then $\exp_{(0,1)}(0,\ln5) = (0,5)$ , and the Riemannian distance is

\boxed{d = \ln 5 \approx 1.609\ \text{(dimensionless metric units)}}

Note: because $K<0$ , $\exp_p$ is a global diffeomorphism (Cartan–Hadamard) — no conjugate points, no cut locus.

⚠️ Common mistake

"Geodesics minimize length." Only locally. A geodesic is a critical point of energy, not necessarily a minimum. On the sphere, the great-circle arc from the North Pole the long way around past the South Pole is a perfectly good geodesic but is not the shortest path. Globally, a geodesic minimizes only up to its first cut point. Beyond the cut locus, it either stops being shortest or stops being unique.

" $\exp_p$ is a global diffeomorphism." False in general. It is only a local diffeomorphism near $0$ . On the sphere $S^n$ , $\exp_p$ wraps all vectors of length $\pi R$ onto the single antipodal point — a catastrophic failure of injectivity. The largest radius on which $\exp_p$ is injective is the injectivity radius.

"The metric is the matrix $g_{ij}$ ." The numbers $g_{ij}$ depend on coordinates; the metric is the coordinate-free tensor. Off-diagonal $g_{ij}$ encode non-orthogonal coordinate axes, not intrinsic curvature.

"Christoffel symbols are tensor components." They are not tensorial — under coordinate change they acquire an inhomogeneous second-derivative term. This is exactly why one can always arrange $\Gamma^k_{ij}(p)=0$ at a chosen point via normal coordinates, even though curvature (built from $\partial\Gamma$ and $\Gamma\Gamma$ ) cannot be made to vanish.

CounterexampleA Geodesic That Fails to Minimize Past Its Conjugate Point

On the unit sphere $S^2$ , fix the North Pole $N$ . The geodesic shooting south along a meridian minimizes distance up to the South Pole $S$ , reached at arclength $\pi$ . The South Pole is the conjugate point of $N$ (all meridians refocus there — a one-parameter family of geodesics of equal length connects $N$ to $S$ ). Continue the geodesic past $S$ : now it is climbing back up another meridian and is strictly longer than the direct route on the other side. Concretely, a geodesic of length $\pi + 0.1$ from $N$ ends at a point whose true distance from $N$ is only $\pi - 0.1$ . Thus $\exp_N$ has a critical point (singular differential) exactly at vectors of length $\pi$ : the Jacobi field measuring geodesic spreading vanishes there, signaling the conjugate point. This is the prototype for how positive curvature focuses geodesics and destroys global minimization.

connection

Riemannian metrics sit at the hub of differential geometry. The Levi-Civita connection [4] is the metric specialization of a general Affine Connection; its non-commutativity yields the Riemann Curvature Tensor, whose traces give Ricci and scalar curvature. The second variation of energy is governed by Jacobi Fields, whose zeros mark conjugate points and detect curvature via the Jacobi equation $\nabla_t^2 J + R(J,\dot\gamma)\dot\gamma = 0$ . The exponential map provides Normal Coordinates, the cleanest local model. Globally, the interplay of curvature and topology is the subject of Comparison Geometry (Bonnet–Myers, Rauch, Toponogov) and the Cut Locus. The energy functional and its critical points generalize to Morse Theory on Path Spaces. Replacing positive-definite $g$ by signature $(-,+,+,+)$ yields a Lorentzian Manifold, the arena of General Relativity, where geodesics are free-fall worldlines. The metric also defines the canonical Riemannian Volume Form $dV_g=\sqrt{\det g_{ij}}\,dx^1\cdots dx^n$ and the Laplace–Beltrami Operator.

🔬 Where this appears

⏳ History

- 1827 — Gauss's Disquisitiones generales circa superficies curvas introduces intrinsic curvature and the Theorema Egregium: curvature depends only on the first fundamental form (the metric), not the embedding.

- 1854 — Riemann's Habilitationsvortrag generalizes Gauss to

n

dimensions, defining manifolds with a quadratic line element and sectional curvature — the conceptual founding of Riemannian geometry.

- 1869–1900 — Christoffel (1869) introduces the symbols

\Gamma^k_{ij}

and covariant differentiation; Ricci-Curbastro and Levi-Civita develop the absolute differential calculus (tensor calculus), publishing the synthesis in 1900.

- 1900s–1915 — Tensor calculus matures into the indispensable computational language; Einstein, guided by Grossmann, adopts it for gravitation.

- 1915–1916 — Einstein's general relativity casts gravity as spacetime curvature; geodesics become free-fall trajectories.

- 1917 — Levi-Civita defines parallel transport, giving the geometric meaning of the connection and clarifying the geodesic equation; independently Schouten and Hessenberg contribute.

- 1925–1928 — Élie Cartan develops the method of moving frames, the modern theory of connections, and symmetric spaces; Cartan–Hadamard theorem on complete nonpositively-curved manifolds.

- 1931 — Hopf and Rinow prove the equivalence of metric completeness, geodesic completeness, and the existence of minimizing geodesics.

- 1930s–1950s — Synge, Myers, Rauch, and Bonnet establish comparison theorems linking curvature to topology; Morse theory of geodesics (Morse, 1934) and de Rham, Hodge develop global analysis on Riemannian manifolds.

open-problem

Cut locus and injectivity radius remain analytically delicate: the cut locus can be non-smooth (a "skeleton" of the manifold), and computing or even estimating it sharply for general metrics is hard. Lower bounds on the injectivity radius in terms of curvature and volume (Cheeger's lemma) are central yet quantitatively subtle.

Convergence and collapse. Gromov–Hausdorff limits of Riemannian manifolds with bounded curvature can degenerate (collapsing), and characterizing limit spaces (e.g., Ricci limit spaces, Cheeger–Colding theory, RCD spaces) is an active frontier; whether all such limits are "rectifiable" in suitable senses is partly open.

Closed geodesics. Whether every compact Riemannian manifold carries infinitely many geometrically distinct closed geodesics is proven for the 2-sphere (Bangert–Franks–Hingston) but open in general for higher-dimensional spheres with arbitrary metrics — a celebrated unresolved problem.

Infinite dimensions. On infinite-dimensional Riemannian manifolds (shape spaces, diffeomorphism groups) the exponential map need not be a local diffeomorphism, Hopf–Rinow can fail, and metric completeness vs. geodesic completeness diverge — requiring substantial functional-analytic care.

Prescribed-curvature and isometric-embedding problems (Nash–Kuiper, the smooth Nash embedding, the regularity of isometric embeddings of negatively curved surfaces) blend geometry with hard nonlinear PDE and retain open regularity questions.

Riemannian Metrics, Geodesics, and the Exponential Map

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also