# Green's identities

Part of a series of articles about :
Calculus :
$\int _{a}^{b}f'(t)\,dt=f(b)-f(a)$ :
Fundamental theorem
Limits
Continuity
Rolle's theorem
Mean value theorem
Inverse function theorem :
Differential
Definitions
Derivative (generalizations)
Differential
infinitesimal
of a function
total
Concepts
Differentiation notation
Second derivative
Implicit differentiation
Logarithmic differentiation
Related rates
Taylor's theorem
Rules and identities
Sum
Product
Chain
Power
Quotient
L'Hôpital's rule
Inverse
General Leibniz
Faà di Bruno's formula
Reynolds :
Integral
Lists of integrals
Integral transform
Leibniz integral rule
Definitions
Antiderivative
Integral (improper)
Riemann integral
Lebesgue integration
Contour integration
Integral of inverse functions
Integration by
Parts
Discs
Cylindrical shells
Substitution (trigonometric, tangent half-angle, Euler)
Euler's formula
Partial fractions (Heaviside's method)
Changing order
Reduction formulae
Differentiating under the integral sign
Risch algorithm :
Series
Geometric (arithmetico-geometric)
Harmonic
Alternating
Power
Binomial
Taylor
Convergence tests
Summand limit (term test)
Ratio
Root
Integral
Direct comparison
Limit comparison
Alternating series
Cauchy condensation
Dirichlet
Abel :
Vector
Gradient
Divergence
Curl
Laplacian
Directional derivative
Identities
Theorems
Gradient
Green's
Stokes'
Divergence
generalized Stokes
Helmholtz decomposition :
Multivariable
Formalisms
Matrix
Tensor
Exterior
Geometric
Definitions
Partial derivative
Multiple integral
Line integral
Surface integral
Volume integral
Jacobian
Hessian :
Advanced
Calculus on Euclidean space
Generalized functions
Limit of distributions :
Specialized
Fractional
Malliavin
Stochastic
Variations :
Miscellanea
Precalculus
History
Glossary
List of topics
Integration Bee
Mathematical analysis
Nonstandard analysis :
vte :
In mathematics, Green's identities are a set of three identities in vector calculus relating the bulk with the boundary of a region on which differential operators act. They are named after the mathematician George Green, who discovered Green's theorem.

## Green's first identity

This identity is derived from the divergence theorem applied to the vector field F = ψ ∇φ while using an extension of the product rule that ∇ ⋅ (ψ X ) = ∇ψ ⋅X + ψ ∇⋅X: Let φ and ψ be scalar functions defined on some region U ⊂ R$^d$, and suppose that φ is twice continuously differentiable, and ψ is once continuously differentiable. Using the product rule above, but letting X = ∇φ, integrate ∇⋅(ψ∇φ) over U. Then $\int _{U}\left(\psi \,\Delta \varphi +\nabla \psi \cdot \nabla \varphi \right)\,dV=\oint _{\partial U}\psi \left(\nabla \varphi \cdot \mathbf {n} \right)\,dS=\oint _{\partial U}\psi \,\nabla \varphi \cdot d\mathbf {S}$ where ∆ ≡ ∇$^2$ is the Laplace operator, ∂U is the boundary of region U, n is the outward pointing unit normal to the surface element dS and dS = ndS is the oriented surface element.
This theorem is a special case of the divergence theorem, and is essentially the higher dimensional equivalent of integration by parts with ψ and the gradient of φ replacing u and v.
Note that Green's first identity above is a special case of the more general identity derived from the divergence theorem by substituting F = ψΓ, $\int _{U}\left(\psi \,\nabla \cdot \mathbf {\Gamma } +\mathbf {\Gamma } \cdot \nabla \psi \right)\,dV=\oint _{\partial U}\psi \left(\mathbf {\Gamma } \cdot \mathbf {n} \right)\,dS=\oint _{\partial U}\psi \mathbf {\Gamma } \cdot d\mathbf {S} ~.$

## Green's second identity

If φ and ψ are both twice continuously differentiable on U ⊂ R$^3$, and ε is once continuously differentiable, one may choose F = ψε ∇φ − φε ∇ψ to obtain $\int _{U}\left[\psi \,\nabla \cdot \left(\varepsilon \,\nabla \varphi \right)-\varphi \,\nabla \cdot \left(\varepsilon \,\nabla \psi \right)\right]\,dV=\oint _{\partial U}\varepsilon \left(\psi {\partial \varphi \over \partial \mathbf {n} }-\varphi {\partial \psi \over \partial \mathbf {n} }\right)\,dS.$
For the special case of ε = 1 all across U ⊂ R$^3$, then, $\int _{U}\left(\psi \,\nabla ^{2}\varphi -\varphi \,\nabla ^{2}\psi \right)\,dV=\oint _{\partial U}\left(\psi {\partial \varphi \over \partial \mathbf {n} }-\varphi {\partial \psi \over \partial \mathbf {n} }\right)\,dS.$
In the equation above, ∂φ/∂n is the directional derivative of φ in the direction of the outward pointing surface normal n of the surface element dS, ${\partial \varphi \over \partial \mathbf {n} }=\nabla \varphi \cdot \mathbf {n} =\nabla _{\mathbf {n} }\varphi .$
Explicitly incorporating this definition in the Green's second identity with ε = 1 results in $\int _{U}\left(\psi \,\nabla ^{2}\varphi -\varphi \,\nabla ^{2}\psi \right)\,dV=\oint _{\partial U}\left(\psi \nabla \varphi -\varphi \nabla \psi \right)\cdot d\mathbf {S} .$
In particular, this demonstrates that the Laplacian is a self-adjoint operator in the L$^2$ inner product for functions vanishing on the boundary so that the right hand side of the above identity is zero.

## Green's third identity

Green's third identity derives from the second identity by choosing φ = G, where the Green's function G is taken to be a fundamental solution of the Laplace operator, ∆. This means that: $\Delta G(\mathbf {x} ,{\boldsymbol {\eta }})=\delta (\mathbf {x} -{\boldsymbol {\eta }})~.$
For example, in R$^3$, a solution has the form $G(\mathbf {x} ,{\boldsymbol {\eta }})={\frac {-1}{4\pi \|\mathbf {x} -{\boldsymbol {\eta }}\|}}~.$
Green's third identity states that if ψ is a function that is twice continuously differentiable on U, then $\int _{U}\left[G(\mathbf {y} ,{\boldsymbol {\eta }})\,\Delta \psi (\mathbf {y} )\right]\,dV_{\mathbf {y} }-\psi ({\boldsymbol {\eta }})=\oint _{\partial U}\left[G(\mathbf {y} ,{\boldsymbol {\eta }}){\partial \psi \over \partial \mathbf {n} }(\mathbf {y} )-\psi (\mathbf {y} ){\partial G(\mathbf {y} ,{\boldsymbol {\eta }}) \over \partial \mathbf {n} }\right]\,dS_{\mathbf {y} }.$
A simplification arises if ψ is itself a harmonic function, i.e. a solution to the Laplace equation. Then ∇$^2$ψ = 0 and the identity simplifies to $\psi ({\boldsymbol {\eta }})=\oint _{\partial U}\left[\psi (\mathbf {y} ){\frac {\partial G(\mathbf {y} ,{\boldsymbol {\eta }})}{\partial \mathbf {n} }}-G(\mathbf {y} ,{\boldsymbol {\eta }}){\frac {\partial \psi }{\partial \mathbf {n} }}(\mathbf {y} )\right]\,dS_{\mathbf {y} }.$
The second term in the integral above can be eliminated if G is chosen to be the Green's function that vanishes on the boundary of U (Dirichlet boundary condition), $\psi ({\boldsymbol {\eta }})=\oint _{\partial U}\psi (\mathbf {y} ){\frac {\partial G(\mathbf {y} ,{\boldsymbol {\eta }})}{\partial \mathbf {n} }}\,dS_{\mathbf {y} }~.$
This form is used to construct solutions to Dirichlet boundary condition problems. Solutions for Neumann boundary condition problems may also be simplified, though the Divergence theorem applied to the differential equation defining Green's functions shows that the Green's function cannot integrate to zero on the boundary, and hence cannot vanish on the boundary. See Green's functions for the Laplacian or for a detailed argument, with an alternative.
It can be further verified that the above identity also applies when ψ is a solution to the Helmholtz equation or wave equation and G is the appropriate Green's function. In such a context, this identity is the mathematical expression of the Huygens principle, and leads to Kirchhoff's diffraction formula and other approximations.

## On manifolds

Green's identities hold on a Riemannian manifold. In this setting, the first two are ${\begin{aligned}\int _{M}u\,\Delta v\,dV+\int _{M}\langle \nabla u,\nabla v\rangle \,dV&=\int _{\partial M}uNv\,d{\widetilde {V}}\\\int _{M}\left(u\,\Delta v-v\,\Delta u\right)\,dV&=\int _{\partial M}(uNv-vNu)\,d{\widetilde {V}}\end{aligned}}$ where u and v are smooth real-valued functions on M, dV is the volume form compatible with the metric, $d{\widetilde {V}}$ is the induced volume form on the boundary of M, N is the outward oriented unit vector field normal to the boundary, and Δu = div(grad u) is the Laplacian.

## Green's vector identity

Green's second identity establishes a relationship between second and (the divergence of) first order derivatives of two scalar functions. In differential form $p_{m}\,\Delta q_{m}-q_{m}\,\Delta p_{m}=\nabla \cdot \left(p_{m}\nabla q_{m}-q_{m}\,\nabla p_{m}\right),$ where p$_m$ and q$_m$ are two arbitrary twice continuously differentiable scalar fields. This identity is of great importance in physics because continuity equations can thus be established for scalar fields such as mass or energy.
In vector diffraction theory, two versions of Green's second identity are introduced.
One variant invokes the divergence of a cross product and states a relationship in terms of the curl-curl of the field $\mathbf {P} \cdot \left(\nabla \times \nabla \times \mathbf {Q} \right)-\mathbf {Q} \cdot \left(\nabla \times \nabla \times \mathbf {P} \right)=\nabla \cdot \left(\mathbf {Q} \times \left(\nabla \times \mathbf {P} \right)-\mathbf {P} \times \left(\nabla \times \mathbf {Q} \right)\right).$
This equation can be written in terms of the Laplacians,
$\mathbf {P} \cdot \Delta \mathbf {Q} -\mathbf {Q} \cdot \Delta \mathbf {P} +\mathbf {Q} \cdot \left[\nabla \left(\nabla \cdot \mathbf {P} \right)\right]-\mathbf {P} \cdot \left[\nabla \left(\nabla \cdot \mathbf {Q} \right)\right]=\nabla \cdot \left(\mathbf {P} \times \left(\nabla \times \mathbf {Q} \right)-\mathbf {Q} \times \left(\nabla \times \mathbf {P} \right)\right).$
However, the terms $\mathbf {Q} \cdot \left[\nabla \left(\nabla \cdot \mathbf {P} \right)\right]-\mathbf {P} \cdot \left[\nabla \left(\nabla \cdot \mathbf {Q} \right)\right],$ could not be readily written in terms of a divergence.
The other approach introduces bi-vectors, this formulation requires a dyadic Green function. The derivation presented here avoids these problems.
Consider that the scalar fields in Green's second identity are the Cartesian components of vector fields, i.e., $\mathbf {P} =\sum _{m}p_{m}{\hat {\mathbf {e} }}_{m},\qquad \mathbf {Q} =\sum _{m}q_{m}{\hat {\mathbf {e} }}_{m}.$
Summing up the equation for each component, we obtain $\sum _{m}\left[p_{m}\Delta q_{m}-q_{m}\Delta p_{m}\right]=\sum _{m}\nabla \cdot \left(p_{m}\nabla q_{m}-q_{m}\nabla p_{m}\right).$
The LHS according to the definition of the dot product may be written in vector form as $\sum _{m}\left[p_{m}\,\Delta q_{m}-q_{m}\,\Delta p_{m}\right]=\mathbf {P} \cdot \Delta \mathbf {Q} -\mathbf {Q} \cdot \Delta \mathbf {P} .$
The RHS is a bit more awkward to express in terms of vector operators. Due to the distributivity of the divergence operator over addition, the sum of the divergence is equal to the divergence of the sum, i.e., $\sum _{m}\nabla \cdot \left(p_{m}\nabla q_{m}-q_{m}\nabla p_{m}\right)=\nabla \cdot \left(\sum _{m}p_{m}\nabla q_{m}-\sum _{m}q_{m}\nabla p_{m}\right).$
Recall the vector identity for the gradient of a dot product, $\nabla \left(\mathbf {P} \cdot \mathbf {Q} \right)=\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} +\mathbf {P} \times \left(\nabla \times \mathbf {Q} \right)+\mathbf {Q} \times \left(\nabla \times \mathbf {P} \right),$ which, written out in vector components is given by
$\nabla \left(\mathbf {P} \cdot \mathbf {Q} \right)=\nabla \sum _{m}p_{m}q_{m}=\sum _{m}p_{m}\nabla q_{m}+\sum _{m}q_{m}\nabla p_{m}.$
This result is similar to what we wish to evince in vector terms 'except' for the minus sign. Since the differential operators in each term act either over one vector (say $p_{m}$’s) or the other ($q_{m}$’s), the contribution to each term must be $\sum _{m}p_{m}\nabla q_{m}=\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \left(\nabla \times \mathbf {Q} \right),$ $\sum _{m}q_{m}\nabla p_{m}=\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} +\mathbf {Q} \times \left(\nabla \times \mathbf {P} \right).$
These results can be rigorously proven to be correct through evaluation of the vector components. Therefore, the RHS can be written in vector form as
$\sum _{m}p_{m}\nabla q_{m}-\sum _{m}q_{m}\nabla p_{m}=\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \left(\nabla \times \mathbf {Q} \right)-\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} -\mathbf {Q} \times \left(\nabla \times \mathbf {P} \right).$
Putting together these two results, a result analogous to Green's theorem for scalar fields is obtained, Theorem for vector fields: $\color {OliveGreen}\mathbf {P} \cdot \Delta \mathbf {Q} -\mathbf {Q} \cdot \Delta \mathbf {P} =\left[\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \times \left(\nabla \times \mathbf {Q} \right)-\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} -\mathbf {Q} \times \left(\nabla \times \mathbf {P} \right)\right].$
The curl of a cross product can be written as $\nabla \times \left(\mathbf {P} \times \mathbf {Q} \right)=\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} -\left(\mathbf {P} \cdot \nabla \right)\mathbf {Q} +\mathbf {P} \left(\nabla \cdot \mathbf {Q} \right)-\mathbf {Q} \left(\nabla \cdot \mathbf {P} \right);$
Green's vector identity can then be rewritten as $\mathbf {P} \cdot \Delta \mathbf {Q} -\mathbf {Q} \cdot \Delta \mathbf {P} =\nabla \cdot \left[\mathbf {P} \left(\nabla \cdot \mathbf {Q} \right)-\mathbf {Q} \left(\nabla \cdot \mathbf {P} \right)-\nabla \times \left(\mathbf {P} \times \mathbf {Q} \right)+\mathbf {P} \times \left(\nabla \times \mathbf {Q} \right)-\mathbf {Q} \times \left(\nabla \times \mathbf {P} \right)\right].$
Since the divergence of a curl is zero, the third term vanishes to yield Green's vector identity: $\color {OliveGreen}\mathbf {P} \cdot \Delta \mathbf {Q} -\mathbf {Q} \cdot \Delta \mathbf {P} =\nabla \cdot \left[\mathbf {P} \left(\nabla \cdot \mathbf {Q} \right)-\mathbf {Q} \left(\nabla \cdot \mathbf {P} \right)+\mathbf {P} \times \left(\nabla \times \mathbf {Q} \right)-\mathbf {Q} \times \left(\nabla \times \mathbf {P} \right)\right].$
With a similar procedure, the Laplacian of the dot product can be expressed in terms of the Laplacians of the factors $\Delta \left(\mathbf {P} \cdot \mathbf {Q} \right)=\mathbf {P} \cdot \Delta \mathbf {Q} -\mathbf {Q} \cdot \Delta \mathbf {P} +2\nabla \cdot \left[\left(\mathbf {Q} \cdot \nabla \right)\mathbf {P} +\mathbf {Q} \times \nabla \times \mathbf {P} \right].$
As a corollary, the awkward terms can now be written in terms of a divergence by comparison with the vector Green equation, $\mathbf {P} \cdot \left[\nabla \left(\nabla \cdot \mathbf {Q} \right)\right]-\mathbf {Q} \cdot \left[\nabla \left(\nabla \cdot \mathbf {P} \right)\right]=\nabla \cdot \left[\mathbf {P} \left(\nabla \cdot \mathbf {Q} \right)-\mathbf {Q} \left(\nabla \cdot \mathbf {P} \right)\right].$
This result can be verified by expanding the divergence of a scalar times a vector on the RHS.