Curves Embedded in a Euclidean Space

In this Chapter, we will return to the analysis of curves. We covered this topic before in Chapter TBD of Introduction to Tensor Calculus. However, at that time, we did so without the benefit of a coordinate system in the ambient space, nor did we have the benefit of a mature tensor framework. As a result, our efforts produced a geometrically insightful theory, but not one capable of producing specific values for specific curves. By contrast, we are now in a strong position to revisit the topic of curves with the full power of the tensor framework behind us.
Our present goal is to adapt the methods that we have developed for higher-dimensional surfaces to curves. Whereas previously we treated curves by taking advantage of their one-dimensional nature, we will now, as much as possible, ignore their one-dimensional nature and treat them as a special case of higher-dimensional surfaces. Since we have already done this for planar curves in Chapter 4, we will focus on curves in a three-dimensional Euclidean space.
The key differential characteristics of a curve are the unit tangent T\mathbf{T}, the curvature normal B\mathbf{B} along with the associated concepts of the absolute curvature σ\sigma and the principal normal P\mathbf{P}, and the binormal Q\mathbf{Q} along with the associated concept of the torsion τ\tau. The vectors T\mathbf{T}, P\mathbf{P}, and Q\mathbf{Q} are related by the Frenet equations given below.
The unit tangent T\mathbf{T} is the derivative of the position vector function R(s)\mathbf{R}\left( s\right) , i.e.
T=Rs.(8.1)\mathbf{T}=\mathbf{R}_{s}.\tag{8.1}
The derivative Ts\mathbf{T}_{s} is the curvature normal B\mathbf{B}, i.e.
B=Ts.(8.2)\mathbf{B}=\mathbf{T}_{s}.\tag{8.2}
Equivalently, B\mathbf{B} is the second derivative of R\mathbf{R}, i.e.
B=Rss.(8.3)\mathbf{B}=\mathbf{R}_{ss}.\tag{8.3}
The curvature normal B\mathbf{B} is orthogonal to T\mathbf{T}. Its magnitude σ\sigma is proportional to the degree to which the curve deviates from a straight line, i.e. curvature. Thus, σ\sigma is known as the absolute curvature or, simply, curvature. The term absolute highlights the fact that σ\sigma is a nonnegative quantity. The unit vector P\mathbf{P} that points in the direction of B\mathbf{B}, i.e.
B=σP,(8.4)\mathbf{B}=\sigma\mathbf{P,}\tag{8.4}
is known as the principal normal. The plane spanned by T\mathbf{T} and P\mathbf{P} is known as the osculating plane. Note that when σ=0\sigma=0, the principal normal is not defined and, therefore, neither is the osculating plane.
The adjective principal in principal normal reminds us of the fact that P\mathbf{P} is not the only normal to a curve in three dimensions since the normal space is two-dimensional. The role of the binormal Q\mathbf{Q} is to supplement the principal normal P\mathbf{P} in providing a basis for the normal space. The binormal Q\mathbf{Q} is chosen to be the unit vector that is orthogonal to P\mathbf{P} in the normal plane. Between the two possible vectors that satisfy this condition, Q\mathbf{Q} is chosen so that the set T,P,Q\mathbf{T,P,Q} is positively oriented. With the help of the cross product, the binormal Q\mathbf{Q} is given by
Q=T×P.(8.5)\mathbf{Q}=\mathbf{T\times P.}\tag{8.5}
The set T,P,Q\mathbf{T,P,Q} is referred to as a local frame or the Frenet--Serret frame.
(8.6)
The derivative P(s)\mathbf{P}^{\prime}\left( s\right) is orthogonal to P\mathbf{P} and is therefore found in the plane spanned by T\mathbf{T} and Q\mathbf{Q}. Its projection onto T\mathbf{T} turns out to be σ-\sigma. Meanwhile its projection onto Q\mathbf{Q} is an additional curve characteristic known as the torsion τ\tau, i.e.
Ps=σT+τQ.(8.7)\mathbf{P}_{s}=-\sigma\mathbf{T}+\tau\mathbf{Q.}\tag{8.7}
Thus, we have
τ=PsQ.(8.8)\tau=\mathbf{P}_{s}\cdot\mathbf{Q.}\tag{8.8}
The geometric interpretation of torsion is the rate at which the curve leaves its instantaneous osculating plane. Finally, the derivative Q(s)\mathbf{Q} ^{\prime}\left( s\right) of the binormal is given by
Qs=τP.(8.9)\mathbf{Q}_{s}=-\tau\mathbf{P}.\tag{8.9}
As a collection, the equations
Ts=B          (8.10)Ps=σT+τQ          (8.11)Qs=τP          (8.12)\begin{aligned}\mathbf{T}_{s} & =\mathbf{B}\ \ \ \ \ \ \ \ \ \ \left(8.10\right)\\\mathbf{P}_{s} & =-\sigma\mathbf{T}+\tau\mathbf{Q}\ \ \ \ \ \ \ \ \ \ \left(8.11\right)\\\mathbf{Q}_{s} & =-\tau\mathbf{P}\ \ \ \ \ \ \ \ \ \ \left(8.12\right)\end{aligned}
are known as the Frenet equations. They have a particularly elegant appearance when written in the matrix form
[TsPsQs]=[σσττ][TPQ].(8.13)\left[ \begin{array} {c} \mathbf{T}_{s}\\ \mathbf{P}_{s}\\ \mathbf{Q}_{s} \end{array} \right] =\left[ \begin{array} {rrr} & \sigma & \\ -\sigma & & \tau\\ & -\tau & \end{array} \right] \left[ \begin{array} {c} \mathbf{T}\\ \mathbf{P}\\ \mathbf{Q} \end{array} \right] .\tag{8.13}
When interpreted as a system of ordinary differential equations, the Frenet equations tell us that the shape of a curve can be reconstructed if σ\sigma and τ\tau are given as functions σ(s)\sigma\left( s\right) and τ(s)\tau\left( s\right) of the arc length.
In the present narrative, where the unit tangent T\mathbf{T} is defined as the derivative of R(s)\mathbf{R}\left( s\right) , it is not an absolute invariant since it changes sign when the orientation of the parameterization is reversed. It can therefore be described either as an invariant with respect to orientation-preserving transformations or as a relative invariant of weight 11. Unlike T\mathbf{T}, the curvature normal B\mathbf{B} is a full-fledged invariant. However, the binormal Q\mathbf{Q}, being the cross-product of T\mathbf{T} and P\mathbf{P}, inherits the unit tangent's lack of absolute invariance. As a result, it can also be described either as an invariant with respect to orientation-preserving transformations or as a relative invariant of weight 11. Finally, the absolute curvature σ\sigma and the torsion τ\tau are absolute tensors and thus characterize the geometry of the curve independent of the parameterization.
Note that we can also take an alternative view on the choice of unit tangent T\mathbf{T}, similar to the view that we took on the choice of the unit normal N\mathbf{N} of a hypersurface. Namely, we could have a priori chosen a particular direction for the unit tangent T\mathbf{T}, which would have made it an invariant along with P\mathbf{P}, Q\mathbf{Q}, σ\sigma, and τ\tau. In that case, however, the derivative R(s)\mathbf{R}^{\prime}\left( s\right) would actually produce T\mathbf{T} only if the orientation of the parameterization is consistent with the a priori choice of T\mathbf{T}.
Now is an appropriate time to present the equations of a helix. This is so for two reasons. First, these equations will be used later in the Chapter for the calculation of the curvature and torsion. Second, these equations immediately will serve to demonstrate that an arc length parameterization is impractical and will thus motivate the development of more robust tensor techniques for analyzing curves.
We have already discussed helices in Chapter TBD of Introduction to Tensor Calculus. If you recall, helices come in one of two orientations: a right-handed helix twists in the counterclockwise direction as it goes up, while a left-handed helix twists in the clockwise direction. If the ambient space is referred to a positively-oriented Cartesian coordinate system where the zz-axis is aligned with the axis of the helix, then a right-handed helix may be given by the equations
x(γ)=rcosγ          (8.14)y(γ)=rsinγ          (8.15)z(γ)=cγ,          (8.16)\begin{aligned}x\left( \gamma\right) & =r\cos\gamma\ \ \ \ \ \ \ \ \ \ \left(8.14\right)\\y\left( \gamma\right) & =r\sin\gamma\ \ \ \ \ \ \ \ \ \ \left(8.15\right)\\z\left( \gamma\right) & =c\gamma,\ \ \ \ \ \ \ \ \ \ \left(8.16\right)\end{aligned}
while a left-handed helix may be given by
x(γ)=rcosγ          (8.17)y(γ)=rsinγ          (8.18)z(γ)=cγ.          (8.19)\begin{aligned}x\left( \gamma\right) & =-r\cos\gamma\ \ \ \ \ \ \ \ \ \ \left(8.17\right)\\y\left( \gamma\right) & =r\sin\gamma\ \ \ \ \ \ \ \ \ \ \left(8.18\right)\\z\left( \gamma\right) & =c\gamma.\ \ \ \ \ \ \ \ \ \ \left(8.19\right)\end{aligned}
The parameter rr corresponds to the radius of the helix while the value 2πc2\pi c, known as the pitch, is the vertical distance between consecutive loops.
  (8.20)
Recall that the arc length ss as a function of the parameter γ\gamma is given by the integral
s(γ)=0γx(t)2+y(t)2+z(t)2dt.(8.21)s\left( \gamma\right) =\int_{0}^{\gamma}\sqrt{x^{\prime}\left( t\right) ^{2}+y^{\prime}\left( t\right) ^{2}+z^{\prime}\left( t\right) ^{2}}dt.\tag{8.21}
For either helix above, this equation yields
s(γ)=r2+c2γ.(8.22)s\left( \gamma\right) =\sqrt{r^{2}+c^{2}}\gamma.\tag{8.22}
Therefore, if we were to re-paramaterize, say, the right-handed helix by ss, the resulting equations would read
x(s)=rcossr2+c2          (8.23)y(s)=rsinsr2+c2          (8.24)z(s)=csr2+c2.          (8.25)\begin{aligned}x\left( s\right) & =r\cos\frac{s}{\sqrt{r^{2}+c^{2}}}\ \ \ \ \ \ \ \ \ \ \left(8.23\right)\\y\left( s\right) & =r\sin\frac{s}{\sqrt{r^{2}+c^{2}}}\ \ \ \ \ \ \ \ \ \ \left(8.24\right)\\z\left( s\right) & =\frac{cs}{\sqrt{r^{2}+c^{2}}}.\ \ \ \ \ \ \ \ \ \ \left(8.25\right)\end{aligned}
Clearly, this parameterization is more cumbersome than the original one. This indicates that, for practical purposes, an arc length parameterization may not be ideal.
However, this does not even begin to illustrate the actual difficulty. Suppose that, instead of a pure circular helix, we considered a slightly more complicated shape, such as an elliptical one, given by the simple equations
x(γ)=acosγ          (8.26)y(γ)=bsinγ          (8.27)z(γ)=cγ.          (8.28)\begin{aligned}x\left( \gamma\right) & =a\cos\gamma\ \ \ \ \ \ \ \ \ \ \left(8.26\right)\\y\left( \gamma\right) & =b\sin\gamma\ \ \ \ \ \ \ \ \ \ \left(8.27\right)\\z\left( \gamma\right) & =c\gamma.\ \ \ \ \ \ \ \ \ \ \left(8.28\right)\end{aligned}
Note that for this shape, which is just as easy to describe with a parameter γ\gamma, the arc length ss is given by the integral
s(γ)=0γa2sin2t+b2cos2t+c2dt,(8.29)s\left( \gamma\right) =\int_{0}^{\gamma}\sqrt{a^{2}\sin^{2}t+b^{2}\cos ^{2}t+c^{2}}dt,\tag{8.29}
which requires special functions in order to be expressed in closed form. As a result, parameterizing an elliptical helix by arc length is, at the very least, impractical. Of course, for more complicated shapes, it is likely to be impossible.
Fortunately, the tensor framework is predicated on the arbitrariness of coordinates. Not only does it serve the practical purpose of analyzing specific shapes but also, as we have learned, of providing an even deeper geometric insight than special coordinates.
We will use the same letter LL for the fundamental objects on curves. Despite the fact that a curve is a one-dimensional object, we will use an index (whose value will always be 11) to "enumerate" the coordinate as well as the elements of all the relevant systems. Having used lowercase Latin letters for ambient indices and lowercase Greek letters for surfaces, we will use capital Greek letters for curves. Thus, for instance, the coordinate L1L^{1} will be denoted by LΦL^{\Phi}.
The benefits of using the indicial notation for one-dimensional objects cannot be overstated. Most importantly, the indicial notation will enable us to carry over all of the methods developed for higher-dimensional objects to curves. The indicial signatures will continue to guide our explorations in the same way they did for higher-dimensional objects and will continue to indicate to us the precise manner by which objects transform under coordinate changes. Thus, one of our first orders of business will be to restate our previous findings summarized above in tensor terms. So effective is the tensor framework in the analysis of curves that we will stick almost exclusively with it.
We are now in a position to march through the fundamental family of objects. The covariant basis textbf{SS}Φ_{\Phi} -- consisting of a single vector textbf{S1S_{1}}, of course -- is defined by
SΦ=R(S)SΦ.(8.30)\mathbf{S}_{\Phi}=\frac{\partial\mathbf{R}\left( S\right) }{\partial S^{\Phi}}.\tag{8.30}
Since R(S)\mathbf{R}\left( S\right) is an ordinary function of one variable S1S^{1}, it would have perhaps been more appropriate to use the ordinary derivative sign, as in
SΦ=dR(S)dSΦ.(8.31)\mathbf{S}_{\Phi}=\frac{d\mathbf{R}\left( S\right) }{dS^{\Phi}}.\tag{8.31}
However, we will stick with the partial derivative to maintain a closer analogy with the multi-dimensional case. Of course, SΦ\mathbf{S}_{\Phi} is tangential to the curve as illustrated in the following figure.
(8.32)
The covariant basis textbf{SS}Φ_{\Phi} immediately demonstrates the benefit of using an index even when an object consists of a single entry. Namely, the index tells us how the object transforms under a change of variables. Introduce the Jacobians
JΦΦ=SΦ(S)SΦ          (8.33)JΦΦ=SΦ(S)SΦ.          (8.34)\begin{aligned}J_{\Phi^{\prime}}^{\Phi} & =\frac{\partial S^{\Phi}\left( S^{\prime }\right) }{\partial S^{\Phi^{^{\prime}}}}\ \ \ \ \ \ \ \ \ \ \left(8.33\right)\\J_{\Phi}^{\Phi^{\prime}} & =\frac{\partial S^{\Phi^{\prime}}\left( S\right) }{\partial S^{\Phi}}.\ \ \ \ \ \ \ \ \ \ \left(8.34\right)\end{aligned}
Of course, each Jacobian is, once again, a single number and the two numbers are the reciprocals of each other at corresponding points. Suppose that SΦS^{\Phi^{\prime}} is an alternative parameterization of the curve and that textbf{SS}Φ_{\Phi^{\prime}} is the corresponding covariant basis. The indicial signature of (the tensor) textbf{SS}Φ_{\Phi} tells us that its relation to textbf{SS}Φ_{\Phi^{\prime}} reads
SΦ=SΦJΦΦ.(8.35)\mathbf{S}_{\Phi^{\prime}}=\mathbf{S}_{\Phi}J_{\Phi^{\prime}}^{\Phi}.\tag{8.35}
Again, in actuality, this relationship represents a simple rescaling by a number, but the indicial signature reveals what that number is.
The metric tensor SΦΨS_{\Phi\Psi} is defined by the familiar formula
SΦΨ=SΦSΨ.(8.36)S_{\Phi\Psi}=\mathbf{S}_{\Phi}\cdot\mathbf{S}_{\Psi}.\tag{8.36}
It has a single element S11S_{11} and thus corresponds to a 1×11\times1 matrix. Once again, the benefit of enumerating the single element by a pair of indices is knowing the exact manner in which SΦΨS_{\Phi\Psi} transforms under a change of variables. Namely, we know that SΦΨS_{\Phi\Psi} is related to SΦΨS_{\Phi ^{\prime}\Psi^{\prime}} by the identity
SΦΨ=SΦΨJΦΦJΨΨ.(8.37)S_{\Phi^{\prime}\Psi^{\prime}}=S_{\Phi\Psi}J_{\Phi^{\prime}}^{\Phi} J_{\Psi^{\prime}}^{\Psi}.\tag{8.37}
Again, it is a single rescaling, but the indicial signature tells us what that rescaling is.
The contravariant metric tensor SΦΨS^{\Phi\Psi} is the "matrix" inverse of SΦΨS_{\Phi\Psi}, i.e.
SΦΨSΨΘ=δΦΘ.(8.38)S_{\Phi\Psi}S^{\Psi\Theta}=\delta_{\Phi}^{\Theta}.\tag{8.38}
In actuality, this identity tells us that the only element of SΨΘS^{\Psi\Theta} is the reciprocal of the only element of SΦΨS_{\Phi\Psi}. That element transforms according to the rule
SΦΨ=SΦΨJΦΦJΨΨ.(8.39)S^{\Phi^{\prime}\Psi^{\prime}}=S^{\Phi\Psi}J_{\Phi}^{\Phi^{\prime}}J_{\Psi }^{\Psi^{\prime}}.\tag{8.39}
The contravariant basis textbf{SS}Φ^{\Phi} is given by
SΦ=SΦΨSΨ.(8.40)\mathbf{S}^{\Phi}=S^{\Phi\Psi}\mathbf{S}_{\Psi}.\tag{8.40}
The vector textbf{SS}Φ^{\Phi} points in the same direction as textbf{SS }Φ_{\Phi} and its length is the reciprocal of the length of textbf{SS }Φ_{\Phi}.
Finally, the length element S\sqrt{S} is the square root of the determinant of the covariant metric tensor SΦΨS_{\Phi\Psi}. Since SΦΨS_{\Phi\Psi} corresponds to a 1×11\times1 matrix, its determinant equals the value of its sole element.
You may be surprised that the Levi-Civita symbols, which we typically associate with skew-symmetry, continue to be highly relevant, even though skew-symmetry is not possible in a one-dimensional space. The Levi-Civita symbols εΦ\varepsilon_{\Phi} and εΦ\varepsilon^{\Phi} are defined by the equations
εΦ=S          (8.41)εΦ=1S.          (8.42)\begin{aligned}\varepsilon_{\Phi} & =\sqrt{S}\ \ \ \ \ \ \ \ \ \ \left(8.41\right)\\\varepsilon^{\Phi} & =\frac{1}{\sqrt{S}}.\ \ \ \ \ \ \ \ \ \ \left(8.42\right)\end{aligned}
As in the higher-dimensional case, the Levi-Civita symbols are tensors with respect to orientation-preserving transformations. They have the additional property that they identically equal 11 when SΦS^{\Phi} is an arc length parameterization. This property can also be used as their definition: εΦ\varepsilon_{\Phi} is the (almost) unique covariant tensor that equals 11 for an arc length parameterization and εΦ\varepsilon^{\Phi} is the (almost) unique contravariant tensor that equals 11 for an arc length parameterization. (The word almost is meant to remind us that εΦ\varepsilon_{\Phi} and εΦ\varepsilon^{\Phi} do not transform by the tensor rule when the orientation of the parameterization changes.) Also note the interesting circumstance that the line element S\sqrt{S} and the covariant Levi-Civita symbol εΦ\varepsilon_{\Phi} have the exact same value. Nevertheless, while S\sqrt{S} (treated as a variant of order 00), is a relative tensor of weight 11, εΦ\varepsilon_{\Phi} (treated as a variant of order 11) is an absolute tensor.
The Christoffel symbol ΓΦΨΘ\Gamma_{\Phi\Psi}^{\Theta} is given by the dot product
ΓΦΨΘ=SΘSΦSΨ(8.43)\Gamma_{\Phi\Psi}^{\Theta}=\mathbf{S}^{\Theta}\cdot\frac{\partial \mathbf{S}_{\Phi}}{\partial S^{\Psi}}\tag{8.43}
In terms of the metric tensors, the Christoffel symbol is given by the familiar identity
ΓΦΨΘ=12SΘΩ(SΩΦSΨ+SΩΨSΦSΦΨSΩ).(8.44)\Gamma_{\Phi\Psi}^{\Theta}=\frac{1}{2}S^{\Theta\Omega}\left( \frac{\partial S_{\Omega\Phi}}{\partial S^{\Psi}}+\frac{\partial S_{\Omega\Psi}}{\partial S^{\Phi}}-\frac{\partial S_{\Phi\Psi}}{\partial S^{\Omega}}\right) .\tag{8.44}
By analogy with equation
Γαβγ=ΓijkZαiZβjZkγ+ZiγZαiSβ,(3.125)\Gamma_{\alpha\beta}^{\gamma}=\Gamma_{ij}^{k}Z_{\alpha}^{i}Z_{\beta}^{j} Z_{k}^{\gamma}+Z_{i}^{\gamma}\frac{\partial Z_{\alpha}^{i}}{\partial S^{\beta }}, \tag{3.125}
the Christoffel symbol ΓΦΨΘ\Gamma_{\Phi\Psi}^{\Theta} is related to its ambient counterpart Γijk\Gamma_{ij}^{k} by the identity
ΓΦΨΘ=ΓijkZΦiZΨjZkΘ+ZiΘZΦiSΨ.(8.45)\Gamma_{\Phi\Psi}^{\Theta}=\Gamma_{ij}^{k}Z_{\Phi}^{i}Z_{\Psi}^{j} Z_{k}^{\Theta}+Z_{i}^{\Theta}\frac{\partial Z_{\Phi}^{i}}{\partial S^{\Psi}}.\tag{8.45}
As was the case for all the preceding objects, the Christoffel symbol ΓΦΨΘ\Gamma_{\Phi\Psi}^{\Theta} has the single element Γ111\Gamma_{11}^{1}. Its rich indicial signature, however, serves the same valuable purpose as before and that is to indicate the precise manner in which the single element of the Christoffel symbol transforms under a change of coordinates. Specifically,
ΓΦΨΘ=ΓΦΨΘJΘΘJΦΦJΨΨ+JΘΘJΦΨΘ,(8.46)\Gamma_{\Phi^{\prime}\Psi^{\prime}}^{\Theta^{\prime}}=\Gamma_{\Phi\Psi }^{\Theta}J_{\Theta}^{\Theta^{\prime}}J_{\Phi^{\prime}}^{\Phi}J_{\Psi^{\prime }}^{\Psi}+J_{\Theta}^{\Theta^{\prime}}J_{\Phi^{\prime}\Psi^{\prime}}^{\Theta},\tag{8.46}
where JΦΨΘJ_{\Phi^{\prime}\Psi^{\prime}}^{\Theta} is given by
JΦΨΘ=2SΘ(S)SΦSΨ.(8.47)J_{\Phi^{\prime}\Psi^{\prime}}^{\Theta}=\frac{\partial^{2}S^{\Theta}\left( S^{\prime}\right) }{\partial S^{\Phi^{\prime}}\partial S^{\Psi^{\prime}}}.\tag{8.47}
The Riemann-Christoffel tensor RΘΦΨΩR_{\cdot\Theta\Phi\Psi}^{\Omega} is given by
RΘΦΨΔ=ΓΨΘΔSΦΓΦΘΔSΨ+ΓΦΩΔΓΨΘΩΓΨΩΔΓΦΘΩ.(8.48)R_{\cdot\Theta\Phi\Psi}^{\Delta}=\frac{\partial\Gamma_{\Psi\Theta}^{\Delta} }{\partial S^{\Phi}}-\frac{\partial\Gamma_{\Phi\Theta}^{\Delta}}{\partial S^{\Psi}}+\Gamma_{\Phi\Omega}^{\Delta}\Gamma_{\Psi\Theta}^{\Omega} -\Gamma_{\Psi\Omega}^{\Delta}\Gamma_{\Phi\Theta}^{\Omega}.\tag{8.48}
However, owing to the one-dimensional nature of curves, the Riemann-Christoffel tensor vanishes identically, i.e.
RΘΦΨΩ=0.(8.49)R_{\cdot\Theta\Phi\Psi}^{\Omega}=0.\tag{8.49}
This can be demonstrated in a number of ways, but the most fundamental way to show this is to note that RΘΦΨΔR_{\cdot\Theta\Phi\Psi}^{\Delta} vanishes under an arc-length parameterization and, thus, being a tensor, under all parameterizations. Of course, the fact that RΘΦΨΔR_{\cdot\Theta\Phi\Psi}^{\Delta} vanishes does not surprise us at all. After all, any curve can be isometrically transformed into a straight line.
The covariant derivative Θ\nabla_{\Theta} applied to a variant TΨΦT_{\Psi }^{\Phi} with a representative collection of curve indices is entirely analogous to the surface covariant derivative, i.e.
ΘUΨΦ=UΨΦSΘ+ΓΘΩΦUΨΩΓΘΨΩUΩΦ.(8.50)\nabla_{\Theta}U_{\Psi}^{\Phi}=\frac{\partial U_{\Psi}^{\Phi}}{\partial S^{\Theta}}+\Gamma_{\Theta\Omega}^{\Phi}U_{\Psi}^{\Omega}-\Gamma_{\Theta\Psi }^{\Omega}U_{\Omega}^{\Phi}.\tag{8.50}
Naturally, the covariant derivative Θ\nabla_{\Theta} possesses all of the familiar properties we expect of it. We have already enumerated those properties on a number of occasions and will not repeat them here.
As was the case for the surface covariant derivative γ\nabla_{\gamma}, the curve covariant derivative Θ\nabla_{\Theta} is not metrinilic with respect to the basis SΦ\mathbf{S}_{\Phi}. This is, of course, due to curvature. Nevertheless, as before, ΦSΨ\nabla_{\Phi}\mathbf{S}_{\Psi} is orthogonal to the tangent space. However, the familiar jump to the would-be curvature tensor BΦΨB_{\Phi\Psi} via the identity
ΦSΨ=NBΦΨ(-)\nabla_{\Phi}\mathbf{S}_{\Psi}=\mathbf{N}B_{\Phi\Psi} \tag{-}
is not possible since there is no such thing as a well-defined normal N\mathbf{N}. Indeed, unlike two-dimensional surfaces, curves embedded in a three-dimensional space are not hypersurfaces as their dimension trails that of the ambient space by 22. As a result, the normal space at each point is two-dimensional, as illustrated in the following figure, and there is no a priori normal direction.
(8.51)
Not all is lost, however. Instead of a scalar curvature tensor BΦΨB_{\Phi\Psi}, we will introduce the vector curvature tensor BΦΨ\mathbf{B}_{\Phi\Psi} given by
BΦΨ=ΦSΨ(8.52)\mathbf{B}_{\Phi\Psi}=\nabla_{\Phi}\mathbf{S}_{\Psi}\tag{8.52}
or, equivalently,
BΦΨ=ΦΨR.(8.53)\mathbf{B}_{\Phi\Psi}=\nabla_{\Phi}\nabla_{\Psi}\mathbf{R.}\tag{8.53}
As we mentioned above, BΦΨ\mathbf{B}_{\Phi\Psi} is orthogonal to the one-dimensional tangent space, i.e.
SΘBΦΨ=0.(8.54)\mathbf{S}^{\Theta}\cdot\mathbf{B}_{\Phi\Psi}=0.\tag{8.54}
Of particular interest is the invariant BΦΦ\mathbf{B}_{\Phi}^{\Phi}, analogous to mean curvature, given by
BΦΦ=ΦSΦ(8.55)\mathbf{B}_{\Phi}^{\Phi}=\nabla_{\Phi}\mathbf{S}^{\Phi}\tag{8.55}
or, equivalently
BΦΦ=ΦΦR.(8.56)\mathbf{B}_{\Phi}^{\Phi}=\nabla_{\Phi}\nabla^{\Phi}\mathbf{R.}\tag{8.56}
The invariant BΦΦ\mathbf{B}_{\Phi}^{\Phi} may be referred to as the vector mean curvature.
Not surprisingly, the invariant BΦΦ\mathbf{B}_{\Phi}^{\Phi} coincides with the curvature normal B\mathbf{B}. This can be easily seen by observing that, on the one hand, BΦΦ\mathbf{B}_{\Phi}^{\Phi} is an invariant and, on the other, it coincides with B\mathbf{B} under an arc length parameterization. Indeed, since the Christoffel symbols vanish and SΦΨ=1S^{\Phi\Psi}=1, we have
ΦΦR=d2Rds2=B.(8.57)\nabla_{\Phi}\nabla^{\Phi}\mathbf{R=}\frac{d^{2}\mathbf{R}}{ds^{2}} =\mathbf{B}.\tag{8.57}
This observation also provides an alternative proof of the fact that the curvature normal B\mathbf{B} is an unqualified invariant, unlike the unit normal T\mathbf{T} which is an invariant only with respect to orientation-preserving coordinate changes.
Since BΦΦ\mathbf{B}_{\Phi}^{\Phi} coincides with the curvature normal B\mathbf{B}, the absolute curvature σ\sigma, being the magnitude of B\mathbf{B}, is given by the equation
σ=BΦΦBΨΨ.(8.58)\sigma=\sqrt{\mathbf{B}_{\Phi}^{\Phi}\cdot\mathbf{B}_{\Psi}^{\Psi}}.\tag{8.58}
Let us now conclude this Section with an interesting observation. While the curvature normal B\mathbf{B} emerged naturally (as the vector mean curvature BΦΦ\mathbf{B}_{\Phi}^{\Phi}) as we followed the blueprint established for surfaces, the unit tangent T\mathbf{T} has eluded our analysis. Of course, since T\mathbf{T} is the unit vector that points in the same direction as SΦ\mathbf{S}_{\Phi}, we could have expressed it as
T=SΦSΦ(8.59)\mathbf{T}=\frac{\mathbf{S}_{\Phi}}{\left\vert \mathbf{S}_{\Phi}\right\vert }\tag{8.59}
or as
T=S1S1(8.60)\mathbf{T}=\frac{\mathbf{S}_{1}}{\left\vert \mathbf{S}_{1}\right\vert }\tag{8.60}
or, since S1=S\left\vert \mathbf{S}_{1}\right\vert =\sqrt{S}, as
T=SΦS.(8.61)\mathbf{T}=\frac{\mathbf{S}_{\Phi}}{\sqrt{S}}.\tag{8.61}
However, none of these expressions live up to the tensor standard as each of them features a hanging index in one form or another. Fortunately, a remedy is available and it is presented in the next Section.
So what is the "tensor" way of converting the covariant basis vector SΦ\mathbf{S}_{\Phi} into the unit tangent T\mathbf{T}? The covariant basis SΦ\mathbf{S}_{\Phi} is a first-order tensor. On a pure logistical level, it could be converted into an invariant if only there existed a universal first-order contravariant tensor with which SΦ\mathbf{S}_{\Phi} could be contracted. Of course, thanks to the one-dimensional nature of curves, such a tensor does exist and it is the Levi-Civita symbol εΦ\varepsilon^{\Phi}. We should once again add the caveat that εΦ\varepsilon^{\Phi} is only a tensor with respect to orientation-preserving coordinate changes but, then again, so is T\mathbf{T}. This caveat applies to most of the statements in this Section and we will therefore simply keep this caveat in mind while using the terms tensor and invariant in the unqualified sense.
And so, consider the combination
εΦSΦ.(8.62)\varepsilon^{\Phi}\mathbf{S}_{\Phi}.\tag{8.62}
On the one hand, it is an invariant since εΦ\varepsilon^{\Phi} is a contravariant tensor while SΦ\mathbf{S}_{\Phi} is a covariant tensor. At the same time, under an arc-length parameterization, εΦ=1\varepsilon^{\Phi}=1 and SΦ\mathbf{S}_{\Phi} is unit length. Consequently, being a vector of length 11 that points in the direction of SΦ\mathbf{S}_{\Phi}, the vector εΦSΦ\varepsilon ^{\Phi}\mathbf{S}_{\Phi} coincides with the unit tangent T\mathbf{T}. Therefore, it coincides with T\mathbf{T} under all parameterizations. In other words, we always have
T=εΦSΦ,(8.63)\mathbf{T}=\varepsilon^{\Phi}\mathbf{S}_{\Phi},\tag{8.63}
which is a fully tensorial expression. We hope that at this point in our narrative, it is no longer necessary to exalt the benefits of such expressions over those that require a special parameterization, such as
T=R(s)(8.64)\mathbf{T}=\mathbf{R}^{\prime}\left( s\right)\tag{8.64}
or those that feature an imbalance of indices.
Let us now apply the same strategy to the differentiation operator
Φ,(8.65)\nabla_{\Phi},\tag{8.65}
i.e. consider the combination
εΦΦ   .(8.66)\varepsilon^{\Phi}\nabla_{\Phi}\ \ \ .\tag{8.66}
It is an invariant operation that coincides with
dds(8.67)\frac{d}{ds}\tag{8.67}
under an arc length parameterization. Thus, as before, this is the case under all parameterizations. In other words, εΦΦ\varepsilon^{\Phi}\nabla_{\Phi} represents the derivative d/dsd/ds in any coordinate system SΦS^{\Phi}. This insight enables us to convert all of our earlier findings to a fully tensorial form.
In particular, the unit normal T\mathbf{T}, which has heretofore eluded our analysis, is given by
T=εΦΦR.(8.68)\mathbf{T}=\varepsilon^{\Phi}\nabla_{\Phi}\mathbf{R}.\tag{8.68}
The curvature normal B\mathbf{B} is given by
B=εΦΦT.(8.69)\mathbf{B}=\varepsilon^{\Phi}\nabla_{\Phi}\mathbf{T}.\tag{8.69}
The expression for the absolute curvature σ\sigma is unchanged, i.e.
σ=BB.(8.70)\sigma=\sqrt{\mathbf{B}\cdot\mathbf{B}}.\tag{8.70}
However, since
BB=εΦΦTεΨΨT(8.71)\mathbf{B}\cdot\mathbf{B}=\varepsilon^{\Phi}\nabla_{\Phi}\mathbf{T} \cdot\varepsilon^{\Psi}\nabla_{\Psi}\mathbf{T}\tag{8.71}
and
εΦεΨ=SΦΨ,(8.72)\varepsilon^{\Phi}\varepsilon^{\Psi}=S^{\Phi\Psi},\tag{8.72}
we have
BB=ΦTΦT(8.73)\mathbf{B}\cdot\mathbf{B}=\nabla_{\Phi}\mathbf{T}\cdot\nabla^{\Phi}\mathbf{T}\tag{8.73}
and, therefore,
σ=ΦTΦT.(8.74)\sigma=\sqrt{\nabla_{\Phi}\mathbf{T}\cdot\nabla^{\Phi}\mathbf{T}}.\tag{8.74}
The definitions of the principal normal P\mathbf{P} and the binormal Q\mathbf{Q} are unchanged, i.e.
σP=B,(8.75)\sigma\mathbf{P}=\mathbf{B},\tag{8.75}
and
Q=T×P.(8.76)\mathbf{Q}=\mathbf{T}\times\mathbf{P.}\tag{8.76}
Finally, the torsion τ\tau is defined by the equation
τQ=σT+εΦΦP(8.77)\tau\mathbf{Q}=\sigma\mathbf{T}+\varepsilon^{\Phi}\nabla_{\Phi}\mathbf{P}\tag{8.77}
and given by tensor expression
τ=εΦQΦP.(8.78)\tau=\varepsilon^{\Phi}\mathbf{Q}\cdot\nabla_{\Phi}\mathbf{P.}\tag{8.78}
We have thus succeeded in removing the first shortcoming of our earlier analysis which was the need for an arc length parameterization. However, we still do not have a practical analytical framework since we still lack a coordinate system in the ambient space. We will remove this second and final shortcoming in the next Section.
Let us now introduce a coordinate system ZiZ^{i} in the ambient space. This will enable us to consider equations of the curve, such as the parametric equations
x(γ)=rcosγ          (8.79)y(γ)=rsinγ          (8.80)z(γ)=cγ          (8.81)\begin{aligned}x\left( \gamma\right) & =r\cos\gamma\ \ \ \ \ \ \ \ \ \ \left(8.79\right)\\y\left( \gamma\right) & =r\sin\gamma\ \ \ \ \ \ \ \ \ \ \left(8.80\right)\\z\left( \gamma\right) & =c\gamma\ \ \ \ \ \ \ \ \ \ \left(8.81\right)\end{aligned}
for a right-handed helix in Cartesian coordinates. In general, the equations of the curve read
Zi=Zi(S),(8.82)Z^{i}=Z^{i}\left( S\right) ,\tag{8.82}
which, in this compact form, coincides with the equations of the surface for multi-dimensional surfaces. However, when we unpack these equations, we observe that there is only one independent variable S1S^{1}, i.e.
Z1=Z1(S1)          (8.83)Z2=Z2(S1)          (8.84)Z3=Z3(S1).          (8.85)\begin{aligned}Z^{1} & =Z^{1}\left( S^{1}\right)\ \ \ \ \ \ \ \ \ \ \left(8.83\right)\\Z^{2} & =Z^{2}\left( S^{1}\right)\ \ \ \ \ \ \ \ \ \ \left(8.84\right)\\Z^{3} & =Z^{3}\left( S^{1}\right) .\ \ \ \ \ \ \ \ \ \ \left(8.85\right)\end{aligned}
The shift tensor ZΦiZ_{\Phi}^{i} is the partial derivative of the equation of the curve with respect to the curve coordinate SΦS^{\Phi}, i.e.
ZΦi=Zi(S)SΦ.(8.86)Z_{\Phi}^{i}=\frac{\partial Z^{i}\left( S\right) }{\partial S^{\Phi}}.\tag{8.86}
The shift tensor relates the ambient and the surface basis vectors
SΦ=ZiZΦi(8.87)\mathbf{S}_{\Phi}=\mathbf{Z}_{i}Z_{\Phi}^{i}\tag{8.87}
and, therefore, the elements of the shift tensor represent the coefficients of the vector SΦ\mathbf{S}_{\Phi} with respect to the ambient covariant basis Zi\mathbf{Z}_{i}. The components TiT^{i} of the unit tangent T\mathbf{T} are given by
Ti=εΦZΦi.(8.88)T^{i}=\varepsilon^{\Phi}Z_{\Phi}^{i}.\tag{8.88}
The curve metric tensor SΦΨS_{\Phi\Psi} can be obtained from the ambient metric tensor according to the equation
SΦΨ=ZijZΦiZΨj.(8.89)S_{\Phi\Psi}=Z_{ij}Z_{\Phi}^{i}Z_{\Psi}^{j}.\tag{8.89}
According to this formula, the determinant SS of SΦΨS_{\Phi\Psi}, which also equals its sole element S11S_{11}, is given by
S=S11=ZijdZidS1dZjdS1.(8.90)S=S_{11}=Z_{ij}\frac{dZ^{i}}{dS^{1}}\frac{dZ^{j}}{dS^{1}}.\tag{8.90}
Thus, the arithmetic expression for the geometric integral
SFdS(8.91)\int_{S}FdS\tag{8.91}
reads
SFdS=abF(S)Zij(S)dZi(S)dS1dZj(S)dS1dS1,(8.92)\int_{S}FdS=\int_{a}^{b}F\left( S\right) \sqrt{Z_{ij}\left( S\right) \frac{dZ^{i}\left( S\right) }{dS^{1}}\frac{dZ^{j}\left( S\right) }{dS^{1} }}dS^{1},\tag{8.92}
which is of course equivalent to the familiar formula
SFdS=γ0γ1F(γ)Zij(γ)dZi(γ)dγdZj(γ)dγdγ(8.93)\int_{S}FdS=\int_{\gamma_{0}}^{\gamma_{1}}F\left( \gamma\right) \sqrt {Z_{ij}\left( \gamma\right) \frac{dZ^{i}\left( \gamma\right) }{d\gamma }\frac{dZ^{j}\left( \gamma\right) }{d\gamma}}d\gamma\tag{8.93}
that we have encountered numerous times.
The covariant derivative Θ\nabla_{\Theta} applied to a tensor UjΨiΦU_{j\Psi }^{i\Phi} with a representative collection of indices is entirely analogous to the surface covariant derivative,
ΘUjΨiΦ=UjΨiΦSΘ+ZΘkΓkmiUjΨmΦZΘkΓkjmUmΨiΦ+ΓΘΩΦUjΨiΩΓΘΨΩUjΩiΦ.(8.94)\nabla_{\Theta}U_{j\Psi}^{i\Phi}=\frac{\partial U_{j\Psi}^{i\Phi}}{\partial S^{\Theta}}+Z_{\Theta}^{k}\Gamma_{km}^{i}U_{j\Psi}^{m\Phi}-Z_{\Theta} ^{k}\Gamma_{kj}^{m}U_{m\Psi}^{i\Phi}+\Gamma_{\Theta\Omega}^{\Phi}U_{j\Psi }^{i\Omega}-\Gamma_{\Theta\Psi}^{\Omega}U_{j\Omega}^{i\Phi}.\tag{8.94}
This differential operator has all of the same properties as the covariant derivative on the surface described in Chapter 5. In particular, it produces tensor outputs for tensor inputs, satisfies the sum and the product rules, commutes with contraction and obeys the chain rule for curve restrictions UjiU_{j}^{i} of ambient variants
ΘUji=ZΘkkUji.(8.95)\nabla_{\Theta}U_{j}^{i}=Z_{\Theta}^{k}\nabla_{k}U_{j}^{i}.\tag{8.95}
Let BΦΨiB_{\Phi\Psi}^{i} be the ambient components of the vector curvature tensor BΦΨ\mathbf{B}_{\Phi\Psi}, i.e.
BΦΨ=BΦΨiZi.(8.96)\mathbf{B}_{\Phi\Psi}=B_{\Phi\Psi}^{i}\mathbf{Z}_{i}.\tag{8.96}
Then BΦΨiB_{\Phi\Psi}^{i} are given by
BΦΨi=ΦZΨi(8.97)B_{\Phi\Psi}^{i}=\nabla_{\Phi}Z_{\Psi}^{i}\tag{8.97}
and the components BΦiΦB_{\Phi}^{i\Phi} of the vector mean curvature BΦΦ\mathbf{B}_{\Phi}^{\Phi}, a.k.a. the curvature normal B\mathbf{B}, are given by
BΦiΦ=ΦZΦi.(8.98)B_{\Phi}^{i\Phi}=\nabla^{\Phi}Z_{\Phi}^{i}.\tag{8.98}
Since the vector mean curvature BΦΦ\mathbf{B}_{\Phi}^{\Phi} coincides with the curvature normal B\mathbf{B} we can denote the same components by BiB^{i}, i.e.
Bi=ΦZΦi.(8.99)B^{i}=\nabla^{\Phi}Z_{\Phi}^{i}.\tag{8.99}
Then the absolute curvature σ\sigma, being the magnitude of B\mathbf{B}, is given by
σ=BiBi.(8.100)\sigma=\sqrt{B_{i}B^{i}}.\tag{8.100}
Let Pi P^{i}\mathbf{\ }be the ambient components of the principal normal P\mathbf{P} and QiQ^{i} be the ambient components of the binormal Q\mathbf{Q}. Then PiP^{i} are given by
Pi=Biσ(8.101)P^{i}=\frac{B^{i}}{\sigma}\tag{8.101}
and QiQ^{i} are given by
Qi=1σεijkεΦZjΦBk(8.102)Q^{i}=\frac{1}{\sigma}\varepsilon^{ijk}\varepsilon^{\Phi}Z_{j\Phi}B_{k}\tag{8.102}
since Q\mathbf{Q} is the cross product of T\mathbf{T} and P\mathbf{P}. Since Bk=ΨZkΨB_{k}=\nabla_{\Psi}Z_{k}^{\Psi}, we find
Qi=1σεijkεΦZjΦΨZkΨ.(8.103)Q^{i}=\frac{1}{\sigma}\varepsilon^{ijk}\varepsilon^{\Phi}Z_{j\Phi}\nabla _{\Psi}Z_{k}^{\Psi}.\tag{8.103}
Finally, since
τ=εΦQΦP,(8.78)\tau=\varepsilon^{\Phi}\mathbf{Q}\cdot\nabla_{\Phi}\mathbf{P,} \tag{8.78}
we have
τ=εΦQiΦPi=εΦQiΦ(1σΨZΨi).(8.104)\tau=\varepsilon^{\Phi}Q_{i}\nabla_{\Phi}P^{i}=\varepsilon^{\Phi}Q_{i} \nabla_{\Phi}\left( \frac{1}{\sigma}\nabla^{\Psi}Z_{\Psi}^{i}\right) \mathbf{.}\tag{8.104}
Thus, we have succeeded in expressing all of the differential characteristics of curves in full tensor form. In the next Section we will apply these formulas to curves specified in Cartesian ambient coordinates.
Refer the ambient space to Cartesian coordinates x,y,zx,y,z and let the equations of the curve read
x=x(γ)          (8.105)y=y(γ)          (8.106)z=z(γ),          (8.107)\begin{aligned}x & =x\left( \gamma\right)\ \ \ \ \ \ \ \ \ \ \left(8.105\right)\\y & =y\left( \gamma\right)\ \ \ \ \ \ \ \ \ \ \left(8.106\right)\\z & =z\left( \gamma\right) ,\ \ \ \ \ \ \ \ \ \ \left(8.107\right)\end{aligned}
where γ\gamma represents the curve coordinate S1S^{1}.
The ambient metric tensors ZijZ_{ij} and ZijZ^{ij} correspond to the identity matrix, i.e.
[Zij]=[Zij]=[100010001].(8.108)\left[ Z_{ij}\right] =\left[ Z^{ij}\right] =\left[ \begin{array} {ccc} 1 & 0 & 0\\ 0 & 1 & 0\\ 0 & 0 & 1 \end{array} \right] .\tag{8.108}
Now, let the symbols xγ,yγ,zγx_{\gamma},y_{\gamma},z_{\gamma}, xγγ,yγγ,zγγx_{\gamma\gamma },y_{\gamma\gamma},z_{\gamma\gamma}, and xγγγ,yγγγ,zγγγx_{\gamma\gamma\gamma} ,y_{\gamma\gamma\gamma},z_{\gamma\gamma\gamma} denote the first, second, and third derivatives of the functions x(γ)x\left( \gamma\right) , y(γ)y\left( \gamma\right) , and z(γ)z\left( \gamma\right) . Below we give the summary of the coordinate representation of all of the differential characteristics of a curve. The derivations are left as an exercise.
The shift tensors are given by
[ZΦi]=[ZiΦ]=[xγyγzγ](8.109)\left[ Z_{\Phi}^{i}\right] =\left[ Z_{i\Phi}\right] =\left[ \begin{array} {c} x_{\gamma}\\ y_{\gamma}\\ z_{\gamma} \end{array} \right]\tag{8.109}
Thus, the curve metrics are given by
[SΦΨ]=[xγ2+yγ2+zγ2]          (8.110)[SΦΨ]=[1xγ2+yγ2+zγ2]          (8.111)S=xγ2+yγ2+zγ2.          (8.112)εΦ=xγ2+yγ2+zγ2          (8.113)εΦ=1xγ2+yγ2+zγ2.          (8.114)\begin{aligned}\left[ S_{\Phi\Psi}\right] & =\left[ x_{\gamma}^{2}+y_{\gamma} ^{2}+z_{\gamma}^{2}\right]\ \ \ \ \ \ \ \ \ \ \left(8.110\right)\\\left[ S^{\Phi\Psi}\right] & =\left[ \frac{1}{x_{\gamma}^{2}+y_{\gamma }^{2}+z_{\gamma}^{2}}\right]\ \ \ \ \ \ \ \ \ \ \left(8.111\right)\\\sqrt{S} & =\sqrt{x_{\gamma}^{2}+y_{\gamma}^{2}+z_{\gamma}^{2}}.\ \ \ \ \ \ \ \ \ \ \left(8.112\right)\\\varepsilon_{\Phi} & =\sqrt{x_{\gamma}^{2}+y_{\gamma}^{2}+z_{\gamma}^{2}}\ \ \ \ \ \ \ \ \ \ \left(8.113\right)\\\varepsilon^{\Phi} & =\frac{1}{\sqrt{x_{\gamma}^{2}+y_{\gamma}^{2} +z_{\gamma}^{2}}}.\ \ \ \ \ \ \ \ \ \ \left(8.114\right)\end{aligned}
Raising the curve index on the shift tensors, we get
[ZiΦ]=[ZiΦ]=1xγ2+yγ2+zγ2[xγyγzγ].(8.115)\left[ Z^{i\Phi}\right] =\left[ Z_{i}^{\Phi}\right] =\frac{1}{x_{\gamma }^{2}+y_{\gamma}^{2}+z_{\gamma}^{2}}\left[ \begin{array} {c} x_{\gamma}\\ y_{\gamma}\\ z_{\gamma} \end{array} \right] .\tag{8.115}
The components TiT^{i} of the unit tangent are
[Ti]=1xγ2+yγ2+zγ2[xγyγzγ].(8.116)\left[ T^{i}\right] =\frac{1}{\sqrt{x_{\gamma}^{2}+y_{\gamma}^{2}+z_{\gamma }^{2}}}\left[ \begin{array} {c} x_{\gamma}\\ y_{\gamma}\\ z_{\gamma} \end{array} \right] .\tag{8.116}
The Christoffel symbol ΓΦΨΘ\Gamma_{\Phi\Psi}^{\Theta} is given by
[ΓΦΨΔ]=xγxγγ+yγyγγ+zγzγγxγ2+yγ2+zγ2(8.117)\left[ \Gamma_{\Phi\Psi}^{\Delta}\right] =\frac{x_{\gamma}x_{\gamma\gamma }+y_{\gamma}y_{\gamma\gamma}+z_{\gamma}z_{\gamma\gamma}}{x_{\gamma} ^{2}+y_{\gamma}^{2}+z_{\gamma}^{2}}\tag{8.117}
For the sake of brevity in the upcoming equations, introduce the symbols
X=yγγzγyγzγγ          (8.118)Y=zγγxγzγxγγ          (8.119)Z=xγγyγxγyγγ          (8.120)\begin{aligned}X & =y_{\gamma\gamma}z_{\gamma}-y_{\gamma}z_{\gamma\gamma}\ \ \ \ \ \ \ \ \ \ \left(8.118\right)\\Y & =z_{\gamma\gamma}x_{\gamma}-z_{\gamma}x_{\gamma\gamma}\ \ \ \ \ \ \ \ \ \ \left(8.119\right)\\Z & =x_{\gamma\gamma}y_{\gamma}-x_{\gamma}y_{\gamma\gamma}\ \ \ \ \ \ \ \ \ \ \left(8.120\right)\end{aligned}
For the components Bi B^{i}\ of the curvature normal, we have
[Bi]=1(xγ2+yγ2+zγ2)2[yγZzγYzγXxγZxγYyγX].(8.121)\left[ B^{i}\right] =\frac{1}{\left( x_{\gamma}^{2}+y_{\gamma} ^{2}+z_{\gamma}^{2}\right) ^{2}}\left[ \begin{array} {c} y_{\gamma}Z-z_{\gamma}Y\\ z_{\gamma}X-x_{\gamma}Z\\ x_{\gamma}Y-y_{\gamma}X \end{array} \right] .\tag{8.121}
As a result, the absolute curvature σ\sigma is given by
σ=X2+Y2+Z2(xγ2+yγ2+zγ2)3/2.(8.122)\sigma=\frac{\sqrt{X^{2}+Y^{2}+Z^{2}}}{\left( x_{\gamma}^{2}+y_{\gamma} ^{2}+z_{\gamma}^{2}\right) ^{3/2}}.\tag{8.122}
Finally the components QiQ^{i} of the binormal Q\mathbf{Q} are given by
[Qi]=1X2+Y2+Z2[XYZ](8.123)\left[ Q^{i}\right] =-\frac{1}{\sqrt{X^{2}+Y^{2}+Z^{2}}}\left[ \begin{array} {c} X\\ Y\\ Z \end{array} \right]\tag{8.123}
and torsion τ\tau is given by
τ=xγγγZ+yγγγY+zγγγZX2+Y2+Z2.(8.124)\tau=-\frac{x_{\gamma\gamma\gamma}Z+y_{\gamma\gamma\gamma}Y+z_{\gamma \gamma\gamma}Z}{X^{2}+Y^{2}+Z^{2}}.\tag{8.124}
The analysis of a helix can be accomplished simply by substituting its equations of the curve into the identities obtained in the previous section. However, it would be more insightful, and perhaps even simpler, to retrace the steps of the previous Section for the specific equations of the curve.
Let us simultaneously analyze a right-handed and a left-handed helix by writing the equations of the curve in the form
x(γ)=Σrcosγ          (8.125)y(γ)=rsinγ          (8.126)z(γ)=cγ,          (8.127)\begin{aligned}x\left( \gamma\right) & =\Sigma r\cos\gamma\ \ \ \ \ \ \ \ \ \ \left(8.125\right)\\y\left( \gamma\right) & =r\sin\gamma\ \ \ \ \ \ \ \ \ \ \left(8.126\right)\\z\left( \gamma\right) & =c\gamma,\ \ \ \ \ \ \ \ \ \ \left(8.127\right)\end{aligned}
where Σ=1\Sigma=1 for a right-handed helix and 1-1 for a left-handed helix.
The shift tensor ZΦiZ_{\Phi}^{i}, which can also be interpreted as the components of the curve basis textbf{SS}Φ_{\Phi} with respect to the ambient basis Zi\mathbf{Z}_{i}, is given by
[ZΦi]=[ZiΦ]=[Σrsinγrcosγc].(8.128)\left[ Z_{\Phi}^{i}\right] =\left[ Z_{i\Phi}\right] =\left[ \begin{array} {r} -\Sigma r\sin\gamma\\ r\cos\gamma\\ c \end{array} \right] .\tag{8.128}
The covariant metric tensor SΦΨ=ZijZΦiZΨjS_{\Phi\Psi}=Z_{ij}Z_{\Phi}^{i}Z_{\Psi}^{j} corresponds to the 1×11\times1 matrix [ZΦi]T[ZΦi]\left[ Z_{\Phi}^{i}\right] ^{T}\left[ Z_{\Phi}^{i}\right] , which yields
[SΦΨ]=[r2+c2].(8.129)\left[ S_{\Phi\Psi}\right] =\left[ r^{2}+c^{2}\right] .\tag{8.129}
Note that it is independent of the coordinate γ\gamma which will therefore be the case for the rest of the metrics. Furthermore, the Christoffel symbol ΓΦΨΘ\Gamma_{\Phi\Psi}^{\Theta} will vanish, i.e.
ΓΦΨΘ=0.(8.130)\Gamma_{\Phi\Psi}^{\Theta}=0.\tag{8.130}
The contravariant metric tensor SΦΨS^{\Phi\Psi} is given by
[SΦΨ]=[1r2+c2].(8.131)\left[ S^{\Phi\Psi}\right] =\left[ \frac{1}{r^{2}+c^{2}}\right] .\tag{8.131}
The length element S\sqrt{S} is the square root of the sole element on the covariant metric tensor, i.e.
S=r2+c2.(8.132)\sqrt{S}=\sqrt{r^{2}+c^{2}}.\tag{8.132}
The Levi-Civita symbol εΦ\varepsilon_{\Phi} and εΦ\varepsilon^{\Phi} are given by
[εΦ]=[r2+c2]          (8.133)[εΦ]=[1r2+c2].          (8.134)\begin{aligned}\left[ \varepsilon_{\Phi}\right] & =\left[ \sqrt{r^{2}+c^{2}}\right]\ \ \ \ \ \ \ \ \ \ \left(8.133\right)\\\left[ \varepsilon^{\Phi}\right] & =\left[ \frac{1}{\sqrt{r^{2}+c^{2}} }\right] .\ \ \ \ \ \ \ \ \ \ \left(8.134\right)\end{aligned}
The components TiT^{i} of the unit tangent are given by εΦZΦi\varepsilon^{\Phi }Z_{\Phi}^{i}, i.e.
[Ti]=1r2+c2[Σrsinγrcosγc].(8.135)\left[ T^{i}\right] =\frac{1}{\sqrt{r^{2}+c^{2}}}\left[ \begin{array} {r} -\Sigma r\sin\gamma\\ r\cos\gamma\\ c \end{array} \right] .\tag{8.135}
Since the Christoffel symbol vanishes, the components of the vector curvature tensor BΦΨiB_{\Phi\Psi}^{i} are given by
BΦΨi=εΦTiγ(8.136)B_{\Phi\Psi}^{i}=\varepsilon^{\Phi}\frac{\partial T^{i}}{\partial\gamma}\tag{8.136}
which yields
[BΦΨi]=r[Σcosγsinγ0].(8.137)\left[ B_{\Phi\Psi}^{i}\right] =-r\left[ \begin{array} {r} \Sigma\cos\gamma\\ \sin\gamma\\ 0 \end{array} \right] .\tag{8.137}
The components BiB^{i} of the curvature normal B\mathbf{B}, where Bi=B^{i}= BΦiΦB_{\Phi}^{i\Phi} is given by
[Bi]=[BΦiΦ]=rr2+c2[Σcosγsinγ0].(8.138)\left[ B^{i}\right] =\left[ B_{\Phi}^{i\Phi}\right] =-\frac{r}{r^{2} +c^{2}}\left[ \begin{array} {r} \Sigma\cos\gamma\\ \sin\gamma\\ 0 \end{array} \right] .\tag{8.138}
Note that the zero third entry confirms the fact that we intuited in Chapter TBD of Introduction to Tensor Calculus -- namely, that the curvature normal points in the direction orthogonal to the axis of rotation.
The absolute curvature σ\sigma is the magnitude of BiB^{i}, i.e.
σ=rr2+c2.(8.139)\sigma=\frac{r}{r^{2}+c^{2}}.\tag{8.139}
The components PiP^{i} of the principal normal P\mathbf{P}, are given by Pi=Bi/σP^{i}=B^{i}/\sigma, i.e.
[Pi]=[Σcosγsinγ0].(8.140)\left[ P^{i}\right] =-\left[ \begin{array} {r} \Sigma\cos\gamma\\ \sin\gamma\\ 0 \end{array} \right] .\tag{8.140}
The binormal Q\mathbf{Q} is the cross product of T\mathbf{T} and P\mathbf{P}. It is left as an exercise to show that the components QiQ^{i} are given by
[Qi]=Σc2+r2[Σcsinγccosγr].(8.141)\left[ Q^{i}\right] =\frac{\Sigma}{\sqrt{c^{2}+r^{2}}}\left[ \begin{array} {r} \Sigma c\sin\gamma\\ -c\cos\gamma\\ r \end{array} \right] .\tag{8.141}
The remaining task is to calculate the torsion τ\tau. We have
τ=QidPidγ(8.142)\tau=Q_{i}\frac{dP^{i}}{d\gamma}\tag{8.142}
and therefore
τ=Σcc2+r2.(8.143)\tau=\frac{\Sigma c}{c^{2}+r^{2}}.\tag{8.143}
Thus, we now have a precise confirmation of the fact that we established qualitatively in Chapter TBD of Introduction to Tensor Calculus. Namely, that a right-handed helix is characterized by positive torsion while a left-handed helix is characterized by negative torsion.
Exercise 8.1Show that ΦSΨ\nabla_{\Phi}\mathbf{S}_{\Psi} is orthogonal to the tangent space, i.e.
SΘΦSΨ=0.(8.144)\mathbf{S}^{\Theta}\cdot\nabla_{\Phi}\mathbf{S}_{\Psi}=0.\tag{8.144}
Exercise 8.2Show that under an arc length parameterization, we have
SΦΨ=1          (8.145)SΦΨ=1          (8.146)S=1          (8.147)εΦ=1          (8.148)εΦ=1          (8.149)ΓΦΨΔ=0          (8.150)\begin{aligned}S_{\Phi\Psi} & =1\ \ \ \ \ \ \ \ \ \ \left(8.145\right)\\S^{\Phi\Psi} & =1\ \ \ \ \ \ \ \ \ \ \left(8.146\right)\\\sqrt{S} & =1\ \ \ \ \ \ \ \ \ \ \left(8.147\right)\\\varepsilon_{\Phi} & =1\ \ \ \ \ \ \ \ \ \ \left(8.148\right)\\\varepsilon^{\Phi} & =1\ \ \ \ \ \ \ \ \ \ \left(8.149\right)\\\Gamma_{\Phi\Psi}^{\Delta} & =0\ \ \ \ \ \ \ \ \ \ \left(8.150\right)\end{aligned}
Exercise 8.3Confirm by a direct calculation under a general coordinate system SΦS^{\Phi} that
ΦεΨ=0.(8.151)\nabla_{\Phi}\varepsilon_{\Psi}=0.\tag{8.151}
Meanwhile
ΦS=ΓΨΦΨS,(8.152)\nabla_{\Phi}\sqrt{S}=\Gamma_{\Psi\Phi}^{\Psi}\sqrt{S},\tag{8.152}
despite the fact that εΦ\varepsilon_{\Phi} and S\sqrt{S} have equal values.
Exercise 8.4Derive all of the results presented in Section 8.6. Using a computer algebra system is recommended!
Exercise 8.5Repeat the analysis of a helix in cylindrical ambient coordinates.
Send feedback to Pavel