4a Serie de Ejercicios de clase - Teorı́a cuántica

Fı́sica moderna II

Ana Maria Alba Villarraga - 20181135059

Proyecto Curricular de Licenciatura en Fı́sica, Universidad Distrital Francisco José de Caldas

amalbav@[Link]

Julio 2022

Aplicaciones de la ecuación de Schrödinger unidimensional a diferentes sistemas. Establezca aproximaciones a

situaciones reales.

1. Considere una partı́cula de masa m que se mueve bajo la acción del potencial cero, esto es, V (x) = 0.
Explique la relación entre la energı́a del punto cero y el principio de incertidumbre.

Debido a que el potencial es nulo y que la fuerza se define como: − ∂V∂x(x) = 0, el momentum y la energı́a es
constante. Ası́ que se toma la ecuación de Schrödinger independiente del tiempo.
ℏ2 ∂ 2 ψ(x)
− + V (x)ψ(x) = Eψ(x) (1)
2m ∂x2
Teniendo en cuenta la condición V (x) = 0 y organizando la expresión anterior.
∂ 2 ψ(x) 2mE
+ 2 ψ(x) = 0 (2)
∂x2 ℏ
Ya que el termino que acompaña la función de onda es una constante, se define como: k 2 = 2mE/ℏ2 . Entonces:
∂ 2 ψ(x)
+ k 2 ψ(x) = 0 (3)
∂x2
La solución a la ecuación (3), tiene la forma de:
ψ(x) = A sin(kx) + B cos(kx) (4)
Se evalúan las condiciones de frontera, pues la función de potencial tiene valor V (x) = 0 en el intervalo 0 ≤ x ≤ L y
V (x) = ∞ en el intervalo x < 0 y x > L. Es decir, que en la zona para la cual se cumple la condición del enunciado:
V (x) = 0 en x = 0.
ψ(0) = A sin(k · 0) + B cos(k · 0) (5)

0=A·0+B·1 ⇒ B=0 (6)

V (x) = 0 en x = L.
ψ(L) = A sin(kL) + B cos(kL) (7)

0 = A sin(kL) + 0 (8)
 Según el resultado de la ecuación (6), B = 0, obteniendo que 0 = A sin(kL) pero el coeficiente A no puede ser
cero, de forma que:
A ̸= 0 ⇒ sin(kL) = 0 (9)
La función seno únicamente es nula para todo múltiplo entero de π. Lo que implica que:
√
2mE
kL = nπ → L = nπ (10)
ℏ
Al reemplazar la constante k, es posible despejar la energı́a de la partı́cula baja la acción del potencial cero.
n2 ℏ2 π 2
E= (11)
2mL2
El estado mı́nimo de energı́a se da respecto al mı́nimo valor de energı́a cinética, ya que la potencial es nula. Esto es:
2
mvmin p2
Emin = Kmin = = min (12)
2 2m
Igualando las expresiones de la energı́a (11) y (12), con n = 1.
ℏ2 π 2 p2min ℏ2 π 2
= → = p2min (13)
2mL2 2m L2
De manera que si ℏ = h/2π, el momentum mı́nimo de la partı́cula está dado por:
ℏπ h
pmin = = (14)
L 2L
Finalmente, este ultimo resultado se puede expresar como una incertidumbre del momentum de la siguiente manera:
h
∆p ∼ (15)
2L
Comparando el resultado con el principio de incertidumbre de Heisenberg, la distancia que dispone la partı́cula para
desplazarse es L. Que corresponde a la incertidumbre de la posición ∆x = L. Se puede inferir que aunque ∆x sea
muy pequeño, ∆x se hace demasiado grande y por ello la partı́cula no está en reposo.

2. Aplique la ecuación de Schrödinger para cada una de las diferentes regiones descritas por el potencial
escalón descrito en la figura adjunta en las que se mueve una partı́cula de masa m.

i) Para el caso del potencial escalón V0 > E, pruebe que la suma de los coeficientes de reflexión y de transmisión
es igual a uno.

Para el potencial escalón que muestra la figura, se soluciona la ecuación de Schrödinger independiente del tiempo
(ecuación 1) en los segmentos en los que está definida la función potencial.

2
• V (x) = 0 en x < 0.
ℏ2 ∂ 2 ψ(x) ∂ 2 ψ(x) 2mE
− 2
= Eψ(x) → + 2 ψ(x) = 0 (16)
2m ∂x ∂x2 ℏ
2 2
Se define la constante de la ecuación como: k1 = 2mE/ℏ .
∂ 2 ψ(x)
+ k12 ψ(x) = 0 (17)
∂x2
La solución de la ecuación se expresa de la siguiente forma.
ψ1 (x) = Aeik1 x + Be−ik1 x (18)

∂ψ(x)
= ik1 Aeik1 x − ik1 Be−ik1 x (19)
∂x

• V (x) = V0 en x > 0.
ℏ2 ∂ 2 ψ(x) ℏ2 ∂ 2 ψ(x)
− + V0 ψ(x) = Eψ(x) → − ψ(x) [V0 − E] = 0 (20)
2m ∂x2 2m ∂x2

∂ 2 ψ(x) 2m
− 2 ψ(x) [V0 − E] = 0 (21)
∂x2 ℏ
Se define la constante de la ecuación como: k2 = (V0 − E)2m/ℏ2 .
2

∂ 2 ψ(x)
− k22 ψ(x) = 0 (22)
∂x2
La solución de la ecuación se expresa de la siguiente forma.
ψ2 (x) = Cek2 x + De−k2 x (23)

∂ψ(x)
= k2 Cek2 x − k2 De−k2 x (24)
∂x

Para que el sistema de ecuaciones se satisfaga, se debe cumplir que la función de onda y su deriva sean finitas,
continuas y univaluadas.

1. Para que sean finitas x → ∞. En la ecuación (23) implica que:

lı́m Cek2 x = ∞ ⇒ C=0 (25)
x→∞

Lo que significa que la exponencial creciente no es aceptable, ya que involucra una probabilidad infinita de
encontrar a la partı́cula en el infinito. Razón por la cual el coeficiente de esa componente debe ser cero.
En consecuencia se tiene:
ψ(x) = De−k2 x (26)

∂ψ1 (x) ∂ψ2 (x)

2. Para que sean continuas |ψ1 (x)|x=0 = |ψ2 (x)|x=0 y ∂x = ∂x .
x=0 x=0

A eik1 x x=0
+ B e−ik1 x x=0
= D e−k2 x x=0
⇒ A+B =D (27)

ik2
ik1 A eik1 x x=0
− ik1 B e−ik1 x x=0
= −k2 D e−k2 x x=0
⇒ A−B =D (28)
k1

Si de (27) se toma B = D − A y se reemplaza en (28), se obtiene:

3
  
D ik2
A= 1+ (29)
2 k1
Por otra parte, si se toma de (27) A = D − B y se considera (28), se conoce:
 
D ik2
B= 1− (30)
2 k1
Se procede a calcular el coeficiente de reflexión con los complejos conjugados de las constantes de las ecuaciones
(29) y (30), definido como:
   
D ik2 D ik2
∗
B B 2 1 + k1 2 1 − k1
R= ∗ → R=     =1 (31)
A A D
1− ik2 D
1+ ik2
2 k1 2 k1

Este resultado indica que cuando las partı́culas inciden con energı́a total menor que la altura del escalón de
potencial, la totalidad de ellas son reflejadas o en otras palabras, rebotan en la región derecha del origen. Por
lo tanto, ningún porcentaje de partı́culas es transmitido, y al realizar la suma del coeficiente de reflexión y de
transmisión se obtiene.
R = 1, T = 0 ⇒ R+T =1+0=1 (32)

ii) Demuestre para el caso en el que qV0 < E, los números de onda k1 y k2 de la onda inicial y final respectivamente
están relacionados por: k2 /k1 = 1 − VE0 , y que los coeficientes de reflexión R y transmisión T están dados


 2
respectivamente por: R = kk11 −k 2
+k2 y T = (k4k+k
1 k2
)2
. Haga el cálculo para E = 1,2V0 y E = 2V0 .
1 2

Evaluando las condiciones de la función de potencial de acuerdo al esquema de la figura. Al considerar la energı́a
total de las partı́culas como la suma de la energı́a cinética y la energı́a potencial.

• V (x) = 0 en x < 0.
mv12 p2 √
E =K +0= = 1 → p1 = 2mE (33)
2 2m

• V (x) = V0 en x > 0.
mv22 p2 p
E = K + V0 = + V0 = 2 + V0 → p2 = 2m(E − V0 ) (34)
2 2m

La ecuación diferencial para x > 0, por la condición de las energı́as, cambia de (21) a lo siguiente.
∂ 2 ψ(x) 2m
+ 2 ψ(x) [E − V0 ] = 0 (35)
∂x2 ℏ

∂ 2 ψ(x)
+ k22 ψ(x) = 0 (36)
∂x2

ψ2 (x) = Ceik2 x + De−ik2 x (37)

∂ψ(x)
= ik2 Ceik2 x − ik2 De−ik2 x (38)
∂x

Pero para x < 0 se mantiene la solución a la ecuación diferencial de la expresión (18). De modo que las constantes
quedan definidas como: k12 = 2mE/ℏ2 y k22 = (E − V0 )2m/ℏ2 . Luego, utilizamos las ecuaciones (33) y (34)
para reescribir las costantes.

4
 √
2mE p1 mv1
k1 = = = (39)
ℏ ℏ ℏ
p
(E − V0 )2m p2 mv2
k2 = = = (40)
ℏ ℏ ℏ
A causa de que las partı́culas inciden sobre el escalón de potencial desde la izquierda, solo se pueden reflejar o
transmitir, no existe ninguna variación de potencial capaz de producir un rebote de partı́culas. Por ello el coeficiente
de reflexión de la región derecha D debe ser nulo. Entonces se debe reescribir la expresión (37) como sigue.
ψ2 (x) = Ceik2 x (41)
Nuevamente, la función de potencial como su derivada deben ser finitas, continuas y univaluadas, para solucionar
el sistema de ecuaciones.
Para que sean continuas |ψ1 (x)|x=0 = |ψ2 (x)|x=0 y ∂ψ∂x1 (x)
= ∂ψ∂x 2 (x)
.
x=0 x=0
ik1 x −ik1 x ik2 x
A e x=0
+B e x=0
=C e x=0
⇒ A+B =C (42)

ik1 A eik1 x x=0

− ik1 B e−ik1 x x=0
= ik2 C e−ik2 x x=0
⇒ k1 (A − B) = k2 C (43)
Reemplazando la expresión (42) en la (43).
k1 (A − B) = k2 (A + B) → k1 A − k1 B = k2 A + k2 B (44)

k1 − k2
(k1 − k2 )A = (k2 + k1 )B ⇒ B= A (45)
k1 + k2
Ahora para encontrar la expresión del coeficiente C, se toma la ecuación (42) y la (45).
k1 − k2 2k1
A+ A=C ⇒ C= A (46)
k1 + k2 k1 + k2
A continuación, se define el coeficiente de reflexión como la fracción entre el número de partı́culas que rebotan y
el número de partı́culas incidentes, relacionando los coeficientes de la ecuación (45).
∗  2
B∗B
  
k1 − k2 k1 − k2 k1 − k2
R= ∗ → R= = (47)
A A k1 + k2 k1 + k2 k1 + k2
El coeficiente de transmisión se expresa como la fracción entre el número de partı́culas que superan el escalón de
potencial y el número de partı́culas incidentes. Solo que considerando la velocidad de la partı́cula cuando incide
y cuando se transmite.
∗  2
v2 C ∗ C
  
v2 2k1 2k1 v2 2k1
T = → T = = (48)
v 1 A∗ A v1 k1 + k2 k1 + k2 v1 k1 + k2
Tomando las relaciones entre los números de onda y las velocidades de las ecuaciones (39) y (40), se reemplaza
de la siguiente manera:
k2 ℏ 

2
m 2k1 k2 4k12 4k1 k2
T = k1 ℏ
= 2 = 2 (49)
m
k1 + k2 k1 (k1 + k2 ) (k1 + k2 )
Sumando los coeficientes, definidos finalmente en (47) y (49).
2 2
k12 − 2k1 k2 + k22 + 4k1 k2

k1 − k2 4k1 k2 (k1 + k2 )
R+T = + 2 = 2 = 2 =1 (50)
k1 + k2 (k1 + k2 ) (k1 + k2 ) (k1 + k2 )
El cociente entre el número de onda de la onda final k2 y el de la onda inicial k1 , de acuerdo a como se definieron
estas dos constantes en las ecuaciones diferenciales, es:

5
 √ 
(E−V0 )2m r r
p
k2 ℏ (E − V0 )2m (E − V0 )2m V0
= √  = √ = = 1− (51)
k1 2mE 2mE 2mE E
ℏ

Ası́ que los coeficientes de reflexión y de transmisión se expresan en términos de la ecuación (51) como se muestra:
!2 !2  q 2

k1 − k2
2 k1
k1 − k2
k1 1 − k2
k1
1 − 1 − VE0
R= = k1 k2
= =  (52)
1 + kk21
q
k1 + k2 k1 + k1 1 + 1 − V0 E

q
k2
k2 4k12 k2 4 k12 k2 4 1 − VE0
4
1
T = = 2 = 2 =  (53)
k1 (k1 + k2 )2
2
k1 k1 k1
  q
k1 + kk12 1 + kk12 1 + 1 − VE0

Cálculo para E = 1, 2V0

 q 2  q 2
V0 1
1− 1− 1,2V0
1− 1− 1,2

1 − 0, 408
2
R= q  = q  = = 0, 176 (54)
1+ 1− V0
1+ 1− 1 1 + 0, 408
1,2V0 1,2

q q
V0 1
4 1 − 1,2V0
4 1 − 1,2 4(0, 408)
T = q 2 =  q 2 = 2 = 0, 823 (55)
V0
1 + 1 − 1,2V0 1 + 1 − 1,21 (1 + 0, 408)

R + T = 0, 176 + 0, 823 = 0, 999 (56)

Cálculo para E = 2V0
 q 2  q 2
V0 1
1− 1− 2V0 1− 1− 2

1 − 0, 707
2
R= q  = q  = = 0, 029 (57)
1+ 1− V0
1+ 1− 1 1 + 0, 707
2V0 2

q q
V0
4 1 − 2V0
4 1 − 21 4(0, 707)
T = q 2 = q 2 = 2 = 0, 970 (58)
(1 + 0, 707)

V0 1
1 + 1 − 2V 0
1 + 1 − 2

R + T = 0, 029 + 0, 970 = 0, 999 (59)

3. Resuelva la ecuación de Schrödinger para una partı́cula de masa m, que se mueve en una caja de
potencial de lados a, b, c.

En la ecuación de Schrödinger independiente del tiempo, se toma la función de onda como una variable de x, y, z.
ψ(r) = X(x)Y (y)Z(z) (60)
Dentro de la caja, la función de potencial es nula V (x, y, z) = 0. Entonces:
ℏ2 ∂ 2 ψ(x) ∂ 2 ψ(y) ∂ 2 ψ(z)
 
− + + = Eψ(r) (61)
2m ∂x2 ∂y 2 ∂z 2
Se establecen las respectivas derivadas.

6
 ∂ψ(x) ∂X(x) ∂ 2 ψ(x) ∂ 2 X(x)
= Y (y)Z(z) ; 2
= Y (y)Z(z) (62)
∂x ∂x ∂x ∂x2

∂ψ(y) ∂Y (y) ∂ 2 ψ(y) ∂ 2 Y (y)

= X(x)Z(z) ; = X(x)Z(z) (63)
∂x ∂y ∂y 2 ∂y 2

∂ψ(z) ∂Z(z) ∂ 2 ψ(z) ∂ 2 Z(z)

= Y (y)X(x) ; = X(x)Y (y) (64)
∂x ∂z ∂z 2 ∂z 2
Reescribiendo la expresión (61), con las segundas derivadas parciales de las ecuaciones (62), (63) y (64), y con (60).
ℏ2 ∂ 2 X(x) ∂ 2 Y (y) ∂ 2 Z(z)
 
− Y (y)Z(z) + X(x)Z(z) + X(x)Y (y) = E [X(x)Y (y)Z(z)] (65)
2m ∂x2 ∂y 2 ∂z 2

1 ∂ 2 X(x) 1 ∂ 2 Y (y) 1 ∂ 2 Z(z) 2m

2
+ 2
+ =− 2 E (66)
X(x) ∂x Y (y) ∂y Z(z) ∂z 2 ℏ
Puesto que la energı́a es constante, y por ende todo el termino derecho lo es, cada termino del lado izquierdo de la
expresión anterior se iguala a una constante que corresponde al vector de onda correspondiente.
1 ∂ 2 X(x) 1 ∂ 2 Y (y) 1 ∂ 2 Z(z)
= −kx2 ; = −ky2 ; = −kz2 (67)
X(x) ∂x2 Y (y) ∂y 2 Z(z) ∂z 2
Entonces la expresión (66) se modifica.
ℏ2
k 2 + ky2 + kz2


E= (68)
2m x
De acuerdo a las expresiones de la ecuación (67), se tienen tres ecuaciones con derivadas parciales, cuyas soluciones
son de la forma:
∂ 2 X(x)
+ kx2 X(x) = 0 ⇒ X(x) = Ax sin(kx x) + Bx cos(kx x) (69)
∂x2

∂ 2 Y (y)
+ ky2 Y (y) = 0 ⇒ Y (y) = Ay sin(ky y) + By cos(ky y) (70)
∂y 2

∂ 2 Z(z)
+ kz2 Z(z) = 0 ⇒ Z(z) = Az sin(kz z) + Bz cos(kz z) (71)
∂z 2
Se procede a evaluar las condiciones de frontera para el movimiento de la partı́cula en la caja respecto a cada
componente, dado por las ecuaciones (69), (70) y (71), cuando su potencial vale cero en el intervalo 0 ≤ x ≤ a.
Para x = 0
X(0) = 0 ⇒ 0 = Ax sin(kx · 0) + Bx cos(kx · 0) (72)

0 = Ax · 0 + B x · 1 ⇒ Bx = 0 (73)
De forma análoga para las componentes y y z:
Y (0) = 0 ⇒ By = 0 (74)

Z(0) = 0 ⇒ Bz = 0 (75)
Para x = 0
X(a) = 0 ⇒ 0 = Ax sin(kx · a) + Bx cos(kx · a) → 0 = Ax sin(kx · a) + 0 (76)
El coeficiente Ax no puede ser nulo, por este motivo, la función seno es la que queda igualada a cero, y esta solo
puede tener este resultado para múltiplos enteros de π.

7
 nx π
sin(kx · a) = 0 ⇒ kx a = nx π → kx = (77)
a
Igualmente, de forma análoga para y y z.
ny π nz π
ky = ; kz = (78)
a a
Donde nx , ny y nz toman valores de 1,2,3,4,...
Se reformula la función de onda (60) con las ecuaciones (69), (70) y (71) al considerar los resultados de (73), (74)
y (75), además de los vectores de onda de las expresiones (77) y (78).
ψ(r) = Ax sin(kx x)Ay sin(ky y)Az sin(kz z) (79)
n π  n π  n π 
x y z
ψ(x, y, z) = Ax Ay Az sin x sin y sin z (80)
a a a
Finalmente también se reformula la energı́a de la partı́cula (68) con los vectores de onda (77) y (78).
ℏ2 π 2 2
nx + n2y + n2z


E= 2
(81)
2m a

4. Resuelva la ecuación de Schrödinger para cada una de las diferentes regiones descritas por el potencial de
barrera descrito en la figura adjunta, en las que se mueve una partı́cula de masa m.

i) Escriba las soluciones generales para las distintas regiones del potencial.

Resolviendo la ecuación de Schrödinger en las tres regiones.

Región 1 y región 3 cuando V (x) = 0

ℏ2 ∂ 2 ψ(x) ∂ 2 ψ(x) 2m
− 2
= Eψ(x) → + 2 Eψ(x) = 0 (82)
2m ∂x ∂x2 ℏ
El termino que acompaña la función de onda es constante, por ello se define como k1 que corresponde al número
de onda, y se expresa de la siguiente manera: k12 = 2mE
ℏ2 .

∂ 2 ψ(x)
+ k12 ψ(x) = 0 (83)
∂x2
Para x < 0 la solución a la ecuación tiene la siguiente forma:
ψ1 (x) = Aeik1 x + Be−ik1 x (84)

Para a < x la solución a la ecuación tiene la siguiente forma:

ψ3 (x) = F eik1 x + Ge−ik1 x (85)

8
 Dado que las partı́culas inciden desde la izquierda sobre la barra de potencial, solo se permite la transmisión en
la región 3, y por ende el segundo termino de la ecuación (85) es nulo.
ψ3 (x) = F eik1 x (86)

∂ψ3 (x)
= ik1 F eik1 x (87)
∂x
Región 2 cuando V (x) = V
Se pueden presentar dos casos, el primero es cuando al energı́a de la partı́cula es menor que la altura de la barrera
de potencial de manera que las partı́culas chocan con la pared, y el otro es cuando la energı́a es mayor que el
potencial y la partı́cula logra superar la barrera, moviéndose por arriba sin chocar.

Para E < V
ℏ2 ∂ 2 ψ(x) ∂ 2 ψ(x) 2m
− + V ψ(x) = Eψ(x) → = 2 (−E + V )ψ(x) (88)
2m ∂x2 ∂x2 ℏ

∂ 2 ψ(x) 2m ∂ 2 ψ(x) 2
− 2 (V − E)ψ(x) = 0 → − k2a ψ(x) = 0 (89)
∂x2 ℏ ∂x2
2
Con el vector de onda como: k2a = 2m(Vℏ2−E)
Debido a la forma de la ecuación (89), su solución es:
ψ2a (x) = Ca ek2a x + Da e−k2a x (90)

∂ψ2a (x)
= k2a Ca ek2a x − k2a Da e−k2a x (91)
∂x

Para E > V
∂ 2 ψ(x) 2m ∂ 2 ψ(x) 2
+ 2 (E − V )ψ(x) = 0 → + k2b ψ(x) = 0 (92)
∂x2 ℏ ∂x2
Con el vector de onda como: k2b 2
= 2m(E−V
ℏ2
)

Cuya solución tiene la siguiente forma:

ψ2b (x) = Cb eik2b x + Db e−ik2b x (93)

∂ψ2b (x)
= ik2b Cb eik2b x − ik2b Db e−ik2b x (94)
∂x

ii) Determine el coeficiente de transmisión T y según sea el caso establezca cuando hay efecto túnel.

Desarrollando el caso en que E < V . Las funciones de onda y las derivadas de la región 1 y la región 2 deben
ser continuas en x = 0, se debe cumplir la condición |ψ1 (x)|x=0 = |ψ2a (x)|x=0 y ∂ψ∂x
1 (x)
= ∂ψ∂x 2a (x)
.
x=0 x=0
ik1 x −ik1 x k2a x −k2a x
A e x=0
+B e x=0
= Ca e x=0
+ Da e x=0
⇒ A+B =C +D (95)

ik1
ik1 A eik1 x x=0
− ik1 B e−ik1 x x=0
= k2a Ca ek2a x x=0
− k2a Da e−k2a x x=0
⇒ (A − B) = Ca − Da
k2a
(96)
Ahora, de igual manera, la funciones de onda y las derivadas de la región 2 y la región 3 deben ser continuas en
x = a, se debe cumplir la condición |ψ2a (x)|x=a = |ψ3 (x)|x=a y ∂ψ∂x
2a (x)
= ∂ψ∂x 3 (x)
.
x=a x=a

9
 Ca ek2a x x=a
+ Da e−k2a x x=a
= F eik1 x x=a
⇒ Ca ek2a a + Da e−k2a a = F eik1 a (97)

ik1 ik1 a
k2a Ca ek2a x x=a
− k2a Da e−k2a x x=a
= ik1 F eik1 x x=a
⇒ Ca ek2a a − Da e−k2a a = Fe (98)
k2a
Al sumar las expresiones (95) y (96).
ik1 ik1 A ik1 B
Ca + Da + Ca − Da = A + B + (A − B) = A + B + − (99)
k2a k2a k2a
   
ik1 ik1
2Ca = A 1 + +B 1− (100)
k2a k2a
Y al restar las mismas ecuaciones.
ik1 ik1 A ik1 B
Ca + Da − Ca + Da = A + B − (A − B) = A + B − + (101)
k2a k2a k2a
   
ik1 ik1
2Da = A 1 − +B 1+ (102)
k2a k2a
Por otra parte, se suman las ecuaciones (97) y (98).
ik1 ik1 a
Ca ek2a a + Da e−k2a a + Ca ek2a a − Da e−k2a a = F eik1 a + Fe (103)
k2a
   
ik1 ik1
2Ca e k2a a
= Fe ik1 a
1+ → 2Ca = F eik1 a
1+ e−k2a a (104)
k2a k2a
También, se restan.
ik1 ik1 a
Ca ek2a a + Da e−k2a a − Ca ek2a a + Da e−k2a a = F eik1 a − Fe (105)
k2a
   
−k2a a ik1 a ik1 ik1 a ik1
2Da e = Fe 1− → 2Da = F e 1− ek2a a (106)
k2a k2a
A partir de las ecuaciones, se iguala la expresión (100) con la (104), y la (102) con la (106).
     
ik1 ik1 ik1
A 1+ +B 1− = F eik1 a 1 + e−k2a a (107)
k2a k2a k2a
     
ik1 ik1 ik1
A 1− +B 1+ = F eik1 a 1 − ek2a a (108)
k2a k2a k2a
Despejando el coeficiente B de ambas ecuaciones.
   
1 + kik2a1 1+ ik1
k2a k2a + ik1 ik1 a −k2a a k2a + ik1
B=  F eik1 a e−k2a a −  A = Fe e − A (109)
1 − kik2a1 1− ik1 k2a − ik1 k2a − ik1
k2a

   
ik1 ik1
1− k2a 1− k2a k2a − ik1 ik1 a k2a a k2a − ik1
B=  F eik1 a ek2a a −  A = Fe e − A (110)
1+ ik1
1+ ik1 k2a + ik1 k2a + ik1
k2a k2a

Por lo que permite igualar las expresiones resultantes y organizar.

   
ik1 a k2a + ik1 −k2a a k2a − ik1 k2a a k2a + ik1 k2a − ik1
Fe e − e =A − (111)
k2a − ik1 k2a + ik1 k2a − ik1 k2a + ik1

10
 Operando los términos entre los corchetes de ambos lados de la igualdad.
 2 2
2ik2a k1 −k2a a k12 − k2a
2

k − k2a −k2a a 2ik2a k1 k2a a 4ik2a k1
F eik1 a 12 2 e + 2 2 e − 2 2 e k2a a
+ 2 2 e =A 2 2 (112)
k1 + k2a k1 + k2a k1 + k2a k1 + k2a k1 + k2a
Se aplican las siguientes funciones hiperbólicas en relación a las funciones exponenciales.
ex + e−x ex − e−x
sinh(x) = ; cosh(x) = (113)
2 2
 2 2

k − k2a 2ik2a k1 −ik1 a 4ik2a k1
F 12 2 (−2 sinh(k2a a)) + 2 2 (2 cosh(k2a a)) = Ae (114)
k1 + k2a k1 + k2a 2
k12 + k2a
Ası́ se obtiene la expresión final del coeficiente F .
Ae−ik1 a 4ik2a k1 2ik2a k1 Ae−ik1 a
F = 2 2 = (115)
−2 sinh(k2a a)(k1 − k2a ) + 4ik1 ka2 cosh(k2a a) 2ik1 ka2 cosh(k2a a) − (k12 − k2a
2 ) sinh(k a)
2a

El cuadrado de su módulo es:

2 2 2 2 2 2
2 4k2a k1 |A| 4k2a k1 |A|
|F | = 2 2 2 2 2 2 2 = 2 2 2 2 )2 sinh2 (k a)
(116)
4k1 ka2 cosh (k2a a) + (k1 − k2a ) sinh (k2a a) 4k1 ka2 + (k1 + k2a 2a

Se denomina coeficiente de transmisión a la proporción de partı́culas incidentes que son transmitidas, entonces se
halla como el cociente entre el cuadrado del modulo de la componente transmitida y la incidente, contemplando
la velocidad de las partı́culas en cada caso.
2
2 4k2a k12 |A|2
2 2
Jt v |F | 4k12 ka2
2 +(k 2 +k 2 )2 sinh2 (k
1 2a 2a a) 4k2a k1
T = = 2 = 2 = 2 )2 sinh2 (k a)
(117)
Ji v |A| |A| 4k12 ka2
2 + (k12 + k2a 2a

Reemplazando los números de onda k1 y k2a respecto a la energı́a y la altura de la barrera.

  
4 2m(Vℏ2−E) 2mE)
ℏ2 4E(V − E)
T =    √ = √ 
2m(V −E) 2m(V −E)
4 2m(Vℏ2−E) 2mE)
ℏ2 + ( 2m(V −E)
ℏ2 + 2mE) 2
ℏ2 ) sinh 2
ℏ a 4E(V − E) + V 2 sinh2
ℏ a
(118)
El efecto tunel se da a través de una barrera cuya altura V0 excede la energı́a total E, con el comportamiento
ondulatorio de la partı́cula, cuando esta atraviesa una barrera que clásicamente no podrı́a superar, es decir, si
E < V o existe cierta probabilidad de que la partı́cula sea transmitida a través de la barrera a la región x > a.
Pero en la región lejana a la barrera se puede detectar una partı́cula tanto antes como después de que penetre la
región sin introducir una incertidumbre significativa en el conocimiento de su energı́a.

iii) Evalúe el coeficiente de transmisión para un electrón de energı́a 2eV que incide sobre una barrera de potencial
rectangular de altura 4eV y espesor 10−10 m.

Se hace uso de la ecuación (118) para hallar el coeficiente de transmisión, donde E = 2eV, V = 4eV y a =
10−10 m. Además, ℏ = h/2π = 4, 136 × 10−15 eV·s/2π = 6, 582 × 10−16 eV·s, y m = 5,1 × 105 eV/c2 .
4(2eV)(4eV − 2eV)
T = √  (119)
2 2 2(5,1×105 eV/c2 )(4eV−2eV) −10
4(2eV)(4eV − 2eV) + (4eV) sinh 6,582×10−16 eV·s 10 m

16eV2
T = = 0, 617 (120)
16eV2 + 16eV2 sinh2 (0, 723)

11
5. Resuelva la ecuación de Schrödinger para cada una de las diferentes regiones descritas por el potencial de
pozo cuadrado finito descrito en la figura adjunta, en las que se mueve una partı́cula de masa m.

Se evalúa la ecuación de Schrödinger independiente del tiempo en cada una de las tres regiones.

En las regiones I, cuando V (x) = 0

ℏ2 ∂ 2 ψ(x) ∂ 2 ψ(x) 2m
− 2
= Eψ(x) → + 2 Eψ(x) = 0 (121)
2m ∂x ∂x2 ℏ
Puesto que la energı́a E es negativa y la constante de la ecuación diferencial se define como una raı́z que contiene
este termino, k1 como número de onda debe ser positivo y por ello se define como: k12 = − 2mEℏ2 . Reorganizando la
expresión anterior.
∂ 2 ψ(x)
− k12 ψ(x) = 0 (122)
∂x2
Para x < 0 la solución a la ecuación tiene la siguiente forma:
ψ1 (x) = Aek1 x + Be−k1 x (123)

Para x > a la solución a la ecuación tiene la siguiente forma:

ψ3 (x) = F e−k1 x + Gek1 x (124)
En vista de que la función de onda debe ser finita, se estudia el comportamiento de la variable en infinito x → ∞
en cada una de las funciones. En el caso de la ecuación (123), x siempre toma valores negativos, es decir, que el
primer término se convierte en una exponencial decreciente.
lı́m Ae−k1 x = 0 ; lı́m Bek1 x = ∞ (125)
x→∞ x→∞

Lo que implica que el segundo termino no es admisible, ya que hace crecer la función de onda sin lı́mite y por tanto,
el coeficiente debe ser nulo, B = 0. Entonces ψ1 (x) queda como:
ψ1 (x) = Aek1 x (126)

∂ψ1 (x)
= k1 Ae−k1 x (127)
∂x
Por otro lado, considerando los valores que toma x en la ecuación (124), esta expresión se mantiene y al valorar el
lı́mite se obtiene:
lı́m F e−k1 x = 0 ; lı́m Gek1 x = ∞ (128)
x→∞ x→∞

Ası́ mismo, como en el caso anterior, el coeficiente G debe ser cero. Considerando este resultado la función de onda
ψ3 (x) es:
ψ3 (x) = F e−k1 x (129)

12
 ∂ψ3 (x)
= −k1 F e−k1 x (130)
∂x

En la región II, cuando V (x) = −V0

ℏ2 ∂ 2 ψ(x) ∂ 2 ψ(x) 2m
− − V 0 ψ(x) = Eψ(x) → = 2 (−E − V0 )ψ(x) (131)
2m ∂x2 ∂x2 ℏ
Recordando que tanto la energı́a como el potencial son constantes y negativos, pero que E sigue siendo mayor que
V0 . Se organiza la expresión y la constante queda determinada como: k22 = 2m(E+V
ℏ2
0)
.
∂ 2 ψ(x)
+ k22 ψ(x) = 0 (132)
∂x2
Para 0 ≤ x ≤ a la solución a la ecuación tiene la siguiente forma:
ψ2 (x) = C sin(k2 x) + D cos(k2 x) (133)

∂ψ2 (x)
= k2 C cos(k2 x) − k2 D sin(k2 x) (134)
∂x
Ahora, se evalúa la continuidad de la función de onda y su derivada. Primero se realiza el procedimiento para las
soluciones pares que se obtienen usando la función coseno de la función de onda dentro del pozo. Relacionando la
ecuación (133) con C = 0 y la ecuación (129).

∂ψ2 (x) ∂ψ3 (x)

Entre ψ2 (x) y ψ3 (x) en x = a, con |ψ2 (x)|x=a = |ψ3 (x)|x=a y ∂x = ∂x .
x=a x=a

D |cos(k2 x)|x=a = F e−k1 x x=a

⇒ D cos(k2 a) = F e−k1 a (135)

−k2 D |sin(k2 x)|x=a = −k1 F e−k1 x x=a

⇒ −k2 D sin(k2 a) = −k1 F e−k1 a (136)
Dividiendo la expresión (136) entre la (135).
−k2 D sin(k2 a) −k1 F e−k1 a
= ⇒ k2 tan(k2 a) = k1 (137)
D cos(k2 a) F e−k1 a
Esta es una ecuación trascendental que no tiene solución analı́tica, y es con la que se obtienen las energı́as permisibles.
Por lo que se necesitan de métodos numéricos o gráficos para resolverla. Sus soluciones se conocen como los valores
propios.
Se adopta la notación siguiente para llegar a una expresión más simplificada.
ap
z = k2 a ; z0 = 2mV0 (138)
ℏ
Si se tienen en cuenta los números de onda respecto a la energı́a y el potencial, se halla la siguiente relación.
−2mE 2m(E + V0 ) 2mV0
q
k12 + k22 = + = → k 1 a = z02 − z 2 (139)
ℏ2 ℏ2 ℏ2
Se escribe la expresión (137) como la igualación de dos funciones respecto a z. Los puntos en los que se intersectan
tales funciones corresponderán a las soluciones de valores propios de la ecuación trascendental.
r
k1 a 1 z02
q
tan(z) = = k1 = 2
z0 − z 2 ⇒ tan(z) = −1 (140)
k2 z z z2
Por ultimo, se evalúa la continuidad de la función de onda y su derivada pero ahora para las soluciones impares, con
la función seno de la ecuación (133) donde D = 0 y la ecuación (129).
C |sin(k2 x)|x=a = F e−k1 x x=a
⇒ C sin(k2 a) = F e−k1 a (141)

13
 k2 C |cos(k2 x)|x=a = −k1 F e−k1 x x=a
⇒ k2 C cos(k2 a) = −k1 F e−k1 a (142)
Dividiendo la expresión (142) entre la (141).
k2 C cos(k2 a) −k1 F e−k1 a
= ⇒ k2 cot(k2 a) = −k1 (143)
C sin(k2 a) F e−k1 a
Al comparar este resultado con el de la ecuación (137) que da lugar a la expresión (140), la expresión resulta.
r
k1 z02
− cot(k2 a) = ⇒ − cot(z) = −1 (144)
k2 z2
Si se le asigna un valor a la anchura del pozo a y a su profundidad V0 . La gráfica muestra las soluciones pares dadas
por la ecuación (140) como las soluciones impares por la ecuación (144), para un pozo de potencial que admite
ciertos estados discretos de energı́a.

6. Verifique la función propia y el valor propio para el estado n = 2 de un oscilador armónico simple por
substitución directa en la ecuación de Schrödinger independiente del tiempo.

La ecuación de Schrödinger unidimensional e independiente del tiempo es:

ℏ2 ∂ 2 ψ(x)
− + V (x)ψ(x) = Eψ(x) (145)
2m ∂x2
Debido a que la energı́a potencial de un oscilador armónico simple corresponde a V (x) = kx2 /2, con k como la
constante elástica, se contempla este valor en la ecuación anterior.
ℏ2 ∂ 2 ψ(x) k 2
− + x ψ(x) = Eψ(x) (146)
2m ∂x2 2
q
k
Si la frecuencia angular del oscilador es ω = m y ası́ mismo se define clásicamente como ω = 2πν, la frecuencia
q
1 k
del sistema termina siendo ν = 2π m y en ese sentido, la constante elástica se reemplaza por k = 4ν 2 π 2 m.

ℏ2 ∂ 2 ψ(x) ∂ 2 ψ(x) 2ν 2 π 2 m(2m) 2 2m

− 2ν 2 π 2 mx2 ψ(x) + Eψ(x) = 0 → − x ψ(x) + 2 Eψ(x) = 0
2m ∂x2 ∂x2 ℏ2 ℏ
(147)

14
 " 2 #
∂ 2 ψ(x)

2m 2νπm
+ E− x2 ψ(x) = 0 (148)
∂x2 ℏ2 ℏ

Se definen dos constantes β = 2mE/ℏ2 y α = 2νπm/ℏ. Y la expresión anterior se expresa:

∂ 2 ψ(x) 
+ β − α2 x2 ψ(x) = 0

2
(149)
∂x
r
√ q
2νπm 2πm
q
k (km)1/4
Ahora se sustituye y = αx, es decir, y = ℏ x = ℏ2π mx =
√
ℏ
x Efectuando las derivadas
correspondientes de y por regla de la cadena y sustituyendo esto en la ecuación (149).
√
dψ dy dψ d( αx) dψ √ dψ d2 ψ d2 ψ
= = = α ; 2
=α 2 (150)
dx dx dy dx dy dy dx dy

d2 ψ  d2 ψ
 
β
+ β − αy 2 ψ = 0 2

α 2
→ + − y ψ=0 (151)
dy dy 2 α
Si se estudia el comportamiento asintótico de la ecuación, cuando y → ∞, el termino β/α se hace despreciable
frente a y 2 , ası́ que la ecuación diferencial se reduce a:
d2 ψ
− y2 ψ = 0 (152)
dy 2
Cuyas solución posee la forma de:
2 2
ψ(y) = Aey /2
+ Be−y /2
(153)
Ahora, la función de onda debe ser finita, y esto se considera al evaluar lo siguiente:
2 2
lı́m Aey /2
=∞ ; lı́m Be−y /2
=0 (154)
z→∞ z→∞

Lo cual siguiere que el coeficiente A sea nulo, pues el termino al que corresponde hace crecer sin lı́mite a la función.
Permitiendo escribir la función de onda como:
2
ψ(y) = Be−y /2
(155)
Para que la función de onda mantenga la caı́da asintótica se debe multiplicar por un polinomio finito. Es decir, un
conjunto de funciones que permiten la expresión correcta de ψ(y).
2
ψ(y) = Be−y /2
H(y) (156)
Se procede a realizar los cálculos de la primera y segunda derivada.
dψ(y) 2 2 dH(y)
= −Bye−y /2 H(y) + Be−y /2 (157)
dy dy

d2 ψ(y) −y 2 /2 2 −y 2 /2 −y 2 /2 dH(y)
2
−y 2 /2 d H(y)
= −Be H(y) + By e H(y) − Bye + Be
dy 2 dy dy 2

d2 ψ(y) dH(y) d2 H(y)

 
−y 2 /2 2
= Be −H(y) + y H(y) − 2y + (158)
dy 2 dy dy 2
Con el fin de sustituir estos resultados en la ecuación (151).
dH(y) d2 H(y)
 
2 β 2 2
Be−y /2 −H(y) + y 2 H(y) − 2y + + Be−y /2 H(y) − y 2 Be−y /2 H(y) = 0 (159)
dy dy 2 α

d2 H(y)
 
dH(y) β
− 2y + − 1 H(y) = 0 (160)
dy 2 dy α

15
Se supone una solución polinómica, en la forma de un polinomio finito.
∞
X
H(y) = al y l = a0 + a1 y + a2 y 2 + a3 y 3 + .... (161)
l=0

Efectuando las respectivas derivadas.

∞
dH(y) X
= lal y l−1 = a1 + 2a2 y + 3a3 y 2 + ... (162)
dy
l=1

∞
d2 H(y) X
= l(l − 1)al y l−2 = 2a2 + 6a3 y + 12a4 y 2 + ... (163)
dy 2
l=2

Como consecuencia la ecuación (160) se convierte en:

 
β
2a2 + 6a3 y + 12a4 y 2 + 20a5 y 3 + ... − 2y(a1 + 2a2 y + 3a3 y 2 + ...) + − 1 (a0 + a1 y + a2 y 2 + a3 y 3 + ....) = 0
α
(164)
Separando los términos con diferente orden de y.
 
β
y0 → 2a2 + − 1 a0 = 0 (165)
α
   
β β
y1 → 6a3 − 2a1 + − 1 a1 = 2 · 3a3 + − 1 − 2 a1 = 0 (166)
α α
   
β β
y2 → 12a4 − 4a2 + − 1 a2 = 3 · 4a4 + − 1 − 2 · 2 a2 = 0 (167)
α α
   
β β
y3 → 20a5 − 6a3 + − 1 a3 = 4 · 5a5 + − 1 − 2 · 3 a3 = 0 (168)
α α
Lo que es suficiente para notar una estructura definida en las expresiones que se dan para cada orden de y, encontrando
una fórmula de recurrencia.
 
β
α − 1 − 2l
al+2 = − al (169)
(l + 1)(l + 2)
Despreciando los términos más pequeños.
 
β
al+2 α − 1 − 2l
2l 2
=− ≈ 2 = (170)
al (l + 1)(l + 2) l l
Para no tener un polinomio con infinitos valores, se establece solo permitir valores enteros positivos, delimitando las
soluciones donde el numerador frente a esta condición se hace cero.
β β
l = n n = 0, 1, 2, 3, ... ⇒ − 1 − 2n = 0 → = 1 + 2n (171)
α α
Y sustituyendo las expresiones de β y α ya definidas.


2mE/ℏ2 E
= = 1 + 2n (172)
(2νπm/ℏ) πℏν
   
hν 1
E = 1 + 2n (πℏν) = 1 + 2n → E= + n hν (173)
2 2

16
Este indica que la energı́a de un oscilador armónico simple no pueden tomar cualquier valor dentro de un rango
continuo, sino que está cuantizada.
Verificando la función propia y el valor propio para el estado n = 2.
 
1 5
E= + 2 hν = hν (174)
2 2
Para el caso, la función de onda de la ecuación (156) es:
2
ψ2 (y) = Be−y /2
(4y 2 − 2) (175)
Con sus derivadas, empleando la regla de la cadena.
d (km)1/4  (km)1/4 
 
dψ dy dψ d  −y2 /2 2 2 2

= = √ x Be (4y − 2) = √ B −4y 3 e−y /2 + 10ye−y /2 (176)
dx dx dy dx ℏ dy ℏ

d2 ψ dy d2 ψ (km)1/4 (km)1/4 
   
d 2 2

2
= = √ √ B −4y 3 e−y /2 + 10ye−y /2
dx dx dydx ℏ dy ℏ

d2 ψ (km)1/2
   
2 2 2
= B 4y 4 e−y /2 − 22y 2 e−y /2 + 10e−y /2 (177)
dx2 ℏ
El resultado de la ecuación (177), (174), (175) y la sustitución de y directamente se sustituye en la ecuación de
Schrödinger independiente del tiempo, expresada para el oscilador armónico simple por (146).
ℏ2 (km)1/2  k  ℏy 2 
  
4 −y 2 /2 2 −y 2 /2 −y 2 /2 2 5 2
− B 4y e − 22y e + 10e + √ Be−y /2 (4y 2 − 2) = hνBe−y /2 (4y 2 − 2)
2m ℏ 2 km 2
(178)
√ r
ℏ2 km 4 2 kℏ 2 2

5 k
4y 2 − 2



− 4y − 22y + 10 + √ y 4y − 2 = ℏ (179)
2m 2 km 2 m

−4y 4 + 22y 2 − 10 + 4y 4 − 2y 2 = 20y 2 − 10 ⇒ 0=0 (180)

7. Determinar los valores ⟨x⟩, ⟨p⟩, x2 y p2 en los estados ψ0 y ψ1 del oscilador armónico. Comprobar
que se cumple la relación de incertidumbre.

Los estados ψ0 y ψ1 del oscilador armónico, se definen con la siguiente ecuación, en los estados n = 0 y n = 1
para el polinomio finito H(y).
1  mω 1/4 −y2 /2
ψn (y) = √ e Hn (y) (181)
2n n! πℏ
1/4
q
Se hace uso de la expresión y = (km)
√
ℏ
x y a la relación ω = k
m.

1  mω 1/4 2
 mω 1/4 √
km/2ℏx2
ψ0 (y) = √ e−y /2
H0 (y) = e− (182)
20 0! πℏ πℏ
 mω 1/4 2
ψ0 (x) = e−mω/2ℏx (183)
πℏ

1  mω 1/4 2 1  mω 1/4 −√km/2ℏx2  mω 1/4 √ 2 2 (km)

1/4
ψ1 (y) = √ e−y /2
H1 (y) = √ e (2y) = e− km/2ℏx √ √ x
21 1! πℏ 2 πℏ πℏ 2 ℏ
(184)

17
 r
 mω 1/4 √ 2mω
− km/2ℏx2
ψ1 (x) = e x (185)
πℏ ℏ
mω 1/4

p mω
Definiendo dos constantes como: α = πℏ yβ= ℏ x. Convirtiendo las expresiones a:
2
ψ0 (x) = αe−β /2
(186)
√ 2
ψ1 (x) = 2αβe−β /2
(187)
Los valores esperados para la posición y el momento de la partı́cula cuando la función de onda es real, son:
Z ∞
2
⟨x⟩n = x |ψn (x)| dx (188)
−∞
Z ∞
2
x2 n
= x2 |ψn (x)| dx (189)
−∞

Z ∞  
ℏ d
⟨p⟩n = |ψn (x)| ψn (x)dx (190)
−∞ i dx
∞
d2
Z  
2 2
p = |ψn (x)| −ℏ ψn (x)dx (191)
n
−∞ dx2
Para n = 0
Z ∞ !r
∞ ∞ ∞
Z Z r Z
2 ℏ ℏ
2 −β 2 /2 −β 2 2
⟨x⟩0 = x |ψ0 (x)| dx = x αe dx = α 2
xe dx = α 2
βe−β dβ
−∞ −∞ −∞ mω mω −∞

  " −β 2 #∞  
2 ℏ e 2 ℏ
=α − =α [0] = 0 (192)
mω 2 mω
−∞

∞ ∞ ∞ ∞
r
β 2 ℏ −β 2
Z Z 2
Z Z
2 2 2 2 −β 2 /2 2 2 −β 2 2 ℏ
x 0
= x |ψ0 (x)| dx = x αe dx = α x e dx = α e dβ
−∞ −∞ −∞ −∞ mω mω
r Z ∞
r √  r r √ 
ℏ ℏ 2 ℏ ℏ π mω ℏ ℏ π ℏ
= α2 β 2 e−β dβ = α2 = = (193)
mω mω −∞ mω mω 2 πℏ mω mω 2 2mω
Z ∞   Z ∞  
ℏ d 2 ℏ d 2
⟨p⟩0 = |ψ0 (x)| ψ0 (x)dx = αe−β /2
αe−β /2
dx
−∞ i dx −∞ i dx
" 2
#∞
∞
ℏ 2 ∞ −β 2
r
e−β
Z  r Z
2ℏ −β 2 /2 mω −β 2 /2 ℏ ℏ 2 ℏ
=α e − βe dβ = − α βe dβ = − α − = − α2 [0] = 0
i −∞ ℏ mω i −∞ i 2 i
−∞
(194)
∞ ∞
d2 d2
Z   Z  
−β 2 /2 2
p2
= |ψ0 (x)| −ℏ 2
ψ0 (x)dx = αe −ℏ 2
αe−β /2
dx
0
−∞ dx2 −∞ dx2

∞  r
mω −β 2 /2 β 2 ℏ
Z  
2 2 −β 2 /2 mω −β 2 /2 ℏ
= −α ℏ e − e − e dβ
−∞ ℏ ℏ mω mω

18
 r ∞ √ 
√
r Z 
mω 2 mω ℏ h
−β 2 2 −β 2
i ℏmω π ℏmω
= ℏ e −β e dβ = √ π− = (195)
πℏ ℏ mω −∞ π 2 2

Para n = 1
Z ∞ !Z
∞ ∞ ∞
r r
2
Z √ −β 2 /2
2
2
Z
2 −β 2 2 ℏ 3 −β 2 ℏ
⟨x⟩1 = x |ψ1 (x)| dx = x 2αβe dx = 2α xβ e dx = 2α β e dβ
−∞ −∞ −∞ mω −∞ mω

" 2
#∞
β 2 + 1 e−β

 
2 ℏ ℏ
= 2α − = α2 [0] = 0 (196)
mω 2 mω
−∞

∞ ∞ ∞ ∞
r
2
Z
2 2
Z
2
√ −β 2 /2
2
2
Z
2 2 −β 2 2
Z
β 2 ℏ 2 −β 2 ℏ
x 1
= x |ψ1 (x)| dx = x 2αβe dx = 2α x β e dx = 2α β e dβ
−∞ −∞ −∞ −∞ mω mω
r Z ∞
r r Z ∞  √ 
ℏ 2 ℏ 4 −β 2 mω ℏ ℏ 4 −β 2 2 ℏ 3 π 3ℏ
= 2α β e dβ = 2 β e dβ = √ = (197)
mω mω −∞ πℏ mω mω −∞ π mω 4 2mω
Z ∞   Z ∞ √
 
√
ℏ d −β 2 /2 ℏ d 2
⟨p⟩1 = |ψ1 (x)| ψ1 (x)dx = 2αβe 2αβe−β /2
dx
−∞ i dx −∞ i dx
Z ∞ r r Z ∞h
−β 2 /2 mω  −β 2 /2 mω −β 2 /2 2  ℏ 2 mω −β 2 3 i
= 2α 2
βe e − e β dβ = 2α2 βe−β − e β dβ
−∞ ℏ ℏ mω −∞ ℏ
" 2 2
#∞
e−β mω (β 2 + 1)e−β
= 2α2 − + = 2α2 [0] = 0 (198)
2 ℏ 2
−∞

Z ∞  2
 Z ∞ √  2

√
2 d −β 2 /2 2 d 2
p 2
1
= |ψ1 (x)| −ℏ 2
ψ1 (x)dx = 2αβe −ℏ 2
2αβe−β /2 dx
−∞ dx −∞ dx
"r 3/2 !# r
∞
r
m2 ω 2 −β 2 /2 3
Z 
2 2 −β 2 /2 mω ℏ 3mω ℏ −β 2 /2 ℏ
= −2α ℏ βe e β − β e dβ
−∞ ℏ ℏ2 mω ℏ mω mω

∞
mω 2 mω ∞ h 4 −β 2
Z r r  r r Z
mω 4 −β 2 mω 2 −β 2 2
i
= −2α ℏ 2 2
β e −3 β e dβ = −2 ℏ β e − 3β 2 e−β dβ
−∞ ℏ ℏ πℏ ℏ −∞
 √ √   
2ℏmω 3 π 3 π 3 3ℏmω
=− √ − = −2ℏmω − =
π 4 2 4 2

19