Stokes' theorem

In vector calculus, and more generally differential geometry, **Stokes' theorem** (sometimes spelled **Stokes's theorem**, and also called the **generalized Stokes theorem** or the **Stokes–Cartan theorem**) is a statement about the integration of differential forms on manifolds, which both simplifies and generalizes several theorems from vector calculus. Stokes' theorem says that the integral of a differential form ω over the boundary of some orientable manifold Ω is equal to the integral of its exterior derivative dω over the whole of Ω, i.e.,

∫

∂
Ω

# ω

∫
Ω

d
ω

.

{\displaystyle \int\_{\partial \Omega}\omega=\int\_\Omega d\omega\,.}

Stokes' theorem was formulated in its modern form by Élie Cartan in 1945, following earlier work on the generalization of the theorems of vector calculus by Vito Volterra, Édouard Goursat, and Henri Poincaré.

This modern form of Stokes' theorem is a vast generalization of a classical result that Lord Kelvin communicated to George Stokes in a letter dated July 2, 1850. Stokes set the theorem as a question on the 1854 Smith's Prize exam, which led to the result bearing his name. It was first published by Hermann Hankel in 1861. This classical Kelvin–Stokes theorem relates the surface integral of the curl of a vector field **F** over a surface (that is, the flux of curl **F**) in Euclidean three-space to the line integral of the vector field over its boundary (also known as the loop integral).

**Mathematical Interpretation**

Let *γ*: [*a*, *b*] → **R**2 be a piecewise smooth Jordan plane curve. The Jordan curve theorem implies that γ divides **R**2 into two components, a compact one and another that is non-compact. Let D denote the compact part that is bounded by γ and suppose *ψ*: *D* → **R**3 is smooth, with *S* := *ψ*(*D*). If Γ is the space curve defined by Γ(*t*) = *ψ*(*γ*(*t*)) and **F** is a smooth vector field on **R**3, then:

∮
Γ

F

⋅

d

Γ

=

∬
S

∇
×

F

⋅

d

S

{\displaystyle \oint\_\Gamma \mathbf{F}\, \cdot\, d{\mathbf{\Gamma}} = \iint\_S \nabla\times\mathbf{F}\, \cdot\, d\mathbf{S} }

This classical statement, along with the classical divergence theorem, the fundamental theorem of calculus, and Green's theorem are simply special cases of the general formulation stated above.

## Introduction

The fundamental theorem of calculus states that the integral of a function f over the interval [*a*, *b*] can be calculated by finding an antiderivative F of f:

∫
a
b

f
(
x
)

d
# x
F
(
b
)
−
F
(
a
)

.

{\displaystyle \int\_a^b f(x)\,dx = F(b) - F(a)\,.}

Stokes' theorem is a vast generalization of this theorem in the following sense.

* By the choice of F, *dF**dx* = *f*(*x*). In the parlance of differential forms, this is saying that *f*(*x*) *dx* is the exterior derivative of the 0-form, i.e. function, F: in other words, that *dF* = *f* *dx*. The general Stokes theorem applies to higher differential forms ω instead of just 0-forms such as F.
* A closed interval [*a*, *b*] is a simple example of a one-dimensional manifold with boundary. Its boundary is the set consisting of the two points a and b. Integrating f over the interval may be generalized to integrating forms on a higher-dimensional manifold. Two technical conditions are needed: the manifold has to be orientable, and the form has to be compactly supported in order to give a well-defined integral.
* The two points a and b form the boundary of the closed interval. More generally, Stokes' theorem applies to oriented manifolds M with boundary. The boundary ∂*M* of M is itself a manifold and inherits a natural orientation from that of M. For example, the natural orientation of the interval gives an orientation of the two boundary points. Intuitively, a inherits the opposite orientation as b, as they are at opposite ends of the interval. So, "integrating" F over two boundary points a, b is taking the difference *F*(*b*) − *F*(*a*).

In even simpler terms, one can consider the points as boundaries of curves, that is as 0-dimensional boundaries of 1-dimensional manifolds. So, just as one can find the value of an integral (*f* *dx* = *dF*) over a 1-dimensional manifold ([*a*, *b*]) by considering the anti-derivative (F) at the 0-dimensional boundaries ({*a*, *b*}), one can generalize the fundamental theorem of calculus, with a few additional caveats, to deal with the value of integrals (dω) over n-dimensional manifolds (Ω) by considering the antiderivative (ω) at the (*n* − 1)-dimensional boundaries (∂Ω) of the manifold.

So the fundamental theorem reads:

∫

[
a
,
b
]

f
(
x
)

d
# x

∫

[
a
,
b
]

d
# F

∫

{
a

}
−

∪
{
b

}

# F
F
(
b
)
−
F
(
a
)

.

{\displaystyle \int\_{[a, b]} f(x)\,dx = \int\_{[a, b]} dF = \int\_{\{a\}^- \cup \{b\}^+} F = F(b) - F(a)\,.}

## Formulation for smooth manifolds with boundary

Let Ω be an oriented smooth manifold with boundary of dimension n and let α be a smooth n-differential form that is compactly supported on Ω. First, suppose that α is compactly supported in the domain of a single, oriented coordinate chart {*U*, *φ*}. In this case, we define the integral of α over Ω as

∫
Ω

# α

∫

φ
(
U
)

(

φ

−
1

)

∗

α

,

{\displaystyle \int\_\Omega \alpha = \int\_{\varphi(U)} \left(\varphi^{-1}\right)^\* \alpha\,,}

i.e., via the pullback of α to **R***n*.

More generally, the integral of α over Ω is defined as follows: Let {*ψi*} be a partition of unity associated with a locally finite cover {*Ui*, *φi*} of (consistently oriented) coordinate charts, then define the integral

∫
Ω

α
≡

∑
i

∫

U
i

ψ
i

α

,

{\displaystyle \int\_\Omega \alpha \equiv \sum\_i \int\_{U\_i} \psi\_i \alpha\,,}

where each term in the sum is evaluated by pulling back to **R***n* as described above. This quantity is well-defined; that is, it does not depend on the choice of the coordinate charts, nor the partition of unity.

The generalized Stokes theorem reads:

**Theorem.** *(Stokes–Cartan) If

ω

{\displaystyle \omega}

 is a smooth

(
n
−
1
)

{\displaystyle (n-1)}

-form with compact support on smooth

n

{\displaystyle n}

-dimensional manifold-with-boundary

Ω

{\displaystyle \Omega}

,

∂
Ω

{\displaystyle \partial\Omega}

 denotes the boundary of

Ω

{\displaystyle \Omega}

 given the induced orientation, and

i
:
∂
Ω
↪
Ω

{\displaystyle i:\partial \Omega\hookrightarrow \Omega}

 is the inclusion map, then*

∫
Ω

d
# ω

∫

∂
Ω

i
∗

ω

{\displaystyle \int\_\Omega d\omega = \int\_{\partial \Omega} i^\*\omega}

.

Conventionally,

∫

∂
Ω

i
∗

ω

{\textstyle \int\_{\partial\Omega} i^\*\omega}

 is abbreviated as

∫

∂
Ω

ω

{\textstyle \int\_{\partial\Omega} \omega}

, since the pullback of a differential form by the inclusion map is simply its restriction to its domain:

i
∗

# ω
ω

∂
Ω

{\displaystyle i^\*\omega=\omega|\_{\partial\Omega}}

. Here

d

{\displaystyle d}

 is the exterior derivative, which is defined using the manifold structure only. The right-hand side is sometimes written as

∮

∂
Ω

ω

{\textstyle \oint\_{\partial\Omega} \omega}

 to stress the fact that the

(
n
−
1
)

{\displaystyle (n-1)}

-manifold

∂
Ω

{\displaystyle \partial\Omega}

 has no boundary. (This fact is also an implication of Stokes' theorem, since for a given smooth

n

{\displaystyle n}

-dimensional manifold

Ω

{\displaystyle \Omega}

, application of the theorem twice gives

∫

∂
(
∂
Ω
)

# ω

∫
Ω

d
(
d
ω
# )
0

{\textstyle \int\_{\partial(\partial \Omega)}\omega=\int\_\Omega d(d\omega)=0}

 for any

(
n
−
2
)

{\displaystyle (n-2)}

-form

ω

{\displaystyle \omega}

, which implies that

∂
(
∂
Ω
# )
∅

{\displaystyle \partial(\partial\Omega)=\emptyset}

.) The right-hand side of the equation is often used to formulate *integral* laws; the left-hand side then leads to equivalent *differential* formulations (see below).

The theorem is often used in situations where

Ω

{\displaystyle \Omega}

 is an embedded oriented submanifold of some bigger manifold, often

R

k

{\displaystyle \mathbf{R}^k}

, on which the form

ω

{\displaystyle \omega}

 is defined.

## Topological preliminaries; integration over chains

Let M be a smooth manifold. A (smooth) singular k-simplex in M is defined as a smooth map from the standard simplex in **R***k* to M. The group *C**k*(*M*, **Z**) of singular k-chains on M is defined to be the free abelian group on the set of singular k-simplices in M. These groups, together with the boundary map, ∂, define a chain complex. The corresponding homology (resp. cohomology) group is isomorphic to the usual singular homology group *H**k*(*M*, **Z**) (resp. the singular cohomology group *H**k*(*M*, **Z**)), defined using continuous rather than smooth simplices in M.

On the other hand, the differential forms, with exterior derivative, d, as the connecting map, form a cochain complex, which defines the de Rham cohomology groups *H**k*
*dR*(*M*, **R**).

Differential k-forms can be integrated over a k-simplex in a natural way, by pulling back to **R***k*. Extending by linearity allows one to integrate over chains. This gives a linear map from the space of k-forms to the kth group of singular cochains, *Ck*(*M*, **Z**), the linear functionals on *Ck*(*M*, **Z**). In other words, a k-form ω defines a functional

I
(
ω
)
(
c
# )

∮
c

ω

{\displaystyle I(\omega)(c) = \oint\_c \omega}

.
on the k-chains. Stokes' theorem says that this is a chain map from de Rham cohomology to singular cohomology with real coefficients; the exterior derivative, d, behaves like the *dual* of ∂ on forms. This gives a homomorphism from de Rham cohomology to singular cohomology. On the level of forms, this means:

1. closed forms, i.e., *dω* = 0, have zero integral over *boundaries*, i.e. over manifolds that can be written as ∂∑*c* *Mc*, and
2. exact forms, i.e., *ω* = *dσ*, have zero integral over *cycles*, i.e. if the boundaries sum up to the empty set: ∑*c* *Mc* = ∅.

De Rham's theorem shows that this homomorphism is in fact an isomorphism. So the converse to 1 and 2 above hold true. In other words, if {*ci*} are cycles generating the kth homology group, then for any corresponding real numbers, {*ai*}, there exist a closed form, ω, such that

∮

c
i

# ω

a
i

,

{\displaystyle \oint\_{c\_i} \omega = a\_i\,,}

and this form is unique up to exact forms.

Stokes' theorem on smooth manifolds can be derived from Stokes' theorem for chains in smooth manifolds, and vice versa. Formally stated, the latter reads:

**Theorem.** *(Stokes' theorem for chains) If

c

{\displaystyle c}

 is a smooth

k

{\displaystyle k}

-chain in a smooth manifold

M

{\displaystyle M}

, and

ω

{\displaystyle \omega}

 is a smooth

(
k
−
1
)

{\displaystyle (k-1)}

-form on

M

{\displaystyle M}

, then*

∫

∂
c

# ω

∫
c

d
ω

{\displaystyle \int\_{\partial c}\omega = \int\_c d\omega}

.

## Underlying principle

To simplify these topological arguments, it is worthwhile to examine the underlying principle by considering an example for *d* = 2 dimensions. The essential idea can be understood by the diagram on the left, which shows that, in an oriented tiling of a manifold, the interior paths are traversed in opposite directions; their contributions to the path integral thus cancel each other pairwise. As a consequence, only the contribution from the boundary remains. It thus suffices to prove Stokes' theorem for sufficiently fine tilings (or, equivalently, simplices), which usually is not difficult.

## Generalization to rough sets

The formulation above, in which Ω is a smooth manifold with boundary, does not suffice in many applications. For example, if the domain of integration is defined as the plane region between two x-coordinates and the graphs of two functions, it will often happen that the domain has corners. In such a case, the corner points mean that Ω is not a smooth manifold with boundary, and so the statement of Stokes' theorem given above does not apply. Nevertheless, it is possible to check that the conclusion of Stokes' theorem is still true. This is because Ω and its boundary are well-behaved away from a small set of points (a measure zero set).

A version of Stokes' theorem that allows for roughness was proved by Whitney. Assume that D is a connected bounded open subset of **R***n*. Call D a *standard domain* if it satisfies the following property: There exists a subset P of ∂*D*, open in ∂*D*, whose complement in ∂*D* has Hausdorff (*n* − 1)-measure zero; and such that every point of P has a *generalized normal vector*. This is a vector **v**(*x*) such that, if a coordinate system is chosen so that **v**(*x*) is the first basis vector, then, in an open neighborhood around x, there exists a smooth function *f*(*x*2, …, *x**n*) such that P is the graph { *x*1 = *f*(*x*2, …, *x**n*) } and D is the region {*x*1 : *x*1 < *f*(*x*2, …, *x**n*) }. Whitney remarks that the boundary of a standard domain is the union of a set of zero Hausdorff (*n* − 1)-measure and a finite or countable union of smooth (*n* − 1)-manifolds, each of which has the domain on only one side. He then proves that if D is a standard domain in **R***n*, ω is an (*n* − 1)-form which is defined, continuous, and bounded on *D* ∪ *P*, smooth on D, integrable on P, and such that dω is integrable on D, then Stokes' theorem holds, that is,

∫
P

# ω

∫
D

d
ω

.

{\displaystyle \int\_P \omega = \int\_D d\omega\,.}

The study of measure-theoretic properties of rough sets leads to geometric measure theory. Even more general versions of Stokes' theorem have been proved by Federer and by Harrison.

## Special cases

The general form of the Stokes theorem using differential forms is more powerful and easier to use than the special cases. The traditional versions can be formulated using Cartesian coordinates without the machinery of differential geometry, and thus are more accessible. Further, they are older and their names are more familiar as a result. The traditional forms are often considered more convenient by practicing scientists and engineers but the non-naturalness of the traditional formulation becomes apparent when using other coordinate systems, even familiar ones like spherical or cylindrical coordinates. There is potential for confusion in the way names are applied, and the use of dual formulations.

### Kelvin–Stokes theorem

Main articles: Kelvin–Stokes theorem

This is a (dualized) (1 + 1)-dimensional case, for a 1-form (dualized because it is a statement about vector fields). This special case is often just referred to as *Stokes' theorem* in many introductory university vector calculus courses and is used in physics and engineering. It is also sometimes known as the **curl** theorem.

The classical Kelvin–Stokes theorem relates the surface integral of the curl of a vector field over a surface Σ in Euclidean three-space to the line integral of the vector field over its boundary. It is a special case of the general Stokes theorem (with *n* = 2) once we identify a vector field with a 1-form using the metric on Euclidean 3-space. The curve of the line integral, ∂Σ, must have positive orientation, meaning that ∂Σ points counterclockwise when the surface normal, **n**, points toward the viewer.

One consequence of the Kelvin–Stokes theorem is that the field lines of a vector field with zero curl cannot be closed contours. The formula can be rewritten as:

Suppose **F**=(*P*(*x*,*y*,*z*),*Q*(*x*,*y*,*z*),*R*(*x*,*y*,*z*)) is defined in region with smooth surface Σ and has first-order continuous partial derivatives. Then

∬
Σ

(

(

∂
R

∂
y

−

∂
Q

∂
z

)

d
y

d
z

(

∂
P

∂
z

−

∂
R

∂
x

)

d
z

d
x

(

∂
Q

∂
x

−

∂
P

∂
y

)

d
x

d
y

)

=

∮

∂
Σ

(

P

d
x
Q

d
y
R

d
z

)

,

{\displaystyle \begin{align}
& \iint\_\Sigma \Bigg(\left(\frac{\partial R}{\partial y}-\frac{\partial Q}{\partial z} \right)\,dy\,dz +\left(\frac{\partial P}{\partial z}-\frac{\partial R}{\partial x} \right) \,dz\,dx +\left (\frac{\partial Q}{\partial x}-\frac{\partial P}{\partial y}\right)\,dx\,dy\Bigg) \\[4pt]
= {} & \oint\_{\partial\Sigma} \Big(P\,dx+Q\,dy+R\,dz\Big)\,,\end{align}}

where P, Q and R are the components of **F**, and ∂Σ is boundary of region with smooth surface Σ.

### Green's theorem

Green's theorem is immediately recognizable as the third integrand of both sides in the integral in terms of P, Q, and R cited above.

#### In electromagnetism

Two of the four Maxwell equations involve curls of 3-D vector fields, and their differential and integral forms are related by the Kelvin–Stokes theorem. Caution must be taken to avoid cases with moving boundaries: the partial time derivatives are intended to exclude such cases. If moving boundaries are included, interchange of integration and differentiation introduces terms related to boundary motion not included in the results below (see Differentiation under the integral sign):

∇
×

E

=
−

∂

B

∂
t

{\displaystyle \quad \nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}} {\partial t} \quad}

∮
C

E

⋅
d

l

=

∬
S

∇
×

E

⋅
d

A

=
−

∬
S

∂

B

∂
t

⋅
d

A

{\displaystyle \begin{align}\quad \oint\_C \mathbf{E} \cdot d\mathbf{l} &= \iint\_S \nabla \times \mathbf{E} \cdot d\mathbf{A} \\ &= -\iint\_S \frac{\partial \mathbf{B}}{\partial t} \cdot d\mathbf{A}\end{align}\quad }

 Ampère's law (with Maxwell's extension):

∇
×

H

=

J

∂

D

∂
t

{\displaystyle \quad \nabla \times \mathbf{H} = \mathbf{J} + \frac{\partial \mathbf{D}} {\partial t} \quad}

∮
C

H

⋅
d

l

=

∬
S

∇
×

H

⋅
d

A

=

∬
S

J

⋅
d

A

∬
S

∂

D

∂
t

⋅
d

A

{\displaystyle \begin{align}\quad \oint\_C \mathbf{H} \cdot d\mathbf{l} &= \iint\_S \nabla \times \mathbf{H} \cdot d\mathbf{A}\\ &= \iint\_S \mathbf{J} \cdot d\mathbf{A} + \iint\_S \frac{\partial \mathbf{D}}{\partial t} \cdot d\mathbf{A}\end{align}\quad }

 (with C and S not necessarily stationary)

The above listed subset of Maxwell's equations are valid for electromagnetic fields expressed in SI units. In other systems of units, such as CGS or Gaussian units, the scaling factors for the terms differ. For example, in Gaussian units, Faraday's law of induction and Ampère's law take the forms:

∇
×

E

=
−

1
c

∂

B

∂
t

,

∇
×

H

=

1
c

∂

D

∂
t

4
π

c

J

,

{\displaystyle \begin{align} \nabla \times \mathbf{E} &= -\frac{1}{c} \frac{\partial \mathbf{B}} {\partial t}\,, \\
\nabla \times \mathbf{H} &= \frac{1}{c} \frac{\partial \mathbf{D}} {\partial t} + \frac{4\pi}{c} \mathbf{J}\,, \end{align}}

respectively, where c is the speed of light in vacuum.

### Divergence theorem

Likewise, the Ostrogradsky–Gauss theorem (also known as the divergence theorem or Gauss's theorem)

∫

V
o
l

∇
⋅

F

d

V
o
l

=

∮

∂

V
o
l

F

⋅
d
Σ

{\displaystyle \int\_\mathrm{Vol} \nabla \cdot \mathbf{F} \, d\_\mathrm{Vol} = \oint\_{\partial \mathrm{Vol}} \mathbf{F} \cdot d\boldsymbol{\Sigma}}

is a special case if we identify a vector field with the (*n* − 1)-form obtained by contracting the vector field with the Euclidean volume form. An application of this is the case **F** = *f***c** where **c** is an arbitrary constant vector. Working out the divergence of the product gives

c

⋅

∫

V
o
l

∇
f

d

V
o
l

=

c

⋅

∮

∂

V
o
l

f

d
Σ

.

{\displaystyle \mathbf c \cdot \int\_\mathrm{Vol} \nabla f \, d\_\mathrm{Vol} = \mathbf c \cdot \oint\_{\partial \mathrm{Vol}} f\, d\boldsymbol{\Sigma}\,.}

Since this holds for all **c** we find

∫

V
o
l

∇
f

d

V
o
l

=

∮

∂

V
o
l

f

d
Σ

.

{\displaystyle \int\_\mathrm{Vol} \nabla f \, d\_\mathrm{Vol} = \oint\_{\partial \mathrm{Vol}} f\, d\boldsymbol{\Sigma}\,.}