Curves Embedded in Surfaces

Curves embedded in surfaces present us with a fascinating interplay among three spaces: the curve itself, its ambient surface, and the overall Euclidean space. There are also three embeddings to consider, i.e. that of the surface within the Euclidean space, of the curve within the Euclidean space, and, finally and most interestingly, of the curve within the surface. Although our interest is in all three embeddings and their interplay, it is the last embedding that is new to us and will therefore receive the greatest amount of attention.
(9.1)
Each embedding has a curvature tensor associated with it. The embedding of the surface within the Euclidean space is characterized by the now-familiar curvature tensor BαβB_{\alpha\beta}. The embedding of the curve within the Euclidean space is characterized by the vector curvature tensor BΦΨ\mathbf{B}_{\Phi\Psi} introduced in the previous Chapter. Finally, the embedding of the curve within the surface is characterized by the geodesic curvature tensor bΦΨb_{\Phi\Psi}, which will be introduced in this Chapter. The geodesic curvature tensor bΦΨb_{\Phi\Psi} is entirely analogous to the curvature tensor BαβB_{\alpha\beta} and the term geodesic simply refers to the fact that bΦΨb_{\Phi\Psi} characterizes the embedding of the curve within the surface.
It is intuitively clear that the three curvature tensors are related. After all, if a surface is curved then it may not be able to accommodate straight curves, at least in some directions. Thus, high values in the curvature tensor BαβB_{\alpha\beta} may induce relatively high values in the vector curvature tensor BΦΨ\mathbf{B}_{\Phi\Psi} depending on the orientation of the curve within the surface. This insight is reflected in the formula
BΦΨ=NBαβSΦαSΨβ+nbΦΨ,(9.72)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi }^{\beta}+\mathbf{n}b_{\Phi\Psi}, \tag{9.72}
where SΦαS_{\Phi}^{\alpha} is the shift tensor of the curve with respect to the surface and n\mathbf{n} is the unit normal to the curve within the surface. This formula is the culmination of this Chapter and may be interpreted as the attribution of the curve's overall curvature characteristics (captured BΦΨ\mathbf{B}_{\Phi\Psi}) to the curvature of the surface (captured by BαβB_{\alpha\beta}) and the curvature of the curve within the surface (captured by bΦΨb_{\Phi\Psi}).
Finally, we will continue to use the tensor notation with uppercase Greek letters for curve indices, despite the fact that curves are one-dimensional objects. In addition to the advantages mentioned in the previous Chapter, the analysis presented in this Chapter will benefit from this approach in one more important way. Namely, in higher-dimensional spaces, whether artificial Euclidean or Riemannian, one can consider similar curve-surface-space embeddings where the "curve", i.e. the innermost object, is of dimension greater than one. Thanks to our present use of the tensor notation to describe embedded curves, most of our analysis will generalize easily to those situations.
In order to describe a curve embedded in a surface which is, in turn, embedded in a Euclidean space we must introduce three coordinate systems: SΦS^{\Phi} on the curve, SαS^{\alpha} on the surface, and ZiZ^{i} in the ambient space. It is somewhat of a problem that the curve and surface coordinates are denoted by the same letter SS. For instance, the symbol S1S^{1} is ambiguous as it simultaneously stands for the curve coordinate and the first surface coordinate. However, we will typically refer to the coordinates collectively either by SΦS^{\Phi} or by SαS^{\alpha} so this potential ambiguity will be avoided. As we are about to see, the only real problem momentarily arises in the vector equations of the surface and the curve.
(9.2)
Denote the dependence of the position vector R\mathbf{R} on the ambient coordinates ZiZ^{i} by the familiar function R(Z)\mathbf{R}\left( Z\right) . Let the vector equations of the surface be denoted by the function R(S)\mathbf{R}\left( S\right) . Finally, the vector equations of the curve will also be denoted by R(S)\mathbf{R}\left( S\right) since both the curve and the surface coordinates are represented by the same letter SS. This is the momentary problem we have just alluded to in the previous paragraph. However, this problem is indeed only momentary since we never discuss the vector equations themselves but, rather, their derivatives with respect to the coordinates which produce the respective covariant bases Zi\mathbf{Z}_{i}, Sα\mathbf{S}_{\alpha}, and SΦ\mathbf{S}_{\Phi}.
Those bases are given by familiar equations
Zi=R(Z)Zi          (9.3)Sα=R(S)Sα          (9.4)SΦ=R(S)SΦ.          (9.5)\begin{aligned}\mathbf{Z}_{i} & =\frac{\partial\mathbf{R}\left( Z\right) }{\partial Z^{i}}\ \ \ \ \ \ \ \ \ \ \left(9.3\right)\\\mathbf{S}_{\alpha} & =\frac{\partial\mathbf{R}\left( S\right) }{\partial S^{\alpha}}\ \ \ \ \ \ \ \ \ \ \left(9.4\right)\\\mathbf{S}_{\Phi} & =\frac{\partial\mathbf{R}\left( S\right) }{\partial S^{\Phi}}.\ \ \ \ \ \ \ \ \ \ \left(9.5\right)\end{aligned}
Note that in the last equation, the superscript of SΦS^{\Phi} in the "denominator" makes it clear R(S)\mathbf{R}\left( S\right) in the "numerator" represents the vector equation of the curve.
Each of the objects above is well familiar to us from the earlier Chapters, along with the other fundamental objects that are constructed from them, such as the metric tensors, the contravariant bases, the volume, area, and length elements, the Levi-Civita symbols, the Christoffel symbols, and the Riemann-Christoffel tensors. Furthermore, we have also established relationships between the ambient and surface objects, as well as between the ambient and curve objects. However, we are yet to discuss the corresponding relationships between the surface and the curve objects, which is precisely the task to which we will now turn our attention.
The shape of a curve embedded in a surface can be described by the relationship between the curve coordinates SΦS^{\Phi} and the surface coordinates SαS^{\alpha}. The resulting equations
Sα=Sα(S).(9.6)S^{\alpha}=S^{\alpha}\left( S\right) .\tag{9.6}
are called the geodesic equations of the curve. In and of themselves, the geodesic equations of the curve do not carry the complete information about the overall shape of the curve in the three-dimensional space. In order to get the full picture, the geodesic equations need to be combined with the equation of the surface
Zi=Zi(S).(9.7)Z^{i}=Z^{i}\left( S\right) .\tag{9.7}
When these equations are composed, the result is, of course, the equations Zi(S)Z^{i}\left( S\right) of the curve in the three-dimensional space, i.e.
Zi(S)=Zi(S(S)).(9.8)Z^{i}\left( S\right) =Z^{i}\left( S\left( S\right) \right) .\tag{9.8}
This identity can be readily differentiated with respect to SΦS^{\Phi} by applying the chain rule on the right. We have
Zi(S)SΦ=Zi(S)SαSα(S)SΦ.(9.9)\frac{\partial Z^{i}\left( S\right) }{\partial S^{\Phi}}=\frac{\partial Z^{i}\left( S\right) }{\partial S^{\alpha}}\frac{\partial S^{\alpha}\left( S\right) }{\partial S^{\Phi}}.\tag{9.9}
Of course, we are already familiar with the shift tensor
ZΦi=Zi(S)SΦ(9.10)Z_{\Phi}^{i}=\frac{\partial Z^{i}\left( S\right) }{\partial S^{\Phi}}\tag{9.10}
that characterizes the embedding of the curve in the ambient space, as well as the shift tensor
Zαi=Zi(S)Sα(9.11)Z_{\alpha}^{i}=\frac{\partial Z^{i}\left( S\right) }{\partial S^{\alpha}}\tag{9.11}
that characterizes the embedding of the surfaces in the ambient space. The only new element in the identity
Zi(S)SΦ=Zi(S)SαSα(S)SΦ.(9.9)\frac{\partial Z^{i}\left( S\right) }{\partial S^{\Phi}}=\frac{\partial Z^{i}\left( S\right) }{\partial S^{\alpha}}\frac{\partial S^{\alpha}\left( S\right) }{\partial S^{\Phi}}. \tag{9.9}
is object Sα(S)/SΦ\partial S^{\alpha}\left( S\right) /\partial S^{\Phi} associated with the embedding of the curve within the surface. This object is denoted by SΦαS_{\Phi}^{\alpha}, i.e.
SΦα=Sα(S)SΦ(9.12)S_{\Phi}^{\alpha}=\frac{\partial S^{\alpha}\left( S\right) }{\partial S^{\Phi}}\tag{9.12}
and is also referred to as a shift tensor, although it may also be referred to as the geodesic shift tensor.
The three shift tensors ZΦiZ_{\Phi}^{i}, ZαiZ_{\alpha}^{i}, and SΦαS_{\Phi }^{\alpha} are related by the identity
ZΦi=ZαiSΦα.(9.13)Z_{\Phi}^{i}=Z_{\alpha}^{i}S_{\Phi}^{\alpha}.\tag{9.13}
This identity does not have an official name, but it is well described as cascading projections. Recall from Chapter 3 that ZαiZ_{\alpha}^{i} represents projection from the ambient three-dimensional space onto the surface. Similarly, ZΦiZ_{\Phi}^{i} represents projection from the ambient space onto the curve and SΦαS_{\Phi}^{\alpha} represents projection from the tangent plane onto the curve. The identity
ZΦi=ZαiSΦα(9.13)Z_{\Phi}^{i}=Z_{\alpha}^{i}S_{\Phi}^{\alpha} \tag{9.13}
tells us that projection from the space onto the curve is equivalent to projection from the space onto the surface followed by projection onto the curve. Our derivation may have been straightforward, but this is actually a nontrivial geometric result.
Next, let us convert the relationship
ZΦi=ZαiSΦα(9.13)Z_{\Phi}^{i}=Z_{\alpha}^{i}S_{\Phi}^{\alpha} \tag{9.13}
among the shift tensors into a relationship between the curve and surface covariant bases. To this end, contract both sides of the above with the ambient basis Zi\mathbf{Z}_{i}, i.e.
ZΦiZi=ZαiSΦαZi.(9.14)Z_{\Phi}^{i}\mathbf{Z}_{i}=Z_{\alpha}^{i}S_{\Phi}^{\alpha}\mathbf{Z}_{i}.\tag{9.14}
Since ZΦiZi=SΦZ_{\Phi}^{i}\mathbf{Z}_{i}=\mathbf{S}_{\Phi} and ZαiZi=SαZ_{\alpha} ^{i}\mathbf{Z}_{i}=\mathbf{S}_{\alpha}, we arrive at the relationship
SΦ=SαSΦα(9.15)\mathbf{S}_{\Phi}=\mathbf{S}_{\alpha}S_{\Phi}^{\alpha}\tag{9.15}
which, given our accumulated experience, we could have all but predicted. And we can also predict that a differentiation of this identity will yield a relationship among the various characteristics of curvature.
Meanwhile, the equation
SΦ=SαSΦα(9.15)\mathbf{S}_{\Phi}=\mathbf{S}_{\alpha}S_{\Phi}^{\alpha} \tag{9.15}
yields the relationship between the curve and the surface metric tensors, i.e.
SΦΨ=SαβSΦαSΨβ.(9.16)S_{\Phi\Psi}=S_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi}^{\beta}.\tag{9.16}
The equivalent forms of this identity are
SΦΨ=SαΦSΨα(9.17)S_{\Phi\Psi}=S_{\alpha\Phi}S_{\Psi}^{\alpha}\tag{9.17}
and
SαΦSΨα=δΨΦ.(9.18)S_{\alpha}^{\Phi}S_{\Psi}^{\alpha}=\delta_{\Psi}^{\Phi}.\tag{9.18}
Since the embedded curve is a hypersurface with respect to the ambient surface, it has a well-defined, to within sign, unit normal n\mathbf{n} that lies in the tangent plane to the surface.
(9.19)
We will call n\mathbf{n} the geodesic normal. It is specified uniquely, to within sign, by the condition that it lies in the tangent plane, is orthogonal to textbf{SS}Φ_{\Phi}, i.e.
nSΦ=0,(9.20)\mathbf{n}\cdot\mathbf{S}_{\Phi}=0,\tag{9.20}
and that it has unit length, i.e.
nn=1.(9.21)\mathbf{n}\cdot\mathbf{n}=1.\tag{9.21}
Denoting the components of n\mathbf{n} with respect to the surface basis Sα\mathbf{S}_{\alpha} by nαn^{\alpha}, i.e.
n=nαSα,(9.22)\mathbf{n}=n^{\alpha}\mathbf{S}_{\alpha},\tag{9.22}
the above conditions read
nαSΦα=0(9.23)n_{\alpha}S_{\Phi}^{\alpha}=0\tag{9.23}
and
nαnα=1.(9.24)n_{\alpha}n^{\alpha}=1.\tag{9.24}
The explicit expression for nαn^{\alpha}, analogous to its classical counterpart
Ni=12εijkεβγZjβZkγ,(3.170)N^{i}=\frac{1}{2}\varepsilon^{ijk}\varepsilon_{\beta\gamma}Z_{j}^{\beta} Z_{k}^{\gamma}, \tag{3.170}
for the components of the surface normal N\mathbf{N}, reads
nα=εαβεΦSβΦ.(9.25)n^{\alpha}=\varepsilon^{\alpha\beta}\varepsilon_{\Phi}S_{\beta}^{\Phi}.\tag{9.25}
Naturally, it exhibits the same limitations: the resulting variant may differ in sign from the components of an a priori choice of normal. Furthermore, the quantity on the right is a tensor only with respect to the orientation-preserving coordinate changes.
The geodesic normal n\mathbf{n} can also be decomposed with respect to the ambient basis Zi\mathbf{Z}_{i}, i.e.
n=niZi.(9.26)\mathbf{n}=n^{i}\mathbf{Z}_{i}.\tag{9.26}
Since the ambient components of a vector are obtained by contracting its surface components with the shift tensor ZαiZ_{\alpha}^{i}, i.e.
Ui=UαZαi.(3.95)U^{i}=U^{\alpha}Z_{\alpha}^{i}. \tag{3.95}
we have
ni=nαZαi(9.27)n^{i}=n^{\alpha}Z_{\alpha}^{i}\tag{9.27}
and therefore
ni=εαβεΦZαiSβΦ.(9.28)n^{i}=\varepsilon^{\alpha\beta}\varepsilon_{\Phi}Z_{\alpha}^{i}S_{\beta} ^{\Phi}.\tag{9.28}
The projection formula
U=(SαU)Sα+(UN)N,(2.55)\mathbf{U}=\left( \mathbf{S}^{\alpha}\cdot\mathbf{U}\right) \mathbf{S} _{\alpha}+\left( \mathbf{U}\cdot\mathbf{N}\right) \mathbf{N,} \tag{2.55}
first derived in Chapter 2, demonstrates that a vector U\mathbf{U} on the surface, can be represented as a sum of its tangential and normal projections, (SαU)Sα\left( \mathbf{S}^{\alpha}\cdot\mathbf{U}\right) \mathbf{S}_{\alpha} and (UN)N\left( \mathbf{U}\cdot\mathbf{N}\right) \mathbf{N}. In component form, the projection formula reads
δji=ZαiZjα+NiNj.(3.159)\delta_{j}^{i}=Z_{\alpha}^{i}Z_{j}^{\alpha}+N^{i}N_{j}. \tag{3.159}
Naturally, an analogous formula exists for curves embedded in surfaces. It features the curve covariant basis textbf{SS}Φ^{\Phi} and the geodesic normal n\mathbf{n} and applies to vectors in the tangent plane. Specifically, for a vector u\mathbf{u} in tangent plane, we have
u=(SΦu)SΦ+(un)n,(9.29)\mathbf{u}=\left( \mathbf{S}^{\Phi}\cdot\mathbf{u}\right) \mathbf{S}_{\Phi }+\left( \mathbf{u}\cdot\mathbf{n}\right) \mathbf{n,}\tag{9.29}
where the terms (SΦu)SΦ\left( \mathbf{S}^{\Phi}\cdot\mathbf{u}\right) \mathbf{S}_{\Phi} and (un)n\left( \mathbf{u}\cdot\mathbf{n}\right) \mathbf{n} are interpreted as the components of u\mathbf{u} that are tangential and orthogonal to the curve. The component form of the same equation reads
δβα=SΦαSβΦ+nαnβ,(9.30)\delta_{\beta}^{\alpha}=S_{\Phi}^{\alpha}S_{\beta}^{\Phi}+n^{\alpha}n_{\beta},\tag{9.30}
which is usually applied in the form
SΦαSβΦ=δβαnαnβ(9.31)S_{\Phi}^{\alpha}S_{\beta}^{\Phi}=\delta_{\beta}^{\alpha}-n^{\alpha}n_{\beta}\tag{9.31}
whenever the combination SΦαSβΦS_{\Phi}^{\alpha}S_{\beta}^{\Phi} is encountered.
Similarly to the object ZαiZjαZ_{\alpha}^{i}Z_{j}^{\alpha} discussed in Chapter 3, SΦαSβΦS_{\Phi}^{\alpha}S_{\beta}^{\Phi} represents orthogonal projection from the tangent space of the surface to the tangent space of the curve. Also recall from Chapter 8, that the shift tensor ZΦiZ_{\Phi }^{i} is proportional to the ambient components TiT^{i} of the unit tangent T\mathbf{T}. In fact,
Ti=εΦZΦi.(9.32)T^{i}=\varepsilon^{\Phi}Z_{\Phi}^{i}.\tag{9.32}
Similarly, the shift tensor SΦαS_{\Phi}^{\alpha} is proportional to the surface components TαT^{\alpha} of T\mathbf{T}, where
T=TαSα.(9.33)\mathbf{T}=T^{\alpha}\mathbf{S}_{\alpha}.\tag{9.33}
Since the combination TαTβT_{\alpha}T^{\beta} also represents orthogonal projection from the tangent space of the surface to the tangent space of the curve, we should expect that
SαΦSΦβ=TαTβ.(9.34)S_{\alpha}^{\Phi}S_{\Phi}^{\beta}=T_{\alpha}T^{\beta}.\tag{9.34}
The simplest way to demonstrate this identity is to note that, similarly to the identity
Ti=εΦZΦi,(9.35)T^{i}=\varepsilon^{\Phi}Z_{\Phi}^{i},\tag{9.35}
we have
Tα=εΦSΦα(9.36)T^{\alpha}=\varepsilon^{\Phi}S_{\Phi}^{\alpha}\tag{9.36}
Thus, recalling that εΦεΨ=δΦΨ\varepsilon_{\Phi}\varepsilon^{\Psi}=\delta_{\Phi} ^{\Psi}, we have
TαTβ=εΦSαΦεΨSΨβ=δΦΨSαΦSΨβ=SαΦSΦβ,(9.37)T_{\alpha}T^{\beta}=\varepsilon_{\Phi}S_{\alpha}^{\Phi}\varepsilon^{\Psi }S_{\Psi}^{\beta}=\delta_{\Phi}^{\Psi}S_{\alpha}^{\Phi}S_{\Psi}^{\beta }=S_{\alpha}^{\Phi}S_{\Phi}^{\beta},\tag{9.37}
as we set out to show.
We have already defined the curve covariant derivative Φ\nabla_{\Phi} for objects and ambient and curve indices. We must now extend it to surface indices. This can be accomplished according to the well established blueprint. However, rather than give an overwhelming expression for a variant TjβΨiαΦT_{j\beta\Psi}^{i\alpha\Phi} with a fully representative collection of indices, we will give three separate expressions for TjiT_{j}^{i} , TβαT_{\beta}^{\alpha}, TΨΦT_{\Psi}^{\Phi} since the tactic of the definition is simply to indicate the proper treatment of each kind of index. We have:
ΘTΨΦ=TΨΦSΘ+ΓΘΩΦTΨΩΓΘΨΩTΩΦ          (9.38)ΘTβα=TβαSΘ+SΘγΓγωαTβωSΘγΓγβωTωα          (9.39)ΘTji=TjiSΘ+ZΘkΓkmiTjmZΘkΓkjmTmi.          (9.40)\begin{aligned}\nabla_{\Theta}T_{\Psi}^{\Phi} & =\frac{\partial T_{\Psi}^{\Phi}}{\partial S^{\Theta}}+\Gamma_{\Theta\Omega}^{\Phi}T_{\Psi}^{\Omega}-\Gamma_{\Theta\Psi }^{\Omega}T_{\Omega}^{\Phi}\ \ \ \ \ \ \ \ \ \ \left(9.38\right)\\\nabla_{\Theta}T_{\beta}^{\alpha} & =\frac{\partial T_{\beta}^{\alpha} }{\partial S^{\Theta}}+S_{\Theta}^{\gamma}\Gamma_{\gamma\omega}^{\alpha }T_{\beta}^{\omega}-S_{\Theta}^{\gamma}\Gamma_{\gamma\beta}^{\omega}T_{\omega }^{\alpha}\ \ \ \ \ \ \ \ \ \ \left(9.39\right)\\\nabla_{\Theta}T_{j}^{i} & =\frac{\partial T_{j}^{i}}{\partial S^{\Theta} }+Z_{\Theta}^{k}\Gamma_{km}^{i}T_{j}^{m}-Z_{\Theta}^{k}\Gamma_{kj}^{m} T_{m}^{i}.\ \ \ \ \ \ \ \ \ \ \left(9.40\right)\end{aligned}
By now, we have significant experience constructing differential operators that possess the tensor property along with a slew of other desirable characteristics. We will therefore limit ourselves to stating the definition and assume that the reader will fill in all of the necessary details.
For a variant TjiT_{j}^{i} defined in the ambient Euclidean space, Φ\nabla_{\Phi} satisfies the chain rule
ΦTji=ZΦkkTji,(9.41)\nabla_{\Phi}T_{j}^{i}=Z_{\Phi}^{k}\nabla_{k}T_{j}^{i},\tag{9.41}
while for a variant TβαT_{\beta}^{\alpha} defined in the wider ambient surface outside, Θ\nabla_{\Theta} satisfies the chain rule
ΦTγβ=SΦααTγβ.(9.42)\nabla_{\Phi}T_{\gamma}^{\beta}=S_{\Phi}^{\alpha}\nabla_{\alpha}T_{\gamma }^{\beta}.\tag{9.42}
As a result, Θ\nabla_{\Theta} is metrinilic with respect to all of the ambient metrics, including the covariant and contravariant bases Zi\mathbf{Z}_{i} and Zi\mathbf{Z}^{i}. However, while Θ\nabla_{\Theta} is metrinilic with respect to the scalar surface metrics, we are not able to conclude the same with respect to the surface covariant and contravariant bases Sα\mathbf{S}_{\alpha} and Sα\mathbf{S}^{\alpha} due to the lack of the metrinilic property of the surface covariant derivative α\nabla_{\alpha} with respect to those objects. Indeed, since
βSα=NBβα,(2.79)\nabla_{\beta}\mathbf{S}_{\alpha}=\mathbf{N}B_{\beta\alpha}, \tag{2.79}
we have, by the chain rule
ΦSα=SΦββSα=NSΦβBβα.(9.43)\nabla_{\Phi}\mathbf{S}_{\alpha}=S_{\Phi}^{\beta}\nabla_{\beta}\mathbf{S} _{\alpha}=\mathbf{N}S_{\Phi}^{\beta}B_{\beta\alpha}.\tag{9.43}
Let us now summarize the results of applying the curve covariant derivative to the ambient, surface, and curve bases:
ΦZi=0          (9.44)ΦSα=NSΦβBβα          (9.45)ΦSΨ=BΦΨ,          (9.46)\begin{aligned}\nabla_{\Phi}\mathbf{Z}_{i} & =\mathbf{0}\ \ \ \ \ \ \ \ \ \ \left(9.44\right)\\\nabla_{\Phi}\mathbf{S}_{\alpha} & =\mathbf{N}S_{\Phi}^{\beta}B_{\beta \alpha}\ \ \ \ \ \ \ \ \ \ \left(9.45\right)\\\nabla_{\Phi}\mathbf{S}_{\Psi} & =\mathbf{B}_{\Phi\Psi},\ \ \ \ \ \ \ \ \ \ \left(9.46\right)\end{aligned}
where BΦΨ\mathbf{B}_{\Phi\Psi} is the vector curvature tensor introduced in the preceding Chapter.
The above expressions bring us face to face with the analysis of curvature which is the task to which we now turn.

9.7.1The intuition behind geodesic curvature

The purpose of geodesic curvature is to characterize the shape of an embedded curve relative to the ambient surface. This concept is exemplified by the difference between roads that we would differentiate as straight and curved.
(9.47)
Naturally, even the "straight" road in the picture above is not straight in the absolute sense since it curves with the curvature of the Earth. However, on an intuitive level, it is "as straight as possible" given that it is constrained to follow the surface of the Earth. We might describe it as straight relative to the surface of the Earth. Such roads also correspond to our intuitive understanding of what it takes to connect two towns by the shortest possible path.
(9.48)
The road in the second picture, on the other hand, would be described as curved relative to the surface of the Earth. We will wish to design the geodesic curvature tensor so that it is low -- ideally, zero -- for the "straight" road and high for the other road. Notice, then, that it is its geodesic curvature that makes a road potentially more dangerous and thus necessitates a lower speed limit -- curvature corresponds to acceleration -- and a double yellow separator line.

9.7.2The geodesic curvature tensor

Let us proceed by analogy with the classical curvature tensor BαβB_{\alpha\beta } which characterizes a surface embedded in a Euclidean space. Recall its definition
αSβ=NBαβ.(2.79)\nabla_{\alpha}\mathbf{S}_{\beta}=\mathbf{N}B_{\alpha\beta}. \tag{2.79}
By dotting both sides with the surface normal N\mathbf{N} we obtain an explicit expression for BαβB_{\alpha\beta}, i.e.
Bαβ=NαSβ.(2.80)B_{\alpha\beta}=\mathbf{N}\cdot\nabla_{\alpha}\mathbf{S}_{\beta}. \tag{2.80}
An analogous definition for the geodesic curvature tensor bΦΨb_{\Phi\Psi} would read
ΦSΨ=nbΦΨ.(-)\nabla_{\Phi}\mathbf{S}_{\Psi}=\mathbf{n}b_{\Phi\Psi}. \tag{-}
However, this definition is flawed. Indeed, while the vector BΦΨ=ΦSΨ\mathbf{B} _{\Phi\Psi}=\nabla_{\Phi}\mathbf{S}_{\Psi} is orthogonal to the embedded curve, it is not necessarily found in the tangent plane and is therefore not necessarily proportional to the geodesic normal n\mathbf{n}. As a result, we will look to the equation
Bαβ=NαSβ(2.80)B_{\alpha\beta}=\mathbf{N}\cdot\nabla_{\alpha}\mathbf{S}_{\beta} \tag{2.80}
for the analogy and define the geodesic curvature tensor bΦΨb_{\Phi\Psi} by the equation
bΦΨ=nΦSΨ(9.49)b_{\Phi\Psi}=\mathbf{n}\cdot\nabla_{\Phi}\mathbf{S}_{\Psi}\tag{9.49}
or, equivalently
bΦΨ=nBΦΨ.(9.50)b_{\Phi\Psi}=\mathbf{n}\cdot\mathbf{B}_{\Phi\Psi}.\tag{9.50}
The tensor bΦΨb_{\Phi\Psi} thus defined is called the geodesic curvature tensor.
Observe, crucially, that the vector curvature tensor BΦΨ\mathbf{B}_{\Phi\Psi} is characteristic strictly of the curve's embedding in the overall Euclidean space and has nothing to do with the ambient surface. Thus, the dependence of the geodesic curvature tensor on the surface enters strictly through the normal n\mathbf{n}, i.e. the direction of the curve relative to the surface.
The invariant
bΦΦ,(9.51)b_{\Phi}^{\Phi},\tag{9.51}
is known as the geodesic curvature although, perhaps, it would have been better described as the geodesic mean curvature. It can be expressed as the dot product of the geodesic normal n\mathbf{n} and the curvature normal BΦΦ\mathbf{B}_{\Phi}^{\Phi}, i.e.
bΦΦ=nBΦΦ.(9.52)b_{\Phi}^{\Phi}=\mathbf{n}\cdot\mathbf{B}_{\Phi}^{\Phi}.\tag{9.52}
In other words, it is the projection of the curvature normal onto the geodesic normal. Thus, the geodesic curvature cannot exceed the absolute curvature σ\sigma, i.e.
bΦΦ<σ.(9.53)\left\vert b_{\Phi}^{\Phi}\right\vert \lt \sigma.\tag{9.53}
Recall from Chapter 8, that σ\sigma is defined as the magnitude of the curvature normal BΦΦ\mathbf{B}_{\Phi}^{\Phi}. Furthermore the geodesic curvature vanishes when the curvature normal BΦΦ\mathbf{B}_{\Phi}^{\Phi} is orthogonal to the geodesic normal n\mathbf{n}.

9.7.3A circle on the surface of a sphere

A great deal of geometric intuition about geodesic curvature and its relationship to the curvature normal can be gained by considering circles on the surface of a sphere. The following figure shows three parallel circles of different radii. For these circles, BΦΦ\mathbf{B}_{\Phi}^{\Phi} is easy to visualize as it lies in the same plane as the circle and has magnitude that is the reciprocal of the radius. Thus, for a small circle, which we would perceive as highly curved, BΦΦ\mathbf{B}_{\Phi}^{\Phi} is closely aligned with n\mathbf{n}. For a large circle BΦΦ\mathbf{B}_{\Phi}^{\Phi} is nearly orthogonal to n\mathbf{n}.
A great circle is a circle on the surface of the sphere with the greatest possible radius which, of course, equals the radius of the sphere itself. It has the property that its center is located at the center of the sphere. In other words, a great circle is the intersection of the sphere with a plane passing through its center. Consequently, its curvature normal BΦΦ\mathbf{B}_{\Phi}^{\Phi} is orthogonal to the surface of the sphere and is therefore orthogonal to n\mathbf{n}. As a result, the geodesic curvature bΦΦb_{\Phi}^{\Phi} vanishes, i.e. bΦΦ=0b_{\Phi}^{\Phi}=0. Note that it is intuitively clear that a "straight" road on the surface of the Earth follows a great circle.
(9.54)
In fact, let us determine the actual value of the geodesic curvature of a circle of radius rr on the surface of a sphere of radius RR. This configuration is illustrated in the following figure.
(9.55)
Since
bΦΦ=nBΦΦ=BΦΦcosθ(9.56)b_{\Phi}^{\Phi}=\mathbf{n}\cdot\mathbf{B}_{\Phi}^{\Phi}=-\left\vert \mathbf{B}_{\Phi}^{\Phi}\right\vert \cos\theta\tag{9.56}
while
BΦΦ=1/r     and    cosθ=R2r2R,(9.57)\left\vert \mathbf{B}_{\Phi}^{\Phi}\right\vert =1/r\text{ \ \ \ \ and \ \ \ }\cos\theta=\frac{\sqrt{R^{2}-r^{2}}}{R},\tag{9.57}
we have
bΦΦ=R2r2Rr,(9.58)b_{\Phi}^{\Phi}=-\frac{\sqrt{R^{2}-r^{2}}}{Rr},\tag{9.58}
or, equivalently,
bΦΦ=1r21R2.(9.59)b_{\Phi}^{\Phi}=-\sqrt{\frac{1}{r^{2}}-\frac{1}{R^{2}}}.\tag{9.59}
Note that this formula was derived without introducing any surface or ambient coordinates, even though the angle θ\theta is reminiscent of the longitudinal angle in spherical coordinates.
According to the above formula bΦΦb_{\Phi}^{\Phi} vanishes when r=Rr=R which confirms our intuition that a great circle is straight relative to the sphere. Also, when rr is small, we have
bΦΦ=1r+O(r)(9.60)b_{\Phi}^{\Phi}=-\frac{1}{r}+O\left( r\right)\tag{9.60}
which tells us that the geodesic curvature of a small circle essentially matches its absolute curvature and is independent of the radius of the sphere. Finally, note that the sign of the geodesic curvature depends on the a priori choice of the geodesic normal n\mathbf{n}. Had we chosen the opposite direction of n\mathbf{n}, all of the values of bΦΦb_{\Phi}^{\Phi} would have had the opposite sign.

9.7.4The coordinate space expression for the geodesic curvature tensor

A shortcoming of the definition
bΦΨ=nBΦΨ(9.50)b_{\Phi\Psi}=\mathbf{n}\cdot\mathbf{B}_{\Phi\Psi} \tag{9.50}
is that it involves a quantity -- namely, the vector curvature tensor BΦΨ\mathbf{B}_{\Phi\Psi} -- which is connected with the embedding of the curve in the overall three-dimensional Euclidean space. Meanwhile, the goal of the geodesic curvature tensor bΦΨb_{\Phi\Psi} is, of course, to characterize the shape of the curve relative to the ambient surface. We would therefore like to obtain an expression for bΦΨb_{\Phi\Psi} strictly in terms of the equations of the curve
Sα=Sα(S)(9.6)S^{\alpha}=S^{\alpha}\left( S\right) \tag{9.6}
and their derivatives.
Let us think back to the classical curvature tensor BαβB_{\alpha\beta} and its coordinate space expression
Bαβ=NiαZβi.(9.61)B_{\alpha\beta}=N_{i}\nabla_{\alpha}Z_{\beta}^{i}.\tag{9.61}
The analogous equation for the geodesic curvature tensor reads
bΦΨ=nαΦSΨα,(9.67)b_{\Phi\Psi}=n_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha}, \tag{9.67}
however, at this point it is not clear whether this definition is equivalent to the one we have already made. In the next Section, we will show that it is indeed equivalent, and will also show the surprising identity
ΦSΨα=nαbΦΨ.(9.68)\nabla_{\Phi}S_{\Psi}^{\alpha}=n^{\alpha}b_{\Phi\Psi}. \tag{9.68}
It is surprising because, as we have already discussed, ΦSΨ=nbΦΨ\nabla_{\Phi }\mathbf{S}_{\Psi}=\mathbf{n}b_{\Phi\Psi} does not hold.
The geodesic curvature bΦΦb_{\Phi}^{\Phi} is given by the equation
bΦΦ=nαΦSαΦ.(9.62)b_{\Phi}^{\Phi}=n^{\alpha}\nabla_{\Phi}S_{\alpha}^{\Phi}.\tag{9.62}
A more general form of this equation is
ΦSαΦ=bΦΦnα.(9.63)\nabla_{\Phi}S_{\alpha}^{\Phi}=b_{\Phi}^{\Phi}n_{\alpha}.\tag{9.63}
Let us go back to the analysis of the equation
BΦΨ=ΦSΨ(9.46)\mathbf{B}_{\Phi\Psi}=\nabla_{\Phi}\mathbf{S}_{\Psi} \tag{9.46}
that defines the vector curvature tensor BΦΨ\mathbf{B}_{\Phi\Psi}. Notice that we can engage the geometry of the surface by substituting the previously obtained expression for textbf{SS}Ψ_{\Psi} in terms of the surface basis Sα\mathbf{S}_{\alpha}, i.e.
SΨ=SαSΨα.(9.15)\mathbf{S}_{\Psi}=\mathbf{S}_{\alpha}S_{\Psi}^{\alpha}. \tag{9.15}
By the product rule, the identity
BΦΨ=Φ(SαSΨα)(9.46)\mathbf{B}_{\Phi\Psi}=\nabla_{\Phi}\left( \mathbf{S}_{\alpha}S_{\Psi} ^{\alpha}\right) \tag{9.46}
becomes
BΦΨ=ΦSα SΨα+SαΦSΨα.(9.64)\mathbf{B}_{\Phi\Psi}=\nabla_{\Phi}\mathbf{S}_{\alpha}~S_{\Psi}^{\alpha }+\mathbf{S}_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha}.\tag{9.64}
As we discovered earlier, ΦSα\nabla_{\Phi}\mathbf{S}_{\alpha} is given by the equation
ΦSα=NSΦβBβα.(9.44)\nabla_{\Phi}\mathbf{S}_{\alpha}=\mathbf{N}S_{\Phi}^{\beta}B_{\beta\alpha}. \tag{9.44}
Thus, we have arrived at the powerful identity
BΦΨ=NBαβSΦαSΨβ+SαΦSΨα(9.65)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi }^{\beta}+\mathbf{S}_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha}\tag{9.65}
which represents a decomposition of the vector curvature tensor BΦΨ\mathbf{B} _{\Phi\Psi} into its normal and tangential components with respect to the surface. The normal component BαβSΦαSΨβB_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi} ^{\beta} will be analyzed shortly. At the present time, let us direct our attention to the tangential components ΦSΨα\nabla_{\Phi}S_{\Psi}^{\alpha}.
As we discussed in the previous Section, we hope that the quantity ΦSΨα\nabla_{\Phi}S_{\Psi}^{\alpha} is related to the geodesic curvature tensor bΦΨb_{\Phi\Psi}. This is indeed easily seen to be the case by dotting both sides of the equation
BΦΨ=NBαβSΦαSΨβ+SαΦSΨα(9.65)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi }^{\beta}+\mathbf{S}_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha} \tag{9.65}
with the geodesic normal n\mathbf{n}, i.e.
nBΦΨ=nNBαβSΦαSΨβ+nSαΦSΨα.(9.66)\mathbf{n}\cdot\mathbf{B}_{\Phi\Psi}=\mathbf{n}\cdot\mathbf{N}B_{\alpha\beta }S_{\Phi}^{\alpha}S_{\Psi}^{\beta}+\mathbf{n}\cdot\mathbf{S}_{\alpha} \nabla_{\Phi}S_{\Psi}^{\alpha}.\tag{9.66}
The dot product nBΦΨ\mathbf{n}\cdot\mathbf{B}_{\Phi\Psi} on the left is the geodesic curvature tensor bΦΨb_{\Phi\Psi}. Additionally, nN=0\mathbf{n} \cdot\mathbf{N}=0, and nSα=nα\mathbf{n}\cdot\mathbf{S}_{\alpha}=n_{\alpha}. Thus, bΦΨb_{\Phi\Psi} is given by
bΦΨ=nαΦSΨα,(9.67)b_{\Phi\Psi}=n_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha},\tag{9.67}
as was our hope.
However, as we have already mentioned, we can demonstrate the stronger relationship
ΦSΨα=nαbΦΨ.(9.68)\nabla_{\Phi}S_{\Psi}^{\alpha}=n^{\alpha}b_{\Phi\Psi}.\tag{9.68}
Indeed, since in the equation
BΦΨ=NBαβSΦαSΨβ+SαΦSΨα,(9.65)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi }^{\beta}+\mathbf{S}_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha}, \tag{9.65}
both BΦΨ\mathbf{B}_{\Phi\Psi} and NBαβSΦαSΨβ\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha }S_{\Psi}^{\beta} are orthogonal to the curve, we can conclude that SαΦSΨα\mathbf{S}_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha} is, too, orthogonal to the curve. Since, at the same time, SαΦSΨα\mathbf{S}_{\alpha}\nabla_{\Phi}S_{\Psi }^{\alpha} lies in the tangent plane, it must be a scalar multiple of the geodesic normal n\mathbf{n}. In other words,
SαΦSΨα=ncΦΨ(9.69)\mathbf{S}_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha}=\mathbf{n}c_{\Phi\Psi}\tag{9.69}
for some cΦΨc_{\Phi\Psi}. Dotting both sides with n\mathbf{n} proves that cΦΨ=bΦΨc_{\Phi\Psi}=b_{\Phi\Psi} and thus
SαΦSΨα=nbΦΨ.(9.70)\mathbf{S}_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha}=\mathbf{n}b_{\Phi\Psi}.\tag{9.70}
From this identity, it immediately follows that
ΦSΨα=nαbΦΨ,(9.68)\nabla_{\Phi}S_{\Psi}^{\alpha}=n^{\alpha}b_{\Phi\Psi}, \tag{9.68}
as we set out to show. A more algebraic proof of this identity can be found in one of the exercises.
With this newly obtained expression for ΦSΨα\nabla_{\Phi}S_{\Psi}^{\alpha}, the equation
BΦΨ=NBαβSΦαSΨβ+SαΦSΨα(9.65)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi }^{\beta}+\mathbf{S}_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha} \tag{9.65}
becomes
BΦΨ=NBαβSΦαSΨβ+SαnαbΦΨ.(9.71)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi }^{\beta}+\mathbf{S}_{\alpha}n^{\alpha}b_{\Phi\Psi}.\tag{9.71}
Since Sαnα=n\mathbf{S}_{\alpha}n^{\alpha}=\mathbf{n}, we arrive at the beautiful identity
BΦΨ=NBαβSΦαSΨβ+nbΦΨ(9.72)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi }^{\beta}+\mathbf{n}b_{\Phi\Psi}\tag{9.72}
which captures the interplay between the three curvature tensors BΦΨ\mathbf{B}_{\Phi\Psi}, BαβB_{\alpha\beta}, and bΦΨb_{\Phi\Psi}. We will, therefore, refer to this identity as the equation of the three curvature tensors. To summarize the roles of these tensors one more time, the vector curvature tensor BΦΨ\mathbf{B}_{\Phi\Psi} characterizes the shape of the curve relative to the Euclidean space, the classical curvature tensor BαβB_{\alpha\beta} characterizes the shape of the surface relative to the Euclidean space, and the geodesic curvature tensor bΦΨb_{\Phi\Psi} characterizes the shape of the curve relative to the surface.
Raising the index Φ\Phi in the equation
BΦΨ=NBαβSΦαSΨβ+nbΦΨ(9.72)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi }^{\beta}+\mathbf{n}b_{\Phi\Psi} \tag{9.72}
and subsequently contracting on Φ\Phi and Ψ\Psi (as well as juggling α\alpha) yields the stunning formula
BΦΦ=NBβαSαΦSΦβ+nbΦΦ(9.73)\mathbf{B}_{\Phi}^{\Phi}=\mathbf{N}B_{\beta}^{\alpha}S_{\alpha}^{\Phi}S_{\Phi }^{\beta}+\mathbf{n}b_{\Phi}^{\Phi}\tag{9.73}
that relates five invariants: the curvature normal B=BΦΦ\mathbf{B} =\mathbf{B}_{\Phi}^{\Phi}, the invariant BβαSαΦSΦβB_{\beta}^{\alpha}S_{\alpha}^{\Phi }S_{\Phi}^{\beta} known as the normal curvature and discussed below, the geodesic curvature bΦΦb_{\Phi}^{\Phi}, the surface normal N\mathbf{N} and the geodesic normal n\mathbf{n}. It is fair to say that it is rare for a fundamental formula to relate as many as five invariants. Interestingly, we can increase that number to six if we replace the curvature normal B\mathbf{B} with the product of the absolute curvature σ\sigma and the principal normal P\mathbf{P}, i.e.
σP=BβαSαΦSΦβN+bΦΦn.(9.74)\sigma\mathbf{P}=B_{\beta}^{\alpha}S_{\alpha}^{\Phi}S_{\Phi}^{\beta} \mathbf{N}+b_{\Phi}^{\Phi}\mathbf{n.}\tag{9.74}

9.9.1The definition

The invariant
BβαSαΦSΦβ(9.75)B_{\beta}^{\alpha}S_{\alpha}^{\Phi}S_{\Phi}^{\beta}\tag{9.75}
in the equation
BΦΦ=NBβαSαΦSΦβ+nbΦΦ(9.73)\mathbf{B}_{\Phi}^{\Phi}=\mathbf{N}B_{\beta}^{\alpha}S_{\alpha}^{\Phi}S_{\Phi }^{\beta}+\mathbf{n}b_{\Phi}^{\Phi} \tag{9.73}
is known as the normal curvature. Since the shift tensor SΦαS_{\Phi }^{\alpha} depends only on the direction of the curve within the surface, the normal curvature is dictated solely by the curvature characteristics of the ambient surface and the direction of the curve. Also, since the sign of the curvature tensor depends on the choice of the surface normal N\mathbf{N}, the same is true for the normal curvature.
There are two alternative expressions for the normal curvature that are free of curve indices. First, with the help of the geodesic projection formula
SΦαSβΦ=δβαnαnβ,(9.31)S_{\Phi}^{\alpha}S_{\beta}^{\Phi}=\delta_{\beta}^{\alpha}-n^{\alpha}n_{\beta}, \tag{9.31}
we obtain an identity that features the mean curvature BααB_{\alpha}^{\alpha}, i.e.
BβαSαΦSΦβ=BααBβαnαnβ.(9.76)B_{\beta}^{\alpha}S_{\alpha}^{\Phi}S_{\Phi}^{\beta}=B_{\alpha}^{\alpha }-B_{\beta}^{\alpha}n_{\alpha}n^{\beta}.\tag{9.76}
Alternatively, with the help of the identity
SαΦSΦβ=TαTβ,(9.34)S_{\alpha}^{\Phi}S_{\Phi}^{\beta}=T_{\alpha}T^{\beta}, \tag{9.34}
where TαT^{\alpha} are the components of the unit tangent T\mathbf{T}, we obtain an expression for the normal curvature
BβαSαΦSΦβ=BβαTαTβ(9.77)B_{\beta}^{\alpha}S_{\alpha}^{\Phi}S_{\Phi}^{\beta}=B_{\beta}^{\alpha }T_{\alpha}T^{\beta}\tag{9.77}
that is a simple quadratic form in TαT^{\alpha}. Recall from Chapter TBD in Introduction to Tensor Calculus, that the largest and the smallest values of BβαTαTβB_{\beta}^{\alpha}T_{\alpha}T^{\beta} equal the eigenvalues of BβαB_{\beta}^{\alpha}, i.e. the principal curvatures κ1\kappa_{1} and κ2\kappa_{2}.
We can now rewrite the equation of the three curvature tensors
BΦΨ=NBαβSΦαSΨβ+nbΦΨ(9.72)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi }^{\beta}+\mathbf{n}b_{\Phi\Psi} \tag{9.72}
and the equation of the three curvatures
BΦΦ=NBβαSαΦSΦβ+nbΦΦ(9.73)\mathbf{B}_{\Phi}^{\Phi}=\mathbf{N}B_{\beta}^{\alpha}S_{\alpha}^{\Phi}S_{\Phi }^{\beta}+\mathbf{n}b_{\Phi}^{\Phi} \tag{9.73}
in the forms
BΦΨ=NBβαTαTβ+nbΦΨ(9.78)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\beta}^{\alpha}T_{\alpha}T^{\beta }+\mathbf{n}b_{\Phi\Psi}\tag{9.78}
and
BΦΦ=NBβαTαTβ+nbΦΦ,(9.79)\mathbf{B}_{\Phi}^{\Phi}=\mathbf{N}B_{\beta}^{\alpha}T_{\alpha}T^{\beta }+\mathbf{n}b_{\Phi}^{\Phi},\tag{9.79}
or, equivalently,
σP=BβαTαTβN+bΦΦn.(9.80)\sigma\mathbf{P}=B_{\beta}^{\alpha}T_{\alpha}T^{\beta}\mathbf{N}+b_{\Phi }^{\Phi}\mathbf{n.}\tag{9.80}
By the Pythagorean theorem, the absolute, normal, and geodesic curvatures are related by the equation
σ2=(BβαTαTβ)2+(bΦΦ)2.(9.81)\sigma^{2}=\left( B_{\beta}^{\alpha}T_{\alpha}T^{\beta}\right) ^{2}+\left( b_{\Phi}^{\Phi}\right) ^{2}.\tag{9.81}
The last equation is illustrated in the following diagram, which shows the unit vectors N\mathbf{N}, n\mathbf{n}, and P\mathbf{P} in the plane orthogonal to the curve. The vectors N\mathbf{N}, n\mathbf{n} are orthogonal to each other, while P\mathbf{P} may be oriented in an arbitrary way with respect to them.
(9.82)

9.9.2On the surface of a sphere

The normal curvature on the surface of a sphere is the same for all curves at all points. If the radius of the sphere is RR, then the common value of normal curvature with respect to the outward normal is 1/R-1/R. To show this, recall that on the surface of a sphere the curvature tensor is related to the metric tensor by the identity
Bβα=1Rδβα.(9.83)B_{\beta}^{\alpha}=-\frac{1}{R}\delta_{\beta}^{\alpha}.\tag{9.83}
Therefore,
BβαTαTβ=1RδβαTαTβ=1RTαTα.(9.84)B_{\beta}^{\alpha}T_{\alpha}T^{\beta}=-\frac{1}{R}\delta_{\beta}^{\alpha }T_{\alpha}T^{\beta}=-\frac{1}{R}T_{\alpha}T^{\alpha}.\tag{9.84}
Finally, since TαT^{\alpha} are the components of a unit vector, we have
BβαTαTβ=1R(9.85)B_{\beta}^{\alpha}T_{\alpha}T^{\beta}=-\frac{1}{R}\tag{9.85}
as we set out to show.
Thus, on the surface of a sphere of radius RR, the equation of three curvatures reads
σP=1RN+bΦΦn.(9.86)\sigma\mathbf{P}=-\frac{1}{R}\mathbf{N}+b_{\Phi}^{\Phi}\mathbf{n}.\tag{9.86}
and therefore the absolute curvature σ\sigma satisfies the identity
σ2=1R2+(bΦΦ)2.(9.87)\sigma^{2}=\frac{1}{R^{2}}+\left( b_{\Phi}^{\Phi}\right) ^{2}.\tag{9.87}
If the curve is a circle of radius rr, then the above equation becomes
1r2=1R2+(bΦΦ)2(9.88)\frac{1}{r^{2}}=\frac{1}{R^{2}}+\left( b_{\Phi}^{\Phi}\right) ^{2}\tag{9.88}
which is consistent with our earlier analysis of geodesic curvature of a circle on the surface of a sphere.
Suppose that a surface is cut by a plane at a point AA and consider a curve that arises at the intersection of the surface and the plane. Meusnier's theorem, named after the French mathematician Jean Baptiste Meusnier, states that the absolute curvature σ\sigma of the resulting curve at the point AA depends only on the curvature tensor BαβB_{\alpha\beta} and the orientation of the plane.
(9.89)
The analytical portion of Meusnier's theorem follows from the equations of the three curvatures
σP=BβαTαTβN+bΦΦn.(9.73)\sigma\mathbf{P}=B_{\beta}^{\alpha}T_{\alpha}T^{\beta}\mathbf{N}+b_{\Phi }^{\Phi}\mathbf{n.} \tag{9.73}
Dotting both sides with the surface normal N\mathbf{N} yields
σPN=BβαTαTβ.(9.90)\sigma\mathbf{P}\cdot\mathbf{N}=B_{\beta}^{\alpha}T_{\alpha}T^{\beta }\mathbf{.}\tag{9.90}
In other words,
σ=BβαTαTβPN.(9.91)\sigma=\frac{B_{\beta}^{\alpha}T_{\alpha}T^{\beta}}{\mathbf{P}\cdot\mathbf{N} }.\tag{9.91}
Note that the curve's unit tangent T\mathbf{T} and its principal normal P\mathbf{P} both lie in the cutting plane and therefore determine its orientation. Therefore, the right side in the above identity indeed depends only on the curvature tensor and the orientation of the surface. Also note that the right side is independent of the choice of the surface normal N\mathbf{N} since both BβαB_{\beta}^{\alpha} and N\mathbf{N} change their sign when the choice is reversed.
If γ\gamma is the angle between the principal normal P\mathbf{P} and the surface normal N\mathbf{N}, then the above formula reads
σ=BβαTαTβcosγ.(9.92)\sigma=\frac{B_{\beta}^{\alpha}T_{\alpha}T^{\beta}}{\cos\gamma}.\tag{9.92}
This is the usual form in which Meusnier's theorem appears. In particular, if the cutting plane contains the surface normal NN then the absolute curvature matches the normal curvature up to sign, i.e.
σ=BβαTαTβ.(9.93)\sigma=\left\vert B_{\beta}^{\alpha}T_{\alpha}T^{\beta}\right\vert .\tag{9.93}
Meusnier's theorem has a very striking geometric interpretation for which we need to introduce the concept of an osculating circle for a planar curve. At a given point AA on a planar curve, the osculating circle is a circle that passes through AA, and has the same tangent and curvature as the curve. Analytically speaking, the curve and the circle agree in the first and second derivatives.
The radius of the circle, i.e. the reciprocal of the absolute curvature, is referred to as the radius of curvature for the curve at the point AA. Meanwhile, the center of the circle is referred to as the center of curvature.
Now, consider a family of cutting planes passing through the point AA that share an axis containing the tangent vector T\mathbf{T}.
(9.95)
All of the resulting curves pass through the point PP and share the tangent T\mathbf{T}. Since the curves share a common tangent direction, they all have the same normal curvature BβαTαTβB_{\beta}^{\alpha}T_{\alpha}T^{\beta}. Let us choose the surface normal N\mathbf{N} so that BβαTαTβB_{\beta}^{\alpha}T_{\alpha }T^{\beta} is negative and denote BβαTαTβ-B_{\beta}^{\alpha}T_{\alpha}T^{\beta} by 1/R1/R, i.e.
BβαTαTβ=1R.(9.96)-B_{\beta}^{\alpha}T_{\alpha}T^{\beta}=\frac{1}{R}.\tag{9.96}
Then the absolute curvatures of the resulting curves are given by
σ=1Rcosγ,(9.97)\sigma=-\frac{1}{R\cos\gamma},\tag{9.97}
where γ\gamma is, once again, the angle between the principal normal P\mathbf{P} and the surface normal N\mathbf{N}.
Observe that these are the very absolute curvatures we would observe if our surface was a sphere of radius RR and, correspondingly, the curves arising as the result of the planar cuts were circles, as in the following figure.
(9.98)
The fact that the absolute curvatures of the circles satisfy the equation
σ=1Rcosγ(9.99)\sigma=-\frac{1}{R\cos\gamma}\tag{9.99}
can be seen in the following diagram which shows that orthogonal cross-section of the sphere.
  (9.100)
Thus, returning to our actual surface, if each of the curves is replaced with its osculating circle the collection of the circles will form a sphere of radius (BβαTαTβ)1\left( B_{\beta}^{\alpha}T_{\alpha}T^{\beta}\right) ^{-1}.
Recall Weingarten's equation
αN=BαβSβ(2.91)\nabla_{\alpha}\mathbf{N}=-B_{\alpha}^{\beta}\mathbf{S}_{\beta} \tag{2.91}
for the surface covariant derivative of the normal N\mathbf{N}. The component form of Weingarten's equation reads
αNi=BαβZβi.(5.2)\nabla_{\alpha}N^{i}=-B_{\alpha}^{\beta}Z_{\beta}^{i}. \tag{5.2}
Let us derive the analogous equations for the geodesic normal n\mathbf{n}.
For reasons that will become apparent shortly, we will first present the component form of Weingarten's equation which reads
Φnα=bΦΨSΨα.(9.101)\nabla_{\Phi}n^{\alpha}=-b_{\Phi}^{\Psi}S_{\Psi}^{\alpha}.\tag{9.101}
Thus, its form is entirely analogous to its classical counterpart. Since the same is true for its derivation, we will leave it as an exercise.
The vector form of Weingarten's equation, on the other hand, is fundamentally different from its classical counterpart, owing to the lack of the metrinilic property of Φ\nabla_{\Phi} with respect to the surface basis Sα\mathbf{S} _{\alpha}. Indeed, differentiating both sides of
n=nαSα,(9.22)\mathbf{n}=n^{\alpha}\mathbf{S}_{\alpha}, \tag{9.22}
we find
Φn=Φnα Sα+nαΦSα.(9.102)\nabla_{\Phi}\mathbf{n}=\nabla_{\Phi}n^{\alpha}~\mathbf{S}_{\alpha}+n^{\alpha }\nabla_{\Phi}\mathbf{S}_{\alpha}.\tag{9.102}
Recall that
ΦSα=NSΦβBβα.(9.44)\nabla_{\Phi}\mathbf{S}_{\alpha}=\mathbf{N}S_{\Phi}^{\beta}B_{\beta\alpha}. \tag{9.44}
Thus,
Φn=bΦΨSΨαSα+SΦβBβαnαN(9.103)\nabla_{\Phi}\mathbf{n}=-b_{\Phi}^{\Psi}S_{\Psi}^{\alpha}\mathbf{S}_{\alpha }+S_{\Phi}^{\beta}B_{\beta\alpha}n^{\alpha}\mathbf{N}\tag{9.103}
and, equivalently,
Φn=bΦΨSΨ+BαβnαSΦβN.(9.104)\nabla_{\Phi}\mathbf{n}=-b_{\Phi}^{\Psi}\mathbf{S}_{\Psi}+B_{\alpha\beta }n^{\alpha}S_{\Phi}^{\beta}\mathbf{N.}\tag{9.104}
The surprising aspect of the formula is that, unlike its component form, it is more complicated than its classical counterpart. Of course, this is an artifact of the curvature of the ambient surface which is captured by the term containing the curvature tensor.
Exercise 9.1From the equation
SΦ=SαSΦα,(9.15)\mathbf{S}_{\Phi}=\mathbf{S}_{\alpha}S_{\Phi}^{\alpha}, \tag{9.15}
show that
SΦΨ=SαβSΦαSΨβ(9.16)S_{\Phi\Psi}=S_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi}^{\beta} \tag{9.16}
as well as the equivalent forms
SΦΨ=SαΦSΨα(9.17)S_{\Phi\Psi}=S_{\alpha\Phi}S_{\Psi}^{\alpha} \tag{9.17}
and
SαΦSΨα=δΨΦ.(9.18)S_{\alpha}^{\Phi}S_{\Psi}^{\alpha}=\delta_{\Psi}^{\Phi}. \tag{9.18}
Exercise 9.2Show that nαn^{\alpha} given by
nα=εαβεΦSβΦ.(9.25)n^{\alpha}=\varepsilon^{\alpha\beta}\varepsilon_{\Phi}S_{\beta}^{\Phi}. \tag{9.25}
satisfies the equations nαSΦα=0n_{\alpha}S_{\Phi}^{\alpha}=0 and nαnα=1n_{\alpha }n^{\alpha}=1.
Exercise 9.3Using the approach introduced in Chapter 2 for the proof of the projection formula
U=(SαU)Sα+(UN)N,(2.55)\mathbf{U}=\left( \mathbf{S}^{\alpha}\cdot\mathbf{U}\right) \mathbf{S} _{\alpha}+\left( \mathbf{U}\cdot\mathbf{N}\right) \mathbf{N,} \tag{2.55}
prove the geodesic projection formula
u=(SΦu)SΦ+(un)n(9.29)\mathbf{u}=\left( \mathbf{S}^{\Phi}\cdot\mathbf{u}\right) \mathbf{S}_{\Phi }+\left( \mathbf{u}\cdot\mathbf{n}\right) \mathbf{n} \tag{9.29}
and establish its component form
δβα=SΦαSβΦ+nαnβ.(9.30)\delta_{\beta}^{\alpha}=S_{\Phi}^{\alpha}S_{\beta}^{\Phi}+n^{\alpha}n_{\beta}. \tag{9.30}
Exercise 9.4Prove the component form
δβα=SΦαSβΦ+nαnβ(9.30)\delta_{\beta}^{\alpha}=S_{\Phi}^{\alpha}S_{\beta}^{\Phi}+n^{\alpha}n_{\beta} \tag{9.30}
of the geodesic projection formula from the explicit expression
nα=εαβεΦSβΦ(9.25)n^{\alpha}=\varepsilon^{\alpha\beta}\varepsilon_{\Phi}S_{\beta}^{\Phi} \tag{9.25}
for the components nαn^{\alpha} of the geodesic normal.
Exercise 9.5Derive the equation
ΦSΨα=nαbΦΨ(9.68)\nabla_{\Phi}S_{\Psi}^{\alpha}=n^{\alpha}b_{\Phi\Psi} \tag{9.68}
from
BΦΨ=NBαβSΦαSΨβ+SαΦSΨα(9.65)\mathbf{B}_{\Phi\Psi}=\mathbf{N}B_{\alpha\beta}S_{\Phi}^{\alpha}S_{\Psi }^{\beta}+\mathbf{S}_{\alpha}\nabla_{\Phi}S_{\Psi}^{\alpha} \tag{9.65}
by dotting both sides with SΔ\mathbf{S}^{\Delta} to obtain
SαΔΦSΨα=0.(9.105)S_{\alpha}^{\Delta}\nabla_{\Phi}S_{\Psi}^{\alpha}=0.\tag{9.105}
Then obtain the desired result, by contracting both sides with SΔβS_{\Delta }^{\beta} and applying the geodesic projection formula.
Exercise 9.6Derive Weingarten's equation for the geodesic normal, i.e.
Φnα=bΦΨSΨα.(9.101)\nabla_{\Phi}n^{\alpha}=-b_{\Phi}^{\Psi}S_{\Psi}^{\alpha}. \tag{9.101}
Exercise 9.7Consider a curve on the surface of a sphere of radius RR referred to spherical coordinates θ,φ\theta,\varphi given by equations
θ=θ(γ)          (9.106)φ=φ(γ).          (9.107)\begin{aligned}\theta & =\theta\left( \gamma\right)\ \ \ \ \ \ \ \ \ \ \left(9.106\right)\\\varphi & =\varphi\left( \gamma\right) .\ \ \ \ \ \ \ \ \ \ \left(9.107\right)\end{aligned}
Show that its geodesic curvature is given by
bΦΦ=1R(θγγφγθγφγγ)sinθφγ(2θγ2+φγ2sin2θ)cosθ(θγ2+φγ2sin2θ)3/2.(9.108)b_{\Phi}^{\Phi}=\frac{1}{R}\frac{\left( \theta_{\gamma\gamma}\varphi_{\gamma }-\theta_{\gamma}\varphi_{\gamma\gamma}\right) \sin\theta-\varphi_{\gamma }\left( 2\theta_{\gamma}^{2}+\varphi_{\gamma}^{2}\sin^{2}\theta\right) \cos\theta}{\left( \theta_{\gamma}^{2}+\varphi_{\gamma}^{2}\sin^{2} \theta\right) ^{3/2}}.\tag{9.108}
Show that in the special case of constant θ\theta, i.e. θ(γ)=θ0\theta\left( \gamma\right) =\theta_{0}, this formula agrees with the equation
bΦΦ=R2r2Rr(9.58)b_{\Phi}^{\Phi}=-\frac{\sqrt{R^{2}-r^{2}}}{Rr} \tag{9.58}
obtained earlier.
Send feedback to Pavel