Mecánica Cuántica 1

Oscilador armónico
Importancia:
• Brook Taylor (1685 – 1731)
• Joseph Fourier (1768 – 1830)

• Potencial: F= –kx V = ½ k x2 = ½ m w2 x2
𝑘
• Se conoce k, la masa m, y la frecuencia clásica 𝜔 =
𝑚
Oscilador armónico
• Ecuación de Schrödinger independiente del tiempo:
ℏ2 2 𝑚𝜔2 𝑥 2
−2𝑚 ∇ 𝜓 + 2
𝜓 = 𝐸𝜓
• Ecuación adimensional:
𝑑2 𝜓
− 2 + 𝑢2 𝜓 = 𝜀𝜓
𝑑𝑢
• Convenciones:
𝑚𝜔 2𝐸
𝑢≡𝑥 𝜀≡
ℏ ℏ𝜔
Solución analítica
𝑑2 𝜓 2 𝜓 =0
+ 𝜀 − 𝑢
𝑑𝑢2 Esto resulta de
buscar la
− 𝑢 2 Τ2
Definición: 𝜓 𝑢 = 𝑒 𝑓(𝑢) tendencia de y
para u grande.

𝑑2 𝑓 𝑑𝑓
2
− 2𝑢 + 𝜀−1 𝑓 =0
𝑑𝑢 𝑑𝑢

Serie de potencias: f(u) = S an un

2𝑛 − 𝜀 + 1 En general,
𝑎𝑛+2 = 𝑎𝑛 esto diverge!
(𝑛 + 1)(𝑛 + 2)
Cuantización de la energía
Si se cumple 2n – e +1 = 0 la serie se trunca
y la función exp(–u2/2) garantiza la convergencia.

Se necesita e = 2n + 1 → E = (n + ½)ℏw (n es entero ≥ 0)

Esto se parece a la propuesta de Planck

Resuelve el problema de la radiación térmica.
Pero la energía mínima no es 0.
Polinomios
E = (n + ½)ℏw (n es entero ≥ 0)
− 𝑢 2 Τ2
Cada estado n tiene función de onda 𝜓𝑛 𝑢 = 𝑒 𝑓𝑛 (𝑢)
Los fn(u) son polinomios con términos 0, 1, 2… n
Enumeramos los términos con el índice “i”
2𝑖 − 𝜀 + 1 e = 2n + 1
𝑎𝑖+2 = 𝑎𝑖 n≥0
(𝑖 + 1)(𝑖 + 2)
2(𝑖 − 𝑛)
𝑎𝑖+2 = 𝑎𝑖 i≤n
(𝑖 + 1)(𝑖 + 2)
Los primeros polinomios
a0 a1 a2 a3 a4 a5
n=0 1
Verificar!
n=1 1

n=2 1 –2 Tradición:
El coeficiente de un
n=3 1 –2/3 es 2n

n=4 1 –4 +4/3

n=5 1 –4/3 +4/15

Se llaman “Polinomios de Hermite”
a0 a1 a2 a3 a4 a5
n=0 1

n=1 2

n=2 –2 4

n=3 –12 8

n=4 12 –48 16

n=5 120 –160 32

Charles Hermite
Matemático francés
1822 – 1901
[Link]
[Link]
[Link]
 n
n
n
n

n
n
n

[Link]
Propiedades de los polinomios de Hermite
𝑑𝐻𝑛
= 2 𝑛 𝐻𝑛−1
𝑑𝑢

𝐻3 = 8𝑢3 − 12𝑢

𝑑𝐻3
= 24𝑢2 − 12 = 6 4𝑢2 − 2 = 2 ∙ 3 ∙ 𝐻2
𝑑𝑢
Propiedades de los polinomios de Hermite
+∞
−𝑢2
න 𝐻𝑖 𝐻𝑗 𝑒 𝑑𝑢 = 2𝑖 𝑖! 𝜋 𝛿𝑖𝑗
−∞
“Delta de Kronecker” dij = 1 si i=j dij = 0 si i≠j

1 𝑢2
−
𝜓𝑛 𝑢 = 𝑒 2 𝐻𝑛 (𝑢)
2𝑛 𝑛! 𝜋

න 𝜓𝑖 ∗ 𝜓𝑗 𝑑𝑢 = 𝛿𝑖𝑗

Las funciones de onda forman un conjunto ortonormal.

 n
n
n
n
n

n
n

[Link]
Lograr lo imposible
• Clásicamente, el oscilador solo llega hasta donde E = V = ½ mw2x2
es decir e = u2

• El límite clásico es 𝑢 = ± 𝜀 = ± 2𝑛 + 1

• Según la gráfica, el oscilador cuántico puede ir más allá!

 n
n
n
n
n

n
n

𝐿í𝑚𝑖𝑡𝑒 𝑐𝑙á𝑠𝑖𝑐𝑜: 𝑢 = ± 𝜀 = ± 2𝑛 + 1

[Link]
Probabilidad de lo clásicamente imposible
• Primer estado excitado: n=1 y = u exp(-u2/2)
• Límite clásico: u2=3
• Probabilidad:
− 3 ∞
‫׬‬−∞ 𝜌 𝑑𝑢 + ‫ ׬‬3 𝜌 𝑑𝑢
𝑃= 0 ∞
‫׬‬−∞ 𝜌 𝑑𝑢 + ‫׬‬0 𝜌 𝑑𝑢

∞ −𝑢2
‫ ׬‬3𝑒 𝑢2 𝑑𝑢 0,050
𝑃= ∞ −𝑢2 = = 0,113
‫׬‬0 𝑒 𝑢2 𝑑𝑢 0,443
Valores esperados: 𝑥
Por lógica, o por simetría, “esperamos” <u>=0

Demostración:
+∞
𝑢 =න 𝜓 2 𝑢 𝑑𝑢
−∞

Como |y|2 es simétrica (par), |y|2u es impar y su integral es nula.

Valores esperados: 2
<u > para n=0
+∞
1 2
𝑢2 = 𝜋
න 𝑒 −𝑢 𝑢2 𝑑𝑢
−∞ +∞
2
𝑢2 = −2 1𝜋 න 𝑢 𝑒 −𝑢 (−2)𝑢 𝑑𝑢
+∞ −∞
2 1
න −𝑢2
𝑢 = 2 𝜋
𝑒 𝑑𝑢
−∞
𝑢2 = 12

𝑥 2 = 12 ℏ
𝑚𝜔
Valores esperados: 2
<u > para n=1
+∞
2
𝑢2 = 1
2 𝜋
න 𝑒 −𝑢 𝐻1 2 𝑢2 𝑑𝑢
−∞ +∞
1 2
𝑢2 = 2 𝜋
න 𝑒 −𝑢 2𝑢 2 𝑢2 𝑑𝑢
+∞ −∞
2 −1
න 3 −𝑢2
𝑢 = 𝜋
𝑢 𝑒 (−2)𝑢 𝑑𝑢
−∞ +∞
3 2
𝑢2 = 𝜋
න 𝑢2 𝑒 −𝑢 𝑑𝑢
−∞
3 ℏ
𝑢2 = 3
2
𝑥2 =
2 𝑚𝜔
Dx
∆𝑥 = 𝑥2 − 𝑥 2

1 ℏ 1
Para n=0: ∆𝑥 = ∆𝑢 = = 0,707
2 𝑚𝜔 2

3 ℏ 3
Para n=1: ∆𝑥 = ∆𝑢 = = 1,225
2 𝑚𝜔 2

Δ𝑥 1
= 3 (guardan la misma proporción que los límites clásicos)
Δ𝑥 0
Valores esperados: 𝑝
Por lógica “esperamos” que <p> = 0 (porque el oscilador está fijo)

Demostración:
El operador “p” es – i ℏ ∂/∂x
Para comenzar, |y|2 es simétrica (par),
La derivada cambia en 1 el orden del polinomio
Luego y* ∂y/∂x es impar y su integral es nula.
 2
Valores esperados: 𝑝 para n=0
El operador “p” es – i ℏ ∂/∂x
𝑑 2 𝜓𝑛
𝑝2 = −ℏ2 න 𝜓𝑛 ∗ 2
𝑑𝑥
𝑑𝑥
1 𝑢2
−4 −2
𝜓0 𝑢 = 𝜋 𝑒 𝑑𝜓0 𝑑 2 𝜓0
= −𝑢 𝜓0 = (𝑢2 − 1)𝜓0
𝑑𝑢 𝑑𝑢2

𝑑 2 𝜓0
𝑝2 0 = −𝑚ℏω න 𝜓0 ∗ 𝑑𝑢 = 𝑚ℏω න 𝜓 0
∗
𝜓 0 1 − 𝑢 2 𝑑𝑢
𝑑𝑢2

𝑝2 0 = 𝑚ℏω 1 − 𝑢2 = 𝑚ℏω 1 − 𝑢2 = 𝑚ℏω 1 − 12 = 1

2 𝑚ℏω
 2
Valores esperados: 𝑝 para n=1
1 − 𝑢 2 Τ2
𝜓1 𝑢 = 𝑢 𝑒
2 𝜋

𝑑𝜓1 2 − 𝑢 2 Τ2 𝑑 2 𝜓1 3 − 𝑢 2 Τ2 𝑑 2 𝜓1
= 1−𝑢 𝑒 = (𝑢 − 3𝑢)𝑒 = (𝑢2 − 3)𝜓1
𝑑𝑢 𝑑𝑢2 𝑑𝑢2

𝑑2 𝜓1
𝑝2 1 = −𝑚ℏω න 𝜓1 ∗ 𝑑𝑢 = 𝑚ℏω න 𝜓 0
∗
𝜓 0 3 − 𝑢 2 𝑑𝑢
𝑑𝑢2

𝑝2 1 = 𝑚ℏω 3 − 𝑢2 = 𝑚ℏω 3 − 𝑢2 = 𝑚ℏω 3 − 32 = 3

2
𝑚ℏω
Dp
1
• Para n=0: Δ𝑝 = 𝑚ℏω
2

3
• Para n=1: Δ𝑝 = 𝑚ℏω
2

Productos:
Para n=0: Dp Dx = ℏ/2
Para n=1: Dp Dx = 3ℏ/2
Energías:
• Para n=0: <V> = ½ k <x2> = ½ k (½ ℏw/k) = ¼ ℏw = ½ E
<K> = <p2>/2m = (½ mℏw)/2m = ¼ ℏw = ½ E

• Para n=1: <V> = ½ k <x2> = ½ k (3/2 ℏw/k) = 3/4 ℏw = ½ E

<K> = <p2>/2m = (3/2 mℏw)/2m = 3/4 ℏw = ½ E

• Conclusión: <x2> = (n + ½) ℏ/mw

<p2> = (n + ½) ℏmw
Estados mixtos
• Las soluciones yn encontradas son estados estacionarios
• Tienen energías En = (n + ½) ℏw
෍ 𝑐𝑛 2 = 1
• Una combinación lineal y = S cn yn 𝑛
satisface la Ecuación de Schrödinger dependiente del tiempo
• Pero NO la ecuación independiente del tiempo (¿cuál E se pondría?)
• El valor esperado de la energía puede tener cualquier valor:
• Pero al medirla, saldrá uno de los En 2
𝐻 = ෍ 𝑐𝑛 𝐸𝑛
con probabilidad |cn|2 𝑛
Estados mixtos
• Una combinación lineal y = S cn yn
satisface la Ecuación de Schrödinger dependiente del tiempo
• Pero NO la ecuación independiente del tiempo (¿cuál E se pondría?)
• No son soluciones estacionarias

• Pero por la combinación de longitudes de onda

podrían ser muy localizadas (paquetes de ondas)
 [Link]
n
n
Ejemplo:
n
Una onda localizada cerca de u=+1
Se puede construir con
y1 + y6 – y3 – y4n

n
n

n
Estados mixtos
• Cualquier función y se puede expandir como y = S cn yn
• Los yn son una “base” ortonormal
• Por ejemplo una onda muy localizada en u ≈ 1
• Eso significa que el oscilador está cerca de u = 1
• <u> ≠ 0
• Entonces el potencial lo empujará hacia la izquierda
• El paquete de ondas se moverá hacia la izquierda
[Link]
En la imagen gif original se ve cómo el oscilador oscila!
Movimiento de una onda localizada
𝜓 𝑥, 𝑡 = Σ 𝑐𝑛 𝜓𝑛 𝑥 𝑒 −𝑖𝜔𝑛 𝑡
𝐸𝑛
𝜔𝑛 = = 𝑛 + 12 𝜔
ℏ

𝜓 𝑥, 𝑡 = 𝑒 −𝑖𝜔𝑡/2 ෍ 𝑐𝑛 𝜓𝑛 𝑥 𝑒 −𝑖𝑛𝜔𝑡
𝑛
Con la y construida calculamos la posición esperada: <x>

𝑥 = න 𝜓 ∗ 𝑥, 𝑡 𝑥 𝜓 𝑥, 𝑡 𝑑𝑥
Movimiento de una onda localizada
𝑥 = න 𝑥 𝜓 𝑥, 𝑡 2 𝑑𝑥

∗ ∗
𝑥 = ෍ ෍ 𝑐𝑝 𝑐𝑞 𝑒 𝑖 𝑝−𝑞 𝜔𝑡 න 𝜓𝑝 𝑥 𝜓𝑞 𝑑𝑥
𝑝 𝑞

Necesitamos los “elementos matriciales de x”: xpq = ∫ yp* x yq dx

Como los yn son ortogonales, solo hay elementos matriciales p = q±1

Elementos matriciales
𝑝 𝑝+1
𝑢𝑝𝑞 = 𝛿𝑞,𝑝−1 + 𝛿𝑞,𝑝+1
2 2

0 1 0 0 0

1 0 2 0 0
1 0 1 0 1 0 0 0 …
𝑢 = 1 0 2 = 2 3
2 2 0 0 0
0 2 0 3 4

0 0 0 4 0

… …
Movimiento de <x>
∗ 𝑖 𝑝−𝑞 𝜔𝑡
𝑝 𝑝+1
𝑢 = ෍ ෍ 𝑐𝑝 𝑐𝑞 𝑒 𝛿𝑞,𝑝−1 + 𝛿𝑞,𝑝+1
2 2
𝑝 𝑞

𝑛
𝑢 =෍ (𝑐𝑛 ∗ 𝑐𝑛−1 𝑒 𝑖𝜔𝑡 + 𝑐𝑛−1 ∗ 𝑐𝑛 𝑒 −𝑖𝜔𝑡 )
2
𝑛

𝑐𝑛 = 𝑐𝑛 𝑒 𝑖𝜙𝑛
𝑛
𝑢 = ෍2 𝑐𝑛 𝑐𝑛−1 cos(𝜔𝑡 + 𝜙𝑛−1 − 𝜙𝑛 )
2
𝑛