Unidad N°3: Potenciales Unidimensionales

Ahora que ya conocemos los postulados, lo que sigue se trata de aplicaciones de la mecanica
cuantica. Comenzamos con la clase mas simple de problemas: una particula en una dimension. Si
bien estos problemas son de alguna manera artificiales, comparten muchas de las caracteristicas
de la mecanica cuantica 3D pero poco de su complejidad.

1 Particula libre
El problema mas simple en esta familia es, por su puesto, el de la particula libre. La ecuacion de
Schrodinger es en este caso

^ j i = p^ j i
2
@ (1)
i~ j i = H
@t 2m

Los estados estacionarios son soluciones de la forma

j i = jE ie¡iEt/~

con lo cual, si reemplazamos en (1) se tiene la ecuacion de Schrodinger independiente del tiempo
para jE i, es decir

^ jE i = p^2
H jE i = E jE i
2m

Este problema puede resolverse sin ir a alguna otra base. Primero notamos que cualquier autoes-
tado de p^ es tambien autoestado de p^2, con lo cual si reemplazamos jE i por jpi se tiene,

p^2 p2
jpi = jpi = E jpi
2m 2m

con lo que los valores permitidos de p son

p = (2mE)1/2

En otras palabras, para cada autovalor E hay dos autoestados de p^ ortogonales, dados como

jE ; +i = jp = +(2mE)1/2i
jE ; ¡i = jp = ¡(2mE)1/2i

es decir que al autovalor E le corresponde un subespacio de degeneracion 2, que puede generarse

con los kets de la ecuacion anterior. Fisicamente esto significa que una particula de energia E puede
estar moviendose a la izquierda o a la derecha con momento jpj = (2mE)1/2. Hasta aqui, uno bien
podria decir que esto es exactamente lo que ocurre en mecanica clasica, de modo que podemos
preguntarnos: ¾Que es lo nuevo aqui? Lo que es nuevo es el hecho de que el estado

jE i = jp = +(2mE)1/2i + jp = ¡(2mE)1/2i

es tambien un autoestado de energia E y representa una particula de energia E que puede encon-
trarse moviendose ya sea a la izquierda o a la derecha con momento (2mE)1/2.

1
Para construir la autobase ortonormal completa de H ^ , debemos escoger de cada autoespacio dege-
nerado cualquier par de vectores ortonormales. La eleccion obvia esta dada por los kets jE ; +i
y jE ; ¡i. En terminos de las ideas discutidas anteriormente, estamos usando el autovalor de una
variable compatible p^ como un rotulo extra dentro del espacio degenerado con respecto a la energia.
Dado que p^ es un operador no-degenerado, el rotulo p por si solo es adecuado. En otras palabras,
no hay necesidad de llamar al estado como jp; E = p2 /2mi, ya que el valor de E = E(p) se obtiene
de p. Podemos entonces tirar este rotulo redundante.
El propagador temporal (operador de evolucion) es entonces

Z 1 Z 1
U^ (t) =
2
dp jpie¡iE(p)t/~ hpj = dp jpie¡ip t/2m~
hpj
¡1 ¡1

y puede evaluarse explicitamente en la base de coordenadas, a partir del elemento de matriz

Z 1
U (x; t; x 0; 0) = hxjU^ (t)jx 0i =
2
dp hxjpie¡ip t/2m~
hpjx 0i
¡1

tal que
Z 1
 1/2
1 ip(x¡x 0)/~ ¡ip2t/2m~ m 0 2
U (x; t; x ; 0) =
0
dpe e = eim(x¡x ) /2~t
2~ ¡1 2~it

Recordemos que, a partir del propagador, podemos resolver cualquier problema de valores iniciales,
es decir
Z 1
j (t)i = U^ (t)j (0)i ! (x; t) = hxj (t)i = dx 0 hxjU^ (t)jx 0ihx 0j (0)i
¡1

con lo que en este caso resulta

Z 1
(x; t) = dx 0 U (x; t; x 0; 0) (x 0; 0)
¡1

Obviamente, podriamos haber asignado algun otro valor al tiempo inicial, de modo que lo anterior
seria de la forma
Z 1
(x; t) = dx 0 U (x; t; x 0; t 0) (x 0; t 0) (2)
¡1

donde U (x; t; x 0; t 0) = hxjU^ (t ¡ t 0)jx 0i, ya que para Hamiltonianos independientes del tiempo U^ solo
depende del intervalo temporal t ¡ t 0. Si t 0 es asignado como cero, entonces escribimos U (x; t; x 0;
0) = U (x; t; x 0).
Podemos dar una interpretacion fisica para el propagador si comenzamos con una particula loca-
lizada en x 0 = x00 , esto es, con (x 0; t 0) = (x 0 ¡ x00 ). Entonces

(x; t) = U (x; t; x00 ; t 0)

En otras palabras, el propagador (en la base x) es la amplitud de probabilidad de que una particula
que, comenzando en el punto x00 al tiempo t 0, termine en el punto x al tiempo t. Obviamente, lo
mismo puede interpretarse en cualquier base:
h!jU^ (t ¡ t 0)j! 0i corresponde a la amplitud de probabilidad de que una particula (o sistema),
comenzando en el estado j! 0i al tiempo t 0, termine en el estado j!i al tiempo t.

2
La Ec. (2) nos dice entonces que la amplitud total de que la particula se encuentre en x al tiempo t
es la suma de las contribuciones de todos los puntos x 0 con peso proporcional a la amplitud inicial
(x 0; t 0) de que la particula haya estado en x 0 al tiempo t 0.

Evolucion temporal del paquete Gaussiano

Hay una tradicion que dice que la derivacion del propagador de la particula libre debe ser seguida
por su aplicacion al paquete Gaussiano. Consideremos entonces un paquete de ondas como funcion
de onda inicial:
02 2
0 e¡x /2

(x 0; 0) = eip0x /~
(
2)1/4

p
Este paquete tiene posicion
p media hx^ i = 0, una incerteza
x =
/ 2 , y momento medio hp^i = p0
con incerteza
p = ~/ 2
. Aplicando la Ec. (2) se obtiene para la amplitud en x al tiempo t lo
siguiente

  ¡1/2   
~t (x ¡ p0t/m)2 p0 p0t
(x; t) =  1/2

+i exp ¡ +i x¡
m
2
2(1 + i~t/m
2) ~ 2m

y probabilidad
      
~t 2 ¡1/2 (x ¡ p0t/m)2
P (x; t) = 
2 + exp ¡ 2
m

+ (~t/m
)2

Las caracteristicas principales de este resultado son las siguientes:

1. La posicion media de la particula es

p0t hp^it
hx^ i = =
m m

En otras palabras, la relacion clasica x = (p/m)t ahora se mantiene para cantidades promedio. Esta
es una de las consecuencias del teorema de Ehrenfest que establece que las ecuaciones clasicas que
obedecen las variables dinamicas tienen contrapartes en mecanica cuantica como relaciones entre
valores de expectacion. Veremos esto en mayor detalle mas adelante.
2. El ancho del paquete crece con el tiempo:

   

~t 2 1/2

x(t) = p 1 +
2 m
2

Esto tiene que ver con que cualquier incerteza en la velocidad inicial se refleja, al pasar el tiempo,
como un crecimiento en la incerteza en la posicion. Imaginando la evolucion del paquete en forma
clasica, a partir de puntos y velocidades, esta consecuencia es facil de ver. De todos modos, el hecho
de que la incerteza en la velocidad inicial sea inevitable (ya que queremos especificar la posicion
con una precision
) es un hecho puramente cuantico.
Si la particula fuese macroscopica, digamos de masa 1g, y quisieramos fijar su posicion inicial dentro
del ancho de un proton, que es aproximadamente 10¡13 cm, la incerteza en velocidad seria

~

v(0) ' p ' 10¡14 cm/s
2 m

3
de modo que tomaria unos 300.000 años antes de que la incerteza en la posicion crezca hasta 1
milimetro! Podemos entonces tratar una particula macroscopica en forma clasica por cualquier
lapso de tiempo razonable.

Algunas caracteristicas generales de autofunciones de energia

Consideremos ahora las autofunciones de energia para algun potencial V (x). A partir de la ecuacion
de Schrodinger estas obedecen
2m[E ¡ V (x)]
(x) = ¡
00
(x)
~2

Preguntemonos ahora que implica la continuidad en V (x). Comencemos en un punto x0 donde

y 0 tienen valores (x0) y 0(x0). Si hacemos de cuenta que x es una variable temporal y es
una coordenada de la particula, el problema de encontrar en todo punto es como encontrar la
trayectoria de una particula (para todo tiempo) dadas su posicion y velocidad en algun tiempo
y su aceleracion como una funcion de la posicion y el tiempo. Es claro que si integramos estas
ecuaciones vamos a obtener 0(x) y (x) continuas. Esta es la situacion tipica (y lo que ocurre en
realidad). Hay, sin embargo, algunas situaciones donde, por simplicidad matematica, consideramos
potenciales que cambian en forma abrupta en algun punto. Esto significa que 00 salta en forma
abrupta en ese punto. Sin embargo, 0 sera continua, ya que el area bajo una funcion es continua
incluso si la funcion salta un poco. ¾Que ocurre si el cambio en V es infinitamente grande?
Significa que 00 es tambien infinitamente grande. Esto, a su vez, implica que 0 puede cambiar
abruptamente a medida que cruzamos este punto, ya que el area bajo 00 puede ser finita en una
region infinitesimal que rodee este punto. Pero ya sea que 0 sea continua o no, , que es el area
bajo 0, sera siempre continua. Ahora veremos algunos casos especificos.

2 Particula en una caja

Ahora vamos a considerar nuestro primer problema con un potencial, si bien a primera vista puede
resultar un tanto artificial (aunque no tanto! Ver por ejemplo sistemas del tipo gases de electrones
bidimensionales)

0 jxj < L/2
V (x) =
1 jxj > L/2

Este potencial (ver Fig. 1) se llama caja ya que hay una barrera de potencial infinita para la
particula si esta intenta salir de la region jxj < L/2.

Figura 1. a) El potencial tipo caja. b) Los primeros dos niveles y funciones de onda de la caja.

La ecuacion de autovalores en la base x es

2m
00
+ (E ¡ V ) = 0
~2

Vamos a comenzar separando el espacio en tres regiones 1, 2, y 3. La solucion es llamada 1,

2, y 3 en cada region.

4
Consideremos primero la region 3, en la cual V = 1. Conviene primero no tomar el caso infinito
sino pensar que V tiene un cierto valor V0 mayor que E. De esta forma, la ecuacion anterior es

2m
00
3 ¡ (V0 ¡ E) 3 = 0
~2
cuya solucion es

3 = Ae¡x + Be+x;  = [2m(V0 ¡ E)/~2]1/2

Si bien A y B son constantes arbitrarias desde el punto de vista matematico, debemos tomar
B = 0 desde el punto de vista fisico ya que Be+x explota exponencialmente cuando x ! 1 y tales
funciones no pertenecen a nuestro espacio de Hilbert. Si ahora tomamos V0 ! 1 y tenemos en
cuenta que en esta region x > 0, entonces

3 = 0

Puede verse lo mismo en la region 1, es decir 1 = 0. En la region 2, dado que V = 0, las soluciones
coinciden exactamente con aquellas de la particula libre, es decir:

2 = Ae+ikx + Be¡ikx; k = (2mE /~2)1/2

Por lo tanto pareciera que los autovalores de energia son continuos nuevamente, como en el caso
de la particula libre. Esto, de todos modos, no es asi, ya que la igualdad entre las autofunciones
de energia y de momento solo se da en la region 2. Debemos pedir que 2 vaya en forma continua
hacia sus contrapartes 1 y 3 a medida que cruzamos hacia las regiones 1 y 3, respectivamente.
En otras palabras, exigimos que
1(¡L/2) = 2(¡L/2)
2(+L/2) = 3(+L/2)

Por otro lado, no podemos pedir que 0 sea continua en las paredes de la caja pues V va a infinito
alli. Si ahora expresamos en forma explicita estas ecuaciones de continuidad, se tiene

0 = Ae¡ikL/2 + Be+ikL/2
Ae+ikL/2 + Be¡ikL/2 = 0

que, en forma matricial, adquiere la forma

!   
e¡ikL/2 e+ikL/2 A 0
=
e+ikL/2 e¡ikL/2 B 0

Esta ecuacion solo tiene soluciones no triviales si el determinante de la matriz se anula, es decir

e¡ikL ¡ e+ikL = ¡2i sin(kL) = 0

lo cual implica que

n
k = ; n = 0; 1; 2; :::
L

Para encontrar las autofunciones correspondientes, volvemos a las ecuaciones de continuidad. Dado
que solo una de estas ecuaciones es independiente, podemos utilizar la primera, es decir

Ae¡in/2 + Be+in/2 = 0 ! A = ¡einB = ¡(¡1)nB

5
por lo que se obtienen dos familias de soluciones:
8 
>  
>
> 2 1/2 nx
< sin n par
L L
n(x) =  1/2  
>
> 2 nx
>
: cos n impar
L L

donde obviamente estamos considerando jxj < L/ 2. Notar que el caso 0 no es solucion, ya que
0 = 0 y esta solucion no corresponde a ningun caso fisico. Por otro lado, ¡n = ¡(¡1) n, y
n

como las funciones de onda (en general estados) que solo difieren en una fase global son fisicamente
equivalentes, debemos restringirnos a n > 0.
Los autovalores de energia pueden obtenerse haciendo el reemplazo en la definicion de k, con lo
que se tiene
 
~2kn2 ~2 n 2
En = = ; n = 1; 2; :::
2m 2m L

Tenemos entonces nuestro primer encuentro con la cuantizacion de una variable dinamica. Las
variables consideradas hasta ahora, x y p, tenian un espectro continuo de autovalores que va desde
¡1 hasta 1, el cual coincide con los valores permitidos en mecanica clasica. De hecho, lo mismo
ocurria para el espectro de energias en el caso de la particula libre. La particula en la caja es el
ejemplo mas simple de una situacion que encontraremos de aqui en adelante donde la ecuacion
de Schrodinger, combinada con condiciones de borde apropiadas, da lugar a la cuantizacion de la
energia.

Energia del punto cero

En base a lo obtenido, puede verse que la menor autoenergia no es cero (como ocurriria en el caso
clasico, correspondiente a la particula en reposo dentro de la caja) sino E1 = ~2 2 /2mL2. La razon
detras de esto es el principio de incerteza, el cual previene que la particula, cuya posicion (y por lo
tanto
x) esta confinada a jxj < L/2, de tener momento nulo bien definido. Esto, a su vez, implica
una cota inferior en la energia, que ahora vamos a deducir. Comencemos con

^ = p^2
H
2m
de modo que

^ i = hp^ i
2
hH
2m

Ahora bien, hp^i = 0 por la siguiente razon: Como un autoestado dentro de la caja es un estado
estacionario, hp^i es independiente del tiempo. Si hp^i =
/ 0 entonces la particula debe (en el sentido
de promedio) moverse a la izquierda o a la derecha y eventualmente escapar de la caja, lo cual no
puede ocurrir fisicamente. Por lo tanto, podemos escribir lo anterior como

^ i = h(p^ ¡ hp^i) i = (
p)

2 2
hH
2m 2m

De modo que si ahora usamos la relacion de incerteza

~

x
p >
2
se tiene

^i > ~2
hH
8m(
x)2

6
Dado que la variable x esta restringida a ¡L/2 6 x 6 L/2, su desviacion estandar no puede exceder
L/2. Consecuentemente

^i > ~2
hH
2mL2

^ i = E, de modo que
En un autoestado de energia, hH

~2
E >
2mL2

Si bien la energia del estado fundamental ocurre ser  2 veces mayor que esta cota, esto nos dice
que nunca puede haber una energia menor que un cierto valor finito. La diferencia entre la menor
energia posible que un sistema cuantico puede tener y la minima energia clasica posible es lo que
se conoce como energia de punto cero, y tiene consecuencias profundas en diversas ramas de la
fisica, incluyendo fisica de altas energias y cosmologia.

Estados ligados
Las soluciones de la caja son tambien ejemplos de estados ligados, a saber, estados en los cuales
un potencial previene que la particula escape al infinito. Los estados ligados estan entonces carac-
terizados por
jxj!1
(x) !
!
!
!!
!
!
!!
!
!
!!
!
!
!! 0
!
!
!

Estos aparecen en mecanica cuantica exactamente donde son esperados clasicamente, es decir, en
situaciones donde V (1) es mayor que E.
Los niveles de energia de los estados ligados estan siempre cuantizados. Veamos porque esto ocurre.
En el problema de la caja, la cuantizacion resulta del requisito de que 2 complete un numero
entero de semi-ciclos dentro de la caja de modo que se conecte suavemente con sus contrapartes 1
y 3 que eran exactamente cero. Consideremos ahora una particula confinada en un pozo finito,
esto es, un potencial que salta de 0 a V0 en jxj = L/2,


0 jxj < L/2
V (x) =
V0 jxj > L/2

Ya hemos visto que en la region clasicamente prohibida (E < V0; jxj > L/2) es la suma de una
exponencial creciente y otra decreciente a medida que jxj ! 1 y que debemos tomar el coeficiente
de la exponencial creciente como cero para obtener una solucion fisicamente admisible. Esto nos da

1 = A1 e+x; 3 = B3 e¡x;  = [2m(V0 ¡ E)/~2]1/2

En la region clasicamente admisible (jxj < L / 2) es la suma de dos exponenciales complejas

(equivalente a un seno y un coseno):

2 = A2 e+ikx + B2 e¡ikx; k = (2mE /~2)1/2

Dado que V es finito en todos lados, podemos pedir que y 0

sean continuas en las interfases.
Tenemos entonces cuatro condiciones sobre ,

1(¡L/2) = 2(¡L/2) 1(¡L/2) = 2(¡L/2)

0 0

2(+L/2) = 3(+L/2) 2(+L/2) = 3(+L/2)

0 0

7
Para puede parecer que tenemos 4 coeficientes a determinar: 2 de las exponenciales decrecientes y
2 de las exponenciales complejas. Sin embargo, la escala de es irrelevante tanto para la ecuacion
de autovalores como para las condiciones de contorno, ya que son lineales en y 0. En todo caso,
la escala de puede tomarse a partir de la condicion de normalizacion, que agrega una ecuacion
mas al conjunto anterior. En este caso se tiene que las soluciones deben satisfacer
8
>
< k tan kL = + solucion par
>
2 (3)
>
> kL
: k cot = ¡ solucion impar
2
8
>  
>
< tan r = + solucion par
r r  r (4)
>
>
: cot r = ¡ solucion impar
r r r

Donde la paridad de la solucion proviene de la forma par del potencial y se vera en breve. Estas
ecuaciones no admiten como solucion cualquier valor para la energia, sino solo ciertos valores
especificos, como se puede ver en la figura 2. Este es el origen de la cuantizacion de la energia aqui.

Figura 2. Soluciones de la ecuacion (3). En esta representacion, los parametros vienen dados como:
= kL/2, = L/2 y r = (mV0L2 /2~2)1/2:

Esta construccion puede generalizarse para un potencial V (x) que tiende a los limites V cuando
x ! 1, respectivamente, y que atrapa a una particula de energia E (menor que V). A partir
de separar el espacio en intervalos, uno puede argumentar que se tiene una condicion mas que el
numero de coeficientes a determinar. La clave esta en notar que los intervalos extremos se vuelven
infinitamente largos, una vez que se alcanzaron los valores limite V. En los intervalos infinitos
se tienen 2 coeficientes (como antes, se pide que la solucion no explote), mientras que en la parte
interna se tienen 2 coeficientes por cada intervalo. De este total debemos restar un coeficiente
porque, como dijimos, la escala global de es irrelevante. Tenemos entonces que si N es el numero
de intervalos, hay 2(N ¡ 1) ¡ 1 coeficientes a determinar y 2(N ¡ 1) ecuaciones. De modo que,
como antes, solo podemos esperar soluciones para valores especiales (discretos) de la energia.
Es claro por qué no hay cuantizacion de energia cuando el estado no es ligado: dado que la particula
es clasicamente permitida en el infinito, oscila alli (es decir, su amplitud es finita) y entonces se
tienen dos parametros mas, uno de cada lado, y por lo tanto dos soluciones para cualquier energia.

8
Paridad definida en autofunciones de energia
No es casual que las soluciones anteriores presenten paridad definida, es decir que sean o bien pares
o bien impares. Esto es una consecuencia de la propiedad par del potencial, es decir V (¡x) = V (x).
Para formalizar esta conclusion, introducimos el operador  ^ conocido como operador paridad, es
decir

^ jxi = j¡xi


cuyo rol es el de invertir el signo de la coordenada. Veamos sus propiedades: este operador es
Hermitico, pues
^ jx 0i
hxj = hxj¡x 0i = (x 0 + x)

^ yjx 0i = (hx 0j

hxj ^ jxi) = (hx 0j¡xi) = (x 0 + x)

Es su propia inversa

^
 ^ jxi = 
^ j¡xi = jxi ! 
^ 2 = ^1

y por lo tanto es unitario. Si queremos ver como actua a izquierda,

^ = (
hxj ^ yjxi) y = (
^ jxi) y = (j¡xi) y = h¡xj

Para un estado cuantico general, tenemos entonces

^ j i ! hxji = hxj
ji =  ^ j i = h¡xj i

es decir

(x) = (¡x)

Sus autoestados pueden obtenerse a partir de la ecuacion de autovalores,

^ j i =  j i ! 
 ^ 2j i =  2j i

^2 = 1
pero como  ^, tenemos que  = 1. Entonces si tomamos

^ j+i = +hxj+i
 = +1 ! hxj
h¡xj+i = +hxj+i ! funciones pares

^ j¡i = ¡hxj¡i
 = ¡1 ! hxj
h¡xj¡i = ¡hxj¡i ! funciones impares

Es decir que las autofunciones del operador paridad son aquellas funciones con paridad definida,
y sus autovalores correspondientes son 1 dependiendo si son pares o impares.
Tambien se puede ver como actua este operador sobre los autoestados de momento:
Z 1 Z 1 Z 1
^ jpi =
 ^ jxihxjpi =
dx  dx j¡xihxjpi = dx jxih¡xjpi
¡1 ¡1 ¡1

9
y como
1
h¡xjpi = e¡ipx/~ = hxj¡pi
(2~)1/2
se tiene entonces
Z 1
^ jpi =
 dx jxihxj¡pi = j¡pi
¡1

Ahora veamos como transforma este operador al operador posicion x^

^ x^
hxj ^ j i = ¡xhxj i = ¡hxjx^ j i ! 
^ x^
^ = ¡x^

a su vez,

^ x^
 ^ 2 = ¡x^
^ ! [
^ ; x^] = 
^ x^ ¡ x^
^ = ¡2x^
^

lo que nos dice que estos operadores no conmutan. Veamos que ocurre con el cuadrado

^ x^2
hxj ^ j i = x2hxj i = hxjx^2j i ! 
^ x^2
^ = x^2

y entonces

^ x^2
 ^2 = 
^ x^2 = x^2
^ ! [
^ ; x^2] = 0

^ actua de igual forma en jpi y en jxi, entonces se tiene

Lo mismo ocurre con p^ y p^2, puesto que 

^ ; p^] = ¡2p^
[ ^
^
[; p^ ] = 0
2

Ahora que conocemos las propiedades de este operador, podemos pensar en problemas donde
el potencial es par en la coordenada. Debido a la paridad, dicha funcion puede ser escrita como
potencias pares de x^, y entonces

^ ; V (x^)] = 0
[

Esto, sumado a que [ ^ ; p^2] = 0, implica que el Hamiltoniano de esta clase de problemas tambien
conmuta con el operador paridad. Por lo tanto, dichos operadores comparten un conjunto de
autovectores y, mas especificamente, un conjunto de autofunciones. Como las autofunciones de  ^
poseen paridad definida, lo mismo ocurre con las autofunciones del Hamiltoniano. Puede ocurrir,
como es el caso de la particula libre, que el espectro de autovalores de energia sea degenerado
(como hemos visto antes, cuando asociamos a las autofunciones con ondas planas). En este caso,
uno puede tener autofunciones del Hamiltoniano que no presenten paridad definida, pero siempre
es posible considerar combinaciones lineales apropiadas que respeten dicha paridad.

Barrera de potencial
Ahora vamos a considerar una situacion donde el potencial no es ligante sino que presenta una
barrera para el paso de la particula. Por simplicidad, vamos a considerar el siguiente potencial
constante a trozos:

V0 jxj 6 L/2
V (x) =
0 jxj > L/2

10
que puede verse en la figura 3.

Figura 3. Esquema del potencial tipo barrera.

Nuevamente, podemos separar el espacio en tres regiones, donde conocemos las soluciones par-
ticulares. A su vez, vamos a considerar en simultaneo los casos E < V0 y E > V0. En la regiones
externas a las barreras, se tiene

i = aie+ikix + bie¡ikix; ki = (2m(E ¡ Vi)/~2)1/2

donde hemos tomado V1 = V3 = 0 y V2 = V0. Claramente, en la region 2 pueden aparecer como solu-
ciones exponenciales reales o complejas dependiendo del valor de E. Como antes, debemos ahora
considerar las condiciones de continuidad tanto en las funciones de onda como en sus derivadas.
Sin embargo, a diferencia del caso del pozo, tenemos ahora un numero mas grande de coeficientes
a determinar que ecuaciones, lo que implica que algunos coeficientes pueden quedar sin determinar.
En esta clase de problemas de dispersion, lo que uno desea describir, tipicamente, son las relaciones
entre los coeficientes asociados a las regiones externas (en este caso 1 y 3). Para ello, podemos
plantear las ecuaciones de continuidad en cada interfase por separado, que matricialmente toma
la forma
M21v1 = M12v2
M32v2 = M23v3
donde
 
ai
vi =
bi

y ahora podemos despejar los coeficientes internos, es decir

v2 = (M12)¡1M21v1 ! v3 = (M23)¡1M32(M12)¡1M21v1

esto nos permite entonces expresar a los coeficientes de la region 3 en terminos de los coeficientes
de la region 1. En nuestro caso se tiene

! !
e¡ik1L/2 +e+ik1L/2 e+ik2L/2 +e¡ik2L/2
M21 = ; M32 =
ik1e¡ik1L/2 ¡ik1e+ik1L/2 ik2e+ik2L/2 ¡ik2e¡ik2L/2
! !
e¡ik2L/2 +e+ik2L/2 e+ik3L/2 +e¡ik3L/2
M12 = ¡ik2L/2
; M23 =
ik2e ¡ik2e+ik2L/2 ik3e +ik3L/2
¡ik3e¡ik3L/2

11
de modo que
   
a3 a1
= M ! M = (M23)¡1M32(M12)¡1M21
b3 b1

Ahora bien, dado que estamos considerando un problema de dispersion, podemos pensar en la
situacion donde hay una onda incidente desde la izquierda con amplitud 1 (sin perdida de gene-
ralidad), la cual se refleja con amplitud r, y se transmite con amplitud t. Esto implica que en la
ecuacion matricial de arriba uno tiene
   
t 1
= M
0 r
de modo que
 (
t = M11 + M12r r = ¡M21 /M22
!
0 = M21 + M22r t = M11 ¡ M12M21 /M22

tomando los reemplazos k1 = k3 = k y k2 = , esto entonces resulta

e¡ikL(1 ¡ e2iL)(k 2 ¡ 2)

r =
(k + )2 ¡ e2iL(k ¡ )2

4ke¡i(k ¡)L
t =
(k + )2 ¡ e2iL(k ¡ )2

Ahora debemos tener cuidado de los casos en energia, ya que si E < V0 entonces  es imaginario,
mientras que si E > V0 entonces  es real. Tenemos pues:
 
2mjE ¡ V0j 1/2
Caso 1: E < V0 !  = i = i0
~2

Entonces las probabilidades de reflexion y transmision son

(1 ¡ e¡20L)2(k 2 + 20)2
R = jrj2 =
(k + 0) + e
2 2 2 ¡40L (k 2 + 20)2 ¡ e¡20L[2k 2 + 220 ¡ 12(k0)2]

16(k0)2e¡20L
T = jtj2 =
(k 2 + 0)2 + e¡40L(k 2 + 20)2 ¡ e¡20L[2k 2 + 220 ¡ 12(k0)2]
2

 1/2
2mjE ¡ V0j
Caso 2: E > V0 !  = = 0
~2

con lo cual las probabilidades de refleccion y transmision son

2[1 ¡ cos(20L)](k 2 ¡ 20)2

R = jrj2 =
(k + 0)4 + (k ¡ 0)4 ¡ 2cos(20L)(k 2 ¡ 20)2

16(k0)2
T = jtj2 =
(k + 0)4 + (k ¡ 0)4 ¡ 2cos(20L)(k 2 ¡ 20)2

y en ambos casos puede verse que

R+T = 1

12
que es una consecuencia de la conservacion de la probabilidad. Esto puede verse en la figura 4,
donde mostramos las probabilidades de reflexion y transmision en funcion de la energia.

Figura 4. Probabilidad de refleccion y transmision en el problema de la barrera.

Es sorprendente notar que, incluso a valores E < V0, existe una probabilidad no nula de que
la particula traspase la barrera. Este efecto, conocido como efecto tunel , es una caracteristica
netamente cuantica, y puede verse tambien en la solucion que hemos construido (ver figura 5).

Figura 5. Partes real (rojo) e imaginaria (azul) de la funcion de onda en el problema de dispersion. El
valor utilizado de la energia es E = 0.8V0.

Aqui podemos ver claramente como la amplitud de hallar a la particula cambia en funcion de la
posicion, al punto que existe una amplitud finita en la region 3, contrario a lo que uno esperaria
en mecanica clasica.

13
3 La ecuacion de continuidad para la probabilidad
Ahora vamos a interrumpir nuestra discusion de problemas unidimensionales para incorporar dos
conceptos que seran usados en futuras discusiones, a saber, el de densidad de corriente de probabi-
lidad y la ecuacion de continuidad que satisface. Dado que el concepto de corriente de probabilidad
tambien sera usado en problemas tridimensionales, discutimos aqui el caso de una particula en tres
dimensiones.
Como preludio de nuestro estudio de la ecuacion de continuidad en mecanica cuantica, recordemos
la ecuacion analoga de electromagnetismo. Sabemos que en ese caso la carga total en el universo
es una constante. Esto es un ejemplo de una ley de conservacion global, ya que se refiere a toda
la carga presente en el universo. Pero la carga tambien se conserva localmente, un hecho que
usualmente se expresa en la forma de una ecuacion de continuidad:

@
(r; t) = ¡r
j
@t

donde  y j son las densidades de carga y corriente, respectivamente. Integrando esta ecuacion
sobre un volumen V rodeado por una superficie SV tenemos, luego de usar la ley de Gauss:
Z Z Z
d
(r; t)d r = ¡
3
r
jd r = ¡
3
j
dS
dt V V SV

Esta ecuacion establece que cualquier disminucion en la carga en el volumen V es tenida en cuenta
por el flujo de la carga hacia fuera de este. En otras palabras, la carga no se crea ni destruye en
ningun volumen.
En mecanica cuantica la cantidad que se conserva globalmente es la probabilidad total de encontrar
a la particula en cualquier lugar del universo. Podemos obtener este resultado expresando la
invarianza de la norma en la base de coordenadas, ya que:

h (t)j (t)i = h (0)jU^ y(t)U^ (t)j (0)i = h (0)j (0)i

entonces
Z Z
const = h (t)j (t)i = h (t)jrihrj (t)id r = 3
P (r; t)d 3r

Para obtener la analogia con el caso local, debemos volver a la ecuacion de Schrodinger:

@ ~2 2
i~ = ¡ r +V
@t 2m
y su conjugada, es decir
@  ~2 2
¡i~ = ¡ r 
+V 
@t 2m

Notar que V debe ser real para que el Hamiltoniano sea Hermitico. Multiplicando la primer
ecuacion por  y la segunda por y tomando la diferencia se tiene

@ ~2
i~ ( 
) = ¡ ( 
r2 ¡ r2 
)
@t 2m

@P ~
! = ¡ r
( 
r ¡ r 
)
@t 2mi

@P
! = ¡r
j
@t

14
donde la densidad de corriente de probabilildad viene dada por

~
j = ( 
r ¡ r 
)
2mi

es decir, el flujo de probabilidad por unidad de tiempo y area perpendicular a j. Para reobtener
a la ecuacion de conservacion global, integramos la ecuacion de continuidad sobre todo el espacio:
Z Z
d
P (r; t)d 3r = ¡ j
dS
dt S1

donde S1 es la esfera en el infinito. Para el caso tipico de funciones

R  2 de onda que estan normalizadas
a la unidad, r3/2 ! 0 a medida que r ! 1 de modo que r d esta acotado, y la integral de
superficie de j en S1 se anula.

Ejemplos:
Para una onda plana que se propaga en el espacio:

(r; t) = hrj (t)i = Aei(k
r ¡!t)

la densidad de probabilidad es constante en todo punto, es decir

P (r; t) = h (t)j(jrihrj)j (t)i = 

(r; t) (r; t) = jAj2

confirmando que las ondas planas son estacionarias. Sin embargo, la densidad de corriente de
probabilidad no es cero,

~ ~ p
j(r; t) = ( 
r ¡ r 
) = jAj2 k =  = v
2mi m m

lo que ilustra que la particula puede estar en movimiento incluso si su densidad de probabilildad
espacial no posee dependencia temporal.
Para el caso de una particula en una caja, confinada en la region 0 < x < L, tenemos que las
autofunciones de energia son de la forma
r  
2 n
n(x) = hxjni = sin x
L L
de modo que

jn = 0

ya que n = n . Notar entonces que este resultado tambien se aplica a cualquier funcion de onda
que sea puramente real.

4 Algunos teoremas a tener en cuenta

Volviendo ahora al caso unidimensional, presentamos en esta seccion algunos teoremas que pueden
resultar de utilidad.
Teorema 1: No hay degeneracion en estados ligados unidimensionales.

15
Prueba: Sean 1 y 2 dos soluciones con igual autovalor E, es decir

~2 d 2
¡ 1+V (x) 1 = E 1
2m dx2
~2 d 2
¡ 2+V (x) 2 = E 2
2m dx2

Multiplicamos entonces la primer ecuacion por 2 y la segunda por 1 y luego restamos ambas
ecuaciones, de modo que se tiene
d2 d2
1 2¡ 2 1 = 0
dx2 dx2
o bien
 
d d d
1 2¡ 2 1 = 0
dx dx dx
tal que
d d
1 2¡ 2 1 = c
dx dx

Para encontrar la constante c, uno puede ir a jxj ! 1, donde 1 y 2 se anulan, de modo que uno
tiene c = 0. Por lo tanto
1 1
d 1 = d 2 ! 1 = ed 2
1 2

donde d es una constante, y entonces las dos autofunciones solo difieren por un factor de escala y
representan el mismo estado.

Teorema 2: Las autofunciones de H pueden siempre elegirse reales puras en la base de coorde-
nadas.
Prueba: Si n es solucion de la ecuacion de Schrodinger:
 
~2 d 2
¡ + V (x) n = En n
2m dx2

entonces podemos tomar su complejo conjugado, es decir

 
~2 d 2
¡ + V (x) 
n = En 
n
2m dx2

de modo que n y n son autofunciones con el mismo autovalor, con lo cual, podemos elegir como
solucion la parte real o la parte imaginaria, pues se escriben como combinacion lineal de estas dos,
es decir
n+
 
n n¡ n
r = ; i =
2 2i

y estas soluciones son tambien autofunciones con energia E.

Es importante notar que este teorema es valido en dimensiones superiores siempre y cuando el
Hamiltoniano, ademas de ser Hermitico, es tambien real.1

1. Notar, por ejemplo, que un campo magnetico introduce terminos no reales en el Hamiltoniano, y por lo tanto el
teorema deja de ser valido.

16
Si ahora aplicamos este teorema a estados ligados unidimensionales, y tenemos en cuenta el teorema
anterior, se tiene que

i = c r

donde c es una constante. Consecuentemente,

= r+i i = (1 + ic) r

Dado que el prefactor es una constante irrelevante, podemos tomar autofunciones reales sin perdida
de generalidad, y entonces podemos decir que bajo esta eleccion los estados ligados no poseen
densidad de corriente.

R
Teorema 3: Si (1 + x2)jV (x)jdx < 1, entoncesRla ecuacion ¡d 2 /dx2 + V (x) admite al menos
un estado ligado para cualquier  > 0 si y solo si V (x)dx 6 0.
La prueba de este teorema se basa en el metodo variacional, que veremos mas adelante, asi que
pospondremos su demostracion hasta ese entonces.
Probemos entonces el ejemplo de la particula en el pozo, donde el potencial viene dado como


¡V0 jxj 6 L/2
V (x) =
0 jxj > L/2

Primero hagamos conversion de unidades para llegar a las formas usadas en el enunciado:

p
V 2 mc
V ! V~ = x ! x~ = x
mc2 ~

de modo que la ecuacion de Schrodinger se escribe

 
d2
¡ + V~ = E~
dx~2

y podemos ahora volver a x~ ! x y V~ ! V teniendo en cuenta que son cantidades adimensionales.

Asi pues, primero debemos ver si
Z
(1 + x2)jV (x)jdx < 1

que, en nuestro caso, es

Z `
2
(1 + x2)v0 dx = `v0 + `3v0 < 1
¡` 3

p
donde ` = mcL/~ 2 y v0 = V0 /mc2. A su vez,
Z 1
V (x)dx = ¡v0` < 0
¡1

y entonces queda claro que en este caso hay al menos un estado ligado. Notar que las unidades
introducidas son simplemente para ser consistentes con las unidades del enunciado, pero no tienen
un correlato fisico.

17
5 Potencial delta de Dirac
Continuando con los problemas unidimensionales, ahora vamos a considerar el caso de un potencial
tipo delta de Dirac, que puede escribirse como

V (x) = (x)

donde  puede ser positivo (barrera delta) o negativo (pozo delta). Vamos primero a considerar el
problema donde el potencial es atractivo, es decir el caso donde  < 0.
Si tomamos E < 0 tenemos posibilidad de hallar estados ligados, mientras que si tomamos E > 0
tendremos un problema de dispersion similar al del pozo.
Caso E < 0: A partir del teorema anterior, es claro que este potencial tiene al menos un estado
ligado. Nuevamente, dividimos el espacio en dos regiones, dadas como x < 0 y x > 0, donde al ser
nulo el potencial podemos utilizar las autofunciones de momento, es decir

r
2mjE j
1 = ae ; x
2 = be ¡x
;  =
~2

El problema aqui, a diferencia de otros potenciales, es que ahora estamos en un caso donde el
potencial es infinito en un punto, lo que implica que la derivada de la funcion de onda puede ser
discontinua en este punto. Tenemos entonces que averiguar cual es el salto en esta derivada. La
continuidad de la funcion de onda, implica a = b. Para ver la derivada, planteamos la ecuacion de
Schrodinger, que en este caso resulta

~2 d 2
¡ (x) + (x) (x) = E (x)
2m dx2

Ahora integramos esta ecuacion en un rango infinitesimal alrededor de x = 0, es decir

Z  Z  Z 
~2 d2
¡ (x)dx + (x) (x)dx = E (x)dx
2m ¡ dx
2
¡ ¡

lo que toma la forma

Z 
~2
¡ [ () ¡
0
(¡)] +  (0) =
0
E (x)dx
2m ¡

Si ahora aplicamos el limite  ! 0 se tiene que el termino a la derecha es nulo (debido a la

continuidad de ) y obtenemos
2m
2(0) ¡ 1(0) = (0)
0 0
~2

Este es justamente el salto en la derivada que queriamos encontrar. Si ahora reemplazamos por
las funciones propuestas, se tiene
mjj
 =
~2

con lo cual se tiene un unico estado ligado, dado por

(x) = ae¡jxj

18
con energia
m2
E = ¡
2~2

Caso E > 0: En esta situacion, podemos preguntarnos por las amplitudes de probabilidad de
reflexion y transmision dada una onda incidente desde la izquierda. Las soluciones entonces son
de la forma
r
2mE
1 = e + re¡ikx; 2 = te ; k =
ikx ikx
~2

La continuidad de la funcion de onda implica que t = 1 + r, y el salto en la derivada implica

2m 1 1
2(0) ¡ 1(0) = (0) ! r = ~2k t =
0 0
;
~2 i m ¡ 1
m
1 + i ~2k

de modo que las probabilidades son

1 1
R = jrj2 = 2~2E
; T = jtj2 = m2
1+ m2
1 + 2~2E

Si consideramos el caso de la barrera, es decir  > 0, esta claro que solo podemos tomar E > 0 y en
este caso tenemos nuevamente un problema de dispersion, donde es facil ver que las probabilidades
de reflexion y transmision no cambian respecto a lo anterior.

Figura 6. Probabilidades de reflexion (rojo) y transmision (azul) para un pozo tipo delta de Dirac. La
energia se da en unidades m2 /2~2.

Puede verse entonces (figura 6) como inicialmente toda la onda se refleja, que corresponde al caso
clasico, pero a medida que incrementamos la energia, una parte comienza a transmitirse, al punto
que a la energia m2 /2~2 ambas probabilidades son iguales, y para valores de E mayores la cantidad
transmitida es mayor a la reflejada.

19