UNIVERSIDAD NACIONAL DE COLOMBIA

Facultad de Ciencias –Departamento de Física

Oscilaciones y Ondas –2016682
TAREA No. 7
GRUPO: No. 4, SECCIÓN: 02, Profesor: Héctor Castro
Integrantes:
Juan Sebastián Castro Millán
Andrés Felipe Escobar de la Hoz
Jaime Felipe Murillo Espitia
Daniel Felipe Peralta Macca
Fecha: 21 de Abril de 2021

1. En la solución hallada para el estado base del oscilador armónico (diapositiva 12)
aparece el constante 𝐶0 , la cual se determina mediante la normalización de la función de
𝑚 𝜔0 1/4
onda, obteniendo el valor: 𝐶0 = ( ) . Realice la normalización y compruebe este
ħ𝜋
resultado.

Como solución de la ecuación se propuso:

2 𝑚𝜔
Ψ(𝑥) = 𝐶0 𝑒 −𝛼𝑥 , 𝑐𝑜𝑛 𝛼 =
2ℏ
Para hallar la constante se normaliza la función de onda, con la siguiente fórmula:
+∞
∫ Ψ2 (𝑥)dx = 1
−∞

Acomodando la expresión para 𝐶0 :

1
𝐶𝑜 =
+∞
√∫−∞ 𝑒 −2𝛼𝑥 2 𝑑𝑥

La integral es un caso especial de la función de error de Gauss, por lo que es igual a:

+∞
2 𝜋
∫ 𝑒 −2𝛼𝑥 𝑑𝑥 = √
−∞ 2𝛼
Entonces la constante es:
1
2𝛼4 𝑚𝜔 4
𝐶𝑜 = √ = ( )
𝜋 𝜋ℏ
 Que coincide con el obtenido en clase.

2. Resuma en pocos pasos la deducción de las soluciones generales del oscilador armónico
mediante series de potencia. No debe hacer todas las deducciones, solo mostrar resultados
intermedios hasta llegar a los polinomios de Hermite.

Partiendo de la ecuación de Schrödinger del oscilador armónico:

ℏ2 𝑑 2 𝜓 1
2
= [ 𝑚𝜔2 𝑥 2 − 𝐸] 𝜓 (1)
2𝑚 𝑑𝑥 2
Para facilitar un poco el cálculo, se introducirá un pequeño cambio de variable
𝑚𝜔
𝜉≡√ 𝑥
ℏ

ℏ
𝑥 ≡ 𝜉√
𝑚𝜔

Donde las derivadas de x con respecto a “𝜉” son:

𝑑𝑥 ℏ 𝑑2 𝑥
=√ =0
𝑑𝜉 𝑚𝜔 𝑑𝜉 2

Aplicando regla de la cadena, obtenemos 𝜓 dependiente de 𝜉 y sus respectivas derivadas

𝑑𝜓 𝑚𝜔 𝑑𝜓
=√
𝑑𝑥 ℏ 𝑑𝜉
𝑑 2 𝜓 𝑚𝜔 𝑑 2 𝜓
=
𝑑𝑥 2 ℏ 𝑑𝜉 2
Y reemplazando en (1):

ℏ2 𝑚𝜔 𝑑 2 𝜓 1 2 √
ℏ 2
= [ 𝑚𝜔 (𝜉 ) − 𝐸] 𝜓
2𝑚 ℏ 𝑑𝜉 2 2 𝑚𝜔

𝑑2𝜓 2
2𝐸
= [𝜉 − ] 𝜓 (2)
𝑑𝜉 2 𝜔ℏ

2𝐸
Para valores de 𝜉 muy grandes podemos despreciar el aporte de la constante 𝜔ℏ, por lo
que una primera aproximación seria la siguiente:
 𝑑2 𝜓
2
≈ 𝜉2𝜓
𝑑𝜉
𝜉2 𝜉2
−( ) ( )
𝜓(𝜉) = 𝐴𝑒 2 + 𝐵𝑒 2
La condición que la función de onda debe cumplir es que sea normalizable, de alguna u
otra forma se hace necesaria su convergencia. El término que acompaña a “B” no
converge, así que la función de onda debe ser de la forma
𝜉2
−( )
𝜓(𝜉) = ℎ(𝜉)𝑒 2 (3)
Donde ℎ(𝜉), se puede hallar usando el método de serie de potencias:
∞

ℎ(𝜉) = ∑ 𝑎𝑛 𝜉 𝑛
𝑛=0
∞

ℎ′(𝜉) = ∑ 𝑛𝑎𝑛 𝜉 𝑛−1

𝑛=1
∞

ℎ′′(𝜉) = ∑ 𝑛(𝑛 − 1)𝑎𝑛 𝜉 𝑛−2

𝑛=2

Hallando las derivadas de (3), obtenemos:

𝜉 2
𝑑𝜓(𝜉) −( )
= [ℎ′ (𝜉) − 𝜉ℎ(𝜉)]𝑒 2
𝑑𝜉
2
𝑑 2 𝜓(𝜉) ′′ (𝜉) ′(𝜉)
𝜉
−( )
= [ℎ − ℎ(𝜉) − 𝜉ℎ′(𝜉) − 𝜉(ℎ − 𝜉ℎ(𝜉))]𝑒 2
𝑑𝜉 2
Reemplazando en dos y después añadiendo las series de potencias se obtiene lo siguiente:
𝜉2 2𝐸 𝜉 2
′′ (𝜉)−2𝜉ℎ′(𝜉) 2 −( ) 2 −( )
[ℎ + (𝜉 − 1)ℎ(𝜉))]𝑒 2 = [𝜉 − ] ℎ(𝜉)𝑒 2
𝜔ℏ
∞ ∞ ∞
𝑛−2 𝑛−1
2𝐸
∑ 𝑛(𝑛 − 1)𝑎𝑛 𝜉 − 2𝜉 ∑ 𝑛𝑎𝑛 𝜉 +( − 1) ∑ 𝑎𝑛 𝜉 𝑛 = 0
𝜔ℏ
𝑛=0 𝑛=1 𝑛=2
∞ ∞ ∞
𝑛−2
2𝐸
𝑛
∑ 𝑛(𝑛 − 1)𝑎𝑛 𝜉 − 2 ∑ 𝑛𝑎𝑛 𝜉 + ( − 1) ∑ 𝑎𝑛 𝜉 𝑛 = 0
𝜔ℏ
𝑛=0 𝑛=1 𝑛=2
Hacemos una sustitución en el contador de la serie de la segunda derivada y un cambio en
el contador de la primera derivada para poder finalmente agruparlo todo en un solo
termino
∞
2𝐸
∑[(𝑛 + 2)(𝑛 + 1)𝑎𝑛+2 − 2𝑛𝑎𝑛 + ( − 1)𝑎𝑛 ]𝜉 𝑛 = 0
𝜔ℏ
𝑛=0

2𝐸
(𝑛 + 2)(𝑛 + 1)𝑎𝑛+2 − 2𝑛𝑎𝑛 + ( − 1) 𝑎𝑛 = 0
𝜔ℏ
2𝐸
𝑈=
𝜔ℏ
[2𝑛+(1−𝑈)]
𝑎𝑛+2 = (𝑛+2)(𝑛+1) 𝑎𝑛 (4)

Para “a” pares:

(1 − 𝑈) (1 − 𝑈)(5 − 𝑈)
𝑎2 = 𝑎 𝑎4 = 𝑎0 …
2∗1 0 4∗3∗2∗1
Para “a” impares:
(3 − 𝑈) (3 − 𝑈)(7 − 𝑈)
𝑎3 = 𝑎1 𝑎5 = 𝑎 …
3∗2∗1 5∗4∗3∗2∗1 1

Y en general:
1 − 𝑈 2 (1 − 𝑈)(5 − 𝑈) 4
ℎ(𝜉) = 𝑎0 [1 + 𝜉 + 𝜉 + ⋯ ] (𝑃𝑎𝑟𝑒𝑠)
2! 4!
3 − 𝑈 3 (3 − 𝑈)(7 − 𝑈) 5
ℎ(𝜉) = 𝑎1 [1 + 𝜉 + 𝜉 + ⋯ ] (𝑃𝑎𝑟𝑒𝑠)
3! 5!

Nótese que para valores de “n” muy grandes, la sucesión toma la siguiente forma:
𝑎𝑛+2 2
≈
𝑎𝑛 𝑛
𝜉2
𝐷 ( )
Con 𝑎𝑛 ≈ 𝑗 . Esto conllevaría a que ℎ(𝜉) = 𝐷𝑒 2 . Que es precisamente lo que se
( )!
2
quiere evitar ya que el objetivo es que la función pueda normalizarse y como se vio
𝜉2
( )
anteriormente, 𝐷𝑒 2 no converge.
 Debe garantizarse que la sucesión se trunque en algún momento. Para ello debe cumplirse
que
2𝑛 + (1 − 𝑈) = 0
𝑈 = 2𝑛 + 1
2𝐸
= 2𝑛 + 1
𝜔ℏ
1
(𝑛 + ) 𝜔ℏ = 𝐸 (𝐸𝑛𝑒𝑟𝑔𝑖𝑎 𝑐𝑢𝑎𝑛𝑡𝑖𝑧𝑎𝑑𝑎)
2
Para garantizar dicho truncamiento, en valores de “j” pares, 𝑎1 = 0, y en valores de “j”
impares, 𝑎0 = 0.

Para n=0, se tiene U= 1, por lo que ℏ0 (𝜉) = 𝑎0 𝐻0 (𝜉) = 1

Para n=1, se tiene U= 3, por lo que ℏ1 (𝜉) = 𝑎1 𝜉 𝐻1 (𝜉) = 2𝜉
Para n=2, se tiene U=5, por lo que ℏ2 (𝜉) = 𝑎0 (1 − 2𝜉 2 ) 𝐻2 (𝜉) = (4𝜉 2 − 2)

Comparando estos tres primeros resultados con los polinomios de Hermite, vemos que
son de la misma forma si se le da el valor adecuado a los 𝑎0 𝑦 𝑎1 .
Así que en general, ℏ𝑛 (𝜉) = 𝐻𝑛 (𝜉)
𝜉2 𝑚𝜔
−( )
𝜓(𝜉) = 𝐻𝑛 (𝜉)𝑒 2 𝑐𝑜𝑛 𝜉 ≡ √ 𝑥
ℏ
Tener en cuenta que dicha función no está normalizada, pero con todo el procedimiento
anterior se garantiza que la función sea normalizable.

3. Estudie la solución de las vibraciones en un cristal con dos átomos distintos en la celda
unitaria, lo cual da origen a las dos ramas, una acústica y una óptica. ver libro: Introduction
to Solid State Physics, C. Kittel, 7th Edition, Cap. 4 pgs. 104-107. Resuma en pocos pasos
los resultados importantes sin deducirlos.

Ecuaciones de movimiento bajo la condición de que cada plano interactúa solo con su plano
vecino mas cercano
𝑑 2 𝑢𝑠
𝑀1 = 𝐶(𝑣𝑠 + 𝑣𝑠−1 − 2𝑢𝑠 )
𝑑𝑡 2
(18)
 𝑑2 𝑣𝑠
𝑀2 2 = 𝐶(𝑢𝑠+1 + 𝑢𝑠 − 2𝑣𝑠 )
𝑑𝑡
Se busca una solución de la forma de una onda viajera, ahora con diferentes amplitudes u
y v en planos alternos
𝑢𝑠 = 𝑢 𝑒𝑥𝑝(𝑖𝑠𝐾𝑎) exp(−𝑖𝑤𝑡)
(19)
𝑣𝑠 = 𝑣 𝑒𝑥𝑝(𝑖𝑠𝐾𝑎) exp(−𝑖𝑤𝑡)
Con a siendo la distancia entre planos idénticos más cercanos, sustituyendo 19 en 18
−𝜔2 𝑀1 𝑢 = 𝐶𝑣[1 + exp(−𝑖𝐾𝑎)] − 2𝐶𝑢
(20)
−𝜔2 𝑀2 𝑣 = 𝐶𝑣[exp(𝑖𝐾𝑎) + 1] − 2𝐶𝑣
Ecuaciones lineales homogéneas que tienen una solución solo si el determinante es 0, o
sea:
𝑀1 𝑀2 𝜔4 − 2𝐶(𝑀1 + 𝑀2 )𝜔2 + 2𝐶 2 (1 − cos(𝐾𝑎)) = 0
1
Para un valor pequeño de Ka tenemos cos (Ka) ≈ 1 - 2 𝐾 2 𝑎2 + ⋯, las dos raíces son

1 1
𝜔2 ≈ 2𝐶( + )
𝑀1 𝑀2
(Rama óptica)
1
𝐶
𝜔2 ≈ 2 𝐾 2 𝑎2
𝑀1 + 𝑀2
(Rama acústica)

4. A partir de la ecuación para la energía total de los fonones U, hallada en la diapositiva 31,
realice la derivada indicada: dU/dT y compruebe la formula hallada para la capacidad
térmica Cv, en la misma diapositiva.

La fórmula dada para el calor especifico o la capacidad térmica en la dispositiva es

𝜕𝑈
𝐶𝑣 = ( )
𝜕𝑇 𝑣
 Donde 𝑈 es el potencial de la red cristalina, tal que
ħ𝜔
𝑈 = ∑ ∫ 𝑑𝜔𝐷𝜆 (𝜔)
𝑒 ħ𝜔⁄𝐾𝐵 𝑇 −1
𝜆

Evaluando la derivada parcial requerida

𝜕𝑈 𝜕 ħ𝜔
𝐶𝑣 = ( ) = ∑ ∫ 𝑑𝜔𝐷𝜆 (𝜔) ħ𝜔⁄𝐾 𝑇
𝜕𝑇 𝑣 𝜕𝑇 𝑒 𝐵 −1
𝜆

Ya que el operador derivada es lineal y la integral no se hace con respecto a 𝑇 ni es parte

de los límites de integración, la derivada parcial pasa directamente al integrando
𝜕 ħ𝜔
𝐶𝑣 = ∑ ∫ 𝑑𝜔𝐷𝜆 (𝜔) ( ħ𝜔⁄𝐾 𝑇 )
𝜕𝑇 𝑒 𝐵 −1
𝜆

𝑓 ′ 𝑓´𝑔−𝑔´𝑓
Por la identidad de la derivada una fracción (𝑔) = 𝑔2

−ħ𝜔 ħ𝜔 2
ħ𝜔(−𝑒 ħ𝜔⁄𝐾𝐵 𝑇 )( ) 𝐾 (𝑒 ħ𝜔⁄𝐾𝐵 𝑇
) (
𝐾𝐵 𝑇 2 𝐵 𝐾𝐵 𝑇 )
𝐶𝑣 = ∑ ∫ 𝑑𝜔𝐷𝜆 (𝜔) = ∑ ∫ 𝑑𝜔𝐷𝜆 (𝜔)
(𝑒 ħ𝜔⁄𝐾𝐵 𝑇 − 1)2 (𝑒 ħ𝜔⁄𝐾𝐵 𝑇 − 1)2
𝜆 𝜆
ħ𝜔
Sea 𝑥 = 𝐾
𝐵𝑇

𝐾𝐵 𝑒 𝑥 𝑥 2
𝐶𝑣 = ∑ ∫ 𝑑𝜔𝐷𝜆 (𝜔) 𝑥
(𝑒 − 1)2
𝜆

Aunque este resultado no sea idéntico al obtenido en clase, solo difieren en un factor 𝐾𝐵 .