Universidad Nacional Mayor de San Marcos
Facultad de Ciencias Fı́sicas Fı́sica del Estado Sólido I
Escuela Académico Profesional de Fı́sica Práctica
1.- Demuestre las ecuaciones de la energı́a potencial, en la aproximación armónica, dadas por las
ecuaciones 22.9 y 22.10 (Ashcroft-Mermin).
2.- Calcule, usando el potencial armónico, la ecuación de autovalores dada por la ecuación 22.53
(Aschroft-Mermin).
3.- Considere una cadena lineal en la que los iones que se alternan tienen masas m1 y m2 y solo
los vecinos más cercanos interactúan a través de un resorte de constante de fuerza K. Encuentre
las relaciones de dispersión para los modos normales. Discuta los casos lı́mites para m1 m2 y
m1 = m2 .
4.-Un ejemplo comúnmente visto de un cristal unidimensional simple es una cadena de masas pun-
tuales, cada una de las cuales tiene dos vecinos más cercanos: uno está a una distancia d y la otra
a una distancia a − d (d a) desde la distancia en el equilibrio. Las dos masas puntuales vecinas
están conectadas por resortes, con constantes de fuerza de los resortes entre las masas puntuales
más próximas K y entre las masas puntuales más distantes G (K > G).
(1) Escriba la energı́a potencial del cristal armónico en términos de los desplazamientos de masas
puntuales desde los sitios equilibrio.
(2) Escriba las ecuaciones de movimiento clásicas para las masas puntuales.
(3) Resuelva las ecuaciones de movimiento y calcule las relaciones de dispersión y los vectores de
polarizacin de los modos normales de las vibraciones de la red.
5.- Un cristal tridimensional tiene una base de dos átomos. Las masas de los dos átomos en la base
son ml y m2 , respectivamente. Sean ~v1 y ~v2 sus velocidades. Demuestre que, para un modo normal
óptico en el centro de la primera zona de Brillouin (~k = 0), m1~v1 + m2~v2 = 0.
Prof. Jaime F. Vento Flores Diciembre del 2019
1
Soluciones de los Ejercicios
1.
a) El teorema de Bloch es independientes de las soluciones de la ecuación de Schrödinger se elige para unposible k
ψ(x + a) = eika ψ(x) (1)
Derivando la ecuación
ψ 0 (x + a) = eika ψ 0 (x) (2)
Para un sólido 1D la función de onda de la estructura de banda ψl (x) es expresado en términos de los coeficientes de
tranmisión y reflexión t y r
a
ψl (x) = eiKx + re−iKx , x≤− (3)
2
La ecuación se modifica como
a
ψl (x) = teiKx , x≥ (4)
2
Para un sólido 1D la función de onda de la estructura de banda ψr (x) es expresado en términos de los coeficientes de
tranmisión y reflexión t y r
a
ψr (x) = te−iKx , x≤− (5)
2
La ecuación se modifica como
a
ψr (x) = e−iKx + reiKx , x≥ (6)
2
Puesto que el hamiltoniano del cristal es idéntico al de un ion en la región −a/2 ≤ x ≤ a/2, cualquier solución de la
ecuación de Schrödinger con energı́a ε debe ser una combinación lineal de ψl y ψr en la región
ψ(x) = Aψl (x) + Bψr (x), −a/2 ≤ x ≤ a/2 (7)
Aplicando las ecuaciones de continuidad a ψ(x) y ψ 0 (x) en x = −a/2, tenemos
ψ(a/2) = eika ψ(−a/2) (8)
h i h i
AteiKa/2 + B e−iKa/2 + reiKa/2 = eika A e−iKa/2 + reiKa/2 + BteiKa/2 (9)
h i h i
A eiKa/2 − eika e−iKa/2 + reiKa/2 + B e−iKa/2 + reiKa/2 − eika teiKa/2 = 0 (10)
Derivando la ecuación
ψ 0 (a/2) = eiKa ψ 0 (−a/2) (11)
h i h i
At(iK)eiKa/2 + B(iK) −e−iKa/2 + reiKa/2 = eika A(iK) e−iKa/2 − reiKa/2 − Bt(iK)eiKa/2 (12)
h i h i
A teiKa/2 − eika e−iKa/2 − reiKa/2 + B −e−iKa/2 + reiKa/2 + eika teiKa/2 = 0 (13)
Introducimos los parámetros α y β,
α = teiKa/2 y β± = e−iKa/2 ± reiKa/2 (14)
Reemplazando
h i h i
A α − eika β+ + B β+ − eika α = 0 (15)
h i h i
A α − eika β− + B −β− + eika α = 0 (16)
2 JFVF-2019
El sistema de ecuaciones tiene soluciones si el determinante es igual a cero
α − eika β+ β+ − eika α
=0 (17)
α − eika β− −β− + eika α
α − eika β+ −β− + eika α − α − eika β− β+ − eika α = 0 (18)
α (β+ + β− ) 1 + 2eika + 2eika β− β+ + 2eika α2 = 0 (19)
(20)
−iKa/2 −iKa 2 iKa
Reemplazando β+ + β= 2e y β+ β= e −r e , tenemos
−α 2e−iKa/2 1 + 2eika + 2eika e−iKa − r2 eiKa + 2eika α2 = 0 (21)
Reemplazando α
2
−teiKa/2 2e−iKa/2 1 + 2eika + 2eika e−iKa − r2 eiKa + 2eika teiKa/2 = 0 (22)
−2t 1 + 2eika + 2eika e−iKa − r2 eiKa 2eika t2 eiKa = 0 (23)
−4t cos ka + 2e−iKa + 2 t2 − r2 eiKa = 0 (24)
t2 − r2 iKa 1
cos ka = e + e−iKa (25)
2t 2t
Por lo tanto, la energı́a del electrón de Bloch está relacionada al vector ~k a través de la ecuación anterior cuando
h̄2 K 2
ε= (26)
2m
b) La ecuación de Schrödinger de una barrera con la misma energı́a es dada por
h̄2 00 h̄2 K 2
− φi + vφi = φi (27)
2m 2m
donde i = 1, 2. La ecuación tiene φ1 y φ2 como soluciones. El Wronskiano w(φ1 , φ2 ) es representadp por
w(φ1 , φ2 ) = φ01 (x)φ2 (x) − φ1 (x)φ02 (x) (28)
Dereivando con respecto a x
d
w(φ1 , φ2 ) = φ001 φ2 + φ01 φ02 − φ01 φ02 − φ1 φ002 (29)
dx
= φ001 φ2 − φ1 φ002 (30)
2mv
= − K 2 (φ1 φ2 (x) − φ1 φ2 ) (31)
h̄2
= 0 (32)
Vemos que w es independiente de x.
c) El wronskiano w(ψl , ψl∗ ) es representado por
w(ψl , ψl∗ ) = iK eiKa − re−iKa e−iKa − r∗ eiKa − eiKa + re−iKa (−iK) e−iKa − r∗ eiKa (33)
= iK eiKa − re−iKa e−iKa − r∗ eiKa + iK eiKa + re−iKa (−iK) e−iKa − r∗ eiKa (34)
= 2iK 1− | r |2 (35)
2
= 2iK | t | (36)
Esta ecuación es válida para la condición x ≥ a/2 y el wronsiano es independiente de x el cual nos da la condición
| r |2 + | r |2 = 1.
d) El wronskiano w(ψl , ψr∗ ) es representado por
w(ψl , ψr∗ ) = iK eiKa − re−iKa t∗ e−iKa − (iK) eiKa − re−iKa = −2iKrt∗ (37)
3 JFVF-2019
Esta ecuación es válida para la condición x ≤ a/2.
para la condición x ≥ +a/2 el wronskiano w(ψl , ψr∗ ) es representado por
w(ψl , ψr∗ ) = iK eiKa + r∗ e−iKa + t∗ (iK) eiKa − r∗ e−iKa t∗ = −2iKr∗ t (38)
Cuando w(ψl , ψr∗ ) es independiente de x se tiene
rt∗ = −r∗ t (39)
(rt ) = −(rt∗ )∗
∗
(40)
es válido solamente cuando rt∗ es puramente imaginario y es de la forma r = ±i | r | eiδ .
e) Asumimos t =| t | eiδ , que nos da r = ±i | r | eiδ .
Ahora, la relación entre t y r es
t2 − r2 | t |2 − | r |2 iδ
= e (41)
t |t|
y el recı́proco de t es
1 1 −iδ
= e (42)
t |t|
Reemplazando en
t2 − r2 iKa 1
cos ka = e + e−iKa (43)
2t 2t
| t |2 + | r |2 iδ iKa 1
= e e + e−iδ e−iKa (44)
2|r| 2|t|
e−i(δ+Ka) 2
= | t | + | r |2 +e−2i(δ+Ka) (45)
2|t|
cos(Ka + δ)
cos ka = (46)
|t|
f) El n-ésimo máximo de | cos(Ka + δ) | para K = Kn es
Kn a + δ = nπ (47)
En el caso de dispersión débil, condición | t |' 1, la banda de enrgı́a prohibida corresponde a un intervalo angosto de
K 0 s de cada valor de Kn .
Por consiguiente, ∆ = K − Kn es la desviación de K de Kn y es representado como
∆ = K − Kn (48)
Hacemos una expasión en series de Taylor cos(Ka + δ) en potencias de ∆
1
cos(Ka + δ) ' 1 − (∆Ka) (49)
2
y cuando la expresión | t | una expasión en series de Taylor en potencias de | r | tenemos
p
|t| = 1− | r |2 (50)
1
= 1 − | r |2 (51)
2
Ahora, el mı́nimo y el máximo de la n-ésima banda prohibida ocurre bajo la condición
| cos((Kn ± ∆K)a + δ) |=| t | (52)
una expasión en series de Taylor
1 1
1− [∆(Kn ± ∆K)a + δ]2 = 1− | r |2 (53)
2 2
∆(Kn ± ∆K)a = |r| (54)
|r|
∆K = (55)
a
4 JFVF-2019
Ahora, la energı́a del gap, εgap , para K = nπ/a y ∆K '| r | /a es representado como
h̄2
εgap = (Kn + ∆K)2 − (Kn − ∆K)2 (56)
m
2h̄2
' Kn ∆K (57)
m
h̄2
εgap ' 2πn |r| (58)
ma2
donde ∆K =| r | /a y K = nπ/a.
g) Ahora, en el lı́mite opuesto es la condición de dispersión fuerte, | t |' 0, la banda de en ergı́a es angosta, y ocurre
próximo de los puntos
1
Ka + δ = n + π (59)
2
Expandiendo en series de Taylor
| Ka + δ | ' ∆Ka (60)
' |t| (61)
El intervalo de la banda de energı́a angosta es representada por
h̄2
εmax − εmin = 2K∆K (62)
2m
= 0|t| (63)
La ecuación de Schrödinger del modelo Kronig Penny, para un potencial tipo delta de Dirac δ(x), es
h̄2 d2 ψ
− + gδ(x)ψ = Eψ (64)
2m dx2
La solución ψ(x) para x < 0 es dado por
ψ(x) = AeiKx + AreiKx , x<0 (65)
iKx
ψ(x) = Ate , x>0 (66)
La función de onda es continua y la condición (1 + r = t) es válida. Integrando la ecuación obtenemos
dψ(0+ ) dψ(0− ) 2mg
− = 2 ψ(0) (67)
dx dx h̄
De la ecuación tenemos
2mgt
iKt − (iK − iKr) = (68)
h̄2
2imgt
t+r−1 = − 2 (69)
h̄ K
img
t 1+ 2 = 1 (70)
h̄ K
1
| t |2 = mg 2
(71)
1+ h̄2 K
1
= (72)
1 + tan2 δ
= cos2 δ (73)
Ahora, la fase cotδ es dado por
<(t)
cotδ = (74)
=(t)
h̄2 K
= − (75)
mg
Para el modelo Kronig Penny