Delta de Dirac
Facultad de Ciencias Fı́sico Matemáticas
Benemérita Universidad Autónoma de Puebla
López Téllez Carlos Alfredo
1. Teorema
La función
1 2 1
δ(~r − r~0 ) = − (∇0 ) (1)
4π |~r − r~0 |
Vale 0 en todas partes, excepto en r~0 = ~r donde tiene un valor infinito. Además
Z Z
δ(~r − r~0 )dv =
0
δ(~r − r~0 )dv 0 = 1 (2)
allspace V0
Si V’ contiene al punto ~r.
2. Demostración
Colocamos el origen en el punto ~r, entonces ~r = 0, deseamos probar que
1 2 1
− (∇0 ) =0 (3)
4π |r~0 |
para r~0 6= ~r. El laplaciano de una función para coordenadas (q1 , q2 , q3 ) se expresa como
" #
2 1 ∂ h2 h3 ∂f ∂ h3 h1 ∂f ∂ h1 h2 ∂f
∇f= + + (4)
h1 h2 h3 ∂q1 h1 ∂q1 ∂q2 h2 ∂q2 ∂q3 h3 ∂q3
Usando coordenadas esféricas (q1 , q2 , q3 ) = (r, θ, φ); los factores de escala son hr = 1, hθ =
r, hφ = rsenθ entonces
1
1 2 1 1 1 1 1
δ(~r − r~0 ) = − (∇0 ) = − ∇2 = − ∇2 (5)
4π ~0
|r | 4π |~r| 4π r
1
aquı́ la función f = r
no depende de θ ni φ, por lo que el laplaciano es
" # " # " #
1 1 ∂ ∂ 1 1 ∂ −1 1 ∂(−1)
∇2 = 2 r2 senθ = 2 senθ r2 2 = 2 = 0 (6)
r r senθ ∂r ∂r r r senθ ∂r r r ∂r
Entonces δ(~r − r~0 ) es 0 para r~0 6= ~r e infinito si r~0 = ~r.
Por otra parte, el Teorema de Gauss afirma
Z I
∇ · Adv = A · ds (7)
V S
para cualquier superficie S que encierre el volumen V, por lo cual
Z Z I
21
1 1
∇ dv = ∇ · ∇ dv = ∇ · ds (8)
V r V r S r
Si elegimos que la superficie S sea una esfera de radio a centrada en el origen entonces r̂ y ds
son paralelos
−1
I I I
1 ∂ 1
∇ · ds = r̂ · ds = ds (9)
S r S ∂r r S r2
El diferencial de superficie para una esfera de radio a es ds = a2 senθdθdφ y como evaluamos
el flujo en la superficie r=a
π 2π
−1 −1 2
I I Z Z
ds = a senθdθdφ = − senθdθ dφ = −4π (10)
S r2 S a2 0 0
De la ecuación (5) y el desarrollo (8)-(10) se obtiene
Z Z
1 1 1
δ(~r − r~0 )dv = − 0
∇2 dv = − (−4π) = 1 (11)
V0 4π V r 4π
Si consideramos ahora a la Delta de Dirac δ(~r − r~0 ) multiplicando a una función vectorial de
varias variables f, obtenemos que esta nueva función es 0 cuando ~r 6= r~0 y solo contribuye f
cuando f(~r) = f(~
r0 ), este valor sale enfrente de la integral, dando
Z
r0 ) =
f(~ δ(~r − r~0 )f(~r)dv (12)
allspace