MathLabs

Topology

Smooth manifolds

Spaces built from an atlas {(Uα,φα)}α∈A\{(U_\alpha, \varphi_\alpha)\}_{\alpha \in A} that locally look like Rn\mathbb{R}^n and support calculus, culminating in Stokes' theorem ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M} \omega, the setting for modern geometry, physics, and robotics.

IntuitionCurved spaces that look flat up close

Standing on Earth's surface, the ground looks flat even though the planet is a sphere. A smooth manifold formalizes this: a space that, zoomed in near any point, looks like ordinary Euclidean space Rn\mathbb{R}^n, even though globally it can be curved, twisted, or have interesting topology (like a sphere, torus, or the space of rotations). This "locally flat, globally curved" idea lets us do calculus — derivatives, integrals, vector fields — on shapes that are not flat. The parametric surface below (sphere, torus, helicoid) shows several manifolds whose local patches are each a distorted copy of the flat plane.

Rotating parametric torus surface.
A torus built by gluing local coordinate patches — a smooth 22-manifold embedded in R3\mathbb{R}^3.

UndergraduateAtlases, tangent spaces, and differentials

Definition: Smooth atlas and manifold

An nn-dimensional smooth manifold MM is a topological space covered by charts {(Uα,φα)}α∈A\{(U_\alpha, \varphi_\alpha)\}_{\alpha \in A}, where each φα:Uα→Rn\varphi_\alpha : U_\alpha \to \mathbb{R}^n is a homeomorphism onto an open subset, and on overlaps the transition maps φβ∘φα−1:φα(Uα∩Uβ)→φβ(Uα∩Uβ)\varphi_\beta \circ \varphi_\alpha^{-1} : \varphi_\alpha(U_\alpha \cap U_\beta) \to \varphi_\beta(U_\alpha \cap U_\beta) are C∞C^\infty (infinitely differentiable). This smoothness of transitions is what lets us transfer calculus from Rn\mathbb{R}^n, done in any one chart, consistently to the whole manifold.

φβ∘φα−1:φα(Uα∩Uβ)→φβ(Uα∩Uβ)\varphi_\beta \circ \varphi_\alpha^{-1} : \varphi_\alpha(U_\alpha \cap U_\beta) \to \varphi_\beta(U_\alpha \cap U_\beta)

At each point p∈Mp \in M, the tangent space TpMT_p M is the nn-dimensional vector space of all velocity vectors of curves through pp — it is the best linear approximation to MM at pp. A smooth map f:M→Nf: M \to N then has a differential (derivative) at each point, a linear map between tangent spaces:

dfp:TpM→Tf(p)Ndf_p : T_p M \to T_{f(p)} N
Manifolds and their dimension / tangent space
ManifoldDimensionTpMT_pM description
Sphere Sn−1⊂RnS^{n-1} \subset \mathbb{R}^nn−1n-1Vectors orthogonal to pp
Rotation group SO(3)SO(3)33Skew-symmetric matrices (angular velocity)
Torus T2=S1×S1T^2 = S^1\times S^122Plane spanned by two angle directions

AdvancedTwo central theorems

Let f:Mm→Nnf: M^m \to N^n be smooth and q∈Nq \in N a regular value (meaning dfpdf_p is surjective for every p∈f−1(q)p \in f^{-1}(q)). Then f−1(q)f^{-1}(q) is a smooth submanifold of MM of dimension m−nm-n.

Why is it true?

This theorem is the primary tool for manufacturing manifolds: instead of exhibiting an atlas by hand, we describe a shape as the zero set of a map and just check a linear-algebra condition (surjectivity of the differential) at each solution point.

Proof

Fix p∈f−1(q)p \in f^{-1}(q). Since dfp:TpM→TqNdf_p: T_pM \to T_qN is surjective and dim⁡TpM=m≥n=dim⁡TqN\dim T_pM = m \ge n = \dim T_qN, its kernel K=ker⁡dfpK = \ker df_p has dimension m−nm-n. Choose a linear complement WW so TpM=K⊕WT_pM = K \oplus W with dim⁡W=n\dim W = n; then dfp∣W:W→TqNdf_p|_W : W \to T_qN is an isomorphism.

Work in local coordinates centered at pp and qq (via charts), so ff becomes a smooth map Rm→Rn\mathbb{R}^m \to \mathbb{R}^n with f(0)=0f(0)=0 and df0df_0 surjective. Reorder coordinates (x,y)∈Rm−n×Rn(x,y) \in \mathbb{R}^{m-n}\times\mathbb{R}^n so that ∂f/∂y\partial f/\partial y at 00 is the invertible n×nn\times n block (possible since df0df_0 has rank nn).

Define Φ(x,y)=(x,f(x,y))\Phi(x,y) = (x, f(x,y)). Then dΦ0=(I0∂f/∂x∂f/∂y)d\Phi_0 = \begin{pmatrix} I & 0 \\ \partial f/\partial x & \partial f/\partial y \end{pmatrix} has det⁡dΦ0=det⁡(∂f/∂y)≠0\det d\Phi_0 = \det(\partial f/\partial y) \ne 0, so Φ\Phi is a local diffeomorphism by the Inverse Function Theorem. In the new coordinates (x,y′)=Φ(x,y)(x,y') = \Phi(x,y), the equation f=qf=q (i.e. f=0f=0) becomes exactly y′=0y'=0.

So near pp, f−1(q)f^{-1}(q) is the set {y′=0}\{y'=0\}, which in these coordinates is literally an (m−n)(m-n)-dimensional coordinate slice — a smooth chart for f−1(q)f^{-1}(q). Since pp was arbitrary, every point of f−1(q)f^{-1}(q) has such a chart, and the transition maps between these charts are restrictions of the (smooth) transition maps of MM, hence smooth. Therefore f−1(q)f^{-1}(q) is a smooth (m−n)(m-n)-dimensional submanifold.

For a compact oriented nn-manifold with boundary MM and a smooth (n−1)(n-1)-form ω\omega on MM, ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M} \omega.

Why is it true?

This single identity unifies the fundamental theorem of calculus, Green's theorem, the divergence theorem, and the classical Stokes' theorem from vector calculus into one statement about differential forms, and is the analytic engine behind de Rham cohomology.

Proof

Step 1 (local case, half-space). First suppose M=Hn={xn≥0}M = \mathbb{H}^n = \{x_n \ge 0\} and ω\omega has compact support in a single chart. Write ω=∑ifi dx1∧⋯dxi^⋯∧dxn\omega = \sum_i f_i\, dx_1\wedge\cdots\widehat{dx_i}\cdots\wedge dx_n. Then dω=∑i(−1)i−1∂fi∂xidx1∧⋯∧dxnd\omega = \sum_i (-1)^{i-1}\frac{\partial f_i}{\partial x_i}dx_1\wedge\cdots\wedge dx_n, and ∫Hndω=∑i(−1)i−1∫∂fi∂xi dx1⋯dxn\int_{\mathbb{H}^n} d\omega = \sum_i (-1)^{i-1}\int \frac{\partial f_i}{\partial x_i}\,dx_1\cdots dx_n.

For i<ni < n, integrating ∂fi/∂xi\partial f_i/\partial x_i over xi∈Rx_i \in \mathbb{R} first and using compact support gives 00 by the ordinary Fundamental Theorem of Calculus (fi→0f_i \to 0 at xi=±∞x_i=\pm\infty). For i=ni=n, integrating over xn∈[0,∞)x_n \in [0,\infty) gives ∫∂fn∂xndxn=[fn]0∞=−fn(x1,…,xn−1,0)\int \frac{\partial f_n}{\partial x_n}dx_n = [f_n]_0^\infty = -f_n(x_1,\dots,x_{n-1},0) (again using compact support at xn=∞x_n=\infty), so only the i=ni=n term survives: ∫Hndω=(−1)n−1∫Rn−1(−fn(x1,…,xn−1,0))dx1⋯dxn−1\int_{\mathbb{H}^n}d\omega = (-1)^{n-1}\int_{\mathbb{R}^{n-1}} \left(-f_n(x_1,\dots,x_{n-1},0)\right)dx_1\cdots dx_{n-1}.

On the boundary ∂Hn={xn=0}\partial\mathbb{H}^n = \{x_n=0\} (oriented so the outward normal −∂n-\partial_n comes last, giving orientation sign (−1)n(-1)^n), the restriction of ω\omega is ω∣∂=fn dx1∧⋯∧dxn−1\omega|_{\partial} = f_n\, dx_1\wedge\cdots\wedge dx_{n-1} (all other terms restrict to 00 since they contain dxndx_n or vanish on the slice). A direct sign check using the standard boundary orientation convention shows ∫∂Hnω=(−1)n∫fn dx1⋯dxn−1\int_{\partial \mathbb{H}^n}\omega = (-1)^n \int f_n\,dx_1\cdots dx_{n-1}, matching the formula above exactly. So ∫Hndω=∫∂Hnω\int_{\mathbb{H}^n}d\omega = \int_{\partial\mathbb{H}^n}\omega in this local model.

Step 2 (partition of unity, globalize). For general MM and general ω\omega, cover MM by finitely many charts {(Uα,φα)}\{(U_\alpha,\varphi_\alpha)\} (using compactness) and choose a smooth partition of unity {ρα}\{\rho_\alpha\} subordinate to this cover, i.e. ∑αρα=1\sum_\alpha \rho_\alpha = 1 with supp⁡ρα⊂Uα\operatorname{supp}\rho_\alpha \subset U_\alpha. Write ω=∑αραω\omega = \sum_\alpha \rho_\alpha\omega; each ραω\rho_\alpha\omega has compact support inside a single chart, where the chart is either entirely interior (in which case ∫∂ραω=0\int_{\partial}\rho_\alpha\omega = 0 trivially and ∫Md(ραω)=0\int_M d(\rho_\alpha\omega)=0 by Step 1 applied to Rn\mathbb{R}^n with no boundary) or meets ∂M\partial M (Step 1 applies directly after transporting via φα\varphi_\alpha, which preserves both dd and orientation).

Since dd is linear, dω=∑αd(ραω)d\omega = \sum_\alpha d(\rho_\alpha\omega) (using ∑αdρα=d(∑αρα)=d(1)=0\sum_\alpha d\rho_\alpha = d(\sum_\alpha \rho_\alpha) = d(1) = 0 to handle the cross terms dρα∧ωd\rho_\alpha \wedge \omega correctly when summed). Integrating and summing the local identities from Step 1 over all α\alpha: ∫Mdω=∑α∫Md(ραω)=∑α∫∂Mραω=∫∂Mω\int_M d\omega = \sum_\alpha \int_M d(\rho_\alpha\omega) = \sum_\alpha \int_{\partial M} \rho_\alpha\omega = \int_{\partial M}\omega, which is exactly ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M} \omega.

AdvancedManufacturing manifolds: the sphere and O(n)O(n)

Apply the Regular Value Theorem to f:Rn→Rf:\mathbb{R}^n \to \mathbb{R}, f(x)=∣x∣2f(x)=|x|^2: since dfx(v)=2⟨x,v⟩df_x(v) = 2\langle x,v\rangle is surjective (nonzero) for every x≠0x\ne 0, q=1q=1 is a regular value, so Sn−1=f−1(1)S^{n-1} = f^{-1}(1) is a smooth (n−1)(n-1)-manifold — no explicit atlas needed! Similarly, for f:Rn×n→Sym⁡(n)f:\mathbb{R}^{n\times n}\to \operatorname{Sym}(n), f(A)=ATAf(A)=A^TA (landing in the (n+12)\binom{n+1}{2}-dimensional space of symmetric matrices), one checks dfAdf_A is surjective at every AA with ATA=IA^TA=I, so O(n)=f−1(I)O(n) = f^{-1}(I) is a smooth manifold of dimension n2−(n+12)=(n2)n^2 - \binom{n+1}{2} = \binom{n}{2} — this single computation instantly gives SO(3)SO(3), the configuration space of a rigid body's rotations, dimension 33.

UndergraduateReal-World Applications and Worked Examples

Robotics models a rigid body's orientation as a point on the manifold SO(3)SO(3), not as three Euler angles (which have coordinate singularities — gimbal lock). Motion planning, control, and state estimation (e.g. the Kalman filter used on satellites and drones) are all done using the tangent space so(3)\mathfrak{so}(3) (angular velocities) and the exponential map, exploiting the smooth manifold structure to avoid singularities. In geometric mechanics, a mechanical system's configuration space is a manifold MM (e.g. a double pendulum's is T2T^2), and Lagrangian/Hamiltonian dynamics live naturally on TMTM and T∗MT^*M; conserved quantities correspond to symmetries of MM via Noether's theorem, a manifold-theoretic upgrade of classical mechanics.

Example: Dimension of the special orthogonal group SO(3)SO(3)

Using f(A)=ATAf(A)=A^TA on 3×33\times 3 real matrices, confirm the dimension of O(3)O(3) (hence SO(3)SO(3), its identity component) predicted by the Regular Value Theorem.

Solution

The domain R3×3\mathbb{R}^{3\times 3} has dimension n2=9n^2=9. The codomain Sym⁡(3)\operatorname{Sym}(3) (symmetric 3×33\times 3 matrices) has dimension (n+12)=(42)=6\binom{n+1}{2} = \binom{4}{2}=6.

One checks dfAdf_A is surjective at every A∈O(3)A\in O(3): differentiating f(A)=ATAf(A)=A^TA gives dfA(H)=HTA+ATHdf_A(H) = H^TA + A^TH; for any symmetric SS, take H=12ASH = \tfrac12 AS, then dfA(H)=12STATA+12ATAS=12S+12S=Sdf_A(H) = \tfrac12 S^TA^TA + \tfrac12 A^TAS = \tfrac12 S + \tfrac12 S = S (using ATA=IA^TA=I), showing surjectivity.

By the Regular Value Theorem, O(3)=f−1(I)O(3) = f^{-1}(I) has dimension 9−6=39 - 6 = 3. Since SO(3)SO(3) is the connected component of the identity, it has the same dimension, 33 — matching the intuitive count of 33 independent rotation axes/angles (e.g. roll, pitch, yaw).

Example: Verifying Stokes' theorem on a disk

Let MM be the unit disk {x2+y2≤1}⊂R2\{x^2+y^2\le 1\}\subset\mathbb{R}^2 with boundary the unit circle, and ω=x dy\omega = x\,dy. Verify ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M}\omega directly.

Solution

Compute the left side: dω=dx∧dyd\omega = dx\wedge dy, so ∫Mdω=∫Mdx dy=π(1)2=π\int_M d\omega = \int_M dx\,dy = \pi(1)^2 = \pi (the area of the unit disk).

Compute the right side: parametrize ∂M\partial M by x=cos⁡t,y=sin⁡tx=\cos t, y=\sin t, t∈[0,2π]t\in[0,2\pi], so dy=cos⁡t dtdy = \cos t\,dt. Then ∫∂Mω=∫02πcos⁡t⋅cos⁡t dt=∫02πcos⁡2t dt\int_{\partial M}\omega = \int_0^{2\pi} \cos t \cdot \cos t\,dt = \int_0^{2\pi}\cos^2t\,dt.

Using cos⁡2t=1+cos⁡2t2\cos^2 t = \tfrac{1+\cos 2t}{2}, ∫02πcos⁡2t dt=12[t+sin⁡2t2]02π=12(2π)=π\int_0^{2\pi}\cos^2t\,dt = \tfrac12\left[t+\tfrac{\sin 2t}{2}\right]_0^{2\pi} = \tfrac12(2\pi) = \pi. Both sides equal π\pi, confirming Stokes' theorem in this case.

What smoothness condition must the transition maps φβ∘φα−1:φα(Uα∩Uβ)→φβ(Uα∩Uβ)\varphi_\beta \circ \varphi_\alpha^{-1} : \varphi_\alpha(U_\alpha \cap U_\beta) \to \varphi_\beta(U_\alpha \cap U_\beta) satisfy for an atlas to define a smooth manifold?

By the Regular Value Theorem applied to f(A)=ATAf(A)=A^TA on 3×33\times 3 matrices, what is dim⁡SO(3)\dim SO(3)?

In robotics, why is SO(3)SO(3) preferred over three Euler angles to represent orientation?

Stokes' theorem ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M} \omega generalizes which classical result when MM is a 11-dimensional interval [a,b][a,b]?

References

  1. John M. Lee (2012). Introduction to Smooth Manifolds
  2. Victor Guillemin, Alan Pollack (1974). Differential Topology
  3. F. Bullo, R. M. Murray (1999). Riemannian Manifolds in Robot Motion Planning and Control