Differential Geometry of Curves in 3D

⚡ Key Takeaways

Differential geometry of curves in 3D studies space curves via the calculus of their parametrizations, capturing shape through two scalar invariants: curvature κ\kappa (how sharply the curve bends) and torsion τ\tau (how it twists out of a plane). The central result is the Frenet–Serret apparatus: a moving orthonormal frame {T,N,B}\{\mathbf{T},\mathbf{N},\mathbf{B}\} whose rotation is governed by a tidy system of ODEs. Applied throughout physics (particle trajectories, optics), engineering (roads, robotics), and graphics. Developed by Frenet (1847) and Serret (1851).

💡 The Core Insight

Attach to each point of a curve a tiny rigid frame of three perpendicular unit vectors that ride along with the curve. As you slide along, this frame rotates. The rate and axes of that rotation are precisely curvature and torsion. Geometry becomes the study of how a frame turns — an idea that scales up to the whole of modern differential geometry.

Formal Statement

DefinitionParametrized curve and regularity

A parametrized differentiable curve is a smooth map γ:I→R3\boldsymbol{\gamma}: I \to \mathbb{R}^3, t↦(x(t),y(t),z(t))t \mapsto (x(t), y(t), z(t)). It is regular if its velocity never vanishes:

1
γ′(t)≠0for all t∈I.\quad \boldsymbol{\gamma}'(t) \neq \mathbf{0} \quad \text{for all } t \in I.

Derivation hint: regularity guarantees a well-defined tangent direction γ′/∥γ′∥\boldsymbol{\gamma}'/\|\boldsymbol{\gamma}'\| at every point; without it cusps and stops are allowed. The arc length measured from t0t_0 is

2
s(t)=∫t0t∥γ′(u)∥ du,dsdt=∥γ′(t)∥=:v(t).\quad s(t) = \int_{t_0}^{t} \|\boldsymbol{\gamma}'(u)\|\, du, \qquad \frac{ds}{dt} = \|\boldsymbol{\gamma}'(t)\| =: v(t).

Derivation hint: the length of an infinitesimal chord is ∥γ′(u)∥ du\|\boldsymbol{\gamma}'(u)\|\,du by the Pythagorean theorem applied to (dx,dy,dz)(dx,dy,dz); integrate. A curve is unit speed (parametrized by arc length) if ∥γ′(s)∥=1\|\boldsymbol{\gamma}'(s)\| = 1 for all ss. Symbols: tt a general parameter, ss arc length, vv speed, prime ′' denotes d/dtd/dt or d/dsd/ds as indicated.

DefinitionThe Frenet frame, curvature, torsion

Let γ(s)\boldsymbol{\gamma}(s) be unit speed with γ′′(s)≠0\boldsymbol{\gamma}''(s) \neq \mathbf{0}. Define:

3
T(s)=γ′(s),κ(s)=∥T′(s)∥=∥γ′′(s)∥,N(s)=T′(s)κ(s),B(s)=T(s)×N(s).\quad \mathbf{T}(s) = \boldsymbol{\gamma}'(s), \qquad \kappa(s) = \|\mathbf{T}'(s)\| = \|\boldsymbol{\gamma}''(s)\|, \qquad \mathbf{N}(s) = \frac{\mathbf{T}'(s)}{\kappa(s)}, \qquad \mathbf{B}(s) = \mathbf{T}(s) \times \mathbf{N}(s).

Derivation hint: differentiating T⋅T=1\mathbf{T}\cdot\mathbf{T}=1 gives T⋅T′=0\mathbf{T}\cdot\mathbf{T}'=0, so T′⊥T\mathbf{T}'\perp\mathbf{T}; its length is κ\kappa and its unit direction is N\mathbf{N}; B\mathbf{B} completes a right-handed frame. T\mathbf{T} is the unit tangent, N\mathbf{N} the principal normal, B\mathbf{B} the binormal; {T,N,B}\{\mathbf{T},\mathbf{N},\mathbf{B}\} is a positively-oriented orthonormal basis. Torsion τ(s)\tau(s) is defined by

4
B′(s)=−τ(s) N(s).\quad \mathbf{B}'(s) = -\tau(s)\,\mathbf{N}(s).

Derivation hint: B′⊥B\mathbf{B}'\perp\mathbf{B} (unit length) and B′⊥T\mathbf{B}'\perp\mathbf{T} (from B⋅T=0\mathbf{B}\cdot\mathbf{T}=0 differentiated, using T′=κN\mathbf{T}'=\kappa\mathbf{N}); hence B′\mathbf{B}' is parallel to N\mathbf{N}, with the sign convention giving torsion.

TheoremFrenet–Serret formulas

For a unit-speed curve with κ>0\kappa > 0, the moving frame satisfies

5
(T′N′B′)=(0κ0−κ0τ0−τ0)(TNB).\quad \begin{pmatrix} \mathbf{T}' \\ \mathbf{N}' \\ \mathbf{B}' \end{pmatrix} = \begin{pmatrix} 0 & \kappa & 0 \\ -\kappa & 0 & \tau \\ 0 & -\tau & 0 \end{pmatrix} \begin{pmatrix} \mathbf{T} \\ \mathbf{N} \\ \mathbf{B} \end{pmatrix}.

Derivation hint: T′=κN\mathbf{T}'=\kappa\mathbf{N} and B′=−τN\mathbf{B}'=-\tau\mathbf{N} are definitions; for N′\mathbf{N}' write N′=aT+bN+cB\mathbf{N}'=a\mathbf{T}+b\mathbf{N}+c\mathbf{B} and use orthonormality to read off a=−κa=-\kappa, b=0b=0, c=τc=\tau. The coefficient matrix is skew-symmetric, reflecting that the frame's evolution is an instantaneous rotation. The rotation axis is the Darboux vector ω=τ T+κ B\boldsymbol{\omega} = \tau\,\mathbf{T} + \kappa\,\mathbf{B}, with T′=ω×T\mathbf{T}'=\boldsymbol{\omega}\times\mathbf{T}, etc.

TheoremFundamental Theorem of Space Curves

Given smooth functions κ(s)>0\kappa(s) > 0 and τ(s)\tau(s) on an interval II, there exists a unit-speed curve γ(s)\boldsymbol{\gamma}(s) with exactly these curvature and torsion. Moreover the curve is unique up to a rigid motion (rotation + translation) of R3\mathbb{R}^3.

ProofExistence and uniqueness

Step 1. Treat the Frenet system [5] as a linear ODE for the nine components of (T,N,B)(\mathbf{T},\mathbf{N},\mathbf{B}). By the Picard–Lindelöf theorem (smooth, in fact linear, coefficients), a unique solution exists given an initial orthonormal frame at s0s_0.

Step 2 (orthonormality preserved). Let G(s)G(s) be the 3×33\times 3 matrix of inner products of the frame vectors. Differentiating and using skew-symmetry of the coefficient matrix AA gives G′=AG+GA ⁣⊤G' = AG + GA^{\!\top} ; since A ⁣⊤=−AA^{\!\top}=-A, G≡IG\equiv I is a solution, and uniqueness forces G(s)=IG(s)=I for all ss. So the frame stays orthonormal.

Step 3 (recover the curve). Define γ(s)=γ(s0)+∫s0sT(u) du\boldsymbol{\gamma}(s) = \boldsymbol{\gamma}(s_0) + \int_{s_0}^{s}\mathbf{T}(u)\,du. Then γ′=T\boldsymbol{\gamma}' = \mathbf{T} is unit length, γ′′=T′=κN\boldsymbol{\gamma}''=\mathbf{T}'=\kappa\mathbf{N} gives curvature κ\kappa, and the Frenet equations reproduce τ\tau.

Step 4 (uniqueness). If two curves share κ,τ\kappa,\tau, apply the rigid motion taking the first frame at s0s_0 onto the second. Both then satisfy the same ODE with the same initial data; by uniqueness their frames coincide for all ss, hence (Step 3) the curves coincide. ■\blacksquare

∎
Figure
The Frenet frame riding along a helix, showing T, N, B at one point
γ(s) T N B Orthonormal frame: T tangent, N points to bend, B = T×N

Why It Holds / Derivation

ProofCurvature and torsion for a general (non-unit-speed) parametrization

Often a curve comes with an arbitrary parameter tt. Using the chain rule dds=1vddt\frac{d}{ds} = \frac{1}{v}\frac{d}{dt} and r′=γ′(t)\mathbf{r}'=\boldsymbol{\gamma}'(t):

Step 1. T=r′/v\mathbf{T} = \mathbf{r}'/v, and r′×r′′=v3 κ B\mathbf{r}' \times \mathbf{r}'' = v^3\,\kappa\,\mathbf{B} after expanding r′′=v′T+v2κN\mathbf{r}'' = v'\mathbf{T} + v^2\kappa\mathbf{N}.

Step 2. Taking magnitudes gives the curvature formula

6
κ=∥r′×r′′∥∥r′∥3.\quad \kappa = \frac{\|\mathbf{r}' \times \mathbf{r}''\|}{\|\mathbf{r}'\|^3}.

Derivation hint: r′=vT\mathbf{r}'=v\mathbf{T} and r′′=v′T+v2κN\mathbf{r}''=v'\mathbf{T}+v^2\kappa\mathbf{N}; cross them, the T×T\mathbf{T}\times\mathbf{T} term dies, leaving v3κ(T×N)=v3κBv^3\kappa(\mathbf{T}\times\mathbf{N})=v^3\kappa\mathbf{B}.

Step 3. Differentiate once more and project onto B\mathbf{B} to isolate torsion:

7
τ=(r′×r′′)⋅r′′′∥r′×r′′∥2.\quad \tau = \frac{(\mathbf{r}' \times \mathbf{r}'') \cdot \mathbf{r}'''}{\|\mathbf{r}' \times \mathbf{r}''\|^2}.

Derivation hint: r′′′\mathbf{r}''' has a component v3κτ Bv^3\kappa\tau\,\mathbf{B}; dotting with r′×r′′=v3κB\mathbf{r}'\times\mathbf{r}''=v^3\kappa\mathbf{B} extracts τ\tau after dividing by ∥r′×r′′∥2=v6κ2\|\mathbf{r}'\times\mathbf{r}''\|^2 = v^6\kappa^2. ■\blacksquare

∎
mathematician · Jean Frédéric Frenet | 1816–1900

French mathematician and astronomer who, in his 1847 doctoral thesis at Toulouse, derived six of the formulas now bearing his name, expressing the derivatives of the moving frame. He later directed the Lyon observatory. Working independently and slightly later (1851), Joseph Alfred Serret published the complete system, which is why the equations carry both names. Their compact "frame rotation" viewpoint became a template for Cartan's method of moving frames decades later.

Building Intuition

🧠 Intuition

Imagine driving a car along the curve at unit speed. The direction you face is T\mathbf{T}. Curvature is how hard you turn the steering wheel: a sharp bend means large κ\kappa, a straight road means κ=0\kappa=0. The vector N\mathbf{N} points toward the inside of the turn — the direction the car accelerates sideways. Now imagine the car's plane (the windshield's plane) can also bank, like an airplane rolling: torsion measures how fast that plane tilts out of the current bending plane. A flat racetrack has τ=0\tau=0 (a planar curve); a corkscrew on-ramp has τ≠0\tau \neq 0. The osculating plane spanned by T,N\mathbf{T},\mathbf{N} is the plane that best hugs the curve locally; torsion is exactly the rate at which the curve lifts off it.

Worked Examples

ExampleCircular helix

Given: γ(t)=(acos⁡t, asin⁡t, bt)\boldsymbol{\gamma}(t) = (a\cos t,\ a\sin t,\ b t) with a=3a = 3, b=4b = 4 (units: metres), t∈Rt\in\mathbb{R}.

Step 1 — derivatives. r′=(−3sin⁡t,3cos⁡t,4)\mathbf{r}' = (-3\sin t, 3\cos t, 4), so v=9+16=5v = \sqrt{9+16}=5. r′′=(−3cos⁡t,−3sin⁡t,0)\mathbf{r}'' = (-3\cos t, -3\sin t, 0), r′′′=(3sin⁡t,−3cos⁡t,0)\mathbf{r}''' = (3\sin t, -3\cos t, 0).

Step 2 — curvature (eq. 6). r′×r′′=(12sin⁡t,−12cos⁡t,9)\mathbf{r}'\times\mathbf{r}'' = (12\sin t, -12\cos t, 9), magnitude 144+81=15\sqrt{144+81}=15. κ=15/53=15/125=0.12 m−1.\kappa = 15 / 5^3 = 15/125 = 0.12\ \text{m}^{-1}. (General formula: κ=a/(a2+b2)=3/25\kappa = a/(a^2+b^2) = 3/25.)

Step 3 — torsion (eq. 7). (r′×r′′)⋅r′′′=(12sin⁡t)(3sin⁡t)+(−12cos⁡t)(−3cos⁡t)+0=36(\mathbf{r}'\times\mathbf{r}'')\cdot\mathbf{r}''' = (12\sin t)(3\sin t)+(-12\cos t)(-3\cos t)+0 = 36. τ=36/152=36/225=0.16 m−1.\tau = 36 / 15^2 = 36/225 = 0.16\ \text{m}^{-1}. (General: τ=b/(a2+b2)=4/25\tau = b/(a^2+b^2)=4/25.)

κ=0.12 m−1,τ=0.16 m−1 (both constant)\boxed{\kappa = 0.12\ \text{m}^{-1},\qquad \tau = 0.16\ \text{m}^{-1}\ \text{(both constant)}}

A helix is the only curve with both κ\kappa and τ\tau constant and nonzero.

ExamplePlane parabola — torsion vanishes

Given: γ(t)=(t, t2, 0)\boldsymbol{\gamma}(t)=(t,\ t^2,\ 0), evaluate at t=1t=1 (dimensionless).

Step 1. r′=(1,2t,0)\mathbf{r}'=(1,2t,0), r′′=(0,2,0)\mathbf{r}''=(0,2,0), r′′′=(0,0,0)\mathbf{r}'''=(0,0,0).

Step 2 — curvature. r′×r′′=(0,0,2)\mathbf{r}'\times\mathbf{r}''=(0,0,2), magnitude 22. At t=1t=1, ∥r′∥=5\|\mathbf{r}'\|=\sqrt{5}, so κ=2/(5)3=2/(55)≈0.179.\kappa = 2/(\sqrt5)^3 = 2/(5\sqrt5) \approx 0.179.

Step 3 — torsion. Numerator (r′×r′′)⋅r′′′=0(\mathbf{r}'\times\mathbf{r}'')\cdot\mathbf{r}''' = 0 since r′′′=0\mathbf{r}'''=\mathbf{0}. τ=0.\tau = 0.

κ≈0.179,τ=0 (curve lies in the xy-plane)\boxed{\kappa \approx 0.179,\qquad \tau = 0\ \text{(curve lies in the }xy\text{-plane)}}

Zero torsion everywhere   ⟺  \iff the curve is planar.

Common Misconceptions

⚠️ Common mistake

Curvature and torsion are NOT parametrization-dependent. A common error is computing ∥γ′′(t)∥\|\boldsymbol{\gamma}''(t)\| and calling it the curvature — this only works for unit-speed curves. For general tt you must use eq. [6]. Likewise, do not confuse the second derivative's magnitude with κ\kappa unless v=1v=1.

Torsion can be negative. Unlike κ≥0\kappa \ge 0, torsion is a signed quantity. Its sign distinguishes a right-handed from a left-handed helix and is physically meaningful (chirality of DNA, screw threads).

The Frenet frame fails where κ=0\kappa = 0. If the curve momentarily straightens, N\mathbf{N} is undefined (division by zero in eq. [3]). Use the Bishop (parallel-transport) frame to avoid this discontinuity.

Counterexample

Consider a curve that is a straight segment, then bends, then a straight segment again, joined smoothly. On each straight piece κ=0\kappa=0, so N\mathbf{N} and B\mathbf{B} are undefined, and torsion is meaningless there. A naive numerical Frenet frame will flip N\mathbf{N} by 180° as κ\kappa passes through zero, producing spurious infinite torsion spikes — a real bug in CAD/animation pipelines. This shows the Frenet frame is a local construction requiring κ>0\kappa>0, not a global one.

Connections

connection

The curve theory here is the one-dimensional ancestor of Differential Geometry of Surfaces, where the single curvature κ\kappa branches into principal, Gaussian, and mean curvatures. The skew-symmetric Frenet matrix is a baby case of a Connection (Differential Geometry) / Cartan's moving frames. Arc length and the variational principle of geodesics connect to Calculus of Variations: straight lines (and geodesics) extremize length. The Darboux vector links to the Angular Velocity of a rigid body, since the frame rotates like a spinning top. In complex analysis and physics, curvature of integral curves underlies Vector Fields and Flows.

Real-World Applications

🔬 Where this appears

Historical Development

⏳ History
- c. 1731 — Alexis Clairaut studies space curves ("courbes à double courbure"), introducing two curvatures.
- 1771 — Leonhard Euler and others develop curvature of plane curves and the radius of curvature.
- 1775 — Gaspard Monge applies analysis to curves and surfaces, founding descriptive/differential geometry.
- 1806 — A.-J. Lancret recognizes that a space curve has two independent invariants (the "first and second curvature").
- 1847 — Jean Frédéric Frenet derives the frame-derivative formulas in his thesis.
- 1851 — Joseph Alfred Serret independently publishes the full system; the equations become "Frenet–Serret."
- 1853 — Bertrand studies curves sharing principal normals (Bertrand mates), enriching the theory.
- 1894 — Gaston Darboux's Leçons systematizes the moving-frame viewpoint (Darboux vector).
- 1900s–1930s — Élie Cartan generalizes to the method of moving frames, the foundation of modern differential geometry.
- 1975 — Richard Bishop introduces the relatively parallel (Bishop) frame, defined even where κ=0\kappa=0.

Subtleties & Limits

open-problem

Several subtleties bound the elementary theory. (1) Frame at inflection points: the Frenet frame is undefined when κ=0\kappa=0; the global existence of a smooth alternative frame relates to topological obstructions and the self-linking / writhe of the curve. (2) Total curvature inequalities: the Fenchel theorem (∫κ ds≥2π\int\kappa\,ds \ge 2\pi for closed curves) and the Fáry–Milnor theorem (∫κ ds>4π\int\kappa\,ds > 4\pi for knotted curves) connect local invariants to global knot topology — quantitative refinements (optimal ropelength, ideal knot shapes) remain active research. (3) Inverse/reconstruction stability: recovering a curve from noisy measured κ(s),τ(s)\kappa(s),\tau(s) is ill-conditioned; robust algorithms are an applied open problem. (4) Generalization to higher dimensions (Frenet frames in Rn\mathbb{R}^n) requires n−1n-1 "generalized curvatures" and degenerates when intermediate derivatives become linearly dependent.

See Also

Curvature · Torsion of a Curve · Frenet–Serret Formulas · Osculating Circle · Bishop Frame · Differential Geometry of Surfaces · Moving Frames · Darboux Vector