Fonones

dinámica de redes
vibraciones cristalinas
Alvaro Morales
Instituto de Fisica
• Un solido se puede considerer como un
Sistema de electrons y iones acoplado
• El Hamiltoniano H (energia mecanica
total) para un solido.
H = He + Hn + He-n

He = energia cinetica energia cinetica de los Interaccion

de los electrons + la iones + la interaccion electron-ion
interaccion e-e n-n
 El Hamiltoniano
Para el solido es:

H = He + Hn + He-n
La aproximación mas importante es, la
aproximacion de un electron, donde se studia un
solo electron moviendose en el potencial efectivo
de todos los demas electrons y iones
 Approximacion #1: separar los electrons en:
Electrones de la coraza & electrons de Valencia
• Ejemplo:
Atomo de Si :

1s22s22p63s23p2
Electrones de la coraza = 1s22s22p6 (capas
cerradas)

Electrones de Valencia: 3s23p2

 El Hamiltoniano
para un crsital perfecto y periodico :
Ne electrones, Ni iones; Ne, Ni ~ 1023

H = He + Hi + He-i

He= ∑i(pi)2/(2mi) + (½)∑i∑i´[e2/|ri - ri´|]

0
(i  i´)

Hi = ∑j(Pj)2/(2Mj) + (½)∑j∑j´[ZjZj´ e2/|Rj - Rj´|] (j  j´)

0

He-i= - ∑i∑j[Zje2/|ri - Rj|]

 Approximacion # 2: aproximacion de Born-
Oppenheimer
• Permite separar el movimiento de electrons y
iones

(me/Mi) ~ 10-3 (<< 1)

 Los iones se mueven mas
lentamente que los electrones
• Frecuencias vibracionales tipicas: υi ~ 1013 s-1
 La escala de tiempo de los iones: ti ~ 10-13 s
• El movimiento electronico ocurre a energias:
Eg= hυe = ħω ~ 1 eV
 υe ~ 1015 s-1  te ~ 10-15 s
• Por lo tanto,
los Electrones responden al mov. ionico
~ Instantaneamente!

Para los electrons los iones estan quietos,

estacionarios
• En la aproximacion de Born-Oppenheimer, el
hamiltonianop de los iones es

• HI = Hi + He-i = j[(Pj)2/(2Mj)] +
(½)jj´[ZjZj´e2/|Rj-Rj´|] - ij[Zje2/|ri-Rj|]
• .
 Relación de dispersión (q)
Para ondas elásticas   vs q
Velocidad depende de la
vs  Y/ dirección, ondas trasversals
Y = Modulo elástico de Bulk viajan a velocidades menores
ρ = Densidad de masa

q→ 0 a longitudes de onda λ
grandes sistema NO
ω Continuo DISPERSIVO (no scattering)
q → ∞ a pequeños λ medio
dispersivo (aparece scattering )
discreto Cuando k aumenta la velocidad
decrece, si q continua
aumentando la dispersión se
q hace importante y el sistema es
altamente dispersivo
 Velocidad del sonido para sólidos típicos
Modulo Velocidad
Distancia Densidad elástico de Velocidad
del
Tipo de a primeros bulk de la onda
Sólido vecinos ρ calculada
sonido
estructura medido
(kg/m3) Y
(A°) (m/s)
(1010 N/m2) (m/s)

Sodio B.C.C 3.71 970 0.52 2320 2250

Cobre F.C.C 2.55 8966 13.4 3880 3830
Aluminio F.C.C 2.86 2700 7.35 5200 5110
Plomo F.C.C 3.49 11340 4.34 1960 1320
Silicio Diamante 2.35 2330 10.1 6600 9150
Germanio Diamante 2.44 5360 7.9 3830 5400
NaCl Rocksalt 2.82 2170 2.5 3400 4730
•VL los valores son comparables con las observaciones directas de la
velocidad del sonido. La velocidad del sonido es del orden de 5000 m/s en
sólidos metálicos, covalentes y ionicos.
Vibraciones Reticulares
Una onda vibracional en una red
cristalina es una secuencia repetitiva
de desplazamientos atómicos,
longitudinales, trasversales o una
combinación de ambos.
La ecuación de movimiento para
cualquier desplazamiento.
puede deducirse de las fuerzas
restauradoras sobre los átomos
producidas por los potenciales
moleculares .

[ F=-V(r) ]
 Vibraciones reticulares de una red 1D con átomos
idénticos
Los átomos interactúan vía un potencial V(r) que puede ser
desarrollado en un serie de Taylor.
 r  a   d 2V 
2
 dV 
V (r )  V (a)  (r  a)     2   ...........
 dr r  a 2  dr r a
0
V(R) Esta ecuación es similar a la energía potencial
de un resorte de constante :

 d 2V 
Repulsivo K   2 
 dr  r a
0 r0=4 Se puede relacionar K con el
módulo elástico de Young Y
Atractivo
min
r  a
R Fuerza  Y 
r a
Fuerza  K (r  a) Y  Ka
Cadena Lineal Monoatómica
La cadena consiste de un número muy grande de átomos idénticos con
idénticas masas.
Los átomos están separados por una distancia “a”.
Los átomos se mueven solo en la dirección paralela a la cadena 1D
Solo se considera interacción con los próximos vecinos.

a a a a a a

Un-2 Un-1 Un Un+1 Un+2

Un+1 > Un > Un-1
a a a a a a

Un-2 Un-1 Un Un+1 Un+2

a a a a a a

Un-2 Un-1 Un Un+1 Un+2

 Cadena Lineal Monoatómica

a a

Un-1 Un Un+1

K (u n1  u n ) Fuerza a la derecha K (u n  u n1 ) Fuerza a la izquierda

Fuerza total = Fuerza a la derecha – Fuerza a la izquierda

..
m u n  K (un1  2un  un1 )
 Cadena Lineal Monoatómica

Se considera una solución en ondas viajeras

un  A exp i  qxn0  t  xn0  na Posición de equilibrio

 i A exp i  qxn  t 

. dun
un   0

dt
d 2un
u n  2   i   2 A exp i  qxn0  t 
..
2

dt

xn  na  un
..
u n   2un Posición instantánea
 Ecuación de movimiento para el enesimo átomo
..
mu n  K (un1  2un  un1 )

2
 m A e

i qxn0 t 
 K Ae 
i qxn10 t 
 2A e

i qxn0 t 
 Ae

i qxn10 t 

qna q ( n  1) a qna q ( n  1) a

2
 m A e
i  qna t 
 K Ae  i  qna  ka t 
 2A e
i  qna t 
 Ae
i  qna  qa t 

2
 m A e
i  qna t 
 K Ae  i  qna t 
e
ika
 2A e
i  qna t 
 Ae
i  qna t 
e
 iqa

Cancelando términos comunes

m 2  K  eiqa  2  eiqa 

 Relación de dispersión (q)
m 2  K  eiqa  2  eiqa  eiqa  eiqa  2cos qa
m 2  K  2cos qa  2 
 2 K (1  cos qa)
1  cos x   2sin 2 
x

2 2  qa  2
m  4 K sin  
 2 
2 4K 2  qa 
  sin  
m  2 
4K  qa 
 sin   Valor máximo 1 cuando q=π/a
m  2 

4K
max   min  2a
m
 Análisis de la relación de dispersión (q)
ω  ωmax Sen(qa/2) 4K ω (q)
ωmax 
m ω  vs q
• Obsérvese la periodicidad:
2 ω (q)
ω(q)  ω(q  ) y ω(-q)  ω(q) 4K
a m
Adicionando un múltiplo de of 2π/a a q no
produce cambio en la frecuencia o en su
velocidad de grupo, así para puntos fuera
de la 1ZB no tienen significado físico
Las soluciones planteadas corresponden
a osciladores desacoplados cada uno
llamado un modo normal. Cada q define
un  oscilando independientemente uno
q
del otro. π π
Hay N modos de vibración en la primera  2π
a a
zona de Brillouin Na
L=Na ( Numero de átomos en la cadena)
1ZB
Análisis de la relación de dispersión (q)
Para longitudes de onda grande, es
decir, q próximo a cero, la relación de
Modelo
dispersión
continuo
4K
ω  Sen(qa/2)
m Modelo
atómico
 K 
Se reduce a : ω   a  q
 m 
 
ω  vs q vs 
K
a
m
Aproximación continua
Vs  Yc 
Vs   d V 
2
La medida de la velocidad del sonido nos
permite determinar las fuerzas
K   2 
 dr r a interatómicas y las constantes elásticas
 Análisis de la relación de dispersión (q)
4K  qa   Ka 2  qa 
 sin   vq   cos  
m  2  q m  2 
Para longitudes de onda grandes v f  v g  vs
Para valores de q en la frontera de zona de Brillouin:
 q=±π/a θ=90o tiene un significado
q=  vg  0 especial: Reflexiones de Bragg se
a obtienen a q= ±nπ/a

Las ondas elásticas en la frontera de zona de Brillouin

experimentan difracción de Bragg
 Modos de vibración en una cadena monoatómica
L
No todos los modos son posibles debido a
el átomo N debe ser el mismo que el (N+n)
cuando la cadena se cierra sobre si misma
Se requiere que haya un número entero de
longitudes de onda en el en el anillo

L  Na  p p entero
Así en la 1ZB en el rango 2π/a de q,
hay N valores permitidos de q .
Na 2 2
L  Na  p      Nk  p
p k a

L 1 dq
g(ω)  g(ω)  2 g(q)
π dω dω
dq
1
2L  qa 
g(ω)  Cos ( 2 ) 
π a max
 Cadena lineal 1D con dos tipos de átomos
 Dos tipos diferentes de átomos estan conectados por resortes
identicos de constante K

(n-2) (n-1) (n) (n+1) (n+2)

K K K K
M m M m M a)

b)

Un-2 Un-1 Un Un+1 Un+2

Este es el modelo más simple para un cristal iónico: la distancia entre

próximos vecinos es a/2 y la constante de red a.

Dos ecuaciones de movimiento: una para el
Relación de dispersión  ( q ) átomo de masa M y otra para el átomo de
masa m
 Cadena diatómica
M m M m M

Un-2 Un-1 Un Un+1 Un+2

Equación de movimiento para la masa M (n):

..
M u n  K (un1  un )  K (un  un1 )
M un  K (un1  2un  un1 )
..

Equación de movimiento para la masa m (n-1):

..
mu n  K (un  un1 )  K (un1  un2 )
-1

..
mu n  K (un  2un1  un2 )
-1
 Cadena diatómica

M m M m M

Un-2 Un-1 Un Un+1 Un+2

Solución para la masa M

un  A exp i  qxn0  t   xn0  na / 2

Para la masa m

un   A exp i  qxn0  t  

-1

α : número complejo que determina la amplitud relativa y la

fase de la onda de vibración.
u n   2 A exp i  qxn0  t  
..
 Cadena diatómica
Para el átomo M (n):
..
M u n  K (un1  2un  un1 )
 q n 1a  q n 1a
 qna
i

t   i

t   qna
i

t  i t  


 2 MAe  2 
 K   Ae  
 2 Ae  2    Ae  
2 2

 
 
 qna
i

t    qna
i

t  i qa
 qna
i

t 
 qna
i t   i qa 


 2 MAe  2 
 K   Ae  2  e 2  2 Ae  2    Ae  2  e 2 
 
 
Cancelando términos

 i
qa
i
qa

 M  K   e  2   e 
2 2 2 eix  eix  2cos x
 
 qa 
 2 M  2 K 1   cos 
 2 
 Cadena diatómica
Para el átomo m (n-1)
..
mu n1  K (un  2un1  un2 )
 q n1a  q n1a  q n  2 a
i

t   i qna t  i

t  i t  


 A 2 me  
 K  Ae  2   2 Ae  
 Ae  
2 2 2

 
 
 qna
i

t   i qa  i qna 
t 
 qna
i

t   i qa
 qna
i t   i 2 qa 


 2 mAe  2 
e 2
 K  Ae  2   2 Ae  2  e 2  Ae  2 e 2 
 
 
Cancelando términos comunes
i 
qa
i
qa
 iqa 
 me
2
 K 1  2 e
2 2
e 
 
 qa 
 i
qa
i
qa
  m  2 K  cos   
2
 m  K  e  2  e 
2 2 2
 2 
 
eix  eix  2cos x
 Cadena diatómica
 qa 
 M  2 K 1   cos 
2
Ecuación Para M
 2 

 qa 
 2 m  K    cos  Ecuación para m
 2 

Un par de ecuaciones para α y ω como función de q. α esta

dada por

2 K cos(qa / 2) 2K   2 M
 
2K   m
2
2 K cos(qa / 2)
Una ecuación cuadrática para ω2 an se puede obtener de la
siguiente manera.
 2 K cos(qa / 2) 2K   2 M
 
2K   m
2
2 K cos(qa / 2)
qa
4 K cos ( )  4 K 2  2 K 2 ( M  m)   4 Mm
2 2

2
2 qa
4 K (1  cos ( ))  2 K 2 (m  M )   4 Mm  0
2

m  M sin 2
(qa / 2)
  2K (
4
)  4 K
2 2
0
mM mM

K (m  M ) m  M 2 4sin 2 (qa / 2) 1/2

m2  mK [( )  ]
mM mM mM
 Análisis de la relación de dispersión (q)
En la región de longitudes de onda larga (qa«1); sin(qa/2)≈ qa/2

K (m  M ) m  M 2 4sin 2 (qa / 2) 1/2

1,2 2  mK [( )  ]
mM mM mM
K m  M 
12
 m  M  4 qa 
2 2 2
 
2
 K    
mM  mM  mM 4 

K m  M    
12

mM
 1  1  q a  
2 2

mM   m  M  2
 
   
1  x 
12
Usando la expansión de Taylor  1  x 2 para x pequeños
K m  M    mM 2 2 
  2
1  1  k a 
  2(m  M )
2
mM 
 Análisis de la relación de dispersión (q)
Con sinqa«1 Valor máximo para la rama optica

2K  m  M 
max opt 
2

mM
Valor Mínimo para la rama acústica

K  m  M   mMq 2 a 2  Kq 2 a 2
min acus.2   2

mM  2(m  M )  2(m  M )

Sustitutyendo los valores de ω en α (amplitud relativa)

y usando la aprox. cos(qa/2) ≈1 for qa«1 se encuentra para α
2K   2 M

2 K cos(qa / 2)
M
  1 vibracion en fase o    vibracion fuera de fase
m
 Ramas acústica y Optica en (q)

 Hay dos valores de ω para cada valor de q la relación de

dispersión tiene dos ramas ;

Rama Optica
Corresponde al signo +
en 2(q)

Rama Acústica
Corresponde al signo -
q en 2(q)
–л/a 0 л/a 2л/a

La relación de dispersión es periódica en q con periodo

2 π /a = 2 π /(Longitud celda unidad).
Los resultados son válidos para una cadena que contiene un
número arbitrario de átomos por celda unitaria.
 Análisis de la relación de dispersión (q)
 2
Substituya min ac
en la amplitude relative α

 
K(q 2a 2 ) 2K   2 M  1

2

2(m  M)
min ac
2 K cos(qa / 2)
Esta solucion representa ondas sonoras de longitud de onda
larga en el cento de la 1ZB (q=0) los dos tipos de átomos
oscilan con al misma amplitud y fase con una velocidad del
sonido dada por. T

Rama Optica   vs q
1/2
w  K 
vs   a  
Rama Acústica
q  2(m  M ) 
q
–π/a 0 π/a 2π/a
1Zona de Brillouin
Representación esquemática de (q)
 Análisis de la relación de dispersión (q)

Substituya max
2
op
en la amplitud relative para obtener,

2K(m  M) 2K   2 M M
 2
max op    
mM 2 K cos(qa / 2) m
  ck  
M
 m
Rama Optica   vs q Esta solución corresponde al máximo
en la rama óptica, el valor α de
muestra que los dos átomos oscilan in
antifase con el centro de masa en
Rama Acústica reposo
q
–π/a 0 π/a 2π/a
1Zona de Brillouin
 Análisis de la relación de dispersión (q)
El otro limite para la solución de ω(q) es para kq/a = π en el límite
de la 1ZB donde sen(qa/2)=1.
1/2
K (m  M )  M  m  2
4  K (m  M ) mK ( M  m)
 2
max ac  mK     
Mm  Mm  Mm  Mm
ω
Rama Acústica M oscila
2K min 
2K
 2
 y m esta en reposo op
m
max ac
M E  h Brecha de Energía
2K
max ac 
M

2K
 2
min op  Rama Optica m oscila M
esta en reposo.
m
q
π π
 2/a
a
a
 Ramas acústica y Optica en (q)
Para q = 0 en el centro de la 1ZB
rama Acústica ωq0  0 y  =1 Los dos átomos oscilan en fase

2K(m  M)
max 
op
mM
ω

Modo trasversal acústico en una cadena diatómica Rama Optica

rama óptica

2K(m  M) M
max   =-
op
mM m
Los dos átomos oscilan en antifase fase
Rama Acústica

π π q
 2/a
a
a
 Modo transversal óptico en una cadena diatómica

Si los iones tienen carga opuestas, se puede excitar una

vibración de este tipo con campos eléctricos (por ejemplo
una onda electromagnética).

Amplitud de la vibracion exagerada

 Ramas acústica y Optica en (q)
Si hay p átomos en la celda
primitiva, habrán 3p ramas
p=2
ω
LO
+: rama óptica longitudinales
TO2
p-1 longitudinales 1 p-1 óptica
2(p-1) transversales 2 1 acústica TO1

LA
TA2
-: rama acústica transversal TA1
1 longitudinales 1 2(p-1) ópticas
2 transversales 2 2 acústicas
q
3 ramas acústicas y 3 ramas ópticas
Dispersión de Neutrones

El método típico para medir las curvas

de dispersión (q)
 GaN

http://www.ioffe.ru/SVA/NSM/Semicond/GaN/figs/siegle971.gif
Dispersiones de fonones en molibdeno a altas presiones

www.esrf.eu/news/spotlight/spotlight36phonon/
Otras cuasiparticulas en solidos
•Electron: consiste en un electron real acompañado de una nube de carga
opuesta debida a la interaccion con los otros electrons.

 F : 5eV, vF  hk F me : 106 ms 1

•Hueco: similar al electron pero de carga opuesta, y es la ausencia de un

electron en la banda de Valencia por debajo del nivel de Fermi

•Polaron: en cristales polares el movimiento de un electron distorsiona la red de

iones a su alrededor. El polaron es el movimiento del electron mas la distorsion
de la red.

PHYS 624: Crystal Lattice Vibrations: Phonons 42

excitaciones colectivas en
solidos
Estas son bosons y no tienen analogos a particulas reales. Son movimientos
colectivos o coherentes de muchas particulas

•fonon: corresponde al movimiento colectivo cuantizado de todos los atomos de

un solido, lleva momento ћK, y energía ћw h  0.1 eV
•Exciton: estado ligado de un electron y un hueco e2  a  0.1eV
•Plasmon: movimiento colectivo del gas de electrons relative a la red, debido a
la interacioon de coulomb e largo alcance h ne 2 m : 5  20eV
e
•Magnon: movimiento colectivo de espines en una red cristalina

h : 0.001  0.1eV

PHYS 624: Crystal Lattice Vibrations: Phonons 43