Differential forms — what integrals are really made of

The grad/curl/div of vector calculus and its four integral theorems collapse into a single line: Mdω=Mω\int_M d\omega = \int_{\partial M} \omega on a manifold.

Opening

Undergraduate vector calculus asks you to memorize three differential operators (grad, curl, div) and four integral theorems (fundamental theorem of calculus, Green’s, Kelvin–Stokes, Gauss divergence) as if they were unrelated facts. This chapter shows that all four are the same one-line statement, Mdω=Mω\int_M d\omega = \int_{\partial M} \omega, written four times with different dimensions. The cast is small: 1-forms, the wedge product, and the exterior derivative dd. By the end of the chapter you should be able to rewrite a vector-calculus identity as a statement about forms, and know which form-language sentence reduces to which classical theorem. Chapter 7’s Lagrangian mechanics will use this language immediately when we take the variation of the action.

Main 1 — 1-forms, the duals of vectors

A 1-form ω\omega at a point is a linear map that eats a single vector and returns a real number. It is therefore the dual of a vector. In coordinates x1,,xnx^1, \ldots, x^n we define the differential dxidx^i to be the machine that extracts the ii-th component of a vector. Any 1-form is then a linear combination of these basis elements,

ω=ωidxi\omega = \omega_i \, dx^i

(Einstein summation), and a 1-form field is one whose coefficients ωi\omega_i depend on position.

The most familiar example is the exact differential of a smooth function f(x,y,z)f(x, y, z):

df=fxdx+fydy+fzdzdf = \frac{\partial f}{\partial x}\, dx + \frac{\partial f}{\partial y}\, dy + \frac{\partial f}{\partial z}\, dz

In freshman physics this was just “the change in ff under a small displacement drd\vec r”. Read again now, dfdf is not a piece of suggestive notation: it is a 1-form, the one whose action on a displacement vector returns the first-order change of ff.

Main 2 — The wedge product and higher-degree forms

To handle oriented areas and volumes we need a way to multiply two 1-forms into a 2-form. The wedge product \wedge is the antisymmetric multiplication defined by

dxdy=dydx,dxdx=0.dx \wedge dy = -\, dy \wedge dx, \qquad dx \wedge dx = 0.

Everything else follows from that one antisymmetry. A generic 2-form in three dimensions is

ω=Pdydz+Qdzdx+Rdxdy,\omega = P\, dy \wedge dz + Q\, dz \wedge dx + R\, dx \wedge dy,

and when fed two vectors u,v\vec u, \vec v it returns the signed area of the parallelogram they span, weighted by P,Q,RP, Q, R. On an nn-dimensional manifold the top-degree form — the nn-form — is exactly a volume element. We define an orientation as the choice of a nowhere-vanishing top form; that is the cleanest way to fix signed volumes without picking coordinates.

Main 3 — The exterior derivative dd unifies grad, curl, div

The exterior derivative dd is a linear operator that turns a kk-form into a (k+1)(k+1)-form. On a 0-form (a function) ff it is the differential we already met: df=ifdxidf = \partial_i f\, dx^i. On a 1-form ω=ωidxi\omega = \omega_i\, dx^i,

dω=(iωj)dxidxj=12(iωjjωi)dxidxj,d\omega = (\partial_i \omega_j)\, dx^i \wedge dx^j = \tfrac{1}{2}\,(\partial_i \omega_j - \partial_j \omega_i)\, dx^i \wedge dx^j,

because antisymmetry of the wedge throws away the symmetric part automatically. The single most important property of dd is

d2=0,d^2 = 0,

which is just the commuting of partial derivatives in disguise. In three dimensions this one operator collapses into all three classical operators:

  • dd of a 0-form ff is the 1-form dfdf — its components are the gradient f\nabla f.
  • dd of a 1-form ω=Axdx+Aydy+Azdz\omega = A_x dx + A_y dy + A_z dz is a 2-form whose components are the curl ×A\nabla \times \vec A.
  • dd of a 2-form is a 3-form whose coefficient is the divergence B\nabla \cdot \vec B.

And the identity d2=0d^2 = 0 is then the two classical identities ×f=0\nabla \times \nabla f = 0 and (×A)=0\nabla \cdot (\nabla \times \vec A) = 0 stated at once.

Main 4 — Stokes’ theorem in one line

Now the payoff. For any compact oriented manifold MM with boundary M\partial M and any (dimM1)(\dim M - 1)-form ω\omega,

Mdω=Mω.\int_M d\omega = \int_{\partial M} \omega.

This single line contains all four of the integral theorems you memorized as separate facts.

  • M=[a,b]M = [a, b] a 1D interval, ω=f\omega = f — the fundamental theorem of calculus, abfdx=f(b)f(a)\int_a^b f'\, dx = f(b) - f(a).
  • MM a planar region, ω=Pdx+Qdy\omega = P\, dx + Q\, dyGreen’s theorem, M(xQyP)dxdy=M(Pdx+Qdy)\iint_M (\partial_x Q - \partial_y P)\, dx\, dy = \oint_{\partial M} (P\, dx + Q\, dy).
  • MM an oriented 2-surface in R3\mathbb R^3, ω=Adr\omega = \vec A \cdot d\vec rKelvin–Stokes, M(×A)dS=MAdr\iint_M (\nabla \times \vec A) \cdot d\vec S = \oint_{\partial M} \vec A \cdot d\vec r.
  • MM a 3D solid, ω\omega a 2-form — Gauss’ divergence theorem, MBdV=MBdS\iiint_M \nabla \cdot \vec B\, dV = \iint_{\partial M} \vec B \cdot d\vec S.

Four classical theorems become four rewordings of one identity. The moment you have an equation that you learn once and apply in four places, vector calculus shrinks dramatically.

In Python

# Verify Green's theorem numerically on the unit disk.
# omega = -y dx + x dy  →  d omega = 2 dx ^ dy
# Both the boundary integral and the area integral should be near 2π.
import numpy as np

# (a) Boundary integral: parametrize the unit circle by t in [0, 2π)
N = 20000
t = np.linspace(0.0, 2.0 * np.pi, N, endpoint=False)
dt = (2.0 * np.pi) / N
x, y = np.cos(t), np.sin(t)
dxdt, dydt = -np.sin(t), np.cos(t)

# omega(γ'(t)) = (-y)·x'(t) + x·y'(t) = sin^2 t + cos^2 t = 1
integrand_boundary = (-y) * dxdt + x * dydt
line_integral = np.sum(integrand_boundary) * dt

# (b) Area integral: d omega = 2 dx ^ dy, so 2 · area(disk) = 2π
M = 2000
xs = np.linspace(-1.0, 1.0, M)
ys = np.linspace(-1.0, 1.0, M)
X, Y = np.meshgrid(xs, ys, indexing="xy")
inside = (X * X + Y * Y) <= 1.0
cell_area = (xs[1] - xs[0]) * (ys[1] - ys[0])
area_integral = 2.0 * np.sum(inside) * cell_area

print(f"∫∂M ω   = {line_integral:.6f}")
print(f"∫M  dω  = {area_integral:.6f}")
print(f"theory 2π = {2.0 * np.pi:.6f}")

When both numbers land near 2π6.28322\pi \approx 6.2832, you have touched the simplest non-trivial case of Stokes’ theorem by hand. The area integral converges more slowly as the grid is refined because the disk’s boundary is approximated by a staircase.

To the next chapter

Chapter 7: Lagrangian mechanics uses the 1-forms and exterior derivative of this chapter directly when taking the variation of the action S=LdtS = \int L\, dt. Even the familiar Euler–Lagrange left-hand side ddtLq˙Lq=0\frac{d}{dt}\frac{\partial L}{\partial \dot q} - \frac{\partial L}{\partial q} = 0 is, in disguise, a component-form statement that a certain 1-form on configuration-velocity space is closed. With the language of forms in place, the next chapter’s jump from Lagrangian to Hamiltonian mechanics will look like a one-line Legendre transform.