Capı́tulo 10

Teorı́a clásica del cristal armónico

Es una suposición artificial considerar que los iones se encuentren fijos en los sitios
de la red de Bravais
~ = n1~a1 + n2~a2 + n3~a3
R (10.1)
con ni entero y ~ai vectores primitivos de la red de Bravais.
En realidad los iones oscilan alrededor de los puntos de la red de Bravais (posiciones
de equilibrio), entonces, la posición instantánea del ión viene dado por
~ = R
~r(R) ~ + ~u(R)
~ (10.2)
~ es el desplazamiento desde su posición de equilinrio del ión en en sitio R.
donde ~u(R) ~
Asuminos que un par de iones separados una distancia ~r contribuyen a la energı́a
potencial con φ(~r) (por ejemplo, φ(~r) puede ser un potencial tipo Lennard-Jones).
Si los iones fueran fijos, entonces
1 X ~ ~0
U= φ(R − R ) (10.3)
2 ~ ~0
R,R

Haciendo un cambio de variable R ~ −R

~ 00 = R ~ 0 , se tiene
N X ~
U= φ(R) (10.4)
2 ~
R6=0

donde N es el número de sitios de la red de Bravais.

~ 6= R,
En general ~r(R) ~ la red no es estática, entonces la energı́a potencial de del cristal
viene dado por
1 X h ~ ~ 0)
i
U= φ ~r(R) − ~r(R (10.5)
2 ~ ~0
R,R

130
Remplazando
1 X h~ ~0 ~ − ~u(R
~ 0)
i
U= φ R − R + ~u(R) (10.6)
2 ~ ~0
R,R

El hamiltoniano del sistema del sistema viene dado por

X ~
P~ 2 (R)
H= +U (10.7)
~
2M
R

donde P~ es el momentum y M es la masa de los átomos respectivamente.

10.1. Aproximación armónica

Usar el potencial de pares φ(~r) es difı́cil. Realizamos una aproximación para ~u(~r)
pequeños (pequeños despazamientos desde la posición de equilibrio).
Hacemos un desarrollo en series de Taylor de U alrededor de su posición de equilibrio
~
R.
El teorema de Taylor en tres dimensiones para una función f (~r)
1 1
f (~r + ~a) = f (~r) + ~a · ∇f (~r) + (~a · ∇)2 f (~r) + (~a · ∇)3 f (~r) + · · · (10.8)
2 3!
Haciendo
~ −R
~r = R ~0 (10.9)
~ − ~u(R
~a = ~u(R) ~ 0) (10.10)

se tiene
N X ~ 1 X nh ~ i o
~ 0 ) · ∇ φ(R
~ −R
~ 0)
U = φ(R) + ~u(R) − ~u(R
2 ~ 2 ~ ~0
R6=0 R,R
1 X nh i o2
+ ~ − ~u(R
~u(R) ~ 0) · ∇ ~ −R
φ(R ~ 0 ) + 0(u3 ) (10.11)
4 ~ R~0
R,

El primer término es la energı́a potencial en el equilibrio, es constante independiente

~ y P~ (R)
de las variables dinámica ~u(R) ~ y puede ser ignorado.
Veamos el término lineal
X X X
~ · ∇φ(R
~u(R) ~ −R
~ 0) = ~ ·
~u(R) ~ −R
∇φ(R ~ 0) (10.12)
~ R
R, ~0 ~
R ~0
R

131
donde
X
~ −R
∇φ(R ~ 0) = 0 (10.13)
~0
R

~ de todos los otros átomos cuando está es su

es la fuerza, sobre el átomo en el sitio R,
posición de equilibrio. Por lo tanto, el término lineal se anula.
La primera corrección es dada por el término cuadrático y correponde a la energı́a
potencial armónica
1 X nh ~ i o
~ 0 ) · ∇ 2 φ(R
~ −R
~ 0)
U armo = ~u(R) − ~u(R (10.14)
4 ~ ~0
R,R

1 X h i h i
U armo = ~ − uµ ( R
uµ (R) ~ 0 ) φµ,ν (R
~ −R
~ 0 ) uν (R)
~ − uν ( R
~ 0) (10.15)
4 ~ R~ 0 ,µ,ν
R,

donde
∂ 2 φ(~r)
φµ,ν (~r) = (10.16)
∂rµ ∂rν

La aproximación armónica es el punto de partida para toda teorı́a dinámica de la red.

También se puede escribir el potencial armónico de manera diferente
1 X
U armo = ~ µ,ν (R
uµ (R)D ~ −R
~ 0 )uν (R
~ 0) (10.17)
2 ~ R~ 0 ,µ,ν
R,

donde
X
~ −R
Dµ,ν (R ~ 0) = δ ~ ~ 0 ~ −R
φµ,ν (R ~ 00 ) − φµ,ν (R
~ −R
~ 0) (10.18)
R,R
~ 00
R

10.2. Calor especı́fico de un cristal clásico

No podemos evaluar las propiedades del equilibrio, donde cada ión está fijo en cada
punto R~ de la red de Bravais. Pero una red no es estática, los iones oscilan alrededor
de su posición de equilibrio. Para calcular las propiedades debemos promediar sobre
todas las posibles configuaciones de los iones dando a cada configuración o estado un
peso

e−βE (10.19)

132
donde E es la energı́a de la configuración y β = 1/kB T .
Si hacemos un tratamiento clásico al cristal, de la mecánica estadı́stica se tiene
R
1 dΓe−βH H
u = R
V dΓe−βH
Z
1 ∂
= − ln dΓe−βH (10.20)
V ∂β
donde dΓ es el elemento de volumen en el espacio de fase.
Y Y
dΓ = ~ P~ (R)
d~u(R)d ~ = ~
duµ (R)dP ~
µ (R) (10.21)
~
R ~
R,ν

En la aproximación armónica

X ~
P~ 2 (R)
H= + U equ + U armo (10.22)
~
2M
R

Remplazando
Z Z P ~ 2 (R)
~

P
−βH −β ~
R 2M
+U equ +U armo
dΓe = dΓe
Z P ~ 2 (R)
~

P armo
−βU equ −β ~ 2M +U
R
= e dΓe
Z P ~ 2 (R)
~ P 
P ~ µ,ν (R−
~ R~ 0 )uν (R
~ 0)
−βU equ −β ~
R 2M
+ 12 ~ R
R, ~ 0 ,µ,ν uµ (R)D
= e dΓe
Z P P 
[β 1/2 P ~ 2
~ (R)]
~ ~ R~ 0 )[β 1/2 uν (R
~ 0 )]
−βU equ − ~
R 2M
+ 12 ~ R
R, ~ 0 ,µ,ν [β
1/2 u
µ (R)]D µ,ν (R−
= e dΓe

Haciendo un cambio de variable

~ = β 1/2 ~u(R)
~u0 (R) ~ ~ = β −1/2 ~u0 (R)
~u(R) ~ (10.23)
~ = β 1/2 P~ (R)
P~ 0 (R) ~ P~ (R) ~
~ = β −1/2 P~ 0 (R) (10.24)

Los elementos de volumen (ejemplo d~u(R) ~

~ = dux (R)du ~ ~
y (R)duz (R)) vienen dado por

~ = β −3/2 d~u0 (R)

d~u(R) ~ (10.25)
dP~ (R)
~ = β −3/2 dP~ 0 (R)
~ (10.26)

Z Z P ~ 02 (R)
~ P 
P ~ µ,ν (R−
~ R~ 0 )u0ν (R
~ 0)
−βH −βU equ −3N − ~
R 2M
+ 12 ~ R
R, ~ 0 ,µ,ν u0µ (R)D
dΓe = e β dΓe (10.27)

133
La integral es independiente de β (o de T ).
Z
equ
dΓe−βH = e−βU β −3N × Cte (10.28)

Remplazando
1 ∂  equ

u = − ln e−βU β −3N × Cte
V ∂β
1 ∂
= − ln (−βU equ − 3N ln β + ln Cte)
V ∂β
1 equ 1 1
= U + 3N (10.29)
V V β

u = uequ + 3nkB T (10.30)

con n = N/V , número de átomos por unidad de volumen.

En T = 0 se tiene u = uequ , corresponde a la teorı́a de la red estática (energı́a del punto
cero, se ignora). En T 6= 0 la energı́a u = uequ es corregida en 3nkB T pequeño, puesto
que kB T es del orden de unos pocos centécimos de 1eV , a temperatura ambiente.
El calor especı́fico, que es mucho más fácil de medir que u, viene dado por
!
∂u
cv =
∂T v
= 3nkB (10.31)

La red estática no contribuye al calor especı́fico. El calor especı́fico por ión, de un

aislante, debido a las vibraciones de la red es 3kB (Ley de Dulong y Petit).
El calor especı́fico molar en un sólido monoatómico viene dado
N N/NA N o moles 3NA kB
cv = 3nkB = 3 kB = 3 kB = 3 NA kB = mol (10.32)
V V /NA V v

donde NA es el número de Avogadro y v mol = V /N o demoles es el volumen molar.

cmol = v mol cv = 3NA kB cmol = 5, 96cal/mol − K (10.33)

El punto 2 puede ser explicado sobre la base puramente clásica como un fracaso de la
aproximación armónica. De acuerdo con la teorı́a clásica, a muy bajas temperaturas
las energı́as térmicas son insufucientes para permitir que el ión se desplace distancias
apreciables desde su posición de equilibrio, por lo tanto la aproximación armónica se
hace buena cuando la temperatura cae. Sin embargo, a temperaturas altas los iones

134
poseen suficiente energı́a para que el ión se mueva lejos de su posición de equilibrio,
tal que no se pueden despreciar los términos anarmónicos ya se vuelven importantes.
Ası́, la mecánica estadı́stica clásica indica que la ley de Dulong y Petit no necesita
ser perfectamente obedecida a temperaturas altas, pero debe ...
La estructura del hamiltoniano del oscilador armónio es cuadrático en ~u(R) ~ y P~ (R)
~
entonces puede ser resuelta exactamente, para pequeñas oscilaciones. En movimiento
general de N iones es la superposición (combinación lineal) de 3N modos normales
de vibración de cada uno con sus propias frecuencias caracterı́sticas.
Los resultados de la teorı́a cuántica, la energı́a permitidas del oscilador armónico es
igual a
 
1
n+ hν (10.34)
2
La generalización para 3N osciladores independientes es obvio. Las energı́as para el
caso de 3N

10.3. Modos normales de una red de Bravais uni-

dimensional monoatómica
Consideremos un conjunto de iones de masa M distribuidos en una red unidimen-
sional.
Sea u(na) es desplazamiento del ión desde su posición de equilibrio, oscilando alrede-
dor de R = na. Si asumimos que solamente iones vecinos interactúan. Entonces, de
() la energı́a potencial armónica es
1 X
U armo = [uµ (R) − uµ (R0 )] φµ,ν (R − R0 ) [uν (R) − uν (R0 )]
4 R,R0 ,µ=x,ν=x
1 X
= [ux (R) − ux (R0 )] φxx (R − R0 ) [ux (R) − ux (R0 )] (10.35)
4 R,R0

Obviando el subı́ndice x y considerando solamente interacciones entre primeros veci-

nos
1X
U armo = [u(R) − u(R + a)] φxx (R − R − a) [u(R) − u(R + a)]
2 R
1X
U armo = [u(R) − u(R + a)] φxx (a) [u(R) − u(R + a)] (10.36)
2 R
Haciendo
φxx (a) = φxx (−a) = K (10.37)

135
 1 X
U armo = K [u(R) − u(R + a)]2 (10.38)
2 R

si R = na
1 X
U armo = K [u(na) − u[(n + 1)a]]2 (10.39)
2 n

Las ecuaciones de movimiento son

∂U armo
M ü = − (10.40)
∂u(na)
1 ∂ n o
= − K · · · + [u((n − 1)a) − u(na)]2 + [u(na) − u((n + 1)a)]2 (10.41)
+···
2 ∂u(na)
Entonces
M ü = −K [2u(na) − u((n − 1)a) + u((n + 1)a)] (10.42)
Las ecuaciones deben deben comportarce como si cada ión estuviese conectado a sus
dos vecinos por un resorte sin masa de constate elástica K.
Si la cadena tiene un número finito de N iones, ¿cómo llevar en cuenta los iones en
los dos extremos? Los iones en los extremos solamente interactúan con sus vecinos
internos, esto complica el análisis (sin alterar el resultado final). Si N es grande y no
estamos interesados en los efectos extremos, entonces aplicamos las condiciones de
frontera de Born-von Karman, unimos los extremos por el mismo resorte que conecta
dos inones internos.
La solución de la ecuación de movimiento es de la forma
u(na, t) ∝ ei(kna−ωt) (10.43)
De la condición de frontera peródica
u(na, t) = u(0, t) (10.44)
Entonces
eikN a = 1 (10.45)
2π n
kNa = 2πn → k= (10.46)
a N
Si consideramos dos vectores permitidos que se diferencian por un vector de la red
recı́proca k 0 = k + 2π/a, representan la misma solución
2π
u(na, t) ∝ ei[(k+ a )na−ωt]
∝ ei2πn ei(kna−ωt)
∝ ei(kna−ωt) (10.47)

136
 h i
Por tanto, hay N valores de k que dan soluciones distintas en en intervalo − πa , πa .
Remplazando
h i
−Mω 2 ei(kna−ωt) = −K 2 − e−ika − eika ei(kna−ωt) (10.48)

h  i
−Mω 2 = −K 2 − eika + e−ika
= −2K [1 − cos ka]
 
2 1
= −4K sen ka (10.49)
2
Entonces
s 
 
K  1 

ω = 2 sen ka  (10.50)
M 2
que depende del momentum de la onda. Esta frecuencia y su energı́a asociada h̄ωk

Figura 10.1: .

es también periódica, con la zona de Brillouin cubriendo un periodo.

Estas soluciones simples, ondas planas oscilantes, son llamadas “modos normales”. Se
puede demostrar que cualquier solución general de la segunda ley de Newton puede

137
escribirse como una suma de tales modos, cada una con su propia amplitud relativa
y fase:
X
un (t) = αk eikna cos (ωk t + φk ) (10.51)
k

donde 2N parámetros αk y φk

10.4. Modos normales de una red de Bravais uni-

dimensional con base
Sea una red de Bravais unidimensional con dos iones por celda primitiva con
posiciones de equilibrio en na y na + d. Consideremos el caso en que los iones son
idénticos y que d ≤ a/2. La fuerza entre vecinos depende de si su posición es d o a − d.
Asumiendo interacciones entre primeros vecinos.
1 X 1 X
U armo = K [u1 (na) − u2 (na)]2 + G [u2 (na) − u1 ((n + 1)a)]2(10.52)
2 n 2 n

donde K y G son los coeficientes de elasticidad, asumimos K ≥ G y d ≤ a/2.

Las ecuaciones de movimiento viene dado por
∂U armo
M ü1 = − (10.53)
∂u1 (na)
1 ∂ n o
= − K · · · + [u1 (na) − u2 (na)]2 + · · ·
2 ∂u1 (na)
1 ∂ n o
− G · · · + [u2 ((n − 1)a) − u1 (na)]2 + · · ·
2 ∂u1 (na)
= −K [u1 (na) − u2 (na)] − G [u2 ((n − 1)a) − u1 (na)] (−1)
= −K [u1 (na) − u2 (na)] − G [u1 (na) − u2 ((n − 1)a)] (10.54)

Análogamente
∂U armo
M ü1 = − (10.55)
∂u1 (na)
= −K [u2 (na) − u1 (na)] − G [u2 (na) − u1 ((n − 1)a)] (10.56)

Cuyas soluciones son de la forma

u1 (na, t) = 1 ei(kna−ωt) (10.57)

u2 (na, t) = 2 ei(kna−ωt) (10.58)

138
donde 1 y 2 son constantes a determinar y 1 /2 es la amplitud de fase relativa de
vibración de los iones dentro de cada celda primitiva. De las condiciones de frontera
periódicas, se encuentran que existen N valores de k no equivalentes dentro de la
celda primitiva de la red recı́proca.
Remplazando en las ecuaciones de movimiento
 
−Mω 2 1 = −K(1 − 2 ) − G 1 − 2 e−ika
despejando
h i  
Mω 2 − (K + G) 1 + K + Ge−ika 2 = 0 (10.59)
Análogamente
 
−Mω 2 2 = −K(2 − 1 ) − G 2 − 1 eika
despejando
  h i
K + Geika 1 + Mω 2 − (K + G) 2 = 0 (10.60)
El sistema de ecuaciones tiene solución si
 
 Mω 2 − (K + G) K + Ge−ika 
 
 ika  =0 (10.61)
 K + Ge Mω 2 − (K + G) 

Resolviendo
h i2     2
Mω 2 − (K + G) = K + Geika K + Ge−ika = K + Geika 
 
= K 2 + G2 + KG eika + e−ika
= K 2 + G2 + 2KG cos ka (10.62)
Entonces
(K + G) 1 h 2 i1/2
ω2 = ± K + G2 + 2KG cos ka (10.63)
M M
s   
K  1 
ω = 2 sen ka  (10.64)

M 2
La amplitud de fase relativa
2 K + Geika K + Geika
= − = −
1 Mω 2 − (K + G) ± |K + Geika |
K + Geika
= ∓ (10.65)
|K + Geika |
Para cada uno de los N valores de k hay dos soluciones, 2N modos normales, 2N
grados de libertad (ya que se tiene dos iones en cada celda primitiva).

139
 Figura 10.2: .

10.5. Modos normales de una red de Bravais uni-

dimensional con base (b)
Sea una cadena lineal con dos tipos de átomos con masas M y M 0 alternandose
regularmente. Tomamos una cadena con 2N átomos, N celdas unitarias cada una de
longitud 2a, y condiciones de frontera periódicas. Cada tipo de átomo esta acoplado
por resortes idéticos a sus vecinos. La energı́a potencial es
N h
X i
U= K(un − vn )2 + K(vn − un+1 )2 (10.66)
n=1

Las ecuaciones de movimiento son

∂ 2 un
M = K(vn + vn−1 − 2un ) (10.67)
∂2
∂ 2 vn
M 0 2 = K(un+1 + un − 2vn ) (10.68)
∂
Nuevamente esperamos que cada átomo semejante tenga ma misma amplitud pero
diferentes tipos de átomos tienen diferentes amplitudes relativas. Asumimos que la
forma en que varian de celda a celda semejante a exp(iqx), llamemos
un = Aei(2nqa−ωt) vn = A0 ei((2n+1)qa−ωt) (10.69)

140
Esta es en efecto una solución de () si
h   i
−Mω 2 A = K A0 einqa + e−inqa − 2A (10.70)
h   i
−M 0 ω 2 A0 = K A e−inqa + einqa − 2A0 (10.71)

es decir, si
  " 2 #1/2
2 1 1 1 1 4sen2 (qa)
ω (q) = K + 0 ±K + 0 − (10.72)
M M M M MM 0
y

A 2K cos(qa) 2K − M 0 ω 2
= = (10.73)
A0 2K − Mω 2 2K cos(qa)

La forma de ω(q) es graficada en la figura. Hay dos soluciones separadas que se

repiten cada π/a. Esto se espera puesto que la celda unitaria tiene una longitud de
2a y la primera zona de Brillouin va de −π/(2a) a π/(2a). Hay N valores de q en
este intervalo y dos ramas son necesarias para proporcionar los 2N modos normales
correspondientes a los 2N grados de libertad de 2N átomos confinados en una linea.
Los dos modos de mismo q tienen diferentes movimientos relativos de los dos átomos
en la celda unitaria. Para pequeños q el signo negativo en () da
 1/2
2K A
ω(q) = K qa; =1 (10.74)
M + M0 A0
y el signo positivo da
" #1/2 !
2K(M + M 0 ) MM 0 q 2 a2 A M0
ω(q) = 1− ; = − (10.75)
MM 0 2(M + M 0 )2 A0 M

La rama de frecuencias bajas tiene la caracterı́stica tı́pica de las ondas de sonido, y

para pequeños q los detalles de la estructura de la cadena no tiene efecto sobre los
modos. Por esta razón la rama más baja es llamada rama acústica. En el modo
superior los dos átomos en la celda se mueven en direcciones opuestas (A/A0 es neg-
ativo en (), ver también figura). Si los dos átomos tienen cargas opuestas como en
un cristal iónico simple tal modo debe crear un conjunto de dipolos oscilantes. Estos
deben estar fuertemente acoplados a un campo eléctrico y debe ser fijado a moverse
por un campo de alta frecuencia. Para muchos sólidos esta frecuencia es similar al
de la luz infraroja y por esta razon esta ramma es llamada rama acústica. Lejos de

141
q = 0 tal distincioón no existe. Por ejemplo en la zon de frontera, q = π/2a, las dos
lociones son:
 
2K 1/2
ω(q) = A0 = 0 (10.76)
M
 
2K 1/2
ω(q) = A=0 (10.77)
M0
con un tipo de átomo estacionario en cada caso.
Si M → M 0 este espectro es semejante al de una cadena monoatómica. El gap en
q = π/2a desaparece y la rama superior debe ser doblado la cubrir la región π/2a a
π/a.

10.6. Modos normales de una red de Bravais tridi-

mensional
Para una red de Bravais tridimensional, la energı́a potencial armónica viene dado
por(). Adoptamos la notación matricial
X
~ µν (R
uµ (R)D ~ −R
~ 0 )uν (R
~ 0 ) ≡ ~u(R)D(
~ ~ −R
R ~ 0 )~u(R
~ 0) (10.78)
µ,ν

Entonces
1 X ~ ~ −R
~ 0 )~u(R
~ 0)
U armo = ~u(R)D(R (10.79)
2 ~ ~0
R,R

donde
!
~ −R
~ 0) = ∂ 2 U arm
Dµν (R (10.80)
~ 0 )∂uν (R)
∂uµ (R ~
0

son las constantes de fuerza y en equilibrio

!
∂U arm
=0 (10.81)
~
∂uν (R) 0

Las formas particulares de las las fuerzas atómicas son llevadas en cuenta en las
~ −R
simetrı́as de las matrices matrices D(R ~ 0 ).

142
a) Simetrı́a 1.

~ −R
Dµν (R ~ 0 ) = Dνµ (R
~ 0 − R)
~ (10.82)

Sabemos que

~ −R
~ 0) = ∂ 2 U arm ∂ 2 U arm
Dµν (R =
~ 0 )∂uν (R)
∂uµ (R ~ ~
∂uν (R)∂u ~0
µ (R )
~ 0 − R)
= Dµν (R ~ (10.83)

independiente del orden de la diferenciación.

b) Simetrı́a 2.
Simetrı́a de invesión.
~ −R
Dµν (R ~ 0 ) = Dµν (R
~ 0 − R)
~ (10.84)

o
~ = D(−R)
D(R) ~ (10.85)

o de la simetrı́a 1 se ve
~ −R
Dµν (R ~ 0 ) = Dνµ (R
~ −R
~ 0) (10.86)

~ con desplazamiento uν (R)

es decir un ión en R ~ y con desplazamiento −uν (R)
~
no cambia.

c) Simetrı́a 3.

X
~ =0
Dµν (R) (10.87)
~
R

o
X
~ =0
D(R) (10.88)
~
R

Si cada ión es desplazado una distancia d~ (constante) tal que ~u(R) ~ esto
~ = d,
representa como si todo el cristal se desplaza como un todo sin cambio en la

143
 energı́a cinética (de los iones que componen el sistema). Si U arm tiene el mismo
vector de desplazamiento entonces
X
0 = ~ −R
dµ Dµν (R ~ 0 )dν
~ R
R, ~ 0 ,µ,ν
X X
= dµ dν ~ −R
Dµν (R ~ 0)
µ,ν ~ R
R, ~0
 
X X
= N dµ dν  ~ 
Dµν (R) (10.89)
µ,ν ~
R

donde se hizo un cambio de variable R ~ −R

~ 00 = R ~0

Tenemos 3N ecuaciones de movimiento (uno para cada componente de desplazamiento

de los N iones).
armo
~ = − ∂U
M üµ (R) (10.90)
~
∂üµ (R)
Entonces
X
~ = −
M üµ (R) ~ −R
Dµν (R ~ 0 )uν (R
~ 0) (10.91)
~ 0 ,ν
R

o
~
d2~u(R) X
M = − ~ −R
D(R ~ 0 )~u(R
~ 0) (10.92)
dt2 ~0 R

La solución armónica de la ecuación de movimiento es

~ ~
~ t) = ~ei(k·R−ωt)
~u(R, (10.93)

donde ~ es un vector que indica la dirección en que los iones se mueven, es denominado
vector de polarización de los modos normales.
De la condición de frontera periódica de Born-von Karman

~u(R ~
~ + Ni~ai ) = ~u(R) (10.94)

donde ~ai es el vector de la celda primitiva y Ni es el número de iones en la dirección

i (i = 1, 2, 3) y N = N1 N2 N3 es el número de iones en la red. Entonces, los valores
permitidos de ~k vienen dados por
~k = m1 ~b1 + m2 ~b2 + m3 ~b3 (10.95)
N1 N2 N3

144
donde ni es un entero y ~bi vectores primitivos de la red reciı́proca. Además, ~k y ~k + K
~
(dos vectores permitidos que se diferencian por un vector de la red recı́proca) dan la
~ t) = ~u~ ~ (R,
misma solución (~u~k (R, ~ t)).
k+K
Remplazando en la ecuación de movimiento

d2 i(~k·R−ωt)
~ X
~ 0 )~ei(~k·R~ 0 −ωt)
~ −R
M 2 ~e = − D(R (10.96)
dt ~0 R

Haciendo un cambio de variable R ~ −R

~ 00 = R ~0
~ ~ X ~ ~ ~ 00 )−ωt]
−Mω 2~ei(k·R−ωt) = −~ ~ 00 )ei[k·(R−R
D(R (10.97)
~ 00
R
 
X
~ 00 )e i~k·R
~ 00
 i[~k·R−ωt]
~
= − ~ D(R e (10.98)
~ 00
R

Entonces

Mω 2~ = D(~k)~ (10.99)

donde
X ~ ~ 00
D(~k) = ~ 00 )eik·R
D(R (10.100)
~ 00
R

~ es simétrica entonces D(~k) también es simétrica, entonces

Como D(R)
X ~ ~
D(~k) = ~ ik·R
D(R)e
~
R
1X h i
= ~ e−i~k·R~ + ei~k·R~ − 2
D(R)
2 ~
R
X h i
= ~ cos ~k · R
D(R) ~ −1
~
R
X  
= −2 ~ sen2 1 ~k · R
D(R) ~ (10.101)
~
2
R

Vemos que D(~k) es una función par y real, es decir, una función simétrica. Teorema:
cualquier matriz simetrı́ca real 3D tiene tres eigenvectores reales ~i (i = 1, 2, 3).

D(~k)~s = λs (~k)~s (10.102)

145
donde los eigenvectores pueden ser normalizados

~s · ~s0 = δss0 (10.103)

Remplazando

Mωs2~s = λs (~k)~s (10.104)

Entonces
v
u
u λ (~
t s k)
ωs = (10.105)
M

Los tres modos normales con vector de onda ~k deben tener ~s vectores de polarización
con frecuencia ωs .
En 3D, análogo para las tres ramas

10.7. Modos normales de la red del diamante

1 X
U armo = ~ + ~u − R
φµ,ν (R ~ 00 ) (10.106)
2 ~ R~ 0 ,µ,ν
R,

donde
X
~ −R
Dµ,ν (R ~ 0) = δ ~ ~ 0 ~ −R
φµ,ν (R ~ 00 ) − φµ,ν (R
~ −R
~ 0) (10.107)
R,R
~ 00
R

146