Estructura de la Materia 2

Clase 4 - Teoría

Docentes
Gustavo Grinblat, Andrea Barral, Tomás Bortolin, Agustina Casafuz

Departamento de Física, FCEN, UBA – 2do Cuatrimestre, 2020

Web: http://materias.df.uba.ar/edlm2a2020c2
 Programa de la materia

• Red cristalina, red recíproca y difracción de rayos X

• Clasificación de los sólidos y energía de cohesión

• Vibraciones, fonones y propiedades térmicas

• Electrones en sólidos (potencial periódico)

• Semiconductores y juntura semiconductora

• Magnetismo en sólidos

• Introducción a los aisladores topológicos

 Modos vibracionales
Posición de los núcleos en un cristal real
Los núcleos en un sólido no se encuentran fijos en el espacio.
ത 𝑡 = 𝑅ത + 𝑢ത 𝑅,
𝑟ҧ 𝑅, ത 𝑡
• Asumimos que los núcleos oscilan alrededor de posiciones
de equilibrio, las cuales determinan una RB.

O O La configuración cristalina varía instante a instante, y

la estructura observada corresponde a un promedio.
Potencial de int. entre núcleos

• Asumimos que el desplazamiento de los núcleos respecto a

las posiciones de equilibrio es mucho menor a la distancia
interatómica.
1 ത
𝑢ത (𝑅)
𝑈 = ෍ 𝜙 𝑟ҧ 𝑅ത − 𝑟ҧ 𝑅′ ത
2 ′
𝑅ത 𝑅ത
~𝑘𝐵 𝑇 (~25 𝑚𝑒𝑉 a 𝑇𝑎 ) 1 O
= ഥ + 𝑢 𝑅ത − 𝑢(𝑅′)
෍ 𝜙 𝑅ത − 𝑅′ ത
2 ′
Distancia entre núcleos 𝑅ത 𝑅ത
 Modos vibracionales: Aproximación armónica
Aproximación armónica en 1D (desarrollo de Taylor a segundo orden en torno al equilibrio)
1 2
𝜙 𝑅 − 𝑅′ + 𝑢 𝑅 − 𝑢(𝑅′) = 𝜙 𝑅 − 𝑅′ + 𝜙′ 𝑅 − 𝑅′ 𝑢 𝑅 − 𝑢(𝑅′) + 𝜙′′ 𝑅 − 𝑅′ 𝑢 𝑅 − 𝑢(𝑅′)
2
𝑓 𝑥 + ℎ ; 𝑥 = 𝑅 − 𝑅′ ; ℎ = 𝑢 𝑅 − 𝑢(𝑅′)
𝑈𝑒𝑞

1 1 1 1
𝑈 = ෍ 𝜙 𝑅 − 𝑅′ + 𝑢 𝑅 − 𝑢(𝑅′) = ෍ 𝜙 𝑅 − 𝑅′ + ෍ 𝜙′ 𝑅 − 𝑅′ 𝑢 𝑅 − 𝑢(𝑅′) + ෍[ ]
2 ′ 2 ′ 2 ′ 2 ′
𝑅𝑅 𝑅𝑅 𝑅𝑅 𝑅𝑅

=0
1 2 ෍ ෍ 𝜙′ 𝑅 − 𝑅′ 𝑢 𝑅 − ෍ ෍ 𝜙′ 𝑅 − 𝑅′ 𝑢 𝑅′
𝑈 = 𝑈𝑒𝑞 + ෍ 𝜙′′ 𝑅 − 𝑅′ 𝑢 𝑅 − 𝑢(𝑅′)
4 ′ 𝑅 𝑅′ 𝑅′ 𝑅
𝑅𝑅
\= 0
= 0 (Fuerza neta sobre el núcleo 𝑅 en el equilibrio)

𝜙′′ 𝑅 − 𝑅′ Constantes de fuerza determinadas por la curvatura del potencial de interacción.

 Modos vibracionales: Cadena unidimensional
Modos normales de una RB unidimensional

(𝑛 − 3)𝑎 (𝑛 − 1)𝑎 (𝑛 + 1)𝑎 (𝑛 + 3)𝑎

(𝑛 − 4)𝑎 (𝑛 − 2)𝑎 𝑅 = 𝑛𝑎 (𝑛 + 2)𝑎 (𝑛 + 4)𝑎 Resortes perfectos de masa nula y constante C

𝑎
𝑢(𝑛𝑎) Masa M
𝜙′′(𝑎)
1 2
1 2
Interacción a primeros vecinos: 𝑈𝑎 = ෍ 𝜙′′ 𝑅 − 𝑅′ 𝑢 𝑅 − 𝑢(𝑅′) = 𝐶 ෍ 𝑢 𝑖𝑎 − 𝑢 𝑖 + 1 𝑎
4 ′ 2
𝑅𝑅 𝑖

Ecuaciones de movimiento
𝑖 =n−1 𝑖=n

𝜕𝑈𝑎 1 𝜕 … + 𝑢 𝑛 − 1 𝑎 − 𝑢 𝑛𝑎 2 + 𝑢 𝑛𝑎 − 𝑢 𝑛 + 1 𝑎 2…
𝑀𝑢ሷ 𝑛𝑎 = − =− 𝐶
𝜕𝑢 𝑛𝑎 2 𝜕𝑢 𝑛𝑎

= −𝐶 2𝑢 𝑛𝑎 − 𝑢 𝑛 − 1 𝑎 − 𝑢 𝑛 + 1 𝑎
 Modos vibracionales: Cadena unidimensional
Modos normales de una RB unidimensional

(𝑛 − 3)𝑎 (𝑛 − 1)𝑎 (𝑛 + 1)𝑎 (𝑛 + 3)𝑎

(𝑛 − 4)𝑎 (𝑛 − 2)𝑎 𝑅 = 𝑛𝑎 (𝑛 + 2)𝑎 (𝑛 + 4)𝑎 Resortes perfectos de masa nula y constante C

𝑎
𝑢(𝑛𝑎) Masa M
Cadena finita
Si la red tiene un número finito N de núcleos con N muy grande, y no nos interesan efectos de borde, entonces
resulta irrelevante cómo tratemos a los extremos, y podemos elegir la condición que nos resulte conveniente.

Condiciones de contorno periódicas de Born-von Karman: 𝑢 𝑁 + 1 𝑎 = 𝑢 𝑎 ; 𝑢 0 = 𝑢 𝑁𝑎

2𝜋 𝑚
Buscamos soluciones de la forma: 𝑢 𝑛𝑎, 𝑡 = 𝜖𝑒 𝑖(𝑘𝑛𝑎−𝜔𝑡) 𝑒 𝑖𝑘𝑁𝑎 =1 𝑘= , −𝑚 entero
𝑎 𝑁
Como desplazar a 𝑘 en 2𝜋/𝑎 no altera el valor de 𝑢 𝑛𝑎 Existen exactamente 𝑁 soluciones
(modos normales) diferentes.
Elegimos tomar 𝑘 entre −𝜋/𝑎 y 𝜋/𝑎.
 Modos vibracionales: Cadena unidimensional
Modos normales de una RB unidimensional: Cadena finita
(𝑛 − 3)𝑎 (𝑛 − 1)𝑎 (𝑛 + 1)𝑎 (𝑛 + 3)𝑎
𝑢 Ambas ondas toman el mismo valor sobre los
distintos puntos de la red, y difieren solo entre ellos.
𝑘
Como la onda física está definida solo sobre
𝑘 + 2𝜋/𝑎
los puntos de la red, entonces ambas ondas
(𝑛 − 2)𝑎 𝑛𝑎 (𝑛 + 2)𝑎
son completamente equivalentes.

Frecuencias de modos normales

2𝐶 1 − cos 𝑘𝑎
𝑀𝑢ሷ 𝑛𝑎 = −𝐶 2𝑢 𝑛𝑎 − 𝑢 𝑛 − 1 𝑎 − 𝑢 𝑛 + 1 𝑎 𝜔 𝑘 =
𝑀
𝑢 𝑛𝑎, 𝑡 = 𝜖𝑒 𝑖(𝑘𝑛𝑎−𝜔𝑡)
𝐶 1
=2 sen 𝑘𝑎
−𝑀𝜔2 𝑒 𝑖 𝑘𝑛𝑎−𝜔𝑡
= −𝐶 2 − 𝑒 −𝑖𝑘𝑎 − 𝑒 𝑖𝑘𝑎 𝑒 𝑖(𝑘𝑛𝑎−𝜔𝑡) 𝑀 2

= −2𝐶 1 − cos(𝑘𝑎) 𝑒 𝑖(𝑘𝑛𝑎−𝜔𝑡) 1 − cos 𝑘𝑎 1

= sen2 𝑘𝑎
2 2
 Modos vibracionales: Cadena unidimensional
Modos normales de una RB unidimensional: Cadena finita
𝜔
4𝐶/𝑀 2da ZB 𝐶 1
𝜔 𝑘 =2 sen 𝑘𝑎
𝑀 2
𝑎
RD
2𝜋/𝑁𝑎
2𝜋/𝑎
RR
𝐾 ∈ 𝑅𝑅

Celda de WZ (1ZB)
𝑘 = −𝜋/𝑎 0 𝑘 = 𝜋/𝑎 2𝜋/𝑎 Velocidad de fase
Casos límites
= 𝑐 = 𝑣𝑔 = 𝜕𝜔/𝜕𝑘
¡Primera ZB!
Velocidad
Cualquier 𝑘 fuera de la 1ZB puede trasladarse a • 𝑘 ≪ 𝜋/𝑎 𝜔 𝑘 = 𝑎 𝐶/𝑀 𝑘 de grupo
un 𝑘’ equivalente dentro de la 1ZB a través de un Relación de dispersión de tipo sonido/luz (ω = 𝑐𝑘).
vector de la RR. • 𝑘 = ±𝜋/𝑎 𝑣𝑔 = 0
 Modos vibracionales: Cadena unidimensional con base
Modos normales de una red unidimensional con una base (átomos iguales)
(𝑛 − 3)𝑎 (𝑛 − 1)𝑎 (𝑛 + 1)𝑎 (𝑛 + 3)𝑎
(𝑛 − 4)𝑎 (𝑛 − 2)𝑎 𝑛𝑎 (𝑛 + 2)𝑎 (𝑛 + 4)𝑎 C >G Masa M

CP 1 2 𝑑 𝑎−𝑑 𝑎

𝐶 2
𝐺 2
Interacción a primeros vecinos: 𝑈𝑎 = ෍ 𝑢1 𝑖𝑎 − 𝑢2 𝑖𝑎 + ෍ 𝑢2 𝑗𝑎 − 𝑢1 𝑗 + 1 𝑎
2 2
𝑖 𝑗

Ecuaciones de movimiento
𝑖=n 𝑗 =n−1

𝜕𝑈𝑎 1 𝜕 … + 𝐶 𝑢1 𝑛𝑎 − 𝑢2 𝑛𝑎 2 + 𝐺 𝑢2 𝑛 − 1 𝑎 − 𝑢1 𝑛𝑎 2
…
𝑀𝑢ሷ 1 𝑛𝑎 = − =−
𝜕𝑢1 𝑛𝑎 2 𝜕𝑢1 𝑛𝑎
= −𝐶 𝑢1 𝑛𝑎 − 𝑢2 𝑛𝑎 − 𝐺 𝑢1 𝑛𝑎 − 𝑢2 𝑛 − 1 𝑎

𝜕𝑈𝑎
𝑀𝑢ሷ 2 𝑛𝑎 = − = −𝐶 𝑢2 𝑛𝑎 − 𝑢1 𝑛𝑎 − 𝐺 𝑢2 𝑛𝑎 − 𝑢1 𝑛 + 1 𝑎
𝜕𝑢2 𝑛𝑎
 Modos vibracionales: Cadena unidimensional con base
Modos normales de una red unidimensional con una base (átomos iguales)
(𝑛 − 3)𝑎 (𝑛 − 1)𝑎 (𝑛 + 1)𝑎 (𝑛 + 3)𝑎
(𝑛 − 4)𝑎 (𝑛 − 2)𝑎 𝑛𝑎 (𝑛 + 2)𝑎 (𝑛 + 4)𝑎 C >G Masa M

CP 1 2
Frecuencias y amplitudes de modos normales
𝑀𝑢ሷ 1 𝑛𝑎 = −𝐶 𝑢1 𝑛𝑎 − 𝑢2 𝑛𝑎 − 𝐺 𝑢1 𝑛𝑎 − 𝑢2 𝑛 − 1 𝑎 𝑢1 𝑛𝑎, 𝑡 = 𝜖1 𝑒 𝑖(𝑘𝑛𝑎−𝜔𝑡)
𝑀𝑢ሷ 2 𝑛𝑎 = −𝐶 𝑢2 𝑛𝑎 − 𝑢1 𝑛𝑎 − 𝐺 𝑢2 𝑛𝑎 − 𝑢1 𝑛 + 1 𝑎 𝑢2 𝑛𝑎, 𝑡 = 𝜖2 𝑒 𝑖(𝑘𝑛𝑎−𝜔𝑡)

−𝑀𝜔2 𝜖1 = −𝐶 𝜖1 − 𝜖2 − 𝐺 𝜖1 − 𝜖2 𝑒 −𝑖𝑘𝑎 𝑀𝜔2 − 𝐶 + 𝐺 𝜖1 + 𝐶 + 𝐺𝑒 −𝑖𝑘𝑎 𝜖2 = 0

−𝑀𝜔2 𝜖2 = −𝐶 𝜖2 − 𝜖1 − 𝐺 𝜖2 − 𝜖1 𝑒 𝑖𝑘𝑎 𝑀𝜔2 − 𝐶 + 𝐺 𝜖2 + 𝐶 + 𝐺𝑒 𝑖𝑘𝑎 𝜖1 = 0

𝑀𝜔2 − 𝐶 + 𝐺 2
= 𝐶 2 + 𝐺 2 + 2𝐶𝐺𝑐𝑜𝑠(𝑘𝑎) Dividiendo entre sí

𝐶+𝐺 1 𝜖2 𝐶 + 𝐺𝑒 𝑖𝑘𝑎
𝜔2 = ± 𝐶 2 + 𝐺 2 + 2𝐶𝐺𝑐𝑜𝑠(𝑘𝑎) =∓
𝑀 𝑀 𝜖1 𝐶 + 𝐺𝑒 𝑖𝑘𝑎
 Modos vibracionales: Cadena unidimensional con base
Modos normales de una red unidimensional con una base (átomos iguales): Cadena finita

2
𝐶+𝐺 1 Casos límites
𝜔 = ± 𝐶 2 + 𝐺 2 + 2𝐶𝐺𝑐𝑜𝑠(𝑘𝑎) 𝜖2
𝑀 𝑀 • 𝑘 ≪ 𝜋/𝑎 = ∓1 (−: RO; +: RA)
𝑢1 = 𝜖1 𝑒 𝑖(𝑘𝑛𝑎−𝜔𝑡) 𝜖2 𝐶 + 𝐺𝑒 𝑖𝑘𝑎 𝜖1
=∓ RA: 𝜔 ⁓ kHz RA: Frecuencias asociadas al espectro audible.
𝑢 = 𝜖 𝑒 𝑖(𝑘𝑛𝑎−𝜔𝑡) 𝜖1 𝐶 + 𝐺𝑒 𝑖𝑘𝑎
2 2

𝜔 C G
2(𝐶 + 𝐺)/𝑀
RO: 𝜔 ⁓ THz (En un sólido iónico podría acoplarse a una OE)
+ − + − + − + − + − + − + − +
2𝐶/𝑀 RO: Frecuencias asociadas a luz infrarroja.
𝜖2
2𝐺/𝑀 • 𝑘 = 𝜋/𝑎 = ∓1 (−: RO; +: RA)
𝜖1
N: número de
RA:
2𝜋/𝑁𝑎
celdas primitivas.
Como hay dos
ramas, hay 2𝑁 modos C G
2N modos normales RO:
normales.
𝑘 = −𝜋/𝑎 0 𝑘 = 𝜋/𝑎
 Modos vibracionales: Planteo general
Vibraciones en RB + Base (3D)
Punto de la RB
ҧ = 𝑅ത𝑚𝛼 + 𝑢ത 𝑚𝛼 𝑖, 𝑗, 𝑙 = Coordenadas cartesianas (𝑥, 𝑦, 𝑧)
𝑟𝑚𝛼 𝑅ത𝑚𝛼 = 𝑅ത𝑚 + 𝑑ҧ𝛼
Índices 𝑚, 𝑛, 𝑝, 𝑞 = Puntos de la RB (1, 2, …, N)
Elemento de la base 𝛼, 𝛽, 𝛾, 𝜎 = Elementos de la base (1, 2, …, P)
𝑖 : Proyección en la dirección 𝑖Ƹ del vector 𝑟ҧ
𝑟𝑚𝛼
𝑎ത1 = 𝑎𝑥ො 𝑎 𝑚𝛼
𝑎ത2 = 𝑎𝑦ො
𝑑1ҧ = 0ത 𝑥ො 1
𝑑ҧ2 = (𝑎/2 )(𝑥ො + 𝑦)
ො RB 𝑈 = ෍ 𝜙 𝑟ҧ 𝑅ത − 𝑟ҧ 𝑅′
ത
2 ′
O ത 𝑅ത
𝑅,
𝑚 = 1, 𝛼 = 2 𝜙: Potencial de interacción entre pares de núcleos

1
RB + Base 𝑈 = ෍ ෍ 𝜙 𝑟𝑝𝛾
ҧ − 𝑟𝑞𝜎
ҧ
2
𝑝,𝑞 𝛾,𝜎
 Modos vibracionales: Planteo general
Aproximación armónica en 3D 𝑖, 𝑗, 𝑙 = Coordenadas cartesianas (𝑥, 𝑦, 𝑧)
𝑚, 𝑛, 𝑝, 𝑞 = Puntos de la RB (1, 2, …, N)
𝑈 𝑟11
ҧ , … , 𝑟1𝑃
ҧ , … , 𝑟𝑁1
ҧ , … , 𝑟𝑁𝑃
ҧ = 𝑈 𝑅ത11 + 𝑢ത 11 , … , 𝑅ത𝑁𝑃 + 𝑢ത 𝑁𝑃 = 𝛼, 𝛽, 𝛾, 𝜎 = Elementos de la base (1, 2, …, P)

𝜕𝑈 𝑅ത11 , … , 𝑅ത𝑁𝑃 𝑖 1 𝜕 2 𝑈 𝑅ത11 , … , 𝑅ത𝑁𝑃 𝑖 𝑗

= 𝑈 𝑅ത11 , … , 𝑅ത𝑁𝑃 + ෍ 𝑖
𝑢𝑛𝛽 + ෍෍෍ 𝑗
𝑢 𝑚𝛼 𝑢 𝑛𝛽
𝜕𝑟𝑛𝛽 2 𝑖
𝜕𝑟𝑚𝛼 𝑟
𝑛,𝛽,𝑖 𝑚,𝑛 𝛼,𝛽 𝑖,𝑗 𝑛𝛽

Taylor a 𝑈𝑒𝑞
= 0 (Componente 𝑖 de la fuerza neta sobre el núcleo 𝑛𝛽 en el equilibrio)
2do orden
𝑗 𝑗
1 𝜕𝑈 1 𝜕𝜙 𝑟𝑝𝛾
ҧ − 𝑟𝑞𝜎
ҧ (𝑟𝑝𝛾 − 𝑟𝑞𝜎 )
𝑈 = ෍ ෍ 𝜙 𝑟𝑝𝛾
ҧ − 𝑟𝑞𝜎
ҧ 𝑗
= ෍෍ 𝑗
𝜕 𝑗
2 𝜕𝑟𝑛𝛽 2 𝜕𝑟 𝜕𝑟
𝑝,𝑞 𝛾,𝜎 𝑝,𝑞 𝛾,𝜎 𝑛𝛽

𝛿𝑝𝑛 𝛿𝛾𝛽 − 𝛿𝑞𝑛 𝛿𝜎𝛽

𝑖 − 𝑟𝑖 )
𝜕2𝑈 1 𝜕𝜙 𝑟𝑝𝛾
ҧ − 𝑟𝑞𝜎 ҧ (𝑟𝑝𝛾 𝑞𝜎
= ෍෍ (𝛿 𝛿
𝑝𝑛 𝛾𝛽 − 𝛿 𝛿
𝑞𝑛 𝜎𝛽 )𝜕
𝑖
𝜕𝑟𝑚𝛼 𝜕𝑟𝑛𝛽
𝑗 2 𝜕𝑟 𝑖 𝜕𝑟 𝑗 𝑖
𝜕𝑟𝑚𝛼
𝑝,𝑞 𝛾,𝜎

𝛿𝑝𝑚 𝛿𝛾𝛼 − 𝛿𝑞𝑚 𝛿𝜎𝛼

 Modos vibracionales: Planteo general
Aproximación armónica en 3D 𝑖, 𝑗, 𝑙 = Coordenadas cartesianas (𝑥, 𝑦, 𝑧)
𝑚, 𝑛, 𝑝, 𝑞 = Puntos de la RB (1, 2, …, N)
1 𝜕 𝑈 𝑅ത11 , … , 𝑅ത𝑁𝑃
2
𝑖 𝑢𝑗 = 𝑈 + 𝑈 𝛼, 𝛽, 𝛾, 𝜎 = Elementos de la base (1, 2, …, P)
𝑈 = 𝑈𝑒𝑞 + ෍ ෍ ෍ 𝑗
𝑢𝑚𝛼 𝑛𝛽 𝑒𝑞 𝑎
2 𝑖
𝜕𝑟𝑚𝛼 𝑟𝑛𝛽
𝑚,𝑛 𝛼,𝛽 𝑖,𝑗

𝜕2𝑈 1 𝜕 2 𝜙 𝑟𝑝𝛾
ҧ − 𝑟𝑞𝜎 ҧ
𝑗
= ෍෍ 𝛿𝑝𝑛 𝛿𝛾𝛽 − 𝛿𝑞𝑛 𝛿𝜎𝛽 𝛿𝑝𝑚 𝛿𝛾𝛼 − 𝛿𝑞𝑚 𝛿𝜎𝛼 =
𝑖
𝜕𝑟𝑚𝛼 𝜕𝑟𝑛𝛽 2 𝜕𝑟 𝑖 𝜕𝑟 𝑗
𝑝,𝑞 𝛾,𝜎
𝛿𝑝𝑛 𝛿𝑝𝑚 𝛿𝛾𝛽 𝛿𝛾𝛼 − 𝛿𝑞𝑛 𝛿𝜎𝛽 𝛿𝑝𝑚 𝛿𝛾𝛼 − 𝛿𝑝𝑛 𝛿𝛾𝛽 𝛿𝑞𝑚 𝛿𝜎𝛼 + 𝛿𝑞𝑛 𝛿𝜎𝛽 𝛿𝑞𝑚 𝛿𝜎𝛼

1 𝜕 2 𝜙 𝑟𝑚𝛼
ҧ − 𝑟𝑞𝜎 ҧ 1 𝜕 2 𝜙 𝑟𝑚𝛼
ҧ − 𝑟ҧ𝑛𝛽 1 𝜕 2 𝜙 𝑟ҧ𝑛𝛽 − 𝑟𝑚𝛼
ҧ 1 𝜕 2 𝜙 𝑟𝑝𝛾
ҧ − 𝑟𝑚𝛼 ҧ
= ෍ 𝛿𝑚𝑛 𝛿𝛼𝛽 − − + ෍ 𝛿𝑚𝑛 𝛿𝛼𝛽
2 𝜕𝑟 𝑖 𝜕𝑟 𝑗 2 𝜕𝑟 𝑖 𝜕𝑟 𝑗 2 𝜕𝑟 𝑖 𝜕𝑟 𝑗 2 𝜕𝑟 𝑖 𝜕𝑟 𝑗
𝑞,𝜎 𝑝𝛾

𝜕 2 𝜙 𝑟𝑚𝛼
ҧ − 𝑟𝑞𝜎 ҧ 𝜕 2 𝜙 𝑟𝑚𝛼
ҧ − 𝑟ҧ𝑛𝛽 1 𝑖𝑗 𝑖𝑗 𝑖 𝑗
=෍ 𝛿𝑚𝑛 𝛿𝛼𝛽 − 𝑈𝑎 = ෍ ෍ ෍ ෍ 𝜙𝑚𝑞 𝛿𝑚𝑛 𝛿𝛼𝛽 − 𝜙𝑚𝑛 𝑢𝑚𝛼 𝑢𝑛𝛽
𝜕𝑟 𝑖 𝜕𝑟 𝑗 𝜕𝑟 𝑖 𝜕𝑟 𝑗 2 𝛼𝜎 𝛼𝛽
𝑞,𝜎 𝑚,𝑛 𝛼,𝛽 𝑖,𝑗 𝑞,𝜎
𝑖𝑗
Evaluado en el equilibrio lo denotamos como: 𝜙𝑚𝑞
𝛼𝜎
 Modos vibracionales: Planteo general
Ecuaciones de movimiento y modos normales 𝑖, 𝑗, 𝑙 = Coordenadas cartesianas (𝑥, 𝑦, 𝑧)
𝑚, 𝑛, 𝑝, 𝑞 = Puntos de la RB (1, 2, …, N)
𝑙
𝜕𝑈𝑎
𝑀𝛾 𝑢ሷ 𝑝𝛾 =− 𝑙 = 𝛼, 𝛽, 𝛾, 𝜎 = Elementos de la base (1, 2, …, P)
𝜕𝑢𝑝𝛾
1 𝑖𝑗 𝑖𝑗 𝑖 𝑢𝑗
𝜕 σ𝑚,𝑛 σ𝛼,𝛽 σ𝑖,𝑗 σ𝑞,𝜎 𝜙𝑚𝑞 𝛿𝑚𝑛 𝛿𝛼𝛽 − 𝜙𝑚𝑛 𝑢𝑚𝛼
2 𝛼𝜎 𝛼𝛽
𝑛𝛽
𝑖𝑙 𝑖𝑙 𝑖
=− 𝑙 = − ෍ ෍ 𝜙𝑚𝑞 𝛿𝑚𝑝 𝛿𝛼𝛾 − 𝜙𝑚𝑝 𝑢𝑚𝛼
𝜕𝑢𝑝𝛾 𝛼𝜎 𝛼𝛾
𝑚,𝛼,𝑖 𝑞,𝜎

Sistema de 3×N×P ecuaciones acopladas (número igual al de grados de libertad del sistema)

Buscamos soluciones de la forma: 𝑢ത 𝑝𝛾 = 𝜖𝛾ҧ 𝑒 𝑖(𝑘ത 𝑅ത𝑝−𝜔𝑡)

തത ത
𝑖𝑙
−𝑀𝛾 𝜔2 𝜖𝛾𝑙 𝑒 𝑖(𝑘𝑅𝑝−𝜔𝑡) = − ෍ ෍ 𝜙𝑚𝑞 𝑖𝑙
𝛿𝑚𝑝 𝛿𝛼𝛾 − 𝜙𝑚𝑝 𝜖𝛼𝑖 𝑒 𝑖(𝑘𝑅ത𝑚 −𝜔𝑡)
𝛼𝜎 𝛼𝛾
𝑚,𝛼,𝑖 𝑞,𝜎

ത ത ത𝑖𝑙 𝑖𝑙 ത 𝑖
𝑀𝛾 𝜔2 𝜖𝛾𝑙 = ෍ ෍ 𝑒 𝑖𝑘(𝑅𝑚 −𝑅𝑝) ෍ 𝜙𝑚𝑞 𝛿𝑚𝑝 𝛿𝛼𝛾 − 𝜙𝑚𝑝 𝜖𝛼𝑖 = ෍ 𝔇𝑖𝑙
𝛼𝛾 (𝑘) 𝜖𝛼
𝛼𝜎 𝛼𝛾
𝛼,𝑖 𝑚 𝑞,𝜎 𝛼,𝑖
 Modos vibracionales: Planteo general
Ecuaciones de movimiento y modos normales 𝑖, 𝑗, 𝑙 = Coordenadas cartesianas (𝑥, 𝑦, 𝑧)
𝑚, 𝑛, 𝑝, 𝑞 = Puntos de la RB (1, 2, …, N)
𝑀𝛾 𝜔2 𝜖𝛾𝑙 = ෍ 𝔇𝑖𝑙 ത 𝑖
𝛼𝛾 (𝑘) 𝜖𝛼 𝛼, 𝛽, 𝛾, 𝜎 = Elementos de la base (1, 2, …, P)
𝛼,𝑖
𝑖𝑙 ത
𝑖𝑙 ത (𝑅ത𝑚 −𝑅ത𝑝 ) 𝑖𝑙
𝐷𝛼𝛾 𝑘 : Elemento genérico de
𝔇𝑖𝑙 ത
𝛼𝛾 𝑘 = ෍ ෍ 𝜙𝑚𝑞 𝛿𝑚𝑝 𝛿𝛼𝛾 − 𝑒 𝑖𝑘
𝜙𝑚𝑝 ന
𝛼𝜎 𝛼𝛾
la matriz dinámica 𝐷
𝑚 𝑞,𝜎
𝜖𝛼ǁ 𝑖 𝔇𝑖𝑙 ത
𝛼𝛾 (𝑘)
Definiendo: 𝜖𝛾ǁ 𝑙 = 𝑀𝛾 𝜖𝛾𝑙 𝑀𝛾 𝜔2 𝜖𝛾ǁ 𝑙 = ෍ 𝔇𝑖𝑙 ത
𝛼𝛾 (𝑘)
෍ − 𝜔2 𝛿𝑖𝑙 𝛿𝛼𝛾 𝜖𝛼ǁ 𝑖 = 0
𝑀𝛼 𝛼,𝑖
𝑀𝛼 𝑀𝛾
𝛼,𝑖

ന − 𝜔2 𝕀 𝜖ǁ = 0
𝐷 Ecuación de autovalores y autovectores

Producto de las proyecciones de los versores

𝑖𝑗 𝜕 2 𝜙 𝑟𝑚𝛼
ҧ − 𝑟𝑞𝜎 ҧ 𝜕 2 𝜙 𝑟ҧ 𝜕 𝑟𝑗 cartesianos de las coordenadas involucradas
𝜙𝑚𝑞 = อ = อ = 𝑖 𝜙′ 𝑟ҧ อ
𝛼𝜎 𝜕𝑟 𝑖 𝜕𝑟 𝑗 𝜕𝑟 𝑖 𝜕𝑟 𝑗 𝜕𝑟 𝑟ҧ sobre la dirección que conecta a los átomos.
𝐸𝑞 𝐸𝑞 𝐸𝑞
=0 =0
𝑖 𝑖 ) (𝑅 𝑗 𝑗
𝑟𝑖 𝑟𝑗 𝛿𝑖𝑗 𝑟𝑖 𝑟𝑗 (𝑅𝑚𝛼 − 𝑅𝑞𝜎 𝑚𝛼 − 𝑅𝑞𝜎 )
= 𝜙 ′′ 𝑟ҧ ′
+ 𝜙 𝑟ҧ − 𝜙′ 𝑟ҧ = 𝐶𝑚𝑞𝛼𝜎 2
𝑟ҧ 𝑟ҧ 𝑟ҧ 𝑟ҧ 3 𝐸𝑞 𝑅ത𝑚𝛼 − 𝑅ത𝑞𝜎
 Modos vibracionales: Planteo general
Condiciones de contorno periódicas de Born-von Karman
ത ത
𝑢 𝑅ത𝑚𝛼 + 𝑁𝑖 𝑎ത𝑖 = 𝑢 𝑅ത𝑚𝛼 , 𝑎ത𝑖 :VP de la RD; 𝑁1 𝑁2 𝑁3 = 𝑁 (𝑁° total de CP en el cristal) 𝑒 𝑖𝑘(𝑅ത𝑚+𝑁𝑖 𝑎ത𝑖) = 𝑒 𝑖𝑘𝑅ത𝑚
𝑚 𝑚 𝑚
ത
𝑒 𝑖𝑘𝑁𝑖 𝑎ത𝑖 =1 ത 𝑖 𝑎ത𝑖 = 2π𝑚, 𝑚 ∈ ℤ
𝑘𝑁 𝑘ത = 1 𝑏ത1 + 2 𝑏ത2 + 3 𝑏ത3 , 𝑚𝑖 ∈ ℤ, 𝑏ത𝑖 : VP de la RR.
𝑁1 𝑁2 𝑁3

ത ഥ 𝑅ത N valores no equivalentes de ത
𝑘.
ഥ
Desplazamientos de 𝑘 en 𝐾 ∈ RR no cambian la solución (𝑒 𝑖𝐾
= 1)
Elegimos tomarlos dentro de la 1ZB.
Ramas acústicas y ópticas
𝜔
Tenemos 3 × 𝑁 × 𝑃 modos normales y
ത
𝑁 valores no equivalentes de 𝑘.
3RO
3 × 𝑃 soluciones para cada valor de 𝑘ത
3 son RA 3RA
3 × 𝑃 ramas
3 × (𝑃 − 1) son RO

d dimensiones 𝑑 RA y 𝑑 × (𝑃 − 1) RO
𝑘
Cristal 3D con 2 elementos en la base
 Resumen

• Posición de núcleos en un cristal real

• Aproximación armónica en 1D

• Modos normales de una cadena lineal (RB)

• Modos normales de una cadena lineal (RB + base)

𝜔
RO
• Sólidos finitos y condiciones de contorno periódicas
RA

• Planteo general en 3 dimensiones

𝑘