Materia Condensada Año 2020

Práctica 4: Vibraciones de la red - fonones

1. Problema:
Estudiar la dinámica de la red de la cadena lineal diatómica de la Figura. Examinar en
particular los modos en el centro y borde de la primera zona de Brillouin. Tener en cuenta
los siguientes casos particulares:

(a) M1 = M2 , G = K (d = a/2)
(b) K  G
(c) M1 = M2 , G > K
(d) G = K, M1 > M2

Solución:
Tenemos que pensar a la cadena lineal diatómica como una red unidimensional de con-
stante de red a con un motivo de dos átomos. Llamando u1,n al desplazamiento del átomo
1 de la celda primitiva n, escribimos la energı́a potencial del sistema,
N N
KX 2 GX
U= (u1,n − u2,n ) + (u2,n − u1,n+1 )2
2 n=1 2 n=1

donde tuvimos en cuenta N celdas primitivas y condiciones periódicas de contorno que

indican:
u1,n+N = u1,n u2,n+N = u2,n
Ahora escribimos la ecuación de movimiento de los átomos de una celda arbitraria n
usando
d2 ux,n dU
M 2
=− x = 1, 2
dt dux,n
Entonces obtenemos,

M1 ü1,n = −K (u1,n − u2,n ) + G (u2,n−1 − u1,n )

M2 ü2,n = K (u1,n − u2,n ) − G (u2,n − u1,n+1 )

Tenemos entonces un sistema de dos ecuaciones diferenciales lineales de segundo orden.
Proponemos una solución que represente la oscilación en modos normales

u1,n = A1 ei(kna−ωt) u2,n = A2 ei(kna−ωt)

1
Materia Condensada Año 2020

Reemplazando la solución propuesta en la ecuación diferencial, se transforma en el sigu-

iente sistema de 2 ecuaciones con dos incognitas

K + G e−ika
 2    
[ω M1 − (K + G)] A1
=0
K + G eika [ω 2 M2 − (K + G)] A2

que tiene solución no trivial cuando la matriz tiene determinante nulo. De plantear la
nulidad del determinante se llega a la condición para ω

ω 4 (M1 M2 ) + ω 2 (−(K + G)(M1 + M2 )) + (2KG(1 − cos ka)) = 0

Ahora, usando la resolvente para ω 2 se obtienen los siguiente valores dependiendo de k,

p
2 (K + G)(M1 + M2 ) ± (K + G)2 (M1 + M2 )2 − 8(M1 M2 )(KG)(1 − cos ka)
ω =
2M1 M2

Para poder hacer un análisis funcional de los modos en el centro y borde de la primer
zona de Brillouin, vamos a tomar M1 = M2 = M para simplificar un poco las expresiones.
Ahora las frecuencias angulares resultan
K +G 1√ 2
ω2 = ± K + G2 + 2KG cos ka
M M
En el centro de la zona de Brillouin k = 0, tenemos

2(K + G)
ωa2 (k = 0) = 0 ωo2 (k = 0) =
M
donde los subı́ndices a y o significan acúsitco y óptico, respectivamente (nombre que
tendrá sentido más adelante). Por otro lado, en el borde de la zona de Brillouin ka = π
tenemos
2G 2K
ωa2 (k = π/a) = ωo2 (k = π/a) =
M M
donde suponemos que K > G, sin perdida de generalidad. Ahora que tenemos los ex-
tremos de la relación de dispersión ω(k), estudiemos mejor el comportamiento en torno a
estos puntos.
Empecemos viendo que pasa cuando ka  π, entonces podemos aproximar cos ka '
1 − 21 a2 k 2 . De esta manera, la relación de dispersión queda
s !
K +G KGa2 2
ω2 = 1± 1− k
M (K + G)2
√
Ahora vamos a usar varias veces la expansión de Taylor 1 + x2 ' 1 + 21 x2 (para cambiar
el signo hay que reemplazar x → ix). Primero obtenemos

KGa2 2(K + G) KGa2

ωa2 (k  π/a) = k2 ωo2 (k  π/a) = − k2
2(K + G)M M 2(K + G)M

2
Materia Condensada Año 2020

Y cuando calculamos la raiz obtenemos (utilizando Taylor nuevamente para ωo )

s r
KGa2
 
KG 2(K + G) 2
ωa (k  π/a) = ak ωo (k  π/a) = 1− k
2(K + G)M M 8(K + G)2

Se puede ver una diferencia fundamental, ωa crece linealmente a partir de 0 en el origen,

mientras que ωo toma un valor finito y decrece como k 2 .
Si ahora analizamos que pasa en torno a k = π/a, utilizando las mismas aproximaciones
y un cambio de variables k 0 = k − π/a para centrarnos en el borde de la zona de Brillouin
obtenemos
r r
Ka2 Ga2
   
0 2G 02 0 2K 02
ωa (k  π/a) = 1− k ωo (k  π/a) = 1+ k
M 8(K − G) M 8(K − G)
Una es creciente y la otra decreciente, pero ambas en formas de cuadrática. Como la
función de k en las relaciones de dispersión es cos ka, que es monótona entre 0 y π/a,
podemos concluir qué puntos del centro se conectan con qué puntos del borde de la zona
(y ası́ es como elegimos cuál corresponde a cuál). Las relaciones de dispersión tienen la
siguiente forma (tomando a = 1, M = 1, G = 1, K = 2):

2.5

1.5
ω

0.5

0
0 π/2 π
a*k

Ahora vamos a analizar por qué pusimos los nombres acústico y óptico a las ramas.
Para eso, estudiemos el comportamiento en el centro de la zona de Brillouin k  π/a
que corresponde a las longitudes de onda largas (varias celdas primitivas). Como ya
conocemos las frecuencias, podemos calcular las amplitudes de oscilación. Y para hacer
más sencillo aún el problema, calculemos en k = 0. Volviendo a nuestro sistema de dos
ecuaciones lineales con dos incógnitas (reemplazando ω = 0 y k = 0), obtenemos que
para la rama acúsitca se cumple que A1 = A2 . Esto implica que todos los átomos se
mueven en la misma dirección con la misma amplitud. Lo que dice el modo k = 0 que
cuesta energı́a ω = 0 es que puedo mover tranquilamente todos los átomos de la cadena
en el mismo sentido y la misma cantidad, sin ningún costo energético (lógico, no?). Si me
alejo un poco del centro de la zona, ahora las amplitudes serán distintas, pero como la
longitud de onda es larga los átomos vecinos se desplazaran en fase. Como van oscilando
todos hacia un lado o hacia el otro, la energı́a que cuestan estos modos es baja, ya que

3
Materia Condensada Año 2020

los resortes no se comprimen ni estiran demasiado. Realmente, el modo k = 0 de energı́a

ω = 0 se conoce como un modo de Goldstone y aparece cuando se rompe una simetrı́a
continua en el problema original. En este caso la simetrı́a continua es la de traslación y
nosotros la rompimos cuando supusimos que los átomos estaban posicionados en na.
Ahora, cuando estamos en k = 0 y ω = 2(K+G) M
tenemos que A1 = −A2 . Esto nos dice
que en este modo, los átomos 1 y 2 se mueven en sentido opuesto. Esto cuesta mucha
energı́a, ya que los resortes se comprimen o estiran al máximo. El nombre de rama óptica
viene de que, si esta fuese una cadena de iones, los iones de distinta carga moviéndose
en direcciones opuestas actuarı́an como pequeños dipolos. Debido a que estos dipolos
interactúan con las ondas electromagnéticas, se da el nombre de rama óptica a esta rama.

Vayamos ahora a los casos particulares:

(a) M1 = M2 , G = K (d = a/2) implica que la cadena es realmente monoatómica y que

contiene 2N átomos idénticos conectados por resortes idénticos. La zona de Brillouin
de esta cadena irı́a desde −2π/a a 2π/a ya que la red de Bravais tiene constante
de red a/2. Pero como nosotros la estamos describiendo a partir de una red con
constante a y motivo de dos átomos, nuestra zona de Brillouin va de −π/a a π/a.
Sin embargo, tenemos que obtener la misma información de ambas maneras. Esto
hace que la relación de dispersión se doble para entrar entera en la zona de Brillouin
reducida. Para M = K = G = a = 1 tenemos:

2.5

1.5
ω

0.5

0
0 π/2 π 3π/2 2π
a*k

La rama violeta es la que falta para completar la rama verde más allá de π/a y
obtener la relación de dispersión de la cadena lineal monoatómica.

4
Materia Condensada Año 2020

(b) Cuando K  G cuesta mucha energı́a comprimir los resortes K y entonces las masas
conectadas por ellos actúan como una sola masa. Esto se ve reflejado en que la rama
ópica se va muy arriba en energı́a y se hace chata, mientras queda una sola rama
acústica a bajas energı́as. Cambiando ahora K = 100, obtenemos que la rama óptica
se fue de la escala y la relación de dispersión es la de una cadena lineal monoatómica
de constante de red a.

2.5

1.5
ω

0.5

0
0 π/2 π
a*k

(c) El caso M1 = M2 , G > K es el que ya analizamos antes, donde hay un gap de

energı́a entre ambas ramas.
(d) Si ahora tenemos G = K, M1 > M2 , tenemos que trabajar nuevamente las expre-
sionesy vamos a encontrar que la situación es cualitativamente igual a cuando G y
K son distintos pero M1 = M2 .

2. Problema:
Calcular la relación de dispersión y analizar los modos de vibración de una red cuadrada
bidimensional con interacciones como se muestra en la Figura (KL : constante de inter-
acción longitudinal; KT : constante de interacción transversal). Analizar en particular los
modos en las direcciones (q, 0) y (q, q).

5
Materia Condensada Año 2020

Solución:
Llamemos ui,j i,j
x al desplazamiento horizontal y uy al desplazamiento vertital (respecto del
equilibrio) del átomo ubicado en la posición (i, j) de la red. El potencial total será

N N h
KL X h i,j i+1,j 2 i,j

i KT X
i,j+1 2 i+1,j 2

i
i,j+1 2
ui,j i,j



U= ux − ux + uy − uy + y − u y + u x − u x
2 i,j=1 2 i,j=1

donde tuvimos en cuenta N 2 celdas primitivas y condiciones periódicas de contorno:

ui+N,j
α = ui,j
α , uαi,j+N = ui,j
α

Tendremos entonces una ecuación diferencial para cada coordenada (i, j) y componente
x, y:
.. i,j
M ux = KL ui+1,j + uxi−1,j − 2 ui,j


x x
+ KT uxi,j+1 + ui,j−1 − 2 ui,j


x x
.. i,j
M uy = KL ui,j+1 + ui,j−1 − 2 ui,j


y y y
+ KT ui+1,j + ui−1,j − 2 ui,j


y y y

Proponemos entonces soluciones de la forma

ui,j = εei(k·r )
i,j −ωt

donde ri,j = a(ix̂ + j ŷ). Reemplazando,

−εx M ω 2 = KL eikx a + e−ikx a − 2 εx

+ KT eiky a + e−iky a − 2 εx

−εy M ω 2 = KL eiky a + e−iky a − 2 εy

+ KT eikx a + e−ikx a − 2 εy

6
Materia Condensada Año 2020

donde vemos que el sistema de ecuaciones ya está diagonalizado, con los autovalores

M ω12 = KL (2 − 2 cos(kx a)) + KT (2 − 2 cos(ky a))

   
2 kx a 2 ky a
= 4KL sin + 4KT sin
2 2
2
M ω2 = KL (2 − 2 cos(ky a)) + KT (2 − 2 cos(kx a))
   
2 ky a 2 kx a
= 4KL sin + 4KT sin
2 2

donde cada una de estas relaciones de dispersión representa una rama. El vector k se
mueve en
la primer zona de Brillouin, que es el cuadrado cuyos vértices se encuentran en
± πa , ± πa .
Para k = (q, 0), tenemos
r
4KL q a
ω1 = sin
r M 2
4KT  q a
ω2 = sin
M 2

Mientras que para k = (q, q), tenemos

r
4KL + 4KT q a
ω12 = sin
r M 2
4KL + 4KT  q a
ω22 = sin
M 2

A la hora de graficar la relación de dispersión es conveniente elegir un camino a lo largo

de la zona de Brillouin bidimensional que pase por los puntos de mayor simetrı́a. En el
caso de la red cuadrada, un camino tı́pico es (0, 0) → ( πa , πa ) → ( πa , 0):

3.5

2.5

2
ω

1.5

0.5

0
(0,0) (π,π) (π,0) (0,0)
a*k

7
Materia Condensada Año 2020

3. Problema:
Considerar una cadena lineal monoatómica como la que se muestra en la Figura. Sobre
la primer masa actúa una fuerza consenoidal F (t) = F0 cos(ω0 t). Describir el movimiento
de los átomos de la cadena.

Solución:
Teniendo en cuenta sólo la interacción entre los átomos, el potencial total se escribe como,
N −1
KX
U= (un − un+1 )2
2 n=1

Las ecuaciones de movimiento seràn

..
m u1 = − K (u1 − u2 ) + F (t)
..
m un = − K (un − un+1 ) + K (un−1 − un ) , n = 2, 3, . . . , N
Para resolver las ecuaciones diferenciales, pongo la fuerza externa como una función com-
pleja,
..
m u1 = − K (u1 − u2 ) + F0 eiω0 t
..
m un = K (un+1 + un−1 − 2 un ) , n = 2, 3, . . . , N
Proponemos soluciones de la forma
un = ε ei (kna+ω0 t) (1)
y
−m ω02 ε = − Kε 1 − ei ka + F0 e−i ka


 
−i ka
−m ω02 ε = Kε e|i ka +
{ze } −2
 , n = 2, 3, . . . , N
2 cos(ka)

Reacomodando,
K − Kei ka − m ω02 ε = F0 e−i ka


(2)
2K [1 − cos(ka)] − m ω02 ε = 0 , n = 2, 3, . . . , N .


(3)

Recordemos que para el problema no forzado, visto en clases, la frecuencia al cuadrado

es ω 2 = 2K
m
(1 − cos(ka)), la cual toma un valor máximo 4K/m. Si ω02 toma algún valor
menor o igual a 4K/m, entonces F estará excitando algún modo normal, ya que habrá
algún k para el cual ω 2 será igual a ω02 .
4K
Separaremos ası́ en tres casos dieferentes: ω02 mayor, menor e igual a m
.

8
Materia Condensada Año 2020

• Caso: ω02 ≤ 4Km

La ecuación (3) se satisface con el k tal que
2K 4K
ω02 = (1 − cos(ka)) = sin2 (ka)
m m

Agregando esta expresión de ω0 a la ecuación (2),

F0 e−i ka = K − Kei ka − m ω02 ε

= K − Kei ka − 2K (1 − cos(ka)) ε

= K 1 − ei ka − 2 + 2 cos(ka) ε

= K (−1 + cos(ka) − i sin(ka)) ε

= K −1 + e−i ka ε

F0 = K 1 − ei ka ε

F0 1
ε= (1 − e−i ka )
K |1 − ei ka |2
F0 1
= [1 − cos(ka) + i sin(ka)]
K (1 − cos(ka))2 + sin2 (ka)
F0 1
= [1 − cos(ka) + i sin(ka)]
K 2 − 2 cos(ka)
Finalmente, para obtener la solución
q final, debemos tomar la parte real de un =
ω 2m
ε ei (kna+ω0 t) , siendo k = arc sin 4K0
y ε y ω0 las calculadas más arriba.
• Caso: ω02 > 4K
m
No podremos satisfacer la ecuación (3) con números k reales, entonces extendemos
al plano complejo: k → k1 + ik2 . Asi, la ecuacion (3) nos pide
K
ω02 = 2 − ei ka − e−i ka


m
K
2 − ei k1 a e−k2 a − e−i k1 a ek2 a


=
m
K
2 − [cos(k1 a) + i sin(k1 a)] e−k2 a − [cos(k1 a) − i sin(k1 a)] ek2 a


=
m
K
2 − cos(k1 a) e−k2 a + ek2 a − i sin(k1 a) e−k2 a − ek2 a




=
m
K
= (2 − 2 [ cos(k1 a) cosh(k2 a) − i sin(k1 a) sinh(k2 a) ] )
m

9
Materia Condensada Año 2020

igualando partes real e imaginaria

2K
(1 − cos(k1 a) cosh(k2 a)) = ω02
m
sin(k1 a) sinh(k2 a) = 0

De la segunda ecuación surgen dos opciones:

– Si k2 = 0
Tenemos entonces que
2K
ω02 = (1 − cos(k1 a))
m
y estamos en el caso anterior, donde ω02 ≤ 4K m
.
– Por lo tanto tiene que ser k2 6= 0.
Entonces tiene que ser k1 = 0 o k1 = ± πa , y será ω02 = 2K
m
(1 ∓ cosh(k2 a)).
Como cosh(k2 a) > 1, la solución correspondiente a k1 = 0, que tiene el signo
menos, no es fı́sica. Por lo tanto será k1 = ± πa , y

2K
ω02 = (1 + cosh(k2 a))
m  
4K 2 k2 a
= cosh
m 2

En la figura graficamos ω0 en función de k1 y k2 (la parte con k1 ≥ 0). Asi, cuando la

q
4K
fuerza externa tenga una frecuencia mayor a m
, la solución (1) estará modulada
−k2 na
por una función de decaimiento exponencial e , con

mω02
 
1
k2 = arc cosh −1
a 2K
"r #
2 mω02
= arc cosh
a 4K

10
Materia Condensada Año 2020

Las soluciones serán de la forma

un = ε e−k2 na ei nπ ei ω0 t = (−1)n e−k2 na ε ei ω0 t
donde tenemos que tomar la parte real. La perturbación que provoca la fuerza
externa sobre el sistema tendrá una longitud de penetración de λ = 1/k2 .

4. Problema:
En los semiconductores covalentes (C, Si, Ge) y también en los parcialmente iónicos
(GaAs), los electrones que forman los enlaces están parcialmente localizados entre los
átomos. La dinámica de la red de estos compuestos es adecuadamente descripta por
el modelo de “carga de enlace” (bond charge model), en el cuál se consideran cargas
electrónicas ubicadas entre los iones, acopladas a ellos y entre sı́. La masas de estas
cargas se desprecian frente a la de los iones, de modo que se reacomodan adiabáticamente
según los iones se desplazan. Este modelo describe una caracterı́stica común de dichos
compuestos: las ramas TA tienen frecuencias muy bajas y son muy chatas hacia el borde
de la primera zona de Brillouin, mientras que sus pendientes cerca de q = 0 son bastante
grandes.
La esencia de este modelo se ilustra en la cadena monoatómica de la Figura, donde se
ubica una carga de enlace entre cada par de iones.

Consideremos un acoplamiento armónico f entre la carga y cada ion, y otro f 0 entre las
cargas a cada lado de un ion. Sea m la masa iónica. La masa de la carga de enlace se
considera nula, de modo que ellas, en lugar de satisfacer una ecuación dinámica, obedecen
una condición de equilibrio (condición adiabática) para una configuración instantánea de
iones.
(a) Obtener las ecuaciones de movimiento y la relación de dispersión ω(q).
(b) Mostrar que para q → 0 se obtiene ω(q) = cq, donde c2 = (f + 2f 0 )/(2m).
(c) Mostrar que para valores crecientes de f 0 la rama se achata hacia el borde de la
primera zona de Brillouin, y que de esta manera, el modelo da cuenta de la carac-
terı́stica común que presentan las ramas TA en los semiconductores covalentes.
Solución:

(a) El primer paso es escribir la energı́a potencial del sistema. Tomemos como celda
unidad un conjunto ión-electrón, entonces
fX 2 2 f0 X
U= (u1,n − u2,n ) + (u1,n − u2,n−1 ) + (u2,n − u2,n+1 )
2 n 2 n

11
Materia Condensada Año 2020

donde u1,n es el desplazamiento del ion de la celda n a partir de su posición de

equilibrio, y u2,n lo mismo pero para el electrón. Derivando con respecto a ux,n para
x = 1, 2 podemos calcular la fuerza y obtener las ecuaciones de movimiento a partir
de la Ley de Newton

mü1,n = −f (u1,n − u2,n ) − f (u1,n − u2,n−1 )

me ü2,n = f (u1,n − u2,n ) + f (u1,n+1 − u2,n ) − f 0 (u2,n − u2,n+1 ) + f 0 (u2,n−1 − u2,n )

Proponiendo soluciones de la forma

u1,n = Aei(kna−ωt) u2,n = Bei(kna−ωt)

donde a es la distancia entre iones, podemos transformar el sistema de ecuaciones

diferenciales en un sistema de ecuaciones lineales. Aún más, vamos a tomar me = 0
ya que es despreciable en comparación a la masa del ion. Esto lleva al sistema

−f (1 + e−ika )
  
2f − mω 2 A
=0
−f (1 + eika ) 2f + 2f 0 (1 − cos ka) B

El determinante de la matriz debe anularse para que exista solución distinta a la

trivial. Esta condición lleva a la ecuación de la relación de dispersión
f 2 (1 + cos ka)
mω 2 = 2f −
f + f 0 (1 − cos ka)
Como cos ka es siempre menos a 1, el denominador no se anula nunca. Evaluemos
los extremos de la relación de dispersión: para k = 0 tengo cos ka = 1 y entonces
resulta ω = 0, como es de esperarse ppara la rama acústica; para k = π/a tengo que
cos ka = −1 y entonces resulta ω = 2f /m, que no depende de f 0 .
(b) Para encontrar la velocidad del sonido tenemos que estudiar el comportamiento de
ω(k) para k  1. Para simplificar los calculos podemos tomar a = 1 o tomar ka
como nuestra nueva variable k 0 (y llamarla k). De cualquier modo, me voy a olvidar
2
de la constante a sin pérdida de generalidad. Empecemos aproximando cos k ' 1− k2
para k chico. Entonces obtengo
k2
2 f 2 (2 − 2
) f 2 (4 − k 2 )
mω = 2f − 2 = 2f − = g(k)
f+ k
f0 2 2f + f 0 k 2

Ahora queremos realizar una expansión de Taylor a orden 2 de esta función

∂g 1 ∂ 2g
g(k) = g(0) + k+ k2
∂k k=0 2 ∂k 2 k=0

ya que esperamos que la relación de dispersión sea lineal en torno al centro de la

zona de Brillouin (por ende, ω 2 ∝ k 2 ). g(0) ya sabemos que es 0. Calculemos la dos
derivadas
∂g 2f 2 k(2f + f 0 k 2 ) − 2f 0 kf 2 (k 2 − 4)
=
∂k (2f + f 0 k 2 )2

12
Materia Condensada Año 2020

que evaluado en k = 0 da 0. Ahora hay que calcular la segunda derivada, lo cual

excede totalmente las ganas de cualquier persona a esta altura de la carrera. En
lugar de eso, sabiendo que queremos evaluar en k = 0, calculemos sólo los términos
que no se anularán.

∂ 2g 4f 3 (2f + f 0 (k = 0)2 )2 + 8f 0 f 2 (2f + f 0 (k = 0)2 )2 16f 5 + 32f 0 f 4

= =
∂k 2 k=0 (2f + f 0 (k = 0)2 )4 16f 4

que finalmente se escribe como (f + 2f 0 ). Finalmente obtenemos

f + 2f 0 2
ω2 = k
2m
que es exactamente lo que pedı́a el enunciado.
(c) Fijando todas las variables libres a 1, y variando f 0 obtenemos las siguientes curvas
donde se puede ver que la banda se va achatando a medida que f 0 crece. Como ya
habı́amos visto, el valor de la relación de dispersión en el borde de zona no depende
de f 0 .

1.6

1.2
ω

0.8

f'=0
0.4 f'=1
f'=10
f'=100
0
0 0.5π/a π/a
k

5. Problema:

(a) Calcular la densidad de modos de vibración g(ω) para una red lineal monoatómica
de N átomos. ¿Por qué se observan divergencias (singularidades de van Hove)?
(b) A partir de la relación de dispersión de una cadena diatómica como la que se muestra
en la Figura, hacer un dibujo cualitativo de g(ω).

Solución:

(a) Para calcular la densidad de modos de vibración g(ω) tenemos que usar la relación

dk
g(ω) dω =
2π

13
Materia Condensada Año 2020

Relación de dispersión de la cadena diatómica. Ejercicio 5-b.

Reescribiendo la expresión,
1 1
g(ω) =
2π dω
dk
vemos que la densidad de modos es inversamente proporcional a la pendiente de la
relación de dispersión. Cómo la relación de dispersión es simétrica con respecto a k,
la densidad de modos total será igual al doble de la densidad de modos de la rama
correspondiente a k > 0. De modo que vamos a calcular
1 1
g(ω) =
π dω
dk

para k > 0. La derivada de la relación de dispersión para la cadena lineal monoatómica

es r  ! r  
dω d K 1 K 1
= 2 sin ka = a cos ka
dk dk M 2 M 2
Ahora tenemos que expresar esta derivada en función de ω, para poder tener la
densidad de modos en función de ω y no de k. De la relación de dispersión sale
r !
2 ω M
k = arcsin
a 2 K
Y reemplazando obtenemos
r r ! ! r r !!
dω K 12 ω M K ω M
= a cos arcsin a = a cos arcsin
dk M 2a 2 K M 2 K
p
Ahora, reescribiendo cos x = 1 − sin2 x llegamos a
v
r u r !2 r
dω Ku t1 − ω M K ω2
=a =a −
dk M 2 K M 4

14
Materia Condensada Año 2020

Finalmente llegamos a
1 1
g(ω) = q
aπ K − ω2
M 4

Esta función toma un valor

p finito cuando ω = 0 y luego es monótona creciente hasta
diverger cuando ω = 2 K/M , correspondiente a los k en el borde de la zona de
Brillouin. Poniendo K = M = a = 1, obtenemos

1.5
g(ω)

0.5

0
0 0.5 1 1.5 2 2.5 3
ω

(b) Para tener una idea cualitativa de la densidad de modos g(ω) de una dada relación
de dispersión hay que tener en cuenta la inversa de las pendientes de la relación de
dispersión. La rama acústica proveerá una curva igual a la del ı́tem anterior. Luego
habrá un gap de energı́a para el cuál no existirá ninguna densidad de modos, y luego
vendrá la densidad asociada a la rama óptica. Como la pendiente de la rama óptica
se anula en el lı́mite inferior y superior, la densidad de modos tendrá dos divergencias
en ambos extremos y un mı́nimo en el medio correspondiente a la energı́a a la cuál
se da la mayor pendiente de la relación de dispersión. Algo ası́:

6. Falta
7. Problema:
Calcular la densidad de modos de vibración y la dependencia con la temperatura del calor
especı́fico de un sólido de Debye en 1, 2, y 3 dimensiones. Repetir lo mismo utilizando el
modelo de Einstein.

15
Materia Condensada Año 2020

Solución:
Existen diferentes formas para calcular la densidad de estados g(ω). Repasemos (para el
caso monoatómico):
Cuando el volumen de un cristal V , donde hay N celdas primitivas, tiende a infinito la
sumatoria en los vectores ~k se puede transformar en una integral sobre la primer ZB:
d Z
1 X X d~k
−→
V s=1 1ZB
(2π)d
~k,s

donde d es la dimensionalidad del cristal y s enumera las diferentes ramas de dispersión

para cada dimensionalidad.
La densidad de estados se obtiene de resolver que
d~k
Z ∞ Z
N
n= = g(ω)dω = d
(4)
V 0 1ZB (2π)

Para relaciones de dispersión isotrópicas (que no dependen de la dirección de ~k) podemos

utilizar coordenadas esféricas/cilı́ndricas:
3D : d~k = k 2 dkdΩ
2D : d~k = k dkdθ
1D : d~k = 2 dk
y será
k2 1
3D : g(ω) = 2
2π dω/dk
k 1
2D : g(ω) = (5)
2π dω/dk
1 1
1D : g(ω) =
π dω/dk
donde k debe expresarse como función de ω.
Otra forma de calcular g(ω) para los casos isotrópicos (y que nos sirve para el modelo de
Debye) es pensando en la cantidad de estados que encierra una esfera de radio k:
Vol. esfera radio k
N (k) =
Vol. que ocupa un sólo estado
El volumen que ocupa un sólo estado es (2π)d /V , y
k3
3D : N (k)/V =
6π 2
k2 1
2D : N (k)/V = (6)
4π dω/dk
k
1D : N (k)/V =
π

16
Materia Condensada Año 2020

Cada una de las ecuaciones previas es un n(k (ω)), el cual puede derivarse para obtener la
densidad de estados
dn(k (ω)) ∂n(k (ω)) dk
g(ω) = =
dω ∂k dω
que reproduce las ecuaciones (5) encontradas antes.
El calor especı́fico por unidad de volumen se obtiene de derivar la energı́a interna respecto
de la temperatura:

∂ 1 X ~ωs (~k)
cv = (7)
∂T V
~k,s
e~ωs (~k)β − 1
 2
~ω ( ~k)β e~ωs (~k)β
1 X s
= kB 2 (8)
V

(~k)β
~k,s e ~ω s −1
 2
~ ~ωs (~k)β
d~k ~ωs (k)β e
X Z
= kB 2 (9)
(2π)d

s e~ωs (~k)β − 1
XZ ∞ e~ωβ
= kB dω g(ω) (~ωβ)2 (10)
s 0 (e~ωβ − 1)2

Modelo de Debye.
Reemplazo las dispersiones ωs (~k) por d ramas con la misma dispersión isotrópica lineal,
ω(k) = c k. El kD y ωD = c kD de Debye, son hasta donde deben hacerse las integrales y
se obtienen haciendo en las ecuaciones (6) que N (kD ) = N (o n(kD ) = n), y entonces
3
3D : kD = 6 n π2
2
2D : kD = 4nπ
1D : kD = n π

Viendo también que dω/dk = c = ωD /kD , obtenemos finalmente:

ω2
3D : gD (ω) = 3 n 3
ωD
ω
2D : gD (ω) = 2 n 2 (11)
ωD
n
1D : gD (ω) =
ωD

Para calcular el calor especı́fico reemplazo estas densidades de estados en la ecuación (10),

17
Materia Condensada Año 2020

y llamando TD = ~ ωD /kB (en los libros aparece como ΘD ),

3 Z TD
x4 ex

T T
3D : cv (ω) = 9 kB n dx
TD 0 (ex − 1)2
2 Z TD
x3 ex

T T
2D : cv (ω) = 4 kB n dx
TD 0 (ex − 1)2
TD
x2 ex
 Z
T T
1D : cv (ω) = kB n dx
TD 0 (ex − 1)2

Modelo de Einstein.
Reemplazo las dispersiones ωs (~k) por d ramas con la misma frecuencia ωE , es decir, sin
dispersión.
La densidad de estados es igual para cualquier dimensión: g(ω) = nδ(ω − ωE ). Lo mismo
sucede con el calor especı́fico,
 2 TE
TE e T
cv = d n k B 2
T
 TE
e T −1

cuyos lı́mites son

cv ∝ e−Te /T para T  TE
cv = d n k B para T  TE

8. Problema:
Calcular el número medio de fonones como función de la temperatura en un cristal
monoatómico construido por NA átomos, utilizando:

(a) El modelo de Debye.

(b) El modelo de Einstein

Analizar en particular el comportamiento para altas y bajas temperaturas. Explicar el

porqué de semejanzas y diferencias.
Solución:
La densidad del número medio de fonones, es decir el número medio de fonones por unidad
de volumen, se escribirá
XZ ∞ 1
ph
n = g(ω) ~ωβ dω
s 0 e −1
Utilizando las densidades de estados del ejercicio anterior tenemos que, para 3D
(a) Modelo de Debye
3 Z TD
x2

ph T T
n = 9n dx
TD 0 ex − 1

18
Materia Condensada Año 2020

A bajas temperaturas el extremo superior de la integral tiende a infinito, y por ende la

integral tiende a un número. Entonces nph ∝ T 3 para bajas temperaturas.
x2
A altas temperaturas podemos expandir el integrando para valor pequeños de x, ex −1
≈
2 3
x − x2 + x12 − . . . .
Entonces, calculando el primer término
 3 Z TD

ph T T
n = 9n xdx + . . .
TD 0
 3 2 TD
T x T
= 9n + ...
TD 2 0
 
T
∝ 9n
TD

(b) Modelo de Einstein

Para 3D es
1 1
nph = 3 n = 3n TE
e~ωE β −1 e T −1
TE
A bajas temperaturas, tiende a cero como e− T , y a altas temperaturas crece linealmente
con T .
De manera similar puede verse el comportamiento para 2 y 3 dimensiones.

9. Problema:
La densidad de estados g(ω) y el calor especı́fico de un sólido de Debye son la siguientes

El tungsteno y el litio tienen las siguientes densidades de estado, dibujadas sobre el fondo
de la de Debye

19
Materia Condensada Año 2020

Las integrales de las tres densidades mostradas son iguales, como también lo son las
velocidades del sonido en los tres sólidos. Dibuje equemáticamente los calores especı́ficos
del W y del Li sobre la curva del de Debye. Explique el porqué de las semejanzas y
diferencias.
Solución:
Debido a que ambas densidades de estados, la del W y la del Li, son iguales a bajas energı́as
(misma velocidad del sonido) las respectivas curvas del Cv a bajas temperaturas serán
parecidas. Sin embargo, las diferencias que aparecen a altas energı́as entre la densidad de
estados de ambos compuestos hará que alcancen a distintas temperaturas el lı́mite clásico.
El lı́mite clásico será igual para ambos compuestos, ya que la integral de la densidad de
estado es igual en ambos casos. Como el W tiene la densidad de estado concentrada a
más bajas energı́as, llegará a más bajas temperaturas al lı́mite clásico. Sucederá al revés
con el Li, que tiene un mayor peso de la densidad de estados a altas energı́as.

10. Falta

11. Problema:
La siguiente Figura representa las interacciones de un material monoatómico unidimen-
sional depositado sobre un sustrato rı́gido.

Los resortes actúan sólo ante desplazamientos longitudinales (esto es, cuando los desplaza-
mientos son en la misma dirección en que está dibujado el resorte).

(a) Calcule la relación de dispersión de los diferentes modos. Grafı́quela.

(b) Calcule y grafique la densidad de estados asociada g(ω).
(c) Calcule la dependencia con la temperatura del calor especı́fico en el lı́mite de bajas
temperaturas.
(d) Suponga que este material es depositado sobre dos sustratos diferentes, de forma tal
que la interacción sustrado-material es en un caso muy fuerte y en el otro muy débil.
Es decir, para un mismo valor de K se tienen los siguientes casos:
Sustrato 1: G1  K.

20
Materia Condensada Año 2020

Sustrato 2: G2  K.
Dibuje esquemáticamente la dependencia con la temperatura del calor especı́fico
para estos dos casos. Explique el porqué de semejanzas y diferencias.

Solución:

(a) La energı́a potencial del sistema es

KX GX
U= (ux,n − ux,n+1 )2 + (uy,n )2
2 n 2 n

donde ux,n es el desplazamiento del átomo n-ésimo en la dirección x a partir de su

posición de equilibrio. Las ecuaciones de movimiento resultan

müx,n = −K (ux,n − ux,n+1 ) + K (ux,n−1 − ux,n )

müy,n = −Guy,n
Como los resortes actúan solo longitudinalmente, ambas ecuaciones diferenciales
están desacolpadas. En la dirección x tenemos la ecuación diferencial de la ca-
dena lineal monoatómica, y en la dirección y tenemos un conjunto de osciladores
armónicos simples desacoplados entre sı́. Entonces una banda es la de la cadena
lineal monoatómica
2K
ωx2 = (1 − cos kx a)
m
y la otra es una banda chata correspondiente a los osciladores desacoplados
G
ωy2 =
m

(b) La densidad de estados seráp

la de una cadena lineal monoatómica sumada a una delta
de Dirac centrada en ω = G/m, que dependiendo del valor de G con respecto a
K puede, o no, estar superpuesta. Entonces
1 1 1  p 
g(ω) = q + δ ω − G/m
aπ K − ω2 a
m 4

21
Materia Condensada Año 2020

(c) Ahora calculemos el calor especı́fico

Z ωM
∂E ∂ ω
Cv = = g(ω)
∂T ∂T 0 eω/T − 1
donde estoy tomando ~ = kB = 1 ya que sólo nos importa la dependencia con T
a bajas temperaturas y ωM es el valor máximo que toma la relación de dispersión.
Reemplazando la relación de dispersión obtenemos
r
∂ G 1
Cv = Cv1D + √G 1
∂T m e m T − 1

donde el primer término es el Cv de una cadena lineal monoatómica que, a bajas

temperaturas, se comportará como un sólido de Debye unidimensional; y el segundo
no es otra cosa que el Cv de Einstein.
En 1D, el Cv de Debye a bajas temperaturas se calcula aproximando ω = ck. En-
1
tonces g(ω) = πc y el calor especı́fico es
ωD ωD ωD /T
ω2 eω/T x2 e x
Z Z Z
1 ∂n 1 T
Cv = ω dω = 2 dω = dx
0 cπ ∂T cπ 0
2
T (e ω/T − 1) cπ 0 (ex − 1)2
proporcional a T , lineal. Mientras que la dependencia del Cv de Einstein con la
temperatura no depende de la dimensión, y es
√
G/m
−
Cv ∝ e T

(d) Dependiendo de si la banda chata está metida dentro de la relación de dispersión

de la cadena lineal monoatómica, o no, el calor especı́fico cambiará cualitativamente
debido a que los estados de la banda chata serán accesibles a más bajas (o más altas)
temperaturas). Esquemáticamente:

Cuando G  K, el Cv llegará al lı́mite clásico a más bajas temperaturas debido a

que los estados de la banda chata están a energı́as bajas, del orden de las energı́as
de la rama acústica. Por otro lado, cuando G  K, la banda chata está más alta
en energı́a que la rama acústica, por lo que el lı́mite clásico se alcanzará a mayores
temperaturas.

12. Falta

22
Materia Condensada Año 2020

13. Falta

14. Falta

15. Falta

23