REVISTA DE LA ESCUELA DE FSICA, UNAH Junio 2014 Vol. II, No.
La Funcin Delta de Dirac
Fernando Andino, Marlon Recarte, Michael Spilsbury
Universidad Nacional Autnoma de Honduras en el Valle de Sula
Abstract
La funcin de Dirac presenta propiedades verdaderamente tiles para modelar problemas de
la fsica matemtica. Se presenta un resumen detallado de las principales caractersticas y
propiedades de dicha funcin, argumentndose la validez de las mismas.
Palabras clave: Funcin de Dirac, Delta de Dirac, Distribuciones
The Dirac function has properties truly useful to model problems in physics and mathematics.
This document includes a detailed overview of the main characteristics and properties of the
function, arguing the validity of them.
Keywords: Dirac function, Delta function, Distributions
I.
cantidad que depende de un parmetro x y
satisface las condiciones:
(
0
six 6= 0
(x) =
six = 0
Z
(x)dx = 1
Definicin
a funcin de Dirac no es propiamente
una funcin, es una distribucin o funcin generalizada (localmente integrable)
introducida por primera vez por el fsico ingles
Paul Dirac quien la denomino como funcin
impropia, desde el punto de vista matemtico,
requiere el uso de la teora de distribuciones
desarrollada por Laurent Schwartz (1940) e
introducidas anteriormente de forma independiente por Sergi Sbolev (1935).
Puede expresarse como limite de una sucesin
de funciones, lo que permite dar una idea intuitiva y operativa de la misma. Esta funcin
constituye una aproximacin muy til para
funciones pulso o funcin de impulso y representa igual tipo de abstraccin matemtica que
el de una carga o masa puntual, permitiendo
definir las derivadas generalizadas de funciones
discontinuas.
En fsica esta funcin es muy til ya que permite expresar magnitudes singulares en un
punto como lmite de magnitudes continuas,
se usa para representar la distribucin de densidad de una masa unidad concentrada en un
punto.
La delta de Dirac es introducida para representar cierto tipo de infinitos y sus argumentos
son variables reales. Para manejar estos infinitos con una notacin rigurosa introducimos la
Figure 1: Representacin de la delta de Dirac
Si queremos tener una imagen rigurosa de ,
consideremos una funcin de variable real x
que sea nula fuera de un pequeo dominio
de amplitud alrededor del origen x = 0 y
que en el interior de este dominio sea igual a
uno. No importa la forma exacta de la funcin
en el interior de este dominio, con tal que
no sufra en l variaciones innecesariamente
bruscas. Tomando el lmite para 0, esta
funcin tender a confundirse con .
55
REVISTA DE LA ESCUELA DE FSICA, UNAH Junio 2014 Vol. II, No. 1
(x) no es una funcin de x segn la definicin
matemtica ordinaria de funcin (que le exige
tener un valor definido para cada punto de su
dominio) sino algo mas general que llamaremos
funcin impropia. Para destacar su diferencia
con las funciones definidas de modo ordinario.
Por tanto, (x) no es un cantidad que pueda
usarse en anlisis matemtico con tanta generalidad como las funciones ordinarias, y su
uso debe restringirse a cierto tipos de expresiones sencillas para las que sea evidente que
no pueda dar lugar a inconsecuencias lgicas.
II.
Figure 3
a sustituir x por a. Este resultado general es
vlido aunque la funcin de x no sea una funcin numrica, sino un vector o un operador
lineal dependiente de x. El intervalo de integracin en las frmulas anteriores no ha de ser
necesariamente de a , sino que puede ser
cualquier dominio que incluya el punto crtico
en el que la funcin no se anula.
Podemos definir la funcin todava de otro
modo, como la derivada u0 (x) de la funcin
u(x) definida as:
(
0 x<0
u(x) =
1 x>0
Propiedades
La propiedad ms importante de (x) puede
expresarse con la siguiente ecuacin,
Z
f (x) (x)dx = f (0)
Para comprobar que esta nueva definicin es
equivalente a la anterior haremos lo siguiente:
Sean a y b dos nmeros positivos cualquiera
Z
Figure 2
f (x) (x)dx =
a
en la que f (x) es cualquier funcin continua
de x. Es fcil ver la validez de esta ecuacin
a partir de la imagen de (x) que acabamos
de dar. El primer miembro de la ecuacin
anterior slo puede depender de los valores
muy prximos al origen, de forma que sin error
notable podemos sustituir f (x) por su valor en
el origen f (0). Esto resulta como consecuencia
de la definicin.
Si hacemos el cambio en el origen en la ecuacin
Z
f (x) (x)dx = f (0)
f (x)u0 (x)dx
Al integrar por partes
Z
b
f (x)u0 (x)dx = f (x)u(x)
f 0 (x)u(x)dx
= f (b)
f 0 (x)dx
= f (b) [f (b) f (0)]
= f (0)
De acuerdo a esto la funcin aparece al
derivar una funcin discontinua.
Algunas ecuaciones elementales se pueden escribir con la funcin siendo esencialmente
reglas de manipulacin para desarrollos algebraicos. Estas ecuaciones dan el mismo resultado en sus dos miembros y se utilizan como
factores en un integrando.
entonces podemos deducir la frmula
Z
f (x) (x a)dx = f (a)
en la que a es un nmero real cualquiera. Luego
el resultado de multiplicar una funcin de x por
(x a) e integrar para todo x es equivalente
56
REVISTA DE LA ESCUELA DE FSICA, UNAH Junio 2014 Vol. II, No. 1
(x) = (x) La funcin delta es una funcin par. Consideremos una funcin f arbitraria, entonces
Z
( x2 a 2 ) =
1
( (x a) + (x a))
2a
(a x) (x b) = (a b)
tado puede demostrarse usando
f (x) (x)dx
Este resul-
f (x) (x)dx
f (x) (x a)dx = f (a), tomando
= f (0)
Z
=
f (x) (x)dx
f (x) = (x b), tendremos :
Z
(x b) (x a)dx = (a b)
por lo que (x) = (x).
Derivada de la funcin delta Para una
funcin continua f y un nmero a 6= 0 tenemos
que
x (x) = 0 Para probar esto sea f (x) una
funcin continua, notemos que
Z
f (x)x (x)dx = 0
entonces como f (x) es arbitraria el factor
x (x) en el integrando debe ser cero.
1
(x2 a2 ) = 2a
( (x a) + (x a)),
para (a > 0) Consideremos una funcin f
arbitraria, entonces
(t + a) (t)
f (t)
dt
a
Z
1
=
f (t) (t + a)dt
a
Z
f (t) (t)dt
1
= [f (a) f (0)]
a
Si a 0 entonces
f (x) (x2 a2 )dx
a0
f (x) (x2 a2 )dx
a0
0
f (x) (x2 a2 )dx
Z 0
= f (a)
(x2 a2 )dx
Z
+ f (a)
(x2 a2 )dx
Por lo que
Z
f (x) 0 (x)dx = f 0 (0)
El resultado anterior tambin puede probarse
usando integracin por partes
Al hacer el cambio z = x2 a2
1
f (x) (x)dx = f (x) (x)
1
f 0 (x) (x)dx
Z
2
(t + a) (t)
dt
a
= lim
f (t)
1
[f (a) f (0)]
a
= f (0)
lim
=
Z
f (x) (x a )dx
Z a2
dz
= f (a)
(z )
2 z + a2
Z
dz
+ f (a)
(z )
2 z + a2
a2
1
1
= f (a) + f (a)
2a
2a
= f (0)
Este resultado tambin puede generalizarse
Z
De esto entonces
57
f (x) (n) (x)dx = (1)n f (n) (0)
REVISTA DE LA ESCUELA DE FSICA, UNAH Junio 2014 Vol. II, No. 1
1
(x xi ) Si f (x)
df
dx (xi )
es una funcin continua con xi races tal que
f (xi ) = 0. Sea g (x) una funcin arbitraria,
entonces
Z
g (x) [f (x)]dx
[f (x)] =
Tomando la serie de Taylor de la funcin f (x)
alrededor de la raz xi tenemos:
f ( x ) = f ( xi ) +
f 0 ( xi )
( x xi )
1!
Figure 4
f 00 (x
i)
( x xi ) 2 +
2!
X
f n ( xi )
=
( x xi ) n
n!
n=0
Donde el termino f (xi ) = 0 y los trminos
(x xi )n de orden n > 1 son despreciables,
quedando nicamente el trmino f 0 (xi )(x
xi ), por lo cual
Z
g (x) [f (x)]dx
XZ
=
g (x) [f 0 (xi )(x xi )]dx
-1.0
g (x)
III.
X
i
-0.5
1
(x xi ) dx
f 0 ( xi )
0.5
1.0
0.5
1.0
Figure 5
Aplicaciones
I. Representaciones
Los siguientes tres ejemplos constituyen representaciones de la delta de Dirac.
3.0
2.5
2 2
n
(x) = lim en x
n
2.0
sin(nx)
n
x
1.5
(x) = lim
1.0
n
1
n 1 + n2 x2
Todas ellas satisfacen
Z
lim
fn ( x ) g ( x ) = g ( 0 )
0.5
(x) = lim
-1.0
-0.5
Figure 6
58
REVISTA DE LA ESCUELA DE FSICA, UNAH Junio 2014 Vol. II, No. 1
Es importante tener presente que en el uso de
las representaciones, las integrales se realizan
primero y luego se toma el lmite n . El
orden no conmuta.
En muchas aplicaciones de la fsica nos encontraremos con las siguientes representaciones
equivalentes de la delta:
sin(nx)
n
x
Z n
1
= lim
eixt dt
n 2 n
Z
1
eixt dt
=
2
(x) = lim
Figure 8
d n
1
lim
n dx 1 + n2 x2
Ver figura 9
1
( x x0 ) =
2
i(xx0 )t
dt
sin(n(x x0 ))
n
( x x0 )
(x x0 ) = lim
(x x0 ) = lim
n
1
2
1 + n ( x x0 ) 2
Figure 9
II. Derivada de la funcin delta
Notemos el comportamiento de las derivadas
de las representaciones de la funcin delta
n n2 x2
d
e
lim
n dx
III. Representaciones utilizando funciones ortogonales
Si ~e(x) es ortogonal a ~e(x0 ) entonces
Ver figura 7
( x x0 ) = e ( x ) e ( x0 )
p
X w ( x ) w ( x0 )
0
=
n (x) n (x )
En2
n
donde
En2 =
d sin(nx)
n dx
x
2n (x)w (x)dx
En la siguiente tabla se presentan algunos polinomios ortogonales y sus caractersticas
Figure 7
(x x ) =
lim
X
k =0
1
k+
Pk ( x ) Pk ( x 0 )
2
Si consideramos x0 = 0 entonces
Ver figura 8
59
REVISTA DE LA ESCUELA DE FSICA, UNAH Junio 2014 Vol. II, No. 1
Table 1: Usando los polinomios de Legendre
Z
n
Hn
Ln
Pn
Nombre
Hermite
Laguerre
Legendre
I
(, )
[0, )
[1, 1]
w (x)
2
ex
ex
1
En2
n
2 n!
1
2
2n+1
R3
Si consideremos un sistema coordenado curvilneo con coordenadas (1 , 2 , 3 ) y factor de
escala
"
hi =
X
1
(x) =
k+
Pk ( x ) Pk ( 0 )
2
f (~r) 3 (~r ~a)d = f (~a)
x
i
2
2
y
i
z
i
2 #1/2
Podemos expresar la funcin delta
k =0
(~r r~0 ) =
(1 10 ) (2 20 ) (3 30 )
h1
h2
h3
V. Coordenadas esfricas
1 = r, 2 = , 3 =
x = r sin( ) cos(),
y = r sin( ) sin(),
Figure 10
z = r cos( )
Usando los esfricos armnicos
h1 = 1, h2 = r, h3 = r sin( ),
(cos(1 ) cos(2 )) (1 2 )
X
l
X
dv = r2 sin( )drdd
Yl,m (1 , 1 )Yl,m
( 2 , 2 )
Por lo que
l=0 m=l
IV. Delta de Dirac tridimensional
(~r r~0 )
1
(r r0 ) ( 0 ) ( 0 )
= 2
r sin
Es simple generalizar la funcin delta en tres
dimensiones
Si se tiene simetra azimutal
3 (~r) = (z ) (y ) (z )
donde ~r = xx + y y + z z es el vector de posicin. La funcin delta tridimensional es cero
en todo punto excepto en (0, 0, 0). Su integral
de volumen es:
Z
R3
3 (~r)d
Z Z
=
(~r r~0 ) =
1
2r2 sin
Si no hay dependencia de ,
(~r r~0 ) =
(r r0 ) ( 0 )
1
(r r0 )
4r
VI. Coordenadas cilndricas
(x) (y ) (z )
1 = , 2 = , 3 = z
=1
Adems generalizando
h1 = 1, h2 = , h3 = 1, dv = dddz
60
REVISTA DE LA ESCUELA DE FSICA, UNAH Junio 2014 Vol. II, No. 1
La correspondiente funcin de Dirac es
(~r r~0 ) =
1
( 0 ) ( 0 ) (z z0 )
(~r r~0 ) =
r
r2
= 4 3 (~r)
Adems
1
( 0 ) (z z0 )
2
1
r
= 2
r
r
Entonces
1
(~r r~0 ) =
( 0 )
2
VII. Divergencia de
r
r2
Consideremos la funcin vectorial ~v =
calcular la divergencia de ~v obtenemos
~v =
1
= 4 3 (~r)
r
2
r
.
r2
References
Al
[1] Arfken and Weber. Mathematical Methods
for Physicists. Elseiver Academic Press,
2005.
[2] Paul Dirac. Principios de mecnica cuntica. Ediciones Ariel, 1958.
1
1
2 1
r
= 2 (1) = 0
r2 r
r2
r r
Por otro lado si calculamos el flujo de V sobre
una esfera de radio R centrada en el origen
tenemos
Z
I
~v d~a =
r
R2
[3] David J. Griffiths. Introduction to Quantum Mechanics. Pearson Prentice Hall,
2005.
[4] Hoskins. Delta Functions Introduction to
Generalised Functions. Woodhead Publishing, 2009.
(R2 sin( )ddr) = 4
Usando lo anterior y el teorema de la divergencia se tiene una paradoja
Z
I
~v d = ~v d~a
[5] Jackson. Classical electrodynamics. John
Wiley and Sons, 1962.
[6] Morse and Feshbach. Methods of Theoretical Physics Vol 1. McGraw Hill, 1953.
Por lo que podemos definir
61
