Green's function

From Citizendium, the Citizens' Compendium
Jump to: navigation, search
This article is developing and not approved.
Main Article
Related Articles  [?]
Bibliography  [?]
External Links  [?]
Citable Version  [?]
This editable Main Article is under development and not meant to be cited; by editing it you can help to improve it towards a future approved, citable version. These unapproved articles are subject to a disclaimer.

In physics and mathematics, Green's functions are auxiliary functions in the solution of linear partial differential equations. Green's function is named for the self-taught English mathematician George Green (1793 – 1841), who investigated electricity and magnetism in a thoroughly mathematical fashion. In 1828 Green published a privately printed booklet, introducing what is now called the Green function. This was ignored until William Thomson (Lord Kelvin) discovered it, recognized its great value and had it published nine years after Green's death. Bernhard Riemann gave it the name "Green function".[1]


Let Lx be a given linear differential operator in n variables x = (x1, x2, ..., xn), then the Green function of Lx is the function G(x,y) defined by

where δ(x-y) is the Dirac delta function. Once G(x,y) is known, any differential equation involving Lx is formally solved. Suppose we want to solve,

for a known right hand side ρ(x). The formal solution is

The proof is by verification,

where in the last step the defining property of the Dirac delta function is used.

The integral operator that has the Green function as kernel may be seen as the inverse of a linear operator,

It is illuminating to make the analogy with matrix equations. Let and be n×n matrices connected by

then the solution of a matrix-vector equation is

Make the correspondence ix, jy, and compare the sum over j with the integral over y, and the correspondence is evident.


We consider a case of three variables, n = 3 with x = (x, y, z).

The Green function of


As an important example of this Green function we mention that the formal solution of the Poisson equation of electrostatics, reading

where ε0 is the electric constant and ρ is a charge distribution, is given by


The integral form of the electrostatic field may be seen as a consequence of Coulomb's law. The field at the point x due to the charge dQ = ρ(y)dy is equal to

The field is additive in the charges, so integration gives the total field at x.

CC Image
Integration regions: yellow domain is Uε, large sphere is VR, small sphere is Vε

Proof of Green function of ∇2

Without loss of generality we take x as the origin (0, 0, 0) and replace y by r = (x, y, z) in the above formulation. The length of r is indicated by r.

The proof uses Green's theorem:

where Va is a sphere with radius a and Sa is the surface of this sphere. The smooth test function φ and its gradient vanish for large r,

Further we notice that

because r ≠ 0 in that region (see the figure, where the region is indicated in yellow). This result is most easily proved if we recall that in spherical polar coordinates

First apply Green's theorem to the large sphere of radius R

because by assumption φ and its gradient vanish on SR.

We consider the integral on the right hand side and we will show that

from which the result to be proved follows directly. The main trick is to write

and to consider first the integral over Uε (the yellow domain in the figure) for non-zero, but small, ε. After the integral has been evaluated, the limit for zero ε is taken.

Since Uε has two surfaces, Green's theorem cannot be applied directly, and therefore we write (see the figure),

and apply Green's theorem to the two terms. Recalling that we saw already the first term, we get

The integral between square brackets is zero because ∇2(1/r) is zero on Uε. The last integral can be shown to vanish for small ε. Because φ and its gradient are smooth and finite, and r is constant (equal to ε) on the surface, we may write for small ε

where we assumed that

and that the value of the gradient averaged over the surface may be taken out of the integral. The remaining surface integral is equal to 4πε2.

In order to evaluate the final integral we use

so that

Here we assumed that the test function φ(r) is constant over the surface (isotropic) for small ε. The limit ε → 0 gives the desired result.




  1. M. Kline, Mathematical Thought from Ancient to Modern Times, Oxford University Press, New York (1972) p. 683
  • P. Roman, Advanced Quantum Theory, Addison-Wesley, Reading, Mass. (1965) Appendix 4.
  • I. M. Gel'fand and G. E. Shilov, Generalized Functions, Vol. 1, Academic Press, New York (1964)