Note: Graph Geodesics Intrinsic

Steve Trettel

July 10, 2024 |

Given a real valued function $f (x, y)$ , its graph defines a surface which inherits an intrinsic geometry from the ambient $E^{3}$ . In this note, we compute the geodesic equations for this surface.

(Geodesics on a Graph)

If $γ (t) = (x (t), y (t))$ is a geodesic on the graph $(x, y, f (x, y))$ of a function $f : R^{2} \to R$ , then

\ddot{x} = - f_{x} \frac{f_{x x} {\dot{x}}^{2} + 2 f_{x y} \dot{x} \dot{y} + f_{y y} {\dot{y}}^{2}}{1 + f_{x}^{2} + f_{y}^{2}}

\ddot{y} = - f_{y} \frac{f_{x x} {\dot{x}}^{2} + 2 f_{x y} \dot{x} \dot{y} + f_{y y} {\dot{y}}^{2}}{1 + f_{x}^{2} + f_{y}^{2}}

In this note we approach the calculation intrinsically, by expanding out the defining equation $\nabla_{\dot{γ}} \dot{γ} = 0$ in coordinates. This involves some preliminary calculation of Christoffel symbols, which are recorded in generality in the note ‘2D Geometry Cheatsheet’. Here we focus on specializing this generality to the case of a surface $(x, y, f (x, y))$ in $R^{3}$ .

The Metric & Christoffel Symbols

From general calculations for a 2D Riemannian metric, we have the following form for the Christoffel symbols.

(General Christoffel Symbols in 2 Dimensions)

The Christoffel symbols of a general metric $g = (\begin{matrix} E & F F & G \end{matrix})$ are

Γ_{x x}^{x} = \frac{G E_{x} - F (2 F_{x} - E_{y})}{2 D} Γ_{x x}^{y} = \frac{E (2 F_{x} - E_{y}) - F E_{x}}{2 D}

Γ_{y y}^{x} = \frac{G (2 F_{y} - G_{x}) - F G_{y}}{2 D} Γ_{y y}^{y} = \frac{E G_{y} - F (2 F_{y} - G_{x})}{2 D}

Γ_{x y}^{x} = Γ_{y x}^{x} = \frac{G E_{y} - F G_{x}}{2 D} Γ_{x y}^{y} = Γ_{y x}^{y} = \frac{E G_{x} - F E_{y}}{2 D}

For $D = det g = E G - F^{2}$

The Metric Components

Specializing to a surface parametrized by the function $F : (x, y) \mapsto (x, y, f (x, y))$ , the metric is computed in coordinates $(x, y)$ by pullback. Specifically, let $T_{x}$ and $T_{y}$ be the $x$ - and $y -$ partial derivatives of the parameterization

T_{x} = ⟨ 1, 0, f_{x} ⟩ T_{y} = ⟨ 0, 1, f_{y} ⟩

Then the components of the metric $g = (\begin{matrix} E & F F & G \end{matrix})$ are the pairwise dot-products of these

E = T_{x} \cdot T_{x} F = T_{x} \cdot T_{y} G = T_{y} \cdot T_{y}

Below we compute these, and their derivatives (which show up in the Christoffel formulas referenced above).

\begin{aligned} E & = 1 + f_{x}^{2} & F & = f_{x} f_{y} & G & = 1 + f_{y}^{2} \\ E_{x} & = 2 f_{x} f_{x x} & F_{x} & = f_{x x} f_{y} + f_{x} f_{x y} & G_{x} & = 2 f_{y} f_{x y} \\ E_{y} & = 2 f_{x} f_{x y} & F_{y} & = f_{x y} f_{y} + f_{x} f_{y y} & G_{y} & = 2 f_{y} f_{y y} \end{aligned}

With these in hand, we compute the three independent Christoffel symbols with upper index $x$ . Under the symmetry exchanging variables $x$ and $y$ (and thus simultaneously exchanging $E$ and $G$ ) these determine all Christoffel symbols, which then determine the geodesic equations. All Christoffel symbols have the same denominator $2 D$ , so its helpful to compute this quickly here:

\begin{aligned} D & = E G - F^{2} \\ = (1 + f_{x}^{2}) (1 + f_{y}^{2}) - (f_{x} f_{y})^{2} \\ = 1 + f_{x}^{2} + f_{y}^{2} + f_{x}^{2} f_{y}^{2} - f_{x}^{2} f_{y}^{2} \\ = 1 + f_{x}^{2} + f_{y}^{2} \end{aligned}

Computing $Γ_{x x}^{x}$

The numerator of this Christoffel symbol is $G E_{x} - F (2 F_{x} - E_{y})$ . Simplifying with what we know,

\begin{aligned} G E_{x} - F (2 F_{x} - E_{y}) & = (1 + f_{y}^{2}) 2 f_{x} f_{x x} - f_{x} f_{y} (2 f_{x x} f_{y} + 2 f_{x} f_{x y} - 2 f_{x} f_{x y}) \\ = (1 + f_{y}^{2}) 2 f_{x} f_{x x} - f_{x} f_{y} (2 f_{x x} f_{y}) \\ = 2 f_{x} f_{x x} + 2 f_{x} f_{y}^{2} f_{x x} - 2 f_{x} f_{y}^{2} f_{x x} \\ = 2 f_{x} f_{x x} \end{aligned}

This gives the full Christoffel symbol upon division by $2 D$ :

Γ_{x x}^{x} = \frac{2 f_{x} f_{x x}}{2 D} = \frac{f_{x} f_{x x}}{1 + f_{x}^{2} + f_{y}^{2}}

Computing $Γ_{x y}^{x}$

The numerator of this Christoffel symbol is $G E_{y} - F G_{x}$ . Simplifying with what we know,

\begin{aligned} G E_{y} - F G_{x} & = (1 + f_{y}^{2}) 2 f_{x} f_{x y} - f_{x} f_{y} \cdot 2 f_{y} f_{x y} \\ = 2 f_{x} f_{x y} + 2 f_{x} f_{y}^{2} f_{x y} - 2 f_{x} f_{y}^{2} f_{x y} \\ = 2 f_{x} f_{x y} \end{aligned}

This gives the full Christoffel symbol upon division by $2 D$ :

Γ_{x y}^{x} = \frac{2 f_{x} f_{x y}}{2 D} = \frac{f_{x} f_{x y}}{1 + f_{x}^{2} + f_{y}^{2}}

Computing $Γ_{x x}^{y}$

The numerator of this Christoffel symbol is $G (2 F_{y} - G_{x}) - F G_{y}$ . Simplifying with what we know,

\begin{aligned} G (2 F_{y} - G_{x}) - F G_{y} & = (1 + f_{y}^{2}) (2 f_{x y} f_{y} + 2 f_{x} f_{y y} - 2 f_{y} f_{x y}) - f_{x} f_{y} \cdot 2 f_{y} f_{y y} \\ = (1 + f_{y}^{2}) (2 f_{x} f_{y y}) - 2 f_{x} f_{y}^{2} f_{y y} \\ = 2 f_{x} f_{y y} + 2 f_{x} f_{y}^{2} f_{y y} - 2 f_{x} f_{y}^{2} f_{y y} \\ = 2 f_{x} f_{y y} \end{aligned}

This gives the full Christoffel symbol upon division by $2 D$ :

Γ_{y y}^{x} = \frac{2 f_{x} f_{y y}}{2 D} = \frac{f_{x} f_{y y}}{1 + f_{x}^{2} + f_{y}^{2}}

The Geodesic Equations

We again reference the 2D geometry ‘cheatsheet’ where the geodesics of a general metric are calculated to satisfy the following ODEs:

(Geodesics)

A curve $γ (t) = (x (t), y (t))$ is a geodesic if it satisfies the following ODEs:

{\begin{matrix} \ddot{x} + {\dot{x}}^{2} Γ_{x x}^{x} + 2 \dot{x} \dot{y} Γ_{x y}^{x} + {\dot{y}}^{2} Γ_{y y}^{x} = 0 \\ \ddot{y} + {\dot{x}}^{2} Γ_{x x}^{y} + 2 \dot{x} \dot{y} Γ_{x y}^{y} + {\dot{y}}^{2} Γ_{y y}^{y} = 0 \end{matrix}}

With our calculation of the Christoffel symbols $Γ_{x x}^{x}, Γ_{x y}^{x}$ and $Γ_{y y}^{x}$ , we have everything to write down the first equation. Noting that all share the same denominator $1 + f_{x}^{2} + f_{y}^{2}$ we combine fractions to

\ddot{x} + \frac{{\dot{x}}^{2} (f_{x} f_{x x}) + 2 \dot{x} \dot{y} (f_{x} f_{x y}) + {\dot{y}}^{2} (f_{x} f_{y y})}{1 + f_{x}^{2} + f_{y}^{2}} = 0

Each term of the fraction contains a multiple of $f_{x}$ , so factoring this out and moving to the other side gives

\ddot{x} = - f_{x} \frac{f_{x x} {\dot{x}}^{2} + 2 f_{x y} \dot{x} \dot{y} + f_{y y} {\dot{y}}^{2}}{1 + f_{x}^{2} + f_{y}^{2}}

This entire situation is symmetric under the exchanges of $x$ and $y$ , and performing this swap gives the second equation. Of interest is that the entire fraction is itself invariant under this interchange, so the only changes are to the second derivative and the single-partial-derivative prefactor.

\ddot{y} = - f_{y} \frac{f_{x x} {\dot{x}}^{2} + 2 f_{x y} \dot{x} \dot{y} + f_{y y} {\dot{y}}^{2}}{1 + f_{x}^{2} + f_{y}^{2}}

◀ Note: Graph Geodesics Extrinsic

Note: Differential Forms Cheatsheet ▶

Note: Graph Geodesics Intrinsic

The Metric & Christoffel Symbols

The Metric Components

Computing Γxxx

Computing Γxyx

Computing Γxxy

The Geodesic Equations

Computing $Γ_{x x}^{x}$

Computing $Γ_{x y}^{x}$

Computing $Γ_{x x}^{y}$