Smooth Manifolds

Six core ideas from a first course on manifolds, each taught in the format that fits it best. Drag points, play an animation, work through an example one step at a time, then check yourself with concept questions and practice problems.

Follows the chapter order of J. M. Lee, Introduction to Smooth Manifolds, 2nd ed. (Springer, 2013). Chapter numbers refer to that book; the questions and examples here are written for this page.

InteractiveLee, Ch. 1: Smooth Manifolds

Charts and transition maps

A manifold is a space that looks like ℝn\R^n up close. Charts are the local pictures. Whether calculus makes sense depends on one thing only: how two overlapping charts translate between each other.

The idea

Topological manifold. A space MM that is Hausdorff, second countable, and locally Euclidean of dimension nn: every point has a neighbourhood UU and a homeomorphism φ:U→U^⊆ℝn\varphi: U \to \hat U \subseteq \R^n onto an open set. The pair (U,φ)(U,\varphi) is a chart.

Two charts (U,φ)(U,\varphi) and (V,ψ)(V,\psi) overlap on U∩VU\cap V. Going from one chart's coordinates to the other's is the transition map

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

a map between open subsets of ℝn\R^n, where we already know what "smooth" means.

Smooth structure. Charts are smoothly compatible if every transition map is smooth with a smooth inverse. A smooth atlas is a collection of compatible charts covering MM; a maximal one is a smooth structure.

Below, the circle S1S^1 gets two atlases. Drag the point and watch both coordinates and the transition map between them.

Explore

Atlas on the circle
Drag the point p around the circle.
The transition map from the first chart's coordinate to the second's.

Worked example

Stereographic projection from the north pole N=(0,…,0,1)N=(0,\dots,0,1) of Sn⊂ℝn+1S^n\subset\R^{n+1} sends xx to the point where the line from NN through xx meets ℝn×{0}\R^n\times\{0\}:

σ(x)=(x1,…,xn)1−xn+1,σ~(x)=(x1,…,xn)1+xn+1 (from the south pole S).\sigma(x) = \frac{(x^1,\dots,x^n)}{1-x^{n+1}}, \qquad \tilde\sigma(x) = \frac{(x^1,\dots,x^n)}{1+x^{n+1}} \ \text{(from the south pole } S).

Show these two charts form a smooth atlas for SnS^n.

  1. What is the overlap of the two chart domains, and what does σ\sigma send it to?

    The domains are Sn∖{N}S^n\smallsetminus\{N\} and Sn∖{S}S^n\smallsetminus\{S\}, so the overlap is Sn∖{N,S}S^n\smallsetminus\{N,S\}. Since σ(S)=0\sigma(S)=0, its image under σ\sigma is ℝn∖{0}\R^n\smallsetminus\{0\}.

  2. Find the inverse σ−1(u)\sigma^{-1}(u) for u∈ℝnu\in\R^n.

    Points on the line from NN through (u,0)(u,0) are (su,1−s)(su,\,1-s). Requiring s2|u|2+(1−s)2=1s^2|u|^2+(1-s)^2=1 gives s=2/(|u|2+1)s = 2/(|u|^2+1), so

    σ−1(u)=(2u, |u|2−1)|u|2+1.\sigma^{-1}(u) = \frac{\big(2u,\ |u|^2-1\big)}{|u|^2+1}.
  3. Compute the transition map σ~∘σ−1\tilde\sigma\circ\sigma^{-1}.

    With x=σ−1(u)x=\sigma^{-1}(u) we have 1+xn+1=2|u|2|u|2+11+x^{n+1} = \dfrac{2|u|^2}{|u|^2+1}, so

    σ~∘σ−1(u)=2u|u|2+1⋅|u|2+12|u|2=u|u|2,u≠0.\tilde\sigma\circ\sigma^{-1}(u) = \frac{2u}{|u|^2+1}\cdot\frac{|u|^2+1}{2|u|^2} = \frac{u}{|u|^2}, \qquad u\neq 0.
  4. Why does this finish the proof?

    u↦u/|u|2u\mapsto u/|u|^2 is smooth on ℝn∖{0}\R^n\smallsetminus\{0\} and is its own inverse, so it is a diffeomorphism. By symmetry the same holds for σ∘σ~−1\sigma\circ\tilde\sigma^{-1}. The two charts are smoothly compatible and cover SnS^n, so they form a smooth atlas. (For n=1n=1 this is the map u↦1/uu\mapsto 1/u graphed above.)

Check your understanding

The "line with two origins" is two copies of ℝ\R glued along ℝ∖{0}\R\smallsetminus\{0\}. It is locally Euclidean and second countable. What goes wrong?

  1. It is not Hausdorff: the two origins cannot be separated by disjoint open sets.
  2. It is not locally Euclidean at the origins.Each origin has a neighbourhood that is a copy of an open interval, so it is locally Euclidean everywhere.
  3. Nothing; it is a 1-manifold.Look at any two neighbourhoods of the two origins: they always meet.
Every neighbourhood of one origin contains points x≠0x\neq 0 near 0, and so does every neighbourhood of the other, so they always intersect. That is exactly why the Hausdorff condition is part of the definition: it makes limits unique.

On ℝ\R, take the charts φ(x)=x\varphi(x)=x and ψ(x)=x3\psi(x)=x^3. Are they smoothly compatible?

  1. Yes: both are homeomorphisms and x3x^3 is smooth.Compatibility needs both transition maps to be smooth, not just one.
  2. No: φ∘ψ−1(y)=y1/3\varphi\circ\psi^{-1}(y)=y^{1/3} is not differentiable at 00.
  3. No: x3x^3 is not a smooth function.x3x^3 is smooth. The problem is its inverse.
ψ∘φ−1(x)=x3\psi\circ\varphi^{-1}(x)=x^3 is smooth, but the other direction is the cube root, which has infinite slope at 0. So ψ\psi defines a different smooth structure on ℝ\R.

Continuing: write ℝψ\R_\psi for ℝ\R with the smooth structure defined by ψ(x)=x3\psi(x)=x^3. Is ℝψ\R_\psi diffeomorphic to the usual ℝ\R?

  1. Yes: F(x)=x3F(x)=x^3 is a diffeomorphism ℝψ→ℝ\R_\psi\to\R.
  2. No: different smooth structures are never diffeomorphic.Different structures can still be diffeomorphic by a map other than the identity.
  3. No: the identity map ℝψ→ℝ\R_\psi\to\R is not smooth.That's true, but it only rules out one particular map.
In the charts ψ\psi and id\mathrm{id}, FF reads F∘ψ−1(y)=(y1/3)3=yF\circ\psi^{-1}(y) = (y^{1/3})^3 = y, the identity, so FF is a diffeomorphism. "Different structure" and "different manifold up to diffeomorphism" are different questions.

What is the smallest number of charts in an atlas for S1S^1?

  1. 1A single chart would make S1S^1 homeomorphic to an open subset of ℝ\R. S1S^1 is compact; no nonempty open subset of ℝ\R is.
  2. 2
  3. 3The stereographic atlas above already does it with fewer.
  4. 4
One chart is impossible because S1S^1 is compact, and two stereographic charts suffice.

In the angle atlas above, the transition map is θ↦θ\theta\mapsto\theta on one part of the overlap and θ↦θ+2π\theta\mapsto\theta+2\pi on the other. Is it smooth?

  1. Yes. The overlap has two components, and on each the map is smooth.
  2. No, it jumps by 2π2\pi.The jump happens between two separate pieces of the domain, never inside one.
Smoothness is a local property. The domain (−π,0)∪(0,π)(-\pi,0)\cup(0,\pi) is disconnected, and on each piece the map is a translation.

Practice

Real projective lineCharts and a transition map

ℝℙ1\RP^1 is the set of lines through 0 in ℝ2\R^2, written [x:y][x:y]. Use the charts u=y/xu=y/x on {x≠0}\{x\neq0\} and v=x/yv=x/y on {y≠0}\{y\neq0\}. Compute the transition map and conclude the charts form a smooth atlas.

Hint

On the overlap both xx and yy are nonzero. Express vv in terms of uu.

Solution

On the overlap u≠0u\neq0 and v=x/y=1/uv = x/y = 1/u. The transition u↦1/uu\mapsto 1/u on ℝ∖{0}\R\smallsetminus\{0\} is smooth and its own inverse. Notice it is the same transition map as stereographic projection of S1S^1. That is a hint that ℝℙ1\RP^1 and S1S^1 are diffeomorphic.

Graphs are manifoldsA single chart

Let f:U→ℝkf:U\to\R^k be continuous on an open U⊆ℝnU\subseteq\R^n. Show its graph Γ(f)={(x,f(x))}⊆ℝn+k\Gamma(f)=\{(x,f(x))\}\subseteq\R^{n+k} is a topological nn-manifold covered by one chart.

Hint

Try the projection (x,y)↦x(x,y)\mapsto x.

Solution

π(x,y)=x\pi(x,y)=x restricted to Γ(f)\Gamma(f) is continuous with continuous inverse x↦(x,f(x))x\mapsto(x,f(x)), so it is a homeomorphism onto UU. Hausdorff and second countable are inherited from ℝn+k\R^{n+k}. With one chart there are no transition maps, so this chart also defines a smooth structure, even when ff is not smooth.

A torus atlasProducts of charts

Using the stereographic charts on S1S^1, write down a smooth atlas on T2=S1×S1T^2=S^1\times S^1. How many charts did you use, and what are the transition maps?

Solution

Take all products σa×σb\sigma_a\times\sigma_b with a,b∈{N,S}a,b\in\{N,S\}: four charts. Each transition map acts coordinate by coordinate as either u↦uu\mapsto u or u↦1/uu\mapsto1/u, so all of them are smooth.

InteractiveLee, Ch. 3: Tangent Vectors

Tangent vectors and the differential

On a manifold there is no ambient space for arrows to live in, so a tangent vector is defined by what it does: it differentiates functions. A smooth map pushes these vectors forward, and in coordinates that push-forward is just the Jacobian matrix.

The idea

Tangent vector. A derivation at pp: a linear map v:C∞(M)→ℝv:C^\infty(M)\to\R with v(fg)=f(p)vg+g(p)vfv(fg)=f(p)\,vg+g(p)\,vf. These form the tangent space TpMT_pM, which has dimension nn and basis ∂∂x1|p,…,∂∂xn|p\left.\pd{}{x^1}\right|_p,\dots,\left.\pd{}{x^n}\right|_p in any chart.
Differential. For smooth F:M→NF:M\to N,  dFp:TpM→TF(p)N\ dF_p:T_pM\to T_{F(p)}N is dFp(v)f=v(f∘F).dF_p(v)f = v(f\circ F). In coordinates it is the Jacobian:
dFp(∂∂xi|p)=∑j∂Fj∂xi(p)∂∂yj|F(p).dF_p\!\left(\left.\pd{}{x^i}\right|_p\right) = \sum_j \pd{F^j}{x^i}(p)\left.\pd{}{y^j}\right|_{F(p)}.
Velocities. If γ\gamma is a curve with γ(0)=p\gamma(0)=p and γ′(0)=v\gamma'(0)=v, then (F∘γ)′(0)=dFp(v)(F\circ\gamma)'(0)=dF_p(v). This is the chain rule, and it is what you can see below.

Explore

The map F
Domain. Drag the point p and the tip of the vector v.
dFp(v)dF_p(v) image of the curve p+tvp+tv

Worked example

On ℝ2∖{0}\R^2\smallsetminus\{0\} use polar coordinates, x=rcos⁡θ, y=rsin⁡θx=r\cos\theta,\ y=r\sin\theta. Express the polar coordinate vectors in Cartesian terms, and a Cartesian vector in polar terms.

  1. Write ∂∂r\pd{}{r} in terms of ∂∂x\pd{}{x} and ∂∂y\pd{}{y}.

    The chain rule, which is the change-of-coordinates formula for tangent vectors, gives

    ∂∂r=∂x∂r∂∂x+∂y∂r∂∂y=cos⁡θ∂∂x+sin⁡θ∂∂y=x∂x+y∂yr.\pd{}{r} = \pd{x}{r}\pd{}{x}+\pd{y}{r}\pd{}{y} = \cos\theta\,\pd{}{x}+\sin\theta\,\pd{}{y} = \frac{x\,\partial_x + y\,\partial_y}{r}.
  2. Now do ∂∂θ\pd{}{\theta}.

    ∂∂θ=−rsin⁡θ∂∂x+rcos⁡θ∂∂y=−y∂∂x+x∂∂y.\pd{}{\theta} = -r\sin\theta\,\pd{}{x}+r\cos\theta\,\pd{}{y} = -y\,\pd{}{x}+x\,\pd{}{y}.

    This is the rotation field. Note it grows with rr: the coordinate vector ∂θ\partial_\theta is not a unit vector.

  3. Sanity check: apply both to f=x2+y2f=x^2+y^2.

    ∂rf=(x⋅2x+y⋅2y)/r=2r\partial_r f = (x\cdot 2x + y\cdot 2y)/r = 2r and ∂θf=−y⋅2x+x⋅2y=0\partial_\theta f = -y\cdot 2x + x\cdot 2y = 0. That matches f=r2f=r^2, which doesn't depend on θ\theta.

  4. Write v=3∂xv = 3\,\partial_x at p=(1,1)p=(1,1) in the polar basis.

    Go the other way with the inverse Jacobian: ∂x=∂r∂x∂r+∂θ∂x∂θ=xr∂r−yr2∂θ\partial_x = \pd{r}{x}\partial_r+\pd{\theta}{x}\partial_\theta = \tfrac{x}{r}\partial_r - \tfrac{y}{r^2}\partial_\theta. At (1,1)(1,1), r=2r=\sqrt2, so

    v=32∂r−32∂θ.v = \frac{3}{\sqrt2}\,\partial_r - \frac32\,\partial_\theta.

    Check: ∂r=(∂x+∂y)/2\partial_r=(\partial_x+\partial_y)/\sqrt2 and ∂θ=−∂x+∂y\partial_\theta=-\partial_x+\partial_y there, and 32⋅12(1,1)−32(−1,1)=(3,0)\tfrac{3}{\sqrt2}\cdot\tfrac{1}{\sqrt2}(1,1)-\tfrac32(-1,1)=(3,0).

Check your understanding

Let v∈TpMv\in T_pM be a derivation and cc a constant function. What is v(c)v(c)?

  1. 00
  2. ccTry the product rule on 1⋅11\cdot1.
  3. It depends on vv.The product rule pins it down for every derivation.
v(1)=v(1⋅1)=1⋅v(1)+1⋅v(1)=2v(1)v(1)=v(1\cdot1)=1\cdot v(1)+1\cdot v(1)=2v(1), so v(1)=0v(1)=0, and by linearity v(c)=cv(1)=0v(c)=c\,v(1)=0.

For smooth F:M→NF:M\to N, v∈TpMv\in T_pM and f∈C∞(N)f\in C^\infty(N), what is dFp(v)fdF_p(v)\,f?

  1. v(f∘F)v(f\circ F)
  2. F(vf)F(vf)ff is a function on NN, so vv can't act on it directly.
  3. (vf)∘F(vf)\circ Fvv only differentiates functions on MM.
Pull ff back to MM by composing with FF, then let vv differentiate it. Everything else in the theory of differentials follows from this formula.

Polar coordinates F(r,θ)=(rcos⁡θ,rsin⁡θ)F(r,\theta)=(r\cos\theta,r\sin\theta): what is the rank of dFdF at a point with r=0r=0?

  1. 2Set r=0r=0 in the Jacobian. One column vanishes.
  2. 1
  3. 0The ∂r\partial_r column is (cos⁡θ,sin⁡θ)(\cos\theta,\sin\theta), which is never zero.
The Jacobian has columns (cos⁡θ,sin⁡θ)(\cos\theta,\sin\theta) and (−rsin⁡θ,rcos⁡θ)(-r\sin\theta,r\cos\theta). At r=0r=0 the second is zero. Moving in θ\theta doesn't move the image point at all, as you can check above.

MM is an nn-manifold. What is the dimension of its tangent bundle TMTM?

  1. nn
  2. 2n2n
  3. n2n^2
A point of TMTM is a base point (nn coordinates) together with a tangent vector there (nn components).

Practice

A paraboloid mapComputing a differential

Let F(u,v)=(u,v,u2+v2)F(u,v)=(u,\,v,\,u^2+v^2). Compute dF(1,2)(∂u)dF_{(1,2)}(\partial_u) and dF(1,2)(∂v)dF_{(1,2)}(\partial_v), and describe the image of dF(1,2)dF_{(1,2)}.

Solution

The Jacobian columns are (1,0,2u)(1,0,2u) and (0,1,2v)(0,1,2v). At (1,2)(1,2): dF(∂u)=∂x+2∂zdF(\partial_u)=\partial_x+2\partial_z and dF(∂v)=∂y+4∂zdF(\partial_v)=\partial_y+4\partial_z. The image is the plane spanned by these, which is the tangent plane to the paraboloid at (1,2,5)(1,2,5). FF is an immersion.

Derivations are localBump functions

Show that if f=gf=g on a neighbourhood of pp, then vf=vgvf=vg for every v∈TpMv\in T_pM.

Hint

Take a bump function ψ\psi that is 11 near pp and supported where f=gf=g. Look at ψ(f−g)\psi\,(f-g).

Solution

Choose ψ\psi supported in the neighbourhood where f=gf=g, with ψ(p)=1\psi(p)=1. Then ψ(f−g)≡0\psi(f-g)\equiv0, so

0=v(ψ(f−g))=ψ(p)v(f−g)+(f−g)(p)vψ=v(f−g).0 = v(\psi(f-g)) = \psi(p)\,v(f-g) + (f-g)(p)\,v\psi = v(f-g).
This is what lets derivations be computed in a chart.

The coordinate vectors are independentUsing coordinate functions

In a chart (x1,…,xn)(x^1,\dots,x^n), show that ∂∂x1|p,…,∂∂xn|p\left.\pd{}{x^1}\right|_p,\dots,\left.\pd{}{x^n}\right|_p are linearly independent.

Solution

If ∑iai∂i|p=0\sum_i a^i\partial_i|_p=0, apply it to the coordinate function xjx^j (extended to a global function by a bump function). This gives ∑iaiδij=aj=0\sum_i a^i\delta_i^j=a^j=0 for each jj.

Interactive galleryLee, Ch. 4: Submersions, Immersions, and Embeddings

Immersions and embeddings

"Smooth", "immersion", "injective immersion" and "embedding" are four different conditions, and the gap between each pair is easiest to see in a picture. Here is a gallery of curves, one for each gap.

The idea

Rank. The rank of FF at pp is the rank of dFpdF_p. FF is an immersion if dFpdF_p is injective everywhere (rank =dim⁡M=\dim M), and a submersion if it is surjective everywhere (rank =dim⁡N=\dim N).
Embedding. An immersion that is also a homeomorphism onto its image with the subspace topology.
Rank theorem. If FF has constant rank rr near pp, there are charts in which F(x1,…,xm)=(x1,…,xr,0,…,0)F(x^1,\dots,x^m) = (x^1,\dots,x^r,0,\dots,0).
Watch out. A smooth curve can have a corner (when its velocity vanishes), and an injective immersion can still fail to be an embedding.

Explore

velocity γ′(t)\gamma'(t) the part traced so far
Curve

Worked example

Show that β(t)=(sin⁡2t, sin⁡t)\beta(t)=(\sin2t,\ \sin t) for t∈(−π,π)t\in(-\pi,\pi) (the figure eight in the gallery) is an injective immersion but not an embedding.

  1. Show β\beta is an immersion.

    β′(t)=(2cos⁡2t, cos⁡t)\beta'(t)=(2\cos2t,\ \cos t). For this to vanish we need cos⁡t=0\cos t=0, so t=±π/2t=\pm\pi/2. But then cos⁡2t=cos⁡(±π)=−1≠0\cos2t=\cos(\pm\pi)=-1\neq0. So β′(t)≠0\beta'(t)\neq0 everywhere, and for a curve that means rank 1.

  2. Show β\beta is injective.

    Suppose β(t)=β(s)\beta(t)=\beta(s). If sin⁡t=sin⁡s≠0\sin t=\sin s\neq0, then sin⁡2t=2sin⁡tcos⁡t\sin2t=2\sin t\cos t forces cos⁡t=cos⁡s\cos t=\cos s, and sine and cosine together determine the angle, so t=st=s. If sin⁡t=0\sin t=0, the only such tt in (−π,π)(-\pi,\pi) is 00.

  3. Find a sequence that shows β−1\beta^{-1} (defined on the image) is not continuous.

    Take tk→πt_k\to\pi. Then β(tk)→(0,0)=β(0)\beta(t_k)\to(0,0)=\beta(0), but tk↛0t_k\not\to0. So β\beta is not a homeomorphism onto its image.

    Another way to see it: the image is the whole figure eight, which is compact, while (−π,π)(-\pi,\pi) is not. A homeomorphism would preserve compactness.

  4. Is the image an embedded submanifold of ℝ2\R^2 at all?

    No. Removing the crossing point from a small neighbourhood of it in the image leaves four pieces, but removing a point from an interval leaves two. So no neighbourhood of the crossing is homeomorphic to an interval.

Check your understanding

γ(t)=(t3,t3)\gamma(t)=(t^3,t^3). Its image is the line y=xy=x, a perfectly good embedded submanifold. Is γ\gamma an immersion?

  1. Yes, because its image is a smooth line.Being an immersion is a property of the map, not of its image.
  2. No, because γ′(0)=0\gamma'(0)=0.
γ′(t)=(3t2,3t2)\gamma'(t)=(3t^2,3t^2) vanishes at t=0t=0. The map pauses there, even though the image looks perfectly smooth. Try it in the gallery and watch the velocity arrow shrink.

If F:Mm→NnF:M^m\to N^n is an immersion, which must hold?

  1. m≤nm\le n
  2. m≥nm\ge nThat's the condition for a submersion.
  3. m=nm=n
An injective linear map TpM→TF(p)NT_pM\to T_{F(p)}N can only exist if dim⁡TpM≤dim⁡TF(p)N\dim T_pM\le\dim T_{F(p)}N.

MM is compact and F:M→NF:M\to N is an injective immersion. Is FF an embedding?

  1. Yes, always.
  2. Not necessarily: the figure eight is a counterexample.The figure eight is parametrised by (−π,π)(-\pi,\pi), which is not compact.
A continuous bijection from a compact space onto a Hausdorff space is a homeomorphism. Closed subsets of MM are compact, their images are compact, hence closed.

The line of irrational slope on the torus T2T^2 is an injective immersion of ℝ\R whose image is dense. Is its image an embedded submanifold?

  1. Yes, it's the image of an injective immersion.That makes it an immersed submanifold. Embedded is stronger.
  2. No. Every small open set of the torus meets it in infinitely many separate segments.
An embedded submanifold looks like a single slice ℝk×{0}\R^k\times\{0\} in some chart around each of its points. Here every neighbourhood contains infinitely many parallel pieces of the line.

The cusp γ(t)=(t2,t3)\gamma(t)=(t^2,t^3). Which is true?

  1. It is not smooth, because the image has a corner.Both components are polynomials, so the map is smooth.
  2. It is smooth and injective, but not an immersion at t=0t=0.
  3. It is an immersion that is not injective.Different tt give different values of t3t^3, so it is injective.
γ′(0)=(0,0)\gamma'(0)=(0,0). A smooth map can have a sharp image exactly where its derivative vanishes.

Practice

The circle, twiceImmersion vs. embedding

Show F(t)=(cos⁡t,sin⁡t)F(t)=(\cos t,\sin t), F:ℝ→ℝ2F:\R\to\R^2, is an immersion but not injective. Then show the inclusion S1↪ℝ2S^1\hookrightarrow\R^2 is an embedding.

Solution

|F′(t)|=1|F'(t)|=1, so FF is an immersion, and F(t+2π)=F(t)F(t+2\pi)=F(t), so it is not injective. The inclusion of S1S^1 is an injective immersion of a compact manifold, so it is an embedding (see the third concept check).

Where a submersion failsCritical points

Where does f(x,y)=x2−y2f(x,y)=x^2-y^2 fail to be a submersion ℝ2→ℝ\R^2\to\R? What does the rank theorem say near any other point?

Solution

df=2xdx−2ydydf=2x\,dx-2y\,dy is zero only at the origin. Everywhere else the rank is 1, so near those points there are coordinates in which f(u,v)=uf(u,v)=u: the level sets near that point are parallel lines.

InteractiveLee, Ch. 5: Submanifolds

Regular level sets

Most manifolds you meet in practice are carved out by equations: spheres, matrix groups, constraint surfaces. The regular level set theorem says when an equation defines a manifold. You can watch it work, and fail, as you move the level.

The idea

Regular value. c∈Nc\in N is a regular value of F:M→NF:M\to N if dFpdF_p is surjective at every p∈F−1(c)p\in F^{-1}(c). Otherwise cc is a critical value.
Regular level set theorem. If cc is a regular value of F:Mm→NnF:M^m\to N^n, then F−1(c)F^{-1}(c) is a properly embedded submanifold of dimension m−nm-n, and TpF−1(c)=ker⁡dFpT_p F^{-1}(c)=\ker dF_p.
Only one direction. If cc is critical the theorem says nothing. The level set may still be a manifold.

For f:ℝ2→ℝf:\R^2\to\R, "regular" just means ∇f≠0\nabla f\neq0 at every point of the curve f=cf=c. Then the curve is smooth, and its tangent line is perpendicular to ∇f\nabla f.

Explore

the level set f=cf=c ∇f\nabla f along it critical points
Function

Worked example

Show that the orthogonal group O(n)={A:A𝖳A=I}O(n)=\{A: A^{\mathsf T}A=I\} is an embedded submanifold of the space M(n,ℝ)≅ℝn2M(n,\R)\cong\R^{n^2} of all n×nn\times n matrices, and find its dimension.

  1. Pick the defining map. What should its target be?

    Φ(A)=A𝖳A\Phi(A)=A^{\mathsf T}A. Its values are always symmetric, so take the target to be the space of symmetric matrices Sym(n)\mathrm{Sym}(n), of dimension n(n+1)/2n(n+1)/2. (Using all of M(n,ℝ)M(n,\R) as the target would make every value critical.)

  2. Compute dΦA(B)d\Phi_A(B).

    Differentiate along A+sBA+sB:

    dΦA(B)=dds|s=0(A+sB)𝖳(A+sB)=B𝖳A+A𝖳B.d\Phi_A(B)=\frac{d}{ds}\Big|_{s=0}(A+sB)^{\mathsf T}(A+sB) = B^{\mathsf T}A + A^{\mathsf T}B.
  3. Show dΦAd\Phi_A is onto Sym(n)\mathrm{Sym}(n) whenever A∈O(n)A\in O(n).

    Given symmetric CC, try B=12ACB=\tfrac12AC. Then

    dΦA(B)=12C𝖳A𝖳A+12A𝖳AC=12C+12C=C,d\Phi_A(B) = \tfrac12 C^{\mathsf T}A^{\mathsf T}A + \tfrac12 A^{\mathsf T}AC = \tfrac12C+\tfrac12C = C,

    using A𝖳A=IA^{\mathsf T}A=I and C𝖳=CC^{\mathsf T}=C. So II is a regular value.

  4. Conclude: what are the dimension and the tangent space at II?

    dim⁡O(n)=n2−n(n+1)2=n(n−1)2,TIO(n)=ker⁡dΦI={B:B𝖳+B=0}.\dim O(n) = n^2-\frac{n(n+1)}2 = \frac{n(n-1)}2, \qquad T_IO(n)=\ker d\Phi_I = \{B: B^{\mathsf T}+B=0\}.

    The tangent space at the identity is the skew-symmetric matrices. This is the Lie algebra 𝔬(n)\mathfrak{o}(n), which comes back in Chapter 8.

Check your understanding

For f(x,y)=x2−y2f(x,y)=x^2-y^2, is 00 a regular value?

  1. Yes, because ff takes the value 0.Every value in the image could be regular or critical. What matters is ∇f\nabla f on the level set.
  2. No: the origin lies on f=0f=0 and ∇f(0,0)=0\nabla f(0,0)=0.
The level set is two crossing lines. Pick "saddle" in the explorer and move cc through 0 to watch the two branches pinch together.

g(x,y)=(x2+y2−1)2g(x,y)=(x^2+y^2-1)^2. The level set g=0g=0 is the unit circle, yet 00 is a critical value. What follows?

  1. The unit circle is not a submanifold.The theorem only goes one way: a critical value proves nothing.
  2. The theorem doesn't apply to this gg. The circle is still a submanifold, as another defining function shows.
  3. The theorem is wrong here.
∇g=4(x2+y2−1)(x,y)\nabla g = 4(x^2+y^2-1)(x,y) vanishes on the whole circle. But h=x2+y2h=x^2+y^2 has 1 as a regular value with the same level set. Being a submanifold is a property of the set; regularity is a property of the function you chose. The "trap" option in the explorer shows this one.

F:ℝ5→ℝ2F:\R^5\to\R^2 is smooth and cc is a regular value. What is the dimension of F−1(c)F^{-1}(c)?

  1. 2
  2. 3
  3. 5
  4. 7
Codimension equals the dimension of the target: 5−2=35-2=3.

At p∈Sn={|x|2=1}p\in S^n=\{|x|^2=1\}, what is TpSnT_pS^n as a subspace of ℝn+1\R^{n+1}?

  1. {v:v⋅p=0}\{v: v\cdot p=0\}
  2. span(p)\mathrm{span}(p)That's the normal line, the direction of the gradient.
  3. All of ℝn+1\R^{n+1}
TpSn=ker⁡d(|x|2)p={v:2p⋅v=0}T_pS^n=\ker d(|x|^2)_p=\{v: 2p\cdot v=0\}. It is the tangent space you'd expect from calculus.

SL(2,ℝ)=det−1⁡(1)SL(2,\R)=\det^{-1}(1). What is its dimension?

  1. 1
  2. 2
  3. 3
  4. 4
11 is a regular value of det⁡:ℝ4→ℝ\det:\R^4\to\R (see the practice problem below), so dim⁡=4−1=3\dim=4-1=3.

Practice

The special linear groupA regular value of the determinant

Show d(det⁡)A(B)=det⁡(A)tr(A−1B)d(\det)_A(B)=\det(A)\,\mathrm{tr}(A^{-1}B) for invertible AA, and conclude that SL(n,ℝ)SL(n,\R) is an embedded submanifold of dimension n2−1n^2-1.

Hint

Write det⁡(A+sB)=det⁡Adet⁡(I+sA−1B)\det(A+sB)=\det A\,\det(I+sA^{-1}B) and expand det⁡(I+sX)\det(I+sX) to first order in ss.

Solution

det⁡(I+sX)=1+strX+O(s2)\det(I+sX)=1+s\,\mathrm{tr}X+O(s^2), which gives the formula. On SL(n,ℝ)SL(n,\R), taking B=cA/nB=cA/n gives d(det⁡)A(B)=cd(\det)_A(B)=c, so the differential is onto ℝ\R and 11 is a regular value. The dimension is n2−1n^2-1.

Critical values of a cubicMatches the "cubic" option in the explorer

Find the critical points and critical values of f(x,y)=y2−x3+xf(x,y)=y^2-x^3+x. For which cc is f−1(c)f^{-1}(c) guaranteed to be a smooth curve?

Solution

∇f=(1−3x2, 2y)=0\nabla f=(1-3x^2,\ 2y)=0 at (±1/3,0)(\pm1/\sqrt3,\,0), with values f=∓233≈∓0.385f=\mp\tfrac{2}{3\sqrt3}\approx\mp0.385. Every other cc is a regular value. At c≈0.385c\approx0.385 the curve has a self-crossing (a node). At c≈−0.385c\approx-0.385 it has an isolated point next to a smooth branch.

Quadric surfacesVarying the level

For which cc is {x2+y2−z2=c}\{x^2+y^2-z^2=c\} a regular level set? Describe the surfaces for c<0c<0, c=0c=0, c>0c>0.

Solution

∇=2(x,y,−z)\nabla=2(x,y,-z) vanishes only at the origin, which lies on the level c=0c=0. So every c≠0c\neq0 is regular. c>0c>0: a one-sheeted hyperboloid. c<0c<0: a two-sheeted hyperboloid. c=0c=0: a double cone, which fails to be a manifold at its tip.

Animated lessonLee, Ch. 8: Vector Fields · Ch. 9: Integral Curves and Flows

Flows and the Lie bracket

The Lie bracket [X,Y][X,Y] is easy to define algebraically and hard to picture. This short animation shows what it measures: follow XX, then YY, then go back along XX and then YY, and see how far you are from where you started.

The idea

Flow. An integral curve of XX satisfies γ′(t)=Xγ(t)\gamma'(t)=X_{\gamma(t)}. The flow θt\theta_t moves every point along its integral curve for time tt.
Lie bracket. [X,Y]f=X(Yf)−Y(Xf)[X,Y]f = X(Yf)-Y(Xf). Second derivatives cancel, so this is again a vector field. In coordinates:
[X,Y]=∑j(XYj−YXj)∂∂xj.[X,Y] = \sum_j \big(X Y^j - Y X^j\big)\pd{}{x^j}.
What it measures. With θ\theta the flow of XX and ψ\psi the flow of YY,
ψ−t∘θ−t∘ψt∘θt(p)=p+t2[X,Y]p+O(t3).\psi_{-t}\circ\theta_{-t}\circ\psi_{t}\circ\theta_{t}(p) = p + t^2\,[X,Y]_p + O(t^3).

So [X,Y]=0[X,Y]=0 exactly when the flows commute.

Watch

Vector fields
XX YY gap t2[X,Y]pt^2[X,Y]_p · drag p
  1. Two vector fields: XX in blue and YY in gold. Pick a start point pp.
  2. Flow along XX for time tt.
  3. Then along YY for time tt.
  4. Then backwards along XX for time tt.
  5. Then backwards along YY. If the flows commuted, we'd be back at pp.
  6. The loop doesn't close. The gap (red) is close to t2[X,Y]pt^2[X,Y]_p (green), and the match improves as tt shrinks.

Worked example

On ℝ2\R^2, let X=−y∂x+x∂yX=-y\,\partial_x+x\,\partial_y (rotation) and Y=x∂x+y∂yY=x\,\partial_x+y\,\partial_y (scaling). Compute [X,Y][X,Y], interpret the answer using the flows, then compute [X,∂x][X,\partial_x].

  1. Compute XYjXY^j and YXjYX^j for each component.

    Y1=x, Y2=yY^1=x,\ Y^2=y, so XY1=X(x)=−yXY^1=X(x)=-y and XY2=X(y)=xXY^2=X(y)=x.

    X1=−y, X2=xX^1=-y,\ X^2=x, so YX1=Y(−y)=−yYX^1=Y(-y)=-y and YX2=Y(x)=xYX^2=Y(x)=x.

  2. Assemble [X,Y][X,Y].

    [X,Y]=(−y−(−y))∂x+(x−x)∂y=0.[X,Y] = (-y-(-y))\,\partial_x + (x-x)\,\partial_y = 0.
  3. Explain the answer using the flows.

    The flow of XX is rotation by angle tt, and the flow of YY is scaling by ese^{s}. Rotating and scaling commute, Rt(esp)=esRtpR_t(e^{s}p)=e^{s}R_tp, so the four-leg loop always closes. Choose "rotation & scaling" above and check that the gap is zero for every tt.

  4. Now replace YY with ∂x\partial_x. Compute [X,∂x][X,\partial_x].

    The components of ∂x\partial_x are constant, so the XYjXY^j terms vanish. And ∂xX1=∂x(−y)=0\partial_x X^1=\partial_x(-y)=0, ∂xX2=∂x(x)=1\partial_x X^2=\partial_x(x)=1. So

    [X,∂x]=−∂y.[X,\partial_x] = -\partial_y.

    In the animation ("rotation & translation"), the gap points straight down with length about t2t^2, at every starting point.

Check your understanding

X=∂xX=\partial_x, Y=x∂yY=x\,\partial_y. What is [X,Y][X,Y]?

  1. ∂y\partial_y
  2. −∂y-\partial_yCheck the order: XYj−YXjXY^j-YX^j.
  3. 00YY depends on xx, which is the direction XX moves in.
  4. x∂xx\,\partial_x
XY2=∂x(x)=1XY^2=\partial_x(x)=1 and every other term is 0. This is the first pair in the animation. Move pp around and the gap stays the same, because the bracket is constant.

Which identity is false in general?

  1. [X,Y]=−[Y,X][X,Y]=-[Y,X]That one is true (antisymmetry).
  2. [fX,Y]=f[X,Y][fX,Y]=f[X,Y]
  3. [X,[Y,Z]]+[Y,[Z,X]]+[Z,[X,Y]]=0[X,[Y,Z]]+[Y,[Z,X]]+[Z,[X,Y]]=0That's the Jacobi identity, which always holds.
The correct rule is [fX,Y]=f[X,Y]−(Yf)X[fX,Y]=f[X,Y]-(Yf)X. The bracket is not C∞C^\infty-linear, so it is not a tensor. It involves derivatives of the coefficients.

On ℝ\R, what is the integral curve of X=x2∂xX=x^2\,\partial_x starting at x=1x=1?

  1. x(t)=etx(t)=e^tThat solves x′=xx'=x, not x′=x2x'=x^2.
  2. x(t)=1/(1−t)x(t)=1/(1-t), defined only for t<1t<1
  3. x(t)=1+t2x(t)=1+t^2
Solving x′=x2x'=x^2 gives x(t)=1/(1−t)x(t)=1/(1-t), which runs off to infinity in finite time. So XX is not complete: its flow is not defined for all tt.

Which vector fields are guaranteed to be complete?

  1. Compactly supported ones (so every vector field on a compact manifold)
  2. Bounded ones on any manifold∂x\partial_x on the open interval (0,1)(0,1) is bounded but runs off the end.
  3. Polynomial ones on ℝn\R^nThe previous question is a counterexample.
Outside the support the flow is stationary, and on the compact support there is a uniform time for which all integral curves exist. That uniform time lets you extend every curve forever.

Practice

A hyperbolic flowSolving for a flow

Find the flow of X=x∂x−y∂yX=x\,\partial_x-y\,\partial_y on ℝ2\R^2 and sketch its integral curves.

Solution

x′=x, y′=−yx'=x,\ y'=-y gives θt(x,y)=(etx, e−ty)\theta_t(x,y)=(e^tx,\ e^{-t}y). The integral curves are the hyperbolas xy=constxy=\text{const}, the four half-axes, and the fixed point at the origin.

The bracket with functionsA product rule

Show [fX,gY]=fg[X,Y]+f(Xg)Y−g(Yf)X[fX,gY]=fg[X,Y]+f(Xg)Y-g(Yf)X.

Solution

Apply both sides to hh:

fX(gYh)−gY(fXh)=f(Xg)(Yh)+fgXYh−g(Yf)(Xh)−gfYXh.\begin{aligned} fX(gYh)-gY(fXh) &= f(Xg)(Yh)+fg\,XYh\\ &\quad -g(Yf)(Xh)-gf\,YXh. \end{aligned}

The two second-order terms combine to fg[X,Y]hfg[X,Y]h.

A shear flowFlows that don't commute

Find the flows of ∂x\partial_x and x∂yx\,\partial_y. Compute the four-leg loop from (0,0)(0,0) exactly and compare it with t2[X,Y]t^2[X,Y].

Solution

θt(x,y)=(x+t,y)\theta_t(x,y)=(x+t,y) and ψt(x,y)=(x,y+tx)\psi_t(x,y)=(x,y+tx). From (0,0)(0,0): (t,0)→(t,t2)→(0,t2)→(0,t2)(t,0)\to(t,t^2)\to(0,t^2)\to(0,t^2). The gap is exactly t2∂y=t2[X,Y]t^2\partial_y=t^2[X,Y], with no O(t3)O(t^3) error, because these fields are so simple.

Draw it yourselfLee, Ch. 11 · 14 · 16 · 17

Forms, Stokes and de Rham

Integrate a 1-form around a loop you draw. For an exact form the answer is always zero. For the area form it is the enclosed area (Stokes). For the angle form dθd\theta it counts how often you wound around the hole, which is the first de Rham cohomology of the punctured plane.

The idea

Line integral. For a 1-form ω=Pdx+Qdy\omega=P\,dx+Q\,dy and a curve γ\gamma, ∫γω=∫abωγ(t)(γ′(t))dt=∫abγ*ω.\displaystyle\int_\gamma\omega=\int_a^b\omega_{\gamma(t)}(\gamma'(t))\,dt=\int_a^b \gamma^*\omega. It doesn't depend on how the curve is parametrised, only on its direction.
Exterior derivative.
d(Pdx+Qdy)=(∂xQ−∂yP)dx∧dy,d(P\,dx+Q\,dy)=\big(\partial_xQ-\partial_yP\big)\,dx\wedge dy,
and d∘d=0d\circ d=0. ω\omega is closed if dω=0d\omega=0 and exact if ω=df\omega=df.
Stokes. ∫Mdω=∫∂Mω\displaystyle\int_M d\omega=\int_{\partial M}\omega. For exact ω=df\omega=df and a closed loop, ∫γdf=0\int_\gamma df = 0.
de Rham. HdR1(M)={closed}/{exact}H^1_{dR}(M)=\{\text{closed}\}/\{\text{exact}\}. For the punctured plane it is one-dimensional, spanned by the class of dθd\theta, and 12π∫γdθ\frac1{2\pi}\int_\gamma d\theta is the winding number.

Explore

Draw a loop with your mouse or finger. It closes itself when you let go. Arrows show the vector field dual to ω\omega.
The 1-form ω
Or try a ready-made loop

Worked example

On ℝ2∖{0}\R^2\smallsetminus\{0\} let ω=−ydx+xdyx2+y2\displaystyle\omega=\frac{-y\,dx+x\,dy}{x^2+y^2}. Show that ω\omega is closed but not exact.

  1. Compute dωd\omega.

    With P=−y/r2P=-y/r^2 and Q=x/r2Q=x/r^2:

    ∂xQ=y2−x2(x2+y2)2,∂yP=−(x2+y2)+2y2(x2+y2)2=y2−x2(x2+y2)2.\partial_xQ = \frac{y^2-x^2}{(x^2+y^2)^2},\qquad \partial_yP = \frac{-(x^2+y^2)+2y^2}{(x^2+y^2)^2}=\frac{y^2-x^2}{(x^2+y^2)^2}.

    They are equal, so dω=0d\omega=0.

  2. Pull ω\omega back along the unit circle γ(t)=(cos⁡t,sin⁡t)\gamma(t)=(\cos t,\sin t) and integrate.

    γ*ω=(−sin⁡t⋅(−sin⁡t)+cos⁡t⋅cos⁡t)dt=dt,∫S1ω=∫02πdt=2π.\gamma^*\omega = \big(-\sin t\cdot(-\sin t)+\cos t\cdot\cos t\big)\,dt = dt, \qquad \int_{S^1}\omega = \int_0^{2\pi}dt = 2\pi.
  3. Why does this show ω\omega is not exact?

    If ω=df\omega=df for a smooth ff on ℝ2∖{0}\R^2\smallsetminus\{0\}, then ∫γdf=f(γ(2π))−f(γ(0))=0\int_\gamma df=f(\gamma(2\pi))-f(\gamma(0))=0 for the closed curve γ\gamma. But the integral is 2π2\pi.

  4. Yet ω\omega "is" dθd\theta. Where does θ\theta live?

    On the plane minus any ray from the origin there is a smooth angle function θ\theta with dθ=ωd\theta=\omega. So ω\omega is locally exact, as the Poincaré lemma promises for every closed form. No single θ\theta works on the whole punctured plane. This obstruction is the nonzero class [ω]∈HdR1(ℝ2∖{0})[\omega]\in H^1_{dR}(\R^2\smallsetminus\{0\}).

Check your understanding

What is d(xdy)d(x\,dy)?

  1. dx∧dydx\wedge dy
  2. dy∧dxdy\wedge dxd(xdy)=dx∧dyd(x\,dy)=dx\wedge dy, and dy∧dxdy\wedge dx is its negative.
  3. 00xdyx\,dy is not closed; only exact forms are guaranteed to be.
  4. dx+dydx+dy
d(xdy)=dx∧dyd(x\,dy)=dx\wedge dy. By Stokes, ∮γxdy\oint_\gamma x\,dy is the area inside γ\gamma when it runs counterclockwise. Choose xdyx\,dy in the explorer to see it.

A loop winds twice counterclockwise around the origin. What is ∮dθ\oint d\theta?

  1. 2π2\pi
  2. 4π4\pi
  3. 00That's what an exact form would give, and dθd\theta isn't exact.
Each full turn contributes 2π2\pi. Use the "twice around" loop in the explorer.

ω\omega is a closed 1-form on all of ℝ2\R^2. Must it be exact?

  1. Yes
  2. Not necessarily: dθd\theta is a counterexample.dθd\theta is not defined at the origin, so it is not a form on all of ℝ2\R^2.
ℝ2\R^2 is star-shaped, so the Poincaré lemma gives HdR1(ℝ2)=0H^1_{dR}(\R^2)=0: every closed form is exact.

For a smooth map FF and a form ω\omega, how do dd and pullback interact?

  1. F*(dω)=d(F*ω)F^*(d\omega)=d(F^*\omega)
  2. F*(dω)=−d(F*ω)F^*(d\omega)=-d(F^*\omega)
  3. Only when FF is a diffeomorphism
dd commutes with pullback by any smooth map. This is why closed forms pull back to closed forms and exact forms to exact forms, so pullback is defined on cohomology.

On the plane with two punctures aa and bb, what is dim⁡HdR1\dim H^1_{dR}?

  1. 1
  2. 2
  3. 0
The angle forms dθad\theta_a and dθbd\theta_b around each puncture are independent: integrals around small loops tell them apart. Mayer–Vietoris confirms the dimension is exactly 2. Try the two-puncture form in the explorer.

Practice

Area of an ellipse by StokesGreen's theorem as a special case

Compute ∮E(xdy−ydx)\oint_E (x\,dy-y\,dx) around the ellipse x2/a2+y2/b2=1x^2/a^2+y^2/b^2=1, counterclockwise, in two ways.

Solution

d(xdy−ydx)=2dx∧dyd(x\,dy-y\,dx)=2\,dx\wedge dy, so by Stokes the integral is twice the area: 2πab2\pi ab. Directly, with x=acos⁡t, y=bsin⁡tx=a\cos t,\ y=b\sin t, the pullback is abdtab\,dt, and integrating over [0,2π][0,2\pi] gives 2πab2\pi ab.

An exterior derivative in ℝ³Computing d

Compute dωd\omega for ω=xydz+exdy\omega = xy\,dz + e^x\,dy, and check that d(dω)=0d(d\omega)=0.

Solution

dω=ydx∧dz+xdy∧dz+exdx∧dyd\omega = y\,dx\wedge dz + x\,dy\wedge dz + e^x\,dx\wedge dy. Then d(dω)=dy∧dx∧dz+dx∧dy∧dz+0=0d(d\omega)=dy\wedge dx\wedge dz+dx\wedge dy\wedge dz+0=0, since dy∧dx=−dx∧dydy\wedge dx=-dx\wedge dy.

Detecting exactnessOne number decides it

Show that a closed 1-form η\eta on ℝ2∖{0}\R^2\smallsetminus\{0\} is exact if and only if ∮S1η=0\oint_{S^1}\eta=0.

Hint

Consider η−c2πdθ\eta-\frac{c}{2\pi}d\theta with c=∮S1ηc=\oint_{S^1}\eta. Show its integral over every closed curve vanishes.

Solution

Exact forms integrate to 0 over closed curves, which gives one direction. Conversely, every closed curve in the punctured plane is homotopic to a multiple of the circle, and integrals of closed forms are homotopy invariant. So if ∮S1η=0\oint_{S^1}\eta=0, every loop integral vanishes, and f(x)=∫x0xηf(x)=\int_{x_0}^x\eta along any path is well defined with df=ηdf=\eta.