Survey of Mathematics

6.7 Differential forms and the one big theorem

A differential form is the thing under an integral sign, a recipe that eats tiny oriented pieces of a curve, surface or solid; the wedge product and the exterior derivative d turn gradient, curl and divergence into one operation, d squared is zero says curl grad and div curl vanish, and a single theorem, the integral of d omega over a region equals the integral of omega over its boundary, contains the fundamental theorem of calculus, Green, Stokes and Gauss.

Beyond the courseRudin, Principles of Mathematical Analysis, ch. 10 (integration of differential forms)Spivak, Calculus on Manifolds, ch. 4Apostol Vol. 2, 12.9 (other notations for surface integrals)MATH 376: the notation P dx + Q dy of 6.1 and 6.2 and the exactness test of Midterm 1Cartan, Sur certaines expressions differentielles et le probleme de Pfaff (1899)

In 6.1 you wrote ∫CP dx+Q dy\int_CP\,dx + Q\,dy and in 6.2 you called P dx+Q dyP\,dx + Q\,dy an exact differential. We never said what dxdx is. This chapter does. The answer unifies the whole unit. An object called a differential form sits under every integral sign of Unit 6; a single operation dd turns gradient, curl and divergence into one rule; the identities curl grad =0= 0 and div curl =0= 0 become one identity, d2=0d^2 = 0; and the four theorems of the unit, the fundamental theorem of calculus, Green, Stokes and Gauss, become one theorem. None of this is needed for the MATH 376 exam. It is here so that the pattern is visible, because it is what a course on manifolds starts from.

Before you start

6.6: curl⁡∇f=0\operatorname{curl}\nabla f = 0 and div⁡curl⁡F=0\operatorname{div}\operatorname{curl}F = 0. 6.5: flux through a parametrized surface is ∬F⋅(ru×rv) du dv\iint F\cdot(r_u\times r_v)\,du\,dv. 6.4: the Jacobian determinant is the area stretching factor. 6.1: ∫CP dx+Q dy=∫(Px′+Qy′) dt\int_CP\,dx + Q\,dy = \int(Px' + Qy')\,dt.

Quick checkFor a path r(t)=(x(t),y(t))r(t) = (x(t), y(t)), the symbol ∫CP dx+Q dy\int_CP\,dx + Q\,dy means which integral?
Quick checkWhat is curl⁡(∇f)\operatorname{curl}(\nabla f) for a C2C^2 function ff on R3\mathbb R^3?
Quick checkA map sends the parameter square with area AA in the uvuv-plane to a region in the xyxy-plane. Its Jacobian determinant is 33 everywhere. The area of the region is

The picture: what a form eats

Think about what an integral sign needs. In ∫CF⋅dr\int_CF\cdot dr it needs, at each point of the curve, a small displacement vector, and it returns a number: the part of FF along the displacement. In ∬SF⋅n dS\iint_SF\cdot n\,dS it needs a tiny piece of surface, a tiny parallelogram with two edge vectors aa and bb in a definite order, and it returns a number: the amount of FF crossing the piece, F⋅(a×b)F\cdot(a\times b). In ∭f dV\iiint f\,dV it needs a tiny parallelepiped with three edges, and returns ff times the signed volume det⁡[a,b,c]\det[a, b, c]. A differential form is exactly that: a rule that eats tiny pieces and returns numbers, linear in each edge, and changing sign when two edges are swapped (because swapping edges reverses orientation).

Here is the whole structure in a table, for three dimensions. The columns are the degree of the form, what it eats, and what vector calculus calls it.

DegreeTypical formEatsWhat it returnsSame as
0ffa pointff at the pointa scalar field
1P dx+Q dy+R dzP\,dx + Q\,dy + R\,dza vector vvPv1+Qv2+Rv3=F⋅vP v_1 + Q v_2 + R v_3 = F\cdot va vector field F=(P,Q,R)F=(P,Q,R), via the line integral
2A dy∧dz+B dz∧dx+C dx∧dyA\,dy\wedge dz + B\,dz\wedge dx + C\,dx\wedge dytwo vectors a,ba, b(A,B,C)⋅(a×b)(A, B, C)\cdot(a\times b)a vector field G=(A,B,C)G = (A,B,C), via the flux
3f dx∧dy∧dzf\,dx\wedge dy\wedge dzthree vectors a,b,ca, b, cf det⁡[a,b,c]f\,\det[a, b, c]a scalar field, via the volume

Read down the ‘same as’ column and notice the pattern: scalars and vector fields appear twice. That is why gradient, curl and divergence look unrelated and are not: they are one operation, dd, acting on degrees 0→10\to1, 1→21\to2, 2→32\to3.

Where this came from· Where forms came from

Hermann Grassmann in 1844 built an algebra of geometric objects in which the product of two vectors depends on their order, with ab=−baab = -ba; it is the algebra of the wedge product below. The step to calculus was Élie Cartan’s, in a long 1899 paper in the Annales scientifiques de l’École Normale Supérieure, ‘Sur certaines expressions différentielles et le problème de Pfaff’. The Pfaff problem asks: given an expression a1 dx1+⋯+an dxna_1\,dx_1 + \dots + a_n\,dx_n, can a change of variables put it in the simplest possible form? Cartan introduced the operation we now call the exterior derivative to answer it, and the paper is where it first appears. The unification of Green, Stokes and Gauss into one statement came later, as the language of forms spread through geometry in the twentieth century.

Forms, and the wedge product

Definition B

One-form

At each point pp, dxdx is the linear function that sends a vector v=(v1,v2,v3)v = (v_1, v_2, v_3) to v1v_1, dydy sends it to v2v_2, and dzdz to v3v_3. A 1-form is ω=P dx+Q dy+R dz\omega = P\,dx + Q\,dy + R\,dz, where P,Q,RP, Q, R are functions of position: at pp it sends vv to P(p)v1+Q(p)v2+R(p)v3=F(p)⋅vP(p)v_1 + Q(p)v_2 + R(p)v_3 = F(p)\cdot v, with F=(P,Q,R)F = (P, Q, R). So 1-forms and vector fields carry the same information. The line integral of the vector field FF along a path is the integral of the 1-form ωF=F⋅dr\omega_F = F\cdot dr.

In one lineIn R3\mathbb R^3 a 1-form is an expression ω=P dx+Q dy+R dz\omega = P\,dx + Q\,dy + R\,dz with functions P,Q,RP, Q, R; at each point it takes a vector vv to the number Pv1+Qv2+Rv3Pv_1 + Qv_2 + Rv_3.

Why it existsIt is the thing under a line integral. dxdx is the 1-form that reads off the xx component of a displacement, so P dx+Q dyP\,dx + Q\,dy is 'the part of the field along the step'.

Definition B

Two-form

The symbol dx∧dydx\wedge dy (’dxdx wedge dydy’) is the 2-form that sends a pair of vectors (a,b)(a, b) to a1b2−a2b1a_1b_2 - a_2b_1, the signed area of the shadow of the parallelogram on the xyxy-plane. Likewise dy∧dzdy\wedge dz sends (a,b)(a, b) to a2b3−a3b2a_2b_3 - a_3b_2 and dz∧dxdz\wedge dx to a3b1−a1b3a_3b_1 - a_1b_3. These are exactly the three components of a×ba\times b. A 2-form is a combination ω=A dy∧dz+B dz∧dx+C dx∧dy\omega = A\,dy\wedge dz + B\,dz\wedge dx + C\,dx\wedge dy, and it sends (a,b)(a, b) to (A,B,C)⋅(a×b)(A, B, C)\cdot(a\times b). The flux of G=(A,B,C)G = (A, B, C) through a surface is the integral of the 2-form A dy∧dz+B dz∧dx+C dx∧dyA\,dy\wedge dz + B\,dz\wedge dx + C\,dx\wedge dy. This is the notation of Apostol 12.9, which the course listed, and it is the one noted in 6.5.

In one lineIn R3\mathbb R^3 a 2-form is A dy∧dz+B dz∧dx+C dx∧dyA\,dy\wedge dz + B\,dz\wedge dx + C\,dx\wedge dy; at each point it takes two vectors a,ba, b to (A,B,C)⋅(a×b)(A, B, C)\cdot(a\times b), the amount of the field (A,B,C)(A, B, C) crossing the parallelogram on aa and bb.

Why it existsIt is the thing under a surface integral: it eats a tiny oriented parallelogram. The sign changes when you swap aa and bb, which is what orientation means.

Definition B

Differential form

A differential form of degree kk (a kk-form) on a region of Rn\mathbb R^n is a sum ω=∑fi1…ik dxi1∧⋯∧dxik\omega = \sum f_{i_1\dots i_k}\,dx_{i_1}\wedge\dots\wedge dx_{i_k} over increasing index lists, with functions fi1…ikf_{i_1\dots i_k}. At each point it eats kk vectors and returns a number, linear in each vector and changing sign when two vectors are swapped. In R3\mathbb R^3 the degrees are 0,1,2,30, 1, 2, 3; in R2\mathbb R^2 they are 0,1,20, 1, 2 with the only 2-form f dx∧dyf\,dx\wedge dy; there are no kk-forms for kk larger than nn.

In one lineA kk-form is a rule that assigns to each point a function of kk vectors that is linear in each and changes sign when two of them are swapped. In coordinates it is a sum of functions times dxi1∧⋯∧dxikdx_{i_1}\wedge\dots\wedge dx_{i_k}. A 0-form is a function.

Why it existsIt is the single object that can be integrated over a kk-dimensional piece, for every kk: points, curves, surfaces, solids.

Definition B

Wedge product

The wedge product ω∧η\omega\wedge\eta of a jj-form and a kk-form is a (j+k)(j + k)-form. The rules: it is bilinear (it distributes over sums and constants and functions come out), it is associative, and the basic rule is that swapping two 1-forms flips the sign: dx∧dy=−dy∧dxdx\wedge dy = -dy\wedge dx, and hence dx∧dx=0dx\wedge dx = 0 (since dx∧dx=−dx∧dxdx\wedge dx = -dx\wedge dx). Two 1-forms always anticommute; a function (0-form) commutes with everything. If you evaluate on vectors, (dx∧dy)(a,b)=dx(a) dy(b)−dx(b) dy(a)=a1b2−a2b1(dx\wedge dy)(a, b) = dx(a)\,dy(b) - dx(b)\,dy(a) = a_1b_2 - a_2b_1.

In one lineThe product of forms, bilinear and associative, with ω∧η=(−1)jkη∧ω\omega\wedge\eta = (-1)^{jk}\eta\wedge\omega for degrees j,kj, k; in particular dxi∧dxj=−dxj∧dxidx_i\wedge dx_j = -dx_j\wedge dx_i and dxi∧dxi=0dx_i\wedge dx_i = 0.

Why it existsIt is how a 1-form eats two vectors at once (it is the signed area of a shadow) and how forms of degrees jj and kk combine into degree j+kj + k.

Everything about the wedge product follows from dxi∧dxj=−dxj∧dxidx_i\wedge dx_j = -dx_j\wedge dx_i. Expand, kill the terms with a repeated dxidx_i, and reorder with a sign.

Worked example
Compute (2 dx−3 dy)∧(dx+4 dy)(2\,dx - 3\,dy)\wedge(dx + 4\,dy) in the plane, and (dx+2 dy+3 dz)∧(dy+dz)(dx + 2\,dy + 3\,dz)\wedge(dy + dz) in space. Compare the second with a cross product.

Solution.

  1. Plane, expand. 2 dx∧dx+8 dx∧dy−3 dy∧dx−12 dy∧dy2\,dx\wedge dx + 8\,dx\wedge dy - 3\,dy\wedge dx - 12\,dy\wedge dy. Distribute: four products.
  2. Plane, simplify. dx∧dx=0=dy∧dydx\wedge dx = 0 = dy\wedge dy and −3 dy∧dx=+3 dx∧dy-3\,dy\wedge dx = +3\,dx\wedge dy, so the result is (8+3) dx∧dy=11 dx∧dy(8 + 3)\,dx\wedge dy = 11\,dx\wedge dy. This is the determinant det⁡[2−314]=11\det\begin{bmatrix}2&-3\\1&4\end{bmatrix} = 11: the wedge of two 1-forms in the plane is a 2×22\times2 determinant.
  3. Space, expand. dx∧dy+dx∧dz+2 dy∧dy+2 dy∧dz+3 dz∧dy+3 dz∧dzdx\wedge dy + dx\wedge dz + 2\,dy\wedge dy + 2\,dy\wedge dz + 3\,dz\wedge dy + 3\,dz\wedge dz.
  4. Space, simplify. Drop dy∧dydy\wedge dy and dz∧dzdz\wedge dz. Rewrite dx∧dz=−dz∧dxdx\wedge dz = -dz\wedge dx and 3 dz∧dy=−3 dy∧dz3\,dz\wedge dy = -3\,dy\wedge dz. The result is (2−3) dy∧dz−dz∧dx+dx∧dy=−dy∧dz−dz∧dx+dx∧dy(2 - 3)\,dy\wedge dz - dz\wedge dx + dx\wedge dy = -dy\wedge dz - dz\wedge dx + dx\wedge dy.
  5. Compare. In the dictionary this is the vector (−1,−1,1)(-1, -1, 1), and (1,2,3)×(0,1,1)=(2−3, 0−1, 1−0)=(−1,−1,1)(1, 2, 3)\times(0, 1, 1) = (2 - 3,\ 0 - 1,\ 1 - 0) = (-1, -1, 1).

Answer. 11 dx∧dy11\,dx\wedge dy and −dy∧dz−dz∧dx+dx∧dy-dy\wedge dz - dz\wedge dx + dx\wedge dy. The wedge product of two 1-forms is the cross product: that is why the cross product is anti-symmetric and why it exists only in three dimensions as a vector.

Quick checkWhat is dy∧dx+dx∧dydy\wedge dx + dx\wedge dy?
Quick checkWhat is (x dx)∧(x dx)(x\,dx)\wedge(x\,dx)?

The exterior derivative

Definition B

Exterior derivative

The exterior derivative dd takes a kk-form to a (k+1)(k+1)-form. For a function ff: df=fx dx+fy dy+fz dzdf = f_x\,dx + f_y\,dy + f_z\,dz (the total differential of 5.3). For a form ω=f dxi1∧⋯∧dxik\omega = f\,dx_{i_1}\wedge\dots\wedge dx_{i_k} with one term: dω=df∧dxi1∧⋯∧dxikd\omega = df\wedge dx_{i_1}\wedge\dots\wedge dx_{i_k}. For a sum, add. Because dxi∧dxi=0dx_i\wedge dx_i = 0, only the derivatives of ff in directions not already present survive.

In one lineThe operation dd that raises the degree of a form by one: on a function df=fx dx+fy dy+fz dzdf = f_x\,dx + f_y\,dy + f_z\,dz, and on f dxIf\,dx_I it is df∧dxIdf\wedge dx_I, extended by linearity.

Why it existsIt is the one derivative that is built to be integrated: ∫Mdω\int_Md\omega over a region equals ∫∂Mω\int_{\partial M}\omega over its boundary. It contains gradient, curl and divergence.

Theorem B Rudin 10.20 and the dictionary of vector calculus; Spivak ch. 4

d is gradient, curl and divergence

Under the dictionary of the table above, dd on 0-forms is the gradient, on 1-forms the curl, and on 2-forms the divergence.

In one lineIn R3\mathbb R^3: df=fx dx+fy dy+fz dzdf = f_x\,dx + f_y\,dy + f_z\,dz (gradient); d(P dx+Q dy+R dz)=(Ry−Qz) dy∧dz+(Pz−Rx) dz∧dx+(Qx−Py) dx∧dyd(P\,dx + Q\,dy + R\,dz) = (R_y - Q_z)\,dy\wedge dz + (P_z - R_x)\,dz\wedge dx + (Q_x - P_y)\,dx\wedge dy (curl); d(A dy∧dz+B dz∧dx+C dx∧dy)=(Ax+By+Cz) dx∧dy∧dzd(A\,dy\wedge dz + B\,dz\wedge dx + C\,dx\wedge dy) = (A_x + B_y + C_z)\,dx\wedge dy\wedge dz (divergence).

Why it is trueIt shows that the three operators of 6.6 are one operator, and it makes their relations, such as the product rules, automatic.

ProofRudin 10.20
The idea: Apply the rule d(f dx_I) = df wedge dx_I term by term, then use dx wedge dx = 0 and the anticommuting rule to put the survivors in the standard order.
df=fx dx+fy dy+fz dzdf = f_x\,dx + f_y\,dy + f_z\,dz, the gradient ∇f=(fx,fy,fz)\nabla f = (f_x, f_y, f_z) in the dictionary.
The 0-form case is the definition.
dω=dP∧dx+dQ∧dy+dR∧dzd\omega = dP\wedge dx + dQ\wedge dy + dR\wedge dz.
Apply the rule to each of the three terms of ω=P dx+Q dy+R dz\omega = P\,dx + Q\,dy + R\,dz.
=(Py dy+Pz dz)∧dx+(Qx dx+Qz dz)∧dy+(Rx dx+Ry dy)∧dz= (P_y\,dy + P_z\,dz)\wedge dx + (Q_x\,dx + Q_z\,dz)\wedge dy + (R_x\,dx + R_y\,dy)\wedge dz.
Expand each differential: dP=Px dx+Py dy+Pz dzdP = P_x\,dx + P_y\,dy + P_z\,dz, and so on.
=(Ry−Qz) dy∧dz+(Pz−Rx) dz∧dx+(Qx−Py) dx∧dy= (R_y - Q_z)\,dy\wedge dz + (P_z - R_x)\,dz\wedge dx + (Q_x - P_y)\,dx\wedge dy.
Terms with a repeated dxidx_i vanished already (Px dx∧dx=0P_x\,dx\wedge dx = 0 and so on). Reorder with dy∧dx=−dx∧dydy\wedge dx = -dx\wedge dy, dz∧dy=−dy∧dzdz\wedge dy = -dy\wedge dz and dx∧dz=−dz∧dxdx\wedge dz = -dz\wedge dx.The three survivors in the standard order dy∧dzdy\wedge dz, dz∧dxdz\wedge dx, dx∧dydx\wedge dy are the cyclic order that matches the cross product.
Add the three: dω=(Ax+By+Cz) dx∧dy∧dzd\omega = (A_x + B_y + C_z)\,dx\wedge dy\wedge dz. ■\blacksquare
For a 2-form, only the derivative in the missing direction survives: d(A dy∧dz)=Ax dx∧dy∧dzd(A\,dy\wedge dz) = A_x\,dx\wedge dy\wedge dz since Ay dy∧dy∧dz=0A_y\,dy\wedge dy\wedge dz = 0 and Az dz∧dy∧dz=0A_z\,dz\wedge dy\wedge dz = 0. Similarly d(B dz∧dx)=By dy∧dz∧dx=By dx∧dy∧dzd(B\,dz\wedge dx) = B_y\,dy\wedge dz\wedge dx = B_y\,dx\wedge dy\wedge dz (a cyclic shift is an even permutation) and d(C dx∧dy)=Cz dz∧dx∧dy=Cz dx∧dy∧dzd(C\,dx\wedge dy) = C_z\,dz\wedge dx\wedge dy = C_z\,dx\wedge dy\wedge dz.

In the plane the same computation gives d(P dx+Q dy)=(Qx−Py) dx∧dyd(P\,dx + Q\,dy) = (Q_x - P_y)\,dx\wedge dy: the integrand of Green’s theorem. In one variable d(f)=f′ dxd(f) = f'\,dx.

Worked example new (compare MATH 376 worksheets 13, 14 and 28)
Compute dd of: (a) f=x2yzf = x^2yz; (b) ω=−y dx+x dy\omega = -y\,dx + x\,dy; (c) ω=yz dx+xz dy+xy dz\omega = yz\,dx + xz\,dy + xy\,dz; (d) ω=x2 dy∧dz+y2 dz∧dx+z2 dx∧dy\omega = x^2\,dy\wedge dz + y^2\,dz\wedge dx + z^2\,dx\wedge dy.

Solution.

  1. (a) df=2xyz dx+x2z dy+x2y dzdf = 2xyz\,dx + x^2z\,dy + x^2y\,dz. The three partial derivatives, as in 5.2.
  2. (b) dω=(∂x(x)−∂y(−y)) dx∧dy=2 dx∧dyd\omega = (\partial_x(x) - \partial_y(-y))\,dx\wedge dy = 2\,dx\wedge dy. Qx−Py=1−(−1)Q_x - P_y = 1 - (-1). Green’s theorem will turn this into twice the area.
  3. (c) P=yzP = yz, Q=xzQ = xz, R=xyR = xy: Ry−Qz=x−x=0R_y - Q_z = x - x = 0, Pz−Rx=y−y=0P_z - R_x = y - y = 0, Qx−Py=z−z=0Q_x - P_y = z - z = 0. So dω=0d\omega = 0. It is the field of worksheet 28: curl zero, and indeed ω=d(xyz)\omega = d(xyz).
  4. (d) A=x2A = x^2, B=y2B = y^2, C=z2C = z^2, so dω=(2x+2y+2z) dx∧dy∧dzd\omega = (2x + 2y + 2z)\,dx\wedge dy\wedge dz. The divergence of (x2,y2,z2)(x^2, y^2, z^2).

Answer. As in steps 1 to 4. Note (c): a form with dω=0d\omega = 0 that is also dd of something is the situation of the next theorem.

Theorem B Rudin 10.20(b); Spivak ch. 4

d squared is zero

If ω\omega has C2C^2 coefficients then d(dω)=0d(d\omega) = 0.

In one lineFor every C2C^2 form ω\omega, d(dω)=0d(d\omega) = 0. For ω=f\omega = f a function in R3\mathbb R^3 this is curl⁡∇f=0\operatorname{curl}\nabla f = 0, and for a 1-form it is div⁡curl⁡F=0\operatorname{div}\operatorname{curl}F = 0.

Why it is trueOne identity, one proof, replaces both identities of 6.6. It is the reason closed and exact forms are different questions, and the starting point of de Rham cohomology.

ProofRudin 10.20
The idea: Differentiating twice gives second partial derivatives f_ij paired with dx_j wedge dx_i. Each pair appears twice with opposite signs: once as f_ij dx_j wedge dx_i and once as f_ji dx_i wedge dx_j, and f_ij = f_ji while the wedge is anti-symmetric.
For a function ff: df=∑ifi dxidf = \sum_if_i\,dx_i, where fi=∂f/∂xif_i = \partial f/\partial x_i.
Definition: dd of a function is the sum of its partials times the basic 1-forms.
d(df)=∑i∑jfij dxj∧dxid(df) = \sum_i\sum_jf_{ij}\,dx_j\wedge dx_i, where fij=∂j∂iff_{ij} = \partial_j\partial_if.
Apply dd to each term fi dxif_i\,dx_i: d(fi dxi)=dfi∧dxid(f_i\,dx_i) = df_i\wedge dx_i.
=∑j<ifij dxj∧dxi+∑j<ifji dxi∧dxj= \sum_{j<i}f_{ij}\,dx_j\wedge dx_i + \sum_{j<i}f_{ji}\,dx_i\wedge dx_j.
Split the double sum into the terms with j<ij<i and j>ij>i (the terms with j=ij = i vanish because dxi∧dxi=0dx_i\wedge dx_i = 0), and rename in the second sum.
=∑j<ifij(dxj∧dxi−dxj∧dxi)=0= \sum_{j<i}f_{ij}\big(dx_j\wedge dx_i - dx_j\wedge dx_i\big) = 0.
Equality of mixed partials for C2C^2 functions (5.5): fij=fjif_{ij} = f_{ji}. And dxi∧dxj=−dxj∧dxidx_i\wedge dx_j = -dx_j\wedge dx_i.These are the same two facts as in the proofs of 6.6, now stated once.
So d(dω)=0d(d\omega) = 0 for every C2C^2 form. ■\blacksquare
For a general form f dxIf\,dx_I: d(f dxI)=df∧dxId(f\,dx_I) = df\wedge dx_I, and applying dd again gives d(df)∧dxI=0d(df)\wedge dx_I = 0 by the previous steps, because d(dxI)=0d(dx_I) = 0 (dxIdx_I has constant coefficients).
Proof workshopRebuild the proof that d(df) = 0 for a function of three variablesRudin 10.20

Tap the lines in the order a proof would write them. Each line can only use what is already above it.

    Why· One identity for two

    Under the dictionary, d∘d=0d\circ d = 0 on a function is curl⁡∘∇=0\operatorname{curl}\circ\nabla = 0, and on a 1-form it is div⁡∘curl⁡=0\operatorname{div}\circ\operatorname{curl} = 0. In the plane, d(df)=(fyx−fxy) dx∧dy=0d(df) = (f_{yx} - f_{xy})\,dx\wedge dy = 0 is the test of 6.2: Py=QxP_y = Q_x for P=fxP = f_x, Q=fyQ = f_y. So the ‘necessary condition for a gradient field’ of 6.2 is now one line: if ω=df\omega = df then dω=d(df)=0d\omega = d(df) = 0.

    Closed and exact

    Definition B

    Closed form

    A form ω\omega is closed if dω=0d\omega = 0. For ω=P dx+Q dy\omega = P\,dx + Q\,dy this is Qx−Py=0Q_x - P_y = 0. For ω=P dx+Q dy+R dz\omega = P\,dx + Q\,dy + R\,dz it is curl⁡F=0\operatorname{curl}F = 0. For a 2-form it is div⁡G=0\operatorname{div}G = 0.

    In one lineA form ω\omega with dω=0d\omega = 0.

    Why it existsThis is the test that can be computed. For a 1-form in the plane it is Py=QxP_y = Q_x, the necessary condition of 6.2.

    Definition B

    Exact form

    A form ω\omega is exact if ω=dη\omega = d\eta for some form η\eta defined on the same region. A 1-form ωF\omega_F is exact exactly when FF is a gradient field, F=∇fF = \nabla f (6.2). The line integral of an exact 1-form between two points depends only on the points. A first order equation P dx+Q dy=0P\,dx + Q\,dy = 0 is an ‘exact differential equation’ (6.2) when the form P dx+Q dyP\,dx + Q\,dy is exact.

    In one lineA form ω\omega that is dηd\eta for some form η\eta (for a 1-form: ω=df\omega = df, a gradient).

    Why it existsThis is what we want: a potential. In 6.2 it meant a gradient field. By d2=0d^2 = 0, every exact form is closed, and the question is when the converse holds.

    By d2=0d^2 = 0 every exact form is closed. The converse is the question of 6.2 in new clothes, and the answer depends on the shape of the region.

    Theorem B Rudin ch. 10; Spivak ch. 4; stated for 1-forms in 6.2

    Poincare lemma

    Let U⊂RnU\subset\mathbb R^n be star-shaped (there is a point pp such that the segment from pp to any point of UU lies in UU) and ω\omega a closed C1C^1 form of degree k≥1k\ge1 on UU. Then ω=dη\omega = d\eta for a form η\eta of degree k−1k-1 on UU.

    In one lineOn a star-shaped region (such as a ball or all of Rn\mathbb R^n) every closed form of degree ≥1\ge1 is exact.

    Why it is trueIt is the converse of d2=0d^2 = 0 on nice regions, and it explains why a gradient field exists wherever curl is zero on R3\mathbb R^3 and why a vector potential exists wherever divergence is zero on R3\mathbb R^3.

    For 1-forms the proof is a formula, the integral along rays from the centre; the stretch problems carry it out. On a region with a hole the lemma fails, and this is exactly where closed and exact part ways.

    Worked example new (compare MATH 376 homework 5 problem 4.1 and 6.2)
    The angle form. Show that ω=−y dx+x dyx2+y2\omega = \dfrac{-y\,dx + x\,dy}{x^2 + y^2} is closed on the punctured plane, and that it is not exact there.

    Solution.

    1. Components. P=−yx2+y2P = \dfrac{-y}{x^2 + y^2}, Q=xx2+y2Q = \dfrac x{x^2 + y^2}.
    2. Partial derivatives. Py=−(x2+y2)+2y2(x2+y2)2=y2−x2(x2+y2)2P_y = \dfrac{-(x^2 + y^2) + 2y^2}{(x^2 + y^2)^2} = \dfrac{y^2 - x^2}{(x^2 + y^2)^2} and Qx=(x2+y2)−2x2(x2+y2)2=y2−x2(x2+y2)2Q_x = \dfrac{(x^2 + y^2) - 2x^2}{(x^2 + y^2)^2} = \dfrac{y^2 - x^2}{(x^2 + y^2)^2}.
    3. Closed. dω=(Qx−Py) dx∧dy=0d\omega = (Q_x - P_y)\,dx\wedge dy = 0 on R2∖{0}\mathbb R^2\setminus\{0\}.
    4. Not exact. On the unit circle r(t)=(cos⁡t,sin⁡t)r(t) = (\cos t, \sin t) we have ∫ω=∫02π(sin⁡2t+cos⁡2t) dt=2π≠0\int\omega = \int_0^{2\pi}(\sin^2t + \cos^2t)\,dt = 2\pi\ne0. If ω=df\omega = df, the line integral around every closed path would be 00 (6.2).
    5. What it is. In polar coordinates ω=dθ\omega = d\theta, the change in angle. A function whose differential is dθd\theta would be θ\theta itself, which increases by 2π2\pi round the origin and cannot be a single-valued function.

    Answer. Closed, not exact on the punctured plane (though exact on any region not circling the origin, such as the right half plane, where ω=darctan⁡(y/x)\omega = d\arctan(y/x)). The hole is what makes ‘closed’ differ from ‘exact’. Counting holes this way is called de Rham cohomology.

    Worked example new (MATH 376 Midterm 1 problem 1 style, 6.2)
    Decide whether ω=(2xy+y3) dx+(x2+3xy2) dy\omega = (2xy + y^3)\,dx + (x^2 + 3xy^2)\,dy is closed and exact on R2\mathbb R^2, and find a potential. Then do the same for y dx−x dyy\,dx - x\,dy.

    Solution.

    1. Closed test. Py=2x+3y2P_y = 2x + 3y^2 and Qx=2x+3y2Q_x = 2x + 3y^2. They are equal, so dω=0d\omega = 0.
    2. Exact. The plane is star-shaped, so by the Poincare lemma ω\omega is exact.
    3. Find the potential. fx=2xy+y3f_x = 2xy + y^3 gives f=x2y+xy3+h(y)f = x^2y + xy^3 + h(y). Then fy=x2+3xy2+h′(y)f_y = x^2 + 3xy^2 + h'(y) must equal Q=x2+3xy2Q = x^2 + 3xy^2, so h′=0h' = 0. The method of 6.2: integrate one equation, differentiate, compare.
    4. Check. df=(2xy+y3) dx+(x2+3xy2) dydf = (2xy + y^3)\,dx + (x^2 + 3xy^2)\,dy.
    5. Second form. y dx−x dyy\,dx - x\,dy has Py=1P_y = 1 and Qx=−1Q_x = -1, so dω=(Qx−Py) dx∧dy=−2 dx∧dy≠0d\omega = (Q_x - P_y)\,dx\wedge dy = -2\,dx\wedge dy\ne0. It is not closed, hence not exact: ∫\int around a counterclockwise circle of radius rr is −2πr2-2\pi r^2.

    Answer. f=x2y+xy3f = x^2y + xy^3 (plus a constant) for the first; the second is not closed.

    Your answerLet ω=(3x2y+z) dx+x3 dy+x dz\omega = (3x^2y + z)\,dx + x^3\,dy + x\,dz. Find ∫Cω\int_C\omega along any path from (0,0,0)(0, 0, 0) to (1,2,3)(1, 2, 3).

    Pullback: how a form is integrated

    To integrate a form over a curve or surface given by a parametrization, we substitute the parametrization into the form. That substitution has a name.

    Definition B

    Pullback

    For φ(u,v)=(x(u,v),y(u,v),z(u,v))\varphi(u, v) = (x(u, v), y(u, v), z(u, v)), the pullback φ∗ω\varphi^*\omega is the form in du,dvdu, dv obtained by replacing x,y,zx, y, z in the coefficients by x(u,v),y(u,v),z(u,v)x(u, v), y(u, v), z(u, v), and each dxdx by xu du+xv dvx_u\,du + x_v\,dv (similarly dydy, dzdz), and expanding with du∧du=0du\wedge du = 0 and dv∧du=−du∧dvdv\wedge du = -du\wedge dv. For a curve r(t)r(t) the same with dx=x′(t) dtdx = x'(t)\,dt.

    In one lineIf φ\varphi is a map from parameters (u,v)(u, v) to (x,y,z)(x, y, z), the pullback φ∗ω\varphi^*\omega of a form ω\omega is obtained by substituting x=x(u,v)x = x(u, v) and so on into the coefficients and dx=xu du+xv dvdx = x_u\,du + x_v\,dv into the basic forms, then expanding with the wedge rules.

    Why it existsIt is the exact meaning of 'substitute the parametrization', and it gives integrals over any parametrized piece in one rule.

    The integral of a form over a parametrized piece is the integral of the pullback over the parameter region, with its orientation: for a curve ∫Cω=∫abr∗ω\int_C\omega = \int_a^b r^*\omega, and for a surface ∬Sω=∬Tr∗ω\iint_S\omega = \iint_Tr^*\omega, where f du∧dvf\,du\wedge dv is integrated as ∬f du dv\iint f\,du\,dv.

    Proposition B Rudin ch. 10 and Spivak ch. 4; computation below

    Forms contain the integrals of 6.1 and 6.5

    Let F=(P,Q,R)F = (P, Q, R) and G=(A,B,C)G = (A, B, C).

    In one lineFor a curve, r∗(P dx+Q dy+R dz)=F(r)⋅r′ dtr^*(P\,dx + Q\,dy + R\,dz) = F(r)\cdot r'\,dt, so ∫Cω\int_C\omega is the line integral. For a surface, r∗(A dy∧dz+B dz∧dx+C dx∧dy)=G(r)⋅(ru×rv) du∧dvr^*(A\,dy\wedge dz + B\,dz\wedge dx + C\,dx\wedge dy) = G(r)\cdot(r_u\times r_v)\,du\wedge dv, so ∬Sω\iint_S\omega is the flux.

    Why it is trueIt shows that the two integrals of 6.1 and 6.5 are the same operation on different forms, and that the components of the fundamental vector product are the three wedge coefficients.

    Proofdirect computation
    The idea: Substitute and expand. The cross product appears because the coefficient of du wedge dv in dy wedge dz is the determinant y_u z_v minus y_v z_u.
    r∗(P dx+Q dy+R dz)=(P(r)x′+Q(r)y′+R(r)z′) dt=F(r)⋅r′ dtr^*(P\,dx + Q\,dy + R\,dz) = (P(r)x' + Q(r)y' + R(r)z')\,dt = F(r)\cdot r'\,dt.
    Pullback of a 1-form along a curve: dx=x′ dtdx = x'\,dt, dy=y′ dtdy = y'\,dt, dz=z′ dtdz = z'\,dt.
    r∗(dy∧dz)=(yu du+yv dv)∧(zu du+zv dv)r^*(dy\wedge dz) = (y_u\,du + y_v\,dv)\wedge(z_u\,du + z_v\,dv).
    Pullback of dy∧dzdy\wedge dz along a surface: dy=yu du+yv dvdy = y_u\,du + y_v\,dv and dz=zu du+zv dvdz = z_u\,du + z_v\,dv.
    =(yuzv−yvzu) du∧dv= (y_uz_v - y_vz_u)\,du\wedge dv.
    Expand; du∧du=dv∧dv=0du\wedge du = dv\wedge dv = 0 and dv∧du=−du∧dvdv\wedge du = -du\wedge dv.
    So r∗(dz∧dx)=(zuxv−zvxu) du∧dvr^*(dz\wedge dx) = (z_ux_v - z_vx_u)\,du\wedge dv and r∗(dx∧dy)=(xuyv−xvyu) du∧dvr^*(dx\wedge dy) = (x_uy_v - x_vy_u)\,du\wedge dv.
    This is the first component of ru×rv=(yuzv−zuyv, zuxv−xuzv, xuyv−yuxv)r_u\times r_v = (y_uz_v - z_uy_v,\ z_ux_v - x_uz_v,\ x_uy_v - y_ux_v). The other two pullbacks are the cyclic copies.
    r∗ω=G(r)⋅(ru×rv) du∧dvr^*\omega = G(r)\cdot(r_u\times r_v)\,du\wedge dv. ■\blacksquare
    Add with coefficients A,B,CA, B, C evaluated at r(u,v)r(u, v).

    So the flux integral of 6.5 is ∬G(r)⋅(ru×rv) du dv\iint G(r)\cdot(r_u\times r_v)\,du\,dv, and the sign: swapping the roles of uu and vv changes du∧dvdu\wedge dv into dv∧du=−du∧dvdv\wedge du = -du\wedge dv, which is why swapping the parameters reverses the orientation (6.5). The fundamental vector product is the pullback made visible.

    Worked example new (compare 6.4, polar and spherical Jacobians)
    Find the pullback of dx∧dydx\wedge dy under polar coordinates x=rcos⁡θx = r\cos\theta, y=rsin⁡θy = r\sin\theta, and of dx∧dy∧dzdx\wedge dy\wedge dz under spherical coordinates. Use it to explain the Jacobian factors.

    Solution.

    1. Differentials. dx=cos⁡θ dr−rsin⁡θ dθdx = \cos\theta\,dr - r\sin\theta\,d\theta and dy=sin⁡θ dr+rcos⁡θ dθdy = \sin\theta\,dr + r\cos\theta\,d\theta.
    2. Wedge. dx∧dy=cos⁡θ⋅rcos⁡θ dr∧dθ+(−rsin⁡θ)sin⁡θ dθ∧dr=rcos⁡2θ dr∧dθ+rsin⁡2θ dr∧dθdx\wedge dy = \cos\theta\cdot r\cos\theta\,dr\wedge d\theta + (-r\sin\theta)\sin\theta\,d\theta\wedge dr = r\cos^2\theta\,dr\wedge d\theta + r\sin^2\theta\,dr\wedge d\theta. The terms dr∧drdr\wedge dr and dθ∧dθd\theta\wedge d\theta vanish, and dθ∧dr=−dr∧dθd\theta\wedge dr = -dr\wedge d\theta.
    3. Simplify. dx∧dy=r dr∧dθdx\wedge dy = r\,dr\wedge d\theta.
    4. Spherical. With x=ρsin⁡φcos⁡θx = \rho\sin\varphi\cos\theta, y=ρsin⁡φsin⁡θy = \rho\sin\varphi\sin\theta, z=ρcos⁡φz = \rho\cos\varphi the same computation gives dx∧dy∧dz=det⁡∂(x,y,z)∂(ρ,φ,θ) dρ∧dφ∧dθ=ρ2sin⁡φ dρ∧dφ∧dθdx\wedge dy\wedge dz = \det\dfrac{\partial(x, y, z)}{\partial(\rho, \varphi, \theta)}\,d\rho\wedge d\varphi\wedge d\theta = \rho^2\sin\varphi\,d\rho\wedge d\varphi\wedge d\theta.

    Answer. dx∧dy=r dr∧dθdx\wedge dy = r\,dr\wedge d\theta and dx∧dy∧dz=ρ2sin⁡φ dρ∧dφ∧dθdx\wedge dy\wedge dz = \rho^2\sin\varphi\,d\rho\wedge d\varphi\wedge d\theta. In general the pullback of dx1∧⋯∧dxndx_1\wedge\dots\wedge dx_n under a map is the Jacobian determinant times du1∧⋯∧dundu_1\wedge\dots\wedge du_n. That is the change of variables formula of 6.4: the ∣J∣|J| factor is the coefficient of the pullback, and the absolute value appears when you turn an oriented integral into an area integral for a map that preserves orientation.

    Your answerPull dx∧dydx\wedge dy back under x=u2−vx = u^2 - v, y=u+v2y = u + v^2. The result is J du∧dvJ\,du\wedge dv. What is JJ at (u,v)=(1,2)(u, v) = (1, 2)?
    Proposition B Rudin ch. 10 and Spivak ch. 4

    Pullback commutes with d

    Let φ\varphi be a C2C^2 map and ω\omega a C1C^1 form. Then φ∗(dω)=d(φ∗ω)\varphi^*(d\omega) = d(\varphi^*\omega).

    In one lineFor a C2C^2 map φ\varphi and a form ω\omega, φ∗(dω)=d(φ∗ω)\varphi^*(d\omega) = d(\varphi^*\omega).

    Why it is trueIt is the chain rule for forms, and it is the step that lets a theorem proved on a rectangle of parameters hold on every parametrized surface.

    The case of a function ff is the chain rule: φ∗(df)=φ∗(fx dx+fy dy+… )=(fx xu+fy yu+… ) du+(… ) dv=d(f∘φ)\varphi^*(df) = \varphi^*(f_x\,dx + f_y\,dy + \dots) = (f_x\,x_u + f_y\,y_u + \dots)\,du + (\dots)\,dv = d(f\circ\varphi). Higher degrees follow because φ∗\varphi^* respects the wedge product and d(d(⋅))=0d(d(\cdot)) = 0. As a check, the form ω=x dy\omega = x\,dy and the map x=uvx = uv, y=u+vy = u + v have φ∗ω=uv (du+dv)\varphi^*\omega = uv\,(du + dv), so d(φ∗ω)=(v du+u dv)∧(du+dv)=(v−u) du∧dvd(\varphi^*\omega) = (v\,du + u\,dv)\wedge(du + dv) = (v - u)\,du\wedge dv after expanding (v du∧dv+u dv∧duv\,du\wedge dv + u\,dv\wedge du), while dω=dx∧dyd\omega = dx\wedge dy and φ∗(dx∧dy)=(v du+u dv)∧(du+dv)=(v−u) du∧dv\varphi^*(dx\wedge dy) = (v\,du + u\,dv)\wedge(du + dv) = (v - u)\,du\wedge dv. They agree.

    The one theorem

    Theorem B Spivak ch. 4 and Rudin 10.33; the special cases are 6.2, 6.4 and 6.6

    Generalized Stokes theorem

    Let MM be a compact oriented kk-dimensional region (an interval, a curve, a plane region, a surface, a solid) with piecewise smooth boundary ∂M\partial M oriented by the rule ‘outward normal first’: a basis of the boundary is positive if, with the outward normal in front, it is positive for MM. Let ω\omega be a C1C^1 form of degree k−1k - 1 on a neighbourhood of MM. Then ∫∂Mω=∫Mdω\int_{\partial M}\omega = \int_Md\omega.

    In one lineIf MM is a compact oriented kk-dimensional region with boundary ∂M\partial M carrying the induced orientation and ω\omega is a C1C^1 (k−1)(k-1)-form on a neighbourhood of MM, then ∫∂Mω=∫Mdω\displaystyle\int_{\partial M}\omega = \int_Md\omega.

    Why it is trueOne statement contains the fundamental theorem of calculus, the fundamental theorem for line integrals, Green's theorem, Stokes' theorem and the divergence theorem. It says the integral of a derivative over a region is the integral of the function over the boundary.

    The special cases are exactly the theorems of the unit. Each line is a theorem you have used.

    MMω\omegadωd\omega∫∂Mω=∫Mdω\int_{\partial M}\omega = \int_Md\omega is
    an interval [a,b][a, b], ∂M={b}−{a}\partial M = \{b\} - \{a\}a function fff′ dxf'\,dxf(b)−f(a)=∫abf′(x) dxf(b) - f(a) = \int_a^bf'(x)\,dx, the fundamental theorem of calculus
    a curve from AA to BBa function φ\varphidφ=∇φ⋅drd\varphi = \nabla\varphi\cdot drφ(B)−φ(A)=∫C∇φ⋅dr\varphi(B) - \varphi(A) = \int_C\nabla\varphi\cdot dr (6.2)
    a region DD in the planeP dx+Q dyP\,dx + Q\,dy(Qx−Py) dx∧dy(Q_x - P_y)\,dx\wedge dy∮∂DP dx+Q dy=∬D(Qx−Py) dA\oint_{\partial D}P\,dx + Q\,dy = \iint_D(Q_x - P_y)\,dA, Green (6.4)
    a surface SS in spaceP dx+Q dy+R dzP\,dx + Q\,dy + R\,dzthe 2-form of curl⁡F\operatorname{curl}F∮∂SF⋅dr=∬Scurl⁡F⋅n dS\oint_{\partial S}F\cdot dr = \iint_S\operatorname{curl}F\cdot n\,dS, Stokes (6.6)
    a solid EEA dy∧dz+B dz∧dx+C dx∧dyA\,dy\wedge dz + B\,dz\wedge dx + C\,dx\wedge dy(div⁡G) dx∧dy∧dz(\operatorname{div}G)\,dx\wedge dy\wedge dz∬∂EG⋅n dS=∭Ediv⁡G dV\iint_{\partial E}G\cdot n\,dS = \iiint_E\operatorname{div}G\,dV, Gauss (6.6)

    The orientation rule ‘outward normal first’ reproduces every rule of 6.4 and 6.6: counterclockwise for a plane region, the right-hand rule for a surface, outward for a solid.

    Why it is true

    ProofSpivak ch. 4
    The idea: Reduce to a cube. On the unit square, the boundary integral is four one-dimensional integrals and the area integral is two, and the fundamental theorem of calculus matches them. Then pull a general piece back to a cube with a parametrization, using that pullback commutes with d.
    ∫∂I2ω=∫01P(x,0) dx+∫01Q(1,y) dy−∫01P(x,1) dx−∫01Q(0,y) dy\displaystyle\int_{\partial I^2}\omega = \int_0^1P(x, 0)\,dx + \int_0^1Q(1, y)\,dy - \int_0^1P(x, 1)\,dx - \int_0^1Q(0, y)\,dy.
    The case of the unit square I2=[0,1]2I^2 = [0, 1]^2 and ω=P dx+Q dy\omega = P\,dx + Q\,dy: ∂I2\partial I^2 is the bottom (y=0y = 0, left to right), the right side (x=1x = 1, up), the top (right to left) and the left side (down).
    ∫I2dω=∫01 ⁣ ⁣[Q(1,y)−Q(0,y)]dy−∫01 ⁣ ⁣[P(x,1)−P(x,0)]dx\displaystyle\int_{I^2}d\omega = \int_0^1\!\!\big[Q(1, y) - Q(0, y)\big]dy - \int_0^1\!\!\big[P(x, 1) - P(x, 0)\big]dx.
    On the other side dω=(Qx−Py) dx∧dyd\omega = (Q_x - P_y)\,dx\wedge dy, and Fubini with the fundamental theorem of calculus in each variable.
    The two sides of the theorem agree on the square. The same argument, with one pair of opposite faces per coordinate, proves it on a kk-cube.
    Compare the two lines: the four terms are the same.In dimension kk the same happens: for each coordinate there are two opposite faces with opposite orientations, and the fundamental theorem of calculus in that coordinate pairs them. The sign pattern is the orientation rule.
    For a piece M=φ(Ik)M = \varphi(I^k): ∫∂Mω=∫∂Ikφ∗ω\int_{\partial M}\omega = \int_{\partial I^k}\varphi^*\omega.
    Definition of the integral over a parametrized piece φ:Ik→M\varphi: I^k\to M and its boundary: ∫∂φω=∫∂Ikφ∗ω\int_{\partial\varphi}\omega = \int_{\partial I^k}\varphi^*\omega.
    =∫Ikd(φ∗ω)= \int_{I^k}d(\varphi^*\omega).
    The cube case applied to the form φ∗ω\varphi^*\omega.
    =∫Ikφ∗(dω)= \int_{I^k}\varphi^*(d\omega).
    Pullback commutes with dd.
    =∫Mdω= \int_Md\omega. ■\blacksquare
    Definition of the integral of dωd\omega over MM.A general region is cut into finitely many such pieces. The integrals over the common faces cancel, because neighbouring pieces traverse a shared face in opposite directions, as in the 'why the inner edges cancel' box of 6.6.
    Proof workshopGive the reason for each line: from the cube to a parametrized pieceSpivak ch. 4

    Each line of the proof is shown in order. Choose the reason that justifies it, then press Check.

    1. ∫∂Mω=∫∂Ikφ∗ω\displaystyle\int_{\partial M}\omega = \int_{\partial I^k}\varphi^*\omega.
    2. =∫Ikd(φ∗ω)\displaystyle= \int_{I^k}d(\varphi^*\omega).
    3. =∫Ikφ∗(dω)\displaystyle= \int_{I^k}\varphi^*(d\omega).
    4. =∫Mdω\displaystyle= \int_Md\omega.
    Worked example new (the proof of Stokes via forms)
    Recover Stokes' theorem of 6.6 from Green's theorem and pullback.

    Solution.

    1. Setup. Let r:T→Sr: T\to S be a smooth parametrization on a rectangle TT, with ru×rvr_u\times r_v pointing the way the normal nn does (so that ∂T\partial T run counterclockwise maps to ∂S\partial S with the induced orientation), and ω=P dx+Q dy+R dz\omega = P\,dx + Q\,dy + R\,dz the 1-form of FF.
    2. Boundary side. ∮∂SF⋅dr=∫∂Tr∗ω\oint_{\partial S}F\cdot dr = \int_{\partial T}r^*\omega, since a line integral is the integral of the pullback.
    3. Green on the rectangle. ∫∂Tr∗ω=∬Td(r∗ω)\int_{\partial T}r^*\omega = \iint_Td(r^*\omega). Green’s theorem for the form r∗ω=p du+q dvr^*\omega = p\,du + q\,dv, which is a planar 1-form on TT.
    4. Pullback commutes with dd. d(r∗ω)=r∗(dω)d(r^*\omega) = r^*(d\omega).
    5. The 2-form dωd\omega is the curl. dω=(Ry−Qz) dy∧dz+…d\omega = (R_y - Q_z)\,dy\wedge dz + \dots, and its pullback is (curl⁡F)(r)⋅(ru×rv) du∧dv(\operatorname{curl}F)(r)\cdot(r_u\times r_v)\,du\wedge dv by the proposition above.

    Answer. ∮∂SF⋅dr=∬T(curl⁡F)(r)⋅(ru×rv) du dv=∬S(curl⁡F)⋅n dS\oint_{\partial S}F\cdot dr = \iint_T(\operatorname{curl}F)(r)\cdot(r_u\times r_v)\,du\,dv = \iint_S(\operatorname{curl}F)\cdot n\,dS: Stokes’ theorem. The proof in 6.6 for a graph was this argument with the chain rule written out.

    Worked example new (compare 6.6 and 6.4)
    Verify the generalized Stokes theorem for ω=x dy∧dz+y dz∧dx+z dx∧dy\omega = x\,dy\wedge dz + y\,dz\wedge dx + z\,dx\wedge dy on the ball of radius RR, and for ω=12(x dy−y dx)\omega = \frac12(x\,dy - y\,dx) on the ellipse x2a2+y2b2≤1\frac{x^2}{a^2} + \frac{y^2}{b^2}\le1.

    Solution.

    1. Ball, the sphere side. r∗ω=F(r)⋅(rφ×rθ) dφ∧dθr^*\omega = F(r)\cdot(r_\varphi\times r_\theta)\,d\varphi\wedge d\theta with F=(x,y,z)=Rρ^F = (x, y, z) = R\hat\rho and rφ×rθ=R2sin⁡φ ρ^r_\varphi\times r_\theta = R^2\sin\varphi\,\hat\rho outward, so F⋅N=R3sin⁡φF\cdot N = R^3\sin\varphi and ∬=2π⋅2⋅R3=4πR3\iint = 2\pi\cdot2\cdot R^3 = 4\pi R^3.
    2. Ball, the volume side. dω=(1+1+1) dx∧dy∧dz=3 dx∧dy∧dzd\omega = (1 + 1 + 1)\,dx\wedge dy\wedge dz = 3\,dx\wedge dy\wedge dz, and ∭3 dV=3⋅43πR3=4πR3\iiint3\,dV = 3\cdot\frac43\pi R^3 = 4\pi R^3. They agree.
    3. Ellipse, the volume (area) side. dω=12(dx∧dy−dy∧dx)=dx∧dyd\omega = \frac12(dx\wedge dy - dy\wedge dx) = dx\wedge dy, and ∬dx dy=πab\iint dx\,dy = \pi ab.
    4. Ellipse, the boundary side. With x=acos⁡tx = a\cos t, y=bsin⁡ty = b\sin t: r∗ω=12(acos⁡t⋅bcos⁡t+bsin⁡t⋅asin⁡t) dt=12ab dtr^*\omega = \frac12(a\cos t\cdot b\cos t + b\sin t\cdot a\sin t)\,dt = \frac12ab\,dt and ∫02π=πab\int_0^{2\pi} = \pi ab. They agree.

    Answer. 4πR34\pi R^3 and πab\pi ab. The second is the area-by-line-integral formula of 6.4; the first is the volume of the ball as a third of the flux of the position vector.

    Your answerUse the generalized Stokes theorem to compute ∬Sz dx∧dy\iint_Sz\,dx\wedge dy over the unit sphere SS with the outward orientation.
    Quick checkIn the generalized Stokes theorem for a solid EE in R3\mathbb R^3, what degree is ω\omega?
    Quick checkFor M=[a,b]M = [a, b] and ω=f\omega = f a function, what is ∫∂Mω\int_{\partial M}\omega?
    Quick checkWhy can you not apply the generalized Stokes theorem to the angle form dθd\theta on the unit disc?
    Python labStokes numerically on a wavy patch, and the sign flip when the parameters are swapped
    import numpy as np
    
    def F(x, y, z):     return (y * z, x**2, np.sin(x) + y * z)            # P, Q, R
    def curlF(x, y, z): return (z, y - np.cos(x), 2 * x - z)               # (R_y - Q_z, P_z - R_x, Q_x - P_y), by hand
    def surf(u, v):     return (u, v, 0.3 * np.sin(2 * u) * np.cos(3 * v) + 0.2 * u * v)
    def d(f, u, v, i, h=1e-6):                                              # partial derivative of the surface, finite differences
        return (np.array(f(u + h * (i == 0), v + h * (i == 1))) - np.array(f(u - h * (i == 0), v - h * (i == 1)))) / (2 * h)
    
    def both_sides(r, u0, u1, v0, v1, n=400):
        u = u0 + (np.arange(n) + .5) * (u1 - u0) / n; v = v0 + (np.arange(n) + .5) * (v1 - v0) / n
        U, V = np.meshgrid(u, v, indexing="ij")
        N = np.cross(d(r, U, V, 0), d(r, U, V, 1), axis=0)                  # r_u x r_v
        c = curlF(*r(U, V)); surface = sum(ci * Ni for ci, Ni in zip(c, N)).sum() * (u1 - u0) * (v1 - v0) / n**2
        def edge(a, t0, t1, along):                  # pull back F . dr to one edge: along = 0, u moves; along = 1, v moves
            t = t0 + (np.arange(n) + .5) * (t1 - t0) / n
            uu, vv = (t, a + 0 * t) if along == 0 else (a + 0 * t, t)
            rr = d(r, uu, vv, along); f = F(*r(uu, vv)); return sum(fi * ri for fi, ri in zip(f, rr)).sum() * (t1 - t0) / n
        line = edge(v0, u0, u1, 0) + edge(u1, v0, v1, 1) - edge(v1, u0, u1, 0) - edge(u0, v0, v1, 1)
        return round(line, 5), round(surface, 5)
    
    print("r(u, v):  (circulation, flux of curl) =", both_sides(surf, 0, 1, 0, 2))
    print("r(v, u):  (circulation, flux of curl) =", both_sides(lambda u, v: surf(v, u), 0, 2, 0, 1))

    Try this. Run it. The surface is the graph z = 0.3 sin(2u) cos(3v) + 0.2uv over the rectangle [0,1] by [0,2], and the field is F = (yz, x^2, sin x + yz). The first number is the line integral of F round the edge of the patch, computed by pulling F . dr back to the four edges of the parameter rectangle. The second is the integral of curl F . (r_u x r_v) over the rectangle, from the cross product of finite-difference tangent vectors. Both are 1.2596: Stokes' theorem on a surface with no symmetry. The second line swaps the parameters (u and v exchange roles, the rectangle becomes [0,2] by [0,1]): the same surface, but both numbers become -1.2596, because the orientation reversed. Change the surface, the field (update the hand-computed curl, P_z, R_x and so on) or the rectangle, and check that the two numbers agree every time.

    Where it is used

    Used inside maths· Inside maths: counting holes, and complex analysis

    The angle form dθd\theta is closed and not exact on the punctured plane. The set of closed 1-forms modulo exact ones is a vector space, the first de Rham cohomology group, and its dimension counts the independent holes of the region: 11 for the punctured plane, 00 for the plane. In complex analysis, dz/z=dr/r+i dθdz/z = dr/r + i\,d\theta and ∮∣z∣=1dz/z=i∫02πdθ=2πi\oint_{\lvert z\rvert=1}dz/z = i\int_0^{2\pi}d\theta = 2\pi i because the integral of the exact part dr/r=d(ln⁡r)dr/r = d(\ln r) around a closed path is 00 and the angle form contributes 2π2\pi. That single computation is the heart of the residue theorem. In the other direction, the divergence-free field G=(x,y,z)/ρ3G = (x, y, z)/\rho^3 of 6.6 is closed and not exact as a 2-form on R3∖{0}\mathbb R^3\setminus\{0\}: its integral over the unit sphere is 4π4\pi, whereas the integral of an exact 2-form over a closed surface is 00.

    Used in the real world· Used in the real world: heat is not exact, entropy is

    For one mole of an ideal gas, the heat absorbed in a small change of temperature TT and volume VV is the 1-form δQ=CV dT+RTV dV\delta Q = C_V\,dT + \dfrac{RT}V\,dV, with R=8.314R = 8.314 J/(mol K) and CV=32RC_V = \frac32R for a monatomic gas. Is δQ\delta Q closed? d(δQ)=RV dT∧dV≠0d(\delta Q) = \dfrac RV\,dT\wedge dV\neq0 (the coefficient of dVdV has a TT-derivative R/VR/V, and the coefficient of dTdT has none). So heat is not exact: there is no function Q(T,V)Q(T, V) whose change is the heat absorbed, and the heat depends on the path, as the line integral of a non-closed form does. Divide by TT: δQT=CVdTT+RdVV\dfrac{\delta Q}T = C_V\dfrac{dT}T + R\dfrac{dV}V has d(δQT)=0d\big(\dfrac{\delta Q}T\big) = 0 because the coefficient of dVdV does not depend on TT and that of dTdT does not depend on VV. It is exact: δQT=d(CVln⁡T+Rln⁡V)\dfrac{\delta Q}T = d(C_V\ln T + R\ln V). That function is the entropy SS, and 1/T1/T is an integrating factor in the sense of the exact differential equations of 6.2. From 300300 K and 11 litre to 600600 K and 22 litres, with CV=32RC_V = \frac32R, the change is ΔS=(32R+R)ln⁡2=2.5⋅8.314⋅0.693≈14.4\Delta S = (\frac32R + R)\ln2 = 2.5\cdot8.314\cdot0.693\approx14.4 J/K, along any path. A heat engine designer uses this: the area enclosed in a temperature-entropy diagram is the net heat of a cycle.

    Your answerOne mole of a monatomic ideal gas (CV=32RC_V = \frac32R, R=8.314R = 8.314 J/(mol K)) goes from (T,V)=(300 K,1 L)(T, V) = (300\text{ K}, 1\text{ L}) to (600 K,2 L)(600\text{ K}, 2\text{ L}). The form δQ/T\delta Q/T is exact. What is the entropy change in J/K?
    A preview· Preview: Maxwell's equations in two lines

    In 6.8 the four Maxwell equations are written with ∇⋅\nabla\cdot and ∇×\nabla\times. In the language of this chapter the electric and magnetic fields together make up one 2-form FF in four-dimensional spacetime, and two of the equations (no magnetic monopoles, and Faraday’s law) say dF=0dF = 0, while the other two (Gauss’s law and Ampere’s law with Maxwell’s correction) say that dd of the Hodge star of FF equals the current 3-form. That dF=0dF = 0 implies that F=dAF = dA locally, with AA the four-potential, by the Poincare lemma. The integral forms of the equations in 6.8 are the generalized Stokes theorem for those forms.

    Traps

    Watch out· The wedge is not a product of numbers

    dx∧dy=−dy∧dxdx\wedge dy = -dy\wedge dx and dx∧dx=0dx\wedge dx = 0. A form dx∧dydx\wedge dy is not the same as the area element dx dydx\,dy: the first has an orientation and the second does not. Integrating dy∧dxdy\wedge dx over a region gives minus its area, unless you swap the orientation of the region at the same time. Reordering without the sign is the most common error.

    Watch out· Substitute dx too, not only x

    To pull back, replace the coefficients and also each differential: dx=xu du+xv dvdx = x_u\,du + x_v\,dv. Leaving dxdx unchanged, or replacing it by x′(t)x'(t) without the dtdt, is the form version of the slip of dotting F(x,y)F(x, y) with r′(t)r'(t) before substituting the path (see the traps of 6.1). The method of 6.1 (parametrize, differentiate, substitute, dot, integrate) is the pullback.

    Watch out· Closed does not mean exact, and d squared needs C two

    dω=0d\omega = 0 is a local test; ω=df\omega = df is global. The angle form is closed everywhere on the punctured plane and exact nowhere around the hole. Also d2=0d^2 = 0 needs continuous second partials: it is the equality of mixed partials in disguise.

    Watch out· Notation of components

    On the 376 midterm, a student wrote the components of FF as Fx,Fy,FzF_x, F_y, F_z, which invites confusion with partial derivatives. In the language of forms the coefficients are P,Q,RP, Q, R (or F1,F2,F3F_1, F_2, F_3), and a subscript xx always means a derivative. Write the component names first and the derivatives after.

    Problems

    check

    Compute dωd\omega for ω=x2y dx+xy2 dy\omega = x^2y\,dx + xy^2\,dy.

    Solution

    P=x2yP = x^2y, Q=xy2Q = xy^2. dω=(Qx−Py) dx∧dy=(y2−x2) dx∧dyd\omega = (Q_x - P_y)\,dx\wedge dy = (y^2 - x^2)\,dx\wedge dy.

    check

    Compute (dx+dy)∧(dx−dy)(dx + dy)\wedge(dx - dy) and (x dy+y dx)∧dx(x\,dy + y\,dx)\wedge dx.

    Solution

    (dx+dy)∧(dx−dy)=dx∧dx−dx∧dy+dy∧dx−dy∧dy=−dx∧dy−dx∧dy=−2 dx∧dy(dx + dy)\wedge(dx - dy) = dx\wedge dx - dx\wedge dy + dy\wedge dx - dy\wedge dy = -dx\wedge dy - dx\wedge dy = -2\,dx\wedge dy. And (x dy+y dx)∧dx=x dy∧dx+y dx∧dx=−x dx∧dy(x\,dy + y\,dx)\wedge dx = x\,dy\wedge dx + y\,dx\wedge dx = -x\,dx\wedge dy.

    checknew (the pullback of 6.1)

    Compute the pullback of ω=y dx+x dy\omega = y\,dx + x\,dy under r(t)=(t2,t3)r(t) = (t^2, t^3) and integrate from t=0t = 0 to 11. Compare with the potential.

    Solution

    dx=2t dtdx = 2t\,dt, dy=3t2 dtdy = 3t^2\,dt, so r∗ω=t3⋅2t dt+t2⋅3t2 dt=5t4 dtr^*\omega = t^3\cdot2t\,dt + t^2\cdot3t^2\,dt = 5t^4\,dt and the integral is 11. The potential is f=xyf = xy with f(1,1)−f(0,0)=1f(1, 1) - f(0, 0) = 1.

    check

    Which of these 1-forms on R2\mathbb R^2 is closed: y dx+x dyy\,dx + x\,dy, y dx−x dyy\,dx - x\,dy, x dx+y dyx\,dx + y\,dy?

    Solution

    dd of P dx+Q dyP\,dx + Q\,dy is (Qx−Py) dx∧dy(Q_x - P_y)\,dx\wedge dy: first 1−1=01 - 1 = 0 (closed, =d(xy)= d(xy)); second −1−1=−2-1 - 1 = -2 (not closed); third 0−0=00 - 0 = 0 (closed, =dx2+y22= d\frac{x^2 + y^2}2).

    usenew (Final 3 of MATH 376 in the language of forms)

    Show that ω=yzcos⁡(xyz) dx+xzcos⁡(xyz) dy+xycos⁡(xyz) dz\omega = yz\cos(xyz)\,dx + xz\cos(xyz)\,dy + xy\cos(xyz)\,dz is exact, and find ∫Cω\int_C\omega for CC the unit circle in the plane z=0z = 0.

    Solution

    ω=d(sin⁡(xyz))\omega = d\big(\sin(xyz)\big) by the chain rule: ∂xsin⁡(xyz)=yzcos⁡(xyz)\partial_x\sin(xyz) = yz\cos(xyz), and the same for yy, zz. So ω\omega is exact and the integral around a closed curve is 00 (the function returns to its value). Also dω=d(dsin⁡(xyz))=0d\omega = d(d\sin(xyz)) = 0 by d2=0d^2 = 0: this is the curl computed term by term in 6.6, in one line.

    usenew (generalized Stokes, a surface integral by boundary)

    Compute ∬Sx dy∧dz\iint_Sx\,dy\wedge dz where SS is the unit sphere with the outward orientation.

    Solution

    ω=x dy∧dz\omega = x\,dy\wedge dz has dω=dx∧dy∧dzd\omega = dx\wedge dy\wedge dz. By the generalized Stokes theorem with MM the unit ball, ∬Sω=∭Bdω=Vol⁡(B)=4π3\iint_S\omega = \iiint_Bd\omega = \operatorname{Vol}(B) = \frac{4\pi}3. (Equivalently the flux of (x,0,0)(x, 0, 0), with divergence 11.)

    usenew (pullback with a Jacobian)

    Pull dx∧dydx\wedge dy back under x=2u+vx = 2u + v, y=u−3vy = u - 3v and compute the area of the image of the square [0,1]2[0, 1]^2.

    Solution

    dx=2 du+dvdx = 2\,du + dv, dy=du−3 dvdy = du - 3\,dv. dx∧dy=2⋅(−3) du∧dv+1⋅1 dv∧du=(−6−1) du∧dv=−7 du∧dvdx\wedge dy = 2\cdot(-3)\,du\wedge dv + 1\cdot1\,dv\wedge du = (-6 - 1)\,du\wedge dv = -7\,du\wedge dv. This is det⁡[211−3]=−7\det\begin{bmatrix}2&1\\1&-3\end{bmatrix} = -7. The map reverses orientation, and the area is ∣−7∣⋅1=7\lvert-7\rvert\cdot1 = 7.

    usenew (thermodynamics, 6.7)

    Show that ω=p dV+V dp\omega = p\,dV + V\,dp is exact and ω′=p dV\omega' = p\,dV is not closed on the (V,p)(V, p)-plane. What does the second say about the work ∫p dV\int p\,dV done by a gas?

    Solution

    Write ω=p dV+V dp\omega = p\,dV + V\,dp: in coordinates (V,p)(V, p) with the form P dV+Q dpP\,dV + Q\,dp, P=pP = p, Q=VQ = V, dω=(QV−Pp) dV∧dp=(1−1) dV∧dp=0d\omega = (Q_V - P_p)\,dV\wedge dp = (1 - 1)\,dV\wedge dp = 0, and ω=d(pV)\omega = d(pV). For ω′=p dV\omega' = p\,dV: P=pP = p, Q=0Q = 0, dω′=(0−1) dV∧dp=−dV∧dp≠0d\omega' = (0 - 1)\,dV\wedge dp = -dV\wedge dp\neq0. So the work ∫p dV\int p\,dV depends on the path in the (V,p)(V, p)-plane: for a cycle it equals the area enclosed, which is why a heat engine can do net work.

    provenew (rule for d of a product)

    Prove that for a function ff and a 1-form ω\omega, d(fω)=df∧ω+f dωd(f\omega) = df\wedge\omega + f\,d\omega.

    Solution

    Write ω=∑iPi dxi\omega = \sum_iP_i\,dx_i. Then fω=∑ifPi dxif\omega = \sum_ifP_i\,dx_i, and by the definition d(fPi dxi)=d(fPi)∧dxi=(f dPi+Pi df)∧dxid(fP_i\,dx_i) = d(fP_i)\wedge dx_i = (f\,dP_i + P_i\,df)\wedge dx_i by the product rule for functions. Sum over ii: d(fω)=f∑idPi∧dxi+df∧∑iPi dxi=f dω+df∧ωd(f\omega) = f\sum_idP_i\wedge dx_i + df\wedge\sum_iP_i\,dx_i = f\,d\omega + df\wedge\omega.

    provenew (exact implies closed, with the line integral)

    Prove that if ω=P dx+Q dy\omega = P\,dx + Q\,dy is exact with C2C^2 potential ff, then ∮Cω=0\oint_C\omega = 0 for every closed piecewise smooth path, and dω=0d\omega = 0.

    Solution

    ω=df\omega = df means P=fxP = f_x and Q=fyQ = f_y. Then dω=(Qx−Py) dx∧dy=(fyx−fxy) dx∧dy=0d\omega = (Q_x - P_y)\,dx\wedge dy = (f_{yx} - f_{xy})\,dx\wedge dy = 0 by equality of mixed partials. For a closed path rr on [a,b][a, b], ∮ω=∫ab(fxx′+fyy′) dt=∫abddtf(r(t)) dt=f(r(b))−f(r(a))=0\oint\omega = \int_a^b(f_xx' + f_yy')\,dt = \int_a^b\frac{d}{dt}f(r(t))\,dt = f(r(b)) - f(r(a)) = 0, by the chain rule and the fundamental theorem of calculus.

    provenew (determinant)

    Prove that for a linear map AA of the plane, A∗(dx∧dy)=det⁡A dx∧dyA^*(dx\wedge dy) = \det A\,dx\wedge dy. Deduce that a linear map multiplies areas by ∣det⁡A∣\lvert\det A\rvert.

    Solution

    Let A=[abcd]A = \begin{bmatrix}a&b\\c&d\end{bmatrix}, so the new coordinates are x′=au+bvx' = au + bv, y′=cu+dvy' = cu + dv in terms of old ones. dx′=a du+b dvdx' = a\,du + b\,dv and dy′=c du+d dvdy' = c\,du + d\,dv. Then dx′∧dy′=ad du∧dv+bc dv∧du=(ad−bc) du∧dvdx'\wedge dy' = ad\,du\wedge dv + bc\,dv\wedge du = (ad - bc)\,du\wedge dv. The area of the image of a region DD is ∬D∣det⁡A∣ du dv=∣det⁡A∣area⁡(D)\iint_D\lvert\det A\rvert\,du\,dv = \lvert\det A\rvert\operatorname{area}(D), since the area counts regardless of orientation.

    stretchnew (the Poincare lemma for 1-forms)

    Let ω=∑iPi dxi\omega = \sum_iP_i\,dx_i be a closed C1C^1 1-form on a ball centred at 00 in Rn\mathbb R^n, meaning ∂jPi=∂iPj\partial_jP_i = \partial_iP_j. Prove that f(x)=∫01∑iPi(tx)xi dtf(x) = \int_0^1\sum_iP_i(tx)x_i\,dt satisfies df=ωdf = \omega.

    Solution

    Differentiate under the integral sign: ∂jf=∫01[Pj(tx)+∑itxi(∂jPi)(tx)]dt\partial_jf = \int_0^1\Big[P_j(tx) + \sum_itx_i(\partial_jP_i)(tx)\Big]dt. By closedness, ∂jPi=∂iPj\partial_jP_i = \partial_iP_j, so the sum is ∑itxi(∂iPj)(tx)=t ddt[Pj(tx)]\sum_itx_i(\partial_iP_j)(tx) = t\,\frac d{dt}\big[P_j(tx)\big] (chain rule: ddtPj(tx)=∑ixi(∂iPj)(tx)\frac d{dt}P_j(tx) = \sum_ix_i(\partial_iP_j)(tx)). So ∂jf=∫01[Pj(tx)+tddtPj(tx)]dt=∫01ddt[tPj(tx)]dt=Pj(x)\partial_jf = \int_0^1\Big[P_j(tx) + t\frac d{dt}P_j(tx)\Big]dt = \int_0^1\frac d{dt}\big[tP_j(tx)\big]dt = P_j(x). Hence df=ωdf = \omega. The integral is the line integral of ω\omega along the ray from 00 to xx, the method of 6.2.

    stretchnew (the cube proof in three variables)

    Verify the generalized Stokes theorem on the unit cube [0,1]3[0,1]^3 for ω=f dy∧dz\omega = f\,dy\wedge dz, and say what the other two kinds of term do.

    Solution

    dω=fx dx∧dy∧dzd\omega = f_x\,dx\wedge dy\wedge dz, so ∫[0,1]3dω=∬[f(1,y,z)−f(0,y,z)]dy dz\int_{[0,1]^3}d\omega = \iint\big[f(1, y, z) - f(0, y, z)\big]dy\,dz by Fubini and the fundamental theorem of calculus in xx. On the boundary, dy∧dzdy\wedge dz restricted to the faces y=y = const or z=z = const is 00 (one of dydy, dzdz is zero there), so only the faces x=1x = 1 (outward normal +x+x, orientation ++) and x=0x = 0 (outward −x-x, orientation −-) contribute, giving ∬f(1,y,z) dy dz−∬f(0,y,z) dy dz\iint f(1, y, z)\,dy\,dz - \iint f(0, y, z)\,dy\,dz, the same. The terms g dz∧dxg\,dz\wedge dx and h dx∧dyh\,dx\wedge dy are the same with the roles of the coordinates cycled, and adding the three gives the divergence theorem on the cube.

    Summary

    • A differential form is the thing under an integral sign: it eats tiny oriented pieces (vectors, parallelograms, boxes) and returns numbers. 1-forms are line-integral integrands P dx+Q dy+R dzP\,dx + Q\,dy + R\,dz (vector fields), 2-forms are flux integrands (vector fields again), 3-forms are volume integrands.
    • The wedge product is bilinear with dxi∧dxj=−dxj∧dxidx_i\wedge dx_j = -dx_j\wedge dx_i and dxi∧dxi=0dx_i\wedge dx_i = 0. On 1-forms in space it is the cross product; in the plane it is the determinant.
    • The exterior derivative dd raises the degree by one: gradient on functions, curl on 1-forms, divergence on 2-forms. In the plane d(P dx+Q dy)=(Qx−Py) dx∧dyd(P\,dx + Q\,dy) = (Q_x - P_y)\,dx\wedge dy.
    • d2=0d^2 = 0, because mixed partials are equal and the wedge is anti-symmetric. This is both curl⁡∇f=0\operatorname{curl}\nabla f = 0 and div⁡curl⁡F=0\operatorname{div}\operatorname{curl}F = 0.
    • A form is closed if dω=0d\omega = 0 and exact if ω=dη\omega = d\eta. Exact implies closed; the converse holds on star-shaped regions (Poincare lemma) and fails around holes: the angle form dθd\theta.
    • The pullback substitutes a parametrization, with dx=xu du+xv dvdx = x_u\,du + x_v\,dv. It turns forms into the integrals of 6.1 and 6.5, the cross product ru×rvr_u\times r_v appears as the coefficient of du∧dvdu\wedge dv, and the Jacobian determinant is the pullback of dx∧dydx\wedge dy.
    • Generalized Stokes: ∫∂Mω=∫Mdω\int_{\partial M}\omega = \int_Md\omega. It is proved on a cube by the fundamental theorem of calculus and extended by pullback. The cases are the fundamental theorem of calculus, the fundamental theorem for line integrals, Green, Stokes and Gauss.

    Watch and read

    The same ideas explained by someone else. Use these when an explanation here does not land.

    Where this sits

    Comes from
    6.1 Paths and line integrals
    6.2 Gradient fields and potentials
    6.4 Green's theorem and change of variables
    6.5 Surfaces, area and flux
    6.6 Curl, divergence, Stokes' theorem and Gauss's theorem
    Leads to
    6.8 Maxwell's equations, conservation laws and heat
    Same idea elsewhere
    MATH 376: Only the notation is 376: the symbol P dx + Q dy as a line integral (6.1) and as an 'exact differential equation' when P_y = Q_x (6.2, Midterm 1 problem 1). This chapter says what that expression is: a 1-form, and the exactness test is d omega = 0. Everything else here is label B, beyond the course.
    MATH 561 and 562: Differential forms are the integrands on manifolds. The generalized Stokes theorem of this chapter, proved on manifolds with boundary, is the central theorem of the course, and de Rham cohomology (closed forms modulo exact forms) is how the course counts holes.
    MATH 521: On the plane or in space the exterior derivative of a 1-form is the curl and of a 2-form the divergence, so d squared = 0 is the pair of identities curl grad = 0 and div curl = 0 of 6.6, proved from the equality of mixed partials in 5.5.
    Physics and geometry: Maxwell's four equations become two, dF = 0 and d of the star of F = J, for one 2-form F. Thermodynamics asks which differential forms are exact: heat is not, but heat divided by temperature is, and that is the entropy.