Capítulo 5.

El oscilador armónico

5.1. El oscilador armónico unidimensional

5.1.1. El cambio de escala
5.1.2. La solución en series
5.1.3. Los valores propios
5.1.4. La normalización
5.1.5. Los elementos de matriz
5.2. Los operadores de creación y de aniquilación
5.2.1. La ecuación de valores propios
5.2.2. La función propia del estado basal
5.3. Problemas

( )
Anexo 5.1. El comportamiento asintótico de u′′ + λ − ξ 2 u = 0

Anexo 5.2. Los polinomios de Hermite, Hn ξ

( )
5. El oscilador armónico
El estudio de las oscilaciones tiene gran importancia en la mecánica cuántica. Las
oscilaciones armónicas se usan como modelo para describir las interacciones que
presentan una posición de equilibrio y, en particular, son el modelo más usado para
estudiar vibraciones.

5.1. El oscilador armónico unidimensional

Un oscilador armónico se caracteriza por un fuerza que es proporcional al
desplazamiento, respecto a la posición de equilibrio. Por lo tanto, para este sistema se tiene

()
un potencial cuadrático, V x = 12 kx . En este caso, la ecuación de valores propios de la
2

energía, ĤuE = EuE , toma la forma

2
!2 d uE 1
− + mω 2 x 2uE = EuE , (5.1)
2m dx 2
2

en donde ω 2 ≡ k m , así,

m2ω 2 2mE
uE′′ − 2 x 2uE = − 2 uE . (5.2)
! !

5.1.1. El cambio de escala

Esta ecuación diferencial de segundo orden con coeficientes polinomiales presenta
varias combinaciones de constantes que pueden eliminarse mediante un cambio de escala

(reescalamiento). Sea ξ ≡ α x la nueva variable independiente, con α > 0 , entonces

d dξ d
dx
=
dx dξ
=α
d
dξ
() ( ) ()
y uE x = uE ξ α = u ξ . Mediante esta transformación la ecuación

diferencial toma la forma

2
⎛ mω ⎞
u′′ − ⎜ 2 ⎟ ξ 2u = − λ u , (5.3)
⎝ !α ⎠

( !α ) > 0 . Si se elige α de tal forma que mω ( !α ) = 1 , entonces

2
en donde λ ≡ 2mE 2

5-2
 u′′ − ξ u = − λ u ,
2
(5.4)

con α = mω ! y E = 12 !ωλ .

5.1.2. La solución en series

Para resolver la ecuación diferencial con coeficientes variables, ecuación (5.4), se usa el
método de Frobenius (método generalizado de serie de potencias). En este método se
propone una solución en series de la forma

∞ ∞

()
u ξ =ξ s
∑C ξ j
j
= ∑ C jξ s+ j , (5.5)
j=0 j=0

en donde C0 ≠ 0 . Al sustituir el serie en la ecuación diferencial se tiene que

( ) ( ) ( ( )( ) )
0 = C0 s s −1 ξ −2 + C1 s +1 sξ −1 + C2 s + 2 s +1 + C0λ
, (5.6)
+ (C ( s + 3)( s + 2) + C λ )ξ + ∑ ⎡⎣C ( s + j + 2)( s + j +1) − C
∞

3 1 j+2 j−2
+ C j λ ⎤⎦ ξ j
j=2

y al igualar los polinomios se obtienen las relaciones siguientes,

(
C s s −1 = 0 ,
0
) (
C s s +1 = 0 ,
1
)
( )(
C s +1 s + 2 + λC0 = 0 ,
2
) ( )(
C s + 2 s + 3 + C1λ = 0 ,
3
) (5.7)

( )(
C s + j +1 s + j + 2 − C j−2 + λC j = 0 ,
j+2
)
( j = 2,3,…) .
A la primera ecuación se le denomina ecuación indicial y permite determinar los
valores de s , en este caso se tiene que s = 0,1 . Debido a que se tiene un potencial
simétrico, las funciones propias deben tener paridad definida, pero, por la forma de la
serie de potencias, el factor que aparece fuera de la suma tiene la misma paridad que el

parámetro s , por lo que la suma siempre debe ser una función par. Por esta razón, C1 = 0 ,

al igual que todos los coeficientes impares. Los coeficientes pares se obtienen a partir de
las relaciones de recurrencia,

λ C0 C j−2 − λC j
C2 = −
(
s + 2 s +1
, C j+2 =
)( )
s + j + 2 s + j +1
,
( )( ) ( j = 2,4,6,…) .

(5.8)

5-3
 Esta ecuación diferencial tiene una relación recursiva de tres términos y no es posible
obtener una solución cerrada de los coeficientes. Por esta razón es necesario transformar la
ecuación diferencial.

Una posibilidad consiste en separar su comportamiento asintótico (Anexo 5.1). Para

2
− αξ
valores grandes de x , u′′ ~ ξ 2u , por lo que la solución debe tener la forma u ~ e . Así,

2
− αξ
u′ ∼ −2αξ e y

0 = u′′ − ξ 2u + λu ∼ e −αξ ⎡⎣4α 2ξ 2 − ξ 2 − 2α + λ ⎤⎦ ~ e −αξ ⎡⎣4α 2 −1⎤⎦ ξ 2 ,

2 2
(5.9)

por lo que 4α 2 = 1 , o bien α = ± 12 . Sólo para α negativa se tiene una solución

normalizable, por tanto u ~ exp − ξ 2 2 .

( )
El resultado anterior sugiere el cambio de variable

()
u ξ = e −ξ 2H ξ .

2

() (5.10)

2
2 2
−ξ 2 2 2 2 2
Así, u′ = −e −ξ 2ξ H + e −ξ 2H ′ , u′′ = −e H + e −ξ 2ξ 2H − 2e −ξ 2
ξ H ′ + e −ξ 2H ′′ y

(
0 = u′′ + λ − ξ 2 u = e −ξ

) 2
2
⎣ ( ) (
⎡ H ′′ − 2ξ H ′ + ξ 2−1 H + λ − ξ 2 H ⎤ ,
⎦ ) (5.11)

()
por lo que H ξ satisface la ecuación diferencial

H ′′ − 2ξ H ′ + λ −1 H = 0 . ( ) (5.12)

Esta ecuación diferencial produce una relación de recurrencia de dos términos, la cual si
permite obtener soluciones cerradas.

Aplicando el método de Frobenius, en donde C0 ≠ 0 , se obtiene

∞ ∞
H = ∑ C jξ j+s
, ξ H ′ = ∑C j s + j ξ ( ) j +s
,
j=0
j=0

∞ ∞

(
H ′′ = ∑ C j s + j s + j −1 ξ )( ) j +s −2
( )(
= ∑ C k+2 s + k + 2 s + k +1 ξ s +k , ) (5.13)
j=0 k=−2

5-4
en donde k = j − 2 . Por lo que

( )
⎡C s −1 sξ −2 + C s +1 sξ −1 +
⎢ 0 1 ( ) ⎤
⎥ s
⎥ξ = 0 .
{ )}
⎢∞ (5.14)
( )( )
⎢ ∑ C k+2 s + k + 2 s + k +1 + C k ⎡⎣ λ −1 − 2 s + k ⎤⎦ ξ ⎥
k
( ) (
⎣ k=0 ⎦

Al igualar a cero cada coeficiente, se tienen las ecuaciones siguientes:

(
s s −1 = 0 , ) s = 0,1 ,

1
(
C s s +1 = 0 , ) C = 0,
1
C
2k+1
=0 , (5.15)

2s + 2k − λ +1
C k+2 =
(
C ,
s + k + 2 s + k +1 k )( )
(k = 0,2,4,…) .


Dado que la solución debe tener paridad, sólo hay términos pares. Así, la solución puede
escribirse en la forma

∞
4m+ 2s +1− λ
H =ξ s
∑A ξ m
2m
, Am+1 =
(
A ,
2m+ 2+ s 2m+1+ s m )( )
(m = 0,1,2,…) ,

(5.16)
m=0

en donde Am = C2m y s = 0,1 .

5.1.3. Los valores propios

Dado que la solución debe ser finita, se compara la serie obtenida con una función
exponencial,

∞
αm Dm+1 α α
e αξ 2
= ∑ Dmξ 2m
, Dm = , = ~ . (5.17)
m! Dm m+1 m
m=0

En nuestro caso,

Am+1 4m 1
~ = , (5.18)
Am 4m2 m

ξ2
y comparando ambos casos se tiene que α = 1 . Así, para valores grandes de ξ , H ~ e y
2
ξ
u ~ e
2
→ ∞ . Entonces, si H es una serie, la función de onda diverge. Por lo tanto, H

5-5
debe ser un polinomio; es decir, existe un entero N , tal que AN+1 = AN+2 =!= 0 , con

A ≠ 0 . Al sustituir esta condición en la relación de recurrencia, ecuación (5.16), se tiene
N
que

4N + 2s +1− λ
0=
( )(
A ,
2N + 2+ s 2N +1+ s N )
( N = 0,1,…) ,

(5.19)

por tanto,

λ = 4N +1+ 2s , ( N = 0,1,…) .

(5.20)

O bien,

⎧4N +1, estado par s = 0

⎪ ( )
λ=⎨ ⇒ λn = 2n +1, (n = 0,1,2,…) . (5.21)

⎪⎩4N + 3, estado impar s = 1 ( )
Esto implica que el espectro es discreto e igualmente espaciado,

1
(
En = !ω 2n +1 = !ω n + 12 ,
2
) ( ) (5.22)

además,

(α x )
2

⎛ξ⎞ 1

() () ( ) ( )
− ξ2 −
un x = un ⎜ ⎟ = u ξ = Nne 2 Hn ξ = Nne 2 Hn α x , (5.23)
⎝α⎠

en donde las funciones Hn son polinomios de grado n , denominados polinomios de

Hermite (Anexo 5.2), y los coeficientes N n garantizan la normalización.

5.1.4. La normalización
Dado que las funciones propias deben estar normalizadas, el coeficiente de
normalización toma la forma

∞ ∞

() ( ) dαξ .
2 2
Nn = ∫ e −ξ2
∫ e Hn ξ
2
−2 −ξ
Hn ξ dx = (5.24)
−∞ −∞

5-6
Este tipo de integrales se puede obtener mediante el uso de la función generadora de los
polinomios de Hermite (Anexo 5.2).1 Considere la integral

∞ ∞
s nt m
( )( )
I = ∫ G s,ξ G t ,ξ e − ξ d ξ = ∑
2

n,m n!m! −∞
∫
2

() ()
e − ξ Hn ξ Hm ξ d ξ
−∞
. (5.25)
∞

= ∫ exp ⎡⎣ −ξ
2
( ⎦) (
+ 2ξ s +t − s +t ⎤ d ξ 2 2
)
−∞

( ) ( )⎦
2
Completando el cuadrado en el exponente, ξ − 2ξ s +t + s +t = ⎡ξ − s +t ⎤ − 2st , se
2 2 2

⎣
tiene que

∞ ∞
( )
2
− ⎡⎣ξ s+t ⎤⎦
∫e dξ = e ∫e
−w
I =e 2st 2st
dw = π e2st , (5.26)
−∞ −∞

( )
en donde w = ξ − s +t . Desarrollando este resultado en una serie de potencias, se tiene

que

∞ ∞
2n s nt n ∞ ∞ s nt m
I= π∑
n!
= ∑ ∑ ∫ e − ξ Hn ξ d ξ
2

n!m!
, () (5.27)
n=0 n=0 m=0 −∞

e igualando coeficientes de ambas series,

⎧⎪ 0 ,n ≠ m α
() ()
∞
∫−∞ n m
−ξ2
e H ξ H ξ dξ = ⎨ = π 2n n!δ mn = δ mn , (5.28)
π 2n n! ,n = m Nn2
⎩⎪
por lo tanto,

1 α α
Nk2 = α = , Nk = . (5.29)
()
∞
∫ e −ξ2
H ξ dξ
2 π 2k k! π 2k k!
−∞ k

Observe que estos coeficientes cumplen las relaciones siguientes

Ni 2 j−i j! Ni Ni
Nj
=
i!
;
Ni+1
= 2 i +1 ;( ) Ni−1
=
1
2i
. (5.30)

1
La función generadora de los polinomios de Hermite (Anexo 5.2) tiene la forma

( )
∞ H ξ
G ξ ,s ≡ e ( ) = e − s +2sξ ≡ ∑ n s n .
ξ 2 − s−ξ
2
2 ()
n=0 n!

5-7
Finalmente, las funciones propias toman su forma final,

α ( )
1

( )
2
− αx
k = e 2
Hk α x , n m = δ nm . (5.31)
π 2k k!

5.1.5. Los elementos de matriz

Para calcular los brakets se pueden utilizar las relaciones de recurrencia de los
polinomios de Hermite (Anexo 5.2).2 Por ejemplo, para los elementos de matriz del
operador diferencial,

1
d ⎡ −12ξ 2 ⎤
d
() ()
− ξ2
n m = ∫ NnNme Hn ξ
2
⎢e Hm ξ ⎥dx
dx dx ⎣ ⎦
∞ − ξ
1 2
dξ d ⎡ − 2ξ
1 2
⎤ dξ
= NnNm ∫ e 2 Hn (ξ ) ⎢ e Hm (ξ )⎥ . (5.32)
−∞ dx dξ ⎣ ⎦α
∞
= NnNm ∫ e − ξ Hn ⎡⎣ −ξ Hm + Hm′ ⎤⎦ dξ
2

−∞

1
Dado que, −ξ Hm + Hm′ = − Hm+1 + mHm−1 , entonces
2

d ∞ ⎡ 1 ⎤
m = NnNm ∫ e − ξ Hn ⎢ − Hm+1 + mHm−1 ⎥ dξ
2
n
dx −∞
⎣ 2 ⎦
⎧ 1 ⎫
= NnNm ⎨− π 2n n!δ n,m+1 + m π 2n n!δ n,m−1 ⎬
⎩ 2 ⎭
. (5.33)
1 α α
= − NnNn−1 2 δ n,m+1 + mNnNn+1 2 δ n,m−1
2 Nn Nn
α Nm N m+1 m
=− δ n,m+1 + mα m δ n,m−1 = −α δ n,m+1 + α δ
2 Nm+1 Nm−1 2 2 n,m−1

En forma similar se pueden calcular el braket

1 m+1 1 m
nxm = δ n,m+1 + δ , (5.34)
α 2 α 2 n,m−1

2
Los polinomios de Hermite tienen relaciones de recurrencia entre polinomios de diferente
grado y con sus derivadas (Anexo 5.2):
()
H ξ = 2ξ Hk ξ − 2kHk−1 ξ ,
k+1
() () () ()
H ′ ξ = 2kHk−1 ξ .
k

5-8

(
y se puede verificar que se cumple la igualdad m E n − E k n x k = −! n ) 2 d
dx
k . Los

promedios de la energía cinética y potencial también pueden calcularse por este método,

!2 2 !2 2 2n +1 !2 mω

T = n−
2m
D n =
2m
α
2
=
2m !
( )
n + 12 = 12 !ω n + 12 = 12 En , ( ) (5.35)

1 ⎛ 2n +1 ⎞ 1
V = n 12 mω 2 x 2 n = 12 mω 2 2⎜
α ⎝ 2 ⎠ ⎟ = 2 mω 2
!
mω
(
n + 12 = 12 En = T .) (5.36)

Esta expresión corresponde a la expresión del teorema virial para el oscilador armónico.

5.2. Los operadores de creación y de aniquilación

La ecuación diferencial del oscilador armónico se puede resolver usando un método
algebraico, aprovechando las propiedades de conmutación de los operadores.

( )
En la ecuación diferencial reescalada (5.4), u′′ − ξ 2u = D2 − ξ 2 u = 0 , el operador de

segundo orden tiene la forma algebraica de una diferencia de cuadrados. Esta observación
sugiere usar los operadores que provendrían de la factorización hipotética,

1 ⎛ ∂ ⎞ 1 ⎛ ∂ ⎞
â ≡ ⎜⎝ ξ + , â+ ≡ ξ − . (5.37)
2 ∂ξ ⎟⎠ 2
⎜⎝ ∂ξ ⎟⎠

En términos de la variables físicas, estos operadores tiene la forma

1 ⎛ 1 d⎞ mω 1
â = ⎜⎝ α x + α dx ⎟⎠ = 2! x̂ + i 2m!ω p̂
2

1 ⎛ 1 d⎞ mω 1
â+ = ⎜⎝ α x − = x̂ − i p̂ , (5.38)
2 α dx ⎟⎠ 2! 2m!ω

que pueden considerarse como parte de la factorización del operador hamiltoniano,

p̂2 1
Ĥ = + mω 2 x̂ 2 . (5.39)
2m 2

5-9
 +
A partir de su definición, es posible demostrar que los operadores â y â no son
hermitianos, pero uno es el adjunto del otro. Dado que estos operadores son
combinaciones de los operadores físicos, están relacionados con ellos de forma simple ,

x̂ =
!
2mω
(
â + â+ , )
p̂ = i
mω ! +
2
(
â − â .) (5.40)

Adicionalmente,

mω 2 i i1 1 1 ⎡ p̂2 mω 2 x̂ 2 ⎤ i
â+ â = x̂ + x̂p̂ − p̂x̂ + p̂2 = ⎢ + ⎥ + ⎡ x̂ , p̂ ⎤⎦
2! 2! 2! 2m!ω !ω ⎣ 2m 2 ⎦ 2! ⎣
, (5.41)
Ĥ 1
= −
!ω 2

por tanto,

(
Ĥ = !ω â+ â + 12 . ) (5.42)

Los nuevos operadores no conmutan,

⎡ mω p̂ mω p̂ ⎤
⎡â, â+ ⎤ = ⎢ x̂ + i , x̂ − i ⎥
⎣ ⎦
⎢⎣ 2! 2mω ! 2! 2mω ! ⎥⎦ , (5.43)
1 1 i
= i ⎡⎣ p̂, x̂ ⎤⎦ − i ⎡⎣ x̂ , p̂ ⎤⎦ = − ⎡⎣ x̂ , p̂ ⎤⎦ = 1
2! 2! !

+ +
así, a partir del conmutador, ââ − â â = 1 , se tiene que

ââ+ = â+ â +1 =
Ĥ 1
+ ,
!ω 2
(
Ĥ = !ω ââ+ − 12 . ) (5.44)

+
Definiendo al operador de número, N̂ ≡ â â , se puede observar que éste está relacionado
con el hamiltoniano,

(
Ĥ = !ω N̂ + 12 , ) ⎡ Ĥ, N̂ ⎤ = 0 .
⎣ ⎦ (5.45)

Analizando las propiedades de conmutación de estos objetos,

⎡ N̂ , â ⎤ = ⎡â+ â, â ⎤ = â+ ⎡â, â ⎤ + ⎡â+ , â ⎤ â = −â , ⎡ Ĥ, â ⎤ = −!ω â ,

⎣ ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦

5-10
 ⎡ N̂ , â+ ⎤ = ⎡â+ â, â+ ⎤ = â+ ⎡â, â+ ⎤ = â+ , ⎡ Ĥ, â+ ⎤ = !ω â+ , (5.46)
⎣ ⎦ ⎣ ⎦ ⎣ ⎦ ⎣ ⎦

se puede observar que los conmutadores involucran a los mismos operadores.

5.2.1. La ecuación de valores propios

Dado que Ĥ y N̂ conmutan, tienen los mismos kets propios,

N̂ n = Mn n , Ĥ n = En n , (5.47)

entonces,

Ĥ 1 E
n N̂ n = n − n = n − 12 = n â+ â n = âun âun ≥ 0 , (5.48)
!ω 2 !ω

por lo tanto, el espectro está acotado, E n ≥ 12 !ω . Aplicando â a la ecuación valores

propios de la energía,

âĤ n = Enâ n , (Ĥâ − ⎡⎣Ĥ,â⎤⎦) n = Ĥâ n + !ω â n . (5.49)

Así,

( ) (
Ĥ â n = En − !ω â n , )( ) (5.50)

es decir, â n es la función propia con valor propio E n − !ω . Por esta razón, al operador â

se le denomina operador de aniquilación. En forma similar,

â+ Ĥ n = Enâ+ n , (Ĥâ − ⎡⎣Ĥâ ⎤⎦) n = Ĥâ
+ + +
n − !ω â+ n , (5.51)

por lo tanto,

( ) (
Ĥ â+ n = En + !ω â+ n , )( ) (5.52)

+
y â n es la función propia con valor propio E n + !ω , por lo que el operador â se le
+

llama operador de creación.

5-11
 Como el espectro está acotado por abajo, denotaremos con 0 a la función propia con

menor valor propio, E0 . Por otro lado, el ket â 0 no puede ser un ket propio, ya que no

existe ninguno con energía menor que E0 , por lo tanto este ket debe ser igual al ket nulo,

â 0 = 0 . Además,

Ĥ 1 E
0 N̂ 0 = 0 â+ a 0 = au0 au0 = 0 = 0 − 2 0 = 0 − 12 , (5.53)
!ω !ω

por lo tanto E0 = 12 !ω . Los valores propios siguientes se obtienen aplicando el operador

de creación,

E = E0 + n!ω = !ω n + 12 .
n
( ) (5.54)

Observe que, N̂ n = n n .

Las funciones propias también se obtiene aplicando los operadores de creación y
aniquilación

â n = cn n −1 , â+ n = dn n +1 , (5.55)

en donde las constantes están fijas por la normalización,

2 2
âun âun = cn n −1 n −1 = cn = n â+ â n = n N̂ n = n ,

2 2
â+un â+un = dn n +1 n +1 = dn = n ââ+ n = n N̂ +1 n = n +1 , (5.56)

entonces, cn = n y dn = n +1 , por lo que

â+ n = n +1 n +1 , â n = n n −1 , (5.57)

y

1
n = â+n 0 . (5.58)
n!

5-12
5.2.2. La función propia del estado basal

El ket propio del estado basal satisface la ecuación â 0 = 0 , por lo tanto,

1⎛ d ⎞
⎜ ξ + ⎟ 0 = 0. (5.59)
2⎝ dξ ⎠

Esta es una ecuación diferencial separable de primer orden,

du0 du0
= −ξu0 , = −ξ dξ . (5.60)
dξ
u0

Integrando,

⎛ ξ2 ⎞
()
u0 ξ = Aexp ⎜ − ⎟ ,
⎝ 2⎠
(5.61)

en donde A se determina por la normalización de la función de onda en el espacio de
coordenadas,

A2 ∞ ξ 2 A2
()
∞ 2
1 = 0 0 = ∫ u0 ξ
α ∫−∞
dx = −e d ξ = π. (5.62)
−∞ α

Así, A2 = α π

α 1 2 − 12(α x )2
()
0 = u0 x =
π1 4
e , (5.63)

y

n n
1 ⎛ d ⎞ α ⎛ d ⎞ − 1ξ 2
n = ⎜⎝ ξ − dξ ⎟⎠ 0 = ⎜ ξ− ⎟ e 2 .
dξ ⎠
(5.64)
2n n! π n!2n ⎝

5.3 Problemas
1. A partir de la fórmula de Rodrigues, obtenga los seis primeros polinomios de Hermite

(n = 0,1,…,5) y calcule las raíces de cada uno de ellos. Haga una gráfica cualitativa de las

( ) ()
2
funciones exp −ξ 2 Hn ξ .

5-13
2. Para el oscilador armónico calcule

d
kxl y k l .
dx

Demuestre que

n T̂ n = 12 En y n V̂ n = 12 En .

3. Una partícula de masa de 100 g está unida a un resorte con constante de fuerza igual a
1000 Jm-2. Calcule: (a) la energía de estado basal, (b) p̂2 , (c) x̂ 2 , (d) la separación

entre los niveles de energía.

5-14
Anexo 5.1. El comportamiento asintótico de u′′ + ( λ − ξ 2 )u = 0

Para este análisis es conveniente transformar la ecuación diferencial en un sistema de
ecuaciones de primer orden. Por tanto se usarán las variables u y y ,

y = u′ , (
y ′ = u′′ = ξ 2 − λ u ,)
que se agrupan en un vector de funciones,

⎛u ⎞
X=⎜ ⎟ .
⎝ y⎠
Así,

⎛ u′ ⎞ ⎛ y ⎞ ⎛ 0 1 ⎞
X′ = ⎜ ⎟ = ⎜ 2 ⎟ =⎜ 2 ⎟ X = BX ,

(
⎝ u′′ ⎠ ⎝ ξ − λ ) u⎠ ⎝ ξ − λ 0 ⎠

en donde

⎛ 0 1 ⎞ ⎛ 0 0 ⎞ ⎛ 0 1 ⎞ ⎡ 1 ⎤
B=⎜ 2 ⎟ =ξ ⎜ +⎜ = ξ 2 ⎢B0 + 2 B2 ⎥ ,
2

ξ −λ 0 ⎠ ⎟ ⎟
⎝ 1 0 ⎠ ⎝ −λ 0 ⎠ ξ
⎝ ⎣ ⎦

B , B son matrices constantes y se ha factorizado la potencia mayor de ξ .

0 2

Para el análisis asintótico se requiere que los valores propios de B0 sean diferentes, por

lo que se propone la transformación

⎛ ξa 0 ⎞ ⎛ ξ au ⎞
Z=⎜ ⎟X=⎜ b ⎟ ,
⎜⎝ 0 ξ b ⎟⎠ ⎜⎝ ξ y ⎟⎠

entonces,

⎛ ξ a−b ⎞ ⎛
aξ aξ ξ a−b ⎞
Z′ = ⎜ b−a 2 ⎟Z=⎜ ⎟Z.

⎜ ξ
⎝
ξ −λ( ) b ξ ⎟
⎠
⎜⎝ ξ b−a+2 − λξ b−a b ξ ⎟⎠

5-15
 Dado que los valores propios de la matriz de mayor potencia tienen que ser diferentes,
entonces, por simplicidad, se puede elegir que la potencia mayor de los elementos no

diagonales sea igual. Entonces, b − a + 2 = a − b , y por tanto 2 b − a +1 = 0 o a = b +1 . Así,

( )
⎛ aξ ξ ⎞ ⎛ 0 ξ ⎞ ⎛ aξ 0⎞
Z′ = ⎜ ⎟ Z = ⎜ ⎟ Z + ⎜ ⎟Z
⎜⎝ ξ − λ ξ b ξ ⎟⎠ ⎜⎝ ξ 0 ⎟⎠ ⎜⎝ − λ ξ b ξ ⎟⎠
.
⎡ ⎛ ⎞ 1⎛ ⎞⎤ ⎡⎛ ⎞ 1⎛ ⎞⎤
= ⎢ξ ⎜ 0 1 ⎟ + ⎜ a 0 ⎟ ⎥ Z = ξ ⎢ ⎜ 0 1 ⎟ + 2 ⎜ a 0 ⎟ ⎥ Z
⎢⎣ ⎝ 1 0 ⎠ ξ ⎝ − λ b ⎠ ⎥⎦ ⎢⎣⎝ 1 0 ⎠ ξ ⎝ − λ b ⎠ ⎥⎦

Para la primera matriz, los valores propios son ±1 . Si se toma la combinación a = 1 y

b = 0 , finalmente se tiene que

⎛ 1 ⎞ ⎛ ⎞ ⎛ ⎞
Z′ = ξ ⎜ A 0 + 2 A 2 ⎟ Z, A0 = ⎜ 0 1 ⎟ , A2 = ⎜ 1 0 ⎟ .
⎝ ξ ⎠ ⎝ 1 0 ⎠ ⎝ −λ 0 ⎠

Para esta ecuación, se propone una solución asintótica de la forma

∞ ∞
ξ µ ∑ Ckξ −k = e 0 ∑C ξ
a0ξ 2 +a1ξ a ξ 2 +a1ξ µ −k
Z=e k
.
k=0 k=0

Entonces

⎡∞ ∞
⎤
Z′ = e ⎢∑
a0ξ 2 +a1ξ
µ − (
k C k
ξ )
µ −k−1
+ 2a0
ξ + a1 ∑ k(C ξ µ−k ⎥ )
⎣ k=0 k=0 ⎦
a ξ 2 +a ξ ⎡ ⎤
∞ ∞ ∞

( )
= e 0 1 ⎢ ∑ µ − l +1 Cl−1ξ −l + 2a0 ∑ Cm+1ξ −m + a1 ∑ Ckξ −k ⎥ ξ µ ,
⎣ l=1 m=−1 k=0 ⎦
⎡ ⎤
{( }
∞
= e 0 1 ξ µ ⎢2a0C0ξ + 2a0C1 + a1C0 + ∑ µ − k +1 Ck−1 + a1Ck + 2a0Ck+1 ξ −k ⎥
a ξ 2 +a ξ
)
⎣ k=1 ⎦

y sustituyendo en la ecuación diferencial se tiene que

⎛ 1 ⎞ a ξ 2+a ξ ∞
( )
∞
Z′ = ξ ⎜ A 0 + 2 A 2 ⎟ e 0 1 ξ µ ∑ Ckξ −k = e 0 0 ξ µ ∑ A 0C kξ 1−k + A 2Ckξ −k−1
a ξ 2 +a ξ

⎝ ξ ⎠ k=0 k=0

⎡ ∞ ∞
⎤
ξ µ ⎢ ∑ A 0Cm+1ξ −m + ∑ A 2Cl−1ξ −l ⎥
a0ξ 2 +a1ξ
=e .
⎣ m=−1 l=1 ⎦
⎡ ∞
⎤
= e 0 1 ξ µ ⎢ A 0C0ξ + A 0C1 + ∑ A 0Ck+1 + A 2Ck−1 ξ −k ⎥
a ξ 2 +a ξ
( )
⎣ k=1 ⎦

5-16
Al igualar los coeficientes de ambos polinomios:

2a C = A 0C0 , 2a C + a C = A 0C1 ,
0 0 0 1 1 0

( )
2a C + a C + µ − k +1 Ck−1 = A 0Ck+1 + A 2Ck−1 ,
0 k+1 1 k
(k = 1,2,…) .

De la primera ecuación, 2a0C0 = A 0C0 , se puede observar que 2a0 corresponde al valor

propio de A 0 , por lo tanto 2a0 = ±1 o a0 = ± 12 . Las otras ecuaciones quedan en la forma

(A
0 )
− 2a0I C1 = a1C0 , (A
0 ) ( )
− 2a0I Ck+1 = a1Ck + ⎡⎣ µ − k +1 I − A 2 ⎤⎦ Ck−1 .

Los vectores propios de A 0 se encuentran resolviendo la ecuación de valores propios:

⎛ −2a 1⎞⎛ C ⎞ ⎛ −2a C + C ⎞

(
0 = A 0 − 2a0I C0 = ⎜)
⎜⎝ 1
0
⎟⎜
0,1

−2a0 ⎟⎠ ⎜⎝ C0,2
⎟ =⎜ 0 0,1 0,2

⎟⎠ ⎜⎝ C0,1 − 2a0C0,2
⎟,
⎟⎠

⎛ 1 ⎞
c ≡C , C = 2a0c0 , C 0 = c0 ⎜ ⎟.
0 0,1 0,2 ⎝ 2a0 ⎠

Entonces,

⎛ −2a0C1,1 + C1,2 ⎞ ⎛ a1c0 ⎞

(A 0
− 2a0
I C1 )
= ⎜ ⎟ = a C
1 0
= ⎜ ⎟,
⎝ C1,1 − 2a0C1,2 ⎠ ⎝ 2a0a1c0 ⎠

( ) = 1 , entonces la ecuación del segundo renglón puede transformarse en

2
Como 2a0

2a C − C = a c y al sumarla a la del primero, obteniendo 2a1c0 = 0 . Por lo tanto a1 = 0 .
0 1,1 1,2 1 0

Para k = 1 ,

⎛ −2a0C2,1 + C2,2 ⎞ ⎛ µ −1 0 ⎞ ⎛ 1 ⎞ ⎛ c0 µ −1 ( ) ⎞
( )
A 0 − 2a0I C2 = ⎜ ⎟ = c0 ⎜ ⎟⎜ ⎟ =⎜ ⎟.

C
⎝ 2,1 − 2a C
0 2,2 ⎠
⎜
⎝ λ µ ⎟⎠ ⎝ 2a0 ⎠ ⎜ c λ + 2a µ
⎝ 0 0 ( ) ⎟
⎠

Multiplicando la ecuación del segundo renglón por 2a0 y sumándola a la del primero, se

obtiene

5-17
 (
0 = λ 2a0 + 2µ −1 c0 ,

)
( )
µ = 1− 2a0λ 2 .

⎛ ξ u⎞ ∞
Como Z = ⎜
⎝y ⎠ ⎟ , de la solución asintótica, u ~ e
a0ξ 2 µ −1
ξ ∑ C k ,1ξ −k , pero u sólo es finita
k=0

para a0 < 0 . Por lo tanto, tomando el valor propio negativo, 2a0 = −1 , la solución

asintótica queda en la forma

1 λ −1

() ( )
− ξ2
u ξ ~e 2
ξ 2
c0 1+! .

5-18
Anexo 5.2. Los polinomios de Hermite, Hn (ξ )

Los polinomios de Hermite son funciones que satisfacen la ecuación diferencial

H ′′ − 2ξ Hn′ + 2nHn = 0 . Ahora se mostrará que estos polinomios también se pueden obtener
n
a partir de una función generadora.

La función generadora de los polinomios de Hermite

( )
Sea G ξ ,s una función de dos variables definida como

( )
H ξ
G ξ ,s ≡ e ( ) = e − s +2sξ ≡ ∑ n s n .
ξ 2 − s−ξ
2
2 ∞
()
n=0 n!

A esta función se le llama función generadora debido a que los polinomios Hn se obtienen

de los coeficientes del desarrollo en series de Taylor en la variable s . Las propiedades de

los polinomios de Hermite Hn también provienen de la función generadora. Por ejemplo,

∂ G ∂ξ = 2sG , así

∞ Hn′ ξ ( )s ∞ Hn ξ ( )s ∞
( )s
Hk−1 ξ
∑ n
= 2∑ n+1
= 2∑ k

n=0 n! n=0 n! k=1 (k −1)!
en donde se usó el cambio de índice k = n +1 . Igualando los coeficientes de ambas series,
se obtiene

H ′ (ξ ) = 0 ,
0 k
()
H ′ ξ = 2kHk−1 ξ , ( ) (k = 1,2,…) .
En forma similar, ∂ G ∂ s = 2 ξ − s G

( )
∞
nHn s n−1 ∞
ξ Hn s n ∞
H s n+1
∑ n!
= 2∑
n!
− 2∑ n
n!
n=1 n=0 n=0
,
∞
Hl+1 s l ∞
ξ Hn s n ∞
H sk
∑ = 2∑ − 2∑ k−1

l=0 l! n=0 n! k=1 k −1 ! ( )
en donde se hicieron las transformaciones de índices l = n −1 , k = n +1 . Comparando los
coeficientes de las series:

5-19
 ()
H ξ = 2ξ H0(ξ ) ,
1
() ()
H ξ = 2ξ Hk ξ − 2kHk−1 ξ ,
k+1
() (k = 1,2,…) .

Las propiedades anteriores permiten obtener la ecuación diferencial de los polinomios Hn ,

′ ′
′ = ⎡⎣2nHn−1 ⎤⎦ = ⎡⎣2ξ Hn − Hn+1 ⎤⎦ = 2Hn + 2ξ Hn′ − Hn+1
Hn′′ = 2nHn−1 ′
,

( )
= 2Hn + 2ξ Hn′ − 2 n +1 Hn = 2ξ Hn′ − 2nHn

así,

H ′′ − 2ξ Hn′ + 2nHn = 0 .
n

( ) ( )
Por lo tanto, la función G ξ ,s ≡ exp −s 2 + 2sξ es la función generadora de los

polinomios de Hermite.

La fórmula de Rodrigues
Dado que los polinomios son los coeficientes del desarrollo en series de Taylor de la
función generadora, G , es posible obtener una expresión explícita para ellos. Al

( )
desarrollar a G ξ ,s con respecto a s se tiene que

∞

( )
G ξ ,s = ∑ an s n ,
n=0

en donde

1 ∂n
an =
n! ∂ s n
G ξ ,s ( ) ,
s=0

por lo tanto

⎛ ⎞
∂n ∂ n ξ 2−( s−ξ )2 ∂ n −t 2 ∂n
() ( ) −( − ξ )
2
ξ2 ξ2 ⎜ ⎟
Hn ξ = n G ξ ,s = ne =e e =e e
∂s ∂s ∂ ⎜
( ) ⎟
n n
t
s=0 s=0 t=− ξ ⎝ ∂ −ξ ⎠,
2
d ne − ξ
( )
n
ξ2
=e −1
dξ n

5-20
 2
d ne − ξ
() ( )
n
ξ2
en donde se usa el cambio de variable t ≡ s − ξ . A la ecuación Hn ξ = e −1 se le
dξ n
denomina la fórmula de Rodrigues de los polinomios de Hermite. Esta fórmula permite
evaluar sucesivamente a todos los polinomios entre 0 y n , por ejemplo, aplicando la
ecuación para n = 2 , se tiene que

() ()
H0 ξ = 1 H1 ξ = 2ξ () ( )
H2 ξ = 2 2ξ 2 −1 .

De igual forma, se pueden usar la relación de recurrencia para evaluar otros polinomios,
por ejemplo,

( )
H = 2ξ H1 − 2H0 = 2ξ 2ξ − 2⋅1 .
2
Es importante mencionar que las raíces de los polinomios de Hermite siempre son reales.

5-21