Skip to content
The Lumenary
First-year graduateGraduate 150 min

Smooth manifolds

Spaces that look like Rn\R^n up close but not from far away — and how to do calculus on them without ever choosing global coordinates.

Formally verified · 15 LeanComputationally tested · 2 PythonAI-drafted · awaiting human reviewMathematically reviewed: Awaiting human review
MotivationH

Calculus without global coordinates

The surface of the Earth is a sphere, and no single map of it is faithful everywhere: Mercator's projection sends the poles to infinity, and every flat map tears the sphere somewhere. Yet navigators have always worked happily with local maps that overlap, and so does every GPS receiver today. Wherever two maps overlap, there is a rule converting one set of coordinates into the other.

The same situation appears all over mathematics and physics. The positions of a double pendulum form a torus S1×S1S^1 \times S^1 (two angles); the orientations of a rigid body — a satellite, a drone, a molecule — form the group SO(3)SO(3) of rotations, which is not a subset of any R3\R^3 in any natural way; the solutions of x2+y2−z2=1x^2 + y^2 - z^2 = 1 form a hyperboloid; and in general relativity spacetime is a four-dimensional manifold with no preferred coordinates at all. In each case we want to differentiate (velocities, forces, gradients) and integrate (work, flux, volume), but there are no global coordinates in which to do it.

A smooth manifold is the minimal structure that makes this possible: a space covered by overlapping coordinate charts whose changes of coordinates are smooth. Everything in this lesson — tangent vectors, vector fields, differential forms, Stokes's theorem and de Rham cohomology — is calculus made independent of the chart. The payoff is two-fold: calculus on curved and abstract spaces, and a bridge from calculus to topology, since de Rham cohomology uses derivatives to count holes.

This lesson builds on point-set topology, multivariable calculus (the inverse function theorem in particular), dual spaces and multilinear algebra, and it upgrades the calculus of differential forms on open subsets of ℝⁿ to manifolds.

IntuitionH

An atlas of the Earth

An atlas of the Earth is a book of maps. Each page shows a region; neighbouring pages overlap; and in the overlap a point has two pairs of coordinates, related by a conversion formula. Three observations drive the whole theory.

  1. Local questions are chart questions. Whether a curve on the Earth is smooth, how fast a ship is moving, whether a function has a maximum — these can be decided on any single page that contains the point.
  2. The answer must not depend on the page. This is guaranteed if the conversion formulas are smooth: a curve that is smooth on one page is then smooth on the other, by the chain rule. That is the entire content of the definition of a smooth atlas.
  3. Global questions are not chart questions. Whether every closed curve can be shrunk to a point, or whether the wind can blow everywhere at once without a calm spot (the hairy ball theorem), depends on how the pages fit together.

Stereographic projection gives an atlas of the sphere with just two pages: project from the North Pole onto the equatorial plane, and from the South Pole.

NSxu
Stereographic projection from the north pole $N$: the point $x$ of the sphere is sent to the point $u$ where the line $Nx$ meets the equatorial plane. Every point except $N$ is covered; a second chart, projecting from the south pole, covers everything except $S$.

On the overlap (everything except the poles) the conversion from north-chart coordinates uu to south-chart coordinates is u↦u/∣u∣2u \mapsto u/|u|^2, an inversion in the unit circle. It is smooth away from u=0u = 0, so the two charts form a smooth atlas. We will prove this, and Lean will check it.

DefinitionsH

Charts, atlases and smooth structures

Definition— topological manifold

A topological nn-manifold is a topological space MM that is Hausdorff, second countable (its topology has a countable base) and locally Euclidean of dimension nn: every point has an open neighbourhood homeomorphic to an open subset of Rn\R^n.

Definition— chart, atlas, smooth structure

A chart on MM is a pair (U,φ)(U, \varphi) with U⊆MU \subseteq M open and φ:U→φ(U)\varphi : U \to \varphi(U) a homeomorphism onto an open subset of Rn\R^n; the components xi=ri∘φx^i = r^i \circ \varphi are local coordinates. Two charts (U,φ)(U, \varphi) and (V,ψ)(V, \psi) are smoothly compatible if the transition map

ψ∘φ−1:φ(U∩V)→ψ(U∩V)\psi \circ \varphi^{-1} : \varphi(U \cap V) \to \psi(U \cap V)

is C∞C^\infty (its inverse φ∘ψ−1\varphi \circ \psi^{-1} is then C∞C^\infty too, by symmetry, so it is a diffeomorphism). A smooth atlas is a collection of pairwise compatible charts whose domains cover MM. Every smooth atlas is contained in a unique maximal one (add every chart compatible with all charts of the atlas). A smooth structure is a maximal smooth atlas, and a smooth manifold is a topological manifold together with a smooth structure.

Definition— smooth map, diffeomorphism

A map F:M→NF : M \to N between smooth manifolds is smooth if for every p∈Mp \in M there are charts (U,φ)(U, \varphi) at pp and (V,ψ)(V, \psi) at F(p)F(p) with F(U)⊆VF(U) \subseteq V such that the local representative ψ∘F∘φ−1\psi \circ F \circ \varphi^{-1} is C∞C^\infty. By the chain rule and the smoothness of transition maps this does not depend on the charts chosen. A diffeomorphism is a smooth bijection with smooth inverse. Write C∞(M)C^\infty(M) for the algebra of smooth functions M→RM \to \R.

Definition— manifold with boundary

Let Hn={x∈Rn:xn≥0}\mathbb{H}^n = \{x \in \R^n : x^n \ge 0\} be the closed upper half-space. A smooth manifold with boundary is defined in the same way, but with charts onto open subsets of Hn\mathbb{H}^n; a map on a subset of Hn\mathbb{H}^n is smooth if it extends locally to a smooth map on an open subset of Rn\R^n. The points sent to ∂Hn={xn=0}\partial \mathbb{H}^n = \{x^n = 0\} by some (equivalently, every) chart form the boundary ∂M\partial M, a smooth (n−1)(n-1)-manifold without boundary. The closed disc D‾n\overline{D}{}^n is a manifold with boundary ∂D‾n=Sn−1\partial \overline{D}{}^n = S^{n-1}.

Why Hausdorff and second countable? Hausdorffness makes limits and integral curves unique; second countability guarantees partitions of unity — smooth functions ψα≥0\psi_\alpha \ge 0 with locally finite supports, each inside a chart, summing to 11. Partitions of unity are how local constructions (integrals, Riemannian metrics, embeddings) are glued into global ones, and they are why both axioms appear in the definition.

Worked examplesH

Examples, non-examples and boundary cases

Example— Euclidean space, open subsets and graphs

Rn\R^n with the single chart id\mathrm{id} is a smooth manifold. Any open subset of a smooth manifold is one, by restricting charts. So GL(n,R)={A:det⁡A≠0}GL(n, \R) = \{A : \det A \ne 0\}, an open subset of the space Mn(R)≅Rn2M_n(\R) \cong \R^{n^2} of matrices, is a smooth manifold of dimension n2n^2. The graph of a smooth function f:U→Rkf : U \to \R^k (U⊆RnU \subseteq \R^n open) is a smooth nn-manifold with the single chart (x,f(x))↦x(x, f(x)) \mapsto x.

Example— The sphere, with two stereographic charts

On Sn={x∈Rn+1:∣x∣=1}S^n = \{x \in \R^{n+1} : |x| = 1\} write points as (x,t)(x, t) with x∈Rnx \in \R^n, t∈Rt \in \R. Stereographic projection from the north pole N=(0,1)N = (0, 1) is σN(x,t)=x/(1−t)\sigma_N(x, t) = x/(1 - t), with inverse

σN−1(u)=(2u∣u∣2+1,  ∣u∣2−1∣u∣2+1),\sigma_N^{-1}(u) = \left( \frac{2u}{|u|^2 + 1}, \; \frac{|u|^2 - 1}{|u|^2 + 1} \right),

and from the south pole it is σS(x,t)=x/(1+t)\sigma_S(x, t) = x/(1 + t). With r=∣u∣2r = |u|^2,

σS∘σN−1(u)=2u/(r+1)1+(r−1)/(r+1)=2u/(r+1)2r/(r+1)=u∣u∣2(u≠0).\sigma_S \circ \sigma_N^{-1}(u) = \frac{2u/(r+1)}{1 + (r-1)/(r+1)} = \frac{2u/(r+1)}{2r/(r+1)} = \frac{u}{|u|^2} \qquad (u \ne 0).

The inversion u↦u/∣u∣2u \mapsto u/|u|^2 is smooth on Rn∖{0}\R^n \setminus \{0\} and is its own inverse, so {σN,σS}\{\sigma_N, \sigma_S\} is a smooth atlas of SnS^n. Its Jacobian determinant is −∣u∣−2n-|u|^{-2n}, which is negative: the two charts induce opposite orientations, a fact we return to below.

Example— Real projective space

RPn\R P^n is the set of lines through 00 in Rn+1\R^{n+1}, i.e. (Rn+1∖{0})/R×(\R^{n+1} \setminus \{0\})/\R^\times, with the quotient topology (it is Hausdorff and second countable). Write [x0:⋯:xn][x^0 : \dots : x^n] for the line through xx. For each ii, the set Ui={xi≠0}U_i = \{x^i \ne 0\} is open and

φi[x0:⋯:xn]=(x0xi,…,xixi^,…,xnxi)∈Rn\varphi_i [x^0 : \dots : x^n] = \left( \tfrac{x^0}{x^i}, \dots, \widehat{\tfrac{x^i}{x^i}}, \dots, \tfrac{x^n}{x^i} \right) \in \R^n

is a chart (the hat means "omit"). The transition maps divide by one coordinate and are smooth where that coordinate is non-zero. For n=1n = 1: φ0[x0:x1]=x1/x0=t\varphi_0[x^0 : x^1] = x^1/x^0 = t and φ1∘φ0−1(t)=1/t\varphi_1 \circ \varphi_0^{-1}(t) = 1/t on R∖{0}\R \setminus \{0\}, and RP1\R P^1 is diffeomorphic to a circle.

Example— Products, tori and configuration spaces

If MM and NN are smooth manifolds, products of charts form an atlas of M×NM \times N, of dimension dim⁡M+dim⁡N\dim M + \dim N. The torus Tn=S1×⋯×S1T^n = S^1 \times \dots \times S^1 is the configuration space of nn independent rotors; the double pendulum lives on T2T^2, and its phase space (positions and momenta) is the 44-manifold T∗T2≅T2×R2T^*T^2 \cong T^2 \times \R^2.

Example— Non-examples

  • A cross. The union of the two coordinate axes in R2\R^2 is not a topological manifold: a small connected neighbourhood of the origin falls into four pieces when the origin is removed, while a connected open subset of R1\R^1 falls into two and one of Rn\R^n (n≥2n \ge 2) stays connected. The double cone x2+y2=z2x^2 + y^2 = z^2 fails at its apex for the same reason (two pieces).
  • The line with two origins. Glue two copies of R\R along R∖{0}\R \setminus \{0\}. Every point has a neighbourhood homeomorphic to R\R, and the space is second countable, but the two origins cannot be separated by open sets: it is not Hausdorff, and the sequence 1/k1/k converges to both origins.
  • The cone z=x2+y2z = \sqrt{x^2 + y^2} (one nappe) is a topological manifold, homeomorphic to R2\R^2 by vertical projection, and it can be given a smooth structure by declaring that projection a chart. But it is not a smooth embedded submanifold of R3\R^3: no plane is tangent at the apex. A manifold is not "a nice subset of RN\R^N"; smoothness is extra structure.

Example— Different smooth structures, same manifold

On R\R, the single chart ψ(x)=x3\psi(x) = x^3 is a smooth atlas. It is not compatible with the identity chart, since the transition map id∘ψ−1(y)=y1/3\mathrm{id} \circ \psi^{-1}(y) = y^{1/3} is not smooth at 00; so it defines a different maximal atlas. Nevertheless F(x)=x3F(x) = x^3 is a diffeomorphism from (R,ψ)(\R, \psi) to (R,id)(\R, \mathrm{id}): its local representative id∘F∘ψ−1(y)=y\mathrm{id} \circ F \circ \psi^{-1}(y) = y is smooth with smooth inverse. "Different smooth structures" and "non-diffeomorphic smooth manifolds" are different statements. Genuinely exotic structures exist — Milnor's exotic 77-spheres (1956), and uncountably many smooth structures on R4\R^4 — but every topological manifold of dimension at most 33 has a smooth structure, unique up to diffeomorphism (Radó, Moise). See differential topology.

DefinitionsH

Tangent vectors, differentials and vector fields

There is no ambient space in which to draw arrows, so we define a tangent vector by what it does: it differentiates functions.

Definition— tangent space

A derivation at p∈Mp \in M is a linear map v:C∞(M)→Rv : C^\infty(M) \to \R satisfying the Leibniz rule v(fg)=f(p) v(g)+g(p) v(f)v(fg) = f(p)\, v(g) + g(p)\, v(f). The derivations at pp form a vector space, the tangent space TpMT_pM.

Two consequences of the Leibniz rule: v(1)=v(1⋅1)=2v(1)v(1) = v(1 \cdot 1) = 2v(1), so vv kills constants; and vv is local — if f=gf = g near pp then v(f)=v(g)v(f) = v(g) (take a bump function χ\chi with χ(p)=1\chi(p) = 1 supported where f=gf = g; then χ(f−g)=0\chi (f - g) = 0, so 0=v(χ(f−g))=χ(p)v(f−g)+(f−g)(p)v(χ)=v(f−g)0 = v(\chi(f-g)) = \chi(p) v(f - g) + (f-g)(p) v(\chi) = v(f - g)). By locality, TpMT_pM depends only on a neighbourhood of pp, and a chart (U,φ)(U, \varphi) with coordinates xix^i gives derivations

∂∂xi∣pf=∂(f∘φ−1)∂ri(φ(p)).\left.\frac{\partial}{\partial x^i}\right|_p f = \frac{\partial (f \circ \varphi^{-1})}{\partial r^i}\big(\varphi(p)\big).

Theorem— Theorem 1 (dimension of the tangent space)

If dim⁡M=n\dim M = n, the derivations ∂/∂x1∣p,…,∂/∂xn∣p\partial/\partial x^1|_p, \dots, \partial/\partial x^n|_p form a basis of TpMT_pM. Every v∈TpMv \in T_pM equals ∑iv(xi) ∂/∂xi∣p\sum_i v(x^i)\, \partial/\partial x^i|_p.

Proof

By locality we may work in the chart and assume p=a∈Rnp = a \in \R^n, with ff smooth on a ball BB around aa. Hadamard's lemma: for x∈Bx \in B,

f(x)−f(a)=∫01ddtf(a+t(x−a)) dt=∑i=1n(xi−ai) gi(x),gi(x)=∫01∂if(a+t(x−a)) dt,f(x) - f(a) = \int_0^1 \frac{d}{dt} f\big(a + t(x - a)\big)\, dt = \sum_{i=1}^n (x^i - a^i)\, g_i(x), \qquad g_i(x) = \int_0^1 \partial_i f\big(a + t(x-a)\big)\, dt,

and the gig_i are smooth with gi(a)=∂if(a)g_i(a) = \partial_i f(a). Apply vv, using that it kills constants and the Leibniz rule:

v(f)=∑i((ai−ai) v(gi)+gi(a) v(xi−ai))=∑i∂if(a) v(xi).v(f) = \sum_i \Big( (a^i - a^i)\, v(g_i) + g_i(a)\, v(x^i - a^i) \Big) = \sum_i \partial_i f(a)\, v(x^i).

So v=∑iv(xi) ∂i∣av = \sum_i v(x^i)\, \partial_i|_a: the ∂i∣a\partial_i|_a span. They are independent because ∂i∣a(xj)=δij\partial_i|_a (x^j) = \delta_i^j. ∎

∎

Definition— differential

The differential of a smooth map F:M→NF : M \to N at pp is the linear map dFp:TpM→TF(p)NdF_p : T_pM \to T_{F(p)}N,   dFp(v)(f)=v(f∘F)\; dF_p(v)(f) = v(f \circ F).

The chain rule is now a one-line computation: d(G∘F)p(v)(f)=v(f∘G∘F)=dFp(v)(f∘G)=dGF(p)(dFp(v))(f)d(G \circ F)_p(v)(f) = v(f \circ G \circ F) = dF_p(v)(f \circ G) = dG_{F(p)}\big(dF_p(v)\big)(f), so d(G∘F)p=dGF(p)∘dFpd(G \circ F)_p = dG_{F(p)} \circ dF_p, and d(id)p=idd(\mathrm{id})_p = \mathrm{id}. In coordinates xx on MM and yy on NN, dFp(∂/∂xi)=∑j(∂F^j/∂xi) ∂/∂yjdF_p(\partial/\partial x^i) = \sum_j (\partial \hat F^j / \partial x^i)\, \partial/\partial y^j: the matrix of dFpdF_p is the Jacobian of the local representative F^\hat F.

Vectors under a change of chart. If xx and x~\tilde x are two coordinate systems at pp, the chain rule gives ∂/∂xi=∑j(∂x~j/∂xi) ∂/∂x~j\partial/\partial x^i = \sum_j (\partial \tilde x^j/\partial x^i)\, \partial/\partial \tilde x^j. So if v=∑ivi∂/∂xi=∑jv~j∂/∂x~jv = \sum_i v^i \partial/\partial x^i = \sum_j \tilde v^j \partial/\partial \tilde x^j, then

v~j=∑i∂x~j∂xivi,i.e.v~=Jv,J=∂x~∂x.\tilde v^j = \sum_i \frac{\partial \tilde x^j}{\partial x^i} v^i, \qquad \text{i.e.} \quad \tilde v = J v, \quad J = \frac{\partial \tilde x}{\partial x}.

The tangent bundle TM=⨆pTpMTM = \bigsqcup_p T_pM is a smooth 2n2n-manifold, with charts (p,v)↦(x1(p),…,xn(p),v1,…,vn)(p, v) \mapsto (x^1(p), \dots, x^n(p), v^1, \dots, v^n); the transition maps are (x,v)↦(x~(x),J(x)v)(x, v) \mapsto (\tilde x(x), J(x) v), which are smooth. A vector field is a smooth section X:M→TMX : M \to TM. Equivalently, it is a derivation of the algebra C∞(M)C^\infty(M), acting by (Xf)(p)=Xpf(Xf)(p) = X_p f. The Lie bracket [X,Y]f=X(Yf)−Y(Xf)[X, Y] f = X(Yf) - Y(Xf) is again a vector field (the second derivatives cancel), with [X,Y]=∑j(X(Yj)−Y(Xj))∂j[X, Y] = \sum_j \big(X(Y^j) - Y(X^j)\big) \partial_j in coordinates; it makes the vector fields a Lie algebra.

Theorem— Theorem 2 (fundamental theorem on flows; stated)

Every smooth vector field XX on MM has a unique maximal flow: a smooth map θ\theta, defined on an open subset of R×M\R \times M containing {0}×M\{0\} \times M, such that t↦θt(p)t \mapsto \theta_t(p) is the maximal integral curve of XX starting at pp. If MM is compact (or XX has compact support), the flow is defined on all of R×M\R \times M, and t↦θtt \mapsto \theta_t is a one-parameter group of diffeomorphisms: θs+t=θs∘θt\theta_{s+t} = \theta_s \circ \theta_t, θ0=id\theta_0 = \mathrm{id}.

This is the existence and uniqueness theorem for ODEs, applied in charts and glued using uniqueness (and Hausdorffness); see Lee, Introduction to Smooth Manifolds, Chapter 9. Example: on S2⊂R3S^2 \subset \R^3 the field X(x,y,z)=(−y,x,0)X(x, y, z) = (-y, x, 0) is tangent to the sphere, its flow rotates the sphere about the vertical axis, and it vanishes at the two poles. On S3⊂C2S^3 \subset \C^2 the field X(z)=izX(z) = iz vanishes nowhere; its flow z↦eitzz \mapsto e^{it}z is the Hopf flow.

DefinitionsH

Covectors and differential forms

Definition— cotangent space, differential of a function

The cotangent space is the dual Tp∗M=(TpM)∗T_p^*M = (T_pM)^*. A smooth function ff has a differential dfp∈Tp∗Mdf_p \in T_p^*M, dfp(v)=v(f)df_p(v) = v(f). In a chart, the dxi∣pdx^i|_p form the basis dual to the ∂/∂xi∣p\partial/\partial x^i|_p, and df=∑i(∂f/∂xi) dxidf = \sum_i (\partial f/\partial x^i)\, dx^i.

Covectors under a change of chart. From dx~j=∑i(∂x~j/∂xi) dxid\tilde x^j = \sum_i (\partial \tilde x^j / \partial x^i)\, dx^i, a covector ω=∑iωi dxi=∑jω~j dx~j\omega = \sum_i \omega_i\, dx^i = \sum_j \tilde \omega_j\, d\tilde x^j has ωi=∑jω~j ∂x~j/∂xi\omega_i = \sum_j \tilde\omega_j\, \partial \tilde x^j/\partial x^i, that is

ω=JTω~,ω~=(JT)−1ω.\omega = J^{\mathsf T} \tilde\omega, \qquad \tilde \omega = (J^{\mathsf T})^{-1} \omega .

Vectors transform by JJ, covectors by the inverse transpose, and the pairing ω(v)=∑iωivi=ω~TJv\omega(v) = \sum_i \omega_i v^i = \tilde\omega^{\mathsf T} J v is the same in both charts. This is why the classical literature calls vectors "contravariant" and covectors "covariant", and it is why the gradient of a function is not a vector until a metric is chosen: dfdf is a covector, and turning it into a vector needs an inner product on each tangent space (see Riemannian geometry).

Definition— differential forms

A differential kk-form is a smooth section of the bundle ΛkT∗M\Lambda^k T^*M of alternating kk-linear forms on tangent spaces. In a chart, ω=∑i1<⋯<ikωi1…ik dxi1∧⋯∧dxik\omega = \sum_{i_1 < \dots < i_k} \omega_{i_1 \dots i_k}\, dx^{i_1} \wedge \dots \wedge dx^{i_k} with smooth coefficients. Write Ωk(M)\Omega^k(M) for the space of kk-forms. The wedge product is graded commutative: α∧β=(−1)klβ∧α\alpha \wedge \beta = (-1)^{kl}\beta \wedge \alpha for α∈Ωk\alpha \in \Omega^k, β∈Ωl\beta \in \Omega^l.

Definition— exterior derivative and pullback

The exterior derivative is the unique family of linear maps d:Ωk(M)→Ωk+1(M)d : \Omega^k(M) \to \Omega^{k+1}(M) such that dfdf is the differential of ff on functions, d(α∧β)=dα∧β+(−1)kα∧dβd(\alpha \wedge \beta) = d\alpha \wedge \beta + (-1)^k \alpha \wedge d\beta, and d∘d=0d \circ d = 0. In a chart, d(∑IωI dxI)=∑IdωI∧dxId\big(\sum_I \omega_I\, dx^I\big) = \sum_I d\omega_I \wedge dx^I. For smooth F:M→NF : M \to N the pullback F∗:Ωk(N)→Ωk(M)F^* : \Omega^k(N) \to \Omega^k(M) is (F∗ω)p(v1,…,vk)=ωF(p)(dFpv1,…,dFpvk)(F^*\omega)_p(v_1, \dots, v_k) = \omega_{F(p)}(dF_p v_1, \dots, dF_p v_k).

Three facts do most of the work: F∗(α∧β)=F∗α∧F∗βF^*(\alpha \wedge \beta) = F^*\alpha \wedge F^*\beta; F∗(df)=d(f∘F)F^*(df) = d(f \circ F); and F∗∘d=d∘F∗F^* \circ d = d \circ F^*. The identity d(df)=∑i,j∂j∂if dxj∧dxi=0d(df) = \sum_{i, j} \partial_j \partial_i f\, dx^j \wedge dx^i = 0 is the symmetry of second derivatives meeting the antisymmetry of ∧\wedge.

Example— Polar coordinates

Let F(r,θ)=(rcos⁡θ,rsin⁡θ)F(r, \theta) = (r\cos\theta, r\sin\theta). Then F∗dx=cos⁡θ dr−rsin⁡θ dθF^*dx = \cos\theta\, dr - r\sin\theta\, d\theta and F∗dy=sin⁡θ dr+rcos⁡θ dθF^*dy = \sin\theta\, dr + r\cos\theta\, d\theta, so F∗(dx∧dy)=r dr∧dθF^*(dx \wedge dy) = r\, dr \wedge d\theta (the Jacobian appears automatically) and

F∗(x dy−y dxx2+y2)=r2 dθr2=dθ.F^*\left( \frac{x\, dy - y\, dx}{x^2 + y^2} \right) = \frac{r^2\, d\theta}{r^2} = d\theta .

The form ω=(x dy−y dx)/(x2+y2)\omega = (x\,dy - y\,dx)/(x^2+y^2) on R2∖{0}\R^2 \setminus \{0\} is closed (dω=0d\omega = 0) and is locally dd of an angle function — but there is no global angle function, and we will see that ω\omega is not exact.

Definition— orientation and integration

An atlas is oriented if all its transition maps have positive Jacobian determinant; MM is orientable if it has one (equivalently, a nowhere-vanishing nn-form). For an nn-form ω\omega with compact support in an oriented chart (U,φ)(U, \varphi), write (φ−1)∗ω=f dx1∧⋯∧dxn(\varphi^{-1})^*\omega = f\, dx^1 \wedge \dots \wedge dx^n and set ∫Mω=∫Rnf dx1⋯dxn\int_M \omega = \int_{\R^n} f\, dx^1 \cdots dx^n. In general, ∫Mω=∑α∫Mψαω\int_M \omega = \sum_\alpha \int_M \psi_\alpha \omega for a partition of unity {ψα}\{\psi_\alpha\} subordinate to an oriented atlas.

This is well defined because, under a change of chart, the form's coefficient changes by det⁡J\det J while the change-of-variables formula for integrals uses ∣det⁡J∣|\det J|; the two agree exactly when det⁡J>0\det J > 0. That is why we need an oriented atlas to integrate forms. The sphere is orientable: the stereographic charts have det⁡J<0\det J < 0, so compose one of them with a reflection. RP2\R P^2 is not orientable, since it contains a Möbius band; in general RPn\R P^n is orientable exactly when nn is odd.

ProofH

The regular value theorem and embeddings

Definition— regular value

Let F:M→NF : M \to N be smooth. A point p∈Mp \in M is a regular point if dFpdF_p is surjective; c∈Nc \in N is a regular value if every point of F−1(c)F^{-1}(c) is regular (in particular if F−1(c)=∅F^{-1}(c) = \varnothing). A subset S⊆MS \subseteq M is an embedded kk-dimensional submanifold if each p∈Sp \in S has a chart (U,φ)(U, \varphi) of MM with φ(U∩S)=φ(U)∩(Rk×{0})\varphi(U \cap S) = \varphi(U) \cap (\R^k \times \{0\}) (a slice chart).

Theorem— Theorem 3 (regular value theorem)

Let F:M→NF : M \to N be smooth, dim⁡M=m\dim M = m, dim⁡N=n\dim N = n, and let cc be a regular value. Then F−1(c)F^{-1}(c) is an embedded submanifold of MM of dimension m−nm - n (or empty), and TpF−1(c)=ker⁡dFpT_p F^{-1}(c) = \ker dF_p.

Proof

The statement is local. Fix p∈F−1(c)p \in F^{-1}(c) and charts centred at pp and at cc; in them FF becomes a smooth map F^:W→Rn\hat F : W \to \R^n on an open W∋0W \ni 0 of Rm\R^m, with F^(0)=0\hat F(0) = 0 and DF^(0)D\hat F(0) of rank nn (so m≥nm \ge n). After reordering coordinates, write x=(x′,x′′)∈Rm−n×Rnx = (x', x'') \in \R^{m-n} \times \R^n with ∂F^/∂x′′(0)\partial \hat F/\partial x''(0) invertible. Define G(x′,x′′)=(x′,F^(x′,x′′))G(x', x'') = \big(x', \hat F(x', x'')\big). Its Jacobian at 00 is block lower-triangular,

DG(0)=(I0∂F^/∂x′∂F^/∂x′′),DG(0) = \begin{pmatrix} I & 0 \\ \partial \hat F/\partial x' & \partial \hat F/\partial x'' \end{pmatrix},

hence invertible, so by the inverse function theorem GG is a diffeomorphism from a neighbourhood W0W_0 of 00 onto a neighbourhood of 00. In the new chart GG (composed with the old one) the map FF reads F^∘G−1(y′,y′′)=y′′\hat F \circ G^{-1}(y', y'') = y'', a linear projection, and F−1(c)∩W0F^{-1}(c) \cap W_0 corresponds to the slice {y′′=0}\{y'' = 0\}. So F−1(c)F^{-1}(c) has a slice chart of dimension m−nm - n at every point. Finally FF is constant on S=F−1(c)S = F^{-1}(c), so dFpdF_p vanishes on TpST_pS; thus TpS⊆ker⁡dFpT_pS \subseteq \ker dF_p, and both have dimension m−nm - n (rank–nullity), so they are equal. ∎

∎

Theorem— Corollary 4 (spheres and matrix groups are manifolds)

  1. SnS^n is an embedded nn-dimensional submanifold of Rn+1\R^{n+1}, and TxSn=x⊥T_xS^n = x^\perp.
  2. O(n)={A:ATA=I}O(n) = \{A : A^{\mathsf T} A = I\} is an embedded submanifold of Mn(R)M_n(\R) of dimension n(n−1)/2n(n-1)/2, and TIO(n)T_I O(n) is the space of skew-symmetric matrices.
  3. SL(n,R)={A:det⁡A=1}SL(n, \R) = \{A : \det A = 1\} is an embedded submanifold of dimension n2−1n^2 - 1, and TISL(n,R)T_I SL(n, \R) is the space of trace-zero matrices.

Proof

(1) F(x)=∣x∣2F(x) = |x|^2 has dFx(h)=2⟨x,h⟩dF_x(h) = 2\langle x, h\rangle, which is onto R\R whenever x≠0x \ne 0; so 11 is a regular value, and TxSn=ker⁡dFx=x⊥T_x S^n = \ker dF_x = x^\perp.

(2) Let Symn\mathrm{Sym}_n be the space of symmetric matrices, of dimension n(n+1)/2n(n+1)/2, and F:Mn(R)→SymnF : M_n(\R) \to \mathrm{Sym}_n, F(A)=ATAF(A) = A^{\mathsf T}A. Expanding,

(A+tH)T(A+tH)=ATA+t(ATH+HTA)+t2HTH,(A + tH)^{\mathsf T}(A + tH) = A^{\mathsf T}A + t\big(A^{\mathsf T}H + H^{\mathsf T}A\big) + t^2 H^{\mathsf T}H,

so dFA(H)=ATH+HTAdF_A(H) = A^{\mathsf T}H + H^{\mathsf T}A. If A∈O(n)A \in O(n) and S∈SymnS \in \mathrm{Sym}_n, then H=12ASH = \tfrac12 AS satisfies ATH+HTA=12S+12STATA=SA^{\mathsf T}H + H^{\mathsf T}A = \tfrac12 S + \tfrac12 S^{\mathsf T}A^{\mathsf T}A = S. So II is a regular value, dim⁡O(n)=n2−n(n+1)/2=n(n−1)/2\dim O(n) = n^2 - n(n+1)/2 = n(n-1)/2, and TIO(n)=ker⁡dFI={H:H+HT=0}T_I O(n) = \ker dF_I = \{H : H + H^{\mathsf T} = 0\}.

(3) For A∈SL(n,R)A \in SL(n, \R), det⁡((1+t)A)=(1+t)n\det\big((1+t)A\big) = (1 + t)^n, so d(det⁡)A(A)=n≠0d(\det)_A(A) = n \ne 0 and d(det⁡)Ad(\det)_A is onto R\R. Hence 11 is a regular value and dim⁡SL(n,R)=n2−1\dim SL(n, \R) = n^2 - 1. Expanding det⁡(I+tH)=1+ttr⁡H+O(t2)\det(I + tH) = 1 + t \operatorname{tr} H + O(t^2) shows d(det⁡)I(H)=tr⁡Hd(\det)_I(H) = \operatorname{tr} H, so TISL(n,R)T_I SL(n, \R) consists of the trace-zero matrices. ∎

∎

Note the codomain in (2): ATAA^{\mathsf T}A is always symmetric, so as a map into all of Mn(R)M_n(\R) the differential is never surjective, and the theorem would not apply. These tangent spaces at the identity are the Lie algebras so(n)\mathfrak{so}(n) and sl(n,R)\mathfrak{sl}(n, \R) of Lie groups.

Theorem— Theorem 5 (Whitney embedding, compact case)

Every compact smooth manifold MM admits a smooth embedding into some RN\R^N.

Proof

For each pp choose a chart (Up,φp)(U_p, \varphi_p) and a smooth bump function ρp:M→[0,1]\rho_p : M \to [0, 1] supported in UpU_p and equal to 11 on an open Vp∋pV_p \ni p. By compactness finitely many V1,…,VkV_1, \dots, V_k cover MM. Let F=(ρ1φ1,…,ρkφk,ρ1,…,ρk):M→Rk(n+1)F = (\rho_1\varphi_1, \dots, \rho_k\varphi_k, \rho_1, \dots, \rho_k) : M \to \R^{k(n+1)}, where ρiφi\rho_i\varphi_i is extended by 00 outside UiU_i (it is smooth). Injective: if F(p)=F(q)F(p) = F(q) and p∈Vip \in V_i, then ρi(q)=ρi(p)=1\rho_i(q) = \rho_i(p) = 1, so q∈Uiq \in U_i, and φi(p)=φi(q)\varphi_i(p) = \varphi_i(q), so p=qp = q. Immersion: near p∈Vip \in V_i the component ρiφi\rho_i\varphi_i equals the chart φi\varphi_i, whose differential is injective, so dFpdF_p is injective. An injective continuous map from a compact space to a Hausdorff space is a homeomorphism onto its image, so FF is an embedding. ∎

∎

Whitney (1936, 1944) proved much more: every smooth nn-manifold (second countable and Hausdorff, compact or not) embeds in R2n\R^{2n}, and 2n2n cannot be lowered in general (RP2\R P^2 does not embed in R3\R^3). But no embedding is canonical, and nothing intrinsic depends on one.

ProofH

Stokes's theorem

Let MM be an oriented nn-manifold with boundary. Give ∂M\partial M the Stokes orientation: a basis (v1,…,vn−1)(v_1, \dots, v_{n-1}) of Tp∂MT_p\partial M is positive when (ν,v1,…,vn−1)(\nu, v_1, \dots, v_{n-1}) is positive in TpMT_pM for an outward-pointing ν\nu. For Hn\mathbb{H}^n the outward vector is −∂n-\partial_n, and since (−∂n,∂1,…,∂n−1)(-\partial_n, \partial_1, \dots, \partial_{n-1}) has orientation sign (−1)⋅(−1)n−1=(−1)n(-1)\cdot(-1)^{n-1} = (-1)^n, the Stokes orientation of ∂Hn=Rn−1\partial\mathbb{H}^n = \R^{n-1} is (−1)n(-1)^n times the standard one.

Theorem— Theorem 6 (Stokes's theorem)

Let MM be an oriented smooth nn-manifold with boundary, with ∂M\partial M carrying the Stokes orientation, and let ω\omega be a compactly supported smooth (n−1)(n-1)-form on MM. Then

∫Mdω=∫∂Mω,\int_M d\omega = \int_{\partial M} \omega ,

where the right-hand side is 00 if ∂M=∅\partial M = \varnothing.

Proof

Step 1: the half-space. Let M=HnM = \mathbb{H}^n and ω=∑i=1nωi dx1∧⋯∧dxi^∧⋯∧dxn\omega = \sum_{i=1}^n \omega_i\, dx^1 \wedge \dots \wedge \widehat{dx^i} \wedge \dots \wedge dx^n, supported in [−R,R]n−1×[0,R)[-R, R]^{n-1} \times [0, R). Then dω=∑i(−1)i−1∂iωi dx1∧⋯∧dxnd\omega = \sum_i (-1)^{i-1} \partial_i\omega_i\, dx^1 \wedge \dots \wedge dx^n, so

∫Hndω=∑i=1n(−1)i−1∫Hn∂iωi dx.\int_{\mathbb{H}^n} d\omega = \sum_{i=1}^n (-1)^{i-1} \int_{\mathbb{H}^n} \partial_i \omega_i \, dx .

For i<ni < n, integrate first in xix^i over [−R,R][-R, R] (Fubini): by the fundamental theorem of calculus the result is ωi(…,R,… )−ωi(…,−R,… )=0\omega_i(\dots, R, \dots) - \omega_i(\dots, -R, \dots) = 0. For i=ni = n, ∫0R∂nωn dxn=−ωn(x′,0)\int_0^R \partial_n\omega_n\, dx^n = -\omega_n(x', 0) since ωn\omega_n vanishes at xn=Rx^n = R. Hence

∫Hndω=(−1)n−1⋅(−1)∫Rn−1ωn(x′,0) dx′=(−1)n∫Rn−1ωn(x′,0) dx′.\int_{\mathbb{H}^n} d\omega = (-1)^{n-1}\cdot(-1) \int_{\R^{n-1}} \omega_n(x', 0)\, dx' = (-1)^n \int_{\R^{n-1}} \omega_n(x', 0)\, dx' .

On ∂Hn\partial \mathbb{H}^n the form dxndx^n restricts to 00, so only the i=ni = n term survives: ω∣∂Hn=ωn(x′,0) dx1∧⋯∧dxn−1\omega|_{\partial\mathbb{H}^n} = \omega_n(x', 0)\, dx^1 \wedge \dots \wedge dx^{n-1}. Integrating with the Stokes orientation, which is (−1)n(-1)^n times the standard one, gives exactly the same number. The same computation with Rn\R^n in place of Hn\mathbb{H}^n gives ∫Rndω=0\int_{\R^n} d\omega = 0, because now every term vanishes.

Step 2: gluing. Cover the compact set supp⁡ω\operatorname{supp}\omega by finitely many oriented charts, each onto an open subset of Rn\R^n or of Hn\mathbb{H}^n, and let {ψα}\{\psi_\alpha\} be a subordinate partition of unity. Since ∑αψα=1\sum_\alpha \psi_\alpha = 1 near supp⁡ω\operatorname{supp}\omega, we have ∑αdψα=0\sum_\alpha d\psi_\alpha = 0 there, so

∑αd(ψαω)=∑αdψα∧ω+∑αψα dω=dω.\sum_\alpha d(\psi_\alpha\omega) = \sum_\alpha d\psi_\alpha \wedge \omega + \sum_\alpha \psi_\alpha\, d\omega = d\omega .

Each ψαω\psi_\alpha\omega is supported in a single chart. Pulling back by that chart preserves integrals (the chart is orientation-preserving), commutes with dd, and carries ∂M\partial M to ∂Hn\partial\mathbb{H}^n with matching Stokes orientations, so Step 1 gives ∫Md(ψαω)=∫∂Mψαω\int_M d(\psi_\alpha \omega) = \int_{\partial M} \psi_\alpha\omega. Summing over α\alpha proves the theorem. ∎

∎

For n=1n = 1 and M=[a,b]M = [a, b] this is the fundamental theorem of calculus (the boundary is {b}\{b\} with sign ++ and {a}\{a\} with sign −-). For domains in R2\R^2 and R3\R^3 it contains Green's theorem, the classical Stokes theorem and the divergence theorem. On a closed manifold (∂M=∅\partial M = \varnothing, MM compact), every exact nn-form integrates to 00; so a volume form on a closed oriented manifold is never exact.

ProofH

de Rham cohomology

Since d∘d=0d \circ d = 0, every exact form (ω=dη\omega = d\eta) is closed (dω=0d\omega = 0). The failure of the converse measures topology.

Definition— de Rham cohomology

HdRk(M)=ker⁡(d:Ωk(M)→Ωk+1(M))im⁡(d:Ωk−1(M)→Ωk(M))H^k_{\mathrm{dR}}(M) = \dfrac{\ker\big(d : \Omega^k(M) \to \Omega^{k+1}(M)\big)}{\operatorname{im}\big(d : \Omega^{k-1}(M) \to \Omega^k(M)\big)}. A smooth F:M→NF : M \to N induces F∗:HdRk(N)→HdRk(M)F^* : H^k_{\mathrm{dR}}(N) \to H^k_{\mathrm{dR}}(M), since F∗F^* commutes with dd; and the wedge product makes HdR∙(M)H^\bullet_{\mathrm{dR}}(M) a graded-commutative algebra.

HdR0(M)H^0_{\mathrm{dR}}(M) is the space of locally constant functions (df=0df = 0 means ff is constant on each component), so dim⁡HdR0(M)\dim H^0_{\mathrm{dR}}(M) is the number of connected components.

Theorem— Theorem 7 (the circle)

Integration over the circle induces an isomorphism HdR1(S1)≅RH^1_{\mathrm{dR}}(S^1) \cong \R, [ω]↦∫S1ω[\omega] \mapsto \int_{S^1}\omega. In particular dθd\theta (the form x dy−y dxx\,dy - y\,dx restricted to S1S^1) is closed but not exact.

Proof

Let π:R→S1\pi : \R \to S^1, π(θ)=(cos⁡θ,sin⁡θ)\pi(\theta) = (\cos\theta, \sin\theta); it is a local diffeomorphism. Every 11-form on the 11-manifold S1S^1 is closed, and π∗ω=f(θ) dθ\pi^*\omega = f(\theta)\, d\theta with ff smooth and 2π2\pi-periodic; ω=dg\omega = dg for some g∈C∞(S1)g \in C^\infty(S^1) exactly when f=G′f = G' for the 2π2\pi-periodic function G=g∘πG = g \circ \pi (a smooth periodic GG descends to a smooth gg on S1S^1, because π\pi is a local diffeomorphism). Also ∫S1ω=∫02πf(θ) dθ\int_{S^1}\omega = \int_0^{2\pi} f(\theta)\, d\theta.

Exact forms integrate to zero: if f=G′f = G' with GG periodic, then ∫02πf=G(2π)−G(0)=0\int_0^{2\pi} f = G(2\pi) - G(0) = 0.

Forms with zero integral are exact: if ∫02πf=0\int_0^{2\pi} f = 0, let G(θ)=∫0θfG(\theta) = \int_0^\theta f. Then G′=fG' = f, and G(θ+2π)−G(θ)=∫θθ+2πf=∫02πf=0G(\theta + 2\pi) - G(\theta) = \int_\theta^{\theta + 2\pi} f = \int_0^{2\pi} f = 0, because the integral of a periodic function over a full period does not depend on where the period starts. So GG is periodic and ω=dg\omega = dg.

Thus the linear map ω↦∫S1ω\omega \mapsto \int_{S^1}\omega has kernel exactly the exact forms. It is onto, since ∫S1dθ=2π≠0\int_{S^1} d\theta = 2\pi \ne 0, so it induces an isomorphism HdR1(S1)→RH^1_{\mathrm{dR}}(S^1) \to \R. ∎

∎

The same argument applied to ω=(x dy−y dx)/(x2+y2)\omega = (x\,dy - y\,dx)/(x^2+y^2) on R2∖{0}\R^2 \setminus \{0\} shows it is not exact: its integral round the unit circle is 2π2\pi, whereas an exact form integrates to 00 round every closed curve. More generally (proofs in Bott–Tu or Lee, Chapters 17–18):

  • Poincaré lemma and homotopy invariance: homotopic maps induce the same map on HdR∙H^\bullet_{\mathrm{dR}}; hence HdRk(Rn)=0H^k_{\mathrm{dR}}(\R^n) = 0 for k≥1k \ge 1, and every closed form is locally exact.
  • Computations (by Mayer–Vietoris): HdRk(Sn)≅RH^k_{\mathrm{dR}}(S^n) \cong \R for k=0,nk = 0, n and 00 otherwise; HdRk(Tn)≅R(nk)H^k_{\mathrm{dR}}(T^n) \cong \R^{\binom nk}; HdR1(R2∖{p1,…,pm})≅RmH^1_{\mathrm{dR}}(\R^2 \setminus \{p_1, \dots, p_m\}) \cong \R^m, detected by the periods around the mm punctures.
  • Top degree: for a compact connected oriented nn-manifold without boundary, ∫M:HdRn(M)→R\int_M : H^n_{\mathrm{dR}}(M) \to \R is an isomorphism (Stokes gives well-definedness; surjectivity is clear; injectivity is the real content).

Theorem— Theorem 8 (de Rham's theorem; stated)

For every smooth manifold MM, integration of forms over smooth singular chains induces an isomorphism HdRk(M)≅Hk(M;R)H^k_{\mathrm{dR}}(M) \cong H^k(M; \R) with singular cohomology with real coefficients, compatible with products.

Stokes's theorem, ∫cdω=∫∂cω\int_c d\omega = \int_{\partial c}\omega, is precisely the statement that integration is a chain map from the de Rham complex to the singular cochain complex; the theorem (de Rham, 1931) says it is a quasi-isomorphism. So de Rham cohomology is a topological invariant, even though it is defined with derivatives. See algebraic topology.

Formal mathematicsF

Smooth manifolds in Lean and Mathlib

Mathlib's manifolds follow the definitions above, with one generalisation: a manifold is modelled on a model with corners I:H→EI : H \to E, which lets the same definitions cover manifolds without boundary (𝓡 n, modelled on Rn\R^n), with boundary (𝓡∂ n, on a half-space) and with corners. ChartedSpace H M is an atlas, IsManifold I n M says the transition maps are CnC^n (with ω meaning analytic), ContMDiff is smoothness of maps, and mfderiv is the differential.

First, the stereographic charts. We define the inverse projection from the north pole and the projection from the south pole on E×RE \times \R for any real inner product space EE, check that the first lands on the sphere, and prove that the transition map is the inversion, which is smooth on E∖{0}E \setminus \{0\} and its own inverse. These are the computations of the sphere example, line by line:

stereoInvNorthLean 4Formally verified
/-- **Inverse stereographic projection from the north pole.** The point `u` of the equatorial
space `E` is sent to the point `(x, t) ∈ E × ℝ` where the line from the north pole `(0, 1)`
through `(u, 0)` meets the unit sphere again:
`x = 2u / (‖u‖² + 1)`, `t = (‖u‖² - 1) / (‖u‖² + 1)`. -/
noncomputable def stereoInvNorth (u : E) : E × ℝ :=
  ((2 / (‖u‖ ^ 2 + 1)) • u, (‖u‖ ^ 2 - 1) / (‖u‖ ^ 2 + 1))
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor
stereoInvNorth_mem_sphereLean 4Proved here from basic factsFormally verified
/-- The inverse stereographic projection really lands on the unit sphere `‖x‖² + t² = 1`.
With `r = ‖u‖²` the claim is `4r + (r - 1)² = (r + 1)²`, an identity of polynomials. -/
theorem stereoInvNorth_mem_sphere (u : E) :
    ‖(stereoInvNorth u).1‖ ^ 2 + (stereoInvNorth u).2 ^ 2 = 1 := by
  have hpos : (0 : ℝ) < ‖u‖ ^ 2 + 1 := by positivity
  simp only [stereoInvNorth, norm_smul, Real.norm_eq_abs, abs_div, abs_two,
    abs_of_pos hpos, mul_pow, div_pow]
  field_simp
  ring
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor
stereoSouthLean 4Formally verified
/-- **Stereographic projection from the south pole** `(0, -1)`: the point `(x, t)` of the sphere
(with `t ≠ -1`) is sent to `x / (1 + t)`, where the line from the south pole through `(x, t)`
meets the equatorial space. -/
noncomputable def stereoSouth (p : E × ℝ) : E :=
  (1 / (1 + p.2)) • p.1
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor
inversionLean 4Formally verified
/-- The **inversion** `u ↦ u / ‖u‖²` in the unit sphere of `E`, defined on `E ∖ {0}`. -/
noncomputable def inversion (u : E) : E :=
  (‖u‖ ^ 2)⁻¹ • u
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor
transition_mapLean 4Proved here from basic factsFormally verified
/-- **The transition map between the two stereographic charts is the inversion.** For `u ≠ 0`,
going up to the sphere by the north chart and down again by the south chart gives
`u / ‖u‖²`. This is the computation in the lesson:
`(2u/(r+1)) / (1 + (r-1)/(r+1)) = (2u/(r+1)) / (2r/(r+1)) = u/r`, with `r = ‖u‖²`. -/
theorem transition_map {u : E} (hu : u ≠ 0) :
    stereoSouth (stereoInvNorth u) = inversion u := by
  have hr : ‖u‖ ^ 2 ≠ 0 := by positivity
  have hpos : ‖u‖ ^ 2 + 1 ≠ 0 := by positivity
  simp only [stereoSouth, stereoInvNorth, inversion, smul_smul]
  congr 1
  field_simp
  ring
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor
inversion_inversionLean 4Proved here from basic factsFormally verified
/-- The inversion is its own inverse on `E ∖ {0}`: `‖u / ‖u‖²‖ = 1/‖u‖`, so inverting again
multiplies by `‖u‖²` and returns `u`. So the transition map is a bijection of `E ∖ {0}`, as the
transition map between two charts must be. -/
theorem inversion_inversion {u : E} (hu : u ≠ 0) : inversion (inversion u) = u := by
  have hn : ‖u‖ ≠ 0 := norm_ne_zero_iff.mpr hu
  have key : ((‖(‖u‖ ^ 2)⁻¹ • u‖ ^ 2)⁻¹ * (‖u‖ ^ 2)⁻¹ : ℝ) = 1 := by
    rw [norm_smul, norm_inv, norm_pow, norm_norm]
    field_simp
  simp only [inversion, smul_smul, key, one_smul]
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor
inversion_smoothLean 4Proved here from basic factsFormally verified
/-- The transition map is **smooth** (`C^∞`) on `E ∖ {0}`: it is the product of the smooth
function `u ↦ (‖u‖²)⁻¹` (smooth where `‖u‖² ≠ 0`) and the identity. Together with
`transition_map` this is what makes the two stereographic charts a smooth atlas. -/
theorem inversion_smooth : ContDiffOn ℝ ∞ (inversion (E := E)) {u | u ≠ 0} := by
  have h1 : ContDiffOn ℝ ∞ (fun u : E => (‖u‖ ^ 2)⁻¹) {u | u ≠ 0} :=
    (contDiff_norm_sq ℝ).contDiffOn.inv fun u hu => by
      have : u ≠ 0 := hu
      positivity
  exact h1.smul contDiffOn_id
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor

Mathlib uses exactly these charts to make the unit sphere of Rn+1\R^{n+1} a manifold, and we cite its instance:

sphere_is_manifoldLean 4Applies a library theoremFormally verified
/-- **The `n`-sphere is a smooth manifold** (indeed analytic: `ω`), modelled on `ℝⁿ`, with the
stereographic charts. This is Mathlib's instance `EuclideanSpace.instIsManifoldSphere`. -/
theorem sphere_is_manifold (n : ℕ) :
    IsManifold (𝓡 n) ω (Metric.sphere (0 : EuclideanSpace ℝ (Fin (n + 1))) 1) :=
  inferInstance
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor

Next, the two algebraic facts behind Corollary 4(2): the expansion of (A+tH)T(A+tH)(A + tH)^{\mathsf T}(A + tH), whose linear term is dFA(H)dF_A(H), and the surjectivity of dFAdF_A onto symmetric matrices at an orthogonal AA:

gram_expansionLean 4Proved here from basic factsFormally verified
/-- Expanding `F(A) = AᵀA` along the line `A + tH`:
`(A + tH)ᵀ(A + tH) = AᵀA + t(AᵀH + HᵀA) + t²HᵀH`. The coefficient of `t` is the derivative
`dF_A(H) = AᵀH + HᵀA`. -/
theorem gram_expansion (A H : Matrix n n ℝ) (t : ℝ) :
    (A + t • H)ᵀ * (A + t • H) = Aᵀ * A + t • (Aᵀ * H + Hᵀ * A) + t ^ 2 • (Hᵀ * H) := by
  simp only [Matrix.transpose_add, Matrix.transpose_smul, Matrix.add_mul, Matrix.mul_add,
    Matrix.smul_mul, Matrix.mul_smul, smul_add, smul_smul, pow_two]
  abel
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor
orthogonal_regular_valueLean 4Proved here from basic factsFormally verified
/-- **The identity is a regular value of `A ↦ AᵀA`.** At an orthogonal matrix `A`, the derivative
`H ↦ AᵀH + HᵀA` maps onto the symmetric matrices: for symmetric `S`, the matrix `H = ½ A S` is a
preimage, because `Aᵀ(½AS) + (½AS)ᵀA = ½S + ½SᵀAᵀA = S`. By the regular value theorem,
`O(n) = F⁻¹(I)` is a manifold of dimension `n² - n(n+1)/2 = n(n-1)/2`. -/
theorem orthogonal_regular_value [DecidableEq n] {A S : Matrix n n ℝ} (hA : Aᵀ * A = 1) (hS : Sᵀ = S) :
    Aᵀ * ((1 / 2 : ℝ) • (A * S)) + ((1 / 2 : ℝ) • (A * S))ᵀ * A = S := by
  rw [Matrix.transpose_smul, Matrix.transpose_mul, hS, Matrix.mul_smul, Matrix.smul_mul,
    ← Matrix.mul_assoc, hA, Matrix.mul_assoc, hA, Matrix.one_mul, Matrix.mul_one, ← add_smul]
  norm_num
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor

The chain rule for manifolds, i.e. functoriality of the tangent map:

chain_ruleLean 4Applies a library theoremFormally verified
/-- **The chain rule on manifolds**: the tangent map of a composite is the composite of the
tangent maps, `T(g ∘ f) = Tg ∘ Tf`. In the language of categories, `T` is a functor. This is
Mathlib's `mfderiv_comp`. -/
theorem chain_rule {E : Type*} [NormedAddCommGroup E] [NormedSpace ℝ E]
    {H : Type*} [TopologicalSpace H] {I : ModelWithCorners ℝ E H}
    {M : Type*} [TopologicalSpace M] [ChartedSpace H M]
    {E' : Type*} [NormedAddCommGroup E'] [NormedSpace ℝ E']
    {H' : Type*} [TopologicalSpace H'] {I' : ModelWithCorners ℝ E' H'}
    {M' : Type*} [TopologicalSpace M'] [ChartedSpace H' M']
    {E'' : Type*} [NormedAddCommGroup E''] [NormedSpace ℝ E'']
    {H'' : Type*} [TopologicalSpace H''] {I'' : ModelWithCorners ℝ E'' H''}
    {M'' : Type*} [TopologicalSpace M''] [ChartedSpace H'' M'']
    {f : M → M'} {g : M' → M''} (x : M)
    (hg : MDifferentiableAt I' I'' g (f x)) (hf : MDifferentiableAt I I' f x) :
    mfderiv I I'' (g ∘ f) x = (mfderiv I' I'' g (f x)).comp (mfderiv I I' f x) :=
  mfderiv_comp x hg hf
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor

Finally Theorem 7, in the coordinate θ\theta: a 11-form on the circle R/TZ\R/T\Z is f(θ) dθf(\theta)\,d\theta with ff continuous and TT-periodic, and it is exact when ff is the derivative of a TT-periodic function. Both halves of the proof are formalised — exact forms have period zero by the fundamental theorem of calculus; forms with period zero have the periodic primitive ∫0xf\int_0^x f — and combined into the "if and only if" that says HdR1(S1)≅RH^1_{\mathrm{dR}}(S^1) \cong \R:

exact_integral_zeroLean 4Proved here from basic factsFormally verified
/-- **Exact forms have zero period.** If `g` is `T`-periodic with continuous derivative `f`, then
`∫₀ᵀ f = g(T) - g(0) = 0`, by the fundamental theorem of calculus. -/
theorem exact_integral_zero {T : ℝ} {f g : ℝ → ℝ} (hg : Function.Periodic g T)
    (hderiv : ∀ x, HasDerivAt g (f x) x) (hf : Continuous f) :
    ∫ x in (0 : ℝ)..T, f x = 0 := by
  rw [intervalIntegral.integral_eq_sub_of_hasDerivAt (fun x _ => hderiv x)
    (hf.intervalIntegrable 0 T)]
  have := hg 0
  rw [zero_add] at this
  rw [this, sub_self]
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor
dtheta_not_exactLean 4Proved here from basic factsFormally verified
/-- **`dθ` is closed but not exact on the circle.** There is no `2π`-periodic function `g` with
`g' = 1`: its integral over a period would be `0`, but `∫₀^{2π} 1 dθ = 2π`. -/
theorem dtheta_not_exact :
    ¬ ∃ g : ℝ → ℝ, Function.Periodic g (2 * Real.pi) ∧ ∀ x, HasDerivAt g 1 x := by
  rintro ⟨g, hg, hderiv⟩
  have h := exact_integral_zero (f := fun _ => (1 : ℝ)) hg hderiv continuous_const
  simp only [intervalIntegral.integral_const, sub_zero, smul_eq_mul, mul_one] at h
  exact Real.pi_ne_zero (by linarith)
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor
primitive_periodicLean 4Proved here from basic factsFormally verified
/-- **Forms with zero period are exact.** If `f` is continuous and `T`-periodic with
`∫₀ᵀ f = 0`, then `G(x) = ∫₀ˣ f` is `T`-periodic (because `∫ₓ^{x+T} f` does not depend on `x`)
and `G' = f` (the fundamental theorem of calculus). -/
theorem primitive_periodic {T : ℝ} {f : ℝ → ℝ} (hf : Continuous f)
    (hper : Function.Periodic f T) (hzero : ∫ x in (0 : ℝ)..T, f x = 0) :
    Function.Periodic (fun x => ∫ t in (0 : ℝ)..x, f t) T ∧
      ∀ x, HasDerivAt (fun x => ∫ t in (0 : ℝ)..x, f t) (f x) x := by
  refine ⟨fun x => ?_, fun x => ?_⟩
  · have hsplit : ∫ t in (0 : ℝ)..x + T, f t = (∫ t in (0 : ℝ)..x, f t) + ∫ t in x..x + T, f t :=
      (intervalIntegral.integral_add_adjacent_intervals (hf.intervalIntegrable _ _)
        (hf.intervalIntegrable _ _)).symm
    have hshift : ∫ t in x..x + T, f t = ∫ t in (0 : ℝ)..0 + T, f t :=
      hper.intervalIntegral_add_eq x 0
    simp only
    rw [hsplit, hshift, zero_add, hzero, add_zero]
  · exact intervalIntegral.integral_hasDerivAt_right (hf.intervalIntegrable _ _)
      hf.aestronglyMeasurable.stronglyMeasurableAtFilter hf.continuousAt
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor
de_rham_circleLean 4Proved here from basic factsFormally verified
/-- **`H¹_dR(S¹) ≅ ℝ`, via the period.** A continuous `T`-periodic `f` (a 1-form `f dθ` on the
circle) is the derivative of a `T`-periodic function (is exact) **if and only if** its period
`∫₀ᵀ f` vanishes. So the linear map "integrate over the circle" from closed 1-forms to `ℝ` has
kernel exactly the exact forms; it is onto (the constant `1` has period `T`), so for `T ≠ 0` it
induces an isomorphism `H¹_dR(S¹) ≅ ℝ`. -/
theorem de_rham_circle {T : ℝ} {f : ℝ → ℝ} (hf : Continuous f) (hper : Function.Periodic f T) :
    (∃ g : ℝ → ℝ, Function.Periodic g T ∧ ∀ x, HasDerivAt g (f x) x) ↔
      ∫ x in (0 : ℝ)..T, f x = 0 := by
  constructor
  · rintro ⟨g, hg, hderiv⟩
    exact exact_integral_zero hg hderiv hf
  · intro hzero
    exact ⟨_, primitive_periodic hf hper hzero⟩
Axioms used: propext, Classical.choice, Quot.soundformal/Lumenary/SmoothManifolds.leanTry it in the Lean web editor

What is and is not formalised. Independently proved: the stereographic computations (the image lies on the sphere, the transition map is the inversion, it is an involution and smooth), the expansion and surjectivity behind O(n)O(n), and the circle computation of Theorem 7, in the coordinate on R\R (continuous rather than smooth ff, which is slightly more general). Cited from Mathlib: the manifold structure on spheres and the chain rule. Not formalised: the regular value theorem itself (Mathlib has the inverse and implicit function theorems on Banach spaces, but not yet the theorem for manifolds), differential forms on manifolds, Stokes's theorem on manifolds with boundary (Mathlib has the divergence theorem on boxes), and de Rham's theorem.

ComputationC

Charts, covectors and periods, computed

The first program computes, exactly, the transition map between the two stereographic charts of S2S^2 and its Jacobian. It then moves a tangent vector and a covector from one chart to the other: vector components by JJ, covector components by (JT)−1(J^{\mathsf T})^{-1}. The pairing ω(v)\omega(v) is the same in both charts; moving the covector by JJ (the classic mistake) changes it.

Stereographic charts and the two transformation lawsPythonSymbolic computationExact computationComputationally tested
"""Two charts on the sphere, and the two transformation laws.

The 2-sphere is covered by the stereographic charts from the north and the
south pole. SymPy computes the transition map between them, its Jacobian
matrix J and its determinant, exactly. Then it takes a tangent vector and a
covector at one point, written in the first chart, and transforms them to the
second chart: vector components by J, covector components by the inverse
transpose of J. The pairing "covector applied to vector" comes out the same in
both charts, as it must; transforming the covector by J instead (the classic
mistake) gives a different number. This is exact symbolic computation for one
explicit chart change, not a proof of the general laws.
"""

import sympy

U, V = sympy.symbols("u v", real=True)

def north_inverse(u, v):
    """Inverse stereographic projection from the north pole: R^2 -> S^2."""
    r = u**2 + v**2
    return (2 * u / (r + 1), 2 * v / (r + 1), (r - 1) / (r + 1))

def south_chart(x, y, z):
    """Stereographic projection from the south pole: S^2 minus (0,0,-1) -> R^2."""
    return (x / (1 + z), y / (1 + z))

def transition():
    """The transition map south o north^-1, simplified, as a pair."""
    x, y, z = north_inverse(U, V)
    return tuple(sympy.simplify(c) for c in south_chart(x, y, z))

def jacobian(pair):
    """Jacobian matrix of a pair of expressions in (u, v)."""
    return sympy.Matrix(pair).jacobian([U, V])

def on_sphere(u, v):
    """x^2 + y^2 + z^2 - 1 for the image of (u, v): it should simplify to 0."""
    x, y, z = north_inverse(u, v)
    return sympy.simplify(x**2 + y**2 + z**2 - 1)

def transform(J, vector, covector):
    """New components: vectors by J, covectors by (J^-1)^T, and the wrong rule."""
    new_vector = J * vector
    new_covector = J.inv().T * covector
    wrong_covector = J * covector
    return new_vector, new_covector, wrong_covector

def pairing(covector, vector):
    """The number covector(vector) = sum of products of components."""
    return sympy.nsimplify((covector.T * vector)[0, 0])

def main():
    t = transition()
    J = jacobian(t)
    det = sympy.factor(J.det())
    print("Charts: north^-1(u, v) on the sphere?", on_sphere(U, V) == 0)
    print(f"Transition south o north^-1: {t}")
    print(f"Jacobian determinant: {det}")
    inv = sympy.simplify(jacobian(t).subs({U: t[0], V: t[1]}, simultaneous=True)
                         * J)
    print(f"J(T(p)) * J(p) = identity (T is its own inverse): {inv == sympy.eye(2)}")
    print()
    point = {U: 1, V: 2}
    Jp = J.subs(point)
    vector = sympy.Matrix([3, -1])
    covector = sympy.Matrix([2, 5])
    new_v, new_w, wrong_w = transform(Jp, vector, covector)
    print("At the point (u, v) = (1, 2) of the north chart:")
    print(f"  J = {Jp.tolist()}")
    print(f"  vector   {vector.T.tolist()[0]} -> {new_v.T.tolist()[0]} (by J)")
    print(f"  covector {covector.T.tolist()[0]} -> {new_w.T.tolist()[0]}"
          " (by inverse transpose)")
    print(f"  pairing in chart 1: {pairing(covector, vector)}")
    print(f"  pairing in chart 2: {pairing(new_w, new_v)}")
    print(f"  pairing if the covector is wrongly moved by J:"
          f" {pairing(wrong_w, new_v)}")
    print()
    print("The determinant is negative everywhere: these two charts reverse")
    print("orientation relative to each other. Composing one with (u, v) -> (u, -v)")
    print("gives an oriented atlas, so the sphere is orientable.")
    print()
    print("Epistemic status: exact symbolic computation (SymPy) for one chart")
    print("change and one point. It illustrates, but does not prove, the general")
    print("transformation laws, which follow from the chain rule.")

if __name__ == "__main__":
    main()
Runs entirely in your browser, in an isolated worker. Nothing is sent to a server.

The determinant −1/(u2+v2)2-1/(u^2+v^2)^2 is negative everywhere, as predicted by the sphere example: the two charts are individually fine but orientation-incompatible. This is exact computation for one chart change, an illustration rather than a proof; the general laws are consequences of the chain rule.

The second program checks Stokes's theorem for x3 dyx^3\,dy on the unit disc and computes periods of the closed form dθd\theta round several loops in the punctured plane, and of a combination of two angle forms in the twice-punctured plane.

Stokes's theorem and the periods of closed formsPythonNumerical approximationComputationally tested
"""Stokes's theorem on a disc, and periods of closed 1-forms.

Part 1 checks Stokes's theorem numerically for the 1-form x^3 dy on the unit
disc: the line integral round the circle and the area integral of
d(x^3 dy) = 3x^2 dx^dy both approximate 3*pi/4.

Part 2 integrates the closed form dtheta = (x dy - y dx)/(x^2 + y^2) on the
punctured plane round several loops. Its periods are 2*pi times winding
numbers, never anything else, and round a loop that misses the puncture the
period is 0, as Stokes predicts. With two punctures, a closed form is
determined up to exact forms by its two periods: H^1 of the twice-punctured
plane is R^2. All integrals are computed by quadrature (numerical evidence).
"""

import math

def line_integral(form, curve, n=4000):
    """Integral of the 1-form P dx + Q dy along a closed curve on [0, 2*pi].

    ``form(x, y)`` returns (P, Q); ``curve(t)`` returns (x, y, dx/dt, dy/dt).
    The trapezoid rule on a periodic integrand is very accurate.
    """
    h = 2 * math.pi / n
    total = 0.0
    for k in range(n):
        x, y, dx, dy = curve(k * h)
        p, q = form(x, y)
        total += p * dx + q * dy
    return total * h

def disc_integral(density, n_r=400, n_t=400):
    """Integral of a function over the unit disc, midpoint rule in polar form."""
    hr, ht = 1.0 / n_r, 2 * math.pi / n_t
    total = 0.0
    for i in range(n_r):
        r = (i + 0.5) * hr
        for j in range(n_t):
            t = (j + 0.5) * ht
            total += density(r * math.cos(t), r * math.sin(t)) * r
    return total * hr * ht

def angle_form(a=0.0, b=0.0):
    """dtheta measured from the point (a, b): ((x-a) dy - (y-b) dx) / r^2."""
    def form(x, y):
        dx, dy = x - a, y - b
        r2 = dx * dx + dy * dy
        return (-dy / r2, dx / r2)
    return form

def circle(cx, cy, radius, turns=1):
    """A circle traversed ``turns`` times anticlockwise, as (x, y, x', y')."""
    def curve(t):
        s = turns * t
        return (cx + radius * math.cos(s), cy + radius * math.sin(s),
                -radius * turns * math.sin(s), radius * turns * math.cos(s))
    return curve

def ellipse(a, b):
    """The ellipse (a cos t, b sin t), anticlockwise."""
    return lambda t: (a * math.cos(t), b * math.sin(t),
                      -a * math.sin(t), b * math.cos(t))

def combine(c1, f1, c2, f2):
    """The 1-form c1 f1 + c2 f2."""
    def form(x, y):
        p1, q1 = f1(x, y)
        p2, q2 = f2(x, y)
        return (c1 * p1 + c2 * p2, c1 * q1 + c2 * q2)
    return form

def main():
    exact = 3 * math.pi / 4
    lhs = line_integral(lambda x, y: (0.0, x**3), circle(0, 0, 1))
    rhs = disc_integral(lambda x, y: 3 * x * x)
    print("1. Stokes for w = x^3 dy on the unit disc (exact value 3*pi/4):")
    print(f"   boundary integral  {lhs:.10f}   error {abs(lhs - exact):.1e}")
    print(f"   integral of dw     {rhs:.10f}   error {abs(rhs - exact):.1e}")
    print()
    print("2. Periods of dtheta = (x dy - y dx)/(x^2 + y^2), divided by 2*pi:")
    loops = [
        ("unit circle round 0", circle(0, 0, 1)),
        ("circle of radius 0.01 round 0", circle(0, 0, 0.01)),
        ("ellipse 3 x 0.2 round 0", ellipse(3, 0.2)),
        ("circle round 0, twice", circle(0, 0, 1, turns=2)),
        ("circle centred at (3, 0)", circle(3, 0, 1)),
    ]
    for name, loop in loops:
        ratio = round(line_integral(angle_form(), loop) / (2 * math.pi), 8)
        print(f"   {name:<32} {ratio + 0.0:+.8f}")
    print("   Annulus 0.01 < r < 1: outer minus inner period = 0 = integral of")
    print("   d(dtheta), so Stokes holds although dtheta is not exact.")
    print()
    form = combine(2.0, angle_form(-1, 0), -0.5, angle_form(1, 0))
    print("3. Two punctures at (-1, 0) and (1, 0); w = 2 dth_left - 0.5 dth_right:")
    for name, loop in [("loop round (-1, 0)", circle(-1, 0, 0.5)),
                       ("loop round (1, 0)", circle(1, 0, 0.5)),
                       ("loop round both", circle(0, 0, 3))]:
        print(f"   {name:<20} period / 2*pi = "
              f"{line_integral(form, loop) / (2 * math.pi):+.6f}")
    print("   The two periods (2, -0.5) pin down the class of w: H^1 = R^2.")
    print()
    print("Epistemic status: numerical quadrature. The agreement in part 1 is")
    print("evidence for one form on one disc; the integer periods in part 2 are")
    print("evidence for the theorem that periods of dtheta are 2*pi times winding")
    print("numbers. The proofs are Stokes's theorem and the lesson's argument.")

if __name__ == "__main__":
    main()
Runs entirely in your browser, in an isolated worker. Nothing is sent to a server.

The periods of dθd\theta come out as 2π2\pi times an integer — the winding number — whatever the size or shape of the loop. With two punctures the two periods can be prescribed independently, illustrating HdR1(R2∖{p,q})≅R2H^1_{\mathrm{dR}}(\R^2 \setminus \{p, q\}) \cong \R^2. Quadrature is numerical evidence; the theorems are proved above.

Explore

Explore: closed but not exact, and a glued manifold

The vector field (−y,x)/(x2+y2)(-y, x)/(x^2 + y^2) corresponds (through the Euclidean dot product) to the 11-form dθ=(x dy−y dx)/(x2+y2)d\theta = (x\,dy - y\,dx)/(x^2+y^2). Drag the loop: when it encircles the origin the circulation ∮ω\oint \omega is 2π2\pi; when it does not, it is 00. The curl is 00 everywhere it is defined (ω\omega is closed), so Green's theorem predicts circulation 00 — and it fails exactly when the loop surrounds the singularity, where ω\omega is not defined. No function has this field as its gradient.

dθ: closed but not exact

Circulation 2π round the origin, 0 elsewhere; curl 0 wherever it is defined.

Examples

F(x, y) = (−y/(x² + y²), x/(x² + y²))

Drag the gold centre (or click anywhere) to move the loop. With the keyboard, focus the centre and use the arrow keys.

Turns so far: 0.09

Circulation
∮ F·dr ≈ 6.2832 | ∬ curl F dA ≈ −0.0322
Flux out
∮ F·n ds ≈ 0 | ∬ div F dA ≈ 0.1043
At the centre
curl F ≈ 0, div F ≈ 0, πr² ≈ 8.042
Index
1

All values are numerical approximations (sums around the loop and over the disc).

The two sides differ: the field must fail to be smooth (or be undefined) somewhere inside the loop, so Green’s theorem does not apply to this loop.

Index 1 ≠ 0: walking once around the loop, the arrow makes 1 net full turns (negative means clockwise). A field that is continuous on the whole disc would have to vanish inside; this one is not (see above), and the point where it breaks down carries the index instead.

A manifold need not live in any RN\R^N to begin with. The torus and the Klein bottle below are squares with their sides glued; each point, including points on the seams and corners, has a neighbourhood that is a disc, so both are smooth 22-manifolds. Walk the creature across the seams: on the Klein bottle it comes back mirror-reversed, so no oriented atlas exists.

Manifolds made by gluing

Charts across the seams make these smooth surfaces; only the torus is orientable.

abab

Keys: ↑ forward, ↓ back, ← and → turn, Space walks or stops. You can also drag it.

Torusa b a⁻¹ b⁻¹
Orientable
Yes
Euler characteristic
χ = 0

Gluing turns the 4 corners into V = 1, the 4 sides into E = 2 and the square into F = 1: V − E + F = 1 − 2 + 1 = 0

Same way round as at the start.

Sides crossed so far: 0

Opposite sides are glued with arrows pointing the same way. Walk off the right side and you come back on the left, the same way round. However it walks, the flatlander is never mirrored: the torus is orientable.

Exercises

Exercises

1. Smooth manifolds

Warm-upHuman

What is the dimension of the orthogonal group O(5)O(5) as a smooth manifold?

2. Smooth manifolds

CoreHuman

Which of these subsets are smooth embedded submanifolds of the ambient Euclidean space? (Choose all that apply.)

Choose all that apply

3. Smooth manifolds

CoreHuman

To show that O(n)O(n) is a manifold we used F(A)=ATAF(A) = A^{\mathsf T}A with codomain Symn\mathrm{Sym}_n, the symmetric matrices. What goes wrong if we take the codomain to be all of Mn(R)M_n(\R)?

Choose one

4. Smooth manifolds

CoreHuman

Coordinates change from xx to x~\tilde x near pp, with Jacobian matrix J=∂x~/∂xJ = \partial \tilde x/\partial x at pp. A tangent vector has components vv in the old chart and JvJv in the new one. How do the components of a covector ω\omega (for example dfdf) change?

Choose one

5. Smooth manifolds

CoreHuman

In the proof that dim⁡TpM=n\dim T_pM = n we wrote f(x)=f(a)+∑i(xi−ai)gi(x)f(x) = f(a) + \sum_i (x^i - a^i) g_i(x) and applied a derivation vv at aa. Why does the term (xi−ai)(a)⋅v(gi)(x^i - a^i)(a) \cdot v(g_i) vanish?

Choose one

6. Smooth manifolds

CoreHuman

Let R=[0,2]×[0,3]R = [0, 2] \times [0, 3] with the anticlockwise boundary. Use Stokes's theorem to compute ∮∂Rx2y dx+(x33+2x)dy\displaystyle\oint_{\partial R} x^2 y\, dx + \Big(\frac{x^3}{3} + 2x\Big) dy.

7. Smooth manifolds

CoreHuman

What is the dimension of HdR1(R2∖{p1,p2,p3})H^1_{\mathrm{dR}}(\R^2 \setminus \{p_1, p_2, p_3\}), the first de Rham cohomology of the plane with three points removed?

8. Smooth manifolds

StretchFormal

Which Lean statement correctly formalises "dθd\theta is not exact on the circle R/2πZ\R/2\pi\Z"?

Choose one

9. Smooth manifolds

StretchComputational

Write period(P, Q, n) that approximates the integral of the 11-form P(x,y) dx+Q(x,y) dyP(x, y)\,dx + Q(x, y)\,dy round the unit circle (cos⁡t,sin⁡t)(\cos t, \sin t), 0≤t≤2π0 \le t \le 2\pi, anticlockwise, by the rectangle rule with n equally spaced values of tt. Then check: the period of dθ=(x dy−y dx)/(x2+y2)d\theta = (x\,dy - y\,dx)/(x^2 + y^2) is 2π2\pi, and the period of the exact form d(x2y)d(x^2 y) is 00.

Runs entirely in your browser, in an isolated worker. Nothing is sent to a server.
Common misconceptionsH

Common misconceptions

Misconception

"A manifold comes with a preferred embedding in some RN\R^N." A manifold is an abstract space with an atlas. Some arise as subsets (Sn⊂Rn+1S^n \subset \R^{n+1}), many do not (RP2\R P^2, the configuration space SO(3)SO(3), spacetime). Whitney's theorem says embeddings exist, but there are many and none is canonical, and intrinsic notions (tangent vectors, forms, cohomology) do not depend on one. Even extra structure can resist embedding: the flat torus R2/Z2\R^2/\Z^2 embeds smoothly in R3\R^3 as a doughnut, but no C2C^2 embedding in R3\R^3 is isometric, because every compact surface in R3\R^3 has a point of positive Gaussian curvature.

Misconception

"Vectors and covectors are the same thing; just write the components in a column." They transform differently under a change of chart: v~=Jv\tilde v = Jv but ω~=(JT)−1ω\tilde\omega = (J^{\mathsf T})^{-1}\omega. Only for orthogonal JJ do the laws coincide, which is why the confusion survives in Euclidean geometry with Cartesian coordinates. Example: in polar coordinates df=fr dr+fθ dθdf = f_r\, dr + f_\theta\, d\theta, but the gradient vector is fr ∂r+r−2fθ ∂θf_r\,\partial_r + r^{-2} f_\theta\, \partial_\theta; converting needs a metric.

Misconception

"A closed form is exact." Only locally (the Poincaré lemma), or on spaces with vanishing cohomology. On R2∖{0}\R^2 \setminus \{0\} the form dθd\theta is closed with non-zero period 2π2\pi, so it is not exact. The obstruction is topological, and measuring it is the point of de Rham cohomology.

Misconception

"A homeomorphism between smooth manifolds is a diffeomorphism," or "the topology determines the smooth structure." The homeomorphism x↦x3x \mapsto x^3 of R\R is smooth but its inverse is not. More seriously, there are manifolds homeomorphic but not diffeomorphic to S7S^7 (Milnor) and to R4\R^4. In dimensions 11, 22 and 33, however, smooth structures exist and are unique up to diffeomorphism.

Misconception

"Hausdorff and second countable are technicalities." Drop Hausdorff and the line with two origins appears, where a sequence has two limits and integral curves of vector fields are not unique. Drop second countability and the long line appears: locally like R\R, Hausdorff, but with no partitions of unity, no Riemannian metric and no embedding in any RN\R^N.

Misconception

"Tangent vectors at different points can be compared or added." TpMT_pM and TqMT_qM are different vector spaces; there is no canonical way to say that a vector at pp "equals" one at qq. Comparing them requires extra structure — a connection and its parallel transport — which is where Riemannian geometry begins.

ViewpointsH

Many ways of seeing

Topology— Derivatives that count holes

de Rham's theorem says that a calculus construction — closed forms modulo exact forms — computes singular cohomology with real coefficients, a topological invariant. Periods of closed forms detect holes: HdR1(R2∖{p1,…,pm})≅RmH^1_{\mathrm{dR}}(\R^2 \setminus \{p_1, \dots, p_m\}) \cong \R^m, HdR∙(Tn)H^\bullet_{\mathrm{dR}}(T^n) is an exterior algebra on nn generators, and on a compact oriented manifold integration identifies HdRnH^n_{\mathrm{dR}} with R\R. What de Rham cohomology cannot see is torsion: RP2\R P^2 has H1(RP2;Z)=0H^1(\R P^2; \Z) = 0 but H2(RP2;Z)≅Z/2H^2(\R P^2; \Z) \cong \Z/2, while all its de Rham groups above degree 00 vanish.

The topologist also asks how far "smooth" is from "topological". In dimensions ≤3\le 3 there is no difference; in dimension 44 there are topological manifolds with no smooth structure at all (Freedman's E8E_8 manifold) and exotic R4\R^4s; in dimensions ≥5\ge 5 surgery theory counts smooth structures, as in Kervaire and Milnor's computation that there are 2828 oriented smooth structures on S7S^7 up to orientation-preserving diffeomorphism. See differential topology and low-dimensional topology.

Geometry— From Riemann's manifolds to metrics and symmetry

Riemann introduced manifolds in order to do geometry in them. A Riemannian metric — an inner product on each tangent space, varying smoothly — exists on every smooth manifold (glue the Euclidean metrics of charts with a partition of unity), and it supplies what bare manifolds lack: lengths, angles, volumes, a canonical isomorphism TpM≅Tp∗MT_pM \cong T_p^*M (so gradients exist), geodesics and curvature. Klein's programme meets Riemann's in homogeneous spaces such as Sn=O(n+1)/O(n)S^n = O(n+1)/O(n) and hyperbolic space, where a Lie group acts transitively by isometries. See Riemannian geometry and Lie groups.

Analysis— The inverse function theorem and elliptic PDE

The analyst sees two engines under the hood. The inverse function theorem produces submanifolds (the regular value theorem), local normal forms of maps and, through the ODE existence theorem, flows. Sard's theorem says almost every value of a smooth map is regular, so regular values are plentiful. And analysis returns at the end: on a compact oriented Riemannian manifold, the Hodge theorem says that each de Rham class contains exactly one harmonic form, Δω=(dd∗+d∗d)ω=0\Delta\omega = (dd^* + d^*d)\omega = 0. So dim⁡HdRk(M)\dim H^k_{\mathrm{dR}}(M), a topological number, equals the dimension of the space of solutions of an elliptic PDE — the prototype of index theory (Atiyah–Singer).

Algebra— The manifold is its algebra of functions

For a smooth manifold MM, the points can be recovered from the R\R-algebra C∞(M)C^\infty(M): every non-zero algebra homomorphism C∞(M)→RC^\infty(M) \to \R is evaluation at a unique point ("Milnor's exercise"). Tangent vectors at pp are derivations of this algebra at pp — our definition — and vector fields are derivations of C∞(M)C^\infty(M), which form a Lie algebra under the commutator. This is the starting point of the algebraic geometer's view (a variety is known by its ring of functions) and of noncommutative geometry, which replaces C∞(M)C^\infty(M) by a noncommutative algebra. The de Rham complex (Ω∙(M),∧,d)(\Omega^\bullet(M), \wedge, d) is a commutative differential graded algebra; Sullivan showed that, for simply connected manifolds, it determines the real homotopy type up to quasi-isomorphism, and a version with polynomial forms captures rational homotopy theory.

Physics— Configuration spaces, phase spaces and fields

Classical mechanics happens on manifolds: a system with configuration space QQ (a torus for a double pendulum, SO(3)SO(3) for a spinning top) has phase space the cotangent bundle T∗QT^*Q, which carries the canonical symplectic form ∑idpi∧dqi\sum_i dp_i \wedge dq^i; momenta are covectors, velocities are vectors, and the distinction of this lesson is the distinction between the Hamiltonian and Lagrangian pictures. In electromagnetism the field strength is a 22-form FF on spacetime and Maxwell's equations read dF=0dF = 0, d⋆F=Jd{\star}F = J. A magnetic monopole would make FF closed but not exact on R3∖{0}\R^3 \setminus \{0\}, with flux through a sphere given by a class in HdR2H^2_{\mathrm{dR}}; Dirac's argument that this flux is quantised if electric charge is (1931) is a physical heuristic that became the mathematics of line bundles and characteristic classes. See symplectic geometry and relativity.

Category theory— The tangent functor

Smooth manifolds and smooth maps form a category Man\cat{Man}. The tangent construction is a functor T:Man→ManT : \cat{Man} \to \cat{Man} (indeed to vector bundles): M↦TMM \mapsto TM, F↦dFF \mapsto dF, and the chain rule d(G∘F)=dG∘dFd(G \circ F) = dG \circ dF together with d(id)=idd(\mathrm{id}) = \mathrm{id} is exactly functoriality. de Rham cohomology is a contravariant functor from Man\cat{Man} to graded-commutative algebras, homotopy invariant, and U↦Ωk(U)U \mapsto \Omega^k(U) is a sheaf. See functors and sheaves.

History & global context

History and global context

  • Intrinsic geometry. Carl Friedrich Gauss's Disquisitiones generales circa superficies curvas (1827), written after his geodetic survey of Hanover, showed that the curvature of a surface can be measured from inside it.
  • Riemann's manifolds. Bernhard Riemann's habilitation lecture Über die Hypothesen, welche der Geometrie zu Grunde liegen (1854, published 1868) introduced nn-dimensional "manifoldnesses" with a metric given by a quadratic form on infinitesimal displacements. His thesis (1851) had already introduced the surfaces now named after him.
  • Topology and forms. Henri Poincaré's Analysis Situs (1895) and its supplements studied manifolds through homology and Betti numbers. Élie Cartan introduced the exterior calculus of differential forms (from 1899). The divergence theorem goes back to Mikhail Ostrogradsky (1826) and Gauss, Green's theorem to George Green (1828), and the classical Stokes theorem to a letter of William Thomson (Kelvin) to George Stokes (1850), which Stokes set as a Cambridge examination question in 1854; higher-dimensional versions were stated by Vito Volterra and Poincaré, and the modern form ∫Mdω=∫∂Mω\int_M d\omega = \int_{\partial M}\omega took shape in the Cartan school.
  • The modern definition. Hermann Weyl's Die Idee der Riemannschen Fläche (1913) defined Riemann surfaces by local coordinate neighbourhoods; Oswald Veblen and J. H. C. Whitehead gave the atlas definition of smooth manifolds (1931–32); Hassler Whitney's Differentiable manifolds (1936) proved that abstract manifolds embed in Euclidean space. Georges de Rham proved his theorem in his 1931 thesis.
  • The global theory. In the 1940s Shiing-Shen Chern's intrinsic proof of the Gauss–Bonnet theorem and his characteristic classes, and W. V. D. Hodge's harmonic integrals (developed further by Kunihiko Kodaira and others), made manifolds central to geometry. In 1956 John Milnor found exotic 77-spheres. Vladimir Arnold's Mathematical Methods of Classical Mechanics (1974) made manifolds and forms the standard language of mechanics.
Further reading

Further reading

  • John M. Lee, Introduction to Smooth Manifolds (2nd edition, 2013) — the standard graduate text; Chapters 1–19 cover this lesson in full detail.
  • Loring W. Tu, An Introduction to Manifolds (2nd edition, 2011) — a gentler first course, through Stokes's theorem and de Rham cohomology.
  • Michael Spivak, Calculus on Manifolds (1965) — a short classic on forms and Stokes's theorem in Rn\R^n; and A Comprehensive Introduction to Differential Geometry, Volume 1 (3rd edition, 1999).
  • Raoul Bott and Loring W. Tu, Differential Forms in Algebraic Topology (1982) — de Rham cohomology, Mayer–Vietoris and beyond.
  • Frank W. Warner, Foundations of Differentiable Manifolds and Lie Groups (1971; Springer reprint 1983) — concise, including a proof of de Rham's theorem and the Hodge theorem.
  • V. I. Arnold, Mathematical Methods of Classical Mechanics (2nd edition, 1989) — manifolds and forms in mechanics.