2  Curvature and Torsion

We have seen how to describe curves and reparametrized them. Now we want to look at local properties of curves:

We will measure two quantities:

For example a 2D spiral is curved, but still lies in a plane. Instead the Helix both deviates from a straight line and pulls away from any fixed plane.

2.1 Curvature

We start with an informal discussion. Suppose γγ is a straight line γγ(t)=a+tv with a,vR3. The tangent vector to γγ is constant γγ˙(t)=v. Whatever the definition of curvature will be, it has to hold that γγ has zero curvature in this case. If we further derive the tangent vector, we obtain γγ¨(t)=00. Thus γγ¨ seems to be a good candidate for the definition of curvature of γγ at the point γγ(t).

Suppose now that γγ is a curve in R2 with unit speed. We have proven that in this case γγ˙γγ¨=0, that is, the vector γγ¨ is orthogonal to the tangent γγ˙ at all times. Now let n(t) be the unit vector orthogonal to γγ˙(t) at the point γγ(t). The amount that the curve γγ deviates from its tangent at γγ(t) after time t0 is (2.1)(γγ(t+t0)γγ(t))n(t), as seen in the figure below.

Amount that γγ deviates from tangent is (γγ(t+t0)γγ(t))n(t)

Equation (2.1) is what we take as measure of curvature. Since γγ˙(t)γγ¨(t)=0 and γγ˙(t)n(t)=0, we conclude that γγ¨(t) is parallel to n(t). Since n(t) is a unit vector, there exists a scalar κ(t) such that γγ¨(t)=κ(t)n(t). As n is unitary, we have κ(t)=γγ¨(t)

Now, approximate γγ at t with its second order Taylor polynomial: γγ(t+t0)=γγ(t)+γγ˙(t)t0+γγ¨(t)2t02+o(t0) where the remainder o(t0) is such that limt00 o(t0)t02=0. Therefore, discarding the remainder, γγ(t+t0)γγ(t)γγ˙(t)t0+γγ¨(t)2t02. Multiplying by n(t) we get (γγ(t+t0)γγ(t))n(t)γγ˙(t)n(t)t0+γγ¨(t)n(t)2t02. Recalling that γγ˙(t)n(t)=0,γγ¨(t)n(t)=κ(t), we then obtain (γγ(t+t0)γγ(t))n(t)12κ(t)t02

Important
The amount that γγ deviates from a straight line is proportional to κ(t)=γγ¨(t).

We take this as definition of curvature for a general unit speed curve in Rn.

Definition 1
Let γγ :(a,b)Rn be a unit speed curve. The curvature of γγ at γγ(t) is κγγ(t):=γγ¨(t).

Note that κ(t) is a function of time. Therefore the curvature of γγ can change from point to point.

We now define curvature for curves which are regular, but not necessarily unit speed.

Definition 2
Let γγ :(a,b)Rn be a regular. The curvature of γγ at γγ(t) is κγγ(t):=γγ~¨(ϕ(t)),t(a,b), where γγ~ is a unit speed reparametrization of γγ, with γγ=γγ~ϕ.

Remark 3

The above definition is well posed:

  • Since γγ is regular, there exist a unit speed reparametrization γγ~ of γγ.
  • If γγ^ is another unit speed reaprametrization of γγ, with γγ=γγ^ϕ^, then κγγ(t)=γγ^¨(ϕ^(t)), showing that there is no ambiguity in the definition of κγγ.

Indeed, since γγ~ and γγ^ are both reparametrizations of γγ, then γγ(t)=γγ~(ϕ~(t)),γγ(t)=γγ^(ϕ^(t)) for some diffeomorphisms ϕ~,ϕ^. Hence (2.2)γγ~(t)=γγ^(ϕ(t)),ϕ:=ϕ^(ϕ~)1, where ϕ is a diffeomorphism, since it is composition of diffeomorphisms. Differentiating (2.2) we get (2.3)γγ~˙(t)=γγ^˙(ϕ(t))ϕ˙(t). Taking the norms of the above, and recalling that γγ~ and γγ^ are unit speed, we get (2.4)|ϕ˙(t)|=1,t. Since ϕ is a diffeomorphism, we already know that |ϕ˙|0. As ϕ˙ is continuous, this means that the sign of ϕ˙ is constant. Thus (2.4) implies ϕ˙(t)1orϕ˙(t)1. In both cases, we have ϕ¨0. Differentiating (2.3) we then obtain γγ~¨(t)=γγ^¨(ϕ(t))ϕ˙2(t)+γγ^˙(ϕ(t))ϕ¨(t)=γγ^¨(ϕ(t))ϕ˙2(t). Taking the norms and using again that |ϕ˙|1, we get that γγ~¨(t)=γγ^¨(ϕ(t)). Recalling that ϕ=ϕ^(ϕ~)1 we get γγ~¨(ϕ~(t))=γγ^¨(ϕ^(t)),t(a,b). Therefore κγγ(t)=γγ~¨(ϕ~(t))=γγ^¨(ϕ^(t)).

Remark 4: Methods for computing curvature
In summary, the curvature of a regular curve γγ :(a,b)Rn is defined via unit speed reparametrizations of γγ. To compute κ we do the following:

  • We find a unit speed reparametrization γγ~ of the regular curve γγ
  • This can be done by computing s the arc-length of γγ, and then defining γγ~:=γγψ,ψ:=s1
  • Then we compute κγγ~(t)=γγ~¨(t)
  • We obtain the curvature of γγ by κγγ(t)=κγγ~(t)

When γγ is regular and has values in R3, there is a way to compute κ without reparametrizing. To do this, we will need the notion of cross product, or vector product. We will see this in the following sections.

Procedure for computing curvature κ

We conclude with two examples in which we compute the curvature κ using unit speed reparametrizations.

Example 5
Consider the circle of radius R>0: γγ(t)=(Rcos(t),Rsin(t)),t[0,2π]. To compute the curvature of γγ we need to find a unit speed reparametrization. We have shown that: γγ regular ϕ=s1 unit speed reparametrization where s is the arc length of γγ: s(t):=t0tγγ˙(τ)dτ. In our case γγ˙(t)=(Rsin(t),Rcos(t))γγ˙(t)=R and so γγ is regular. However γγ is not unit speed, therefore we need to find a unit speed reparametrization. The arc length starting at t0=0 is s(t)=0tRdτ=tR. The inverse of s is ϕ(t):=s1(t)=tR. Therefore a unit speed reparametrization of γγ is γγ~:=γγϕ which reads γγ~(t):=(Rcos(tR),Rsin(tR)). We have γγ~˙(t)=(sin(tR),cos(tR))γγ~¨(t)=(1Rcos(tR),1Rsin(tR)) Therefore the curvature of γγ is κ(t)=γγ~¨(t)=1R. In this case κ(t) is constant. The curvature also tells us that the smaller the circle, the higher the curvature. For a large circle, like the Earth, the curvature is barely noticeable.

Before proceeding with the next example, let us give a short overview of the Hyperbolic functions.

Remark 6: Hyperbolic functions

The Hyperbloic functions are the analogous of the trigonometric functions, but defined using the hyperbola rather than the circle. Their formulas can be obtained by means of the exponential function et. We have:

  • Hyperbolic cosine: The even part of the function et, that is, cosh(t)=et+et2=e2t+12et=1+e2t2et.

  • Hyperbolic sine: The odd part of the function et, that is, sinh(t)=etet2=e2t12et=1e2t2et.

  • Hyperbolic tangent: Defined by tanh(t)=sinhtcosht=etetet+et=e2t1e2t+1.

  • Hyperbolic cotangent: The reciprocal of tanh for t0, cotht=coshtsinht=et+etetet=e2t+1e2t1.

  • Hyperbolic secant: The reciprocal of cosh sech(t)=1cosht=2et+et=2ete2t+1.

  • Hyperbolic cosecant: The reciprocal of sinh for t0, csch(t)=1sinht=2etet=2ete2t1.

For a plot cosh,sinh,tanh see Figure 2.1 below. The properties of the hyperbolic functions which are of interest to us are:

  1. Identities: cosh(t)+sinh(t)=etcosh(t)sinh(t)=etcosh2(t)sinh2(t)=1sech2(t)tanh2(t)=1

  2. Derivatives: ddt[sinh(t)]=cosh(t)ddt[cosh(t)]=sinh(t)ddt[tanh(t)]=1tanh2(t)=csch2(t)

  3. Integrals: t0tsinh(u)du=cosh(t)cosh(t0)t0tcosh(u)du=sinh(t)sinh(t0)t0ttanh(u)du=log(cosh(t))log(cosh(t0))

Figure 2.1: Plot of cosh,sinh,tanh.

Example 7: The Catenary
The catenary is the shape of a heavy chain suspended at its ends. The chain is only subjected to gravity, see Figure 2.2. This shape looks similar to a parabola, but it is not a parabola. This was first noted by Galilei, see this Wikipedia page. The profile of the hanging chain can be obtained via a minimization problem, and one can show it is of the form γγ(t)=(t,cosh(t)),tR. See Figure 2.3 for a plot of γγ. Let us check if γγ is regular. We have γγ˙(t)=(1,sinh(t)) so that γγ˙2=1+sinh2(t)=cosh2(t)γγ˙=cosh(t). Note that cosh(t)1 showing that γγ is regular. However γγ˙(1)=cosh(1)=e+e121.54, proving that γγ is not unit speed. Let us then compute the arc length of γγ starting at t0=0 s(t)=0tγγ˙(u)du=0tcosh(u)du=sinh(t) since sinh(0)=0. We need to invert s. We have s=sinh(t)s=etet2e2t2set1=0, where the last equation was obtained multuplying both sides by et. Now we substitute y=et and obtain e2t2set1=0y22sy1=0y=s±1+s2. Recalling that y=et, we only consider the positive solution, and obtain that et=s+1+s2t=log(s+1+s2). We have proven that the inverse of the arc length s(t) is ψ(t):=s1(t)=log(t+1+t2). Therefore γγ~(t):=γγ(ψ(t)) is a unit speed reparametrization of γγ. Substituting ψ and using the definition of γγ we have γγ~(t)=(log(t+1+t2),1+t2). We can now compute the curvature. We have: γγ~˙(t)=(11+t2,t1+t2)γγ~¨(t)=(t(1+t2)3/2,1(1+t2)3/2) Moreover γγ~¨(t)2=t2(1+t2)3+1(1+t2)3=1(1+t2)2. Therefore the curvature is κ(t)=γγ~¨(t)=11+t2.

Figure 2.2: The catenary is the shape of a heavy chain suspended at its ends. Image from Wikipedia.

Figure 2.3: Plot of the catenary curve γγ(t)=(t,cosh(t)).

2.2 Vector product in R3

The discussion in this section follows (). We start by defining orientation for a vector space.

Definition 8: Same orientation
Consider two ordered basis of R3 B=(b1,b2,b3),B~=(b~1,b~2,b~3). We say that B and B~ have the same orientation if the matrix of change of basis has positive determinant.

When two basis B and B~ have the same orientation, we write bb~. The above is clearly an equivalence relation on the set of ordered basis. Therefore the set of ordered basis of R3 can be decomposed into equivalence classes. Since the determinant of the matrix of change of basis can only be positive or negative, there are only two equivalence classes.

Definition 9: Orientation
The two equivalence classes determined by on the set of ordered basis are called orientations.

Definition 10: Positive orientation

Consider the standard basis of R3 E=(e1,e2,e3) where we set e1=(1,0,0),e2=(0,1,0),e3=(0,0,1). Then:

  • The orientation corresponding to E is called positive orientation of R3.
  • The orientation corresponding to the other equivalence class is called negative orientation of R3.

For a basis B of R3 we say that:

  • B is a positive basis if it belongs to the class of e.
  • B is a negative basis if it does not belong to the class of e.

Example 11
Since we are dealing with ordered basis, the order in which vectors appear is fundamental. For example, we defined the equivalence class of E=(e1,e2,e3), to be the positive orientation of R3. In particular e is a positive basis.
Consider instead E~=(e2,e1,e3). The matrix of change of variables between E~ and E is (e2|e1|e3)=(010100001) and the latter has negative determinant. Thus E~ does not belong to the class of E, and is therefore a negative basis.

We are now ready to define the vector product in R3.

Definition 12: Vector product in R3
Let u,vR3. The vector product of u and v is the unique vector u×vR3 which satisfies the property: (2.5)(u×v)w=|u1u2u3v1v2v3w1w2w3|,wR3. Here |aij| denotes the determinant of the matrix (aij), and u=i=13uiei,v=i=13viei,w=i=13wiei, with (e1,e2,e3) standard basis of R3.

The following proposition gives an explicit formula for computing u×v.

Proposition 13
Let u,vR3. Then (2.6)u×v=|u2u3v2v3|e1|u1u3v1v3|e2+|u1u2v1v2|e3.

Proof
Denote by (u×v)i the i-th component of u×v with respect to the standard basis, that is, u×v=i=13(u×v)iei. We can use (2.5) with w=e1 to obtain (u×v)e1=|u1u2u3v1v2v3100|=|u2u3v2v3| where we used the Laplace expansion for computing the determinant of the 3×3 matrix. As the standard basis is orthonormal, by bilinearity of the scalar product we get (u×v)e1=i=13(u×v)ieie1=(u×v)i. Therefore we have shown (u×v)1=|u2u3v2v3|. Similarly we obtain (u×v)2=|u1u2u3v1v2v3010|=|u1u3v1v3| and (u×v)3=|u1u2u3v1v2v3001|=|u1u2v1v2|, from which we conclude.

Sometimes we will denote formula (2.6) by u×v=|ijku1u2u2v1v2v3|.

Let us collect some crucial properties of the vector product.

Proposition 14

The vector product in R3 satisfies the following properties: For all u,vR3

  1. u×v=v×u
  2. u×v=00 if and only if u and v are linearly dependent
  3. (u×v)u=0, (u×v)v=0
  4. For all wR3, a,bR (au+bw)×v=au×v+bw×w

The proof, which is based on the properties of determinants, is omitted.

Remark 15: Geometric interpretation of vector product

Let u,vR3 be linearly independent. We make some observations:

  1. Property 3 in Proposition 14 says that (u×v)u=0,(u×v)v=0. Therefore u×v is orthogonal to both u and v.

  2. In particular u×v is orthogonal to the plane generated by u and v.

  3. Since u and v are linearly independent, Property 2 in Proposition 14 says that u×v00

  4. Therefore we have (u×v)(u×v)=u×v2>0

  5. On the other hand, using the definition of u×v with w=v×w yields (u×v)(u×v)=|u1u2u3v1v2v3(u×v)1(u×v)2(u×v)3|

  6. Therefore the determinant of the matrix (u|v|u×v) is positive. This shows that (u,v,u×v) is a positive basis of R3.

  7. For all u,v,x,yR3 it holds (2.7)(u×v)(x×y)=|uxvxuyvy|. Indeed, one can check that the above formula holds for the standard vectors ei, and thus the general formula follows by linearity.

  8. Using (2.7) we get u×v2=(u×v)(u×v)=|uuvuuvvv|=u2v2|uv|2=u2v2u2v2cos2(θ)=u2v2(1cos2(θ))=u2v2sin2(θ)=A2 where A is the area of the parallelogram with sides u and v.

For u,v linearly independent, u×v is orthogonal to the plane generated by u,v. Moreover |u×v| is the area of the parallelogram with sides u,v, and (u,v,u×v) is a positive basis of R3

Let us summarize the above remark.

Remark 16: Summary: Properties of u×v

Let u,vR3 be linearly independent. Then

  • u×v is orthogonal to the plane spanned by u,v
  • u×v is equal to the area of the parallelogram with sides u,v
  • u×v is such that (u,v,u×v) is a positive basis of R3.

We conclude with noting that the cross product is not associative, and with a useful proposition for differentiating the cross product of curves in R3.

Proposition 17
The vector product is not associative. In particular, for all u,v,wR3 it holds: (2.8)(u×v)×w=(uw)v(vw)u.

The proof is omitted. It follows by observing that both sides of (2.8) are linear in u,v,w. Therefore it is sufficient to verify (2.8) for the standard basis vectors ei. This is left as an exercise.

Proposition 18
Suppose γγ,ηη :(a,b)R3 are parametrized curves. Then the curve γγ×ηη :(a,b)R3 is smooth, and (2.9)ddt(γγ×ηη)=γγ˙×ηη+γγ×ηη˙.

The proof is omitted. It follows immediately from formula (2.6).

2.3 Curvature formula in R3

Given a unit speed curve γγ :(a,b)Rn we defined its curvature as κ(t)=γγ¨(t). If γγ is not unit speed then the curvature is not defined. However, when γγ is regular, then we can find a unit-speed reparametrization γγ~ of γγ, and compute κ as κ(t)=γγ~¨(t). If γγ is a regular curve in R3, there is a way to compute κ without passing through γγ~. The formula for computing κ is as follows.

Proposition 19: Curvature formula
Let γγ:(a,b)R3 be a regular curve. The curvature κ(t) of γγ at γγ(t) is given by (2.10)κ(t)=γγ˙×γγ¨γγ˙3.

We delay the proof of the above Proposition, as this will get easier when the Frenet frame is introduced. For a proof which does not make use of the Frenet frame, see the proof of Proposition 2.1.2 in ().

For now we use (2.10) the above proposition to compute the curvature on specific curves.

Example 20
Consider the straight line γγ(t)=a+tv for some a,vR3 fixed, with v00. Then γγ˙(t)=v,γγ¨(t)=00. Therefore γγ˙(t)=v0 showing that γγ is regular. We have γγ˙×γγ¨=v×00=00. Therefore the curvature is κ=γγ˙×γγ¨γγ˙3=0, as expected.

Example 21

Consider the Helix of radius R>0 and rise H>0 γγ(t)=(Rcos(t),Rsin(t),Ht),tR. Then γγ˙(t)=(Rsin(t),Rcos(t),H)γγ¨(t)=(Rcos(t),Rsin(t),0) From this we deduce that γγ˙(t)=R2+H2, showing that γγ is regular. Finally γγ˙×γγ¨=|γ˙2γ˙3γ¨2γ¨3|e1|γ˙1γ˙3γ¨1γ¨3|e2+|γ˙1γ˙2γ¨1γ¨2|e3=|Rcos(t)HRsin(t)0|e1|Rsin(t)HRcos(t)0|e2+|Rsin(t)Rcos(t)Rcos(t)Rsin(t)|e3=(RHsin(t),RHcos(t),R2cos2(t)+R2sin2(t))=(RHsin(t),RHcos(t),R2) and therefore γγ˙×γγ¨=RR2+H2. By the general formula we have κ=γγ˙×γγ¨γγ˙3=R(R2+H2)12(R2+H2)32=RR2+H2

We notice the following:

  • If H=0 then the Helix is just a circle of radius R. In this case the curvature is κ=1R which agrees with the curvature computed for the circle of radius R.

  • If R=0 then the Helix is just parametrizing the z-axis. In this case the curvature is κ=0, which agrees with the curvature of a straight line.

2.4 Signed curvature of plane curves

In this section we assume to have plane curves, that is, curves with values in R2. In this case we can give a geometric interpretation for the sign of the curvature. This cannot be done in higher dimension.

Definition 22
Let γγ :(a,b)R2 be unit speed. We define the signed unit normal to γγ at γγ(t) as the unit vector n(t) obtained by rotating γγ˙(t) anti-clockwise by an angle of π/2.

Definition 23
Let γγ :(a,b)R2 be unit speed. The signed curvature of γγ at γγ(t) is the scalar κs(t) such that γγ¨(t)=ks(t)n(t)

Remark 24
Notice that since n is a unit vector and γγ is unit speed, then |κs(t)|=γγ¨(t)=κ(t). Thus the signed curvature is related to the curvature by κs(t)=±κ(t).

Remark 25
It can be shown that the signed curvature is the rate at which the tangent vector γγ˙ of the curve γγ rotates. The signed curvature is:

  • positive if γγ˙ is rotating anti-clockwise
  • negative if γγ˙ is rotating clockwise

In other words,

  • ks>0 means the curve is turning left,
  • ks<0 means the curve is turning right.

A rigorous justification of the above statement is found in Proposition 2.2.3 in ().

For curves which are not unit speed, we define the signed curvature as the signed curvature of the unit speed reparametrization.

Definition 26
Let γγ :(a,b)R2 be regular and let γγ~ be a unit speed reparametrization of γγ. The signed curvature of γγ at γγ(t) is the scalar κs(t) such that γγ~¨(t)=ks(t)n(t), where n(t) is the unit vector obtained by rotating γγ~˙(t) anti-clockwise by an angle π/2.

The signed curvature completely characterizes plane curves, in the sense of the following theorem.

Theorem 27: Characterization of plane curves

Let ϕ :RR be smooth. Then:

  1. There exists a unit speed curve γγ :RR2 such that its signed curvature κs satisfies κs(t)=ϕ(t),tR.

  2. Suppose that γγ~ :RR2 is a unit speed curve such that its signed curvature κ~s satisfies κ~s(t)=ϕ(t),tR. Then γγ~=γγ up to rotations and translations.

We do not prove the above theorem. For a proof, see Theorem 2.2.6 in ().

2.5 Space curves

In this section we deal with space curves, that is, curves with values in R3. There are several issues compare to the plane case:

  • A 3D counterpart of the signed curvature does not exist, since there is no notion of turning left or turning right.

  • We have seen in the previous section that the signed curvature completely characterizes plane curves. In 3D however curvature is not enough to characterize curves: there exist γγ and ηη space curves such that κγγ=κηη,γγηη, that is, γγ and ηη have same curvature but are different curves.

Example 28
Let γγ be a circle of radius R>0 γγ(t)=(Rcos(t),Rsin(t),0), and ηη be a helix of radius S>0 and rise H>0 ηη(t)=(Scos(t),Ssin(t),Ht). We have computed that κγγ=1R,κηη=SS2+H2. If we now choose R=2 and we impose that κγγ=κηη we get 1R=SS2+H2H2=2SS2 Therefore choosing S=1 and H=1 yields κγγ=κηη,γγηη..

Therefore curvature is not enough for characterizing space curves, and we need a new quantity. As we did with curvature, we start by considering the simpler case of unit speed curves. We will also need to assume that the curvature is never zero.

Definition 29: Principal normal vector
Let γγ:(a,b)R3 be a unit speed curve with κ(t)0,t(a,b). The principal normal vector to γγ at γγ(t) is n(t):=1κ(t)γγ¨(t).

Remark 30

Since for γγ unit speed we defined κ(t):=γγ¨(t), we have that n(t)=1, thus n is a unit vector. Moreover n is orthogonal to γγ˙, that is, γγ˙n=0.

This is because γγ˙n=1κγγ˙γγ¨=0, where the last equality follows from γγ˙γγ¨=0, being γγ unit speed.

Principal normal vector n(t) to γγ at γγ(t).

Question 31
Why is the principal normal interesting? Because it can tell the difference beween a plane curve and a space curve. See picture below.

Left: Principal normal to a circle. Note that n always points towards the origin 00. Right: Principal normal to a helix. Note that n points towards the z-axis, but never towards the same point.

Definition 32: Binormal vector
Let γγ:(a,b)R3 be a unit speed curve with κ(t)0,t(a,b). The binormal vector to γγ at γγ(t) is b(t):=γγ˙(t)×n(t).

Definition 33: Orthonormal basis
Let v1,v2,v3 be vectors in R3. We say that the triple {v1,v2,v3} is orthonormal if vi=1,vivj=0, for ij.

Proposition 34
Let γγ:(a,b)R3 be a unit speed curve with κ(t)0,t(a,b). Then the triple B=(γγ˙(t),n(t),b(t)) is a positive orthonormal basis of R3 for all t(a,b).

Proof
Since γγ is unit speed we have γγ˙(t)1. Moreover we have already observed that n(t)1,γγ˙(t)n(t)0. As b is defined by b:=γγ˙×n, by the properties of the vector product, see Proposition 14, it follows that bγγ˙=0,bn=0. By the calculation in Remark 15 Point 8, we have that b2=γγ˙2n2|γγ˙n|2=1. This shows that the vectors {γγ˙,n,b} are orthonormal. By the properties of the vector product, see Remark 15 Point 6, we also know that (γγ˙,n,b) is a positive basis of R3.

Proposition 35
Let γγ:(a,b)R3 be a unit speed curve with κ0. Then (2.11)b=γγ˙×n,n=b×γγ˙,γγ˙=n×b.

Proof
The first equality in (2.11) is true by definition of b. For the other 2 equalities, recall formula (2.8): (2.12)(u×v)×w=(uw)v(vw)u, for all u,v,wR3. Applying the above with u=γγ˙,v=n,w=γγ˙, yields (γγ˙×n)×γγ˙=(γγ˙γγ˙)n(nγγ˙)γγ˙=γγ˙2n0=n, where we used that γγ˙ is a unit vector and nγγ˙=0. Therefore, by definition of b, we have b×γγ˙=(γγ˙×n)×γγ˙=n showing the second equality in (2.11). For showing the third equality in (2.11), we apply (2.12) with u=γγ˙,v=n,w=n, to get (γγ˙×n)×n=(γγ˙n)n(nn)γγ˙=0n2γγ˙=γγ˙ where we used that n is a unit vector and γγ˙n=0. Therefore, by definition of b and anti-commutativity of the vector product, we have n×b=b×n=(γγ˙×n)×n=γγ˙, showing the last equality in (2.11).

Proposition 36
Let γγ:(a,b)R3 be a unit speed curve with κ0. Then (2.13)b˙(t)=τ(t)n(t), for some τ(t)R.

Proof
By definition of b and the formula of derivation of the cross product (2.9) we have b˙=ddt(γγ˙×n)=γγ¨×n+γγ˙×n˙=γγ˙×n˙, where we used that γγ¨×n=0, since n is defined by n:=γγ¨/κ, and therefore n and γγ¨ are parallel. Hence, we have proven that (2.14)b˙=γγ˙×n˙. By the properties of the cross product we have that u×v is orthogonal to both u and v. Thus (2.14) implies that b˙γγ˙=0. Further, observe that ddt(bb)=b˙b+bb˙=2b˙b. On the other hand, since b is a unit vector, we have ddt(bb)=ddt(b2)=ddt(1)=0 Therefore b˙b=0. To summarize, we have shown that b˙ is orthogonal to b and γγ˙. Since (γγ˙,n,b) is an orthonormal basis of R3 we conclude that b˙ is parallel to n. Therefore there exists τ(t)R such that b˙=τ(t)n(t), concluding the proof.

The scalar τ in equation (2.13) is called the torsion of γγ.

Definition 37: Torsion of unit speed curve
Let γγ:(a,b)R3 be a unit speed curve, with κ0. The torsion of γγ at γγ(t) is the unique scalar τ(t)R such that b˙(t)=τ(t)n(t).

Remark 38

In particular the torsion satisfies: τ(t)=b˙(t)n(t).

The above can be immediately obtained by multiplying (2.13) by n. Indeed, b˙=τnb˙n=τnn=τ, since n is a unit vector.

Warning
We defined the torsion only for space curves γγ:(a,b)R3 which are unit speed and have non-vanishing curvature, that is, such that γγ˙(t)=1,κ(t)=γγ¨(t)0, for all t(a,b).

We can extend the definition of torsion to regular curves γγ with non-vanishing curvature. In this case the torsion of γγ is defined as the torsion of a unit speed reparametrization of γγ.

Definition 39
Let γγ:(a,b)R3 be a regular curve with non-vanishing curvature. Let γγ~ be a unit speed reparametrization of γγ, with γγ=γγ~ϕ,ϕ:(a,b)(a~,b~). We define the torsion of γγ at γγ(t) as τγγ(t):=τγγ~(ϕ(t)), where τγγ~(s) denotes the torsion of γγ~ at γγ~(s).

As usual, it is possible to check that the above definition of torsion does not depend on the choice of unit speed reparametrization γγ~. As with curvature, there is a general formula to compute the torsion without having to reparametrize.

Proposition 40: Torsion formula
Let γγ:(a,b)R3 be a regular curve with non-vanishing curvature. The torsion τ(t) of γγ at γγ(t) is given by τ(t)=(γγ˙×γγ¨)γγγγ˙×γγ¨2.

We delay the proof of the above proposition for a bit. In the meantime, let us look at examples.

Example 41: Torsion Helix
Consider the Helix of radius R>0 and rise H>0 γγ(t)=(Rcos(t),Rsin(t),Ht),tR. We have already shown that γγ˙(t)=R2+H2,κ=RR2+H2. Therefore the Helix is regular with non-vanishing curvature. The torsion can be then computed via the formula τ(t)=(γγ˙×γγ¨)γγγγ˙×γγ¨2. Let us compute the quantities appearing in the formula for τ γγ˙(t)=(Rsin(t),Rcos(t),H)γγ¨(t)=(Rcos(t),Rsin(t),0)γγ(t)=(Rsin(t),Rcos(t),0) Moreover we had already computed that γγ˙×γγ¨=(RHsin(t),RHcos(t),R2)γγ˙×γγ¨=RR2+H2. Finally we compute (γγ˙×γγ¨)γγ=R2H. We are ready to find the torsion: τ=(γγ˙×γγ¨)γγγγ˙×γγ¨2=HR2+H2.

Example 42: Curvature and Torsion of Circle
The Circle of radius R>0 is γγ(t):=(Rcos(t),Rsin(t),0). The curvature and torsion of the Helix of radius R and rise H>0 are κ=RR2+H2,τ=HR2+H2. For H=0 the Helix coincides with the Circle γγ. Therefore we can set H=0 in the above formulas to obtain the curvature and torsion of the Circle κ=1R,τ=0.

From the above example we notice that the torsion of the circle is 0. This is true in general for space curves which are contained in a plane: we will prove this result in general. For the moment, let us give an example for which this happens, that is, an example of space curve γγ which is contained in a plane.

Example 43
Define the space curve γγ(t):=(45cos(t),1sin(t),35cos(t)), for tR. As seen in the plot below, γγ is just a Circle which has been rotated an translated. Therefore γγ is contained in a plane, and we expect curvature and torsion to be κ=1R,τ=0, for some R>0, radius of the Circle γγ. Let us proceed with the calculations: γγ˙=(45sin(t),cos(t),35sin(t)) so that γγ˙2=1625sin2(t)+cos2(t)+925sin2(t)=1, showing that γγ is regular and unit speed. Further γγ¨=(45cos(t),sin(t),35cos(t)). As γγ is unit speed, we have κ=γγ¨=1625cos2(t)+sin2(t)+925cos2(t)=1. As γγ is unit speed, the normal vector is n=1κγγ¨=(45cos(t),sin(t),35cos(t)). We can then compute the binormal b=γγ˙×n=|ijk45sin(t)cos(t)35sin(t)45cos(t)sin(t)35cos(t)|=(35cos2(t)35sin2(t),1225cos(t)sin(t)+1225cos(t)sin(t),45sin2(t)45cos2(t))=(35,0,45). Therefore b˙=0, and we obtain that the torsion is τ=b˙n=0.

Plot of the curve in example above

2.6 Frenet frame

For a unit speed curve γγ:(a,b)R3 with non-vanishing curvature we computed the triple {γγ˙,n,b}. We saw that the above is a positive orthonormal basis of R3. We also used this triple to compute curvature κ and torsion τ of γγ: κ=γγ¨,τ=b˙n. This triple is so important that it has a name.

Definition 44: Frenet frame
Let γγ:(a,b)R3 be unit speed with κ0. The positive orthonormal basis {γγ˙,n,b} is called Frenet frame of γγ.

We can also define the Frenet frame for regular curves with non-vanishing curvature.

Definition 45
Let γγ:(a,b)R3 be regular with κ0. The Frenet frame of γγ is defined as the Frenet frame of a unit speed reparametrization γγ~ of γγ.

Remark 46

We should check that the above definition is well-posed:

  • Note that γγ~ is unit speed. Moreover the curvature of κγγ~ is given by κγγ~(t)=κγγ(ϕ(t)) for some ϕ diffeomorphism. Therefore κγγ~0 as we are assuming κγγ0. Therefore the Frenet-Frame of γγ~ is well defined.

  • If γγ^ is another unit speed reparametrization of γγ, then the Frenet frame generated by γγ^ coincides with the one generated by γγ~. The proof is left as an exercise.

From the Frenet frame we can define the Frenet-Serret equations.

Theorem 47: Frenet-Serret equations
Let γγ:(a,b)R3 be unit speed with κ0. The Frenet-Serret equations are γγ¨=κnn˙=κγγ˙+τbb˙=τn

Proof
The first Frenet-Serret equation (2.15)γγ¨=κn holds by definition of n and κ. The third Frenet-Serret equation (2.16)b˙=τn holds by Proposition 36. Now, recall that in Proposition 35 we have proven (2.17)b=γγ˙×n,n=b×γγ˙,γγ˙=n×b. Differentiating the second equation in (2.17) and using (2.15)-(2.16) we get n˙=b˙×γγ˙+b×γγ¨=(τn×γγ˙)+b×κn=τ(γγ˙×n)κ(n×b)=τbκγγ˙, where in the last equality we used the first and third equations in (2.17). The above is exactly the second Frenet-Serret equation.

Remark 48
We can write the Frenet-Serret ODE sysyem in vectorial form. To this end, introduce the matrix F:=(0κ0κ0τ0τ0). It is immediate to check that the Frenet-Serret equations are equivalent to (γγ¨n˙b˙)=F(γγ˙nb).

Important: Summary
Recall that:

  1. Curvature κ is defined only for regular curves.
  2. Torsion τ is defined only for regular curves with non-vanishing κ.

The two strategies for computing κ and τ are discussed in the diagram in Figure 2.4 below.

Figure 2.4: Summary for computing κ and τ for regular curve γγ.

Let us conclude the section with an example. We compute the Frenet frame of the helix. As a consequence we obtain curvature and torsion.

Example 49: Frenet frame of helix
Consider the helix of radius 1 and rise 1 given by γγ(t)=(cos(t),sin(t),t), for tR. We now proceed following the diagram at Figure 2.4. We ask the first question: Is γγ unit speed? We have that γγ˙(t)=(sin(t),cos(t),1), and therefore γγ˙=2. This shows that γγ is regular but not unit speed. We ask the second question in the diagram: Can we find a unit speed reparametrization of γγ? Let us try. We compute the arc length of γγ starting at t0=0 s(t):=0tγγ˙(u)du=2t. The arc length is invertible with ψ(t):=s1(t)=t2. Therefore a unit speed reparametrization of γγ is given by γγ~(t):=γγ(ψ(t))=(cos(t2),sin(t2),t2). The next step in the diagram is Compute Frenet frame {γγ˙,n,b} and curvature κ, torsion τ We compute γγ~˙(t)=12(sin(t2),cos(t2),1)γγ~¨(t)=12(cos(t2),sin(t2),0) Therefore the curvature is κ(t)=γγ~¨(t)=12. From the curvature we obtain the principal normal vector n(t)=1κ(t)γγ~¨(t)=(cos(t2),sin(t2),0). We can now compute the binormal b(t)=γγ~˙×n=12|ijksin(t2)cos(t2)1cos(t2)sin(t2)0|=12(sin(t2),cos(t2),1). We have therefore computed the Frenet frame of γγ. This is given by γγ~˙(t)=12(sin(t2),cos(t2),1)n(t)=(cos(t2),sin(t2),0)b(t)=12(sin(t2),cos(t2),1). See below for a picture of the Frenet frame of the helix. Given the Frenet frame, we can compute the torsion via the formula τ(t)=b˙n. Indeed, we have b˙=12(cos(t2),sin(t2),0) and therefore b˙n=12(cos2(t2)sin2(t2))=12. The torsion is then τ(t)=b˙n=12. The Frenet-Frame of the unit-speed Helix is plotted in Figure 2.5.

Figure 2.5: Frenet frame of the helix of radius 1 and rise 1.

2.7 Consequences of Frenet-Serret

The most important consequence of the Frenet-Serret equations is that they allow to fully characterize space curves in terms of curvature and torsion. Precisely, the following theorem holds.

Theorem 50: Characterization of space curves

Let κ,τ :RR be smooth functions, with κ>0. Then:

  1. There exists aunit speed curve γγ :RR3 such that its curvature κγγ and torsion τγγ satisfy κγγ(t)=κ(t),τγγ(t)=τ(t),tR.

  2. Suppose that γγ~ :RR3 is a unit speed curve such that its curvature κ~γγ~ and torsion τγγ~ satisfy κγγ~(t)=κ(t),τγγ~(t)=τ(t),tR. Then γγ~=γγ up to rotations and translations.

The proof of Theorem 50 is omitted, and it can be found in Theorem 2.3.6 in ().

Theorem 50 is a very strong result. It is saying two things:

  1. If we prescribe curvature and torsion, then there exists a unit speed curve which has such curvature and torsion.

  2. If two unit speed curves have same curvature and torsion, then they must be the same curve, up to translations and rotations.

In other words, curvature and torsion fully characterize space curves. This result is the 3D counterpart of Theorem 27, which said that signed curvature characterizes 2D curves.

Example 51
In Example 43 we have considered the unit speed curve γγ(t):=(45cos(t),1sin(t),35cos(t)), for t[0,2π]. We have computed that κγγ=1,τγγ=0. If we plot γγ, we clearly see that γγ is just obtained by translating and rotating a unit circle, see plot below. Theorem 50 enables us to rigorously prove this claim. Indeed, consider the unit speed circle γγ~(t):=(cos(t),sin(t),0), for t[0,2π]. In Example 42 we have proven that curvature and torsion are
κγγ~=1,τγγ~=1. Therefore κγγ=κγγ~,τγγ=τγγ~, and by Theorem 50 we conclude that γγ is equal to γγ~ up to rotations and translations.

Plot of the curve in example above

Another consequence of the Frenet-Serret equations is that they allow us to finally prove the curvature and torsion formulas given in Proposition 19 and Proposition 40. For reader’s convenience we recall these two results.

Proposition 52: Curvature and torsion formulas
Let γγ:(a,b)R3 be a regular curve. The curvature κ(t) of γγ at γγ(t) is given by κ(t)=γγ˙×γγ¨γγ˙3. Suppose in addition that γγ has non-vanishing curvature. The torsion τ(t) of γγ at γγ(t) is given by τ(t)=(γγ˙×γγ¨)γγγγ˙×γγ¨2.

Before proceeding with the proof, we need to establish some notation.

Notation: Compact notation for arc length reparametrization

Suppose γγ:(a,b)Rn is regular and denote by s:(a,b)(a~,b~),ts(t) its arc length. We already know that in this case s invertible, with inverse s1 giving a unit speed reparametrization γγ~:(a~,b~)Rn of γγ, defined by γγ~=γγψ,ψ:=s1:(a~,b~)(a,b) Sometimes it is more convenient to adopt more compact notation. In the new notation the unit speed reparametrization γγ~ is denoted by γγ(s): tγγ~(t)sγγ(s). Thus, the reparametrization is denoted with the same symbol γγ, but this time γγ is considered as a function of the arc length parameter s(a~,b~). We will denote:

  • The derivative of s by dsdt

  • The derivative of ψ=s1 by dtds.

Moreover:

  • The derivative of γγ(t) is denoted by dγγdt(t)=γγ˙(t),t(a,b)

  • The derivative of γγ(s) is denoted by dγγds(s)=γγ˙(s),s(a~,b~).

We also have new notations for the chain rule:

  • The chain rule for γγ is the old notations is: γγ(t)=γγ~(s(t))γγ˙(t)=γγ~˙(s(t))s˙(t),t(a,b). In the new notations the above chain rule is written dγγdt(t)=dγγds(s(t))dsdt(t),t(a,b). We will often omit the dependence on the point t by writing dγγds=dγγdtdtds.

  • The chain rule for the reparametrization γγ~ in the old notation is: γγ~(t)=γγ(ψ(t))γγ~˙(t)=γγ˙(ψ(t))ψ˙(t),t(a~,b~). In the new notations the above chain rule is written dγγds(s)=dγγdt(ψ(s))dtds(s),s(a~,b~), since ψ˙ is written dt/ds in the new notations. Without dependence on the point s, the above reads dγγds=dγγdtdtds.

Example 53: How to use the new notations
Let γγ and γγ~ be as above. We know that γγ~ is unit speed. Thus γγ(s) is unit speed with respect to s, that is, (2.18)γγ˙(s)=1,s(a~,b~). As an exercise, let us check that (2.18) holds, using the new notations. By chain rule we have γγ˙(s)=dγγds(s)=dγγdt(ψ(s))|dtds(s)|=γγ˙(ψ(s))|dtds(s)|. Now, recall that (2.19)dsdt(t)=s˙(t)=γγ˙(t),t(a,b). According to the new notations and the inverse function theorem, dtds(s)=1(dsdt(ψ(s)))=1γγ˙(ψ(s)),s(a~,b~), where we used (2.19) evaluated at t=ψ(s). Thus γγ˙(s)=γγ˙(ψ(s))|dtds(s)|=γγ˙(ψ(s))1γγ˙(ψ(s))=1, concluding (2.18).

Let us highlight the main feature of the above notation.

Important: New Notation!

Let γγ:(a,b)Rn be a regular curve:

  1. We denote by tγγ(t),t(a,b) the given curve γγ.

  2. We denote by sγγ(s),s(a~,b~) the arc length reparametrization of the curve γγ. The parameter s is the arc length parameter. In particular γγ(s) is unit speed with respect to s.

We will heavily rely on the new notations for proving Proposition 52.

Proof: Proof of Proposition 52
We only prove the formula for κ, as the one for τ can be obtained similarly, just with more calculations. For a proof see Proposition 2.3.1 in ().

Since γγ is regular, we can reparametrize γγ by arc length s(t). We denote the arc lenght reparametrization by γγ(s). We know that γγ(s) is unit speed, that is, dγγds=1. Therefore is well define the Frenet frame {t(s),n(s),b(s)},t(s):=γγ˙(s)=dγγds(s). The Frenet-Serret equations are t˙(s)=κ(s)n(s)n˙(s)=κ(s)t(s)+τ(s)b(s)b˙(s)=τ(s)n(s) By chain rule dγγdt=dγγdsdsdt=(dsdt)t. Differentiating the above we infer d2γγdt2=ddt[(dsdt)t]=d2sdt2t+(dsdt)dtdt. By chain rule we have dtdt=dtdsdtds, and therefore d2γγdt2=d2sdt2t+(dsdt)dtdt=d2sdt2t+(dsdt)2dtds. Hence γγ˙(t)×γγ¨(t)=dγγdt×d2γγdt2=(dsdt)t×[d2sdt2t+(dsdt)2dtds]=[(dsdt)(d2sdt2)t×t]+[(dsdt)3t×dtds]=(dsdt)3t×dtds, since t×t=0 by the properties of the cross product. Now we recall that dtds=κ(s)n(s) by the first Frenet-Serret equation. Moreover dsdt(t)=γγ˙(t)2. Therefore γγ˙(t)×γγ¨(t)=(dsdt)3t×dtds=γγ˙(t)3κ(s(t))t×n=γγ˙(t)3κ(s(t))b, where in the last line we used the definition of b b(s)=γγ˙(s)×n(s)=t(s)×n(s). We can now take the norms and obtain γγ˙(t)×γγ¨(t)=γγ˙(t)3κ(s(t))b=γγ˙(t)3κ(s(t)) using that b=1. As γγ is regular, we can divide by γγ˙(t)3 and obtain κ(s(t))=γγ˙(t)×γγ¨(t)γγ˙(t)3. Recalling that the curvature of γγ at t is defined as the curvature of γγ(s) at s(t), we conclude that the above is the desired formula.

We now state and prove two more results which directly follow from the Frenet-Serret equations. They state, respectivley:

  1. A curve has torsion τ=0 if and only if it is contained in a plane.

  2. A curve has constant curvature and zero torsion if and only if it is part of a circle.

Before proceeding, we recall the following.

Remark 54: Equation of a plane
The general equation of a plane ππd in R3 is given by

ππd={xR3:xP=d}, for some vector PR3 and scalar dR. Note that:

  • If d=0, the condition xP=0 is saying that the plane ππ0 contains all the points x in R3 which are orthogonal to P. In particular ππ0 contains the origin 00.

  • If d0, then ππd is the translation of ππ0 by the quantity d in direction P.

In both cases, P is the normal vector to the plane, as shown in Figure 2.6 below.

Figure 2.6: The plane ππ0 is the set of points of R3 orthogonal to P. The plane ππd is obtained by translating ππ0 by a quantity d in direction P.
Proposition 55

Let γγ:(a,b)R3 be regular and such that κ0. They are equivalent:

  1. The torsion of γγ satisfies τ(t)=0 for all t(a,b).

  2. The image of γγ is contained in a plane, that is, there exists a vector PR3 and a scalar dR such that γγ(t)P=d,t(a,b).

Proof
Without loss of generality we can assume that γγ is unit speed. Indeed, if we were to consider γγ~ a unit speed reparametrization of γγ, then

  • γγ~ would still be contained in the same plane in which γγ is contained.
  • The torsion of γγ~ would not change, i.e., it would still be identically zero.

Thefore the Frenet frame of γγ exists. We denote it by {γγ˙(t),n(t),b(t)}.

Step 1. Suppose that τ=0 for all t. By the Frenet-Serret equations we have b˙=τ(t)n=00, so that b(t) is constant. As by definition b=γγ˙×n, we conclude that the vectors γγ˙(t) and n(t) always span the same plane, which has constant normal vector b. Intuition suggests that γγ should be contained in such plane, see Figure Figure 2.7 below. Indeed, recall that the Frenet frame is orthonormal. Hence γγ˙b=0,t(a,b). Then ddt(γγb)=γγ˙b+γγb˙=0,t(a,b), since b˙=0. Thus γγb is a constant scalar function, meaning that there exists costant dR such that γγ(t)b=d,t(a,b). The says that γγ is contained in a plane.

Step 2. Suppose that γγ is contained in a plane. Hence there exists PR3 and dR such that γγ(t)P=d,t(a,b). We can differentiate the above equation twice to obtain γγ˙P=0,γγ¨P=0, where we used that P and d are constant. By Frenet-Serret we have γγ¨(t)=κ(t)n(t). Therefore the already proven relation γγ¨P=0 implies κ(t)n(t)P=0. As we are assuming κ0, we deduce that n(t)P=0,t(a,b). We have shown that γγ˙(t) and n(t) are both orthogonal to P. Since b(t) is orthogonal to γγ˙(t) and n(t), we conclude that b(t) is parallel to P. Hence, there exists λ(t)R such that (2.20)b(t)=λ(t)Pt(a,b). Since b=1 and P is constant, from (2.20) we conclude that λ(t) is constant. Differentiating (2.20) we obtain b˙(t)=0,t(a,b). By definition of torsion we thus have τ(t)=b˙n(t)=0,t(a,b).

Figure 2.7: If b is constant, then γγ lies in the plane spanned by γγ˙ and n.
Proposition 56

Let γγ:(a,b)R3 be a unit speed curve. They are equivalent:

  1. The image of γγ is contained in a circle of radius 1/c.

  2. The curvature and torsion of γγ satisfy κ(t)=c,τ(t)=0,t(a,b), for some constant cR.

Proposition 56 is actually a consequence of Theorem 50, and of the fact that we have computed that for a circle of radius R one has κ=1R,τ=0. Therefore, by Theorem 50, every unit speed curve γγ with constant curvature and torsion must be equal to a circle, up to rigid motions.

Nevertheless, we still give a proof of Proposition 56, to show yet another application of the Frenet-Serret equations.

Proof
Step 1. Suppose the image of γγ is contained in a circle of radius 1/c. Then, up to a translation, γγ is parametrized by γγ(t)=(1ccos(t),1csin(t),0) for t in some interval (a~,b~). We have already seen that in this case κ=c,τ=0, concluding the proof.

Step 2. Suppose that κ(t)=c,τ(t)=0,t(a,b), for some constant cR. Since γγ is unit speed, its Frenet-Serret equations are: γγ¨=κn=cnn˙=κγγ˙+τb=cγγ˙b˙=τn=0 In particular b˙=0 and so b is a constant vector. As seen in the proof Proposition 55, this implies that γγ is contained in a plane ππ orthogonal to b, see Figure 2.7. As c is constant we get ddt(γγ+1cn)=γγ˙+1cn˙=γγ˙1ccγγ˙=0, where we used the second Frenet-Serret equation. Therefore γγ(t)+1cn(t)=p,t(a,b), for some constant point pR3. In particular γγ(t)p=1cn(t)=1c, since n is a unit vector. The above shows that γγ is contained in a sphere of radius 1/c and center p. In formulas: γγ((a,b))S:={xR3:xp=1/c}. The intersection of S with the plane ππ is a circle C with some radius R. Since γγ((a,b))ππ,γγ((a,b))S, this implies (2.21)γγ((a,b))ππS=C. Thus γγ parametrizes part of C. From Step 1 it follows that the curvature and torsion of γγ must satisfy κ=1R,τ=0. Since we already know that κ=c, we conclude that R=1/c. Therefore the circle C has radius 1/c and the thesis follows by (2.21).