Proofs Chapter 5: Second-Order Differential Identities

Laplacian, curl-of-grad, div-of-curl, and product rules

5. Second-Order Differential Identities

This chapter proves vector calculus identities (36)-(44): the definition of the Laplacian, the fundamental second-order identities (curl of gradient = 0, divergence of curl = 0, curl of curl), and the Laplacian product rules. The central tools are Schwarz's theorem on the symmetry of mixed partial derivatives and the antisymmetry of the Levi-Civita symbol $\varepsilon_{ijk}$ together with the $\varepsilon$-$\delta$ identity.

Assumption: throughout this chapter scalar fields $f$ and vector fields $\mathbf{A}, \mathbf{B}$ are assumed to be of class $C^2$ (continuous second partial derivatives). Schwarz's theorem $\dfrac{\partial^2 f}{\partial x_j \partial x_k} = \dfrac{\partial^2 f}{\partial x_k \partial x_j}$ applies. Einstein's summation convention is used.

5.1 Laplacian

(36) Definition of the Laplacian

Identity: $\nabla \cdot (\nabla f) = \nabla^2 f = \Delta f$
Conditions: $f \in C^2$
Proof

By the definition of gradient, the $i$-th component of $\nabla f$ is

\begin{equation}[\nabla f]_i = \dfrac{\partial f}{\partial x_i} \label{eq:5-36-1}\end{equation}

By the definition of divergence, for a vector field $\mathbf{V}$,

\begin{equation}\nabla \cdot \mathbf{V} = \displaystyle\sum_{i} \dfrac{\partial V_i}{\partial x_i} \label{eq:5-36-2}\end{equation}

Substituting the vector field $\nabla f$ from $\eqref{eq:5-36-1}$ into $\eqref{eq:5-36-2}$,

\begin{equation}\nabla \cdot (\nabla f) = \displaystyle\sum_{i} \dfrac{\partial}{\partial x_i}\!\left(\dfrac{\partial f}{\partial x_i}\right) = \displaystyle\sum_{i} \dfrac{\partial^2 f}{\partial x_i^2} \label{eq:5-36-3}\end{equation}

The operator $\displaystyle\sum_i \dfrac{\partial^2}{\partial x_i^2}$ on the right is called the Laplacian and is denoted by $\nabla^2$ or $\Delta$.

\begin{equation}\nabla \cdot (\nabla f) = \nabla^2 f = \Delta f \label{eq:5-36-4}\end{equation}

Remark: the symbol $\nabla^2$ (read "nabla squared") encodes the formal inner product $\nabla \cdot \nabla$. The notation $\Delta$ emphasizes the Laplacian as a differential operator. Physically the Laplacian appears as the spatial part of the heat equation, the Poisson equation for the electric potential, and the wave equation.

(37) Cartesian component form of the Laplacian

Identity: $\nabla^2 f = \dfrac{\partial^2 f}{\partial x^2} + \dfrac{\partial^2 f}{\partial y^2} + \dfrac{\partial^2 f}{\partial z^2}$
Conditions: 3D Cartesian coordinates, $f \in C^2$
Proof

Take the dimension-independent expression from (36),

\begin{equation}\nabla^2 f = \displaystyle\sum_{i} \dfrac{\partial^2 f}{\partial x_i^2} \label{eq:5-37-1}\end{equation}

and substitute $(x_1, x_2, x_3) = (x, y, z)$.

\begin{equation}\nabla^2 f = \dfrac{\partial^2 f}{\partial x^2} + \dfrac{\partial^2 f}{\partial y^2} + \dfrac{\partial^2 f}{\partial z^2} \label{eq:5-37-2}\end{equation}

Remark: in curvilinear coordinates (cylindrical, spherical) extra metric-dependent terms appear. The Cartesian form is simple precisely because the metric components are constant and every Christoffel symbol vanishes.

5.2 Fundamental Identities

(38) Curl of a gradient is zero

Identity: $\nabla \times (\nabla f) = \mathbf{0}$
Conditions: $f \in C^2$
Proof

Write the $i$-th component of the curl with the Levi-Civita symbol $\varepsilon_{ijk}$:

\begin{equation}[\nabla \times \mathbf{V}]_i = \varepsilon_{ijk}\,\dfrac{\partial V_k}{\partial x_j} \label{eq:5-38-1}\end{equation}

Set $\mathbf{V} = \nabla f$, i.e. $V_k = \dfrac{\partial f}{\partial x_k}$:

\begin{equation}[\nabla \times (\nabla f)]_i = \varepsilon_{ijk}\,\dfrac{\partial}{\partial x_j}\!\left(\dfrac{\partial f}{\partial x_k}\right) = \varepsilon_{ijk}\,\dfrac{\partial^2 f}{\partial x_j \partial x_k} \label{eq:5-38-2}\end{equation}

Since $f \in C^2$, Schwarz's theorem (equality of mixed partials) gives

\begin{equation}\dfrac{\partial^2 f}{\partial x_j \partial x_k} = \dfrac{\partial^2 f}{\partial x_k \partial x_j} \label{eq:5-38-3}\end{equation}

so the second-derivative tensor is symmetric in $(j, k)$. The Levi-Civita symbol is antisymmetric in $(j, k)$:

\begin{equation}\varepsilon_{ijk} = -\varepsilon_{ikj} \label{eq:5-38-4}\end{equation}

The contraction of a symmetric tensor $S_{jk} = S_{kj}$ with an antisymmetric tensor $A_{jk} = -A_{kj}$ vanishes. Pairing the $(j, k)$ and $(k, j)$ terms in the sum,

\begin{equation}\varepsilon_{ijk}\,\dfrac{\partial^2 f}{\partial x_j \partial x_k} = \dfrac{1}{2}\!\left(\varepsilon_{ijk} + \varepsilon_{ikj}\right)\dfrac{\partial^2 f}{\partial x_j \partial x_k} = 0 \label{eq:5-38-5}\end{equation}

Every component vanishes, so

\begin{equation}\nabla \times (\nabla f) = \mathbf{0} \label{eq:5-38-6}\end{equation}

Remark: physically this says "a potential field has zero vorticity". An electrostatic field $\mathbf{E} = -\nabla\varphi$ satisfies $\nabla \times \mathbf{E} = \mathbf{0}$. Conversely, on a simply connected domain, $\nabla \times \mathbf{F} = \mathbf{0}$ implies the existence of a potential $\mathbf{F} = \nabla\varphi$ (Poincaré's lemma).
The gradient field ∇f is orthogonal to the level curves, and the circulation around any closed curve vanishes, i.e. the curl (vorticity) is zero f minimum level curves f = const f grows outward ∇f (gradient) closed curve C The gradient field ∇f is orthogonal to the level curves, and its circulation around any closed curve C vanishes — hence ∇ ✕ (∇f) = 0 everywhere.
Figure 1. Geometric meaning of the vanishing curl of a gradient. In this figure the centre is a minimum of $f$ and the outer level curves carry larger values. Since $\nabla f$ points in the direction in which $f$ increases, the arrows point outward and are orthogonal to the level curves. The circulation $\oint_C \nabla f\cdot d\mathbf{r}$ along any closed curve $C$ is then always $0$, because $f$ returns to its original value after one loop. For a gradient field this global property (vanishing circulation) and the local property $\nabla\times(\nabla f)=\mathbf{0}$ are two sides of the same fact, and the contraction of a symmetric with an antisymmetric tensor in the proof is its algebraic counterpart. If the centre were a maximum instead, all arrows would simply point inward and the conclusion would be unchanged.

(39) Divergence of a curl is zero

Identity: $\nabla \cdot (\nabla \times \mathbf{A}) = 0$
Conditions: $\mathbf{A} \in C^2$
Proof

The $i$-th component of $\nabla \times \mathbf{A}$ is

\begin{equation}[\nabla \times \mathbf{A}]_i = \varepsilon_{ijk}\,\dfrac{\partial A_k}{\partial x_j} \label{eq:5-39-1}\end{equation}

Apply the definition of divergence $\nabla \cdot \mathbf{V} = \dfrac{\partial V_i}{\partial x_i}$ (summation convention).

\begin{equation}\nabla \cdot (\nabla \times \mathbf{A}) = \dfrac{\partial}{\partial x_i}\!\left(\varepsilon_{ijk}\,\dfrac{\partial A_k}{\partial x_j}\right) = \varepsilon_{ijk}\,\dfrac{\partial^2 A_k}{\partial x_i \partial x_j} \label{eq:5-39-2}\end{equation}

Since $\mathbf{A} \in C^2$, Schwarz's theorem gives

\begin{equation}\dfrac{\partial^2 A_k}{\partial x_i \partial x_j} = \dfrac{\partial^2 A_k}{\partial x_j \partial x_i} \label{eq:5-39-3}\end{equation}

so the second-derivative tensor is symmetric in $(i, j)$. The Levi-Civita symbol is antisymmetric in $(i, j)$ ($\varepsilon_{ijk} = -\varepsilon_{jik}$). By the same "symmetric × antisymmetric = 0" argument as in (38),

\begin{equation}\varepsilon_{ijk}\,\dfrac{\partial^2 A_k}{\partial x_i \partial x_j} = 0 \label{eq:5-39-4}\end{equation}

Therefore

\begin{equation}\nabla \cdot (\nabla \times \mathbf{A}) = 0 \label{eq:5-39-5}\end{equation}

Remark: in Maxwell's theory $\nabla \cdot \mathbf{B} = 0$ (no magnetic monopoles) is a physical law; the identity above does not derive it. The logic runs the other way: a magnetic field satisfying $\nabla \cdot \mathbf{B} = 0$ can be written, locally in general and globally under suitable topological assumptions, as $\mathbf{B} = \nabla \times \mathbf{A}$ with a vector potential $\mathbf{A}$. Once written in that form, the identity guarantees $\nabla \cdot \mathbf{B} = 0$ automatically.

(40) Curl of a curl

Identity: $\nabla \times (\nabla \times \mathbf{A}) = \nabla(\nabla \cdot \mathbf{A}) - \nabla^2 \mathbf{A}$
Conditions: $\mathbf{A} \in C^2$ (the identity itself is coordinate-free; the componentwise proof below is carried out in Cartesian coordinates)
Proof

Write the $i$-th component of the double curl in components:

\begin{equation}[\nabla \times (\nabla \times \mathbf{A})]_i = \varepsilon_{ijk}\,\dfrac{\partial}{\partial x_j}\,[\nabla \times \mathbf{A}]_k \label{eq:5-40-1}\end{equation}

Expand the inner curl with another Levi-Civita symbol:

\begin{equation}[\nabla \times \mathbf{A}]_k = \varepsilon_{klm}\,\dfrac{\partial A_m}{\partial x_l} \label{eq:5-40-2}\end{equation}

Substitute $\eqref{eq:5-40-2}$ into $\eqref{eq:5-40-1}$:

\begin{equation}[\nabla \times (\nabla \times \mathbf{A})]_i = \varepsilon_{ijk}\,\varepsilon_{klm}\,\dfrac{\partial^2 A_m}{\partial x_j \partial x_l} \label{eq:5-40-3}\end{equation}

Apply the contraction identity ($\varepsilon$-$\delta$ identity):

\begin{equation}\varepsilon_{ijk}\,\varepsilon_{klm} = \varepsilon_{kij}\,\varepsilon_{klm} = \delta_{il}\,\delta_{jm} - \delta_{im}\,\delta_{jl} \label{eq:5-40-4}\end{equation}

Substitute $\eqref{eq:5-40-4}$ into $\eqref{eq:5-40-3}$ and contract with the Kronecker deltas:

\begin{equation}[\nabla \times (\nabla \times \mathbf{A})]_i = \delta_{il}\,\delta_{jm}\,\dfrac{\partial^2 A_m}{\partial x_j \partial x_l} - \delta_{im}\,\delta_{jl}\,\dfrac{\partial^2 A_m}{\partial x_j \partial x_l} \label{eq:5-40-5}\end{equation}

The first term sets $l = i$, $m = j$; the second sets $m = i$, $l = j$.

\begin{equation}[\nabla \times (\nabla \times \mathbf{A})]_i = \dfrac{\partial^2 A_j}{\partial x_j \partial x_i} - \dfrac{\partial^2 A_i}{\partial x_j \partial x_j} \label{eq:5-40-6}\end{equation}

By Schwarz's theorem the order of partial derivatives in the first term may be exchanged, giving

\begin{equation}\dfrac{\partial^2 A_j}{\partial x_j \partial x_i} = \dfrac{\partial}{\partial x_i}\!\left(\dfrac{\partial A_j}{\partial x_j}\right) = \dfrac{\partial}{\partial x_i}(\nabla \cdot \mathbf{A}) = [\nabla(\nabla \cdot \mathbf{A})]_i \label{eq:5-40-7}\end{equation}

The second term is the component of the (Cartesian) Laplacian:

\begin{equation}\dfrac{\partial^2 A_i}{\partial x_j \partial x_j} = \displaystyle\sum_{j} \dfrac{\partial^2 A_i}{\partial x_j^2} = \nabla^2 A_i = [\nabla^2 \mathbf{A}]_i \label{eq:5-40-8}\end{equation}

Substituting $\eqref{eq:5-40-7}$ and $\eqref{eq:5-40-8}$ into $\eqref{eq:5-40-6}$:

\begin{equation}[\nabla \times (\nabla \times \mathbf{A})]_i = [\nabla(\nabla \cdot \mathbf{A})]_i - [\nabla^2 \mathbf{A}]_i \label{eq:5-40-9}\end{equation}

True for every $i$, so

\begin{equation}\nabla \times (\nabla \times \mathbf{A}) = \nabla(\nabla \cdot \mathbf{A}) - \nabla^2 \mathbf{A} \label{eq:5-40-10}\end{equation}

Remark: this is the central identity used to derive the wave equation in electromagnetism. Maxwell's equations give $\nabla \times (\nabla \times \mathbf{E}) = -\partial_t (\nabla \times \mathbf{B})$, and this identity rewrites the left-hand side as $\nabla(\nabla \cdot \mathbf{E}) - \nabla^2 \mathbf{E}$. In free space ($\nabla \cdot \mathbf{E} = 0$) this yields the wave equation $\nabla^2 \mathbf{E} = \mu_0 \varepsilon_0 \,\partial_t^2 \mathbf{E}$.

5.3 Laplacian Product Rules

(41) Laplacian of a scalar product

Identity: $\nabla^2(fg) = f\,\nabla^2 g + 2(\nabla f \cdot \nabla g) + g\,\nabla^2 f$
Conditions: $f, g \in C^2$
Proof

By (36) the Laplacian is the divergence of the gradient:

\begin{equation}\nabla^2(fg) = \nabla \cdot \bigl(\nabla(fg)\bigr) \label{eq:5-41-1}\end{equation}

The product rule for the gradient (identity (15), Leibniz rule) gives

\begin{equation}\nabla(fg) = f\,\nabla g + g\,\nabla f \label{eq:5-41-2}\end{equation}

Substitute $\eqref{eq:5-41-2}$ into $\eqref{eq:5-41-1}$:

\begin{equation}\nabla^2(fg) = \nabla \cdot (f\,\nabla g) + \nabla \cdot (g\,\nabla f) \label{eq:5-41-3}\end{equation}

Apply the divergence product rule (identity (23): $\nabla \cdot (h\mathbf{V}) = (\nabla h) \cdot \mathbf{V} + h\,(\nabla \cdot \mathbf{V})$) to each term.

\begin{equation}\nabla \cdot (f\,\nabla g) = (\nabla f) \cdot (\nabla g) + f\,\nabla \cdot (\nabla g) = \nabla f \cdot \nabla g + f\,\nabla^2 g \label{eq:5-41-4}\end{equation}

\begin{equation}\nabla \cdot (g\,\nabla f) = (\nabla g) \cdot (\nabla f) + g\,\nabla \cdot (\nabla f) = \nabla f \cdot \nabla g + g\,\nabla^2 f \label{eq:5-41-5}\end{equation}

Substitute $\eqref{eq:5-41-4}$ and $\eqref{eq:5-41-5}$ into $\eqref{eq:5-41-3}$:

\begin{equation}\nabla^2(fg) = f\,\nabla^2 g + 2(\nabla f \cdot \nabla g) + g\,\nabla^2 f \label{eq:5-41-6}\end{equation}

Remark: a direct component computation gives the same result. Expanding $\dfrac{\partial^2(fg)}{\partial x_i^2} = \dfrac{\partial}{\partial x_i}\!\left(\dfrac{\partial f}{\partial x_i} g + f \dfrac{\partial g}{\partial x_i}\right) = \dfrac{\partial^2 f}{\partial x_i^2} g + 2 \dfrac{\partial f}{\partial x_i}\dfrac{\partial g}{\partial x_i} + f \dfrac{\partial^2 g}{\partial x_i^2}$ and summing over $i$ recovers $\eqref{eq:5-41-6}$.

(42) Definition of the vector Laplacian

Identity: $\nabla^2 \mathbf{A} = (\nabla^2 A_x,\ \nabla^2 A_y,\ \nabla^2 A_z)$
Conditions: Cartesian coordinates, $\mathbf{A} \in C^2$
Proof

In Cartesian coordinates the basis vectors $\mathbf{e}_x, \mathbf{e}_y, \mathbf{e}_z$ are constant (independent of position), so any differential operator passes through the basis and acts only on the components.

The vector Laplacian is the quantity that appears as the second term on the right of (40); from $\eqref{eq:5-40-8}$ in that proof,

\begin{equation}[\nabla^2 \mathbf{A}]_i = \dfrac{\partial^2 A_i}{\partial x_j \partial x_j} = \displaystyle\sum_{j} \dfrac{\partial^2 A_i}{\partial x_j^2} = \nabla^2 A_i \label{eq:5-42-1}\end{equation}

For $i = 1, 2, 3$ (i.e. $i = x, y, z$),

\begin{equation}\nabla^2 \mathbf{A} = (\nabla^2 A_x,\ \nabla^2 A_y,\ \nabla^2 A_z) \label{eq:5-42-2}\end{equation}

Remark: in curvilinear coordinates (cylindrical, spherical) the basis vectors depend on position, so $\nabla^2 \mathbf{A}$ is no longer the componentwise scalar Laplacian. In a general coordinate system one usually adopts identity (40), $\nabla^2 \mathbf{A} = \nabla(\nabla \cdot \mathbf{A}) - \nabla \times (\nabla \times \mathbf{A})$, as the definition of the vector Laplacian.

(43) Laplacian of a scalar times a vector

Identity: $\nabla^2(f\mathbf{A}) = (\nabla^2 f)\mathbf{A} + 2(\nabla f \cdot \nabla)\mathbf{A} + f\,\nabla^2\mathbf{A}$
Conditions: Cartesian coordinates, $f \in C^2$, $\mathbf{A} \in C^2$
Proof

By (42), in Cartesian coordinates the $i$-th component of $\nabla^2(f\mathbf{A})$ is $\nabla^2(fA_i)$:

\begin{equation}[\nabla^2(f\mathbf{A})]_i = \nabla^2(fA_i) \label{eq:5-43-1}\end{equation}

Apply (41) to the scalar fields $f$ and $A_i$:

\begin{equation}\nabla^2(fA_i) = (\nabla^2 f)\,A_i + 2(\nabla f \cdot \nabla A_i) + f\,\nabla^2 A_i \label{eq:5-43-2}\end{equation}

The term $\nabla f \cdot \nabla A_i$ is the $i$-th component of the directional-derivative operator $(\nabla f \cdot \nabla)$ applied to $\mathbf{A}$:

\begin{equation}\nabla f \cdot \nabla A_i = \displaystyle\sum_{j} \dfrac{\partial f}{\partial x_j}\,\dfrac{\partial A_i}{\partial x_j} = [(\nabla f \cdot \nabla)\mathbf{A}]_i \label{eq:5-43-3}\end{equation}

Substitute $\eqref{eq:5-43-3}$ into $\eqref{eq:5-43-2}$ and collect components:

\begin{equation}[\nabla^2(f\mathbf{A})]_i = (\nabla^2 f)\,A_i + 2[(\nabla f \cdot \nabla)\mathbf{A}]_i + f\,[\nabla^2 \mathbf{A}]_i \label{eq:5-43-4}\end{equation}

True for every $i$, hence

\begin{equation}\nabla^2(f\mathbf{A}) = (\nabla^2 f)\mathbf{A} + 2(\nabla f \cdot \nabla)\mathbf{A} + f\,\nabla^2\mathbf{A} \label{eq:5-43-5}\end{equation}

Remark: the operator $(\nabla f \cdot \nabla)$ is the unnormalised directional derivative along the vector $\nabla f$ (the usual directional derivative scaled by $\lVert\nabla f\rVert$); acting componentwise on a vector field, it returns a vector field.

(44) Laplacian of a dot product

Identity: $\nabla^2(\mathbf{A} \cdot \mathbf{B}) = \mathbf{A} \cdot \nabla^2\mathbf{B} + \mathbf{B} \cdot \nabla^2\mathbf{A} + 2\,\displaystyle\sum_{i} (\nabla A_i \cdot \nabla B_i)$
Conditions: Cartesian coordinates, $\mathbf{A}, \mathbf{B} \in C^2$
Proof

By the definition of the dot product,

\begin{equation}\mathbf{A} \cdot \mathbf{B} = \displaystyle\sum_{i} A_i\,B_i \label{eq:5-44-1}\end{equation}

By linearity of the Laplacian,

\begin{equation}\nabla^2(\mathbf{A} \cdot \mathbf{B}) = \displaystyle\sum_{i} \nabla^2(A_i\,B_i) \label{eq:5-44-2}\end{equation}

For each $i$, apply (41) to the scalar fields $A_i$ and $B_i$:

\begin{equation}\nabla^2(A_i\,B_i) = A_i\,\nabla^2 B_i + 2\,(\nabla A_i \cdot \nabla B_i) + B_i\,\nabla^2 A_i \label{eq:5-44-3}\end{equation}

Sum over $i$:

\begin{equation}\nabla^2(\mathbf{A} \cdot \mathbf{B}) = \displaystyle\sum_{i} A_i\,\nabla^2 B_i + 2\,\displaystyle\sum_{i} (\nabla A_i \cdot \nabla B_i) + \displaystyle\sum_{i} B_i\,\nabla^2 A_i \label{eq:5-44-4}\end{equation}

By (42), in Cartesian coordinates $[\nabla^2 \mathbf{B}]_i = \nabla^2 B_i$ and $[\nabla^2 \mathbf{A}]_i = \nabla^2 A_i$, so the first and third sums collapse into dot products:

\begin{equation}\displaystyle\sum_{i} A_i\,\nabla^2 B_i = \mathbf{A} \cdot \nabla^2 \mathbf{B},\quad \displaystyle\sum_{i} B_i\,\nabla^2 A_i = \mathbf{B} \cdot \nabla^2 \mathbf{A} \label{eq:5-44-5}\end{equation}

Substituting $\eqref{eq:5-44-5}$ into $\eqref{eq:5-44-4}$,

\begin{equation}\nabla^2(\mathbf{A} \cdot \mathbf{B}) = \mathbf{A} \cdot \nabla^2 \mathbf{B} + \mathbf{B} \cdot \nabla^2 \mathbf{A} + 2\,\displaystyle\sum_{i} (\nabla A_i \cdot \nabla B_i) \label{eq:5-44-6}\end{equation}

Remark: the third term $\displaystyle\sum_i \nabla A_i \cdot \nabla B_i$ is the sum of componentwise inner products of gradients. In terms of the Jacobian matrices $D\mathbf{A} = (\partial_j A_i)_{ij}$ and $D\mathbf{B} = (\partial_j B_i)_{ij}$ it equals $\displaystyle\sum_{i,j} \partial_j A_i\, \partial_j B_i = D\mathbf{A} : D\mathbf{B}$, the Frobenius inner product (the Hilbert-Schmidt inner product in finite dimensions). It appears, for example, when computing the Laplacian of the kinetic energy density $|\mathbf{v}|^2 / 2 = \mathbf{v} \cdot \mathbf{v} / 2$ in fluid dynamics.

5.4 Further topics: harmonic functions and Green's identities

The identities of this chapter lead directly into the theory of elliptic partial differential equations. Two topics that follow immediately from (36) and the divergence theorem are collected here.

Harmonic functions: a function $f$ satisfying $\nabla^2 f = 0$ is called harmonic. Typical examples are the Coulomb potential $\phi = 1/r$ in free space ($r \neq 0$), the real and imaginary parts of a complex analytic function, and steady-state temperature fields. Harmonic functions enjoy a rich analytic theory, starting with the maximum principle: a nonconstant function that is harmonic on a bounded connected domain and continuous up to the boundary attains its maximum and minimum on the boundary. Adding a source term to the Laplace equation $\nabla^2 f = 0$ gives the Poisson equation $\nabla^2 f = -\rho/\varepsilon_0$.

Suppl. Green's identities

Formula: $\displaystyle\iiint_V f\,\nabla^2 g\,dV = -\iiint_V \nabla f \cdot \nabla g\,dV + \oiint_S f\,\dfrac{\partial g}{\partial n}\,dS$ (first identity)
Conditions: $V$ is a bounded domain with piecewise $C^1$ boundary $S = \partial V$. For the first identity, $f \in C^1(\overline{V})$ and $g \in C^2(\overline{V})$; for the second identity, $f, g \in C^2(\overline{V})$ (the version with $f$ and $g$ interchanged involves $\nabla^2 f$)
Proof

Apply the product rule for the divergence (identity (23)) to the vector field $f\,\nabla g$.

\begin{equation}\nabla \cdot (f\,\nabla g) = \nabla f \cdot \nabla g + f\,\nabla \cdot (\nabla g) = \nabla f \cdot \nabla g + f\,\nabla^2 g \label{eq:5-45-1}\end{equation}

The last equality uses (36), $\nabla \cdot (\nabla g) = \nabla^2 g$. Integrate $\eqref{eq:5-45-1}$ over $V$ and apply the divergence theorem to the left-hand side.

\begin{equation}\iiint_V \nabla \cdot (f\,\nabla g)\,dV = \oiint_S (f\,\nabla g) \cdot \mathbf{n}\,dS = \oiint_S f\,\dfrac{\partial g}{\partial n}\,dS \label{eq:5-45-2}\end{equation}

Here $\dfrac{\partial g}{\partial n} = \nabla g \cdot \mathbf{n}$ is the directional derivative along the outward normal. Equating with the integral of the right-hand side of $\eqref{eq:5-45-1}$ gives Green's first identity.

\begin{equation}\iiint_V f\,\nabla^2 g\,dV = -\iiint_V \nabla f \cdot \nabla g\,dV + \oiint_S f\,\dfrac{\partial g}{\partial n}\,dS \label{eq:5-45-3}\end{equation}

From here on assume in addition that $f \in C^2(\overline{V})$. Writing $\eqref{eq:5-45-3}$ again with $f$ and $g$ interchanged and subtracting, the term $\nabla f \cdot \nabla g$ cancels by symmetry and Green's second identity remains.

\begin{equation}\iiint_V \left(f\,\nabla^2 g - g\,\nabla^2 f\right)dV = \oiint_S \left(f\,\dfrac{\partial g}{\partial n} - g\,\dfrac{\partial f}{\partial n}\right)dS \label{eq:5-45-4}\end{equation}

Remark: the second identity moves "a difference of second derivatives over the volume" to "boundary values and normal derivatives", and it is the starting point for uniqueness theorems and Green's function representations in elliptic PDE theory. Setting $f = g$ gives $\displaystyle\iiint_V f\,\nabla^2 f\,dV = -\iiint_V |\nabla f|^2 dV + \oiint_S f\,\dfrac{\partial f}{\partial n}\,dS$; if $V$ is connected, $f = 0$ on the boundary and $\nabla^2 f = 0$ in $V$, then $\nabla f \equiv \mathbf{0}$, so $f$ is constant and the boundary condition forces $f \equiv 0$ (uniqueness for the Dirichlet problem).

Chapter summary

  • (36)(37) Definition of the Laplacian: $\nabla \cdot (\nabla f) = \nabla^2 f$ is the sum of second partial derivatives.
  • (38) curl-of-grad = 0, (39) div-of-curl = 0: both follow from the structural identity "symmetric (Schwarz) × antisymmetric (Levi-Civita) = 0".
  • (40) curl of curl: the $\varepsilon$-$\delta$ identity decomposes it into $\nabla(\nabla \cdot \mathbf{A}) - \nabla^2 \mathbf{A}$. Key step in deriving the wave equation.
  • (41)-(44) Laplacian product rules: all follow from applying (41) componentwise.
  • §5.4 further topics: applying the divergence theorem to $\nabla \cdot (f\nabla g)$ yields Green's identities and leads to the analysis of harmonic functions ($\nabla^2 f = 0$).

References

Frequently Asked Questions

What are the key formulas involving the Laplacian ∇²?

For a scalar field: $\nabla^2 f=\partial^2f/\partial x^2+\partial^2f/\partial y^2+\partial^2f/\partial z^2$. The product rule gives $\nabla^2(fg)=f\nabla^2g+2(\nabla f)\cdot(\nabla g)+g\nabla^2f$. For vector fields, the vector Laplacian is $\nabla^2\mathbf{F}=\nabla(\nabla\cdot\mathbf{F})-\nabla\times(\nabla\times\mathbf{F})$.

What is a harmonic function?

A function satisfying $\nabla^2 f=0$ is called harmonic. Examples include the Coulomb potential $\phi=1/r$ (for $r>0$), real and imaginary parts of complex analytic functions, and steady-state temperature distributions. By the maximum principle, a nonconstant function that is harmonic on a bounded connected domain and continuous up to the boundary attains its maximum and minimum on the boundary.

What are Green's identities?

Derived from the divergence theorem: Green's first identity is $\iiint_V f\nabla^2g\,dV=-\iiint_V\nabla f\cdot\nabla g\,dV+\oiint_S f\frac{\partial g}{\partial n}\,dS$; the second identity is $\iiint_V(f\nabla^2g-g\nabla^2f)\,dV=\oiint_S(f\frac{\partial g}{\partial n}-g\frac{\partial f}{\partial n})\,dS$. These are fundamental in the theory of elliptic PDEs and potential theory.