Notas sobre Mecánica Cuántica

Mauricio Velásquez
 Índice general

1. El oscilador armónico 1
1.1. Las funciones de onda . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.1.1. Función de onda del estado fundamental en la representación posición . . . . . . . . . . . . . . . . . . . 4
1.1.2. Segundo método: resolver la ecuación de Schrödinger . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
1.2. Funciones hipergeométricas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
1.3. Pozo de potencial infinito . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
1.4. Escalón de potencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
1.5. La barrera de potencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22
1.6. La barrera de potencial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 27
1.7. Casos límites . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
1.8. Un electrón en un pozo de potencial con barrera finitas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32
1.8.1. Pozo de potencial con una de las paredes finita . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32

2. Ondas, operadores y postulados de la mecánica cuántica 37

2.1. Vista clásica del experimento de la doble rendija . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
2.2. Vista cuántica de partículas y ondas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 40
2.3. La función de onda . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 43
2.4. Principio de superposición lineal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
2.5. Paquetes de onda . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
2.6. Paquetes de onda localizado . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45
2.7. Paquetes de onda y las relaciones de incertidumbre . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2.8. Espacios de Hilbert . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 49
2.9. Dimensión y base de un espacio vectorial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
2.10. Funciones de cuadrado integrable: funciones de onda . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 50
2.11. Notación de Dirac . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51
2.12. Operadores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 53
2.12.1. Propiedades de los conmutadores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54
2.12.2. Adjunto de un operador . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
2.12.3. Operadores hermíticos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
2.13. Autofunciones y autovalores de operadores . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
 ii

;A<
2.13.1. Operadores hermíticos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
2.13.2. Operador proyección . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
2.13.3. Operador unitario . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
2.14. Algunos teoremas importantes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57
2.15. Mecánica ondulatoria y matricial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
2.15.1. Mecánica matricial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 59
2.16. La función delta de Dirac . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60
2.16.1. Propiedades de la función delta . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 62
2.17. La función delta tridimensional . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 63
2.18. Representación posición . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 64
2.19. Representación momentum . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 65
2.20. Mecánica ondulatoria . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
2.21. Postulados de la mecánica cuántica . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 67
2.22. Operador evolución . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
2.22.1. Propiedades del operador evolución . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68

Apéndices 70
 Capítulo

El oscilador armónico

E
l oscilador armónico es un problema muy importante en todas las ramas de la Física. Suministra un útil modelo
para una variedad de fenómenos vibracionales que se presentan en la naturaleza, por ejemplo: en mecánica
clásica, electrodinámica, mecánica estadística, estado sólido, física de partículas, física nuclear y física atómica.
En mecánica cuántica ilustra los conceptos básicos y el formalismo teórico. El hamiltoniano de una partícula de
masa m la cual oscila con frecuencia angular ω bajo la influencia de un potencial armónico unidimensional es:
2
b = pb + 1 mω 2 xb2
H (1.1)
2m 2
El problema es como encontrar los autovalores de energía y autoestados de este hamiltoniano. Existen varios métodos para
resolver este problema: el método de resolver la ecuación de Schrödinger independiente del tiempo.
Otro método sería el método algebraico. Este método trata con el álgebra de operadores, que involucra operadores algebraicos
conocidos como operadores de creación y destrución. Este método de operadores será el que se desarrollará.
Se escribe el hamiltoniano en (1.1) en términos de dos operadores adimensionales hermíticos:

Pb √
pb = √ y qb = X
b mω/h̄
mh̄ω

Hb = h̄ω ( pb2 + qb2 ). (1.2)

2
Se definen dos operadores adimensionales no hermíticos:

1 1
a = √ (qb + i pb), b
b a+ = √ (qb − i pb)
2 2
+ 1
a = (qb − i pb)(qb + i pb)
a +b
b
2
1
= (qb2 + i qbpb − i pbqb + pb2 )
2
1 i
= (qb2 + pb2 ) + (qbpb − pbqb)
2 2
1 i
= (qb2 + pb2 ) + [qb, pb] (1.3)
2 2
 2
;A<

Ahora, se desarrolla el conmutador entre qb y pb.

" r #
mω Pb
[qb, pb] = X
b ,√
h̄ mωh̄
" #
b Pb PbX
X b 1 b b bb
= − = [X P − PX]
h̄ h̄ h̄

=i
Aquí utilizamos el resultado [ X,
b Pb] = ih̄
Entonces:
1 2 1
a+ + b
b a= (qb + pb2 ) −
2 2
o
1 2 1
a+ + b
(qb + pb2 ) = b a+ (1.4)
2 2
Insertando (1.4) en (1.2) se tiene:
   
+ 1 1
H = h̄ω b
b a +b
a+ = h̄ω N +
b , con N a+ + b
b =b a (1.5)
2 2

Donde Nb es conocido como el operador número, el cual es hermítico.

Ahora se halla el conmutador [b a+ ]
a, b
 
1 1
[b a + ] = [b
a, b a+ − b
ab a+ b
a] = (qb + i pb)(qb − i pb) − (qb − i pb)(qb + i pb)
2 2
1 2
[b a+ ] =
a, b pb2 − 
[qb − i qbpb + i pbqb +  qb2 − i qbpb + i pbqb − 
pb2 ]
2
1
= [−2i qbpb + 2i pbqb]
2
= i [ pbqb − qbpb] = i [ pb, qb] = −i [qb, pb] = 1

[b a+ ] = 1,
a, b (1.6)

[ H,
b Nb] = 0

esto implica que poseen un conjunto de autoestados comunes definidos como |ni, tal que:

b | n i = n | n i;
N
(1.7)
b |ni = En |ni.
H

Los autoestados |ni son llamados autoestados de energía. Combinando (1.5) y (1.7), se obtiene:
   
1 1
H |ni = h̄ω N +
b b |ni = h̄ω n + |ni
2 2
 
1
En = h̄ω n +
2

Sean los conmutadores:

[b b ] = h̄ω b
a, H a+ , H
a, [b a+
b ] = −h̄ω b

Hb (ba|ni) = (b b − ωb
aH a)|ni = ( En − h̄ω )(b
a|ni)
+ + +
a |ni) = (b
H (b
b a H + ωb
b a+ |ni)
a )|ni = ( En + h̄ω )(b
 3
;A<

Así b a+ |ni son autoestados de H

a|ni y b b con autovalores ( En − h̄ω ) y ( En + h̄ω ). De esta manera las acciones de b
ayba+ sobre
|ni generan nuevos estados de energía que son más bajos y más altos por una unidad de h̄ω. Entonces b ayb +
a son conocidos
como los operadores de creación y destrucción u operadores escalera.
Puesto que b ayba+ no conmutan con N, b los estados |ni no son estados ni de b a ni de ba+ usando (1.6) junto con la siguiente
identidad

[A
bB,b Cb] = A b[ B,
b C b] + [ A,
b C b] B,
b se obtiene:

[ N,
b b a+ b
a ] = [b a] = b
a, b a + [b a+ , b
a ] + [b
a, b a ]b
a
a+ , b
= [b a = −b
a ]b a (1.8)
[ N,
b ba + ] = [b
a+ b a+ ] = b
a, b a + [b a + ] + [b
a, b a+ , b
a + ]b
a
a+
=b (1.9)

De (1.8) y (1.9) se obtienen los siguientes términos:

N
bba=b b − 1) y N
a( N bba+ = b
a+ ( N
b + 1),

combinando estas relaciones con (1.7), se obtiene:

N a|ni) = b
b (b b − 1)|ni = (n − 1)(b
a( N a|ni)
(1.10)
+ +
a |ni) = b
b (b
N a (N a+ |ni)
b + 1)|ni = (n + 1)(b

a+ actúan sobre |ni ellos decresen e incrementan n por una unidad. De (1.10) se puede inferir que:
ayb
Cuando los operadores b

a | n i = c n | n − 1i,
b (1.11)

Donde cn es una constante que se puede hallar a partir de las condiciones de normalización para los autoestados |ni.
De (1.11):

a + ) · (b
(hn|b a|ni = |cn |2 hn − 1||n − 1i
a+ b
a|ni) = hn|b
= | c n |2 , (1.12)

de igual modo y utilizando N a+ , b

b =b a y (1.7):

a + ) · (b
(hn|b a|ni) = nhn|ni = n. (1.13)

Combinando (1.12) y (1.13), se obtiene:

|cn |2 = n.
a|ni, no puede ser negativo: n > 0. La norma no puede ser negativa.
Esto implica que n el cual es igual a la norma de b
De (1.11), √
a | n i = n | n − 1i
b (1.14)
Esta ecuación muestra que aplicaciones repetidas del operador ba sobre |ni generan una secuencia de autovectores |n − 1i, |n −
2i, | n − 3i, ·, Y b
a|0i = 0.
De igual forma, √
a + | n i = n + 1| n + 1i.
b (1.15)
Esto implica que aplicaciones repetidas de b a+ sobre |ni generan una secuencia infinita de autovectores |n + 1i, |n + 2i, |n +
3i, · · · Y el espectro de energía de un oscilador armónico es descrito como:
 
1
En = n + h̄ω, n = 0, 1, 2, 3, · · ·
2

El espectro de energía de un oscilador armónico consiste de niveles de energía igualmente espaciados,

En+1 − En = h̄ω.

Esto está basado en la idea de Planck. La energía de la radiación emitida por las cargas oscilantes (desde las paredes de la cavidad)
deben estar solo en paquetes que son múltiplos enteros de h̄ω, lo cual guía al nacimiento de la mecánica cuántica.
1.1 Las funciones de onda 4
;A<

1.1 Las funciones de onda

E
l método de operadores puede también ser empleado para determinar las funciones de onda del oscilador armónico.
Empleando (1.15), se observan diferentes vectores escritos en términos del estado fundamental |0i.

a + |0i,
|1i = b
1 + 1
|2i = √ b a |1i = √ (b a + )2 |0i
2 2!
1 + 1 + 2 1 2 1 3
|3i = √ b a |2i = √ (b a ) √ (b a + ) |0i = √ (b a + ) |0i
3 3 2! 3!
..
.
1 + 1 n
|ni = √ b a | n − 1i = √ (b a + ) |0i.
n n!

Para encontrar cualquier estado excitado |ni simplemente se opera b a+ sobre |0i n veces. Cualquier conjunto de Kets |ni y |n0 i,
0
corresponden a diferentes autovalores, y deben ser ortogonales hn |ni = δn0 ,n . Puesto que H b es hermítico y ninguno de sus
estados es degenerado.
Los estados |0i, |1i, · · · , |ni son estados simultáneos de H b El conjunto {|ni} constituye una base completa ortonormal:
b y N.

+∞
hn0 |ni = δn0 ,n ; |nihn0 | = 1.
X

n =0

Ahora se determinarán las funciones de onda del oscilador armónico en la representación posición.

1.1.1. Función de onda del estado fundamental en la representación posición

El operador
Pb −ih̄ d d
pb = √ = √ = −ix0 ,
mh̄ω mh̄ω dx dx
h̄
con x0 = , constante con dimensiones de longitud. Los operadores de creación y destrucción se pueden reescribir como:
mω
1
a = √ (qb + i pb)
b
2
 r 
1 mω d
a= √
b X
b + x0
2 h̄ dx
   
1 x d 1 2 d
a= √
b + x0 = √ X + x0
b (1.16)
2 x0 dx 2x0 dx
 
1 d
a+ = √
b b − x0 2
X
2x0 dx

a|0i = |0i en el espacio de posición como:

Empleando (1.16) se puede escribir la ecuación b
   
a |0i = √
h x |b
1
hx| Xb + x0 2 d |0i = √ 1 d
x h x |0i + x0 2 h x |0i = 0
2x0 dx 2x0 dx
utilizando ψ0 ( x ) = h x |0i,
 
1 d
√ xψ0 ( x ) + x0 2 ψ0 ( x ) =0
2x0 dx

dψ0 ( x ) x
= − 2 ψ0 ( x )
dx x0
1.1 Las funciones de onda 5
;A<

2 /2x 2
ψ0 ( x ) = Ae− x 0 : función de onda del estado fundamental y A es una constante que se puede hallar con la condición de
normalización:
Z∞ Z∞
2 2 /x 2 √
1= dx |ψ0 ( x )| = A 2
dxe− x 0 = A2 πx0
−∞ −∞
donde
 mω 1/4 1
A= = p√ .
πh̄ πx0
La función de onda normalizada es:
1 2 2
ψ0 ( x ) = p√ e− x /2x0
πx0
Esta es una fuanción Gaussiana.
El primer estado excitado es:
 
h x |1i = b+
h x | a |0i = √
1
hx| X b − x0 2 d |0i
2x0 dx
 
1 d
ψ1 ( x ) = √ xψ0 ( x ) − x0 2 ψ0 ( x )
2x0 dx
 
1 d
ψ1 ( x ) = √ x − x0 2 ψ0 ( x )
2x0 dx

−x
  
1
ψ1 ( x ) = √ x − x0 2 ψ0 ( x )
2x0 x0 2
2 x
= √ ψ0 ( x )
2 x0
s
√ x 2 2 2
ψ1 ( x ) = 2 ψ0 ( x ) = √ 3
xe− x /2x0
x0 πx0
para los estados excitados segundo y tercero:
  2
1
h x |2i = √ h x | √
1 b − x0 2 d
X |0i
2! 2x0 dx
 2   
= √
1
√
1
hx| X b − x0 2 d Xb − x0 2 d |0i
2! 2x0 dx dx
 2   
1 1 2 d 2 d
= √ √ x h x | − x0 X − x0
b |0i
2! 2x0 dx dx
 2   
1 1 2 d 2 d
= √ √ x − x0 x h x | − x0 |0i
2! 2x0 dx dx
 2  2
1 1 2 d
ψ2 ( x ) = √ √ x − x0 ψ0 ( x ),
2! 2x0 dx

del mismo modo,

3 
d 3
 
1 1
h x |3i = ψ3 ( x ) = √ √ x − x0 2 ψ0 ( x )
3! 2x0 dx
Siguiendo la secuencia para ψ0 ( x ), ψ1 ( x ), ψ2 ( x ), ψ3 ( x ) se puede escribir como:
 n
1 1 2 d 2 2
ψn ( x ) = p√ n + 1/2
x − x0 e− x /2x0
π2 n! x0
n dx
que es la relación de recurrencia para la función de onda en el n + h estado excitado del oscilador armónico.
1.1 Las funciones de onda 6
;A<

1.1.2. Segundo método: resolver la ecuación de Schrödinger

Escribiendo el correspondiente hamiltoniano mecánico-cuántico de la forma:
2 2
b = −h̄ d + m ω 2 x2
H
2m dx2 2
Puesto que el potencial es constante en el tiempo, la ecuación de Schrödinger independiente del tiempo determina la solución
estacionaria ψn y las correspondientes energías o autovalores En . La ecuación de Schrödinger independiente del tiempo toma
la forma:
−h̄2 d2 m
ψ( x ) + ω 2 x2 ψ( x ) = Eψ( x ). (1.17)
2m dx2 2
Debido a la importancia del oscilador armónico y sus soluciones para la mecánica cuántica, ahora se considera el método para
resolver (1.17) en detalle. Empleando las siguientes sustituciones:

2mE mω
K2 = 2
,λ= ,
h̄ h̄

se puede reescribir la ecuación (1.17) como:

d2 ψ
+ (K2 − λ2 x2 )ψ = 0. (1.18)
dx2
La ecuación (1.18) se conoce como la ecuación diferencial de Weber. Se introduce la transformación y = λx2 en (1.18):

dψ dψ dy dy d2 y
= , = 2λx, = 2λ,
dx dy dx dx dx2

d2 ψ dψ d2 y
 
d dψ dy d dψ dy
= = + ,
dx2 dx dy dx dx dy dx dy dx2

d2 ψ d dy dψ dy dψ
= 2λ
dx2 dy dx dy dx dy

d2 ψ d2 ψ dψ
2
= 2 4λ2 x2 + 2λ ,
dx dy dy
y
pero x2 =
λ
d2 ψ d2 ψ dψ
= 4y + 2λ ,
dx2 dy2 dy

reemplazando en (1.18),
d2 ψ dψ
4yλ 2
+ 2λ + ( K 2 − λ2 x 2 ) ψ = 0
dy dy

Dividiendo esta última expresión por 4λ:

d2 ψ 1 dψ K2 λx2
 
y + + − ψ=0
dy2 2 dy 4λ 4

K2 2Emh̄ E
haciendo b = , b= 2
=
2λ 2h̄ mω h̄ω
d2 ψ 1 dψ
 
b y
y 2 + + − Ψ=0 (1.19)
dy 2 dy 2 4

Se propone una solución de la forma:

ψ(y) = e−y/2 φ(y)
1.2 Funciones hipergeométricas 7
;A<

(1.19) se puede escribir en función de φ(y) como:

 
dψ 1 dφ(y) 1 dφ(y) −y/2
= − e−y/2 ψ(y) + e−y/2 = − φ(y) + e
dy 2 dy 2 dy

d2 ψ 1 −y/2 1 −y/2 dφ(y) 1 −y/2 dφ(y) −y/2 d

2
= e φ ( y ) − e − e + e φ(y)
dy2 4 2 dy 2 dy dy2

d2 ψ dφ(y) d2 φ(y) −y/2

 
1
= φ(y) − + e , (1.20)
dy2 4 dy dy2

introduciendo la ecuación (1.20) en (1.19):

dφ(y) d2 φ(y) −y/2 1
     
1 1 dφ(y) −y/2 b y −y/2
y φ(y) − + e + − φ(y) + e + − e φ(y) = 0
4 dy dy2 2 2 dy 2 4

d2 φ ( y )
   
1 dφ(y) b 1
y + − y + − φ(y) = 0 (1.21)
dy2 2 dy 2 4

1.2 Funciones hipergeométricas

En 1769 Euler construyó la ecuación diferencial hipergeométrica:

d2 φ dφ
z (1 − z ) + [ c − ( a + b + 1) z ] − abφ = 0 (1.22)
dz2 dz
contiene los tres parámetros libres a, b, c y posee una gran variedad de soluciones. Tiene tres puntos singulares 0, 1, ∞. La
ecuación diferencial hipergeométrica es un prototipo: cada ecuación diferencial de segundo orden con al menos tres puntos
singulares regulares puede ser llevada a la ecuación diferencial hipergeométrica por medio de un adecuado cambio de variables.
Las soluciones de la ecuación diferencial hipergeométrica incluyen muchos de las más interesantes funciones especiales de la
física matemática. Las soluciones a la ecuación diferencial hipergeométrica son construidas por las series hipergeométricas. La
solución de la ecuación diferencial hipergeométrica de Euler se llama función hipergeométrica o función Gaussiana 2 F1 introducida
por Gauss.
Para resolver (1.22) se sustituye la serie de potencias
∞
X
φ(z) = z σ
ck zk
k =0

Se encuentra la relación de recurrencia:

∞ ∞ ∞
ck (k + σ )(k + σ − 1)zk−2 + [c − ( a + b + 1)z]zσ ck (k + σ )zk−1 − abzσ
X X X
z (1 − z ) z σ ck zk = 0
k =0 k =0 k =0

reordenando los términos

∞
c 0 σ ( c + σ − 1 ) z σ −1 +
X
[ck+1 (k + σ + 1)(k + c + σ) − ck (k + a + σ)(k + b + σ)] zk+σ = 0.
k =0

Para que esta expresión se haga cero, todos los coeficientes tiene que ser iguales a cero:

σ(c − 1 + σ) = 0 (1.23)

(k + a + σ)(k + b + σ)
σk+1 = c .
(k + 1 + σ)(k + c + σ) k

Una solución de (1.22) con c0 = 1 es:

∞
X ( a + σ )k (b + σ )k k
φ(z) = zσ z
(1 + σ ) k ( c + σ ) k
k =0
1.2 Funciones hipergeométricas 8
;A<

Utilizando las abreviaturas (símbolos de Pochhammer):

( a ) k = a ( a + 1) · · · ( a + k − 1),
( a)0 = 1.
El radio de convergencia se calcula de la prueba para la razón de la convergencia:
 
 ck 
k = lı́m    = 1.
k → ∞ c k +1 

En (1.23) hay dos posibles valores para σ:

(i) σ = 0. La solución en este caso son las series hipergeométricas:

∞
X ( a)k (b)k zk
φ1 (z) = 2 F1 ( a, b, c; z) = . (1.24)
(c)k k!
k =0

Los índices en la notación 2 F1 provienen de la generalización de la serie hipergeométrica en la forma:

∞
X ( α1 ) k ( α2 ) k · · · ( α p ) k z k
p Fq ( α1 , · · · , α p ; β 1 , · · · , β p ; z) = .
( β 1 )k · · · ( β q )k k!
k =0

La solución (1.24) sólo tiene sentido si, en la serie 2 F1 , ninguno de los denominadores de los varios términos en la serie des-
aparecen. La existencia de 2 F1 implica la condición c 6= n, donde n = 0, 1, · · · . Entonces la serie es holomórfica en el círculo
unitario. Cuando a = −n y b = −n, la serie termina y define un polinomio de n + h grado. Por ejemplo:

2 F1 (− n, n + 1; 1; x ) = Pn (1 − 2x )
es un polinomio de Legendre.
Casos especiales adicionales son los polinomios de Gegenbauer y Tschebycheff, entre otros.

(ii) σ = 1 − c. En (1.24), la solución para este caso es:

φ2 (z) = z1−c 2 F1 ( a + 1 − c, b + 1 − c; 2 − c; z)

Note el factor z1−c , la solución φ2 solo existe si c 6= 2, 3, · · · . La solución general de la ecuación diferencial hipergeométrica es
por tanto:
φ(z) = A 2 F1 ( a, b; c; z) + Bz1−c 2 F1 ( a + 1; c, b + 1 − c; 2 − c; z). (1.25)
Se puede deducir otra importante ecuación diferencial de (1.22). Se hace la sustitución de la transformación lineal
dx x
x = bz, = b, z=
dz b
dφ dφ dx dφ
= =b
dz dx dz dx
2
d φ

d dφ dx

d dx dφ dx dφ d2 x
= = +
dz2 dz dx dz dx dz dx dz dx dz2

d2 φ 2
2d φ
= b .
dz2 dx2
Se hace la sustitución y (1.22) se escribe:

x x  2 d2 φ h x i dφ
1− b 2
+ c − ( a + b + 1) b − abφ = 0
b b dx b dx
 x  d2 φ h x i dφ
x 1− 2
+ c − ( a + b + 1) − aφ = 0
b dx b dx
 x  d2 φ h x i dφ
x 1− + c − ( a + 1 ) − x − aφ = 0 (1.26)
b dx2 b dx
1.2 Funciones hipergeométricas 9
;A<

tomando el límite cuando b → ∞, se obtiene la ecuación diferencial de Kummer:

d2 φ dφ
+ [c − x ] − aφ = 0 (1.27)
dx2 dx
La ecuación tiene una singularidad no esencial en x = 0 y una esencial cuando x → ∞, la cual resurge a través de la confluencia
de z = 1 y z → ∞.
La solución general de (1.27) se obtiene expandiendo una serie de potencias alrededor de x = 0,
φ( x ) = A 1 F1 ( a; c, x ) + Bx1−c 1 F1 ( a − c + 1; 2 − c; x ),
con la función hipergeométrica confluente:
∞
X ( a)k ( x )k a x a ( a + 1) x 2
1 F1 ( a; c; x ) = = 1+ + +··· (1.28)
(c)k k! c 1! c(c + 1)2!
k =0

La solución de (1.28) se origina de (1.25) tomando el límite cuando b → ∞, con x = bz.

Para a = −n, polinomios de n + h grado aparecen, en particular los polinomios de Laguerre:

(m) (n + m)!
Ln (z) = F (−n; m + 1; z)
n!m! 1 1
y los polinomios de Hermite:
 
(2n)! 1 2
H2n (z) = (−1)n F − n; ; z
n! 1 1 2
y
 
(2n + 1)! 3
H2n+1 (z) = (−1)n 2z 1 F1 −n; ; z2 .
n! 2
Después de bosquejar un suplemento matemático para enunciar las soluciones de las funciones hipergeométricas, se pasa a
analizar la solución del oscilador armónico.
Comparando (1.21) con (1.27) se identifica a (1.21) como una ecuación diferencial de Kummer. La solución general es:
   
1 1 3
φ(y) = A 1 F1 a; ; y + By1/2 1 F1 a + ; ; y (1.29)
2 2 2

donde  
b 1
a=− − .
2 4
La solución del problema físico se determina por la función de onda

ψ(y) = e−y/2 φ(y)

La función ψ es de cuadrado integrable, i.e., ψ desaparece en el infinito.
De (1.28) el comportamiento límite para | x | −→ ∞ es:
Γ(c) −iaπ − a Γ(c) x a−c
1 F1 ( a; c; x ) −→ e x + e x (1.30)
Γ(c − a) Γ( a)
Reemplazando (1.30) en (1.29) se obtiene el siguiente resultado:
Γ(1/2) −iaπ −a Γ(1/2) y a−1/2 Γ(3/2) −i(a+1/2)π − a−1/2 1/2 Γ(3/2) y a−1
φ(y) = A e y +A e y +B e y y + By1/2 e y
Γ(1/2 − a) Γ( a) Γ (1 − a ) Γ( a + 1/2)

φ(y) = C1 e−iaπ y− a + C2 ey y a−1/2 + C3 e−i(a+1/2)π y− a + C4 ey y a−1/2

en y −→ ∞,
φ(y) −→ (C2 + C4 )ey y a−1/2 = constante ey y a−1/2

ψ(y) = constante ey/2 y a−1/2 .

Para la cuantización de la energía, se consideran los dos posibles casos:
1.2 Funciones hipergeométricas 10
;A<

(i) a = −n yB = 0 con

 n =0, 1, 2, ·
· · ,.
b 1 b 1
−n = − − −→ − = n, con el autovector
2 4 2 4

2 1 2
ψ( x ) = Nn e−(λ/2) x 1 F1 (− n, , x λ) (1.31)
2
E
donde b =
h̄ω
b 1 1 E 1
= n + , b = 2n + , = 2n + ,
2 4 2 h̄ω 2
y la energía:
1
En = (2n + )h̄ω (1.32)
2
1
(ii) Haciendo a + = −n, A = 0 en (1.29), se obtiene:
2

1
a = −n −
2
 
1 b 1
−n − =− − ,
2 2 4

resolviendo para b,
1 b 1
n+ = − ,
2 2 4
b 1 1
= n+ + ,
2 2 4
1
b = 2n + 1 + ,
2
3 En 3
b = 2n + −→ = 2n +
2 h̄ω 2
 
3
En = 2n + h̄ω (1.33)
2

El autovector correspondiente es:

 
2 /2 3
ψn ( x ) = Nn e−λx x 1 F1 −n, , λx2 (1.34)
2

Combinando las ecuaciones (1.32) y (1.33), se obtiene el espectro de energía discreto total i.e., los autovalores permitidos en
(1.32) y (1.33). Se pueden concebir conjuntamente de la forma:
 
1
En = n + h̄ω, con n = 0, 1, 2, · · · ,
2

y este es el espectro de energía del oscilador armónico.

Este espectro de energía es equidistante con el espaciamiento h̄ω y tiene un valor finito en el estado fundamental (n = 0). La
h̄ω
energía de punto cero es .
2
Los polinomios que ocurren en (1.31) y (1.34) se conocen como los polinomios de Hermite y se definen así:
 
(2n)! 1 2
H2n ( E) = (−1)n F − n, , E
n! 1 1 2
 
n 2(2n + 1) ! 3 2
H2n+1 ( E) = (−1) E 1 F1 −n, , E .
n! 2
1.2 Funciones hipergeométricas 11
;A<

Los autovectores y energías pueden ser escritas como:

−λx2 /2
√ 
1

ψn = Nn e H2n ( λx ), En = 2n + h̄ω (1.35)
2

−λx2 /2
√ 
1

ψn = Nn e H2n+1 ( λx ), En = (2n + 1) + h̄ω (1.36)
2
con
n = 0, 1, 2, 3, · · · ,
(1.35) y (1.36) se pueden escribir conjuntamente como:

−λx2 /2
√ 
1

ψn = Nn e Hn ( λx ), En = n+ h̄ω (1.37)
2
con
n = 0, 1, 2, 3, · · · ,
Para los polinomios de Hermite, se tiene la relación
2 d n − E2
Hn ( E) = (−1)n e E e (1.38)
dEn
La ecuación diferencial de Hermite, es dada por:
d2 y dy
2
− 2E + 2ny = 0 (1.39)
dE dE
Y para diferentes aplicaciones se requiere encontrar soluciones las cuales son finitas para todos los valores finitos de E y son
de tal forma que e1/2 , y( E) −→ 0 cuando E −→ ∞.
La solución general para (1.39) se denota por Hn ( E) y se llama el polinomio de Hermite de orden n:
[1/2n]
n!
(−1)r (2E)n−2r
X
Hn ( E) = (1.40)
r!(n − 2r )!
r =0

2
∞ tn
Teorema 1.1 (Función generatriz). e2tE−t =
P
Hn ( E)
n=0 n!

La demostración es la siguiente:

Demostración.
2 2
e2tE−t = e2tE e−t
∞ ∞ s
X (2tE)r X (−t2 )
=
r! s!
r =0 s =0

(2E)r r+2s
(−1)s
X
= t (1.41)
r!s!
r,s=0

Para un valor fijo de s se obtiene tn , tomando r + 2s = n, i.e., r = n − 2s, para un valor fijo de s el coeficiente de tn en (1.41)
es:
(2E)n−2s
(−1)s (1.42)
(n − 2s)!s!
Los coeficientes totales de tn se obtienen al sumar sobre todos los valores permitidos de s y n − 2s > 0.
1 1 1
s 6 n. Así si n es par, s va de 0 a n, mientras que si n es impar, s va de 0 a (n − 1).
2 2     2
1 1
Es decir, en todos los casos, s va de 0 a n , con n definido por:
2 2
1

  
1  n si n es par
n = 2
2  1 (n − 1) si n es impar

2
1.2 Funciones hipergeométricas 12
;A<

(1.41) queda como:

[n/2] [1/2n]
X (−1)s (2E)n−2s (−1)s (2E)n−2s tn n!
2tE−t2
X
e = tn =
(n − 2s)!s! (n − 2s)!s!n!
s =0 s=0,n=0

∞ [X
n/2]
2tE−t2
X (−1)s (2E)n−2s tn
e = n!
(n − 2s)!s! n!
n =0 s =0

∞
X tn
= Hn ( E) (1.43)
n!
n =0

Ahora, se demostrará la ecuación (1.38).

El teorema de Taylor es:

dF (0) t d2 F (0) t
F ( t ) = F (0) + + + · · · +, (1.44)
dt 1! dt2 2!
∞  n 
X d F tn
F (t) = n
(1.45)
dt t=0 n!
n =0

Se tiene entonces:
∞
X dn Hn tn
Hn ( E) =
∂tn n!
n =0

∂n 2tE−t2
 
= e
∂tn t =0
 n 
∂ E2 −(E−t)2
= e
∂tn t =0
 n 
E2 ∂ −( E−t)2
=e e
∂tn t =0

∂ ∂
F( x − t ) = − F( x − t )
∂t ∂x
a si que:
∂n n ∂
n
F ( x − t ) = (− 1 ) F( x − t )
∂tn ∂x n
y se tiene que:
∂n −(E−t)2
 
2
Hn ( E) = (−1)n e E e
∂En t =0

2 ∂ n − E2
= (−1)n e E e
∂En
que es (1.38)
1.2 Funciones hipergeométricas 13
;A<

Los Hn ( E) son polinomios de grado n en E con el término dominante 2n En .

Los primeros cinco polinomios, se pueden derivar de (1.38):

H0 ( E) = 1
H1 ( E) = 2E
H2 ( E) = 4E2 − 1 (1.46)
3
H3 ( E) = 8E − 12E
H4 ( E) = 16E4 − 48E2 + 12,

cada uno de estos polinomios es una solución de la ecuación (1.39). La solución de (1.37) se combinó con (1.46) y se introdujo
√
E = λx
2 /2
ψn ( x ) = Nn e− E Hn ( E) (1.47)

La constante Nn , la cual depende de los n índices, se determina por la condición de normalización:

Z∞
|ψn ( x )|2 dx = 1, (1.48)
−∞

el 1 en la ecuación (1.48) da la probabilidad total para encontrar la partícula en todo el espacio.

Z∞ Z∞
2 1 2
|ψn ( x )| dx = √ Nn2 e− E Hn ( E)2 dE = 1
λ
−∞ −∞

se emplea (1.38)
Z∞ Z∞
2 1 d n − E2
|ψn ( x )| dx = √ Nn 2 (−1)n Hn ( E) e dE
λ dEn
−∞ −∞

Z∞ Z∞
dn 2
h 2
i∞ dHn dn−1 − E2
Hn ( E) n e− E dE = e− E Hn ( E) − e dE
dE −∞ dEn dEn−1
−∞ −∞

El primer término desaparece.

Y después se integra n veces el segundo término:

Z∞ Z∞
dn 2 dn 2
Hn ( E) n e− E dE = (−1)n Hn e− E dE,
dE dEn
−∞ −∞

por una propiedad de los polinomios de Hermite:

dn
Hn ( E) = 2n n!
dEn
Z∞ Z∞
dn 2 2 √
Hn ( E) n e−E dE = (−1)n 2n n! e− E dE = (−1)n 2n n! π
dE
−∞ −∞

y la constante de normalización es:

Nn 2 (−1)n √
√ (−1)n 2n n! π = 1
λ
s√
λ 1
Nn = √
π 2n n!
1.2 Funciones hipergeométricas 14
;A<

Los estados estacionarios del oscilador armónico en mecánica cuántica son:

s
√
r
1 λ −1/2λx2
ψn ( x ) = n
e Hn ( λx ) (1.49)
2 n! π
para n = 0 y de la ecuación (1.46)
r
4 λ −1/2λx2
ψ0 ( x ) = e
π
para n = 1:
s r
1 λ −1/2λx2 √
ψ1 ( x ) = 2 e λx
2 π
para n = 2:
s r
1 λ −1/2λx2
ψ2 ( x ) = e (4λx2 − 2)
8 π
x0 en (1.16) es:
r
h̄ mω
x0 = y λ=
mω h̄
h̄ 1 1
x02 = = λ= 2
mω λ x0
v s !
u
u1 1 −1/2λx2 /x2 x2
ψ2 ( x ) = t e 0 4 2 −2
8 x02 π x0

(1.49) es la relación de recurrencia que se encontró con el método de operadores.

Las autofunciones de (1.49) tienen la siguiente propiedad de simetría:

ψn (− x ) = (−1)n ψn ( x )

Lo que significa que si

n es par:
ψn (− x ) = ψn ( x ) −→ paridad +1
n es impar:
ψn (− x ) = −ψn ( x ) −→ paridad −1.

Clásicamente, una partícula puede estar en el estado de equilibrio en x = 0, P = 0, E = 0. En Mecánica Quántica, la energía
h̄ω
posible más pequeña tiene el valor E = : la energía de punto cero.
2
Esta energía de punto cero es una consecuencia directa de la relación de incertidumbre:

h̄
∆X 2 ∆P2 = (1.50)
4
Se analizará con más detalle la (1.50).
2 2 2
∆P2 = ( P − P) = P2 − 2PP + P = P2 − 2PP + P
2
= P2 − P .
Análogamente
2
∆X 2 = X 2 − X .
1.3 Pozo de potencial infinito 15
;A<

En un estado con un valor de energía fijo, los valores medios P y X son iguales a cero, puesto que el integrando es una función
impar:
Z∞ Z∞
X= ψn∗ ( x ) xψn ( x )dx = |ψn ( x )|2 xdx = 0,
−∞ −∞

y
Z∞ Z∞
d
Px = ψn∗ ( x ) Pbx ψn ( x )dx = −ih̄ Ψ∗n ( x ) ψn ( x )dx
dx
−∞ −∞

Z∞ ∞
−ih̄ d −ih̄
|ψn ( x )|2 dx = |ψn ( x )|2 

= =0
2 dx 2 −∞
−∞

por tanto
∆P2 = P2 , ∆X 2 = X 2 (1.51)

Con (1.51) se puede escribir la relación de incertibumbra en la forma:

h̄
P2 X 2 > (1.52)
4
La energía promedio del oscilador es:
P2 mω 2 2
H= + X . (1.53)
2m 2
Comparando (1.52) y (1.53), se nota que un incremento en la energía potencial conduce a una disminución en la energía cinética
y viceversa.
Combinando ambas ecuaciones, se obtiene:
P2 mω 2 h̄2
H= +
2m 8P2
1
La función H = H ( P2 ) posee un mínimo en P2 = mωh̄. Puesto que un estado de energía fija se caracteriza por H = E, se
2
obtiene como valor mínimo de los posibles autovalores de energía

h̄ω
mı́n E > = E0 .
2
La energía de punto cero E0 por tanto es el valor de energía más pequeño que es compatible con la relación de incertidumbre.
Es posible indicar el movimiento de punto cero (lo cual guía a la energía de punto cero) por observar la dispersión de la luz
en cristales. En cuerpos sólidos, los átomos y moléculas ejecutan pequeñas oscilaciones; de acuerdo a la teoría clásica, ellos
deberían disminuir con la temperatura decreciente. Estas oscilaciones atómicas son la causa de la dispersión de la luz, lo cual,
por tanto, desapareciera con la disminución de la temperatura. Sin embargo, los experimentos demuestran que la intensidad de
la luz dispersada converge hacia un límite finito, demostrando que aún en el punto del cero absoluto, las oscilaciones atómicas
ocurren.

1.3 Pozo de potencial infinito

Siguiendo con las aplicaciones de la ecuación de Schrödinger, se encuentra el pozo infinito de potencial. Se resolverá el pozo
infinito de potencial unidimensional.
Se considera una partícula de maasa m confinada a moverse dentro de un pozo infinito de potencial, (véase Fig. 2.1).

+∞ si x < 0

V ( x ) = 0 si 0 6 x 6 L (1.54)
+∞ si x > L


1.3 Pozo de potencial infinito 16
;A<

O L
Figura 1.1. Pozo infinito de potencial

L es el ancho del pozo.

Clásicamente la partícula permanece confinada dentro del pozo con energía total:

1 2 P2
E=K= mv =
2 2m
√
P = ± 2mE
con movimiento hacia adelante y hacia atrás como resultado de constantes reflexiones contra las paredes del pozo.
Cuánticamente, la solución está condicionada por la ecuación de Schrödinger independiente del tiempo:

−h̄2 d2
ψ( x ) + V ( x )ψ( x ) = Eψ( x )
2m dx2
De acuerdo a (1.54), V ( x ) = 0 entre x = 0 y x = L

−h̄2 d2 d2 Ψ ( x ) −2mE
ψ ( x ) = Eψ ( x ) −→ = Ψ( x ) (1.55)
2m dx2 dx 2
h̄2
2mE
Se hace K2 = en (1.55):
h̄2
d2 ψ ( x )
= −h̄2 ψ( x ),
dx2
las raíces de la ecuación diferencial son: (
2 2 λ1 = iK
λ = −K
λ2 = −iK
La solución general de (1.55) es:
ψ( x ) = A1 eihx + B1 e−ihx (1.56)
Las funciones de onda desaparecen en las paredes del pozo, de modo que en x = 0, la solución (1.56) es:
ψ(0) = A1 + B1 = 0 −→ B1 = − A1
(1.56) se transforma en:
h i
ψ( x ) = A1 eiKx − ie−iKx

ψ( x ) = A1 [  + i sin Kx − 
Kx
cos  + i sin Kx ]
hx
cos

ψ( x ) = 2A1 i sin Kx,

se hace A = 2A1 i
ψ( x ) = A sin Kx (1.57)
En la pared del pozo también ψ = 0 si x = L. De (1.57)
0 = A sin KL −→ KL = nπ
nπ
Kn = , n = 0, 1, 2, · · · ,
2
1.4 Escalón de potencial 17
;A<

nπ
(1.57) queda como ψ( x ) = A sin x
L
y
2mE n2 π 2 h̄2 π 2 2
K2 = = −→ En = n (1.58)
h̄2 L2 2mL2

Estos son los valores permitidos que puede tomar la partícula dentro del pozo, i.e., la energía es cuantizada.
Para hallar A en (1.57), se normaliza la función de onda:

ZL
|ψ( x )|2 dx = 1
0

ZL
nπx
A2 sin2 dx = 1
L
0

ZL ZL  
2 nπx
2 1 1 2nπx
A sin dx = 1 −→ A2 − cos dx = 1
L 2 2 L
0 0
"  #
1  L 2nπx  L

L
A2 x − sin =1
2 0 4nπ L 0

1  L
 r
2 1 2
A x = 1, LA2 = 1 −→ A = .
2 0 2 L

La función de onda para la partícula de masa m en el pozo de potencial infinito se escribe de la siguiente forma:
r
2 nπx
ψ( x ) = sin
L L
para n = 1, 2, 3, 4, · · · ,

1.4 Escalón de potencial

Otro problema simple, donde se aplica la ecuación de Schrödinger, consiste de una partícula que es libre de un potencial para
( pero para x > 0 el potencial se incrementa. Un potencial de esta forma se llama escalón de potencial, (véase Fig.).
cierta región;
0, x < 0,
V (x) =
V0 , x > 0.
En este problema se intenta analizar la dinámica de un flujo de partículas, con la misma masa m y la misma velocidad v, mo-
viéndose de izquierda a derecha, (véase Fig.)
Se consideran dos casos:

(i) Caso E > V0 .

Las partículas son libres para la región x < 0, y sienten un potencial repulsivo V0 que comienza en x = 0 y se mantiene
constante para x > 0.
Se analizan los casos clásicos y cuánticos de este ejercicio.
√
Clásicamente: las partículas se aproximan al escalón de potencial desdep la izquierda con momentum constante P = 2mE.
Cuando las partículas entran en la región x > 0, el momentum es P = 2m( E − V0 ), este es el momentum que la partícula
experimenta cuando se mueven a la derecha. Como E > V0 , las partículas tendrán la energía sufuciente para entrar en la región
x > 0.
Estas partículas tendrian transmisión total.
Esto es un problema simple de dispersión en una dimensión.
Cuánticamente: la dinámica de estas partículas es regulada por la ecuación de Schrödinger, la cual para este problema, se aplica
1.4 Escalón de potencial 18
;A<

de la siguiente forma:
para x < 0, V = 0
h̄2 d2 ψ I ( x )
− = Eψ I ( x ),
2m dx2
d2 ψ I ( x ) −2mE 2mE
= ψ I ( x ) −→ K12 = 2
dx2 h̄ 2
h̄
d2 ψ I ( x )
+ K12 ψ I ( x ) = 0 para x < 0 (1.59)
dx2
La ecuación de Schrödinger para la región x > 0 es:
−h̄2 d2 ψ I I ( x )
+ V0 ψ I I ( x ) = Eψ I I ( x ) (1.60)
2m dx2
d2 ψ I I 2m
= −( E − V0 ) 2 ψ I I ( x ),
dx2 h̄
2m
se hace K22 = ( E − V0 )
h̄2

d2 ψ I I ( x )
+ K22 ψ I I ( x ) = 0 para x > 0 (1.61)
dx2
Las soluciones de las ecuaciones (1.59) y (1.61) son soluciones de onda plana.
ψ I ( x ) = AeiK1 x + Be−iK1 x ,

x < 0,
(1.62)
ψ I I ( x ) = CeiK2 x + De−iK2 x , x>0


en donde AeiK1 x y Ce−iK2 x representan ondas que se mueven en la dirección x positiva, y Be−iK1 x y De−iK2 x corresponden a
ondas que se mueven en la dirección x negativa, (véase Fig.)
Las ondas pueden ser transmitidas o reflejadas en x = 0, puesto que ninguna onda se refleja despues de entrar en la región
x > 0, entonces D = 0.
La función de onda completa Ψ( x, t) para estas soluciones es:

ψ I ( x )e−iωt = Aei(K1 x−ωt) + Be−i(K1 x+ωt) , x < 0
Ψ( x, t) = (1.63)
ψ ( x )e−iωt = Cei(K2 x−ωt) , x > 0
II

Donde Aei(K1 x−ωt) , Be−i(K1 x+ωt) , Cei(K2 x−ωt) representan las ondas incidente, reflejada y transmitida.
La densidad de probabilidad es una constante para x > 0 ya que:
2
|ψ I I ( x )|2 = |Cei(K2 x−ωt) | = CC ∗ ei(K2 x−ωt) e−i(K2 x−ωt)
= | C |2 .
Ahora se evaluaran los coeficientes de reflexión y transmisión, R y T. Es decir, las probabilidades de que las partículas se reflejen
o se transmitan.
Estos se definen:
   
 Densidad de corriente reflejada   JR 
R=   =  ,
J 
Densidad de corriente incidente  I
   
J  J 
T =  R  , R =  R  (1.64)
JI JI
La cantidad de corriente se obtiene de la siguiente manera. Se parte de la ecuación de Schrödinger dependiente del tiempo en
una dimensión:
∂Ψ −h̄2 ∂2 Ψ
ih̄ = + VΨ, (1.65)
∂t 2m ∂x2
1.4 Escalón de potencial 19
;A<

donde Ψ = Ψ( x, t)
Se conjuga (1.65),
∂Ψ∗ −h̄2 ∂2 ∗
−ih̄ = Ψ + VΨ∗ (1.66)
∂t 2m ∂x2
Se multiplica (1.65) por Ψ∗ :
∂Ψ ∗ −h̄2 ∂2 Ψ∗
ih̄ Ψ = Ψ + VΨΨ∗ , (1.67)
∂t 2m ∂x2
se multiplica (1.66) por Ψ:
∂Ψ∗ −h̄2 ∂2 ∗
−ih̄ Ψ= Ψ Ψ + VΨΨ∗ (1.68)
∂t 2m ∂x2
se resta (1.68) de (1.67):

∂Ψ ∗ ∂Ψ∗ −h̄2 ∂2 Ψ ∗ h̄2 ∂2 Ψ∗

ih̄ Ψ + ih̄ Ψ= Ψ + Ψ
∂t ∂t 2m ∂x2 2m ∂x2
−h̄2 ∂2 Ψ ∂2 Ψ ∗
 
∂
ih̄ [Ψ∗ Ψ] = Ψ∗ 2 − Ψ
∂t 2m ∂x ∂x2
∂Ψ∗
 
∂ ih̄ ∂ ∂Ψ
[Ψ∗ Ψ] = Ψ∗ −Ψ
∂t 2m ∂x ∂x ∂x

∂Ψ∗
 
∂ ∗ ih̄ ∂ ∗ ∂Ψ
[Ψ Ψ] = − Ψ −Ψ
∂t 2m ∂x ∂x ∂x

∂Ψ∗
  
∂ ∂ ih̄ ∂Ψ
[Ψ∗ Ψ] = − Ψ − Ψ∗ (1.69)
∂t ∂x 2m ∂x ∂x

∂
haciendo ∇ x = y ρ = Ψ∗ Ψ (densidad de probabilidad) en (1.69):
∂x
∂ρ
= ∇x J
∂t
∂ρ → →
+ ∇x · J = 0 (1.70)
∂t
esta ecuación representa la conservación del flujo de probabilidad. Se denomina ecuación de continuidad, y representa la den-
sidad de corriente de probabilidad en (1.69). Se comprueba que:

∂Ψ∗
 
ih̄ ∗ ∂Ψ
J= Ψ −Ψ . (1.71)
2m ∂x ∂x

Se calcula la densidad de corriente incidente.

De (1.62):

Ψ Ii = AeiK1 x onda incidente en el escalón

ih̄  iK1 x 
Jincidente = Ae A(−iK1 )e−iK1 x − Ae−iK1 x A(iK1 )eiK1 x
2m
ih̄  
Jincidente = − A2 iK1 − iK1 A2
2m
h̄
Jincidente = K A2 . (1.72)
m 1
1.4 Escalón de potencial 20
;A<

Luego se calcula la densidad de corriente reflejada:

Ψ Ir = Be−iK1 x onda reflejada

ih̄  −iK1 x 
Jreflejada = Be B(iK1 )eiK1 x − BeiK1 x B(−iK1 )e−iK1 x
2m
ih̄
Jreflejada = 2B2 iK1
2m
h̄
Jreflejada = − K B2 . (1.73)
m 1
Finalmente, se calcula la densidad de corriente transmitida. Por (1.62):

Ψ I It = CeiK2 x

ih̄  iK2 x 
Jtransmitida = Ce C (−iK2 )e−iK2 x − Ce−iK2 x C (iK2 )eiK2 x
2m
ih̄  
Jtransmitida = −C2 iK2 − iK2 C2
2m
h̄
Jtransmitida = K2 C 2 (1.74)
m
Se reemplaza a la ecuación (1.64). Y se calculan los coeficientes de reflexión y transmisión.

 −h̄K1 B2 /m  2
 
R=  = B ,
h̄K1 A2 /m  A2
 h̄K2 C2 /m 
 
T =  ,
h̄K1 A2 /m 

K2 C 2
T= . (1.75)
K1 A 2

Para determinar las constantes B y C, se utilizan las condiciones de frontera de las funciones de onda en x = 0. La función de
onda y su primera derivada son continuas en x = 0.
De (1.62):

en x = 0
Ψ I = Ψ I I −→ A + B = C. (1.76)
 
dΨ I  dΨ I I  iK1 x

−iK1 x 

iK2 x 
= −→ AiK 1 e − BiK 1 e = CiK 2 e ,
dx x=0 dx x=0
 
x =0 x =0

AiK1 − BiK1 = CiK2 ,

K1 ( A − B) = K2 C. (1.77)

Combinando (1.76) y (1.77) se obtiene:

K1 K K
A+B = ( A − B) −→ B + 1 B = 1 A − A
K2 K2 K2
K1 − K2
   
K2 + K1
B =A
K2 K2
K1 − K2
B= A (1.78)
K2 + K1
1.4 Escalón de potencial 21
;A<

del mismo modo se tiene:

A − C = −B

AK1 − BK1 = K2 C

AK1 − K2 C = BK1
K2
A− C = B =C−A
K1
 
K2 2K1
2A = C 1 + −→ C = A (1.79)
K1 K1 + K2

La constante A puede también ser calculada de la condición de normalización para la función de onda, pero no es necesaria
aquí.
De (1.75):

B2 B K − K2
R= , de (1.78) = 1
A2 A K1 + K2

( K1 − K2 ) 2 K2 4K12
R= 2
, de (1.75) y (1.79) T=
( K1 + K2 ) K1 ( K1 + K2 ) 2

4K1 K2
T= ,
( K1 + K2 ) 2
haciendo la siguiente sustitución:
s r
K2 2m( E − V0 ) V0
F= = 2 2
= 1−
K1 h̄ 2m/h̄ E

Así,
(1 − F )2
R= ,
(1 + F )2
4F
T= .
(1 + F )2

La probabilidad total es R + T = 1 como debe ser.

En contraste con la mecánica clásica, en donde ninguna de las partículas se reflejaría, en mecánica cuántica el coeficiente de
reflexión es diferente de cero a pesar de que E > V0 . Este efecto se debe al comportamiento ondulatorio de las partículas.

(ii) Caso E < V0 .

√
Desde el punto de vista clásico, las partículas que llegan al escalón de potencial desde la izquierda (con momentum P = 2mE)
se detendrán en x = 0 y regresarán con las magnitudes de sus momentos sin cambiar. Ninguna de las partículas transpasarán
la barrera en x = 0. Existe reflexión total de las partículas.
Desde la cuántica, la ecuación de Schrödinger y la función de onda en la región x < 0 son dadas por la ecuación (1.60) y el
primer término de la (1.62) respectivamente. Para x > 0, la ecuación de Schrödinger es:
 2 
d 2
− K 3 Ψ2 ( x ) = 0, ( x > 0) (1.80)
dx2
Donde,
2m(V0 − E)
K32 =
h̄2
La solución de (1.80) es:
Ψ2 ( x ) = Ce−K3 x + DeK3 x , ( x > 0).
1.5 La barrera de potencial 22
;A<

Ya que la función de onda debe ser finita en todas partes y el término eK3 x diverge cuando x → ∞, la constante D → 0. La
función de onda completa es:
(
Aei(K1 x−ωt) + Be−i(K1 x+ωt) , x < 0, (1.81)
Ψ( x, t) =
Ce−K3 x eiωt , x > 0 (1.82)

Como se hizo en el caso (i ), se evalúan los coeficientes de reflexión y transmisión.

El coeficiente de transmisión T, el cual corresponde a la función de onda ψT = Ce−K3 x , es cero ya que Ψ T es una función
puramente real (ψT∗ = ψT ).
 
h̄ d d
JTransmitido = ψT ψT − ψT ψT = 0
2im dx dx
Por lo tanto, el coeficiente de reflexión R debe ser igual a 1. Este resultado se puede derivar, aplicando las condiciones de
continuidad en x = 0 de (1.81) y (1.82) las funciones de onda son:

Ψ1 ( x ) = AeiK1 x + Be−K1 x
Ψ2 ( x ) = Ce−K3 x
Ψ1 (0) = Ψ2 (0)
A+B = C
 
dΨ1  dΨ2 
=
dx x=0 dx x=0
iK1 A − BiK1 = −K3 C
iK1 A − BiK1 = −K3 A − K3 B
B(K3 − iK1 ) = A(−K3 − iK1 )
B(K1 + iK3 ) = A(−K1 − iK3 )
K − iK3
B= 1 A,
K1 + iK3
de igual forma,
2K1
C= A
K1 + iK3
anteriormente se había hallado que,
B2 | B2 |
R= =
A2 | A |2
 K1 − iK3 2
 
R=
 
K1 + iK3 

K12 + K32
R= = 1, yT =0
K12 + K32

se tiene reflexión total como en el caso clásico, pero existe una diferencia con este caso, clásicamente ninguna de las partículas
pueden ser halladas en la región x > 0, cuánticamente existe una probabilidad diferente de cero en que la función de onda
penetre la región clásicamente prohibida. Se halla la densidad de probabilidad relativa. De (1.82)

P( x ) = |Ψ( x, t)|2 = |ψ2 ( x )|2 = |C |2 e−2K3 x

4K12 | A|2 −2K3 x

P( x ) = e
K12 + K32

De la Fig. 7 se ve que la probabilidad es apreciable cerca a x = 0.

1.5 La barrera de potencial

Considere un haz de partículas de masa m que son enviadas desde la izquierda sobre una barra de potencial:
1.5 La barrera de potencial 23
;A<

0, x < 0

V ( x ) = V0 , 0 6 x 6 L

0, x > L,


Este, como en el caso anterior, representa un caso de dispersión.

(i) Caso E > V0 .

Clásicamente
√ pizquierda (véase Fig. 8) con momentum P1 =
las partículas que se aproximan a la barrera provenientes de la
2mE, cuando ellas entran a la región 0 6 x 6 L con momentum P2 = 2m( E − V0 ), mantendrán este momentum P2
hasta que alcancen el punto x = L. √
Entonces, tan pronto como pasan más allá de x > L poseen momentum P3 = 2mE y mantienen este valor para toda
la región x > L, puesto que las partículas poseen E > V0 ninguna de las partículas serán reflejadas de regreso. Todas las
partículas pasarán al lado derecho de x = L, se tendrá transmisión total.
Cuánticamente la función de onda tiene un comportamiento oscilatorio en las tres regiones. Su amplitud se reduce cada vez
que las partículas pasan a una nueva región.
La ecuación de Schrödinger para la región x < 0 es:

h̄2 d2 Ψ1 d2 ψ1 2mE
− = Eψ1 −→ = − 2 ψ1
2m dx2 dx2 h̄
haciendo,
2mE d2 ψ1
K12 = 2
−→ = −K12 ψ1
h̄ dx2

La solución de esta ecuación diferencial es:

ψ1 ( x ) = AeiK1 x + Be−iK1 x x<0 (1.83)

Para la región 0 6 x 6 L la ecuación de Schrödinger es:

h̄2 d2 ψ2
− + V0 ψ2 = Eψ2
2m dx2
d2 ψ2 2m
= − 2 ( E − V0 )ψ2 ,
dx2 h̄
se hace,
2m
K22 = ( E − V0 )
h̄2
d2 ψ2
= −K22 ψ2 (1.84)
dx2
La solución de (1.84) es:
ψ2 ( x ) = CeiK2 x + De−iK2 x , 0 6 0 6 L, (1.85)
Y finalmente para la región x > L

d2 ψ3
= −K12 ψ3 , (1.86)
dx2
la solución para (1.86) es:
ψ3 ( x ) = EiK1 x , x>L (1.87)

Las constantes A, B, C, D, E se pueden obtener en términos de A utilizando las condiciones de frontera.

dΨ( x )
ψ( x ) y deben ser continuas en x = 0 y x = L.
dx
Entonces, se aplican las condiciones de frontera en x = 0.

x=0
ψ1 (0) = ψ2 (0)
1.5 La barrera de potencial 24
;A<

Igualando (1.83) y (1.85) en x = 0,

A+B = C+D (1.88)

dψ1   
= AiK1 eiK1 x − BiK1 e−iK1 x  AiK1 − BiK1

dx 
x =0 x =0

dψ1  
iK2 x −iK2 x

= CiK 2 e − DiK 2 e = CiK2 − DiK2
dx x=0

igualando,
 
dψ1  dψ2 
=
dx x=0 dx x=0

−→ AiK1 − BiK1 = CiK2 − DiK2

iK1 ( A − B) = iK2 (C − D ) (1.89)

igualando (1.85) y (1.87) en x = L

CeiK2 L + De−iK2 L = EeiK1 L (1.90)

 
dψ2  dψ3 
= −→ CiK2 eiK2 L − DiK2 e−iK2 L = EiK1 eiK1 L
dx x= L dx x= L
 
iK2 CeiK2 L − De−K2 L = iK1 EeiK1 L . (1.91)
 
J 
Se resuelve E en términos de A. De (1.64), el coeficiente de transmisión es T =  T , la densidad de corriente se obtiene de
JI
(1.71):

∂Ψ∗
 
ih̄ ∗ dΨ
J= Ψ −Ψ
2m dx dx
JT = densidad de corriente transmitida
J I = densidad de corriente incidente

De la Fig. 8,

Ψ I = AeiK1 x (onda incidente)

Ψ T = EeiK1 x (onda transmitida a la región x > L)

ih̄  
JI = − AeiK1 x AiK1 e−iK1 x − Ae−iK1 x AiK1 eiK1 x
2m
ih̄ h̄
JI = − 2A2 iK1 = A2 K1
2m m
ih̄  
JT = − EeiK1 x EiK1 e−iK1 x − Ee−iK1 x EiK1 eiK1 x
2m
ih̄ h̄
JT = − 2E2 iK1 = E2 K1
2m m
h̄ 2
E K1 E2 | E |2
T= m = 2 = . (1.92)
h̄ 2
A K1
A | A |2
m
1.5 La barrera de potencial 25
;A<

De (1.90) y (1.91):

CeiK2 L = EeiK1 L − De−iK2 L

K2 − K1
 
i ( K1 − K2 ) L −2iK2 L i ( K1 − K2 ) L E
C = Ee − De = Ee − e i ( K1 − K2 ) L
2 K2
K1 iK1 L
EeiK1 L − De−iK2 L − De−iK2 L = Ee ,
K2
 
K
2De−iK2 L = E eiK1 L − 1 eiK1 L
K2

E K 2 − K 1 i ( K1 + K2 ) L
 
D= e .
2 K2

De (1.88) y (1.89):

B =C+D−A (1.93)

E K 2 − K 1 i ( K1 + K2 ) L
   
K2
A−C−D+A = C− e
K1 2 K2

E K 2 − K 1 i ( K1 + K2 ) L K E K 2 − K 1 i ( K1 + K2 ) L
     
K
2A − e = C 1+ 2 − 2 e
2 K2 K1 K1 2 K2

E K 2 − K 1 i ( K1 + K2 ) L E K 2 − K 1 i ( K1 + K2 ) L E K 2 − K 1 i ( K1 − K2 ) L
        
i ( K1 − K2 ) L K1 + K2
2A = e − e + Ee − e
2 K2 2 K1 2 K2 K1
" ! #
K2 − K1 iK2 L K2 − K1 iK2 L K22 − K12
    
E K2 + K1 −iK2 L −iK2 L
e − e +2 e − e eiK1 L = 2A
2 K2 K1 K1 K1 K2

Eh i eiK1 L
K1 (K2 − K1 )eiK2 L − K2 (K2 − K1 )eiK2 L + 2K2 (K1 + K2 )e−iK2 L − K22 e−iK2 L + K12 e−iK2 L = 2A
2 K1 K2

Eh i eiK1 L
(K1 K2 − K12 − K22 + K1 K2 )eiK2 L + (2K1 K2 + 2K22 − K22 + K12 )e−iK2 L = 2A
2 K1 K2
Eh i
(K1 + K2 )2 e−iK2 L − (K1 − K2 )2 eiK2 L = 2K2 K2 Ae−iK1 L
2
4K1 K2 Ae−iK1 L
E= h i
(K1 + K2 )2 e−iK2 L − (K1 − K2 )2 eiK2 L

4K1 K2 Ae−iK1 L
E= h i
(K1 + K2 )2 (cos(K2 L) − i sin(K2 L)) − (K1 − K2 )2 (cos(K2 L) + i sin(K2 L))

4K1 K2 Ae−iK1 L
E=  . (1.94)
4K1 K2 cos(K2 L) − 2i (K12 + K22 ) sin(K2 L)
1.5 La barrera de potencial 26
;A<

De (1.92), el coeficiente de transmisión es:

| E |2 16K12 K22
T= =
| A |2 |4K1 K2 cos(K2 L) − 2i (K12 + K22 ) sin(K2 L)||4K1 K2 cos(K2 L) + 2i (K12 + K22 ) sin(K2 L)|

16K12 K22
T= h 2
i
16K12 K22 cos2 (K2 L) + 4(K12 + K22 ) sin2 (K2 L)

16K12 K22
T= 
16K12 K22 cos2 (K2 L) + (4K14 + 4K24 + 8K12 K22 ) sin2 (K2 L)


1
T= " ! #
K12 K2 1
cos2 (K 2 L) + + 22 + 2
sin (K2 L)
4K22 4K1 2

1
T= " #
K12 K2 sin2 (K2 L)
2
1 − sin (K2 L) + sin (K2 L) + 22 sin2 (K2 L) +
2
4K22 4K1 2

1
T= " !#
sin (K2 L) sin2 (K2 L)
2 K12 K22
1− + +
2 4 K22 K12

1
T= " ! #
1 K14
+ K24 2
1+ −2 + sin (K2 L)
4 K12 K22

1
T=  !2  (1.95)
K12 − K22
1 + 1 sin2 (K2 L)
4 K1 K2

!2
K12 − K22 V02
= . (1.96)
K1 K2 E( E − V0 )

Reemplazando (1.96) en (1.95):

1
T= " # (1.97)
V02 2√
 q
2
1+ sin L 2mV0 /h̄ E/V0 − 1
4E( E − V0 )
q
2mV0 E
usando la sustitución λ = L ye = ,
h̄2 V0
1
T=  (1.98)
√

1 2
1+ sin (λ e − 1)
4e(e − 1)

De igual forma:
−1
4e(e − 1)

R = 1+ √ (1.99)
sin2 (λ e − 1)
Se sabe que R + T = 1, R = 1 − T.

• Si E > V0 , e >> 1, T =∼ 1 y R = 0. A muy altas energías y potenciales débiles, las partículas no sienten el efecto de la
barrera, hay transmisión total.
1.6 La barrera de potencial 27
;A<

√ √ En
• Se tiene transmisión total cuando sin(λ e − 1) = 0 implica que λ e − 1 = nπ y T = 1 ocurre cuando en = V0 =
n2 π 2 h̄2 n2 π 2 h̄2
2mL2 V0
+ 1 o siempre que la energía incidente de la partícula sea En = V0 + 2mL2
con n = 1, 2, 3, · · · . El máximo de
el coeficiente de transmisión coincide con los autovalores de energía del pozo infinito de potencial. Estos son conocidos
como resonancias. Este fenómeno de resonancia, el cual no ocurre en física clásica, resulta de interferencia constructiva
entre las ondas incidente y reflejada (términos CeiK2 x y De−iK2 x en la Fig. 8). Este fenómeno se observa en un número
de casos tales como la dispersión a bajas energías (E ∼ 0,1ev) de electrones de átomos nobles (conocido como el efecto
Ramsauer-Townsend) y también en dispersión de neutrones fuera del núcleo.

√ √
• En el límite e → 1, sin(λ e − 1) ∼ λ e − 1 que da de (1.98) y (1.99) que:
−1
λ2 ( e − 1)

T = 1+
4e(e − 1)

E  −1
  
L 2 2mV − 1
0 
T=
 V 0 
1 + 2 E  E   
h̄ 4 − 1
V0 V 0
−1
L2 mV0

T = 1+ 2 ,
2h̄ E/V0
pero e → 1 → E ∼ V0 ,
−1
mL2 V0

T = 1+ (1.100)
2h̄2
y
!−1
2h̄2
R= 1+
mL2 V0

(ii) Caso E < V0 : efecto tunel.

1.6 La barrera de potencial

Clásicamente, se espera reflexión total, cada partícula que arriba a la barrera ( x = 0) se reflejará de nuevo, ninguna partícula
puede penetrar la barrera ya que tendrían una energía cinética negativa.
Las predicciones cuánticas son muy diferentes de las predicciones clásicas, ya que la función de onda no es cero más allá de la
barrera.
Las soluciones de la ecuación de Schrödinger producen soluciones similares a (1.83), (1.85) y (1.87) excepto que en (1.85), la
solución ψ2 ( x ) = CeiK2 x + De−iK2 x se cambia por ψ2 ( x ) = CeK2 x + De−K2 x .
Explícitamente se pueden escribir: 
ψ1 ( x ) = AeiK1 x + Be−iK1 x , x < 0, 


K2 x − K2 x
ψ2 ( x ) = Ce + De , 0 < x < L, (1.101)

iK1 x

ψ3 ( x ) = Ee , x>L


2mE 2m(V0 − E)
donde K12 = y K22 = .
h̄2 h̄2
El comportamiento de la densidad de probabilidad correspondiente a esta función de onda es la que se espera según la Fig. 11.
Es oscilatoria en las regiones x < 0 y x > L y decae exponencialmente en la región 0 6 x 6 L.
Se encontrarán los coeficientes de transmisión y reflexión:

| R |2 | E |2
R= 2
, T= (1.102)
| A| | A |2
1.6 La barrera de potencial 28
;A<

Aplicando las condiciones de continuidad para las funciones de onda en las regiones x = 0 y x = L, se sigue de (1.101):
en x = 0,

ψ1 (0) = ψ2 (0) −→ A + B = C + B (1.103)

 
dψ1  dψ2 
=
dx x=0 dx x=0
   
AiK1 eiK1 x − BiK1 e−iK1 x  = CiK2 eK2 x − DK2 eK2 x 
 
x =0 x =0

AiK1 − BiK1 = CK2 − DK2

iK1 ( A − B) = K2 (C − D ) (1.104)

En X = L,
ψ2 ( L) = ψ3 ( L) −→ CeK2 L + De−K2 L = EeiK1 L , (1.105)
 
dψ2  dψ3 
=
dx x= L dx x= L

CK2 eK2 L − DK2 e−K2 L = iK1 EeiK1 L

 
K2 CeK2 L − De−K2 L = iK1 EeiK1 L , (1.106)

Se combinan las ecuaciones (1.105) y (1.106) de la siguiente forma:

De−K2 L = EeiK1 L − CeK2 L

D = Ee(iK1 +K2 ) L − Ce2K2 L

K2 (CeK2 L − EeiK1 L + CeK2 L ) = iK1 EeiK1 L

K1 iK1 L
2C K2 L = i Ee + EeiK1 L
K2
 
E K
C= 1 + i 1 e(iK1 −K2 ) L (1.107)
2 K2
 
(iK1 +K2 ) L E K1 (iK1 +K2 ) L
D= Ee − 1+i e
2 K2
 
E K
D= 2 − 1 − i 1 e(iK1 +K2 ) L
2 K2
 
E K
D= 1 − i 1 e(iK1 +K2 ) L (1.108)
2 K2

Se combinan (1.107) y (1.108) con (1.103) y (1.104):

   
E K1 (iK1 −K2 ) L E K1 (iK1 +K2 ) L
A+B = 1+i e + 1−i e
2 K2 2 K2
 
B E 1 (iK1 −K2 ) L i K1 (iK1 −K2 ) L 1 (iK1 +K2 ) L i K1 (iK1 +K2 ) L
1+ = e + e + e − e
A A 2 2 K2 2 2 K2

E iK1 L eK2 L + e−K2 L K 1 e K2 L − e − K2 L

  
B
1+ = e −i
A A 2 K2 2
 
B E iK1 L K1
1+ = e cosh(K2 L) − i sinh(K2 L) (1.109)
A A K2
1.6 La barrera de potencial 29
;A<

K2
A − B = −i (C − D )
K1
    
K E K K
A−B = −i 2 1 + i 1 e(iK1 −K2 ) L − 1 − i 1 e(iK1 +K2 ) L
K1 2 K2 K2
     
iK1 L i K2 K 1 − K2 L i K2 K 1 K2 L
A−B = Ee − 1+i e + 1−i e
2 K1 K2 2 K1 K2

e − K2 L i K 2 K2 L e K2 L
 
B E iK1 L iK
1− = e − 2 e − K2 L + + e +
A A 2 K1 2 2 K1 2

E iK1 L eK2 L + e−K2 L e 2 − e − K2 L

  KL 
B K
1− = e +i 2
A A 2 K1 2
 
B E iK1 L K2
1− = e cosh(K2 L) + i sinh(K2 L) (1.110)
A A K1
B E
Se resuelven (1.109) y (1.110) para A yA .
Primero se suman (1.109) y (1.110):
 
E iK1 L K1 K2
2= e cosh(K2 L) − i sinh(K2 L) + cosh(K2 L) + i sinh(K2 L)
A K2 K1
   
E K2 K
2 = eiK1 L 2 cosh(K2 L) + i − i 1 sinh(K2 L)
A K1 K2
" #
E iK1 L K22 − K12
2= e 2 cosh(K2 L) + i sinh(K2 L)
A K1 K2
" #−1
E −iK1 L K22 − K12
= 2e 2 cosh(K2 L) + i sinh(K2 L) (1.111)
A K1 K2

De (1.109) y (1.111):
 
K1
2 cosh(K2 L) − i sinh(K2 L)
B K2
= −1
A K22 − K12
2 cosh(K2 L) + i sinh(K2 L)
K1 K2
" #" #−1
B K1 K22 − K12 K22 − K12
= ( (((K2 L) − 2i
2 cosh sinh(K2 L) − ( (((K2 L) − i sinh(K2 L) 2 cosh(K2 L) + i sinh(K2 L)
(
2 cosh
( ( ( (
(
A K2 K1 K2 K1 K2
!" #−1
B K K2 − K12 K2 − K12
= sinh(K2 L) −2i 1 − i 2 2 cosh(K2 L) + i 2 sinh(K2 L)
A K2 K1 K2 K1 K2
!" #−1
B 2K12 K2 + K2 (K22 − K12 ) K2 − K12
= −i sinh(K2 L) 2 cosh(K2 L) + i 2 sinh(K2 L)
A K1 K22 K1 K2
" #−1
B K12 + K22 K22 − K12
= −i sinh(K2 L) 2 cosh(K2 L) + i sinh(K2 L) (1.112)
A K1 K2 K1 K2
1.6 La barrera de potencial 30
;A<

Utilizando (1.111) y (1.112) en (1.102), se hallan los coeficientes de transmisión y reflexión

!2  
2 − K2 2
!
| B |2 K12 + K12 K
R= = sinh2 (K2 L) 4 cosh2 (K2 L) 2 1
sinh2 (K2 L)
| A |2 K1 K2 K1 K2

 !2 −1
| B |2 2 K22 − K12 2
T= = 44 cosh (K2 L) + sinh (K2 L) (1.113)
| A |2 K1 K2

Se utiliza la identidad cosh2 (K2 L) = 1 + sinh2 (K2 L) y se puede escribir (1.113) como:
 !2 −1
K22 − K12
T = 44 + 4 sinh2 (K2 L) + sinh2 (K2 L)
K1 K2

  !2  −1
K22 − K12 1
T = 1 + 1 + sinh2 (K2 L)
K1 K2 4

" " # #−1

K12 K22 + K24 + K14 − 2K12 K22 2
T = 1+ 4 sinh (K2 L)
4K12 K22
" ! #−1
1 K12 + K22 2
T = 1+ sinh (K2 L) (1.114)
4 K1 K2
!2
1 K12 + K22
R= T sinh2 (K2 L) (1.115)
4 K1 K2

T es finito, significa que la probabilidad para la transmisión de las partículas en la región x > L no es cero, clásicamente las
partículas no podrían estar en la región x > 0.
Este efecto es solo cuántico y es debido al aspecto ondulatorio de los objetos microscópicos.
Este efecto se conoce como efecto de tunelamiento: objetos mecánico-cuánticos pueden pasar a través de barreras clásicamente
impenetrables.
El efecto de tunelamiento tiene importantes aplicaciones en varias ramas de la física moderna como son: fusión nuclear, mi-
croscopio de efecto tunel, decaimiento alfa transporte de carga en mecanismos electrónicos.
2
2m(V0 − E)

!2 2mE
K12 + K22  h̄2 +
= h̄2 

K1 K2  2m p 
E(V0 − E)
h̄2
!2 !2
2E + 2(V0 − E) V0
= p = p
2 E(V0 − E) E(V0 − E)

V02
=
E(V0 − E)

Se pueden reescribir (1.114) y (1.115) como:

V02 T
 q 
L
R= sinh2 2m(V0 − E)
4E(V0 − E) h̄
" #−1
V02
 q
1 L
T = 1+ sinh2 2m(V0 − E)
4 E(V0 − E) h̄
1.7 Casos límites 31
;A<

Se hace la siguiente sustitución:

s
2mV0 E
λ=L , e=
h̄2 V0

T √
R= sinh2 (λ 1 − e), (1.116)
4e(1 − e)
−1
√

1
T = 1+ sinh2 (λ 1 − e) (1.117)
4e(1 − e)

1.7 Casos límites

Si E << V0 , por lo tanto e << 1, λ >> 1.
√ 1 √
Se hace la aproximación sinh(λ 1 − e) ' eλ 1−e , con esta aproximación (1.117) queda:
2
"  #−1
1 λ √1− e 2 √

1
T' e = 16e(1 − e)e−2λ 1−e
4e(1 − e) 2

  2L √
E E − 2m(V0 − E)
T ' 16 1− e h̄ (1.118)
V0 V0

Esto muestra que el coeficiente de transmisión no es cero, como sería clásicamente, y tiene un valor finito.
La mecánica cuántica predice que existe un tunelamiento más allá de la barrera para x > L.

Cuando E ' V0 , e = 1, (1.116) y (1.117) guian a la ecuación (1.100).

En el límite clásico, h̄ −→ 0, las ecuaciones (1.116) y (1.117) resultan: R −→ T, T = 0.

Ejemplo 1.1. Un electrón con energía de 30ev incide sobre una barrera cuadrada de altura 50ev.
Hallar el coeficiente de transmisión T, es decir, la probabilidad de encontrar al electrón en la región x > L. El ancho de la barrera es de
0.5nm.

Solución.
E = 30ev = 4,806 × 10−18 j
V0 = 50ev = 8,01 × 10−18 j, m = 9,11 × 10−31 kg
L = 0,1nm = 0,1 × 10−9 m

Se utiliza la ecuación (1.118):

2L
q
− 2m(V0 − E) = −4,56
h̄
El coeficiente de transmisión T es:

T∼
= 16(0,6)(1 − 0,6)e−4,56
T∼
= 0,040 = 4

Este resultado indica que de cada 100 partículas que inciden en la barrera 4 la tunelan, i.e., pasan a la región x > L.
¡¡¡Clásicamente esto es imposible!!!
Ejemplo 1.2. Los fenómenos cuánticos son frecuentemente despreciables en el mundo macroscópico.
Muestre esto numéricamente para los siguientes casos:

(a) La amplitud de oscilación de punto cero para un péndulo de longitud L = 1m y masa m = 1kg.
1.8 Un electrón en un pozo de potencial con barrera finitas 32
;A<

(b) La probabilidad de tunelamiento para una canica de masa m = 5g que se mueve a velocidad de 10cm/s contra un obstáculo rígido
de altura h = 5cm y ancho L = 1cm.

(c) La difracción de una pelota de tenis de masa m = 0,1kg que se mueve a velocidad v = 0,5m/s por una ventana de área 1 × 1,5m2 .
Solución.

(a) La energía cinética promedio K = 21 E.

Pero
1 1
E= KA2 = ω 2 A2
2 2
1 2 h̄ω
ω mA2 =
p
, ω = g/L = 31rad/s,
2 r 2
h̄
A= = 5,8 × 10−18 m.
mω
Esta oscilación es despreciable.

(b) Si se considera el ancho y altura del obstáculo rígido como el ancho y altura de una barrera potencial gravitacional, la probabilidad de
tunelamiento es:
De la ecuación (1.118),

E

E
 √
T∼ = 16 1− e−2L/h̄ 2m(V0 − E)
V0 V0
E: energía de la canica es cinética.
1
E = mv2 = 2,5 × 10−5 j
2
V0 = mgh = 2,45 × 10−3 j
29
T∼
= 0,1616e−9,3×10 ≈ 0,

la probabilidad de tunelamiento para la canica es cero.

(c) La longitud de onda de Broglie para la bola de tenis es:

h h
λ= = = 1,32 × 10−32 m,
p mv

y los ángulos de difracción en la dirección horizontal y vertical son:

λ
θ1 ≈ = 1,3 × 10−32 rad
D
λ
θ2 ≈ = 9 × 10−33 rad.
L
No hay difracción.

1.8 Un electrón en un pozo de potencial con barrera finitas

La tecnología en la fabricación de nanoestracturas con propiedades físicas dadas permite la realización experimental del confi-
namiento de electrones dentro del pozo de potencial de varios tipos. Los más importantes en efectos prácticos son potenciales
unidimensionales de pozos rectangulares con una de las paredes de potencial finita o con ambas paredes finitas.

1.8.1. Pozo de potencial con una de las paredes finita

Sea la energía potencial dada por:
1.8 Un electrón en un pozo de potencial con barrera finitas 33
;A<

∞, x 6 0,


V ( x ) = 0, 0 < x < L, (1.119)

V0 , x > L


tal potencial se muestra en la Fig. 12.

Es ahora necesario distinguir entre dos diferentes tipos de movimiento electrónico en el pozo.
El primer tipo corresponde a E 6 V0 , y el electrón en el pozo de potencial puede tener estados ligados, los cuales corresponden
al conjunto discreto de núcleos de energía. El segundo tipo de movimiento corresponde a E > V0 y el electrón no tiene estados
ligados y puede ser encontrado en algún punto con coordenadas x > 0. Este tipo de movimiento corresponde a un espectro de
energía continuo. En esta sección se considera sólo el caso cuando E 6 V0 y el electrón es confinado en el pozo de potencial.
Para la región 0 6 x 6 L, la ecuación de Schrödinger es:

d2 ψ1 ( x ) 2mE
+ K12 ψ1 ( x ) = 0, donde K12 = , (1.120)
dx2 h̄2
para la región x > L,
d2 ψ2 ( x ) 2m(V0 − E)
− K22 ψ2 ( x ) = 0, donde K22 = (1.121)
dx2 h̄
m representa la masa del electrón.
Las soluciones de (1.120) y (1.121) son:

ψ1 ( x ) = AeiK1 x + Be−iK1 x 0<x<L

(1.122)
ψ2 ( x ) = CeK2 x + De−K2 x x>L

Donde A, B, C, D son constantes, las cuales pueden ser encontradas de las condiciones de frontera. En este caso las condiciones
de frontera son la continuidad de las funciones de onda y de sus derivadas evaluadas en las fronteras del pozo de potencial.
Aplicando las condiciones de frontera en x = 0 0 = A + B, puesto que en la región x 6 0 la función de onda ψ( x ) = 0. La
primera ecuación en (1.122) se puede escribir como:

ψ1 ( x ) = 2iA sin(K1 x ) (1.123)

La función de onda ψ2 ( x ) cuando x −→ ∞ debe ser finita, esto implica que C = 0.

Luego ψ2 ( x ) = De−K2 x .
En la frontera x = L las siguientes condiciones para la continuidad de las funciones de onda y sus derivadas y sus deben
satisfacer:

ψ1 ( L) = ψ2 ( L)
 
dψ1 ( x )  dψ2 ( x ) 
= ,
dx x= L dx x= L
esto produce,
2iA sin(K1 L) = De−K2 L , (1.124)
− K2 L
2iAK1 cos(K1 L) = − DK2 e (1.125)

Ahora se divide la ecuación (1.124) por la ecuación (1.125)

K1
tan(K1 L) = − , (1.126)
K2

(1.126) define el espectro de energía del electrón en el pozo de potencial.

Esta es una ecuación transcendental, cuya solución no puede ser obtenida analíticamente.
Para el análisis des estas ecuaciones se utilizan métodos gráficos junto con métodos númericos.
1.8 Un electrón en un pozo de potencial con barrera finitas 34
;A<

Se tratará de transformar (1.126), para ello se utiliza la identidad trigonométrica:

tan φ
sin φ = ± q , (1.127)
1 + tan2 φ

tan(K1 L) −K /K
sin(K1 L) = ± p = ±s 1 2
1 + tan2 (K1 L) K2
1 + 12
K2

K1
sin(K2 L) = ± q (1.128)
K12 + K22

Se sabe que
q
K1 = 2mE/h̄2
s
2m(V0 − E)
K2 = ,
h̄2

entonces el término en la raíz cuadrada de la ecuación (2.129) se puede escribir como:

2mE 2m(V0 − E) 2mV0

K12 + K22 = + = ,
h̄2 h̄2 h̄2
luego se puede escribir la ecuación (2.129) de la forma:

K1
sin(K1 L) = ± r = ±γK1 L, (1.129)
2mV0
h̄2
donde,
h̄
γ= p ,
2mL2 V0

la cual define la pendiente

F = (K1 L) = ±γK1 L
se hace la gráfica del lado derecho e izquierdo de la ecuación (1.129) como funciones del parámetro KL (véase Fig. 13). Los
puntos de intersección de la función seno con las líneas de f = ±γKL definen las raíces de la ecuación (1.129). De todas las
raíces de la ecuación (1.129) se escogen las raíces las cuales satisfacen la condición dada por la ecuación (1.126). Las raíces
elegidas deben satisfacer la condición tan(KL) < 0.
La función tangente es negativa solo en los siguientes intervalos de valores para KL donde las raíces pueden existir:
π
+ nπ ≤ KL ≤ π + nπ (1.130)
2
Se eligen las raíces de la ecuación (1.129) las cuales existen en el intervalo (1.130). En la Fig. 13 los intervalos de KL permitidos
son mostrados en líneas de color púrpura. De las expresiones halladas se puede identificar muy bien para cuales soluciones de
la ecuación (1.126) ella no existe. No habrían raíces ni estados de energía ligados del electrón en esos pozos de potencial. Se
toma una solución particular para la ecuación (1.129), de la Fig. 13 se ve la línea punteada que corresponde a la solución máxima
con γ = π2 cuando las raíces de la ecuación (1.126) existen para valores más altos de la pendiente γ, no hay estados ligados.
Así la condición para la cual hay, al menos, un estado de energía con E ≤ V0 se reduce a la igualdad:

2
γ≤ (1.131)
π
se sabe que
h̄ 2
γ= p 6 (1.132)
2mL2 V0 π
1.8 Un electrón en un pozo de potencial con barrera finitas 35
;A<

Esta condición se puede escribir más convenientemente como

π 2 h̄2
V0 L2 > , (1.133)
8m
donde el lado izquierdo de la desigualdad se define solo por dos parámetros del pozo de potencial, y el derecho por el tipo
de partícula confinada en el pozo de potencial. Relativamente en pozos de potencial estrechos y superficiales, puede que no
2
existan estados ligados para el electrón. El segundo nivel de energía en el pozo puede existir en γ 6 3π i.e.,

9π 2 h̄2
V0 L2 > (1.134)
8m
sean las expresiones para las funciones de onda en cada región (0 < x < L) y ( x > L) dadas por:

ψ1 ( x ) = 2iA sin(K1 L)
ψ1 ( x ) = A1 sin(K1 L), (1.135)
donde A1 = 2iA
y
ψ2 ( x ) = De−K2 x (1.136)

De la ecuación (1.124) D = A1 sin(K1 L)eK2 L .

Se reemplaza esta expresión para D en la ecuación (1.136)

ψ2 ( x ) = A1 sin(K1 L)eK2 L e−K2 L (1.137)

De (1.135) y (1.136) se sigue que dentro del pozo de potencial la función de onda es una función oscilatoria como lo fue en el
caso del pozo infinito de potencial. Afuera del pozo, el electrón se aparta de la barrera de potencial, su fucnción de onda decre-
se exponencialmente. Esto significa que en un estado ligado hay una probabilidad diferente de cero de encontrar al electrón
afuera del pozo donde su energía es E < V0 , y por tanto la energía cinética K sería negativa.
Esto puede ser explicado como sigue: en mecánica cuántica la energía total E, no puede ser presentada como una suma de
energía cinética y potencial, porque, de acuerdo al principio de incertidumbre de Heisenberg, las energías cinética y potencial
(como también el momentum y la coordenada) no pueden ser ambos medidos simultáneamente. Las relaciones de incerti-
dumbre relacionan las coordenadas de las partículas y sus correspondientes momentum, los cuales definen energías potencial
(coordenadas) y cinética (momentum).
Así, la situación con E < V0 es posible para una partícula cuantizada, y es muy común en varios problemas teóricos y aplicados.
E1
Ejemplo 1.3. Encuentre la relación V0 para un electrón en un pozo con el perfil de potencial descrito por la ecuación (1.119) cuando solo
3π
un estado ligado existe el cual corresponde al parámetro K1 L = 4 .
q
Solución. De la ecuación (1.120) K1 = 2mE 2 .
1
h̄
3π
Se emplea la relación K1 L = 4 .

3π
K1 =
4L
s
3π 2mE1 9π 2 2mE1
= 2
−→ =
4L h̄ 16L2 h̄2

9π 2 h̄2
E1 = (1.138)
32mL2
El parámetro γ, el cual define la pendiente de la función F = γK1 L como es mostrado en la Fig. 13, se encuentra utilizando la ecuación
(1.129):
√ √
sin(K1 L) 2/2 2 2
γ= = = (1.139)
K1 L 3π/4 3π
1.8 Un electrón en un pozo de potencial con barrera finitas 36
;A<

Pero √
h̄
2 2
γ= p =
2mL2 V0 3π

h̄2 4(2) 9π 2 h̄2

= −→ V0 = (1.140)
2mL2 V0 9π 2 16mL2
Se toma en cuenta la ecuación (1.138),
E1 1
= = 0,5
V0 2

Ejemplo 1.4. Encuentre la constante de normalización A1 que se muestra en las ecuaciones (1.135) y (1.137).
Solución. Puesto que el comportamiento del electrón en las regiones 1 y 2 es descrita por sus funciones de onda, la condición de normaliza-
ción debe ser de la forma:
Z∞ ZL Z∞
|ψ|2 dx = |ψ1 |2 dx + |ψ2 |2 dx = 1 (1.141)
0 0 L

Se sustituyen en esta condición las expresiones para las funciones de onda ψ1 y ψ2 :

ZL Z∞
2
| A1 | sin(K1 x )dx + | A1 |2 sin2 (K1 L)e2K2 L e−2K2 x dx = 1
0 L

Z∞
 L 
Z
| A1 |2  sin(K1 x )dx + sin2 (K1 L)e2K2 L e−2K2 x dx = 1,
0 L

resolviendo las integrales, se obtiene:

" #
2 L sin(2K1 L) sin2 (K1 L)
| A1 | − + =1
2 4K1 2K2
v
u 2
A1 = u
sin(2K1 L) sin2 (K1 L)
u
t
L− +
2K1 K2
 Capítulo

2 Ondas, operadores y postulados de la me-

cánica cuántica

2.1 Vista clásica del experimento de la doble rendija

E
n física clásica, las partículas y ondas exhiben comportamientos completamente diferentes. Mientras la descrip-
ción completa de una partícula requiere sólo un parámetro, el vector posición ~r (t), la descripción completa de
una onda requiere dos, la amplitud y la fase.
ˆ
Una onda plana tridimensional puede ser descrita por la función de onda Ψ(~r, t) = Aei(K0~r−ωt) = Aeiφ , donde
A es la amplitud y φ la fase de la onda K̂ es el vector de onda y ω la frecuencia angular.
El significado de Ψ es dado por la intensidad de la onda: I = |Ψ|2 .

(a) Una fuente de partículas clásicas.

Sea F una fuente de partículas clásicas, indestructibles que se acerca a la pantalla, en la Fig. 2.1, sólo la fuente F1 está
abierta. Sólo la rendija F1 está abierta y I1 (y) es la correspondiente intensidad observada en la pantalla (el número de
partículas por segundo que llegan al punto y sobre la pantalla).

Fuente: F

Pantalla
Fuente de
partículas clásicas
Figura 2.1. Sólo la rendija F1 está abierta.

En el segundo experimento, sea I2 (y) la intensidad que se recoge sobre la pantalla cuando sólo la rendija F2 se abre. (Véase
Fig. 2.2).
Si las fuentes F1 y F2 se abren simultáneamente, la intensidad total que se recoge sobre la pantalla detrás de las dos
2.1 Vista clásica del experimento de la doble rendija 38
;A<

Fuente: F

Fuente de
partículas clásicas

Pantalla

Figura 2.2. Sólo la rendija F2 está abierta.

rendijas debe ser igual a la suma de I1 y I2 , (véase Fig. 2.3).

I (y) = I1 (y) + I2 (y) (2.1)

Fuente: F

Pantalla

Figura 2.3. Ambas rendijas, F1 y F2 , son abiertas.

2.1 Vista clásica del experimento de la doble rendija 39
;A<

(b) Una fuente de ondas.

Sea F una fuente de ondas (ondas de luz o sobre el agua).
Sea I1 la intensidad recogida sobre la pantalla cuando sólo F1 se abre e I2 la intensidad cuando sólo F2 se abre (véanse Fig.
2.4 y 2.5).

Fuente
F

Pantalla

Figura 2.4. Sólo la rendija F1 está abierta.

F
Fuente

Pantalla

Figura 2.5. Sólo la rendija F2 está abierta.

Una onda es representada por una función compleja ψ y su intensidad es proporcional a su amplitud (por ejemplo en las
ondas que se producen en un estanque de agua, la amplitud será la altura del agua) al cuadrado I1 = |ψ1 |2 , I2 = |ψ2 |2 .
Cuando ambas rendijas son abiertas (véase Fig. 2.3) la intensidad total recogida sobre la pantalla muestra un patrón de
interferencia, por lo tanto no puede ser la suma de I1 y I2 . Las amplitudes, no las intensidades, deben adicionarse: la
2.2 Vista cuántica de partículas y ondas 40
;A<

amplitud total ψ es la suma de ψ1 y ψ2 ; y la intensidad es:

I = |ψ1 + ψ2 |2 = |ψ|2 + |ψ2 |2 + ψ1∗ ψ2 + ψ2∗ ψ1

I = I1 + I2 + 2<(ψ1∗ ψ2 )
p
I = I1 + I2 + 2 I1 I2 cos ϑ,
√
donde ϑ es la diferencia de fase entre ψ1 y ψ2 , y 2 I1 I2 es un término oscilante el cual es responsable por el patrón de
interferencia, (véase Fig. 2.6).

Pantalla

Figura 2.6. Ambas rendijas, F1 y F2 están abiertas.

Así la distribución de intensidad resultante no puede ser obtenida de I1 o de I2 solamente, ya que depende de una
fase ϑ la cual no puede ser medida cuando solo una rendija es abierta.
Conclusión. Clásicamente, las ondas presentan patrones de interferencia mientras las partículas no lo hacen. Cuando dos corrientes de
partículas no interactuantes se combinan en la misma región del espacio, sus amplitudes se suman pero sus intensidades no.

2.2 Vista cuántica de partículas y ondas

Se discute el experimento de la doble rendija con partículas materiales de carácter cuántico tales como electrones. Aquí la
fuente F arroja un chorro de electrones (véase Fig. 2.7). Sólo F1 está abierta y luego F2 se abre (véase Fig. 2.8).
Cuando ambas rendijas se abren se presenta un patrón de interferencia, los electrones parecen interferir consigo mismos, esto
muestra como si los electrones pudieran ir a través de cada rendija, F1 y F2 al mismo tiempo a pesar de ser discretos.
Alguien se podrá preguntar: si un electrón no se partir o dividir ¿Cómo puede viajar a través de las dos rendijas al mismo tiem-
po? Estos experimentos se llevaron a cabo también con chorros débiles de electrones, i.e., el patrón de interferencia observado
no depende de la intensidad (véase Fig. 2.9).
Cada electrón se envió uno por uno hasta que alcanzaban la pantalla de observación, i.e., si ambas rendijas están abiertas y
esperamos cierto tiempo para lograr impactos suficientes sobre la pantalla, el patrón de interferencia aparece de nuevo.
La pregunta crucial es: ¿Por cuál rendija pasó el electrón?.
Para responder a esto: un experimento puede ser realizado para ver por cual rendija entra el electrón. Para saber esto se instala
una fuente de luz detrás de la pared que contiene las rendijas como muestra la Fig. 2.10. Se acomoda un contador Geiger tal que
un electrón cuando golpea la pantalla se escucha un «clik» sobre el contador. Puesto que las cargas eléctricas dispersan la luz,
siempre que un electrón pasa por alguna de las dos rendijas, emitirá luz hacia nuestros ojos, siempre que se escuche un «clik»
2.2 Vista cuántica de partículas y ondas 41
;A<

Fuente de electrones

Pantalla

Figura 2.7. Sólo la rendija F1 está abierta.

Fuente de electrones

Pantalla

Figura 2.8. Sólo la rendija F2 está abierta.

sobre el contador, se verá un destello cerca de las redijas F1 o F2 ; pero nunca cerca de ambos a la vez. Después de recolectar los
varios conteos con ambas rendijas abiertas, se encuentra que la distribución es similar a la de las partículas clásicas en la Fig.
3. ¡¡¡El patrón de interferencia ha desaparecido!!! Pero si se apaga la fuente de luz, el patrón de interferencia aparece de nuevo.
De este experimento se concluye que el mero acto de medir, i.e., de saber por cual rendija pasan los electrones, esto afecta la
distribución sobre la pantalla de observación.
Claramente, los electrones son muy delicados: su movimiento se modifica cuando son observados. Esto es el principio de la
mecánica cuántica el cual dice que «la medida interfiere con los estados de objetos microscópicos». Se puede pensar en bajar
la intensidad de la fuente de luz, para no interferir «demasiado» con los electrones que pasan por las rendijas, en este caso,
algunos electrones pasaron a través de las rendijas sin ser vistos, porque no habían fotones disponibles en el momento correcto
para identificarlos. Este proceso es importante porque confirma que la luz tiene propiedades de partícula.
Dos perfiles de distribución son considerados de este experimento con fuentes de luz tenue, uno que corresponde a electrones
que fueron vistos y el otro a electrones que fueron vistos (pero que fueron detectados por el detector Geiger), la primera distri-
bución no presenta interferencias (i.e., es similar al experimento clásico), pero la segunda distribución muestra un patrón de
interferencia. Esto resulta del hecho que cuando los electrones no son vistos, ellos presentan interferencia. Para los electrones
2.2 Vista cuántica de partículas y ondas 42
;A<

que producen interferencia, es imposible identificar la rendija por la cual ha pasado cada electrón. Este hallazgo experimental
introduce un nuevo concepto fundamental: el mundo microfísico es indeterminado.
A diferencia de la física clásica, en donde se puede rastrear la trayectoria de las partículas, no se puede seguir las trayectorias
de las partículas microscópicas. Es técnicamente imposible hacer tal trazado, con detalle, de sus trayectorias. Tal resultado
inspiró a Heisenberg para postular sus famosas relaciones de incertidumbre, el cual dice: «es imposible diseñar un aparato el
cual permita determinar la rendija por la cual el electrón atraviesa sin perturbar bastante al electrón y destruir el patrón de
interferencia». El patrón de interferencia que se obtiene del experimento de la doble rendija indica que los electrones muestran

Pantalla

Figura 2.9. Ambas rendijas, F1 y F2 están abiertas.

un comportamiento como partícula y como onda.

Cuando los electrones son observados o detectados (por la fuente de luz) uno por uno, ellos se comportan como partículas,
pero cuando ellos son detectados después de muchas medidas, ellos se comportan como ondas de longitud de onda λ = hp y
muestran un patrón de interferencia.
Los diferentes hallazgos experimentales discutidos hasta acá: radiacción de cuerpo negro, efecto fotoeléctrico, efecto compton,
experimento de Davisson-Germer, experimento de la doble rendija, revelan que los fotones, electrones y cualquier otra partí-
cula microscópica se comportan diferente a partículas clásicas diferente a ondas clásicas. Estos descubrimientos indican que
en el mundo microscópico la naturaleza puede mostrar comportamiento de partícula como también un comportamiento on-
dulatorio.
La pregunta ahora es: ¿Cómo puede algo comportarse como partícula y como onda al mismo tiempo?.
Son estas ideas diferentes en el dominio de la clásica la respuestas es si, pero no en la mecánica cuántica. Este comportamiento
dual no puede ser reconciliado en la física clásica de ningún modo, son conceptos totalmente diferentes. Pero los sistemas mi-
croscópicos pueden abordar ambos comportamientos. Las manisfestaciones de onda y de partícula no contradicen ni excluyen
una a la otra, pero como es sugerido por Bhor, ellas son sólo complementarias. Ambos conceptos (descripción de partícula o de
onda de un fenómeno) son complementarios en describir la verdadera naturaleza de los sistemas microscópicos siendo rasgos
complementarios del mundo microscópico. Partículas y ondas son igualmente importantes para una descripción completa de
los sistemas cuánticos. De aquí proviene la esencia del principio de complementariedad.
Se ha visto que cuando el concepto rígido de partícula o onda es indiscriminadamente aplicado o impuesto sobre la física cuán-
tica, entramos en conflicto con la realidad, sin el principio de complementariedad, la mecánica cuántica no habría estado en
2.3 La función de onda 43
;A<

una posición para producir los resultados correctos que ha realizado.

Fuente de luz
para detectar
por cual rendija pasan Pantalla
los electrones

Figura 2.10. Ambas rendijas, F1 y F2 están abiertas; pero existe una fuente luminosa detrás de las rendijas para identificar por cual pasa el electrón prove-
niente de la fuente F de electrones.

2.3 La función de onda

E
s posible mejorar la comprensión de la física cuántica si se hace un enlance conceptual entre partículas y ondas
con la noción del uso probabilidad.
Al reconocer la dualidad entre onda y partícula de la radiación electromagnética y la materia, se dice que la
probabilidad de hallar la partícula es proporcional al cuadrado de la amplitud de una onda que representa la
partícula. En el capítulo anterior se aprendió que existe una onda de De Broglie asociada con cualquier partícula. La amplitud
de la onda de De Broglie asociada con una partícula no es una cantidad que se pueda medir, porque la función de onda que
representa una partícula es, generalmente, una función compleja y se denota con la letra griega Ψ y si la función depende de
la posición y el tiempo y ψ si solo depende de la posición. A esta función se le conoce como función de onda.
La función de onda ψ es con frecuencia una función compleja. El cuadrado absoluto |ψ|2 = ψ∗ ψ donde ψ∗ es el complejo
conjugado de ψ, es siempre real y positivo, y es proporcional a la probabilidad por unidad de volumen de hallar una partícula
en un punto determinado en algún instante. La función de onda contiene en su interior toda la información que se pueda
conocer de una partícula. Aunque ψ no se observar, se puede medir la cantidad real |ψ|2 . Si ψ representa una sola partícula,
entonces |ψ|2 , se considera como densidad de probabilidad, y si se tiene una partícula en el espacio, es la probabilidad relativa
en cada punto del espacio de encontrar la partícula.
Matemáticamente se puede interpretar esto como:

P( X, Y, Z )dV = |ψ|2 dV, (2.2)

2.4 Principio de superposición lineal 44
;A<

que significa la probabilidad de encontrar a la partícula en un elemento de volumen que rodea algún punto.1
De Broglie relaciona el momento de una partícula con su longitud de onda λ = Ph . Si una partícula libre ideal posee un momen-
tum Px conocido, su función de onda es una onda senosoidad infinitamente larga de longitud de onda λ = Phx y la partícula
tiene la misma probabilidad de estar en cualquier punto a lo largo del eje x, (véase Fig. 2.11).

Figura 2.11. Distribución de probabilidad de De Broglie para una partícula de materia.

La función de onda para tal partícula libre que se mueve a lo largo del eje x se puede escribir como:

ψ( x ) = AeiKx , (2.3)
donde A es la amplitud constante y K = 2π λ es el número de onda angular.
Se expande la zona con sistemas unidimensionales, tal que (2.2) queda P( x ) = |ψ|2 dx: la probabilidad de que la partícula se
encuentre en el intervalo infinitesimal alrededor del punto x. No es posible especificar la posición de una partícula con certeza
completa, pero por medio de |ψ|2 se deduce la probabilidad de observarla en una región que rodea un punto x determinado.
La probabilidad de hallar la partícula en un intervalo arbitrario a ≤ x ≤ b es:
Zb
P= |ψ|2 dx (2.4)
a
2
La probabilidad es el área bajo la curva de |ψ| en función de x entre los puntos x = a y x = b.
Debido a que la partícula debe estar en algún lugar a lo largo del eje x, la suma de las probabilidades en todos los valores de x
debe ser 1:
Z∞
|ψ|2 dx = 1 (2.5)
−∞
Cualquier función de onda que satisfaga (2.5) se dice que está normalizada, la normalización quiere decir, simplemente, que la
partícula existe en algún punto en el espacio.
Una vez conocida la función de onda para una partícula, es posible calcular la posición promedio a la que se opera hallanda
después de muchas mediciones. Esta posición promedio se denomina valor esperado de x y es definido como:
Z∞
hxi = ψ∗ xψdx, (2.6)
−∞

además se puede extender esta definición para hallar el valor esperado de cualquier función F ( x ) asociado con la partícula:
Z∞
hF ( x )i = ψ∗ F ( x )ψdx (2.7)
−∞

2.4 Principio de superposición lineal

¿Cómo se calcula matemáticamente la existencia del patrón de interferencia en el experimento de la doble rendija con partí-
culas materiales tales como electrones? Una respuesta se ofrece por el principio de superposición. La interferencia resulta de
1 Esta interpretación probabilística de la función de onda fue sugerida primero por Max Born (1882-1970) en 1928. Erwin Schrödinger (1887-1961) propuso

una ecuación de onda que describe la forma en la que la función de onda evoluciona en el espacio y tiempo.
2.5 Paquetes de onda 45
;A<

la superposición de las ondas emitidas por las rendijas F1 y F2 . S las funciones ψ1 (~r, t) y ψ2 (~r, t) las cuales denotan las ondas
que emiten las rendijas F1 y F2 representan dos estados físicos posibles del sistema, entonces la superposición lineal es:
ψ(~r, t) = α1 ψ1 (~r, t) + α2 ψ2 (~r, t) (2.8)
También representan una posible descripción física del sistema α1 y α2 son constantes complejas. Este es el principio de super-
posición. La intensidad producida sobre la pantalla de observación debido a la rendija F1 es |ψ1 (~r, t)|2 y la debido solo a F2 es
|ψ2 (~r, t)|2 , cuando ambas rendijas son abiertas, la intensidad es:
|ψ(~r, t)|2 = |ψ1 (~r, t) + ψ2 (~r, t)|2 (2.9)
2 2
= |ψ1 (~r, t)| + |ψ2 (~r, t)| + ψ1∗ (~r, t)ψ2 (~r, t) + ψ1 (~r, t)ψ2∗ (~r, t),
donde es asterisco denota el conjugado complejo.
Note que (2.9) no es igual a la suma de |ψ1 (~r, t)|2 y |ψ2 (~r, t)|2 , esta ecuación contiene un término adicional ψ1∗ (~r, t)ψ2 (~r, t) +
ψ1 (~r, t)ψ2∗ (~r, t), este término es el que da lugar al patrón de interferencia.
El patrón de interferencia por tanto resulta de la existencia de un cambio de fase entre ψ1 (~r, t) y ψ2 (~r, t).
Se hace un resumen sobre los resultados del experimento de la doble rendija.
1. Las intensidades se suman para partículas clásicas: I = I1 + I2 .

2. Las amplitudes y no las intensidades se suman para partículas cuánticas: ψ(~r, t) = ψ1 (~r, t) + ψ2 (~r, t).
Esto da resurgimiento al patrón de interferncia.

3. Siempre que se intenta determinar experimentalmente el resultado de eventos individuales para partículas materiales
microscópicas (tal como intentar especificar la rendija por la cual un electrón pasa), el patrón de interferencia se destruye.
En este caso las intensidades se suman de la misma forma que se hizo con las partículas clásicas: I = I1 + I2 .

2.5 Paquetes de onda

¿Cómo describir una partícula dentro del contexto de la mecánica cuántica?
Una partícula la cual se localiza dentro de una cierta región del espacio puede ser denotada por una onda de materia cuya
amplitud es grande en esa región y cero fuera de ella.
Esta onda de materia debe entonces ser localizada alrededor de la región del espacio en el cual la partícula es confinada.
Una función de onda localizada se llama un paquete de onda.
Un paquete de onda consiste de un grupo de ondas de longitudes de ondas ligeramente diferentes con fases y amplitudes que
interfieren destructivamente en las otras regiones. Estos paquetes de onda son útiles en la descripción de partículas aisladas que
son confinadas a una cierta región espacial y son clave en entender la conexión entre mecánica cuántica y mecánica clásica.
Los paquetes de onda son una herramienta unificadora que pueden servir como explicación para la naturaleza dual onda y
partícula en maecánica cuántica.

2.6 Paquetes de onda localizado

Los paquetes de onda localizados se pueden construir superponiendo, en la misma región del espacio ondas con longitudes de
onda diferente, pero con faces y amplitudes tal que la superposición sea constructiva en la región deseada y destructiva fuera
de ella.
Esta superposición se lleva a cabo por medio de transformadas de Fourier. Se considera el paquete de onda unidimensional.
Este paquete describe una partícula clásica confinada en una región unidimensional. una partícula que se mueve en el eje x. Se
construye el paquete Ψ( x, t) por superponer ondas planas (sobre x) de diferentes longitudes de onda:
Z∞
1
Ψ( x, t) = √ φ(K )ei(Kx−ωt) dk, (2.10)
2π
−∞

φ(K ) es la amplitud del paquete de onda.

Se mira la evaluación del paquete de onda en un tiempo particular, se toma por simplicidad z = 0.
Z∞
1
ψ0 ( x ) = √ φ(K )eiKx dK, (2.11)
2π
−∞
2.6 Paquetes de onda localizado 46
;A<

donde φ(K ) es la transformada de Fourier de ψ0 ( x )

Z∞
1
φ(K ) = √ ψ0 ( x )e−iKx dx (2.12)
2π
−∞

Este paquete tiene la propiedad que se requiere de localización: el pico de |ψ0 ( x )| en x = 0 y desaparece lejos de x = 0.
Esto implica en el lenguaje de la interpretación probabilista de Born que la partícula tiene una probabilidad mayor de ser
encontrada cerca de x = 0 y una escasa oportunidad de ser encontrada lejos de x = 0. Lo mismo se aplica para el pico φ(K )
en K = 0, desaparece fuera de K = 0.

x K
O
 1/4  1/4 2
2
/L2 eiK0 x L2 2
Figura 2.12. Dos paquetes de onda localizados: ψ0 ( x ) = 2
πL2
e− x y φ( x ) = 2π e− L (K −K0 ) /4 . Ellas tienen un pico en x = 0 y
K = K0 respectivamente y desaparecen lejos de acá.

La Fig. 2.12 muestra un paquete de onda típico que tiene las propiedades de localización discutidas.
En resumen, la partícula se representa no por una única onda de De Broglie de frecuencia bien definida, sino por un paquete
de onda que se obtiene sumando un gran número de diferentes frecuencias.
La interpretación física del paquete de onda o amplitud se hace clara: ψ0 ( x ) es la función de onda o amplitud de probabilidad
para encontrar la partícula en la posición x, donde |ψ0 ( x )|2 da la densidad de probabilidad de encontrar la partícula en x y
P( x )dx = |ψ0 ( x )|2 dx da la probabilidad de encontrar la partícula entre x y x + dx.
De las ecuaciones (2.11) y (2.12), se ve que:

Z∞ Z∞
2
|ψ0 ( x )| dx = |φ(K )|2 dK, (2.13)
−∞ −∞

entonces si ψ0 ( x ) es normalizada lo mismo pasa con φ(K ). La función φ(K ) puede ser interpretada como una amplitud de pro-
babilidad en el espacio de los vectores K, para una partícula en el estado φ(K ) y |φ(K )|2 representa la densidad de probabilidad
para medir K que es el vector de onda de la partícula. La cantidad P(K )dK = |φ(K )|2 dK da la probabilidad de encontrar el
vector de onda de la partícula entre K y K + dK.
Se puede extraer información acerca del movimiento de la partícula expresando su correspondiente onda de materia en tér-
minos de la energía de la partícula, E, y de su momentum P.
φ(K )
Se usa K = Ph̄ , dK = dP h̄ , E = h̄ω e introduciendo la definición φ ( P ) =
e √ .
h̄
2.6 Paquetes de onda localizado 47
;A<

Se escriben (2.10), (2.11) y (2.12) como:

Z∞
1 e( P)ei( Px− Et)/h̄ dP,
Ψ( x, t) = √ φ (2.14)
2πh̄
−∞

Z∞
1 e( P)eiPx/h̄ dP,
ψ0 ( x ) = √ φ (2.15)
2πh̄
−∞

Z∞
1
e( P) = √
φ ψ0 ( x )e−iPx/h̄ dx, (2.16)
2πh̄
−∞

donde E( P) es la energía total de la partícula descrita por el paquete de onda Ψ( x, t) y φ

e( P) es la amplitud de momentum del
paquete.
Ejemplo 2.1.
2 2
(a) Encuentre Ψ( x, 0) para un paquete de onda Gaussiano φ(K ) = Ae− L (K −K0 ) /4 , donde A es el funtor de normalización que será
encontrado.
Calcule la probabilidad de encontrar la partícula en la región − L2 ≤ x ≤ L2 .

(b) Encuentre φ(K ) para un paquete de onda cuadrado


 AeiK0 x , | x | ≤ L,
ψ0 ( x ) =
0, | x | > L

encuentre el funtor A tal que ψ( x ) sea normalizado.

Solución.
(a) El funtor de normalización A se obtiene de
Z∞ Z∞
2 2 (K −K 2
1= |φ(K )| dK = A 2
e− L 0) /2
dK
−∞ −∞

Se hace cambio de variable,

z = K − K0
Z∞ √
− L2 z2 /2 2π
e dZ = ,
2
−∞

esto da como resultado que

1/4
L2

A= ,
2π
el paquete de onda se escribe de la forma siguiente:
1/4
L2

2 /4( K − K 2
φ(K ) = e− L 0)
2π
De (2.11),
Z∞
1
ψ0 ( x ) = √ φ(K )eiKx dK
2π
−∞
1/4 Z∞
L2

1 2 /4( K − K 2
ψ0 ( x ) = √ e− L 0) +iKx
dK (2.17)
2π 2π
−∞
2.6 Paquetes de onda localizado 48
;A<

Se hace un arreglo en el exponente de (2.17):

L2 ix 2 x2
 
2 L
− (K − K0 ) + iKx = − (K − K0 ) − − 2 + iK0 x
4 2 L L

y se introduce una nueva variable de sustitución:

L x 2
y= ( K − K0 ) − i y dK = dy
2 L L
Entonces (2.17) se escribe commo:
1/4 Z∞
L2

1 2 /L2 2 2
ψ0 ( x ) = √ e− x eiK0 x e−y dy
2π 2π L
−∞

Z∞
2 1/4 − x2 /L2 iK0 x
 
1 2
ψ0 ( x ) = √ 2
e e e−y dy (2.18)
2π πL
−∞

R∞ 2 √
Se utiliza la integral e−y dy = π, entonces (2.18) es:
−∞
 1/4
2 2 /L2
ψ0 ( x ) = e− x eiK0 x (2.19)
πL2

donde eiK0 x es la fase de ψ0 ( x ); ψ0 ( x ) es una onda oscilante con número de onda K0 modulado por un envolvente Gaussiano centrado
en el origen (véase Fig. 12).
Este funtor de fase eiK0 x tiene significado físico real.
La función de onda ψ0 ( x ) es compleja, como lo necesita la mecánica cuántica. Tanto que ψ0 ( x ) como φ(K ) son normalizados la
L
probabilidad de encontrar la partícula en la región − L2 ≤ x ≤ se obtiene de (2.4)
2
ZL/2  1/2 ZL/2
2 2 2 /L2
P= |ψ0 ( x )| dx = e−2x dx
πL2
− L/2 − L/2

haciendo el cambio de variable z = 2x

= L2 dz,
L , dx
Z1 Z1
2 1/2 L
 
−z2 /2 1 2 2
P= 2
e dz = √ e−z /2 dz '
πL 2 2π 3
−1 −1

(b) La normalización para ψ0 ( x ) es:

Z∞ Z∞ ZL
2 2 −iK0 x iK0 x 2
1= |ψ0 ( x )| dx = A e e dx = A dx = 2LA2
−∞ −∞ −L

1
A= √
2L
Entonces la transformada de Fourier es:
Z∞
1
φ(K ) = √ ψ0 ( x )e−iKx dx
2π
−∞

ZL
1 1 sin [(K − K0 ) L]
φ(K ) = √ eiK0 x e−iKx dx = √
2 πL πL K − K0
−L
2.7 Paquetes de onda y las relaciones de incertidumbre 49
;A<

2.7 Paquetes de onda y las relaciones de incertidumbre

Se quiere mostrar aquí que el ancho de un paquete de onda ψ0 ( x ) y el ancho de su amplitud φ( x ) no son independientes sino
que están correlacionados. El ancho en los espacios de x y K tienen una conexión directa con el principio de incertidumbre de
Heisenberg.
Los paquetes de onda que se trataron en el ejemplo previo fueron:
1/4 1/4
L2
 
2 /4( K − K 2 2 2 /L2
φ(K ) = e− L 0) y ψ0 ( x ) = e− x eiK0 x (2.20)
2π πL2

como se muestra en la (Fig. 2.12), |ψ0 ( x )|2 y |φ(K )|2 son centrados en x = 0 y K = K0 .
Se definen los anchos medios ∆x y ∆K correspondientes al máximo medio de |ψ0 ( x )|2 y |φ(K )|2 . De esta forma cuando se varía
de 0 a ±∆x y K de K0 a K0 ± ∆K, las funciones |ψ0 ( x )|2 y |φ(K )|2 descienden a e−1/2 :

|ψ0 (±∆x )|2

= e−1/2 ,
|ψ0 (0)|2
(2.21)
|φ(K0 ± ∆K )|2
= e−1/2
|φ(K0 )|2

Estas ecuaciones combinadas con (2.20), conduden a:

2 /L2 2 /2∆K 2
e−2∆x = e−1/2 y e− L = e−1/2

esto conduce a
∆x = L
2 y ∆K = 1
L
∆P
por lo tanto ∆x∆K = 12 , puesto que ∆K =
h̄
h̄
∆x∆P = (2.22)
2
Esta relación demuestra que si el ancho del paquete de onda es estrecho en el espacio de las x su ancho en el espacio de los
momentos debe incrementar y viceversa.

2.8 Espacios de Hilbert

El formalismo de la mecánica cuántica trata con operadores que son lineales y funciones de onda que pertenecen a un espacio
de Hilbert abstracto. Las propiedades matemáticas y estructuras de espacio de Hilbert son esenciales persé un entendimiento
propio del formalismo de la mecánica cuántica.
Sea L un espacio lineal sobre R o C. Una aplicación L × L −→ R, C denotada por φ, que cumple:

(i ) (φ, φ) > 0 y (φ, φ) = 0 si y solo si φ = 0.

(ii ) (φ, φ1 + φ2 ) = (φ, φ1 ) + (φ, φ2 ),
(2.23)
(iii ) (φ, λψ) = λ(φ, ψ).
(iv) (ψ, φ)∗ , ∀φ, ψ, φ1 , φ2 ∈ L y ∀λ ∈ R o C.

en el caso de R sobran las conjugaciones. A L se le llama espacio con producto escalar. Un espacio de Hilbert (H ) es un
espacio con producto escalar, completo con la norma k φ k= (φ, φ)1/2 .
Desde un punto de vista geométrico, la gran ventaja de la estructura de producto escalar radica en la posibilidad que ofrece de
introducir generalizaciones naturales del concepto de ortogonalidad o perpendicularidad de la geometría griega clásica.
Dos vectores φ, φ2 ∈ H son ortogonales si (φ1 , φ2 ) = 0 (su producto escalae es cero), del mismo modo un conjunto de vectores
B = {φ1 } ⊂ H se dice que son ortogonales si (φi , φj ) = 0, ∀i 6= j, si además (φi , φi ) = 1, ∀i se dice que el conjunto B es
ortonormal.
Si dos conjuntos B1 y B2 de vectores en H se dice que B1 ⊥ B2 si (φ1 , φ2 ) = 0, ∀φ1 ∈ B1 y ∀φ2 ∈ B2 .
Se resumirá que es un espacio de Hilbert (H ).
Este espacio satisface las cuatro propiedades:
2.9 Dimensión y base de un espacio vectorial 50
;A<

i H es un espacio lineal.

ii H tiene definido un producto escalar que es estrictamente positivo, i.e., (ψ, ψ) > 0 cumple las propiedades de (2.23).

iii H es separable. Existe una sucesión de Cauchy con ψn ∈ H (n = 1, 2, 3, · · · ) tal que para ψ ∈ H y e > 0 existe al
menos un ψn de la sucesión para la cual k ψ − ψn k< e.

iv H es completo.
Cada sucesión de Cauchy φn ∈ H converge a un elemento de H . Para cualquier ψn , la relación lı́m k ψn − ψm k= 0
n,m→∞

En los espacios de Hilbert puede ser útil el siguiente esquema gráfico para recordar con facilidad el contexto en que se desen-
vuelven:

Figura 2.13

se puede decir que todo espacio de Hilbert es un espacio de Banach; pero solo algunos espacios de Banach son espacios de Hil-
bert.
Ejemplos de espacios de H son Rn y Cn .
En Rn y Cn el estudiante ha efectuado ya infinidad de cálculos y de ejercicios de geometría griega antigua con el producto
n
ai∗ bi
P
escalar habitual i.e., ψ = ( a1 , a2 , a3 , · · · , an ) y φ = (b1 , b2 , b3 , · · · , bn ) el producto escalar (ψ, φ) =
i =1

2.9 Dimensión y base de un espacio vectorial

Un conjunto de N vectores no cero φ1 , φ2 , φ3 , · · · , φN se llaman linealmente independientes (L.I) si y solo si la solución de la
N
P
ecuación: ai φi = 0 y se cumple que a1 = a2 = · · · = a N = 0 pero si existe un conjunto de escales no todos cero entonces
i =1
N
P
podemos expandir un vector ψ1 en la base {φi }, de la forma: ψn = ai φi .
i =1
Generalmente se trabajará en una base ortonormal como se definió anteriormente, esta base cumple con:
(
δij = 0, si i 6= j
(φi φj ) = δij en donde
δ1j = 1, si i = j

y esta base es completa, i.e., genera el espacio en donde está.

2.10 Funciones de cuadrado integrable: funciones de onda

En el caso de un espacio de funciones un elemento vector se da por una función compleja y el producto escalar por integrales,
el producto escalar de las funciones ψ( x ) y φ( x ) es:
Z
(ψ( x ), φ( x )) = ψ∗ ( x )φ( x )dx (2.24)
2.11 Notación de Dirac 51
;A<

Si esta integral diverge, el producto escalar no existe.

Y son de cuadrado integrable si: Z
(ψ( x ), ψ( x )) = |ψ( x )|2 dx, (2.25)

es finita.
La dimensión de un espacio de Hilbert H de funciones de cuadrado integrable es infinita, puesto que cada función de onda
puede ser expandida en términos de un número infinito de funciones linealmente independientes.
La dimensión del espacio se da por el máximo número de vectores base requeridos para expandir el espacio.
Como ejemplos de funciones de onda de cuadrado integrable están las funciones de onda de la mecánica cuántica Ψ(~r, t).
De acuerdo a la interpretación probabilística de Born para (~r, t), la cantidad | (~r, t) |2 dv representa la probabilidad de encontrar,
en un tiempo t, la partícula en un elemento de volumen dv. Pero la ondulatoriedad en una región es la probabilidad de encontrar
el objeto en esa región. Cuidado: no es la probabilidad de que el objeto esté ahí. ¡Hay una gran diferencia! El objeto no estaba
allí antes de que lo encontráramos. Es nuestra localización del objeto lo que causó que estuviera allí.
Siguiendo esto, la probabilidad de encontrar la partícula en un lugar en el espacio debe ser igual a 1.
Z Z∞ Z∞ Z∞
2
| (~r, t) | dv = dx dy dz| (~r, t) |2 dv = 1 (2.26)
−∞ −∞ −∞

Entonces las funciones de onda de la mecánica cuántica son de cuadrado integrable o normalizadas. Cualquier función de onda
que no sea de cuadrado integrable no tiene sentido físico en la mecánica cuántica.

2.11 Notación de Dirac

El estudio físico de un sistema es representado en mecánica cuántica por elementos de un espacio de Hilbert (H ), estos ele-
mentos son llamados vectores estado.
Dirac introdujo lo que se consolidó como una notación invaluable en mecánica cuántica.
Permite manipular el formalismo de la mecánica cuántica con facilidad y claridad. Dirac introdujo los conceptos de kets (| i)
y de bras (h |) y el de brakets (h | i).

kets: son elementos de un espacio vectorial, Dirac denotó el vector estado ψ por el símbolo |ψi.
kets∈ H , al espacio ket.

bras: son elementos del espacio dual de Hilbert Hd para cada ket existe un bra∈ Hd .
El producto escalar se define como:
(φ, ψ) = hφ|ψi
Z
hφ|ψi = φ∗ ( x, t)ψ( x, t)dx (2.27)

Si hφ|ψi∗ = hψ|φi, si |ψi y |φi son complejos; pero si son reales, hφ|ψi = hψ|φi.
Para cualquier par de estados |ψi y |φi del espacio de H se tienen las siguientes propiedades:
i La norma se define como hψ|ψi = 1.

ii Desigualdad de Schwarz:
|hψ|φi|2 6 hψ|ψihφ|φi (2.28)

iii Desigualdad triangular: para estos estados se tiene que:

q q q
hψ + φ|ψ + φi 6 hψ|ψi + hφ|φi (2.29)

iv Estados ortogonales: dos kets |ψi y |φi son ortonormales si cumplen hψ|φi = 0, propiedad de ortogonalidad, y
hψ|ψi = 1, hφ|φi = 1
2.11 Notación de Dirac 52
;A<

v Conmutatividad respecto a la suma:

|ψi + |φi = |φi + |ψi

vi Vector nulo: |ψi + |0i = |ψi.

vii El vector opuesto se define como |ψi = −|ψi tal que:

|ψi + |ψi = 0

viii Multiplicación por un escalar: α|ψi = |αψi

ix (|ψi)† = hψ| y (hψ|)† = |ψi

x |ψ + φi = |ψi + |φi
Ejemplo 2.2. (En notación de matrices los kets son representados por columna y los bras por matrices fila).
Sean los dos kets:    
−3i 2
| ψ i =  2 + i  y | φ i =  −i 
4 2 − 3i

(a) Encuentre el bra hφ|.

(b) Evalúe el producto escalar hφ|ψi

Solución.


(a) hφ| = 2 i 2 + 3i .
 

−3i
(b) hφ|ψi = 2i 2 + 3i  2 + i  = −6i + 2i − 1 + 8 + 12i = 7 + 8i
4

Los productos hψ|φi o hφ|hψ| no están definidos.

Ejemplo 2.3. Considere los estados:
|ψi = 3i |φ1 i − 7i |φ2 i y |χi = −|φ1 i + 2i |φ2 i
donde |φ1 i y |φ2 i son ortonormales.

(a) Calcule |ψ + χi y hψ + χ|.

(b) Calcule el producto escalar:

hψ|χi y hχ|ψi, ¿son iguales?

Solución.

(a)

|ψ + χi = |ψi + |χi = 3i − 1 |φ1 i − 5i |φ2 i


hψ + χ| = −3i − 1 hφ1 | + 5i hφ2 |

(b)
hψ|χi = −14 + 3i y hχ|ψi = −14 − 3i

Ejercicio 2.1. Muestre que los estados |ψi y |χi satisfacen la desigualdad triangular y la desigualdad de Schwarz.
2.12 Operadores 53
;A<

2.12 Operadores
b es una regla matemática que cuando se aplica a un estado ket |ψi, la transforma en otro estado ket |ψ0 i del
Un operador A
mismo espacio
b|ψi = |ψ0 i, hψ| A
A b = hψ0 | (2.30)
Ejemplos de operadores.

I |ψi = |ψi
Operador unitario: b (2.31)
→ ∂ ∂ ∂
Operador gradiente: ∇ ψ(~r ) = bi ψ + bj ψ + b
k ψ (2.32)
∂x ∂y ∂z
∂
Operador momentum: Pb|ψi = −ih̄ |ψi (2.33)
∂x
Operador puridad: ℘
b ψ( x ) = ψ(− x ) (2.34)
b |ψi = x |ψi
Operador posición: X (2.35)

Los operadores fundamentales, por decir algo, son:

→ ∂
b y Pb = −ih̄ ∇
X o en una dimensión Pb = −ih̄
∂x

muchos operadores se construyen de X

b y de P,
b por ejemplo el operador hamiltoniano para una sola partícula se define:

b2
b = P +V
H b (X
b ), (2.36)
2m
b2
P
donde 2m es el operador energía cinética K
byV
b es el operador energía potencial.
Esto demuestra que se pueden sumar operadores para obtener otros operadores.
El producto de dos operadores se define cuando operan sobre una función, por ejemplo, sean los operadores A
byB
b y que operen
sobre la función F unidimensional.
Entonces la expresión:
A
bBbF ( x )

es una función, i.e., A

bBb es un nuevo operador C.b
El significado de A BF ( x ) quiere decir que B opera primero sobre F ( x ), dando una nueva función y A
b b b b luego opera sobre esta
nueva función.
Ejemplo 2.4. Sean

b=X
A b
y
b = Pb = −ih̄ d ,
B
dx
entonces,
A
bBbF ( x ) −→ X
b PbF ( x ),

se puede construir un nuevo operador PbX b F ( x ), como se mencionó arriba primero opera X b y luego P.
b Se compara el resultdo
de operar con el producto P X y X P sobre F ( x ):
b b b b
 
b ∂ F (x) − ∂ X
 
b Pb − PbX
X b F ( x ) = −ih̄ X b F (x)
∂x ∂
 
b (x) − ∂ X
= −ih̄ Xg b F (x) (2.37)
∂x
2.12 Operadores 54
;A<

b F ( x ) = x F ( x ) reemplazando estos resultados en (2.37):

b ( x ) = xg( x ) y X
donde Xg
 
b F ( x ) = −ih̄ xg( x ) − ∂ x F ( x )
 
b Pb − PbX
X
∂x
 
∂ ∂
= ih̄ X F ( x ) − x F ( x )
∂x ∂x

∂F ( x )
 
∂
= ih̄ X F ( x ) − x − F (x)
∂x ∂x

= ih̄F ( x ) (2.38)

y puesto que esto debe conservarse para cualquier función diferenciable F ( x ), se puede deducir que:

b Pb − PbX
X b = ih̄,

así se ha demostrado que el producto entre los operadores X b 6= X

b y Pb no conmuta ya que PbX b P.
b
Ya que operadores que combinan de la forma A A − B A 6= 0 y frecuentemente resurgen en cálculos de la mecánica cuántica,
b b b b
se acostumbra escribir esta resta entre multiplicación de operadores como:

[ A,
b Bb] = A
bBb−B
bAb (2.39)

y este término se define como el conmutador de A b con B.

b Se debe conservar el orden.
Si [ A, B] 6= 0, entonces se dice que A y B no conmutan.
b b b b
Si [ A,
b Bb] = 0 entonces A byB b se dice que conmutan.
Una operación de la forma [ A, b Bb] 6= 0, se llama relación de conmutación.
Por ejemplo,
[ X,
b Pb] = ih̄ (2.40)
se denomina relación de conmutación fundamental.
También se puede definir una relación de anticonmutación como

[ A,
b Bb] = A
bBb+B
bAb (2.41)

2.12.1. Propiedades de los conmutadores

(i ) [ A,
b Bb] = −[ B,
b Ab]

(ii ) [ A,
b bBb + cC
b] = b[ A,
b Cb] + c[ B,
b Cb]

(iii ) [ a A
b + b B,
b Cb] = a[ A, b] + b[ B,
b C b Cb]

(iv) [ A,b BbCb] = Bb[ A,

b Cb] + [ A,
b Bb]C
b (2.42)
(v) [ AbB,
b Cb] = A b[ B,
b Cb] + [ A,
b Cb] B
b

(vi ) [ A,
b [ B,
b Cb]] + [C,
b [ A,
b Bb]] + [ B,
b [ A,
b Ab]] = 0

esta última expresión se denomina identidad de Jacobi.

2.13 Autofunciones y autovalores de operadores 55
;A<

2.12.2. Adjunto de un operador

(i ) (B)† = B ∗ .Si B es una cantidad escalar. V. gr., (5i )† = −5i = (5i )∗

b† )† = A
(ii ) ( A b

(iii ) ( A
bBb)† = B
b† A
b†
b)† = A
(iv) (B A b† B ∗ = B ∗ A
b†

(v) ( Ab+B b)† = A

b† + Bb†

(vi ) hφ| A b† ψi = hφ A
b|ψi = hφ| A b† |ψi (2.43)
∗
(vii ) hφ| A
b|ψi = hψ| A |φi b†

(viii ) |B Aψ
b i = BA
b|ψi
b | = hψ|B ∗ A
(ix ) hB Aψ b†

2.12.3. Operadores hermíticos

(i ) Ab† = Ab operador hermítico (2.44)

(ii ) Ab† = − A
b operador antihermítico (2.45)
b es hermítico y |φi es un estado arbitrario, entonces:
(iii ) Si A


b2 |φi > 0
h φ | A
(2.46)
 b|φi∗ = hφ| A
hφ| A b|φi

2.13 Autofunciones y autovalores de operadores

Se ha definido que un operador es un símbolo matemático que cuando es aplicado a una función da una nueva función. Así si
se tiene una función F ( x ) y un operador A,
b entonces:
bF ( x ) = nueva función
A (2.47)
La nueva función puede ser proporcional a F ( x ) en ese caso, se puede escribir:
bF ( x ) = a f ( x )
A (2.48)
donde a es una constante de proporcionalidad si se da (2.48) F ( x ) se llama una autofunción del operador A
b y a el correspondiente
autovalor.
Ejemplo 2.5. Considere la función F ( x ) = ei(kx−ωt) esto corresponde a una onda viajera en la dirección x que opera sobre F ( x ) con el
operador momentum.
Solución.
∂
PbF ( x ) = −ih̄ F ( x ) = −ih̄(ik )ei(kx−ωt)
∂x
PbF ( x ) = h̄k F ( x )
y puesto que por de Broglie, la relación h̄k es el momentum P de la partícula, se tiene:
PbF ( x ) = PF ( x ) (2.49)
La ecuación (2.49) explica la elección de signo en la definición del operador momentum, cf. (2.34).
El valor esperado de un operador A b (se escribe como h Abi) con respecto a un estado |ψi del sistema se define como:

hψ| A
b|ψi
hA
bi = (2.50)
hψ|ψi
2.13 Autofunciones y autovalores de operadores 56
;A<

2.13.1. Operadores hermíticos

b se llama conjugado hermítico2 de A
El operador A b si:
Z Z
(Ab† ψ)ψdx = ψ∗ Aψdx
b (2.51)

Un operador es hermítico si cumple: Z Z

∗
( Aψ
b ) ψdx = ψ∗ Aψdx,
b (2.52)

En notación de Dirac, se escribe:

b|ψi∗ = hψ| A
hψ| A b|ψi, (2.53)
Ejemplo 2.6.

a Verifique que el operador X

b es hermítico.
De (2.52): Z Z Z Z
b )∗ ψdx = ( xψ)∗ ψdx = ψ∗ xψdx = ψ∗ Xψdx
( Xψ b

∂
b Verifique que el operador Pb = −ih̄ es hermítico.
∂x
De (2.52):
Z∞ Z∞
b )∗ ψdx = ih̄
( Pψ ∂
( ∂x
∗
ψ) ψdx (2.54)
−∞ −∞

La última expresión en (2.54) se resuelve por integración por partes:

dψ

U = ψ, dU =
 dx
dx
∗
dv = dψ dx, v = ψ∗

dx
Z∞ Z∞   Z∞
∗ ∗ ∞ ∂ψ
( Pψ) ψdx = ih̄|ψ ψ||−∞ +
b −ih̄ dx = ψ∗ Pψdx
b
∂x
−∞ −∞ −∞

El segundo término en la última expresión es cero por la condición de cuadrado integrable que debe cumplir la función de onda.

Pb2
c verifique que el operador energía cinética K
b= es hermítico.
2m
De De (2.52):
Z Z Z
b )∗ ψdx = 1
(0Kψ 2
( P ψdx ) =
b 1 b )∗ Pψdx
( Pψ b
2m 2m
Z Z
1
= ψ P ψdx = ψ∗ Kψdx
∗ b2 b
2m

d El hamiltoniano (2.36) es hermítico:

b=K
H b+V
b (X
b)

Aquí V
b es un operador hermítico por virtud de ser una función del operador hermítico X,
b y puesto que K
b se demostró en el
numeral ( c ) que es hermítico, entonces H
b es hermítico, i.e.,

b† = H
H b (2.55)
2 Otro nombre para conjugado hermítico es adjunto.
2.14 Algunos teoremas importantes 57
;A<

Se recuerdan varias propiedades importantes de los operadores:

• (A
b+B
b+C b + · · · +)† = A
b+ D b† + B
b† + C
b† + D
b† + · · · (2.56)
†
• (A
bBbC
bD b †C
b) = D b† B
b† A
b† (2.57)
†
• (|ψihφ|) = |φihψ| (2.58)

Para un operador antihermitiano se cumple:

b† = − B
B b o hψ| B
b|φi = −hφ| B
b|ψi

2.13.2. Operador proyección

Un operador Pb se dice que es el operador proyección si es hermítico y es igual a su propio cuadrado:

Pb† = P,
b Pb2 = Pb (2.59)

2.13.3. Operador unitario

Un operador U b −1 es igual a su adjunto U

b es unitario si su inverso, U b † , i.e.,

U b −1 o U
b† = U bUb† = U
b †U
b=b
I (2.60)

Donde bI es el operador identidad.

Se han comentado algunas propiedades de operadores y estados. Se está listo para discutir como encontrar los autovalores y
autovectores de un operador.
un vector estado |ψi es un autovector de un operador A
b si:

b|ψi = a|ψi,
A (2.61)

a es un número complejo llamado el autovalor de A. b

La ecuación (2.61) se llama la ecuación de autovalores.

2.14 Algunos teoremas importantes

Teorema 2.1. Para un operador hermítico, todos sus autovalores son reales y los autovectores que corresponden a diferentes autovalores
son ortogonales.
b† = A,
Es decir, A b las a en (2.61) son reales y hφn |φm i = δnm

Demostración. Si
b|φn i = an |φn i,
A (2.62)
multiplicando (2.62) por el bra hφm |:

hφm | A
b|φn i = an hφm |φn i (2.63)
b† = a∗ hφm |,
hφm | A (2.64)
m

multiplicando (2.64) por el ket |φn i

b† |φn = a∗ hφm |φn i
hφm | A (2.65)
Como A
b es hermítico, entonces b†
A = A,
b y se puede escribir la ecuación (2.65) como:

b|φn i = a∗ hφm |φn i

hφm | A (2.66)

Restando la ecuación (2.66) menos la ecuación (2.63):

( a∗ − an )hφm |φn i = 0 (2.67)


2.15 Mecánica ondulatoria y matricial 58
;A<

De (2.67) se debe tener en cuenta dos casos:

(i) Si m = n entonces hφm |φn i > 0 y se debe tener que an = a∗m . Aque los autovalores son reales.

(ii) Si m 6= n, en general se tiene que an 6= a∗m y se debe tener hφm |φn i = 0, i.e., |φm i y |φn i deben ser ortogonales.

La doble implicación se mantiene: un operador es hermítico si sus autovalores son reales. Se deja esta prueba como ejercicio.
Corolario 2.1. Los autovalores del operador hamiltoniano son reales. 

Puesto que por definición el hamiltoniano tiene las dimensiones de energía, por tanto los autovalores del hamiltoniano son los
las energías del sistema descrito por la función de onda.
Teorema 2.2. Los autoestados de un operador hermítico definen un conjunto completo de estados base mutuamente ortonormales. El
operador es diagonal en esta base con sus elementos diagonales iguales a sus autovalores. Este conjunto base es único si el operador no tiene
autovalores degenerados3 y es no único (de hecho infinito) si existen degeneraciones.

Demostración. Se deja como ejercicio la demostración. 

Teorema 2.3. Si dos operadores hermíticos, A b y B,
b conmutan y si Ab no tiene autovalores degenerados, entonces cada autovector de A
b es
también un autovector de B.
b
En resumen, se puede construir una base ortonormal común que se compone de los autovectores conjuntos de Ab y B.
b

Demostración. Puesto que Ab es hermítico sin autovalores degenerados a cada autovalor de A

b le corresponde un único autovector
(por teorema (2.2)). Sea:
Ab|φn i = an |φn i (2.68)
y como A
byB
b conmutan, entonces:

B
bAb|φn i = A
bBb|φn i o A b|φn i) = an ( B
b( B b|φn i) (2.69)

entonces ( Bb|φn i) es un autovector de Ab (cf. (2.61)) con autovalor an . Este autovector es único (los autovalores no son degene-
rados). El ket |φn i debe ser un autovector de B:b
b|φn i = bn |φn i
B (2.70)
Puesto que cada autovector de Ab es también autovector de B, b y viceversa, ambos operadores deben tener una base común. Esta
base es única, y está conformada de los autovectores conjugados de A b y B. b|φn i
b Ahora si an es degenerada, se puede decir que B
es un autovector de A con autovalor an . |φn i no es necesariamente un autovector de B. Si uno de los operadores es degenerado,
b b
existe un número infinito de conjuntos base ortonormales que son comunes a estos dos operadores, i.e., la base común existe
y no es única. 
Teorema 2.4. Los autovalores de un operador antihermítico son solo imaginarios o cero.

Demostración. Se deja como ejercicio para el lector. 

Teorema 2.5. Los autovalores de un operador unitario son números complejos de módulo igual a uno. Los autovectores de un operador
unitario que no tienen autovalores degenerados son mutuamente ortogonales.

Demostración. Se deja como ejercicio para el lector. 

 
2 i 0
Ejercicio 2.2.  3 1 5  , ¿Es esta matriz unitaria?
0 − i −2

2.15 Mecánica ondulatoria y matricial

La teoría de la mecánica cuántica trata en esencia de resolver el siguiente problema de autovalores:
b | ψ i = E | ψ i,
H (2.71)

donde H
b es el hamiltoniano del sistema. Esta ecuación, en general, no depende del sistema coordenado o representación.
3 Degenerados significa que tienen autovalores repetidos.
2.15 Mecánica ondulatoria y matricial 59
;A<

2.15.1. Mecánica matricial

La representación de la mecánica cuántica en una base discreta es igual al problema de hallar los autovalores de una matriz. La
representación de (2.71) en una base discreta {|φn i} conlleva al siguiente problema de autovalores matricial:

=0

(2.72)
o en forma simplificada, se puede escribir (2.72) como:
 
 Hn×n − E b
In×n  = 0
b 
(2.73)

Esta es una ecuación de orden n en E. La solución de (2.72) cede al espectro4 de energía del sistema: E1 , E2 , E3 , · · · , En .
Si se conoce este conjunto de autovalores E1 , E2 , E3 , · · · , En , se puede determinar los correspondientes autovectores,
|φ1 i, |φ2 i, |φ3 i, · · · , |φn i.
La disgonalización de la matriz H b de un sistema da el espectro de energía como también a los vectores estado del sistema.
Este procedimiento fue trabajado por Heisenberg e involucra solo cantidades matriciales. Esta formulación de la mecánica
cuántica se conoce como mecánica matricial.
El punto de partida de Heisenberg, en su intento para encontrar una fundamentación teórica a las ideas de Bohr fueron las
transiciones atómicas, donde Fmn = ( Em − En )/h̄ la cual da las frecuencias de la radiación asociada con las transiciones del
electrón de una órbita m a una órbita n (los saltos cuánticos de Schrödinger). Estas frecuencias Fmn pueden ser acomodadas
en una matriz cuadrada donde el elemento mn corresponde a la transición del estado cuántico m al estado n.
Otros ejemplos en donde se podría hacer representaciones matriciales, son el operador posición por una matriz posición; el
operador momentum por una matriz momentum; el operador energía por una matriz energía y así otras. Así, dentro del con-
texto de la mecánica matricial, se trata con cantidades que no conmutan, como predice la mecánica cuántica.
Esto es un rasgo esencial que distingue la mecánica matricial de la mecánica clásica, donde todas las canridades conmutan (se
pueden medir simultáneamente). La mecánica matricial de Heisenberg requiere la introducción de alguna maquinaria mate-
mática: espacios vectoriales lineales, espacios de Hilbert, álgebra de conmutadores y álgebra matricial. Esto es completamente
diferente de la maquinaria matemática de la mecánica clásica.
Ejercicio 2.3. Considere un sistema cuyo estado inicial |ψ(0)i0 y hamiltoniano son dados por:
   
3 3 0 0
1
|ψ(0)i =  0  , H =  0 0 5 
5
4 0 5 0
a Si una medida de la energía es realizada, ¿Qué valores se obtendrán y con qué probabilidad?

b Encuentre el estado del sistema |ψ(0)i en términos de la base {|φn i}.

c Encuentre la energía total del sistema en t = 0.

Ejercicio 2.4. Realice el ejercicio anterior con el estado |ψ(0)i y hamiltoniano dados por:

a    
3 1 0 0 1
1 0   0 0 1 0 
|ψ(0)i =  , H= 
5 4   0 1 0 0 
0 2 0 0 1

b    
1 1 0 0 i
 1   0 0 1 0 
|ψ(0)i = 
 −1  ,
 b=
H  0

1 0 0 
0 −i 0 0 1
4 Conjunto de autovalores del sistema.
2.16 La función delta de Dirac 60
;A<

2.16 La función delta de Dirac

La función delta de Dirac fue introducida por P. Dirac al final de la década de 1920 en un esfuerzo por crear una herramienta
matemática para desarrollar la teoría cuántica de campos. Puede ser considerada como una función de probabilidad o función
generalizada. En dinámica se conoce como una función impulso. Es una herramienta muy útil que aparece en muchas ramas
de la física: mecánica cuántica, electromagnetismo, óptica, problemas de ingeniería.
La ecuación básica asociada con la función delta de Dirac δ( x ) es:
Z∞
δ( x ) f ( x )dx = f (0), (2.74)
−∞

donde f ( x ) es alguna función que es continua en x = 0.

La ecuación (2.74) parece extraña: se tiene una integral que solo depende sobre el valor de la función f ( x ) en x = 0, ya que una
integral es una área bajo la curva, entoncces se espera que su valor no dependa de un valor particular de x. No hay función δ( x )
que satisfaga (2.74); pero existe otra clase de objetos matemáticos, conocidos como funciones generalizadas o distribuciones
que satisfacen (2.74).
Una función generalizada puede ser definida como el límite de una sucesión de funciones.
Se verá como trabaja esto en el caso de δ( x ).
Se comienza con la función Gaussiana normalizada:
r
n −nx2
gn ( x ) = e , (2.75)
π
1
aquí n = δ2
, con δ siendo el parámetro Gaussiano estandar. Esta clase de funciones son normalizadas ya que cumplen:

Z∞
gn ( x )dx = 1, (2.76)
−∞

para cualquier valor con n > 0, la sucesión de funciones para n = 1, 2, 3, · · · , ver v. gr., la (Fig. 2.14) son:
r r r
1 − x2 2 −2x2 100 −100x2
g1 ( x ) = e , g2 ( x ) = e , · · · , g100 ( x ) = e ,···
π π π

15

n=1
n=2
10
n=100
n=200
5 n=1000

-3 -2 -1 1 2 3

Figura 2.14. Función generalizada para n = 1, n = 2, n = 100, n = 200, n = 1000

Se puede afirmar que si n aumenta, se pueden escribir las cuatro propiedades siguientes:

(i) gn (0) es grande,

(ii) gn ( x 6= 0) eventualmente es pequeño,

(iii) El ancho del pico central es pequeño,

2.16 La función delta de Dirac 61
;A<

(iv)
Z∞
gn ( x )dx = 1
−∞

permanece constante.
La sucesión (2.75) es útil si aparece como parte de una integral, v. gr.,:
Z∞ Z∞ r
n −nx2
lı́m gn ( x ) f ( x )dx = lı́m e f ( x )dx, (2.77)
n→∞ n→∞ π
−∞ −∞

cuando n → ∞, gn ( x ) se estrecha en valores próximos a x = 0. Entonces para n → ∞ sólo f ( x ) en x = 0, es importante. Así

se puede reemplazar f ( x ) por f (0) en la integral, lo cual produce: (véase Fig. para un ejemplo).
Z∞ r Z∞ r
n −nx2 n −nx2
lı́m e f ( x )dx = f (0) lı́m e dx = f (0) lı́m 1 = f (0) (2.78)
n→∞ π n→∞ π n→∞
−∞ −∞

en (2.78) se utilizó la propiedad (iv).

15

10
g(x)
f(x)
5

-4 -2 2 4

1000 −1000x2
Figura 2.15. f ( x ) = cos( x ) y g( x ) = e
π
Entonces la función delta de Dirac es definida por la ecuación:
Z∞ Z∞ r
n −nx2
δ( x ) f ( x )dx = lı́m e f ( x )dx, (2.79)
n→∞ π
−∞ −∞

se podría decir, evitando formalismos, que: r

n −nx2
δ( x ) = lı́m e , (2.80)
n→∞ π
Esta es la delta de Dirac.
n −nx2
La ecuación (2.80) no tiene significado a menos que la función πe aparezca dentro de una integral y el límite lı́mn→∞
aparezca afuera de la misma integral.
Por ejemplo la (Fig. 2.15) representa la integral
Z∞ Z∞
δ( x ) f ( x )dx = δ( x ) cos( x )dx = 1 (2.81)
−∞ −∞

como se puede ver en dicha figura.

Sin embargo, cuando se trabaja con la función delta, raramente se necesita el límite que la define.
Esta sucesión particular de funciones definida en (2.75) no es única. Existen muchas sucesiones que pueden ser usadas para
definir la función delta. Por ejemplo también se podría haber definido δ( x ) como:
1 sin(nx )
δ( x ) = lı́m (2.82)
n→∞ π x
2.16 La función delta de Dirac 62
;A<

30

20
n=10
n=20
10 n=30
n=100

- 1.0 - 0.5 0.5 1.0

1 sin(nx )
Figura 2.16. Sucesiones de π x para n = 10, 20, 30, 100

sin(nx )
La sucesión de funciones x se ilustra en la (Fig. ). Las propiedades clave de las sucesiones en (2.82) fueron expresadas
previamente de la (i ) a la (iv).
Para sintetizar:

2.16.1. Propiedades de la función delta

(i)
f ( x ) · δ ( x ) = f (0) δ ( x )
f ( x )δ( x − a) = f ( a)δ( x − a) (2.83)
δ(− x ) = δ( x ), δ( x − a) = δ( a − x )

(ii)
ZB (
f (0), A < 0 < B
f ( x )δ( x )dx = (2.84)
0, 0 ∈ / ( A, B)
A

ZB (
f ( a ), A < a < B
f ( x )δ( x − a)dx =
0, a ∈ / ( A, B)
A

(iii)
ZB
1
f ( x )d(kx )dx = f (0) (2.85)
|k|
A

(iv)
n
X δ ( x − xi )
δ( f ( x )) = (2.86)
| f 0 ( xi )|
i =1

(v)
1
δ(kx ) = δ( x ) (2.87)
|k|

(vi)
1
δ ( x 2 − a2 ) = [δ( x − a) + δ( x + a)] (2.88)
2| a |

(vii)
1
δ [k( x − x0 )] = δ ( x − x0 ) (2.89)
|k|
2.17 La función delta tridimensional 63
;A<

(viii)  
1 b
δ [(kx ± b)] = δ x± (2.90)
|k| k

(ix)     
1 b b
δ(kx2 − b) = r δ x+ +δ x− (2.91)
b a a
2k
a

(x)
1
δ [( x − a)( x − b)] = [δ( x − a) + δ( x − b)] (2.92)
| a − b|
Z∞
1
δ( x ) = eikx dk transformada de Fourier
2π
−∞

(xi) Usando la transformada de Fourier de la función delta se escribe:

Z∞
dδ( x ) i
= δ0 ( x ) = keikx dk (2.93)
dx 2π
−∞

(xii)
Z∞
f ( x )δ0 ( x − a)dx = − f 0 ( a), (2.94)
−∞

Z∞
f ( x )δ00 ( x − a)dx = (−1)2 f 00 ( a) = f 00 ( a), (2.95)
−∞

Z∞
f ( x )δn ( x − a)dx = (−1)n f n ( a), (2.96)
−∞

Z∞
f ( x )δ0 ( x − a)dx = 0 (2.97)
−∞

(xiii)
δ0 ( x ) = −δ0 (− x ),
xδ0 ( x ) = −δ( x ),
(2.98)
x2 δ0 ( x ) = 0,
x2 δ00 ( x ) = 2δ( x )

2.17 La función delta tridimensional

La forma tridimensional de la función delta de Dirac es:
δ(~r − ~r 0 ) = δ( x − x 0 )δ(y − y0 )δ(z − z0 ), (2.99)
y en coordenadas esféricas es:
1
δ(~r − ~r 0 ) = 2 δ(r − r 0 )δ(cos θ − cos θ 0 )δ(φ − φ0 )
r
1
= δ (r − r 0 ) δ ( θ − θ 0 ) δ ( φ − φ 0 ) (2.100)
r2 sin θ
2.18 Representación posición 64
;A<

donde,
δ(θ − θ 0 )
δ(cos θ − cos θ 0 ) = (2.101)
sin θ
La transformada de Fourier de la función delta es:
Z
1 ~ ~0
δ(~r − ~r 0 ) = d3 keik·(~r−r ) , (2.102)
2π 3
Z
d3 r f (~r )δ(~r ) = f (0)
y
Z
d3 r f (~r )δ(~r − ~
r0 ) = f (~
r0 ) (2.103)

Las siguientes relaciones son frecuentemente encontradas:

 
−→ r
∇ · 2 = 4πδ
b
r
y
 
2 1
∇ = −4πδ(~r )
r
(2.104)

r es el vector unitario.
donde b
Se menciona que la dimensión física de la función delta de Dirac es uno sobre la dimensión del argumento. Por ejemplo, si x
1 1
es una distancia, la dimensión física de δ( x ) está dada por [δ( x )] = = , donde L es una longitud. De igual manera, la
[x] L
1
dimensión física de δ(~r ) es 3 puesto que:
L
1 1 1 1
[δ(~r )] = [δ( x )δ(y)δ(z)] = = 3
[ x ] [y] [z] L

2.18 Representación posición

En la representación posición, la base consiste de un número infinito de vectores {|~r i}, los cuales son kets del operador posición
~R,
b
~R
b
|~r i = ~r |~r i, (2.105)

donde ~r es el vector posición y es al mismo tiempo el autovalor del operador ~R.

b
La condición de completez para la base es:
Z
d3 r |~r ih~r | = b
I, (2.106)

y la condición de ortonormalidad es:

h~r |~r 0 i = δ(~r − ~r 0 ) como se muestra en (2.99) (2.107)
cada vector estado |ψi puede expandirse como:
Z Z
|ψi = d r |~r ih~r |ψi = d3 rψ(~r )|~r i,
3
(2.108)

en (2.108) se ha utilizado (2.106), y se tiene en cuenta que

h~r |ψi = ψ(~r ), (2.109)

2.19 Representación momentum 65
;A<

(2.109) se conoce como la función de onda para el vector estado |ψi.

El producto escalar entre dos vectores estado |ψi y |φi, puede ser expresado en la representación posición como:
Z 
3
hφ|ψi = hφ| d r |~r ih~r | |ψi
Z
= φ∗ (~r )ψ(~r )d3 r, (2.110)

De (2.105),
bn
~R |~r i = ~r n |~r i,
multiplicando por el bra h~r 0 |
h~r 0 |~R|~r i = ~r n h~r 0 |~r i = ~r n δ(~r 0 −~r ),
b
(2.111)

b†
en (2.111) se utilizó (2.105) el operador posición en tres dimensiones ~R es hermítico, i.e., ~R = ~R .
b b

2.19 Representación momentum

La base {|~
Pi} de la representación momentum se obtiene del operador momentum ~
P:
b

~P
b~
| Pi = ~P|~Pi, (2.112)

donde ~
P es el vector momentum.
La condición de ortonormalidad es:
~ 0 i = δ(~P − P
h~P| P ~ 0 ), (2.113)
y la condición de completez es: Z
d3 P | ~
Pih~
P| = b
I (2.114)

Expandiendo |ψi en esta base, se obtiene:

Z Z
3
|ψi = d P|~
Pih~
P||ψi = d3 Pψ(~
P)|~
P i, (2.115)

donde el coeficiente de expansión ψ(~P) representa la función de onda en el espacio de los momentos.
El producto escalar entre dos estados es:
Z  Z
hφ|ψi = hφ| d P| Pih P| |ψi = d3 Pφ∗ (~
3 ~ ~
P ) ψ ( P ), (2.116)

bn ~
~ 0 |~P ~ 0 − ~P)
hP | Pi = ~Pn δ( P (2.117)

para cambiar de la base {|~r i} a la base {|~

Pi} se introduce la función de transformación h~r |~
P i.
Z  Z
h~r |ψi = h~r | d3 P | ~ P | | ψ i = d3 P hr | ~
Pih~ Piψ(~
P ), (2.118)

i.e.,
Z
ψ(~r ) = d3 Ph~r |~
Piψ(~
P ), (2.119)

similarmente,
Z  Z
3
ψ(~
P) = h~
P|ψi = h~
P| d r |~r ih~r | |ψi = d3 r h ~
P|~r iψ(~r ) (2.120)

Las relaciones (2.119) y (2.120), ψ(~r ) y ψ(~

P), son la transformada de Fourier una de la otra.
En mecánica cuántica la transformada de Fourier de la función F (~r ) viene dada por:
Z
1 ~
F (~r ) = 3/2
d3 Pei P·~r/h̄ g(~
P ), (2.121)
(2πh̄)
2.19 Representación momentum 66
;A<

se compara (2.119) y (2.121), obteniéndose:

1 ~
h~r |~Pi = 3/2
ei P·~r/h̄ , (2.122)
(2πh̄)
esta función transforma de la representación momentum a la representación posición.
La función inversa h~ P|~r i es:
∗ 1 ~
h~P|~r i = h~r |~Pi = 3/2
e−i P·~r/h̄ (2.123)
(2πh̄)
 2
La cantidad h~r |~
Pi representa la densidad de probabilidad para encontrar la partícula en una región alrededor de ~r cuando
 

su momento es igual a ~ P.
Se utiliza (2.119) y (2.122) y se escribe la función de onda en la representación posición:
Z
1 ~
ψ(~r ) = 3/2
d3 Pei P·~r/h̄ ψ(~
P ), (2.124)
(2πh̄)
si Z
d3 rψ∗ (~r )ψ(~r ) = 1, es normalizada, entonces
Z
d3 Pψ∗ (~
P)ψ(~
P) = 1

este resultado se conoce como teorema de Parseval (su demostración se deja como ejercicio para el lector).
Para determinar la forma del operador momentum ~ P en la representación posición, se calcula:
b
Z Z
h~r | P|ψi = h~r | P| Pih P|ψid P = ~Ph~r |~Pih~P|ψid3 P
~ ~ 3
b b~ ~

Z
1 ~Pei~P·~r/h̄ ψ(~P)d3 P
= 3/2
(2.125)
(2πh̄)

donde se ha utilizado la condición de completez (2.114) y (2.122).

Ahora de (2.125),
−→
~Pei~P·~r/h̄ = −ih̄ ∇ ei~P·~r/h̄ (2.126)
y (2.125) se escribe como:
Z !
−→ 1 i~
P·~r/h̄
h~r |~P|ψi = −ih̄ ∇ e P ) d3 P
ψ(~
b
(2πh̄)3/2
Z 
−→
= −ih̄ ∇ h~r |~Pih~P|ψid3 P

−→
= −ih̄ ∇ h~r |ψi (2.127)

entonces ~
P en la representación posición es:
b
−→
~P
b
= −ih̄ ∇ (2.128)
En componentes cartesianas:
∂ ∂ ∂
Pbx = −ih̄ , Pby = −ih̄ , Pbz = −ih̄ (2.129)
∂x ∂y ∂z
ya se sabe que Pb es hermítico, i.e., Pb† = P.
b
b2
~P
El operador Hamiltoniano H b=
b se escribe como H +V~ en la representación posición:
2m

h̄2 2 b h̄2
 2
∂2 ∂2

∂
H=−
b ∇ + V (~r ) = − + 2 + 2 +V b (~r ), (2.130)
2m 2m ∂x2 ∂y ∂z
2.20 Mecánica ondulatoria 67
;A<

donde ∇2 es el operador Laplaciano.

En el espacio de los momentos el operador posición es:

b j = ih̄ ∂ ,
R ( j = x, y, z) (2.131)
∂Pj

b = ih̄ ∂ ,
X b = ih̄ ∂ ,
Y b = ih̄ ∂
Z (2.132)
∂Px ∂Py ∂Pz

2.20 Mecánica ondulatoria

Representando el formalismo de la mecánica cuántica en una base continua, se llega a un problema de autovalores, no en la
forma de la ecuación matricial de Heisenberg, sino en la forma de una ecuación diferencial. La representación de una ecuación
de autovalores como en (2.71) en el espacio da:

h~r | H
b |ψi = Eh~r |ψi (2.133)
se utiliza (2.130),
!
h̄2 2 ~
h~r | − ∇ + V (~r ) ψ = Eψ(~r )
2m

h̄2 2
− ∇ ψ(~r ) + V
b (~r )ψ(~r ) = Eψ(~r ) (2.134)
2m
Esta ecuación es conocida como la ecuación de Schrödinger, sus soluciones dan el espectro de energía del sistema como también
la función de onda. La formulación de la mecánica cuántica en la representación posición se llama mecánica ondulatoria.
A diferencia de Heisenberg, Schrödinger tomó un camino diferente al tratar de explicar teóricamente las ideas de Bohr. Él
comenzó con la hipótesis de la onda-partícula de De Broglie y las extendió a los electrones que orbitan el núcleo del átomo.

2.21 Postulados de la mecánica cuántica

Se analizarán los postulados sobre los cuales la comprensión cuántica de un sistema físico está basado. Estos postulados no se
derivan, sino que resultan del experimento.
Postulado 2.1. En un tiempo fijo t0 , el estado de un sistema físico se define al especificar un ket |ψ(t0 )i que pertenece al espacio de Hilbert
H.
b que actúa sobre H ; este operador es un observable.
Postulado 2.2. Toda cantidad física medible A se describe por operador A
Postulado 2.3. El único posible resultado de la medida de una cantidad física A es uno de los autovalores del correspondiente observable
A.
b

De este postulado se debe tener en cuenta que:

i Una medida de A siempre da una valor real ya que el operador asociado con A es hermítico.

ii Si el espectro de A
b es discreto, el resultado que se obtiene por medir A es cuantizado.

Postulado 2.4. Cuando la cantidad física A se mide sobre un sistema en el estado normalizado |ψi, la probabilidad Pn de obtener el
autovalor no degenerado an del correspondiente observable A
b es:

Pn = |hφn |ψi|5 , (2.135)

donde |φn i es un autovector normalizado de A asociado con el autovalor an .

5 Si |ψi es no normalizado, entonces: Pn = |hφn |ψi|2 /hψ|ψi
2.22 Operador evolución 68
;A<

Para autovalores an degenerados, la probabilidad se escribe como:

n
|hφn |ψi|2 ,
X
Pn = (2.136)
i =1

donde n es el número de degeneración para an y {|φn i} con n = 1, 2, 3, · · · , es un conjunto ortonormal de vectores los cuales
forman una base en el autovalor an de A.
b
Para el caso de un espectro no degenerado continuo, cuando la cantidad física A es medida sobre un sistema en el estado
normalizado |ψi, la probabilidad dPα de obtener un resultado incluido entre 2 y 2 + dα es igual a:

dPα = |hφα |ψi|2 dα , (2.137)

donde |φα i del observable A

b asociado con A.

Postulado 2.5. Si la medida de la cantidad física A sobre el sistema en el estado |ψi da el resultado an , el estado del sistema inmediatamente
después de la medida es la proyección normalizada
P |ψi
p n ,
hψ| Pn |ψi
de |ψi sobre el espacio asociado con an .
Postulado 2.6. La evolución temporal del vector estado |ψ(t)i es gobernada por la ecuación de Schrödinger:

d
ih̄ |ψ(t)i = H
b (t)|ψ(t)i, (2.138)
dt

donde H
b (t) es el observable asociado con la energía total del sistema.

2.22 Operador evolución

La transformación del vector estado |ψ(t0 )i en un tiempo inicial t0 a un vector |ψ(t)i en un tiempo arbitrario t es lineal.
b (t, t0 ) tal que:
Entonces existe un operador lineal U

|ψ(t)i = U
b (t, t0 )|ψ(t0 )i, (2.139)

U (t, t0 ) se define como el operador evolución del sistema.

2.22.1. Propiedades del operador evolución

i Ya que |ψ(t0 )i es arbitrario, de (2.139) se observa que:

b (t0 , t0 ) = 1,
U (2.140)

Si se sustituye (2.139) dentro de la ecuación de Schrödinger (2.138), se obtiene:

∂ b
ih̄ U (t, t0 ) = H (t)U
b (t, t0 ), (2.141)
∂t
Zt Zt
b (t, t0 ) = − i
∂U H ( t 0 )U
b (t0 , t0 )dt0 , (2.142)
h̄
t0 t0

t Zt
0 i
U ( t , t0 ) = − H ( t 0 )U
b (t0 , t0 )dt0 ,

b (2.143)
t0 h̄
t0

Zt
b (t, t0 ) = 1 − i
U H ( t 0 )U
b (t0 , t0 )dt0 , (2.144)
h̄
t0
2.22 Operador evolución 69
;A<

(2.139) se puede escribir en la forma:

|ψ(t)i = Ub (t, t0 )|ψ(t0 )i

(2.145)
|ψ(t0 )i = U
b (t0 , t00 )|ψ(t00 )i

Combinando las expresiones en (2.145):

b (t, t0 )U
|ψ(t)i = U b (t0 , t00 )|ψ(t00 )i, (2.146)
Además se puede escribir:
b (t, t00 )|ψ(t00 )i,
|ψ(t)i = U (2.147)
Comparando (2.146) y (2.147)
b (t, t00 ) = U
U b (t, t0 )U
b (t0 , t00 ), (2.148)

Generalizando este procedimiento:

U b ( t n , t n −1 ) · · · U
b ( t n , t1 ) = U b ( t 3 , t 2 )U
b ( t2 , t1 ) (2.149)

Si se asume que t1 < t2 < t3 < · · · < tn el sistema progresa de un tiempo menor a un tiempo mayor.
En (2.148) si t = t00 , se obtiene:
b (t, t0 )U
1=U b ( t 0 , t ), (2.150)
por tanto se obtiene que:
b (t0 , t) = U
U b −1 (t, t0 ) (2.151)

Ahora se escribe la ecuación de Schrödinger (2.138) en la forma:

d|ψ(t)i = |ψ(t + dt)i − |ψ(t)i, (2.152)

i b
=− H (t)|ψ(t)idt
 h̄ 
i b
|ψ(t + dt)i = 1 − H (t)dt |ψ(t)i, (2.153)
h̄
Se emplea (2.139), en donde:
|ψ(t + dt)i = U
b (t + dt, t)|ψ(t)i, (2.154)
y (2.153) se escribe como:
b (t + dt, t) = 1 − i Hdt,
U b (2.155)
h̄

este operador se llama, el operador evolución infinitesimal. Ya que H b (t) es hermítico, entonces U
b (t + dt, t) también lo es.
Y de (2.151) se concluye:
Ub (t0 , t) = U
b −1 (t, t0 ) = Ub † (t, t0 )

es un operador unitario.
Apéndices

70