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 \(\mathbf{T}\), the curvature normal \(\mathbf{B}\) along with the associated concepts of the absolute curvature \(\sigma\) and the principal normal \(\mathbf{P}\), and the binormal \(\mathbf{Q}\) along with the associated concept of the torsion \(\tau\). The vectors \(\mathbf{T}\), \(\mathbf{P}\), and \(\mathbf{Q}\) are related by the Frenet equations given below.
The unit tangent \(\mathbf{T}\) is the derivative of the position vector function \(\mathbf{R}\left( s\right) \), i.e.
\[ \mathbf{T}=\mathbf{R}_{s}.\tag{8.1} \]
The derivative \(\mathbf{T}_{s}\) is the curvature normal \(\mathbf{B}\), i.e.
\[ \mathbf{B}=\mathbf{T}_{s}.\tag{8.2} \]
Equivalently, \(\mathbf{B}\) is the second derivative of \(\mathbf{R}\), i.e.
\[ \mathbf{B}=\mathbf{R}_{ss}.\tag{8.3} \]
The curvature normal \(\mathbf{B}\) is orthogonal to \(\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 \(\mathbf{P}\) that points in the direction of \(\mathbf{B}\), i.e.
\[ \mathbf{B}=\sigma\mathbf{P,}\tag{8.4} \]
is known as the principal normal. The plane spanned by \(\mathbf{T}\) and \(\mathbf{P}\) is known as the osculating plane. Note that when \(\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 \(\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 \(\mathbf{Q}\) is to supplement the principal normal \(\mathbf{P}\) in providing a basis for the normal space. The binormal \(\mathbf{Q}\) is chosen to be the unit vector that is orthogonal to \(\mathbf{P}\) in the normal plane. Between the two possible vectors that satisfy this condition, \(\mathbf{Q}\) is chosen so that the set \(\mathbf{T,P,Q} \) is positively oriented. With the help of the cross product, the binormal \(\mathbf{Q}\) is given by
\[ \mathbf{Q}=\mathbf{T\times P.}\tag{8.5} \]
The set \(\mathbf{T,P,Q}\) is referred to as a local frame or the Frenet--Serret frame.
(8.6)
The derivative \(\mathbf{P}^{\prime}\left( s\right) \) is orthogonal to \(\mathbf{P}\) and is therefore found in the plane spanned by \(\mathbf{T}\) and \(\mathbf{Q}\). Its projection onto \(\mathbf{T}\) turns out to be \(-\sigma\). Meanwhile its projection onto \(\mathbf{Q}\) is an additional curve characteristic known as the torsion \(\tau\), i.e.
\[ \mathbf{P}_{s}=-\sigma\mathbf{T}+\tau\mathbf{Q.}\tag{8.7} \]
Thus, we have
\[ \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 \(\mathbf{Q} ^{\prime}\left( s\right) \) of the binormal is given by
\[ \mathbf{Q}_{s}=-\tau\mathbf{P}.\tag{8.9} \]
As a collection, the equations
\[ \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
\[ \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 \(\sigma\left( s\right) \) and \(\tau\left( s\right) \) of the arc length.
In the present narrative, where the unit tangent \(\mathbf{T}\) is defined as the derivative of \(\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 \(1\). Unlike \(\mathbf{T}\), the curvature normal \(\mathbf{B}\) is a full-fledged invariant. However, the binormal \(\mathbf{Q}\), being the cross-product of \(\mathbf{T}\) and \(\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 \(1\). 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 \(\mathbf{T}\), similar to the view that we took on the choice of the unit normal \(\mathbf{N}\) of a hypersurface. Namely, we could have a priori chosen a particular direction for the unit tangent \(\mathbf{T}\), which would have made it an invariant along with \(\mathbf{P}\), \(\mathbf{Q}\), \(\sigma\), and \(\tau\). In that case, however, the derivative \(\mathbf{R}^{\prime}\left( s\right) \) would actually produce \(\mathbf{T}\) only if the orientation of the parameterization is consistent with the a priori choice of \(\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 \(z\)-axis is aligned with the axis of the helix, then a right-handed helix may be given by the equations
\[ \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
\[ \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 \(r\) corresponds to the radius of the helix while the value \(2\pi c\), known as the pitch, is the vertical distance between consecutive loops.
  (8.20)
Recall that the arc length \(s\) as a function of the parameter \(\gamma\) is given by the integral
\[ 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\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 \(s\), the resulting equations would read
\[ \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
\[ \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 \(s\) is given by the integral
\[ 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 \(L\) 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 \(1\)) 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 \(L^{1}\) will be denoted by \(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{\(S\)}\(_{\Phi}\) -- consisting of a single vector textbf{\(S_{1}\)}, of course -- is defined by
\[ \mathbf{S}_{\Phi}=\frac{\partial\mathbf{R}\left( S\right) }{\partial S^{\Phi}}.\tag{8.30} \]
Since \(\mathbf{R}\left( S\right) \) is an ordinary function of one variable \(S^{1}\), it would have perhaps been more appropriate to use the ordinary derivative sign, as in
\[ \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, \(\mathbf{S}_{\Phi}\) is tangential to the curve as illustrated in the following figure.
(8.32)
The covariant basis textbf{\(S\)}\(_{\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
\[ \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^{\Phi^{\prime}}\) is an alternative parameterization of the curve and that textbf{\(S\)}\(_{\Phi^{\prime}}\) is the corresponding covariant basis. The indicial signature of (the tensor) textbf{\(S\)}\(_{\Phi}\) tells us that its relation to textbf{\(S\)}\(_{\Phi^{\prime}}\) reads
\[ \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_{\Phi\Psi}\) is defined by the familiar formula
\[ S_{\Phi\Psi}=\mathbf{S}_{\Phi}\cdot\mathbf{S}_{\Psi}.\tag{8.36} \]
It has a single element \(S_{11}\) and thus corresponds to a \(1\times1\) matrix. Once again, the benefit of enumerating the single element by a pair of indices is knowing the exact manner in which \(S_{\Phi\Psi}\) transforms under a change of variables. Namely, we know that \(S_{\Phi\Psi}\) is related to \(S_{\Phi ^{\prime}\Psi^{\prime}}\) by the identity
\[ 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^{\Phi\Psi}\) is the "matrix" inverse of \(S_{\Phi\Psi}\), i.e.
\[ S_{\Phi\Psi}S^{\Psi\Theta}=\delta_{\Phi}^{\Theta}.\tag{8.38} \]
In actuality, this identity tells us that the only element of \(S^{\Psi\Theta}\) is the reciprocal of the only element of \(S_{\Phi\Psi}\). That element transforms according to the rule
\[ S^{\Phi^{\prime}\Psi^{\prime}}=S^{\Phi\Psi}J_{\Phi}^{\Phi^{\prime}}J_{\Psi }^{\Psi^{\prime}}.\tag{8.39} \]
The contravariant basis textbf{\(S\)}\(^{\Phi}\) is given by
\[ \mathbf{S}^{\Phi}=S^{\Phi\Psi}\mathbf{S}_{\Psi}.\tag{8.40} \]
The vector textbf{\(S\)}\(^{\Phi}\) points in the same direction as textbf{\(S\) }\(_{\Phi}\) and its length is the reciprocal of the length of textbf{\(S\) }\(_{\Phi}\).
Finally, the length element \(\sqrt{S}\) is the square root of the determinant of the covariant metric tensor \(S_{\Phi\Psi}\). Since \(S_{\Phi\Psi}\) corresponds to a \(1\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
\[ \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 \(1\) when \(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 \(1\) for an arc length parameterization and \(\varepsilon^{\Phi}\) is the (almost) unique contravariant tensor that equals \(1\) 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 \(\sqrt{S}\) and the covariant Levi-Civita symbol \(\varepsilon_{\Phi}\) have the exact same value. Nevertheless, while \(\sqrt{S}\) (treated as a variant of order \(0\)), is a relative tensor of weight \(1\), \(\varepsilon_{\Phi}\) (treated as a variant of order \(1\)) is an absolute tensor.
The Christoffel symbol \(\Gamma_{\Phi\Psi}^{\Theta}\) is given by the dot product
\[ \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
\[ \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
\[ \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 \(\Gamma_{ij}^{k}\) by the identity
\[ \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 \(\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,
\[ \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_{\Phi^{\prime}\Psi^{\prime}}^{\Theta}\) is given by
\[ 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_{\cdot\Theta\Phi\Psi}^{\Omega}\) is given by
\[ 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_{\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_{\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_{\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_{\Psi }^{\Phi}\) with a representative collection of curve indices is entirely analogous to the surface covariant derivative, i.e.
\[ \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 \(\mathbf{S}_{\Phi}\). This is, of course, due to curvature. Nevertheless, as before, \(\nabla_{\Phi}\mathbf{S}_{\Psi}\) is orthogonal to the tangent space. However, the familiar jump to the would-be curvature tensor \(B_{\Phi\Psi}\) via the identity
\[ \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 \(\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 \(2\). 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_{\Phi\Psi}\), we will introduce the vector curvature tensor \(\mathbf{B}_{\Phi\Psi}\) given by
\[ \mathbf{B}_{\Phi\Psi}=\nabla_{\Phi}\mathbf{S}_{\Psi}\tag{8.52} \]
or, equivalently,
\[ \mathbf{B}_{\Phi\Psi}=\nabla_{\Phi}\nabla_{\Psi}\mathbf{R.}\tag{8.53} \]
As we mentioned above, \(\mathbf{B}_{\Phi\Psi}\) is orthogonal to the one-dimensional tangent space, i.e.
\[ \mathbf{S}^{\Theta}\cdot\mathbf{B}_{\Phi\Psi}=0.\tag{8.54} \]
Of particular interest is the invariant \(\mathbf{B}_{\Phi}^{\Phi}\), analogous to mean curvature, given by
\[ \mathbf{B}_{\Phi}^{\Phi}=\nabla_{\Phi}\mathbf{S}^{\Phi}\tag{8.55} \]
or, equivalently
\[ \mathbf{B}_{\Phi}^{\Phi}=\nabla_{\Phi}\nabla^{\Phi}\mathbf{R.}\tag{8.56} \]
The invariant \(\mathbf{B}_{\Phi}^{\Phi}\) may be referred to as the vector mean curvature.
Not surprisingly, the invariant \(\mathbf{B}_{\Phi}^{\Phi}\) coincides with the curvature normal \(\mathbf{B}\). This can be easily seen by observing that, on the one hand, \(\mathbf{B}_{\Phi}^{\Phi}\) is an invariant and, on the other, it coincides with \(\mathbf{B}\) under an arc length parameterization. Indeed, since the Christoffel symbols vanish and \(S^{\Phi\Psi}=1\), we have
\[ \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 \(\mathbf{B}\) is an unqualified invariant, unlike the unit normal \(\mathbf{T}\) which is an invariant only with respect to orientation-preserving coordinate changes.
Since \(\mathbf{B}_{\Phi}^{\Phi}\) coincides with the curvature normal \(\mathbf{B}\), the absolute curvature \(\sigma\), being the magnitude of \(\mathbf{B}\), is given by the equation
\[ \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 \(\mathbf{B}\) emerged naturally (as the vector mean curvature \(\mathbf{B}_{\Phi}^{\Phi}\)) as we followed the blueprint established for surfaces, the unit tangent \(\mathbf{T}\) has eluded our analysis. Of course, since \(\mathbf{T}\) is the unit vector that points in the same direction as \(\mathbf{S}_{\Phi}\), we could have expressed it as
\[ \mathbf{T}=\frac{\mathbf{S}_{\Phi}}{\left\vert \mathbf{S}_{\Phi}\right\vert }\tag{8.59} \]
or as
\[ \mathbf{T}=\frac{\mathbf{S}_{1}}{\left\vert \mathbf{S}_{1}\right\vert }\tag{8.60} \]
or, since \(\left\vert \mathbf{S}_{1}\right\vert =\sqrt{S}\), as
\[ \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 \(\mathbf{S}_{\Phi}\) into the unit tangent \(\mathbf{T}\)? The covariant basis \(\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 \(\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 \(\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
\[ \varepsilon^{\Phi}\mathbf{S}_{\Phi}.\tag{8.62} \]
On the one hand, it is an invariant since \(\varepsilon^{\Phi}\) is a contravariant tensor while \(\mathbf{S}_{\Phi}\) is a covariant tensor. At the same time, under an arc-length parameterization, \(\varepsilon^{\Phi}=1\) and \(\mathbf{S}_{\Phi}\) is unit length. Consequently, being a vector of length \(1\) that points in the direction of \(\mathbf{S}_{\Phi}\), the vector \(\varepsilon ^{\Phi}\mathbf{S}_{\Phi}\) coincides with the unit tangent \(\mathbf{T}\). Therefore, it coincides with \(\mathbf{T}\) under all parameterizations. In other words, we always have
\[ \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
\[ \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
\[ \nabla_{\Phi},\tag{8.65} \]
i.e. consider the combination
\[ \varepsilon^{\Phi}\nabla_{\Phi}\ \ \ .\tag{8.66} \]
It is an invariant operation that coincides with
\[ \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/ds\) in any coordinate system \(S^{\Phi}\). This insight enables us to convert all of our earlier findings to a fully tensorial form.
In particular, the unit normal \(\mathbf{T}\), which has heretofore eluded our analysis, is given by
\[ \mathbf{T}=\varepsilon^{\Phi}\nabla_{\Phi}\mathbf{R}.\tag{8.68} \]
The curvature normal \(\mathbf{B}\) is given by
\[ \mathbf{B}=\varepsilon^{\Phi}\nabla_{\Phi}\mathbf{T}.\tag{8.69} \]
The expression for the absolute curvature \(\sigma\) is unchanged, i.e.
\[ \sigma=\sqrt{\mathbf{B}\cdot\mathbf{B}}.\tag{8.70} \]
However, since
\[ \mathbf{B}\cdot\mathbf{B}=\varepsilon^{\Phi}\nabla_{\Phi}\mathbf{T} \cdot\varepsilon^{\Psi}\nabla_{\Psi}\mathbf{T}\tag{8.71} \]
and
\[ \varepsilon^{\Phi}\varepsilon^{\Psi}=S^{\Phi\Psi},\tag{8.72} \]
we have
\[ \mathbf{B}\cdot\mathbf{B}=\nabla_{\Phi}\mathbf{T}\cdot\nabla^{\Phi}\mathbf{T}\tag{8.73} \]
and, therefore,
\[ \sigma=\sqrt{\nabla_{\Phi}\mathbf{T}\cdot\nabla^{\Phi}\mathbf{T}}.\tag{8.74} \]
The definitions of the principal normal \(\mathbf{P}\) and the binormal \(\mathbf{Q}\) are unchanged, i.e.
\[ \sigma\mathbf{P}=\mathbf{B},\tag{8.75} \]
and
\[ \mathbf{Q}=\mathbf{T}\times\mathbf{P.}\tag{8.76} \]
Finally, the torsion \(\tau\) is defined by the equation
\[ \tau\mathbf{Q}=\sigma\mathbf{T}+\varepsilon^{\Phi}\nabla_{\Phi}\mathbf{P}\tag{8.77} \]
and given by tensor expression
\[ \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 \(Z^{i}\) in the ambient space. This will enable us to consider equations of the curve, such as the parametric equations
\[ \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
\[ 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 \(S^{1}\), i.e.
\[ \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_{\Phi}^{i}\) is the partial derivative of the equation of the curve with respect to the curve coordinate \(S^{\Phi}\), i.e.
\[ 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
\[ \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 \(\mathbf{S}_{\Phi}\) with respect to the ambient covariant basis \(\mathbf{Z}_{i}\). The components \(T^{i}\) of the unit tangent \(\mathbf{T}\) are given by
\[ T^{i}=\varepsilon^{\Phi}Z_{\Phi}^{i}.\tag{8.88} \]
The curve metric tensor \(S_{\Phi\Psi}\) can be obtained from the ambient metric tensor according to the equation
\[ S_{\Phi\Psi}=Z_{ij}Z_{\Phi}^{i}Z_{\Psi}^{j}.\tag{8.89} \]
According to this formula, the determinant \(S\) of \(S_{\Phi\Psi}\), which also equals its sole element \(S_{11}\), is given by
\[ 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
\[ \int_{S}FdS\tag{8.91} \]
reads
\[ \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
\[ \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 \(U_{j\Psi }^{i\Phi}\) with a representative collection of indices is entirely analogous to the surface covariant derivative,
\[ \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 \(U_{j}^{i}\) of ambient variants
\[ \nabla_{\Theta}U_{j}^{i}=Z_{\Theta}^{k}\nabla_{k}U_{j}^{i}.\tag{8.95} \]
Let \(B_{\Phi\Psi}^{i}\) be the ambient components of the vector curvature tensor \(\mathbf{B}_{\Phi\Psi}\), i.e.
\[ \mathbf{B}_{\Phi\Psi}=B_{\Phi\Psi}^{i}\mathbf{Z}_{i}.\tag{8.96} \]
Then \(B_{\Phi\Psi}^{i}\) are given by
\[ B_{\Phi\Psi}^{i}=\nabla_{\Phi}Z_{\Psi}^{i}\tag{8.97} \]
and the components \(B_{\Phi}^{i\Phi}\) of the vector mean curvature \(\mathbf{B}_{\Phi}^{\Phi}\), a.k.a. the curvature normal \(\mathbf{B}\), are given by
\[ B_{\Phi}^{i\Phi}=\nabla^{\Phi}Z_{\Phi}^{i}.\tag{8.98} \]
Since the vector mean curvature \(\mathbf{B}_{\Phi}^{\Phi}\) coincides with the curvature normal \(\mathbf{B}\) we can denote the same components by \(B^{i}\), i.e.
\[ B^{i}=\nabla^{\Phi}Z_{\Phi}^{i}.\tag{8.99} \]
Then the absolute curvature \(\sigma\), being the magnitude of \(\mathbf{B}\), is given by
\[ \sigma=\sqrt{B_{i}B^{i}}.\tag{8.100} \]
Let \(P^{i}\mathbf{\ }\)be the ambient components of the principal normal \(\mathbf{P}\) and \(Q^{i}\) be the ambient components of the binormal \(\mathbf{Q}\). Then \(P^{i}\) are given by
\[ P^{i}=\frac{B^{i}}{\sigma}\tag{8.101} \]
and \(Q^{i}\) are given by
\[ Q^{i}=\frac{1}{\sigma}\varepsilon^{ijk}\varepsilon^{\Phi}Z_{j\Phi}B_{k}\tag{8.102} \]
since \(\mathbf{Q}\) is the cross product of \(\mathbf{T}\) and \(\mathbf{P}\). Since \(B_{k}=\nabla_{\Psi}Z_{k}^{\Psi}\), we find
\[ Q^{i}=\frac{1}{\sigma}\varepsilon^{ijk}\varepsilon^{\Phi}Z_{j\Phi}\nabla _{\Psi}Z_{k}^{\Psi}.\tag{8.103} \]
Finally, since
\[ \tau=\varepsilon^{\Phi}\mathbf{Q}\cdot\nabla_{\Phi}\mathbf{P,} \tag{8.78} \]
we have
\[ \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,z\) and let the equations of the curve read
\[ \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 \(S^{1}\).
The ambient metric tensors \(Z_{ij}\) and \(Z^{ij}\) correspond to the identity matrix, i.e.
\[ \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_{\gamma},y_{\gamma},z_{\gamma}\), \(x_{\gamma\gamma },y_{\gamma\gamma},z_{\gamma\gamma}\), and \(x_{\gamma\gamma\gamma} ,y_{\gamma\gamma\gamma},z_{\gamma\gamma\gamma}\) denote the first, second, and third derivatives of the functions \(x\left( \gamma\right) \), \(y\left( \gamma\right) \), and \(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
\[ \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
\[ \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
\[ \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 \(T^{i}\) of the unit tangent are
\[ \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
\[ \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
\[ \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 \(B^{i}\ \)of the curvature normal, we have
\[ \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
\[ \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 \(Q^{i}\) of the binormal \(\mathbf{Q}\) are given by
\[ \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
\[ \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
\[ \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 \(\Sigma=1\) for a right-handed helix and \(-1\) for a left-handed helix.
The shift tensor \(Z_{\Phi}^{i}\), which can also be interpreted as the components of the curve basis textbf{\(S\)}\(_{\Phi}\) with respect to the ambient basis \(\mathbf{Z}_{i}\), is given by
\[ \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_{\Phi\Psi}=Z_{ij}Z_{\Phi}^{i}Z_{\Psi}^{j}\) corresponds to the \(1\times1\) matrix \(\left[ Z_{\Phi}^{i}\right] ^{T}\left[ Z_{\Phi}^{i}\right] \), which yields
\[ \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.
\[ \Gamma_{\Phi\Psi}^{\Theta}=0.\tag{8.130} \]
The contravariant metric tensor \(S^{\Phi\Psi}\) is given by
\[ \left[ S^{\Phi\Psi}\right] =\left[ \frac{1}{r^{2}+c^{2}}\right] .\tag{8.131} \]
The length element \(\sqrt{S}\) is the square root of the sole element on the covariant metric tensor, i.e.
\[ \sqrt{S}=\sqrt{r^{2}+c^{2}}.\tag{8.132} \]
The Levi-Civita symbol \(\varepsilon_{\Phi}\) and \(\varepsilon^{\Phi}\) are given by
\[ \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 \(T^{i}\) of the unit tangent are given by \(\varepsilon^{\Phi }Z_{\Phi}^{i}\), i.e.
\[ \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_{\Phi\Psi}^{i}\) are given by
\[ B_{\Phi\Psi}^{i}=\varepsilon^{\Phi}\frac{\partial T^{i}}{\partial\gamma}\tag{8.136} \]
which yields
\[ \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 \(B^{i}\) of the curvature normal \(\mathbf{B}\), where \(B^{i}=\) \(B_{\Phi}^{i\Phi}\) is given by
\[ \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 \(B^{i}\), i.e.
\[ \sigma=\frac{r}{r^{2}+c^{2}}.\tag{8.139} \]
The components \(P^{i}\) of the principal normal \(\mathbf{P}\), are given by \(P^{i}=B^{i}/\sigma\), i.e.
\[ \left[ P^{i}\right] =-\left[ \begin{array} {r} \Sigma\cos\gamma\\ \sin\gamma\\ 0 \end{array} \right] .\tag{8.140} \]
The binormal \(\mathbf{Q}\) is the cross product of \(\mathbf{T}\) and \(\mathbf{P}\). It is left as an exercise to show that the components \(Q^{i}\) are given by
\[ \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
\[ \tau=Q_{i}\frac{dP^{i}}{d\gamma}\tag{8.142} \]
and therefore
\[ \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 \(\nabla_{\Phi}\mathbf{S}_{\Psi}\) is orthogonal to the tangent space, i.e.
\[ \mathbf{S}^{\Theta}\cdot\nabla_{\Phi}\mathbf{S}_{\Psi}=0.\tag{8.144} \]
Exercise 8.2Show that under an arc length parameterization, we have
\[ \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^{\Phi}\) that
\[ \nabla_{\Phi}\varepsilon_{\Psi}=0.\tag{8.151} \]
Meanwhile
\[ \nabla_{\Phi}\sqrt{S}=\Gamma_{\Psi\Phi}^{\Psi}\sqrt{S},\tag{8.152} \]
despite the fact that \(\varepsilon_{\Phi}\) and \(\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