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 ∫CPdx+Qdy and in 6.2 you called Pdx+Qdy an exact differential. We never said what dx 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 d turns gradient, curl and divergence into one rule; the identities curl grad =0 and div curl =0 become one identity, d2=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 and divcurlF=0. 6.5: flux through a parametrized surface is ∬F⋅(ru×rv)dudv. 6.4: the Jacobian determinant is the area stretching factor. 6.1: ∫CPdx+Qdy=∫(Px′+Qy′)dt.
Quick checkFor a path r(t)=(x(t),y(t)), the symbol ∫CPdx+Qdy means which integral?
dx stands for x′(t)dt and dy for y′(t)dt. This is the definition of the line integral in 6.1, and the pullback of this chapter makes the substitution a rule.Review: the line integral
Quick checkWhat is curl(∇f) for a C2 function f on R3?
Each component of the curl is a difference of two equal mixed partials. This chapter proves it in one line.Review: curl grad = 0
Quick checkA map sends the parameter square with area A in the uv-plane to a region in the xy-plane. Its Jacobian determinant is 3 everywhere. The area of the region is
The Jacobian determinant is the factor by which the map stretches areas. This is the change of variables formula of 6.4, which this chapter explains as a rule for the wedge product.Review: the Jacobian as an area factor
The picture: what a form eats
Think about what an integral sign needs. In ∫CF⋅dr it needs, at each point of the curve, a small displacement vector, and it returns a number: the part of F along the displacement. In ∬SF⋅ndS it needs a tiny piece of surface, a tiny parallelogram with two edge vectors a and b in a definite order, and it returns a number: the amount of F crossing the piece, F⋅(a×b). In ∭fdV it needs a tiny parallelepiped with three edges, and returns f times the signed volume 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.
Degree
Typical form
Eats
What it returns
Same as
0
f
a point
f at the point
a scalar field
1
Pdx+Qdy+Rdz
a vector v
Pv1+Qv2+Rv3=F⋅v
a vector field F=(P,Q,R), via the line integral
2
Ady∧dz+Bdz∧dx+Cdx∧dy
two vectors a,b
(A,B,C)⋅(a×b)
a vector field G=(A,B,C), via the flux
3
fdx∧dy∧dz
three vectors a,b,c
fdet[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, d, acting on degrees 0→1, 1→2, 2→3.
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=−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 a1dx1+⋯+andxn, 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 p, dx is the linear function that sends a vector v=(v1,v2,v3) to v1, dy sends it to v2, and dz to v3. A 1-form is ω=Pdx+Qdy+Rdz, where P,Q,R are functions of position: at p it sends v to P(p)v1+Q(p)v2+R(p)v3=F(p)⋅v, with F=(P,Q,R). So 1-forms and vector fields carry the same information. The line integral of the vector field F along a path is the integral of the 1-form ωF=F⋅dr.
In one lineIn R3 a 1-form is an expression ω=Pdx+Qdy+Rdz with functions P,Q,R; at each point it takes a vector v to the number Pv1+Qv2+Rv3.
Why it existsIt is the thing under a line integral. dx is the 1-form that reads off the x component of a displacement, so Pdx+Qdy is 'the part of the field along the step'.
Definition B
Two-form
The symbol dx∧dy (’dx wedge dy’) is the 2-form that sends a pair of vectors (a,b) to a1b2−a2b1, the signed area of the shadow of the parallelogram on the xy-plane. Likewise dy∧dz sends (a,b) to a2b3−a3b2 and dz∧dx to a3b1−a1b3. These are exactly the three components of a×b. A 2-form is a combination ω=Ady∧dz+Bdz∧dx+Cdx∧dy, and it sends (a,b) to (A,B,C)⋅(a×b). The flux of G=(A,B,C) through a surface is the integral of the 2-form Ady∧dz+Bdz∧dx+Cdx∧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 a 2-form is Ady∧dz+Bdz∧dx+Cdx∧dy; at each point it takes two vectors a,b to (A,B,C)⋅(a×b), the amount of the field (A,B,C) crossing the parallelogram on a and b.
Why it existsIt is the thing under a surface integral: it eats a tiny oriented parallelogram. The sign changes when you swap a and b, which is what orientation means.
Definition B
Differential form
A differential form of degree k (a k-form) on a region of Rn is a sum ω=∑fi1…ikdxi1∧⋯∧dxik over increasing index lists, with functions fi1…ik. At each point it eats k vectors and returns a number, linear in each vector and changing sign when two vectors are swapped. In R3 the degrees are 0,1,2,3; in R2 they are 0,1,2 with the only 2-form fdx∧dy; there are no k-forms for k larger than n.
In one lineA k-form is a rule that assigns to each point a function of k 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∧⋯∧dxik. A 0-form is a function.
Why it existsIt is the single object that can be integrated over a k-dimensional piece, for every k: points, curves, surfaces, solids.
Definition B
Wedge product
The wedge productω∧η of a j-form and a k-form is a (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∧dx, and hence dx∧dx=0 (since dx∧dx=−dx∧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.
In one lineThe product of forms, bilinear and associative, with ω∧η=(−1)jkη∧ω for degrees j,k; in particular dxi∧dxj=−dxj∧dxi and dxi∧dxi=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 j and k combine into degree j+k.
Everything about the wedge product follows from dxi∧dxj=−dxj∧dxi. Expand, kill the terms with a repeated dxi, and reorder with a sign.
Worked example
Compute (2dx−3dy)∧(dx+4dy) in the plane, and (dx+2dy+3dz)∧(dy+dz) in space. Compare the second with a cross product.
Solution.
Plane, expand.2dx∧dx+8dx∧dy−3dy∧dx−12dy∧dy. Distribute: four products.
Plane, simplify.dx∧dx=0=dy∧dy and −3dy∧dx=+3dx∧dy, so the result is (8+3)dx∧dy=11dx∧dy. This is the determinant det[21−34]=11: the wedge of two 1-forms in the plane is a 2×2 determinant.
Space, simplify. Drop dy∧dy and dz∧dz. Rewrite dx∧dz=−dz∧dx and 3dz∧dy=−3dy∧dz. The result is (2−3)dy∧dz−dz∧dx+dx∧dy=−dy∧dz−dz∧dx+dx∧dy.
Compare. In the dictionary this is the vector (−1,−1,1), and (1,2,3)×(0,1,1)=(2−3,0−1,1−0)=(−1,−1,1).
Answer.11dx∧dy and −dy∧dz−dz∧dx+dx∧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∧dy?
dy∧dx=−dx∧dy, so the sum is 0. This sign rule is the whole content of the wedge product.
Quick checkWhat is (xdx)∧(xdx)?
Any 1-form wedged with a multiple of itself is zero, because dx∧dx=0. In general ω∧ω=0 for a 1-form.
The exterior derivative
Definition B
Exterior derivative
The exterior derivatived takes a k-form to a (k+1)-form. For a function f: df=fxdx+fydy+fzdz (the total differential of 5.3). For a form ω=fdxi1∧⋯∧dxik with one term: dω=df∧dxi1∧⋯∧dxik. For a sum, add. Because dxi∧dxi=0, only the derivatives of f in directions not already present survive.
In one lineThe operation d that raises the degree of a form by one: on a function df=fxdx+fydy+fzdz, and on fdxI it is df∧dxI, extended by linearity.
Why it existsIt is the one derivative that is built to be integrated: ∫Mdω over a region equals ∫∂Mω over its boundary. It contains gradient, curl and divergence.
Theorem BRudin 10.20 and the dictionary of vector calculus; Spivak ch. 4
d is gradient, curl and divergence
Under the dictionary of the table above, d on 0-forms is the gradient, on 1-forms the curl, and on 2-forms the divergence.
In one lineIn R3: df=fxdx+fydy+fzdz (gradient); d(Pdx+Qdy+Rdz)=(Ry−Qz)dy∧dz+(Pz−Rx)dz∧dx+(Qx−Py)dx∧dy (curl); d(Ady∧dz+Bdz∧dx+Cdx∧dy)=(Ax+By+Cz)dx∧dy∧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=fxdx+fydy+fzdz, the gradient ∇f=(fx,fy,fz) in the dictionary.
The 0-form case is the definition.
dω=dP∧dx+dQ∧dy+dR∧dz.
Apply the rule to each of the three terms of ω=Pdx+Qdy+Rdz.
Expand each differential: dP=Pxdx+Pydy+Pzdz, and so on.
=(Ry−Qz)dy∧dz+(Pz−Rx)dz∧dx+(Qx−Py)dx∧dy.
Terms with a repeated dxi vanished already (Pxdx∧dx=0 and so on). Reorder with dy∧dx=−dx∧dy, dz∧dy=−dy∧dz and dx∧dz=−dz∧dx.The three survivors in the standard order dy∧dz, dz∧dx, dx∧dy are the cyclic order that matches the cross product.
Add the three: dω=(Ax+By+Cz)dx∧dy∧dz. ■
For a 2-form, only the derivative in the missing direction survives: d(Ady∧dz)=Axdx∧dy∧dz since Aydy∧dy∧dz=0 and Azdz∧dy∧dz=0. Similarly d(Bdz∧dx)=Bydy∧dz∧dx=Bydx∧dy∧dz (a cyclic shift is an even permutation) and d(Cdx∧dy)=Czdz∧dx∧dy=Czdx∧dy∧dz.
In the plane the same computation gives d(Pdx+Qdy)=(Qx−Py)dx∧dy: the integrand of Green’s theorem. In one variable d(f)=f′dx.
Worked example new (compare MATH 376 worksheets 13, 14 and 28)
(a)df=2xyzdx+x2zdy+x2ydz. The three partial derivatives, as in 5.2.
(b)dω=(∂x(x)−∂y(−y))dx∧dy=2dx∧dy. Qx−Py=1−(−1). Green’s theorem will turn this into twice the area.
(c)P=yz, Q=xz, R=xy: Ry−Qz=x−x=0, Pz−Rx=y−y=0, Qx−Py=z−z=0. So dω=0. It is the field of worksheet 28: curl zero, and indeed ω=d(xyz).
(d)A=x2, B=y2, C=z2, so dω=(2x+2y+2z)dx∧dy∧dz. The divergence of (x2,y2,z2).
Answer. As in steps 1 to 4. Note (c): a form with dω=0 that is also d of something is the situation of the next theorem.
Theorem BRudin 10.20(b); Spivak ch. 4
d squared is zero
If ω has C2 coefficients then d(dω)=0.
In one lineFor every C2 form ω, d(dω)=0. For ω=f a function in R3 this is curl∇f=0, and for a 1-form it is divcurlF=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 f: df=∑ifidxi, where fi=∂f/∂xi.
Definition: d of a function is the sum of its partials times the basic 1-forms.
d(df)=∑i∑jfijdxj∧dxi, where fij=∂j∂if.
Apply d to each term fidxi: d(fidxi)=dfi∧dxi.
=∑j<ifijdxj∧dxi+∑j<ifjidxi∧dxj.
Split the double sum into the terms with j<i and j>i (the terms with j=i vanish because dxi∧dxi=0), and rename in the second sum.
=∑j<ifij(dxj∧dxi−dxj∧dxi)=0.
Equality of mixed partials for C2 functions (5.5): fij=fji. And dxi∧dxj=−dxj∧dxi.These are the same two facts as in the proofs of 6.6, now stated once.
So d(dω)=0 for every C2 form. ■
For a general form fdxI: d(fdxI)=df∧dxI, and applying d again gives d(df)∧dxI=0 by the previous steps, because d(dxI)=0 (dxI 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.
The definition of the exterior derivative on a function.
Start by writing df: you cannot differentiate it before you have it.
The rule that d of a function times a basic form is the differential of the function wedged with that form.
You need df first.
Expand each differential of a first partial; the terms with a repeated dx or dy or dz are zero because a form wedged with itself is zero.
Expand the differentials of the first partials after the rule has been applied.
Reorder with the sign rule dy wedge dx equals minus dx wedge dy, and collect the terms with the same basic form.
Do the expansion before collecting.
The mixed partials of a function with continuous second derivatives are equal, so every bracket is zero.
Collect the brackets first.
Why· One identity for two
Under the dictionary, d∘d=0 on a function is curl∘∇=0, and on a 1-form it is div∘curl=0. In the plane, d(df)=(fyx−fxy)dx∧dy=0 is the test of 6.2: Py=Qx for P=fx, Q=fy. So the ‘necessary condition for a gradient field’ of 6.2 is now one line: if ω=df then dω=d(df)=0.
Closed and exact
Definition B
Closed form
A form ω is closed if dω=0. For ω=Pdx+Qdy this is Qx−Py=0. For ω=Pdx+Qdy+Rdz it is curlF=0. For a 2-form it is divG=0.
In one lineA form ω with dω=0.
Why it existsThis is the test that can be computed. For a 1-form in the plane it is Py=Qx, the necessary condition of 6.2.
Definition B
Exact form
A form ω is exact if ω=dη for some form η defined on the same region. A 1-form ωF is exact exactly when F is a gradient field, F=∇f (6.2). The line integral of an exact 1-form between two points depends only on the points. A first order equation Pdx+Qdy=0 is an ‘exact differential equation’ (6.2) when the form Pdx+Qdy is exact.
In one lineA form ω that is dη for some form η (for a 1-form: ω=df, a gradient).
Why it existsThis is what we want: a potential. In 6.2 it meant a gradient field. By d2=0, every exact form is closed, and the question is when the converse holds.
By d2=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 BRudin ch. 10; Spivak ch. 4; stated for 1-forms in 6.2
Poincare lemma
Let U⊂Rn be star-shaped (there is a point p such that the segment from p to any point of U lies in U) and ω a closed C1 form of degree k≥1 on U. Then ω=dη for a form η of degree k−1 on U.
In one lineOn a star-shaped region (such as a ball or all of Rn) every closed form of degree ≥1 is exact.
Why it is trueIt is the converse of d2=0 on nice regions, and it explains why a gradient field exists wherever curl is zero on R3 and why a vector potential exists wherever divergence is zero on R3.
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 ω=x2+y2−ydx+xdy is closed on the punctured plane, and that it is not exact there.
Solution.
Components.P=x2+y2−y, Q=x2+y2x.
Partial derivatives.Py=(x2+y2)2−(x2+y2)+2y2=(x2+y2)2y2−x2 and Qx=(x2+y2)2(x2+y2)−2x2=(x2+y2)2y2−x2.
Closed.dω=(Qx−Py)dx∧dy=0 on R2∖{0}.
Not exact. On the unit circle r(t)=(cost,sint) we have ∫ω=∫02π(sin2t+cos2t)dt=2π=0. If ω=df, the line integral around every closed path would be 0 (6.2).
What it is. In polar coordinates ω=dθ, the change in angle. A function whose differential is dθ would be θ itself, which increases by 2π 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)). 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 is closed and exact on R2, and find a potential. Then do the same for ydx−xdy.
Solution.
Closed test.Py=2x+3y2 and Qx=2x+3y2. They are equal, so dω=0.
Exact. The plane is star-shaped, so by the Poincare lemma ω is exact.
Find the potential.fx=2xy+y3 gives f=x2y+xy3+h(y). Then fy=x2+3xy2+h′(y) must equal Q=x2+3xy2, so h′=0. The method of 6.2: integrate one equation, differentiate, compare.
Check.df=(2xy+y3)dx+(x2+3xy2)dy.
Second form.ydx−xdy has Py=1 and Qx=−1, so dω=(Qx−Py)dx∧dy=−2dx∧dy=0. It is not closed, hence not exact: ∫ around a counterclockwise circle of radius r is −2πr2.
Answer.f=x2y+xy3 (plus a constant) for the first; the second is not closed.
Your answerLet ω=(3x2y+z)dx+x3dy+xdz. Find ∫Cω along any path from (0,0,0) to (1,2,3).
Test closedness with the curl, then find f by integrating fx and comparing.
ω=df for f=x3y+xz (check: fx=3x2y+z, fy=x3, fz=x), so the integral is f(1,2,3)−f(0,0,0)=2+3=5.
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)), the pullbackφ∗ω is the form in du,dv obtained by replacing x,y,z in the coefficients by x(u,v),y(u,v),z(u,v), and each dx by xudu+xvdv (similarly dy, dz), and expanding with du∧du=0 and dv∧du=−du∧dv. For a curve r(t) the same with dx=x′(t)dt.
In one lineIf φ is a map from parameters (u,v) to (x,y,z), the pullback φ∗ω of a form ω is obtained by substituting x=x(u,v) and so on into the coefficients and dx=xudu+xvdv 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∗ω, and for a surface ∬Sω=∬Tr∗ω, where fdu∧dv is integrated as ∬fdudv.
Proposition BRudin ch. 10 and Spivak ch. 4; computation below
Forms contain the integrals of 6.1 and 6.5
Let F=(P,Q,R) and G=(A,B,C).
In one lineFor a curve, r∗(Pdx+Qdy+Rdz)=F(r)⋅r′dt, so ∫Cω is the line integral. For a surface, r∗(Ady∧dz+Bdz∧dx+Cdx∧dy)=G(r)⋅(ru×rv)du∧dv, so ∬Sω 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.
Pullback of a 1-form along a curve: dx=x′dt, dy=y′dt, dz=z′dt.
r∗(dy∧dz)=(yudu+yvdv)∧(zudu+zvdv).
Pullback of dy∧dz along a surface: dy=yudu+yvdv and dz=zudu+zvdv.
=(yuzv−yvzu)du∧dv.
Expand; du∧du=dv∧dv=0 and dv∧du=−du∧dv.
So r∗(dz∧dx)=(zuxv−zvxu)du∧dv and r∗(dx∧dy)=(xuyv−xvyu)du∧dv.
This is the first component of ru×rv=(yuzv−zuyv,zuxv−xuzv,xuyv−yuxv). The other two pullbacks are the cyclic copies.
r∗ω=G(r)⋅(ru×rv)du∧dv. ■
Add with coefficients A,B,C evaluated at r(u,v).
So the flux integral of 6.5 is ∬G(r)⋅(ru×rv)dudv, and the sign: swapping the roles of u and v changes du∧dv into dv∧du=−du∧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∧dy under polar coordinates x=rcosθ, y=rsinθ, and of dx∧dy∧dz under spherical coordinates. Use it to explain the Jacobian factors.
Solution.
Differentials.dx=cosθdr−rsinθdθ and dy=sinθdr+rcosθdθ.
Wedge.dx∧dy=cosθ⋅rcosθdr∧dθ+(−rsinθ)sinθdθ∧dr=rcos2θdr∧dθ+rsin2θdr∧dθ. The terms dr∧dr and dθ∧dθ vanish, and dθ∧dr=−dr∧dθ.
Simplify.dx∧dy=rdr∧dθ.
Spherical. With x=ρsinφcosθ, y=ρsinφsinθ, z=ρcosφ the same computation gives dx∧dy∧dz=det∂(ρ,φ,θ)∂(x,y,z)dρ∧dφ∧dθ=ρ2sinφdρ∧dφ∧dθ.
Answer.dx∧dy=rdr∧dθ and dx∧dy∧dz=ρ2sinφdρ∧dφ∧dθ. In general the pullback of dx1∧⋯∧dxn under a map is the Jacobian determinant times du1∧⋯∧dun. That is the change of variables formula of 6.4: the ∣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∧dy back under x=u2−v, y=u+v2. The result is Jdu∧dv. What is J at (u,v)=(1,2)?
The coefficient is the determinant of the matrix of partial derivatives, xuyv−xvyu.
dx=2udu−dv and dy=du+2vdv, so dx∧dy=(2u⋅2v−(−1)(1))du∧dv=(4uv+1)du∧dv. At (1,2) this is 9, the Jacobian determinant.
Proposition BRudin ch. 10 and Spivak ch. 4
Pullback commutes with d
Let φ be a C2 map and ω a C1 form. Then φ∗(dω)=d(φ∗ω).
In one lineFor a C2 map φ and a form ω, φ∗(dω)=d(φ∗ω).
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 f is the chain rule: φ∗(df)=φ∗(fxdx+fydy+…)=(fxxu+fyyu+…)du+(…)dv=d(f∘φ). Higher degrees follow because φ∗ respects the wedge product and d(d(⋅))=0. As a check, the form ω=xdy and the map x=uv, y=u+v have φ∗ω=uv(du+dv), so d(φ∗ω)=(vdu+udv)∧(du+dv)=(v−u)du∧dv after expanding (vdu∧dv+udv∧du), while dω=dx∧dy and φ∗(dx∧dy)=(vdu+udv)∧(du+dv)=(v−u)du∧dv. They agree.
The one theorem
Theorem BSpivak ch. 4 and Rudin 10.33; the special cases are 6.2, 6.4 and 6.6
Generalized Stokes theorem
Let M be a compact oriented k-dimensional region (an interval, a curve, a plane region, a surface, a solid) with piecewise smooth boundary ∂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 M. Let ω be a C1 form of degree k−1 on a neighbourhood of M. Then ∫∂Mω=∫Mdω.
In one lineIf M is a compact oriented k-dimensional region with boundary ∂M carrying the induced orientation and ω is a C1(k−1)-form on a neighbourhood of M, then ∫∂Mω=∫Mdω.
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.
M
ω
dω
∫∂Mω=∫Mdω is
an interval [a,b], ∂M={b}−{a}
a function f
f′dx
f(b)−f(a)=∫abf′(x)dx, the fundamental theorem of calculus
a curve from A to B
a function φ
dφ=∇φ⋅dr
φ(B)−φ(A)=∫C∇φ⋅dr (6.2)
a region D in the plane
Pdx+Qdy
(Qx−Py)dx∧dy
∮∂DPdx+Qdy=∬D(Qx−Py)dA, Green (6.4)
a surface S in space
Pdx+Qdy+Rdz
the 2-form of curlF
∮∂SF⋅dr=∬ScurlF⋅ndS, Stokes (6.6)
a solid E
Ady∧dz+Bdz∧dx+Cdx∧dy
(divG)dx∧dy∧dz
∬∂EG⋅ndS=∭EdivGdV, 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.
The case of the unit square I2=[0,1]2 and ω=Pdx+Qdy: ∂I2 is the bottom (y=0, left to right), the right side (x=1, up), the top (right to left) and the left side (down).
On the other side dω=(Qx−Py)dx∧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 k-cube.
Compare the two lines: the four terms are the same.In dimension k 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ω=∫∂Ikφ∗ω.
Definition of the integral over a parametrized piece φ:Ik→M and its boundary: ∫∂φω=∫∂Ikφ∗ω.
=∫Ikd(φ∗ω).
The cube case applied to the form φ∗ω.
=∫Ikφ∗(dω).
Pullback commutes with d.
=∫Mdω. ■
Definition of the integral of dω over M.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.
∫∂Mω=∫∂Ikφ∗ω.
=∫Ikd(φ∗ω).
=∫Ikφ∗(dω).
=∫Mdω.
Worked example new (the proof of Stokes via forms)
Recover Stokes' theorem of 6.6 from Green's theorem and pullback.
Solution.
Setup. Let r:T→S be a smooth parametrization on a rectangle T, with ru×rv pointing the way the normal n does (so that ∂T run counterclockwise maps to ∂S with the induced orientation), and ω=Pdx+Qdy+Rdz the 1-form of F.
Boundary side.∮∂SF⋅dr=∫∂Tr∗ω, since a line integral is the integral of the pullback.
Green on the rectangle.∫∂Tr∗ω=∬Td(r∗ω). Green’s theorem for the form r∗ω=pdu+qdv, which is a planar 1-form on T.
Pullback commutes with d.d(r∗ω)=r∗(dω).
The 2-form dω is the curl.dω=(Ry−Qz)dy∧dz+…, and its pullback is (curlF)(r)⋅(ru×rv)du∧dv by the proposition above.
Answer.∮∂SF⋅dr=∬T(curlF)(r)⋅(ru×rv)dudv=∬S(curlF)⋅ndS: 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 ω=xdy∧dz+ydz∧dx+zdx∧dy on the ball of radius R, and for ω=21(xdy−ydx) on the ellipse a2x2+b2y2≤1.
Solution.
Ball, the sphere side.r∗ω=F(r)⋅(rφ×rθ)dφ∧dθ with F=(x,y,z)=Rρ^ and rφ×rθ=R2sinφρ^ outward, so F⋅N=R3sinφ and ∬=2π⋅2⋅R3=4πR3.
Ball, the volume side.dω=(1+1+1)dx∧dy∧dz=3dx∧dy∧dz, and ∭3dV=3⋅34πR3=4πR3. They agree.
Ellipse, the volume (area) side.dω=21(dx∧dy−dy∧dx)=dx∧dy, and ∬dxdy=πab.
Ellipse, the boundary side. With x=acost, y=bsint: r∗ω=21(acost⋅bcost+bsint⋅asint)dt=21abdt and ∫02π=πab. They agree.
Answer.4πR3 and π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 ∬Szdx∧dy over the unit sphere S with the outward orientation.
Compute dω first. A cyclic shift of three basic 1-forms does not change the sign.
d(zdx∧dy)=dz∧dx∧dy=dx∧dy∧dz (a cyclic shift is an even permutation), so the integral is the volume of the unit ball, 34π. This is the flux of the field (0,0,z), whose divergence is 1.
Quick checkIn the generalized Stokes theorem for a solid E in R3, what degree is ω?
M=E is 3-dimensional, so dω must be a 3-form, and ω has degree 3−1=2: it is the 2-form of the field whose flux is on the boundary. That gives the divergence theorem.
Quick checkFor M=[a,b] and ω=f a function, what is ∫∂Mω?
The boundary of the interval is the two end points, with + at b (outward normal points along the interval) and − at a. The theorem says f(b)−f(a)=∫abf′(x)dx.
Quick checkWhy can you not apply the generalized Stokes theorem to the angle form dθ on the unit disc?
The theorem needs ω to be C1 on the whole region. The angle form is defined away from the origin only, which is the same hole that breaks Green's theorem in 6.4.
Python labStokes numerically on a wavy patch, and the sign flip when the parameters are swapped
import numpy as npdef F(x, y, z): return (y * z, x**2, np.sin(x) + y * z) # P, Q, Rdef 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 handdef 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θ 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: 1 for the punctured plane, 0 for the plane. In complex analysis, dz/z=dr/r+idθ and ∮∣z∣=1dz/z=i∫02πdθ=2πi because the integral of the exact part dr/r=d(lnr) around a closed path is 0 and the angle form contributes 2π. That single computation is the heart of the residue theorem. In the other direction, the divergence-free field G=(x,y,z)/ρ3 of 6.6 is closed and not exact as a 2-form on R3∖{0}: its integral over the unit sphere is 4π, whereas the integral of an exact 2-form over a closed surface is 0.
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 T and volume V is the 1-form δQ=CVdT+VRTdV, with R=8.314 J/(mol K) and CV=23R for a monatomic gas. Is δQ closed? d(δQ)=VRdT∧dV=0 (the coefficient of dV has a T-derivative R/V, and the coefficient of dT has none). So heat is not exact: there is no function 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 T: TδQ=CVTdT+RVdV has d(TδQ)=0 because the coefficient of dV does not depend on T and that of dT does not depend on V. It is exact: TδQ=d(CVlnT+RlnV). That function is the entropy S, and 1/T is an integrating factor in the sense of the exact differential equations of 6.2. From 300 K and 1 litre to 600 K and 2 litres, with CV=23R, the change is ΔS=(23R+R)ln2=2.5⋅8.314⋅0.693≈14.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=23R, R=8.314 J/(mol K)) goes from (T,V)=(300 K,1 L) to (600 K,2 L). The form δQ/T is exact. What is the entropy change in J/K?
Integrate each term: ∫dT/T=ln(T2/T1) and ∫dV/V=ln(V2/V1).
δQ/T=CVdT/T+RdV/V=d(CVlnT+RlnV), so ΔS=CVln2+Rln2=2.5Rln2=2.5⋅8.314⋅0.6931≈14.41 J/K.
A preview· Preview: Maxwell's equations in two lines
In 6.8 the four Maxwell equations are written with ∇⋅ and ∇×. In the language of this chapter the electric and magnetic fields together make up one 2-form F in four-dimensional spacetime, and two of the equations (no magnetic monopoles, and Faraday’s law) say dF=0, while the other two (Gauss’s law and Ampere’s law with Maxwell’s correction) say that d of the Hodge star of F equals the current 3-form. That dF=0 implies that F=dA locally, with A 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∧dx and dx∧dx=0. A form dx∧dy is not the same as the area element dxdy: the first has an orientation and the second does not. Integrating dy∧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=xudu+xvdv. Leaving dx unchanged, or replacing it by x′(t) without the dt, is the form version of the slip of dotting F(x,y) with 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ω=0 is a local test; ω=df is global. The angle form is closed everywhere on the punctured plane and exact nowhere around the hole. Also d2=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 F as Fx,Fy,Fz, which invites confusion with partial derivatives. In the language of forms the coefficients are P,Q,R (or F1,F2,F3), and a subscript x always means a derivative. Write the component names first and the derivatives after.
Problems
check
Compute dω for ω=x2ydx+xy2dy.
Solution
P=x2y, Q=xy2. dω=(Qx−Py)dx∧dy=(y2−x2)dx∧dy.
check
Compute (dx+dy)∧(dx−dy) and (xdy+ydx)∧dx.
Solution
(dx+dy)∧(dx−dy)=dx∧dx−dx∧dy+dy∧dx−dy∧dy=−dx∧dy−dx∧dy=−2dx∧dy. And (xdy+ydx)∧dx=xdy∧dx+ydx∧dx=−xdx∧dy.
checknew (the pullback of 6.1)
Compute the pullback of ω=ydx+xdy under r(t)=(t2,t3) and integrate from t=0 to 1. Compare with the potential.
Solution
dx=2tdt, dy=3t2dt, so r∗ω=t3⋅2tdt+t2⋅3t2dt=5t4dt and the integral is 1. The potential is f=xy with f(1,1)−f(0,0)=1.
check
Which of these 1-forms on R2 is closed: ydx+xdy, ydx−xdy, xdx+ydy?
Solution
d of Pdx+Qdy is (Qx−Py)dx∧dy: first 1−1=0 (closed, =d(xy)); second −1−1=−2 (not closed); third 0−0=0 (closed, =d2x2+y2).
usenew (Final 3 of MATH 376 in the language of forms)
Show that ω=yzcos(xyz)dx+xzcos(xyz)dy+xycos(xyz)dz is exact, and find ∫Cω for C the unit circle in the plane z=0.
Solution
ω=d(sin(xyz)) by the chain rule: ∂xsin(xyz)=yzcos(xyz), and the same for y, z. So ω is exact and the integral around a closed curve is 0 (the function returns to its value). Also dω=d(dsin(xyz))=0 by d2=0: this is the curl computed term by term in 6.6, in one line.
usenew (generalized Stokes, a surface integral by boundary)
Compute ∬Sxdy∧dz where S is the unit sphere with the outward orientation.
Solution
ω=xdy∧dz has dω=dx∧dy∧dz. By the generalized Stokes theorem with M the unit ball, ∬Sω=∭Bdω=Vol(B)=34π. (Equivalently the flux of (x,0,0), with divergence 1.)
usenew (pullback with a Jacobian)
Pull dx∧dy back under x=2u+v, y=u−3v and compute the area of the image of the square [0,1]2.
Solution
dx=2du+dv, dy=du−3dv. dx∧dy=2⋅(−3)du∧dv+1⋅1dv∧du=(−6−1)du∧dv=−7du∧dv. This is det[211−3]=−7. The map reverses orientation, and the area is ∣−7∣⋅1=7.
usenew (thermodynamics, 6.7)
Show that ω=pdV+Vdp is exact and ω′=pdV is not closed on the (V,p)-plane. What does the second say about the work ∫pdV done by a gas?
Solution
Write ω=pdV+Vdp: in coordinates (V,p) with the form PdV+Qdp, P=p, Q=V, dω=(QV−Pp)dV∧dp=(1−1)dV∧dp=0, and ω=d(pV). For ω′=pdV: P=p, Q=0, dω′=(0−1)dV∧dp=−dV∧dp=0. So the work ∫pdV depends on the path in the (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 f and a 1-form ω, d(fω)=df∧ω+fdω.
Solution
Write ω=∑iPidxi. Then fω=∑ifPidxi, and by the definition d(fPidxi)=d(fPi)∧dxi=(fdPi+Pidf)∧dxi by the product rule for functions. Sum over i: d(fω)=f∑idPi∧dxi+df∧∑iPidxi=fdω+df∧ω.
provenew (exact implies closed, with the line integral)
Prove that if ω=Pdx+Qdy is exact with C2 potential f, then ∮Cω=0 for every closed piecewise smooth path, and dω=0.
Solution
ω=df means P=fx and Q=fy. Then dω=(Qx−Py)dx∧dy=(fyx−fxy)dx∧dy=0 by equality of mixed partials. For a closed path r on [a,b], ∮ω=∫ab(fxx′+fyy′)dt=∫abdtdf(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 A of the plane, A∗(dx∧dy)=detAdx∧dy. Deduce that a linear map multiplies areas by ∣detA∣.
Solution
Let A=[acbd], so the new coordinates are x′=au+bv, y′=cu+dv in terms of old ones. dx′=adu+bdv and dy′=cdu+ddv. Then dx′∧dy′=addu∧dv+bcdv∧du=(ad−bc)du∧dv. The area of the image of a region D is ∬D∣detA∣dudv=∣detA∣area(D), since the area counts regardless of orientation.
stretchnew (the Poincare lemma for 1-forms)
Let ω=∑iPidxi be a closed C1 1-form on a ball centred at 0 in Rn, meaning ∂jPi=∂iPj. Prove that f(x)=∫01∑iPi(tx)xidt satisfies df=ω.
Solution
Differentiate under the integral sign: ∂jf=∫01[Pj(tx)+∑itxi(∂jPi)(tx)]dt. By closedness, ∂jPi=∂iPj, so the sum is ∑itxi(∂iPj)(tx)=tdtd[Pj(tx)] (chain rule: dtdPj(tx)=∑ixi(∂iPj)(tx)). So ∂jf=∫01[Pj(tx)+tdtdPj(tx)]dt=∫01dtd[tPj(tx)]dt=Pj(x). Hence df=ω. The integral is the line integral of ω along the ray from 0 to x, the method of 6.2.
stretchnew (the cube proof in three variables)
Verify the generalized Stokes theorem on the unit cube [0,1]3 for ω=fdy∧dz, and say what the other two kinds of term do.
Solution
dω=fxdx∧dy∧dz, so ∫[0,1]3dω=∬[f(1,y,z)−f(0,y,z)]dydz by Fubini and the fundamental theorem of calculus in x. On the boundary, dy∧dz restricted to the faces y= const or z= const is 0 (one of dy, dz is zero there), so only the faces x=1 (outward normal +x, orientation +) and x=0 (outward −x, orientation −) contribute, giving ∬f(1,y,z)dydz−∬f(0,y,z)dydz, the same. The terms gdz∧dx and hdx∧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 Pdx+Qdy+Rdz (vector fields), 2-forms are flux integrands (vector fields again), 3-forms are volume integrands.
The wedge product is bilinear with dxi∧dxj=−dxj∧dxi and dxi∧dxi=0. On 1-forms in space it is the cross product; in the plane it is the determinant.
The exterior derivatived raises the degree by one: gradient on functions, curl on 1-forms, divergence on 2-forms. In the plane d(Pdx+Qdy)=(Qx−Py)dx∧dy.
d2=0, because mixed partials are equal and the wedge is anti-symmetric. This is both curl∇f=0 and divcurlF=0.
A form is closed if dω=0 and exact if ω=dη. Exact implies closed; the converse holds on star-shaped regions (Poincare lemma) and fails around holes: the angle form dθ.
The pullback substitutes a parametrization, with dx=xudu+xvdv. It turns forms into the integrals of 6.1 and 6.5, the cross product ru×rv appears as the coefficient of du∧dv, and the Jacobian determinant is the pullback of dx∧dy.
Generalized Stokes:∫∂Mω=∫Mdω. 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.
Wikipedia, Differential form · the definitions in this chapter in full generality, with the dictionary to vector calculus in three dimensions
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.