Mecánica Cuántica

12/05/2020

Parte I
Solución de la ecuación radial para partículas
en potenciales centrales
1. La ecuación radial
Recordemos de las clases anteriores que quedó pendiente la solución radial R(r) de la ecuación
de Schrödinger para potenciales V (r) con simetría esférica
✓ ◆
@ 2 @R 2µr2
r [V (r) E] R = R. (1)
@r @r ~2
Del análisis de la parte angular del problema se demostró que = l(l + 1), por lo que la ecuación (1)
después de realizar la simplificación (@/@r ! d/dr) se convierte en
✓ ◆
d 2 dR 2µr2
r [V (r) E] R =l(l + 1)R;
dr dr ~2
 ✓ ◆
~2 d 2 dR ~2
r + l(l + 1) + V (r) R(r) =ER(r). (2)
2µr2 dr dr 2µr2

El caso más simple que podemos estudiar es el de una partícula libre en coordenadas esféricas.

2. La partícula libre (V (r) = 0) en coordenadas esféricas

El caso de una partícula libre corresponde a V (r) = 0, por lo que la ecuación radial (2) se reduce
a:
 ✓ ◆
~2 d 2 dR ~2 2µ
r + l(l + 1) R(r) =ER(r); multiplicamos por ,
2µr2 dr dr 2µr2 ~2
 ✓ ◆
1 d dR l(l + 1) 2µE
r2 + R(r) = 2 R(r); (3)
r2 dr dr r2 | ~{z }
⌘k2

2µ
si además definimos k 2 = ~2 E, lo que implica que la energía es

~2 k 2
E= ,
2µ
y la ecuación diferencial (3), después de desarrollar el término d
dr dr , se convierte en
r2 dR

d2 2 d l (l + 1)
+ R (r) = k 2 R (r) . (4)
dr2 r dr r2

1
Considere a continuación el cambio de variable x = kr, con los siguientes operadores diferenciales
d d dx d
= =k ,
dr dx dr dx
y
d2 2 d
2
= k .
dr2 dx2
La ecuación (4) se reduce a:

d2 R(r) 2 dR(r) l (l + 1)
+ + 1 R(r) = 0, (5)
dx2 x dx r2
la cual podemos identificar con ayuda de la siguiente tabla.
-Ecuación esférica de Bessel


d2 f (x) 2 df (x) l (l + 1)
+ + 1 f (x) = 0, (6)
dx2 x dx x2

-Solución general de la ecuación esférica de Bessel

Las soluciones f (x) son combinaciones lineales de funciones esféricas de Bessel jl (x) y funciones
esféricas de Neumann nl (x)

f (x) = A jl (x) + B nl (x),

-Funciones esféricas de Bessel jl (x)
l=0 1 2
3.14 4.49 5.76
Las funciones jl (x) son regulares en x = 0, esto Ceros de jl (x) :
es, son finitas en el origen. Algunas funciones 6.28 7.73 9.10
esféricas de Bessel : 9.42 10.90 12.32

senx Gráficas de jl (x):

j0 (x) = ; (7)
x
senx cosx
j1 (x) = 2 ;
✓x x
◆
3 1 3
j1 (x) = senx cosx.
x3 x x2

Ortogonalidad:
ˆ 1
⇡
jl (kr) jl0 (k 0 r) r2 dr = 2 (k k0 ) .
0 2k
(8)
-Funciones esféricas de Neumann nl (x)
Las funciones nl (x) son singulares en x = 0, es- Gráficas de nl (x):
to es, divergen en el origen. Algunas funciones
esféricas de Neumann:
cosx
n0 (x) = ;
x
cosx senx
n1 (x) = 2
;
✓x ◆x
3 1 3
n1 (x) = cosx senx.
x3 x x2

De la tabla anterior vemos que (5) es la ecuación esférica de Bessel, por lo que las solución general
es:
R (x) = C1 jl (x) + C2 nl (x) .

2
Además, como estamos considerando el origen (r = 0) y sabemos que nl (x) es singular en ese punto,
hacemos C2 = 0 en la solución anterior para solo quedarno con jl (x) (omitiendo la constante C1 ),

R (x) = jl (x) .

Por lo tanto la solución para el caso libre en tres dimensiones en coordenadas esféricas es:

klm = klm (r, ✓, ) = R (r) Ylm (✓, ) = jl (kr) Ylm (✓, ) .

La ortogonalidad de la klm libre en coordenadas esféricas se sigue de las ortogonalidad de las

funciones esféricas de Bessel jl (kr) dada por la ecuación (8), de tal forma que :

˚ ˆ 1 ˆ 2⇡ ˆ ⇡
0
⇤ 0 2
klm (r, ✓, ) k0 ,l0 ,m0 (r, ✓, ) dV = jl (kr) j (k r) r dr
l0 Ylm⇤ (✓, ) Yl”m (✓, ) sen (✓) d✓ d ;
0
V |0 0
{z }
m,m0 l,l0
˚ ˆ 1
⇤
klm (r, ✓, ) k0 ,l0 ,m0 (r, ✓, ) dV = m,m0 l,l0 jl (kr) jl0 (k 0 r) r2 dr;
0
V
⇡
˚
⇤
klm (r, ✓, ) k0 ,l0 ,m0 (r, ✓, ) dV = (k k0 ) m,m0 l,l0 .
2k 2
V

3. El punto cuántico esférico

El potencial para una partícula de masa µ en un punto cuántico esférico (o caja esférica impene-
trable) esta dado por (
0, r  a;
V (r) =
+1, r > a,
que satisface la ecuación radial (5). En este caso el origen (r = 0) esta incluido, por lo que utilizamos
la solución regular en la región r  a de la esfera

R (r) = jl (kr) .
La condición de que la solución se anule en la frontera r = a nos conduce a

R (a) =jl (ka) = 0;

r !
2µE
jl a =0. (9)
~2

Observe que:

La ecuación (9) es una ecuación trascendental. Para cada valor de l tendremos un número
infinito de valores de E.
La energía está cuantizada, pues solo ciertos valores (los ceros de jl ) están permitidos.
Para el caso particular de los estados s (l = 0), podemos calcular analíticamente la energía ya que
podemos utilizar la ecuación (7) y expresar a (9) de la siguiente forma:
r !
2µE
J0 a =0;
~2
✓q ◆
2µ
sen ~ 2 E a
✓q ◆ =0,
2µ
~2 E a

3
de donde vemos que debe cumplirse la condición
r
2µE
a = n⇡, n = 1, 2, · · · ,
~2
que nos conduce a la energía del pozo cuántico esférico

~2 n2 ⇡ 2
E = En0 = ; n = 1, 2, · · · ,
2µ a2
la expresión anterior coincide con el valor de la energía de una caja de paredes infinitas en una
dimensión.

Parte II
Soluciones para potenciales centrales
(V (r) 6= 0) y las propiedades de la solución
radial
Partimos de la ecuación (2)
2 3
6 ~2 d ✓ d ◆ ~2 7
6 7
6 2
r2 + 2
l (l + 1) + V (r)7 R (r) = ER (r) . (10)
4 2µr dr dr 2µr 5
| {z }
término centrífugo

Considere el cambio de variable u = rR en la ecuación anterior, y después de algunos pasos algebraicos

llegamos a
 2 2
~ d ~2
+ l (l + 1) + V (r) u (r) = Eu (r) . (11)
2µ dr2 2µr2
A continuación estudiaremos las propiedades generales de la solución radial u(r) mediante un análisis
del comportamiento asintótico de dicha función, sin especificar el valor de V (r).
La idea es conocer las propiedades generales de u (r) en los limites r ! 0 y r ! 1. La únicas
condiciones que impondremos al potencial son las siguientes:

lı́m r2 V (r) = 0, (12)

r!0
y
lı́m V (r) = 0. (13)
r!1

Esto es, estamos pidiendo que V (r) sea menos singular1 que 1/r2 (que es la misma dependencia en
r que tiene el término centrífugo 2 en (10) ). Dicho en otras palabras, pedimos que V (r) no se vaya
más rápido a cero que 1/r2 .

3.1. Comportamiento de u(r) cuando r ! 0

Vemos que al multiplicar la ecuación (11) por r2
1 Los potenciales de interés como el coulombiano (V (r) ⇠ r 1 ) o el potencial de Yukawa (V (r) ⇠ e µ0 r g r 1 , con
s
µ0 = mc/~) caen en esta categoría.
2 Se le llama término centrífugo pues dicho término actúa como un potencial V
cen = ~ l(l + 1)/2µr , que produce
2 2

una fuerza centrífuga.

4
 2 3
6 ~ 2 d
2 2
@
r@~ 2 2
2⇠⇠⇠ : 07 2 ⇠⇠ :0
⇠
4 r 2
+ l (l + 1) + ⇠
r V (r) 5 u (r) = ⇠
r ⇠
Eu (r), (14)
2µ dr r@
2µ@ 2 | {z }
ecuación (12)

y tomando el límite r ! 0, resulta que la ecuación diferencial que gobierna este régimen es

~2 2 d2 u (r) ~2
r 2
+ l (l + 1) u (r) =0;
2µ dr 2µ
d2 u (r) l (l + 1)
u (r) =0. (15)
dr2 r2
Proponemos una solución de la forma u(r) = r↵ y procedemos a calcular el valor de ↵ en el siguiente
ejercicio.

Ejercicio 3.1. Determine el valor del parámetro ↵ de la solución u = r↵ para la ecuación (15).

Demostración.
d2 r ↵ l (l + 1) ↵
r =0;
dr2 r2
l (l + 1) ↵
↵ (↵ 1) r↵ 2
r =0;
r2
↵ (↵ 1) l (l + 1) =0;
↵2 ↵ l (l + 1) =0;
[↵ + l] [↵ (l + 1)] =0,

y por tanto
(
l + 1;
↵=
l,

que conduce a dos soluciones

(
rl+1 , regular en r 0;
u (r) =
r l, singular en r = 0.
⌅
Del resultado del ejercicio anterior, la solución físicamente aceptable es la solución regular, por lo
que concluimos que en el límite r ! 0

u (r) ⇠ rl+1 . (16)

3.2. Comportamiento de u(r) cuando r ! 1

Cuando r es muy grande, la ecuación (11) se simplifica con ayuda de (13)
" #
~2 d 2 ~2 *0
+ l (l + 1) + V (r) u(r) = Eu(r)
2µ dr2 2µr2
y llegamos a
~2 d2 u(r)
= Eu(r).
2µ dr2

5
Si consideramos energías negativas3 (E < 0) y definimos la energía de amarre W = E > 0, la
ecuación anterior se expresa como
d2 u(r) 2µ
= 2 W u(r),
dr2 ~
y las soluciones físicamente aceptables cuando r ! 1 son de la forma:
q
2µ
W r
u(r) ⇠e ~2 ;
kr
u(r) ⇠e , (17)

en donde definimos r
2µ
k= W.
~2

Una vez conocido el comportamiento de las funciones u(r) en los límites r ! 0 y r ! 1, podemos
abordar problemas con potenciales V (r) específicos. Un problema típico es el del átomo hidrógeno
que estudiaremos en la siguiente sección.
3 En el caso de energías positiva (E > 0) tendríamos una situación de «colisión» (tipo dispersión de Rutherford)

que nos conduce a otro tipo de soluciones.