Potenciales delta de Dirac

P. H. Rivera*
Facultad de Ciencias Físicas
Universidad Nacional Mayor de San Marcos, Lima, Perú

Ciudad Universitaria, 22 de diciembre del 2021

*
priverar@unmsm.edu.pe

1
1. Potencial delta de Dirac
Un potencial de interés para analizar el efecto de las impurezas en el trans-
porte electrónico en metales y semiconductores es el potencial definido por una
función delta de Dirac.
La función delta de Dirac tiene como dimensiones el recíproco de la longitud.
Para los potenciales atractivos generado por las impurezas sobre los electrones
consideramos los potenciales definido por
~2λ
V (x) = − δ(x) (1)
2ma
donde a es una constante que posee unidades de longitud y representa una dis-
tancia arbitraria y de interés. λ es una constante adimensional que representa
la profundidad del potencial.

2
 6 6
x=a
-
x -
x

x=0

La ecuación a resolver está dado por

~2 d2
− φ(x) + V (x)φ(x) = Eφ(x) (2)
2m dx2
reacomodando tenemos
d2 λ 2mE 2
φ(x) + δ(x)φ(x) = − φ(x) = κ φ(x) (3)
dx2 a ~2

3
aquí consideramos que la energía E < 0, luego
2m|E| 2mE
κ2 = = − (4)
~2 ~
luego reescribimos la ecuación de la forma
d2φ(x) 2 λ
− κ φ(x) = − δ(x)φ(x) . (5)
dx2 a
Observamos que esta ecuación es válida sólo para x = 0, mientras que para
x 6= 0 se cumple
d2φ(x) 2
− κ φ(x) = 0 (6)
dx2
y se observa que la autofunción es cero para x → ±∞ resultando que
e−κx , x > 0



φ(x) = 
 eκx
. (7)
, x<0

4
Figura 1: Las autofunciones de energía para un potencial δ ubicado en x = 0.

5
 Estamos considerando que la amplitud de la autofunción es la unidad. Asu-
mimos que la autofunción es continua en x = 0 pero al analizar la derivada de
la autofunción en general observamos de que alrededor de x = a tenemos que
dφ  dφ  2m
   
Z a+
 −  = a− dx 2 [V (x) − E]φ(x) (8)
dx a+ dx a− ~
y en nuestro caso para x = 0 tenemos
2
dφ dφ 2m λ λ
   
Z  ~
  −  = − dx 2 [− δ(x) − E]φ(x) = − φ(0) . (9)
dx x=0+ dx x=0− ~ 2ma a
De la Ec.(7) tenemos que
λ λ
−κ − κ = − =⇒ κ= . (10)
a 2a
Para calcular el autovalor de energía tenemos que
2
λ 2m|E| ~2  λ 2
   
2
κ =  = =⇒ |E| = (11)
2a ~2 2m 2a
6
2. Dos potenciales delta de Dirac
Consideramos dos pozos de potenciales delta de Dirac de la forma
~2λ
V (x) = − [δ(x − a) + δ(x + a)] (12)
2ma
que están separados uno del otro por una distancia igual a 2a y es un potencial
simétrico bajo el intercambio de x por −x.
6
−a a -
x

Este es el caso límite de dos pozos que están separados una distancia igual
a 2a desde sus centros y el estado fundamental muestra dos ondas tipo s que

7
se acoplan, formando un estado ligante y simétrico. Los estados de los pozos
deben ser caracterizados por su paridad.
Las soluciones de paridad par son descritas por







e−κx , x>a
φ(x) = 

A cosh κx , −a < x < a . (13)
eκx , x < −a





Las condiciones de matching entre las autofunciones y también entre sus

respectivas derivadas. Aquí la masa efectiva de los portadores de carga es el
mismo aún en el potencial delta de Dirac.

8
Para x = −a
e−κa = A cosh κa (14)
λ
−κe−κa − κA senh κa = − e−κa
a
λ −κa
e − e−κa = A senh κa (15)
κa
λ
− 1 = tanh κa (16)
κa
λ eκa − e−κa 1 − e−2κa
− 1 = κa = (17)
κa e + e−κa 1 + e−2κa
haciendo
λ
 
 
−2κa
 −1 1+e
 = 1 − e−2κa
κa
λ
(1 + e−2κa) − (1 + e−2κa) = (1 − e−2κa)
κa
2κa 2κa
(1 + e−2κa) = =⇒ e−2κa = −1
λ λ
9
 De este útlimo término podemos asegurar de que al menos existe una solu-
ción para un estado ligado con esta simetría.
Para λ grande, 2κa debe ser un poco más grande que λ. Si escogemos que
2κa = λ + , entonces
λ+  
e−(λ+) = −1=1+ −1=
λ λ λ

e−λe =
λ

e−λ ≈ =⇒  ≈ λe−λ
λ
de modo que
2κa = λ + λe−λ
La intersección de las funciones
λ
− 1 = tanh κa
κa
10
presenta al menos una solución para esta paridad.
Las soluciones impares están dados por







e−κx , x>a
φ(x) = 

A senh κx , −a < x < a . (18)
−eκx , x < −a





Debido a la antisimetría, aplicamos las condiciones de matching para x =

a.

A senh κa = e−κa (19)

λ
−κe−κa − κA cosh κa = − e−κa (20)
a
λ
 
 − 1 e−κa = A cosh κa
κa
λ eκa + e−κa 1 + e−2κa
− 1 = coth κa = κa =
κa e − e−κa 1 − e−2κa

11
 De la misma forma para encontrar al menos una solución impar para λ  1
tenemos que

2κa = λ − λe−λ (21)

la solución debe ser la intersección entre las funciones

−1
λ
 

tanh κa =  − 1 (22)
κa
se observa que la solución impar es menos ligada que la solución par. Esto veri-
fica que los estados acoplados de dos pozos finitos con un autoestado cada uno
se acoplan formando un estado ligante (simetría par) y otro estado antiligante
(simetría impar).
Una combinación lineal de los estados ligantes (par) y antiligantes (impar)
me permite obtener estados que dependen del tiempo y que presentan a los
electrones en el pozo de la izquierda durante un semiperiodo y en el pozo de
la derecha durante el semiperiodo siguiente. Este comportamiento oscilatorio

12
Figura 2: Autofunciones de energía para potenciales tipo delta de Dirac geme-
los, ubicados a x = −a y x = a.

13
de los electrones entre los dos pozos permite entender que los electrones están
tunelando a través de la barrera intermediaria.
Suponiendo que una combinación lineal tal de las soluciones pares e impares
está dada de la forma
φpar(x) + αφimpar(x)
ψ(x) = √ (23)
1+α 2

se escoge α de modo tal que los electrones se localizen en el pozo de la derecha.

Después de un tiempo t el sistema evoluciona
e−iEpart/~φpar(x) + e−iEimpart/~αφimpar(x)
ψ(x, t) = √
1 + α2
−i(Eimpar −Epar )t/~
−iEpar t/~ φpar (x) + e √
αφimpar(x)
=e (24)
1 + α2
El sistema va a evolucionar de modo que necesariamente
e−i(Eimpar−Epar)t/~ = −1 (25)

14
donde el electrón se localiza en el lado izquierdo de manera similar a como el
electrón estuvo localizado en el pozo de la derecha para t = 0.
La frecuencia de oscilación entre el pozo de la derecha y el pozo de la
izquierda, donde los electrones tunelan a través de la barrera está dado por
Eimpar − Epar
ω = ωip = (26)
~
Para los pozos delta bastante profundos, λ es muy grande, de la Ec.(11)
tenemos que
~2
E=− 2
(κa)2
2ma
y
2κa = λ ± λe−λ
~2 2
E=− (2κa)
8ma2
~2 −λ 2 ~2 −λ 2
Eimpar = − (λ − λe ) E par = − (λ + λe )
8ma2 8ma2
15
 ~
ω=− 2
[(λ − λe−λ)2 − (λ + λe−λ)2]
8ma
~
= 2
4λ2e−λ
8ma
~ −λ+2 ln λ
= e (27)
2ma2
La definición de frecuencia angular
2π
ω=
T
nos permite calcular el tiempo de tunelamiento entre las barreras
4πma2
T = −λ+2 ln λ (28)
~e
Ver los videos:

simple https://www.youtube.com/watch?v=1dW_izzvfOk.

16
gemelos https://www.youtube.com/watch?v=E38myvlfyUQ.

wikipedia https://en.wikipedia.org/wiki/Delta_potential

17