green's function

In mathematics, a Green's function is a type of function used to solve inhomogeneous differential equations subject to boundary conditions. Technically, a Green's function of a linear operator L acting on distributions over a manifold M, at a point x₀, is any solution of (Lf)(x) = δ(x − x₀), where δ is the Dirac delta function. If the kernel of L is nontrivial, then the Green's function is not unique. However, in practice, some combination of symmetry, boundary conditions and/or other externally imposed criteria would give us a unique Green's function. Also, please note Green's functions are distributions in general, not functions, meaning they can have discontinuities. Not every operator L admits a Green's function. A Green's function can also be thought of as a one-sided inverse of L. Green functions are also a useful tool in condensed matter theory - the Green function of the Hamiltonian is a key concept, with important links to the concept of density of states. The Green's function was named for British mathematician George Green, who first developed the concept in the 1830s.

Motivation

Convolving with a Green's function gives solutions to inhomogeneous differentio-integral equations, most commonly a Sturm-Liouville problem . If g is the Green's function of an operator L, then the solution for f of the equation Lf = h is given by

f(x) = \int{ h(s) g(x,s) ds}.

This can be thought as an expansion of f according to Dirac delta function basis (projecting f over δ(x − s)) and a superposition of the solution on each projection.

Green's Function for solving inhomogeneous Boundary value problem

The primary use of Green's functions in mathematics is to solve inhomogeneous boundary value problems. In particle physics, Green's functions are also usually used as propagators in Feynman diagrams (and the phrase "Green's function" is often used for any correlation function).

Working Frame

Let

L

be a linear differential operator in the form of

L = {d \over dx}\left( p(x) {d \over dx} \right) + q(x)

and let D be the boundary conditions operator

Du = \left\{\begin{matrix} \alpha _1 u'(0) + \beta _1 u(0) = 0 \\ \alpha _2 u'(l) + \beta _2 u(l) = 0 \end{matrix}\right.

Let

f(x)

be a continuous function in

0,1

. We shall also suppose that the problem

\begin{matrix}Lu = f \\ Du = 0 \end{matrix}

is regular, i.e. only the trivial solution exists for the homogenous problem.

Theorem

Then there is one and only one solution u(x) which satisfies

\begin{matrix}Lu = f \\ Du = 0 \end{matrix}

and it is given by

u(x) = \int_{0}^{l}{ f(s) g(x,s) ds}

where g(x,s) is Green's function and satisfies the following demands:

g(x,s) is continuous in x and s.
For $x \ne s$ , $L g ( x, s ) = 0$ .
For $s \ne 0, l$ , $D g ( x, s ) = 0$ .
Derivative "jump": $g ' ( s_{ + 0}, s ) - g ' (s_{ - 0}, s ) = 1 / p(s)$ .
Symmetry: g(x,s) = g(s,x).

Example

Given the problem

\begin{matrix}Lu\end{matrix} = u ' ' + u = f( x )

Du =  u(0) = 0 \quad, \quad u( \pi / 2) = 0

Find Green's function. First step: From demand-2 we see that

g(x,s) = c_1 (s) \cdot \cos(x) + c_2 (s) \cdot \sin(x)

For x < s we see from demand-3 that the

c_1 (s) = 0

, while for x > s we see from demand-3 that the

c_2 (s) = 0

(we leave it to the reader to fill in the in-between steps). Summarize the results:

g(x,s)=\left\{\begin{matrix} a(s) \sin(x), & x < s \\ b(s) \cos(x), & s < x \end{matrix}\right.

Second step: Now we shall determine a(s) and b(s). Using demand-1 we get

a(s) \sin(s) = b(s) \cos(s)

Using demand-4 we get

b(s) \cdot ( - \sin(s) ) - a(s) \cdot cos (s) = 1/1 = 1

Using Cramer's rule or by intelligent guess solve for a(s) and b(s) and obtain that

a(s) = - \cos(s) \quad ; \quad b(s) = - \sin(s)

. Check that this automatically satisfies demand-5. So our Green's function for this problem is:

$g(x,s)=\left\{\begin{matrix} -1 \cdot \cos(s) \cdot \sin(x), & x < s \\ -1 \cdot \sin(s) \cdot \cos(x), & s < x \end{matrix}\right.$

Futher examples

Let the manifold be R and L be d/dx. Then, the Heaviside step function H(x − x₀) is a Green's function of L at x₀.
Let the manifold be the quarter-plane { (x, y) : x, y ≥ 0 } and L be the Laplacian. Also, assume a Dirichlet boundary condition is imposed at x=0 and a Neumann boundary condition is imposed at y=0. Then the Green's function is G(x, y ; x₀, y₀)

=\frac{1}{2\pi}\left \ln\sqrt{(x-x_0)^2+(y-y_0)^2}-\ln\sqrt{(x+x_0)^2+(y-y_0)^2}+\ln\sqrt{(x-x_0)^2+(y+y_0)^2}-\ln\sqrt{(x+x_0)^2+(y+y_0)^2}\right