This identity is derived from the divergence theorem applied to the vector field F = ψ∇φ: Let φ and ψ be scalar functions defined on some region U ⊂ Rd, and suppose that φ is twice continuously differentiable, and ψ is once continuously differentiable. Then
where ∆ is the Laplace operator, ∂U is the boundary of region U, n is the outward pointing unit normal of surface element dS and dS is the oriented surface element.
This form is used to construct solutions to Dirichlet boundary condition problems. To find solutions for Neumann boundary condition problems, the Green's function with vanishing normal gradient on the boundary is used instead.
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.
Green's identities hold on a Riemannian manifold, In this setting, the first two are
where u and v are smooth real-valued functions on M, dV is the volume form compatible with the metric, is the induced volume form on the boundary of M, N is oriented unit vector field normal to the boundary, and Δu = div(grad u) is the Laplacian.
Green’s second identity establishes a relationship between second and (the divergence of) first order derivatives of two scalar functions. In differential form
where pm and qm 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
This equation can be written in terms of the Laplacians,
However, the terms
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.
Summing up the equation for each component, we obtain
The LHS according to the definition of the dot product may be written in vector form as
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.
Recall the vector identity for the gradient of a dot product,
which, written out in vector components is given by
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 ’s) or the other (’s), the contribution to each term must be
^Strauss, Walter. Partial Differential Equations: An Introduction. Wiley.
^M. Fernández-Guasti. Complementary fields conservation equation derived from the scalar wave equation. J. Phys. A: Math. Gen., 37:4107–4121, 2004.
^A. E. H. Love. The Integration of the Equations of Propagation of Electric Waves. Philosophical Transactions of the Royal Society of London. Series A, Containing Papers of a Mathematical or Physical Character, 197:pp. 1–45, 1901.
^J. A. Stratton and L. J. Chu. Diffraction Theory of Electromagnetic Waves. Phys. Rev., 56(1):99–107, Jul 1939.
^N. C. Bruce. Double scatter vector-wave Kirchhoff scattering from perfectly conducting surfaces with infinite slopes. Journal of Optics, 12(8):085701, 2010.
^W. Franz, On the Theory of Diffraction. Proceedings of the Physical Society. Section A, 63(9):925, 1950.
^Chen-To Tai. Kirchhoff theory: Scalar, vector, or dyadic? Antennas and Propagation, IEEE Transactions on, 20(1):114–115, jan 1972.
^M. Fernández-Guasti. Green's second identity for vector fields. ISRN Mathematical Physics, 2012:7, 2012. Article ID: 973968.