Función de Green

En matemáticas, una función de Green es una función matemática usada como núcleo de un operador lineal integral y usada en la resolución de ecuaciones diferenciales inhomogéneas con condiciones de contorno especificadas. La función de Green recibe ese nombre por el matemático británico George Green, que desarrolló el concepto hacia 1830.

El término también aparece en física, particularmente en teoría cuántica de campos, para referirse a varios tipos de funciones de correlación y operadores integrales para ciertas magnitudes calculables a partir del operador de campo.

Motivación intuitiva

El término función de Green se usa para designar a un operador lineal K que tiene forma de integral, siendo el núcleo de este operador integral la función de Green propiamente dicha. Para explicar que es la función de Green consideremos un operador diferencial lineal L que actúa sobre cierto espacio de funciones definidas sobre una variedad diferenciable M, y pongamos que pretendemos resolver la ecuación diferencial:

(1) $L[u(x)] = f(x) \qquad x\in\Omega\subset M$

La idea del método basado en la función de Green es encontrar una función de dos variables G(x, s) continua y diferenciable en el sentido de la teoría de distribuciones que cumpla:

(2) $L[G(x,s)] = \delta(x-s)\,$

Donde $\delta()\;$ es la distribución delta de Dirac. Si se puede hallar una función G que cumpla la ecuación (2) entonces la solución de la ecuación (1) sea cual sea la función f puede escribirse en la forma:

(3) $u(x) = K[f(x)] := \int G(x,s)f(s)\quad ds$

Puede verse informalmente que la solución así calculada es solución de la ecuación (1) ya que:

$L[u(x)] = \int L[G(x,s)] f(s) ds = \int \delta(x-s)f(s) ds = f(x)$

Por tanto, tenemos la siguiente relación entre el operador integral dado por la función de Green y el operador diferencial asociado a la ecuación diferencial:

$L\circ K = K\circ L = Id \qquad \Rightarrow \qquad L^{-1} = K$

Conviene añadir algunas precisiones al planteamiento informal que hemos presentado:

1. Si el núcleo de L no es trivial, entonces la función de Green no es única, aunque en la práctica una combinación de las simetrías del problema, las condiciones de contorno y otros criterios prácticos externos nos proporcionan un única función de Green.
2. La función de Green G usualmente no es una Función matemática ordinaria sino que puede ser una distribución o función generalizada.
3. No cualquier operador diferencial lineal L admite función de Green. En el caso más general K es sólo un inverso por la derecha de L.

Las funciones de Green son muy útiles en teoría de la materia condensada donde permiten resolver la ecuación de difusión y también en mecánica cuántica donde la función de Green del hamiltoniano es un concepto clave, para el desarrollo de la teoría cuántica de campos.

Definición formal

Para definir la función de Green que hace de núcleo integral del operador que resuelve cierta ecuación diferencial inhomogénea es necesario introducir algunos conceptos. Empezando con un operador diferencial de Sturm-Liouville $\scriptstyle L$ de la forma:

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

Y expresando mediante operador D las condiciones de frontera de Dirichlet:

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

Sea $\scriptstyle f(x)$ una función continua en [0,l], con la cual planteamos el siguiente problema:

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

Este es un problema regular, lo cual significa, que para la ecuación homogénea la única solución existente es la solución trivial.

Teorema Solamente existe una solución u (x) que satisface

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

Dicha solución viene además dada por la siguiente expresión:

$u(x) = \int_0^\ell f(s) G(x,s) \, ds$

En la cual G (x, s) es la función de Green que satisface las siguientes condiciones:

1. G (x, s) es continua x y s.
2. Para $x \ne s$, $L G( x, s ) = 0 \,$.
3. Para $s \ne 0, l$, $D G( x, s ) = 0 \,$.
4. Salto en la derivada: $G'(s_{ + 0}, s ) - G'(s_{ - 0}, s ) = 1 / p(s). \,$
5. Simetría: G(x, s) = G(s, x).

Ejemplos

Ejemplo introductorio

Dado el problema

$\begin{cases} \cfrac{d}{dx}\left[ \cfrac{d}{dx} u(x) \right] + u(x) = f(x) \\ u(0) = 0, \quad \quad u\left(\frac{\pi}{2}\right) = 0 \end{cases}$

Donde la última línea representa las condiciones de contorno o frontera. Para encontrar la función de Green del problema anterior se siguen los siguienes pasos:

• Primer paso. La función de Green para el operador lineal es definida como la solución para_

$g'' + g=\delta(x-s).\,$

Si $\scriptstyle x\ \ne\ s\,\!$, entonces, la distribución delta asume una valor nulo y la solución general para el problema es

$g(x,s) = A \cos x + B \sin x.\,$

Para $\scriptstyle x\ <\ s$, la condición de frontera en $\scriptstyle x=0$ significa que:

$g(0,s) = c_1 \cos 0 + c_2 \sin 0 = c_1 \cdot 1 + c_2 \cdot 0 = 0, \quad c_1 = 0.$

La ecuación para $\scriptstyle g(\pi/2,s) = 0$ se omite pues $\scriptstyle x \ne \pi/2$ si $\scriptstyle x < s$ y $\scriptstyle s \ne \pi/2$. Para $\scriptstyle x > s$ la condición de frontera en $\scriptstyle x=\pi/2$ implica que:

$g\left(\frac{\pi}{2},s\right) = c_3 \cdot 0 + c_4 \cdot 1 = 0, \quad c_4 = 0.$

La ecuación $\quad g(0,s) = 0$ es omitida por similares razones. Combinando ambos resultados anteriores, obtenemos, finalmente:

$g(x,s)=\left\{\begin{matrix} c_2 \sin x, \;\; x < s \\ c_3 \cos x, \;\; s < x \end{matrix}\right.$

• Segundo paso. A continuación, vamos a encontrar c₂ y c₃. Debemos asegurar la continuidad de la función de green para el intervalo escogido. Cuando $x=s\,\!$ se tiene que:

$c_2 \sin s = c_3 \cos s.\,$

También debemos asegurar la discontinuidad de la primera derivada por integración de la ecuación diferencial de $\scriptstyle x=s-\epsilon$ a $\scriptstyle x=s+\epsilon$ y tomando el límite cuando $\scriptstyle \epsilon$ tiende a cero. Por lo cual, derivando la igualdad anterior y garantizando la discontinuidad de esta, tenemos:

$c_3 \cdot [ - \sin s ] - c_2 \cdot \cos s = 1\,$

En la cual se iguala a 1 pues $\scriptstyle p(x)$= 1. Resolvemos para las constantes. $\scriptstyle c_2$ y $\scriptstyle c_3$ obteniendo:

$c_2 = - \cos s \quad ; \quad c_3 = - \sin s$

Entonces, la función de Green es:

$g(x,s) = \begin{cases} -\cos(s) \sin(x) & x < s, \\ -\sin(s) \cos(x) & s < x \end{cases}$

• Solución final, recopilando los resultados anteriores tenemos que la solución final al problema planteado es:

$u(x) = \int_0^{\pi/2} f(s)g(x,s)\ ds = -\int_0^x f(s)\sin s\cos x\ ds - \int_x^{\pi/2} f(s)\cos s \sin x\ ds$

Dicha solución existe para cualquier función $\scriptstyle f \in L^1([0,\pi/2])$ integrable en el intervalo $\scriptstyle [0,\pi/2]$.

Oscilador armónico forzado

En el caso de un oscilador armónico tenemos la siguiente ecuación diferencial

$mx''(t)+bx'(t)+kx(t)=F(t) \,$

Siendo $\scriptstyle F(t)$ la fuerza que provoca la oscilación. Supondremos que la fuerza comienza actuar en t=0 de modo que:

$F(t)=\begin{cases} 0&t<0 \\ f(t) &t \geq 0 \end{cases}$

Asumiendo como condiciones de contorno x(0)=x'(0)=0 la solución de la ecuación de movimiento es:

(*) $x(t)=\int_{-\infty}^{\infty}G(t,s)F(s)ds$

con

$G(t,s) = \begin{cases} \frac{1}{m\omega_1}e^{-\beta (t-s)} \sin(\omega_1(t-s)) &t\geq s\\ 0 & t<s\end{cases} \qquad\qquad \omega_1^2=\frac{k}{m}-\beta^2 \qquad \beta=\frac{b}{2m}$

Eliminando la parte nula de la función de Green resulta:

$x(t)=\int_{0}^{t}\frac{1}{m\omega_1}e^{-\beta (t-s)} \sin(\omega_1(t-s))f(s)ds \qquad\qquad\forall t\geq 0$

La integral (*) se conoce como integral de Duhamel.

Ecuaciones diferenciales con coeficientes constantes

Una ecuación diferencial lineal con coeficientes constantes de orden n-ésimo se caracteriza por ser de la forma:

$\sum_{k=0}^{n}a_k\frac{d^kx(t)}{dt^k}=F(t)$

Supondremos que F(t) es de la forma:

$F(t)=\begin{cases} 0&t<0 \\ f(t) &t \geq 0 \end{cases}$

La solución que cumple las condiciones de contorno:

$\frac{d^{n-1}x(t)}{dt^{n-1}}|_{t=0}=\frac{d^{n-2}x(t)}{dt^{n-2}}|_{t=0}=\cdot\cdot\cdot=\frac{dx(t)}{dt}|_{t=0}=x(0)=0$

viene dada por:

$x(t)=\int_{-\infty}^{\infty}G(t,s)F(s)ds$

con

$G(t,s) = \begin{cases} \frac{dx_h(u)}{du}|_{t-s} &t\geq s\\ 0 & t<s\end{cases}$

siendo x_h la solución particular de la ecuación homogénea que verifica:

$\begin{cases} \sum_{k=0}^{n}a_k\cfrac{d^kx_{h}(t)}{dt^k}=0\\ \cfrac{d^{n-1}x_h(t)}{dt^{n-1}}|_{t=0}=\cfrac{d^{n-2}x_h(t)}{dt^{n-2}}|_{t=0}=\dots = \cfrac{dx_h(t)}{dt}|_{t=0}=0\\ x_h(0)=-\cfrac{1}{a_0} \end{cases}$

La solución de la ecuación inhogénea viene dada por tanto:

$x(t)=\int_{0}^{t}\frac{dx_h(u)}{du}\bigg|_{t-s}f(s)ds, \qquad \forall t\geq 0$

Aplicaciones

El uso principal del formalismo de la función de Green en matemáticas y física es la resolución de ecuaciones diferenciales inhomogéneas con condiciones de contorno dadas. En física las funciones de Green además son usadas como propagadores en el cálculo de diagramas de Feynman.

Bibliografía

• Eyges, Leonard, The Classical Electromagnetic Field, Dover Publications, New York, 1972. ISBN 0-486-63947-9. (Chapter 5 contains a very readable account of using Green's functions to solve boundary value problems in electrostatics.)
• A. D. Polyanin and V. F. Zaitsev, Handbook of Exact Solutions for Ordinary Differential Equations (2nd edition), Chapman & Hall/CRC Press, 2003. ISBN 1-58488-297-2
• A. D. Polyanin, Handbook of Linear Partial Differential Equations for Engineers and Scientists, Chapman & Hall/CRC Press, Boca Raton, 2002. ISBN 1-58488-299-9
• Jerry B. Marion Dinámica clásica de las partículas y sistemas

Enlaces externos

• Weisstein, Eric W. «Green's Function» (en inglés). MathWorld. Wolfram Research.
• Green's function for differential operator en PlanetMath
• Green's function en PlanetMath
• Tutorial on Green's function