TEORÍA DE BANDAS

FÍSICA DEL ESTADO SÓLIDO

FC-UNI 2021-2
OCTUBRE 2021

ARTURO TALLEDO
DOCTOR EN FÍSICA
 ESTADOS ELECTRÓNICOS EN UN
CRISTAL (TEORÍA DE BANDAS)
• INTRODUCCIÓN
• GAS DE FERMI DE ELECTRONES LIBRES (Kittel)
• TEOREMA DE BLOCH (Aschcroft)
• MODELO DE ELECTRONES CASI LIBRES (Ashcroft)
• MODELO DE ELECTRONES FUERTEMENTE
LIGADOS (Omar)
• MODELO KRONIG PENNEY (Mckelvey)
Valores posibles de energía de un electrón

E (eV) -E (eV)
-13.6

-3.4
-13.6/n2

Valores posibles de energía de un Valores posibles de energía de un

electrón libre (como en un tubo de electrón ligado (en un átomo) Sólo
rayos catódicos o en un TEM) hay valores discretos de energía
constituyen un continuum
 Valores posibles de energía de un electrón en un sólido cristalino

E (eV) -E (eV) -E (eV)

-13.6

-3.4
-13.6/n2

Se sabe que en un sólido cristalino,

un electrón puede tomar valores de
energía contenidos sólo en ciertos
intervalos. Estos intervalos se llaman
Bandas de energía
 Modelos básicos para explicar la existencia de bandas de energía

-E (eV) -E (eV) -E (eV)

-13.6

NFE TBE

Modelo de Modelo de
electrones electrones
casi libres
-3.4
fuertemente
ligados -13.6/n2

Electrones en caja vacía Se sabe que en un sólido cristalino,

(quasicontinuum) un electrón puede tomar valores de
energía contenidos sólo en ciertos
intervalos. Estos intervalos se llaman
Bandas de energía
 ELECTRONES EN UNA CAJA VACÍA

Gas de Fermi de electrones libres, Kittel

Modelo de Sommerfeld de un metal,.. Ashcroft and Mermin,

Física de Sólidos de Brown
 ELECTRONES EN UNA CAJA VACÍA 1D
Cristal 1D, ignorando en primera aproximación el efecto de los átomos en los puntos de red

Ecuación de Schrodinger:
x
V ( x )  0, si x   0, L 
-e V=0
x=0 x=L=Na
 2 d 2
 E
2m dx 2
Soluciones:
 ( x)  ?
 ( x)  A ei q x
E ? Relación de dispersión:
 2q 2
E
2m
 Condiciones de frontera periódicas

x x+L

x=-2L x=-L x=0 x=L x=2L

 ( x)   ( x  L)

iq ( x  L )
Ae iqx
 Ae

2
e i qL 1 q  n, n  0,  1,  2, ....
L
 ELECTRONES EN UNA CAJA 1D, E vs q
Relación de dispersión. Esquema de Zona Extendida  2q 2 2
E q n, n
2m L
E

L  N at a

Nfdo número de
funciones de
onda
disponibles en
cada zona de
.... q
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . Brillouin

4  0 2
 3 2 a 2 3 4
a   L a
a a a a a
SZB PZB SZB
N fdo en PZB  N at en la muestra
EQUIVALENCIA ENTRE REPRESENTACIÓN EN ZONA EXTENDIDA Y EN ZONA
REDUCIDA
Segmentos del gráfico E vs q de todas las ZBs se desplazan a la Primera Zona de Brillouin

E L  N at a

4
a

2
a

2
a

. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

4 3 
2 
a
0 2
a
2 3 4 q
a  L
a a a a a
CZB TZB SZB PZB SZB TZB CZB
ELECTRONES EN CAJA VACÍA 1D

Esquema de Zona Reducida

 ELECTRONES EN CAJA VACÍA 1D
Esquema de Zona Periódica.
E

Tercera banda

Segunda banda
Primera banda


π π k
a a
 ELECTRONES EN CAJA VACÍA 3D
Ecuación de Schrodinger:

V ( x, y, z )  0, si ( x, y, z ) dentro de la caja
y
2 2
    E L2
(x,y,z)

2m V=0
Soluciones: x
L1
 
 (r )  A e iq• r

Relación de dispersión:
2
 q 2
E
2m
 Condiciones de frontera periódicas
(x+L1,y+L2,z+L3)

(x,y,z)

La densidad de nube electrónica en (x,y,z) es igual a la densidad en (x+L1,y+L2,z+L3)

 ( x, y, z )   ( x  L1 , y  L2 , z  L3 )

ei ( q1 L1 )  ei ( q2 L2 )  ei ( q3 L3 )  1

2 2 2
 
q1 , q2 , q3  (
L1
m,
L2
n,
L3
p)
Representación gráfica de E vs q por superficies isoenergéticas
2
2 L1
2 q  2q 2 qy
E 
2m 2m ....................
....................
....................
. . . . . . . . . . . .q. . . . . . . .
 q , q , q   ( 2L m, 2
n,
2
p) . . . . . . . . . . . . .3 . . . . . . .
.. .. .. .. .. .. .. .. .. .. .. .. .. ..E..3 .. .. .. .. ..
1 2 3
1
L2 L3

.................... 2
. . . . . . . . .b . . . . . q. . . . . . L2
L1L2 L3  N at abc . . . . . . . . . . . . . . .2 . . . . .
. . . . . . . . . . . . . . .E.2 . . . .
. . . . . . . . . . . . . .q . . . . . .
8 3 . . . . . . . . . . . . . . 1. . . . . .
. . .. . . . . . . . . . .E1. .a. . . .
N fdo  abc3 . . . a. . . . . . . . . . . . . . . . . qx
8 ....................
N at abc ....................
. . . . . . . . . . . . . . . . . . . .
....................
Nfdo funciones de
. . . . . . . . .b . . . . . . . . . . .
....................
onda espaciales ....................
....................
disponibles en cada ....................
zona de Brillouin
N fdoen PZB  N at en la muestra
 Número de estados cuánticos de un electrón con vector
de onda entre 0 y q :
2 Cada punto representa una función
L1
qx de onda posible para un electrón en
.................... una caja vacía 3D.
.................... Por cada punto hay un paralepípedo
....................
.................... de volumen  2  3 / L1L2 L3
....................
. . . . . . . . . . . . . . . E. . . . . El número de funciones de onda
.. .. .. .. .. .. .. .. .. .. .. .. .. q.. .. .. .. .. .. .. posibles con valores de q entre 0 y
.................... 2
q es:
.................... L2
.................... 4 3
.................... q
.................... qx N func. de onda  3
....................  2 
3

....................
.................... L1L 2 L3
....................
....................
.................... Por cada función de onda posible hay dos
....................
.................... estados cuánticos de un electrón
....................
2 V 4 3
N estados cuánti cos  q
 2  3
3
 ENERGÍA DE FERMI DE UN SISTEMA DE N ELECTRONES
Si de un sistema de N electrones, entran uno por uno todos los electrones a una caja vacía de
volumen V, a 0 K, estos van ocupando sucesivamente los estados cuánticos de menor energía,
dos electrones por cada función de onda espacial. La energía del estado cuántico ocupado por
el último electrón (el de máxima energía) se llama energía de Fermi. La respectiva superficie
isoenergética se llama Superficie de Fermi.
2 En un metal cada átomo contribuye con 1, 2,
qx L1
ó 3 electrones de valencia al sistema de N
.................... electrones.
....................
....................
.................... 2 V 4 3
.................... N qF
..............E . .F . . . .  2 
3
3
.. .. .. .. .. .. .. .. .. .. .. .. ..qF.. .. .. .. .. .. ..
.................... 2
.................... L2
.................... qF  (3  2 n)1/3
....................
.................... qx
.................... 2
.................... EF  (3  2 n) 2/3
....................
.................... 2m
....................
....................
....................
.................... EF es la energía de Fermi y
....................
n es la concentración de electrones
 Densidad de Estados (1)
Número de funciones de onda espaciales de un electrón con vector de onda
entre q y q + dq :
4 q 2 dq
g (q ) dq 
 2 
2 3
L1
qx L1L2 L3
....................
....................
.................... V
. . . . . . . . . . q.+. . . . . . . . . g (q) dq  4  q  dq
2
. . . . . . . d.E. . . . d. q. . . . . . . 2 
3
. . . . . . +. . . . . . . . . . . . . .
E
.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..
. . . . . . . . . . . . . .q. . . . . . 2
.................... L2 Número de estados cuánticos con energía
..............E ...... entre E y E + dE
....................
....................
....................
. . . . . . . . . . . . . . . . . . . . qx
.................... (2)V
.................... g ( E ) dE  4  q  dq
2
....................
2 
3
....................
....................
....................
....................
 Densidad de Estados (2)

De diapositiva anterior

2V 2
g ( E ) dE  q dq
2 2

De la relación E vs q
EF
2m 1/2 1/2
q ( ) E
 2 3/2
V  2m 
g ( EF )  EF    2  EF3/2
2 2  
 2m 
2q dq   dE
2  3/2
   2  N
E 
 3/2
   (3  2 )
F
2 m V
3/2
V  2m  3N
g (E)    E 1/2 g ( EF ) 
2 2   2  2 EF
 TEOREMA DE BLOCH

Las funciones de onda de un electrón sometido a un

potencial periódico se identifican con un vector de onda k en
la primera zona de Brillouin
 Teorema de Bloch
• Las funciones de onda espaciales de un electrón en un
cristal (potencial periódico) son de la forma:
   
  (r )  e i k • r u  (r )
nk nk
   
donde: u  ( r  R )  u  (r ) y k  PZB
nk nk

y
(x,y,z)
L2
V=0
x
L1
 Teorema de Bloch
• Las funciones de onda espaciales de un electrón en un
cristal (potencial periódico) son de la forma:
   
  (r )  e i k • r u  (r )
nk nk
   
donde: u  ( r  R )  u  (r ) y k  PZB
nk nk

PRUEBA:
2 2 
En la ecuación de Schrodinger     V (r )  E  ...(1)
2m
   
V ( r  R )  V ( r )  V ( r )   Vg e • r
 
i g ...(2)

g


CFP   (r )   cqe • r
 
i q
...(3)

q
 Teorema de Bloch
2 2 
    V (r )  E (1)
2m 2 2 2


   c q 2e i q .r
(1a)
2m 2m q
 q


 (r )   cqei q • r
 


2

q 
E c  E e i q .r
(1b)
2m q
 q

V ( r )   Vg e • R
 
i g

g
       
V (r ) ( r )   Vg e   cq e . r 
i g . r i q

 g  q 

 
V (r ) (r )  Vgcq e V
    
i ( g  q ). r
 
g cq' g e i q '. r
 
gq g q'

 
V (r ) (r )  Vg' cq g' e
 
i q. r

g 'q (1c)
 Teorema de Bloch
Reemplazando (1a) y (1b) y (1c) en (1) tendremos entonces la ecuación de
Schrodinger expresada en la forma de un sistema infinito de ecuaciones finitas:

 2 2 
e  V  0

iq .r
 q  E  cq   cq g '

q  2m  
g'
g'

(Una ecuación para cada vector q posible)

 
Puesto que las funcionese i q .r son linealmente independientes, tendremos para cada
q una ecuación así:

 2 2 
 q  E  cq  V 
g' cq g'  0 Sistema de Ecuaciones A
 2m  g'
 Teorema de Bloch
Cualquier q siempre puede ser escrito como k - g donde k está en la PZB
y g es un vector de la red recíproca. Luego:

 2   2 

 2m
( k  g )  E  ck g 

V 
g'

g'
c  
k g g '
0

  
Cambiando el nombre de los índices:   g   g
g

 2   2 

 2m
( k  g )  E  ck g 

V

g'
 
g ' g
c 
k g '
0
Subsistema de Ecuaciones B

Enfatizamos que ésta es simplemente la ecuación de Schrodinger en el

espacio recíproco.
Hay N valores posibles de k en la PZB
Por cada k en la PZB podemos escribir un subsistema de ecuaciones del tipo
(subsistema de ecuaciones B)
Las ecuaciones con un k determinado no se mezclan con las ecuaciones para otros
k
 2   2 
 ( k  g )  E  ck g   Vg'  g c    0
 2m 1
 1 
g' 1 k g '

 2   2 

 2m
( k  g 2
)  E  ck g 
 2
V

g'
 
g '  g2
c 
k g '
0

Subsistema de Ecuaciones B
   2
2


 2m
( k  g 3
)  E  ck g 
 3
V

g'
 
g '  g3
c 
k g '
0

. . . . . . . . . . . . . .0

 2   2 

 2m
( k  g n
)  E  ck g 
 n
V

g'
 
g 'gn
c 
k g '
0
   2   2  
 2m (k  g1 )  E   Vg  g  Vg  g  ...Vg  g 
  
1 1 2 1 n 1

  ck g1 
 2   
  (k  g ) 2  E    
Vg  g  Vg  g  ...Vg  g  
 2m 2 
 2 
c

  k  g2 
1 2 2 2 n
 

  
  2     
 (k  g3 ) 2  E   Vg  g  Vg  g  ...Vg  g  c     0
 Subsistema de Ecuaciones B
 2m  1 3 2 3 n 3
  k  g2 
  
   
. . . . .  ... 
  
  
  2   2   ck  gn 
 (k  g n )  E   Vg  g  Vg  g  ...Vg  g   
 2m  
1 n 2 n n n

 
Haciendo el Determinante = 0 tendremos los valores propios En y para cada valor propio un vector
propio: (ck-g1, ck-g2, ck-g3,…….. ck-gn). Estas componentes del vector propio son al mismo tiempo los
coeficientes de la función de onda, es decir:

 
  (r )   c  e  e  c  e
    
 e u  (r ). Lqqd

i ( k  g )• r i k . r  i g • r i k . r
nk  k g  k g nk
g g
 Teorema de Bloch: CASO DEL CRISTAL 1D
q1 q2 qN qN 1 q2 N 1

-4π -3π -2π -π π 2π 3π 4π q

a a a a 0 a a a a

 2 2 
 q1  E  c q1  V g' c q1g '  0
 2m  g'

Resolver la ecuación de Schrodinger  2 2 

se convierte en el problema de resolver un  q 2  E  c q2  V g' cq2 g '  0
sistema de infinitas ecuaciones lineales  2m  g'

finitas, una por cada q posible entre 

 2 2 
y  .  q N  E  c qN  V g' c q N g '  0
 2m  g'

 2 2 

2 m
q N 1  E  cqN 1  V g' cq N 1g '  0
  g'

 2 2 

2 m
q2 N 1  E  c q2 N 1  V g' c q2 N 1 g '  0
  g'
 Teorema de Bloch: CASO DEL CRISTAL 1D
Re enumerando los valores de los q posibles: q = k+g donde k es un q posible en la PZB
Y g es un vector de la red recíproca.
q2 qN qN 1
q1
k1  0 k1  g1 k1  2 g2

-4π -3π -2π -π 0 π 2π 3π 4π

q
a a a a a
a a a

q1  k1  0

q2  k2  0

2
qN 1  k1  g1  k1 
a

4
q2 N 1  k2  2 g1  k1 
a
 Teorema de Bloch: CASO DEL CRISTAL 1D
Con este cambio de nomenclatura, el sistema de ecuaciones lineales equivalente a la Ecuación
de Schrodinger puede reescribirse en la forma:
q1 q2 qN qN 1
k1  0 k1  g1 k1  2 g2

-π π 2π 3π 4π
-4π -3π -2π 0 a a q
a a a a a a

 2 
 (k1  0) 2  E  c k1 0  V g' c k1 0g '  0
 2m  g'

 2 
 (k2  0) 2  E  ck2 0  V g' c k 2  0 g '  0
 2m  g'

 2 

2 m
( k1  g1 ) 2
 E  c k1  g1  V g' c k1  g1 g '  0
  g'

 2 
 ( k2  g1 ) 2  E  c k2  g1  V g' c k2  g1 g '  0
 2m  g'

 2 

2 m
( k1  2 g1 ) 2
 E  c k1 2 g1  V g' c k12 g1g '  0
  g'
 Teorema de Bloch: CASO DEL CRISTAL 1D
Hecho esto podemos ver que el conjunto infinito de ecuaciones lineales se puede
reescribir como un conjunto de subsistemas, uno para cada k posible de la PZB. Los
subsistemas son independientes entre sí. Son de la forma:
k1  0 k1  g1 k1  2 g2

-π π 2π 3π 4π
-4π -3π -2π 0 a a q
a a a a a a

 2   2 
 ( k1  0) 2
 E  c k1 0  V g' ck10g '  0  ( k  0) 2
 E  c k 2 0  V c k 2 0  g '  0
 2m
2 g'
 g'  2m  g'

 2   2 
 ( k1  g1 ) 2
 E  c k1  g1  V g' c k1 g1g '  0 
2 m
( k 2  g1 ) 2
 E  c k2  g1  V g' ck2  g1 g '  0
 2m  g'   g'

 2 
 2
 (k2  2 g1 ) 2  E  ck2 2 g1  V ck2 2 g1 g '  0
 ( k1  2 g1 ) 2
 E  c k12 g1   Vg' ck12 g1g'  0 
 2m  g'
g'

 2m  g'

El potencial es dato, es decir todas las Vg son datos. Las incógnitas son las E y las f.d.o., es decir, los ck+g.
 Cada sistema lineal de ecuaciones asociado a un k de la PZB se puede
escribir en forma matricial como se muestra aquí:
  2   Ck1 0 
 ( k1  0)  E  Vg1 Vg 2 . . .  
2

C
 2m    k1  g1 
  2   Ck  g 
(k1  g1 ) 2  E  Vg 2 . . .  
2 
 Vg1  1

  2m   .  0
  . 
 Vg Vg . . .   (k  g ) 2  E  . . .  
2

 1 2  1 n   C 
 2m   
k1  g n

. . .  . 

valores de energía se obtienen igualando a cero el determinante:
  2  
 ( k1  0) 2
 E  Vg 1 Vg 2 .. . 
 2 m  
  2  
 Vg1  ( k1  g1 )  E  Vg 2 .. . 
2

  2m   0
 
 Vg Vg . . .  2

( k1  g n ) 2  E  . . .

 1 2  
 2m 
 
.. . 
Y las funciones de onda asociadas con cada valor de energía son combinaciones lineales
de la forma:

 nk ( x)   ck  g ei ( k  g ) x  ei k x  ck  g e  i g x  ei k xuk ( x )

g g

n es el número de la zona de Brillouin o número de

banda
 DEL MODELO DE ELECTRONES CASI LIBRES AL DE CAJA VACÍA

 2 
 (k  g) 2  E  ck g  V g ' g ck g '  0 Subsistema de Ecuaciones B
 2m  g'

Potencial periódico nulo, todos los coeficientes Fourier de potencial son Vg = 0

 2 
 (k  g)  E  c k g  0
2
Sistema trivial de ecuaciones, cada ecuación tiene un solo
 2m  término.

 2 

 2m
(k ) 2  E  ck  0
 c k -g  1, si q  k - g
 2  c k -g '  0, si q  k - g'
 (k g ) 2  E  ck g  0
 2m 
2 2 2 2
E  E (0) , E k -g , E q
 2  k-g
2m 2m
 (k 2g ) 2  E  ck g  0
 2m 

 2

(k 3g ) 2  E  ck g  0
q (r )  Aeiq.r   (r )  Ae
 
e
i k .r ig .r
 nk
 2m 
 MODELO DE ELECTRONES CASI LIBRES:
Al potencial V=0 de un electrón en una caja vacía se le agrega un potencial
periódico débil
En este modelo se resuelve la ecuación de Shrodinger para un electrón en un
potencial periódico, resolviendo el subsistema de ecuaciones B pero con la
condición de que el potencial periódico es débil.

El modelo asume que las nuevas funciones de onda y valores de energía van a
ser las mismas funciones de onda y valores de energía de caja vacía sólo un
poquito perturbadas por el potencial periódico. Es decir, una función de onda que
en caja vacía era
     2   2
 (r )  A e i ( k  g1 )• r
con energía E   
(0) k  g1
k  g1 2m
Se convertirá debido al potencial periódico en una función de onda

  (r )   c  e  g )•r con c 
  
i ( k
  1 y todos los otros c   0
nk  k g k  g1 k  g1
g

y energía incógnita E
 El modelo distingue dos tipos de soluciones:

1) Para k-gi individuales que están lejos de una frontera de zona de

Brillouin (CASO NO DEGENERADO)
2) Para conjuntos de k-gi que están cerca a una frontera de zona de
Brillouin CASO CASI DEGENERADO.

En el caso no degenerado se comienza a resolver dividiendo el subsistema de

ecuaciones B en dos:  2 2
(i)la ecuación individual del sistema que contiene la energía a modificar E (0)  k  g1
k  g1 2m

(ii) Todos los demás términos correspondientes del subsistema de ecuaciones B

En el caso casi degenerado se comienza a resolver dividiendo el subsistema de

ecuaciones B en dos:
(i)Las m ecuaciones del subsistema que contienen valores iguales o casi iguales de
energía y que estando separados por un vector g de la red recíproca, están ambos muy
cerca de m fronteras de zona de Brillouin, m = 2 en el caso 1D
puede ser m=4 si el cristal 2D de red cuadrada o un número mayor en cristales
3D
(ii) Todos los demás términos correspondientes del subsistema de ecuaciones B
 MODELO DE ELECTRONES CASI LIBRES
Caso no degenerado ó k-g Lejos de las fronteras de zonas de Brillouin
(0)
E k-g 1
- E (0)
k -g  V, g  g1

E  E  c (0)
k - g1 k g 1   Vg g 1 ck g Ecuación B1
g

E  E  c
(0)
k -g k g  Vg1 -g ck - g1  V g 'g c k g ' Ss.Ecs B2
g '  g1

Vg1 -g ck - g1 Vg ' g ck g '

ck g 
E  E   E  E 
(0)
k -g

g '  g1
(0)
k -g

Puesto que V es pequeño, esperamos: ck -g1  1, y ck -g  0, si g  g1

Vg1 -g c k - g1
c k g   0(V 2 )
E  E  (0)
k -g
Ecuación C
 MODELO DE ELECTRONES CASI LIBRES

Caso no degenerado:
(0)
E k-g 1
- E (0)
k -g  V, g  g1

Reemplazo ecuación C en ecuación B1:

Vg g 1 Vg1 g
E  E  c
(0)
k - g1 k g 1 
E  E k(0)-g
ck g1  0(V 3 )
g

Vg g 1 Vg1 g
EE (0)
k - g1   0(V 3 )
g E  E k(0)-g

Los niveles no degenerados no cambian apreciablemente, no forman parte del

fenómeno de aparición de bandas de energía que queremos explicar.
 MODELO DE ELECTRONES CASI LIBRES
Caso Casi Degenerado ó k-g cerca de la frontera de zona de Brillouin:
Primero escribimos las ecuaciones del subsistema que
involucran a los k-gj cuya modificación queremos (0)
E k-g i
- E (0)
k-g j  V, si i, j  1, m.
estudiar, que están cerca a la FZB y que satisfacen:

E  E  c
m
(0)
k -g i k g i   Vg j g i ck g j  V g g i c k g , i  1, 2, ... m Ss.ecuaciones B1
j1 g  g1 , g 2 ,..g m

Luego escribimos las ecuaciones del subsistema que

involucran a los k-g cuyos valores de energía de caja vacía
(0)
E k-g i
- E (0)
k-g  V, si g  g1,g 2 ,...g m
están muy lejos de los gi anteriores,es decir, satisfacen:

E  E  c
m
(0)
k -g k g   Vg j g c k g j  V g ' g ck g ' , g  g1 , g 2 ,... g m Ss.ecuaciones B2
j1 g ' g1 , g 2 ,..g m

Hay que resolver B1 y B2 considerando ck-gi >> ck-g

m
1
ck g 
 E  E k(0)-g   g1 -g k - g1
V
j 1
c  0(V 2
)
Ecuación D
 MODELO DE ELECTRONES CASI LIBRES

Caso Casi Degenerado:

E  E  c
m
(0)
k -g i k gi   Vg j g i c k g j  V g g i c k g , i  1, 2, ... m Ecuación B1
j1 g  g1 , g 2 ,..g m

m
1
c k g 

E  E k(0)-g   V
j 1
c
g1 -g k - g1  0(V 2 ) Ecuación D

Reemplazo la ecuación D en B1:

Vg g i Vg j g
E  E  c
m m
  Vg j g i c k g j   (  E  E  )ck g j  0(V 3 ) , i  1, 2, ... m
(0)
k -g i k g i (0)
j1 j1 g  g1 , g 2 ,..g m k -g
 MODELO DE ELECTRONES CASI LIBRES
Caso Casi Degenerado:
Las nuevas energías de las funciones de onda con k-gi cerca a las FZB, respecto al
potencial cero, se calculan entonces resolviendo el siguiente «sistema de ecuaciones
del caso casi degenerado» :
E  E  c
m
(0)
k -g i k g i   Vg j g i ck g j , i  1, 2, ... m SECCD
j1

Matricialmente, por ejemplo, para m = 4:

E  E k(0)-g 1 - Vg 2 g 1 - Vg 3 g 1 - Vg 4 g 1  ck -g 
  1 
- Vg1 g 2 E  E k(0)-g 2 - Vg 3 g 2 - Vg 4 g 2  ck -g 2 
    0
- V
 g1 g 3 - Vg 2 g 3 E  E (0)
k -g 3 - Vg 4 g 3  ck -g 3 
 (0)  c
- Vg1 g 4 - Vg 2 g 4 - Vg 3 g 4 E  E k -g 4   k -g 4 

El resultado importante es que se rompe la cuasicontinuidad de la relación de

dispersión de caja vacía originándose «intervalos de energía prohibidos».
El orden de degeneración m depende del tipo de estructura (de red)
y también de la posición de un k dentro de la zona de Brillouin.

qy PZB

k1

 k2
a k2
 
 k1
a a
PZB qx


a

Para una red recíproca cuadrada

k1 es doble degenerado Para una red recíproca sc
k1 es 8- degenerado
k2 es 4-degenerado
k2 es 2-degenerado
MODELO DE ELECTRONES CASI LIBRES: El Cristal 1D
Caso Casi Degenerado: frontera entre PZB y SZB

De SECCD:
E
EE  c
2
(0)
k -g i kg i
 V j 1
g j g i
ck  g , i  1, 2
j

EE  c
(0)
k -g 1 k  g1
 Vg  g ck  g  Vg  g ck  g
1 1 1 2 1 2

EE  c(0)
k -g 2 k  g2
 Vg  g ck  g  Vg
1 2 1 2
g 2
ck  g
2

2
g1  0, g2  g0 
a

Eg
EE  c (0)
k k
 V0 ck  Vg ck  g
0 0



a
k  2
a a
kk EE  c (0)
k g 0 k  g0
 V g ck  V0 ck  g
0 0
 EE  c (0)
k k
 V0 ck  Vg ck  g
0 0

EE  c (0)
k g 0 k  g0
 V g ck  V0 ck  g
0 0

Consideramos V0 = 0, nivel de
referencia de energías  E  Ek(0) - Vg0  ck 
 
c   0
-V g0 E  Ek(0)
 g0   k  g0 

Ek(0)  Ek(0)
 g0

 E0  E0 2  4 V 2 
 k  g0 k g0 


E 
2 2

Cuando k=g0/2 E  Ek(0)  Vg0 El gap en k=g0/2

E  Ek(0)  V 2 E g  2 V 2
( ) ( )
a a
MODELO DE ELECTRONES CASI LIBRES:El Cristal 1D
Caso Casi Degenerado: frontera entre SZB y TZB
E De SECCD: cuando las f.d.o casi degeneradas son k  g0 y k  g0

EE  c
2
(0)
k -g i k g i
 V
j 1
g j g i
ck  g , i  1, 2
j

EE  c
(0)
k -g 1 k  g1
 Vg  g ck  g  Vg  g ck  g
1 1 1 2 1 2

Gap 2
EE  c(0)
k -g 2 k  g2
 Vg  g ck  g  Vg  g ck  g
1 2 1 2 2 2

2 2
k  g0 g1  g0  , g2   g0  
k  g0 a a
-2π
a 
π
k
π
2π k
a
k EE  c (0)
k -g 0 k  g0
 V0 ck  g  V2 g ck  g
0 0 0

EE  c
a a (0)
kg 0 k  g0
 V2 g ck  g  V0 ck  g
0 0 0
EE  c  V c
(0)
k -g 0 k  g0 0 k  g0
 V2 g ck  g
0 0

EE  c  V
(0)
kg 0 k  g0 2 g0
ck  g  V0 ck  g
0 0

 E  Ek(0)
 g0 - V2 g0  ck  g0 
    0
-V2 g0 E  Ek(0)
 g0
c
  k  g0 


 E0  E0   4 V2 g0 
2 2

E (0)
k  g0 E (0)
k  g0  k  g0 k  g0

E 
2 2
Cuando k= 0
E  Eg(0)0  V2 g0 El gap en k=g0

E  E g(0)  V 4 E g  2 V 4
0 ( ) ( )
a a
BANDAS DE ENERGÍA DE UN CRISTAL 1D SEGÚN MODELO NFE

L  Na

E
4
a

Eg3=2Vg3
2
a

2
a
Eg2=2Vg2

Eg1=2Vg1
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .

4
3 
2 
a
0 2
a
2 3 4 q
a  L
a a a a a
CZB TZB SZB PZB SZB TZB
CZB
 Interpretación de Bragg

Cerca de la frontera de zona de Brillouin: λ  2a

2dsenθ  n

La superposición de ondas incidentes y reflejadas produce ondas

estacionarias seno y coseno con diferentes valores de energía


 ( x)  B Cos x
a


 ( x)  BSen x
a
 MODELO DE ELECTRONES
FUERTEMENTE LIGADOS

Los niveles de energía de átomos libres se convierten en

bandas de energía cuando un sistema de N átomos
distribuidos en una red cristalina de parámetros de red muy
grandes se aproxima a un estado de equilibrio
Evolución de la función de onda de un electrón en el potencial de dos iones
positivos
La molécula H 2 : -
+ Niveles de energía
+

degeneradas

no degeneradas
 N núcleos y un solo electrón
-
+ + + + + + + + +

N funciones de onda para N

niveles de energía

N funciones de onda para

un solo nivel de energía
 MODELO DE ELECTRONES
FUERTEMENTE LIGADOS
Cuando hay N iones y un solo electrón cada nivel atómico da
lugar a N niveles de energía, los cuales constituyen una banda
de energía.
MODELO DE ELECTRONES
FUERTEMENTE LIGADOS

a) Potencial periódico
b) Funciones de onda atómicas
c) Funciones de onda atómicas multiplicadas por onda plana
 MODELO DE ELECTRONES
FUERTEMENTE LIGADOS
Escojamos una función de onda que satisfaga el teorema de Bloch y que
se reduzca a una función de onda atómica para valores de x cerca de cada ion

N
1
e
i kX j
 k (x)  1/2  ν (x  X j ) Ec (1)
N j1

Satisface el teorema de Bloch

1 ikx N -i k( x -X j )
 k (x)  1/2 e  e  ν (x  X j ) Ec (1a)
N j1

Cerca del j-ésimo ion se comporta como un orbital atómico:

i kX j
 k (x)  e  ν (x  X j )   ν (x  X j ) Ec (1b)
 MODELO DE ELECTRONES FUERTEMENTE LIGADOS
Cálculo del valor esperado de la energía
E (k )   k
 * H  dx
k

E (k )   k H  k Ec (2)
Reemplazando Ec (1) en Ec (2):
1 N i k ( X j  X 'j )
E (k )   e  ( x  X j ' ) H  ( x  X j )
N j, j '

E (k ) 
N
e 
1 i k ( X j  X 'j )
 ( x  X j ' ) H  ( x  X 1 )   ( x  X 2 )   ( x  X 3 )  ...  (S)
Asumimos, despreciando efectos de frontera, que para cualquier j´ la suma S da lo mismo
N 1
j
2
e
i k Xj
E (k)   ν (x) H  ν (x  X j ) Ec (3)
N
j 
2
 MODELO DE ELECTRONES FUERTEMENTE LIGADOS
Cálculo de la energía
 N 1 /2
E (k )  
j  N /2
ei k X j  ( x) H  ( x  X j ) Ec (3)

E (k )   ( x) H  ( x)   ' ei k X j  ( x) H  ( x  X j ) Ec (4)
j

2 d 2
H  -  V ( x)
2 m0 dx 2

V ( x)   v (x - X
j
j
)

V ( x)  v( x)  V ' ( x)
 MODELO DE ELECTRONES
FUERTEMENTE LIGADOS
Primer término de la ecuación (4)

 2 d 2 
 ν (x) H  ν (x)   ν (x)  2
 v( x)  ν (x)   ν (x) V' ( x)  ν (x)
 2m0 dx 

β  -   (x) V' ( x)  ν (x)


*

 ν (x) H  ν (x)  E   Ec (4a)

 MODELO DE ELECTRONES
FUERTEMENTE LIGADOS
Segundo término de la ecuación (4), considerando sólo vecinos más cercanos


' i k Xj
e  ν (x) H  ν (x  X j )  eika  ν (x) H  ν (x  a)  e -ika  ν (x) H  ν (x  a)
j

 2 d 2 
 ν (x) H  ν (x - a)   ν (x)  2
 v ( x  a )   ν (x - a)   ν (x) V' ( x  a)  ν (x)
 2m0 dx 

  -  * (x) V' ( x  a)  ν (x - a) Integral de traslape

E (k )  E    2  Cos ka Ec (5)
 MODELO DE ELECTRONES
FUERTEMENTE LIGADOS

ka
E (k )  E 0  4  sen ( )
2

E0  E    2 

Ek  E0  γ a k 2 2
 MODELO DE ELECTRONES
FUERTEMENTE LIGADOS
Caso 3 dimensional: banda s de un sólido con estructura fcc monoatómica

e
ik  R 0j
E s (k)  E s   - γ
j, vecinos más cercanos

a
 2 (i  j)

a
k  (k x , k y , k z ) y R 0   (k  j)
2
a
 2 (i  k )

  k xa   k ya   k ya   k za   k za   k x a 
E s (k )  E s    4 Cos  Cos   Cos  Cos   Cos  Cos 
  2   2   2   2   2   2 

E s (k )  E s    12  γ k a 2 2
 MODELO KRONIG PENNEY
0, si x   0, b v

V(x)  - V0 , si x   b, a 
V(x)  V(x  n a)

V
 2 d 2
bb
-  V(x)  E 
2 m dx 2 x
V0
ac
a

d 2 d 2
  2   0 , si x   0, b   2   0 , si x   b, a 
dx 2 dx 2

1/2 1/2
2mE  2 m (E  V0 ) 
α  β 
 2   
    2 
 V
b

V0
a
 MODELO KRONIG PENNEY
d 2 d 2
  2   0 , si x   0, b   2   0 , si x   b, a 
dx 2 dx 2

Ψ k (x)  e ikx u k (x)

Ψ k (x)  e u1k (x)

ikx
Ψ k (x)  e ikx u 2k (x)

d 2u du
1
 2i k 1
 (k 2  α 2 ) u 1  0 , si x   0, b
dx 2 dx

d 2u du
2
 2i k 2
 (k 2  β 2 ) u 2  0 , si x   b, a 
dx 2 dx
 MODELO KRONIG PENNEY
d 2u du
1
 2i k 1
 (k 2  α 2 ) u 1  0 , si x   0, b
dx 2 dx

d 2u du
2
 2ik 2
 (k 2  β 2 ) u 2  0 , si x   b, a 
dx 2 dx
Caso 1: Si E>0

u1 (x)  A e i(α - k)x  B e -i(α  k)x , si x   0, b

u 2 (x)  C e i(β-k)x
 De -i(β  k)x
, si x   b,a 
 MODELO KRONIG PENNEY
Usando como condiciones de frontera la continuidad de la función y
su derivada en x = 0 y en x = b

ABCD  0

i α - k  A  i  α  k  B  i  β - k  C  i  β  k  D  0

A e i(α - k)b  B e -i(α  k)b  C e -i(β - k)c  De i(β  k)c  0

i α - k  A e i(α - k)b  i  α  k  B e -i(α  k)b  i β - k  C e -i(β - k)c  i β  k  De i(β  k)c  0

 MODELO KRONIG PENNEY
Los coeficientes A,B, C y D se encuentran resolviendo estas ecuaciones
simultáneas y la solución existe sólo si el determinante de sus
coeficientes es nulo, lo cual implica:

 α 2  β 2 
 
 Sen α b Sen β c  Cos  b Cos  c  Cos k (b  c) (Sol.1)
2α β
 MODELO KRONIG PENNEY
Caso 2: -V0 < E < 0

El procedimiento es igual pero α es imaginario y podemos escribirlo así:

α  iγ
Y teniendo en cuenta que:

Cos (i x)  Cos h x y Sen (i x)  i Senh x

Puede deducirse que habrá solución no trivial para A,B,C y D sólo si:

 γ 2  β 2 
 
Senh γ b Sen β c  Cosh γ b Cos β c  Cos k a .... (Sol. 2)
2γ β
 MODELO KRONIG PENNEY
Para resolver exactamente E vs k puede procederse del siguiente modo:

1) Asegurarse de tener los datos V0 (en eV), a y c (en angstron) .

2) Tomar 2000 valores de energía E en el intervalo (- V0 , V0) y para cada

y(α ó
β γ
uno de ellos escribir los valores de ).
β α γ
3) Reemplazar y( ó ), así como a, b y c en la solución (1) o

solución (2). De este modo se obtiene f (E).

4) Para los valores de f(E) en el intervalo (-1,1) tome arco cos f ( E ) y

así podrá obtener k. Con lo cual habrá obtenido un punto de la relación de
Dispersión E vs k).
MODELO KRONIG PENNEY
Otros métodos de cálculo de bandas

El método celular

El método de ondas planas aumentadas

El método del seudo potencial

 EL MÉTODO CELULAR

Resuelve la ecuación de schrodinger en una celda Wigner Seitz

con condiciones de frontera:

 (r)  e -ik.R  (r  R)

 
nˆ (r )   (r)  -e -ik.R nˆ (r  R)   (r  R)

En la publicación original para calcular el estado de más

baja energía de la banda 3s del sodio se aproximó la celda
WS a una esfera y se resolvió el problema como si fuera un
problema atómico pero con condiciones de frontera :

 0, (r0 )  0
SUPERFICES DE FERMI
 SUPERFICIES DE FERMI DE ELECTRONES en caja vacía
Si de un sistema de N electrones, entran uno por uno todos los electrones a una caja vacía de
volumen V, a 0 K, estos van ocupando sucesivamente los estados cuánticos de menor energía,
dos electrones por cada función de onda espacial. La energía del estado cuántico ocupado por
el último electrón (el de máxima energía) se llama energía de Fermi. La respectiva superficie
isoenergética se llama Superficie de Fermi.
2 En un metal cada átomo contribuye con 1, 2,
qx L1
ó 3 electrones de valencia al sistema de N
.................... electrones. N = Nat, N = 2Nat ó N = 3 Nat
....................
....................
.................... 2 V 4 3
.................... N qF
..............E . .F . . . .  2 
3
3
.. .. .. .. .. .. .. .. .. .. .. .. ..qF.. .. .. .. .. .. ..
.................... 2
.................... L2
.................... qF  (3  2 n)1/3
....................
.................... qx
.................... 2
.................... EF  (3  2 n) 2/3
....................
.................... 2m
....................
....................
....................
.................... EF es la energía de Fermi y
....................
n es la concentración de electrones
SUPERFICIES DE FERMI DE ELECTRONES en caja vacía

2 VPZB  VPr imrec

qx L1

....................
.................... 8 3
.................... VPZB 
.................... VPr imdirecta (1)
....................
....................
.. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. .. ..
.................... 2  2 
3

.................... L2 Vcubitorec 
.................... L1L2 L3
....................
.................... qx
.................... 8 3
.................... Vcubitorec 
.................... Vmuestra
....................
....................
....................
.................... Vmuestra  N ´VPr imdirecta(2)
....................
....................

(2) en (1) VPZB  N ´Vcubitosenredrecíproc

#puntos en PZB = N´= #celdas primitivas en red directa

 SUPERFICIES DE FERMI DE ELECTRONES en caja vacía

2
Si el material es monoatómico
qy L1

....................
....................
les
eta

.................... #puntos en PZB = N’= Nat = # átomos en la

m
ra

....................
ES p a

muestra
. . . . . . . SZB .............
NT mi
LE er
A eF

....................
IV a d

.. .. .. .. .. Esfera
.. .. de.. Fermi
.. .. para
.. ..metales
.. .. .. .. .. .. .. ..
TR sfer

#Estados cuánticos en PZB = 2Nat

E

....................
DIVALENTES
2
.................... L2
.................... N= número de electrones
. . . . . . .Esfera
SZB . .de.Fermi. .para. . . . . . . .
. . . . . . .MONOVALENTES
.............
metales
qx
. . . . . . . . . . . . . . . . SZB ....
....................
.. .. .. .. .. .. .. .. .. .. PZB.. .. .. .. .. .. .. .. .. .. Para átomos monovalentes:
.................... N = Nat
.. .. .. .. .. .. .. .. .. .. SZB.. .. .. .. .. .. .. .. .. ..
.................... Para átomos divalentes:
....................
N = 2Nat

En cada punto «entran» 2 electrones Para átomos trivalentes:

N = 3Nat
 SUPERFICIES DE FERMI
La superficie de Fermi es definida como la superficie en el espacio k dentro de la
cual todos los estados están ocupados por electrones de valencia a 0K. Todos los
estados fuera de la esfera de Fermi están vacíos.

ky

a


a kx

En el modelo de electrones La superficie de Fermi real

en caja vacía, la superficie puede tener diferentes formas
de Fermi es la esfera de manteniendo las propiedades
de simetría de la red.
Fermi
SUPERFICIES DE FERMI DE ELECTRONES en caja vacía
Método de Harrison: Una técnica para construir superficies de Fermi de metales.
En primera aproximación consideramos que los electrones son libres. Las superficies
De Fermi son esferas. Si salen de la PZB usamos dos o mas zonas de Brillouin.
SUPERFICIES DE FERMI DE ELECTRONES en caja vacía

3ZB

1ZB
2ZB
 SUPERFICIES DE FERMI DE ELECTRONES en caja vacía
Método de Harrison para materiales trivalentes
 SUPERFICIES DE FERMI DE ELECTRONES en caja vacía

PROBLEMA: Considere un sólido cristalino monoatómico 2D con red cuadrada de

parámetro de red a y átomos divalentes. Usando el modelo de electrones en caja
vacía encuentre la fracción de electrones que ocupan funciones de onda con valores
de q en la segunda zona de Brillouin.

Relación entre q de Fermi y concentración de electrones:

 qF2
2 2
N
 2 
 
 L 

qF2  2 n 2 2 
q  2
2
qF 
F
a2 a
 SUPERFICIES DE FERMI DE ELECTRONES en caja vacía
2
a

2  cos  0 
qF  2
a qF
De figura:  0  27.5971

 /a
cos  0   0  0.4816 rad
qF
Sen2 0  0.82109

qF

4 Asomb
f SZB  Asomb   0 (rad )  qF2  q F2 sin  0 cos  0 f SZB  0.09
 qF2
 SUPERFICIES DE FERMI DE ELECTRONES en caja vacía

Problema 5.14 de Omar: Recordando que para un gas de electrones libres:

qF  (3  2 n)1/3

a) Cuando la concentración de electrones aumenta, la esfera de Fermi se expande.

Demostrar que esta esfera empieza a tocar las caras de la PZB en una red fcc cuando
La relación entre concentración de electrones y de átomos es n/n a = 1,36
b) Suponga que algunos de los átomos en un cristal de cobre que tiene estructura
fcc monoatómica son gradualmente reemplazados por átomos de zinc.
Considerando que el zinc es divalente y el cobre monovalente, calcular la
razón de
átomos de cobre a átomos de zinc en una aleación CuZn (bronce) a la cual la
esfera de Fermi toca las caras de la zona. ( El bronce experimenta un cambio
estructural a esta razón de concentraciones)
La recíproca de una fcc es una bcc 4
a
a

Cuando la esfera de Fermi toca a la PZB es porque

su radio es
3 4 3
qF aleación  
4 a a

 naleación
qF -aleación  (3  2 n)1/3 n aleación  3 nCu 
4
 1.36
a 3 a3 nCu
 SUPERFICIES DE FERMI en el modelo de NFE
Método de Harrison para materiales trivalentes
¿ Cómo puede pasarse de las superficies de Fermi para electrones libres a
las de electrones casi libres?

a) El potencial del cristal redondea las esquinas afiladas de las

superficies de Fermi.

(b) Las superficies de Fermi cortan perpendicularmente a las fronteras

de las zonas de Brillouin.

(c) La interacción del electrón con el potencial periódico del cristal

origina la aparición de bandas prohibidas en los límites de la zona.

d) El volumen total encerrado por la superficie de Fermi depende

solamente de la concentración electrónica y es independiente de los
detalles de la interacción con la red.
 SUPERFICIES DE FERMI DE ALGUNOS
METALES

Ver tabla de superficies de Fermi en archivo de internet

L = (1,1,1), X = (0,0,1), W = (1,1,0)
 CONDICIÓN DE LAUE

k´
g/2 g
O
(Origen de red
frontera de zona
recíproca) de Brillouin

Los electrones cuyas funciones de onda

k´  k  g  2 k . g
2 2 2

tienen su k cerca a la frontera de zona de

Brillouin satisfacen la condición de Laue y
experimentan reflexión de Bragg. Eso
2 k . g  g2
explica que los niveles de energía de los
electrones a uno y otro lado de una frontera
g 1
de zona de Brillouin experimenten un salto,
produciéndose un gap de energía. k.  g
g 2
 Densidad de estados

Densidad de estados para electrones libres en caja vacía 3D

3/2
V  2m 
g (E)   2  E 1/2
2 2  

Densidad de estados para electrones en sólido cristalino 3D

V dS
g (E ) 
4 3   E (k)
SE k
 Bandas Típicas

EF
eV


a

metales semiconductores

EF(meV)

Superconductores
FUNCIONES DE DISTRIBUCIÓN

• En el estado fundamental, es decir a 0K los

N electrones en un cristal ocupan los N
estados de energía más baja, pero ¿qué pasa
a una temperatura T ?
 Estadística Fermi Dirac
• La probabilidad de que un estado de energía E sea ocupado por un
electrón a una temperatura T es:
1
f (E) 
e (E-μ)/K BT  1
donde μ es el potencial químico (depende de cada material) y KB la constante
universal de Boltzmann.
Si la densidad de estados es g(E), entonces:

Número total de partículas N  f(E) g(E) dE

Energía total del sistema de partículas E total   E f(E) g(E) dE

Es válida para N partículas indistinguibles que obedecen principio de exclusión de

Pauli
 Estadística Fermi Dirac

f (E)
   EF EF   EF


a

metales semiconductores
 La Distribución Maxwell Boltzmann(*)
• Dado un sistema de N partículas como las moléculas de un gas en un recipiente,
el número de partículas en un nivel de energía Ei, con degeneración gi es:

N i  A g i e - E i /K BT
donde A es una constante de normalización y KB es la constante de Boltzmann
Si la densidad de estados es g(E), entonces:

Número total de partículas N  A  e -E/K BT g(E) dE

E total  A  E e -E/K BT g(E) dE

Energía total del sistema de partículas
Es válida para N partículas distinguibles que no obedecen principio de exclusión de
Pauli

(*) Ver: Solid -State and Semiconductor Physics por J.P. McKelvey
Modelo Semiclásico para la dinámica de los electrones en
un cristal

¿Cómo se comportan los electrones de un sólido

cristalino en presencia de campos externos?

1. Permanecen en su banda, n es constante

2. Un electrón con vector de onda k y energía En(k) tendrá

1
v  k E

2. Obedecen la segunda ley de Newton en la forma:
dk
F
dt
 Límite de validez del modelo semiclásico

Los electrones obedecen la dinámica semiclásica con la condición

de que los campos E y B no sean muy intensos. En la práctica esto
significa:
2
 E gap (k) 
e  a <<
EF
2
 E gap (k) 
c <<
EF

eB
c 
m
a es el parámetro de red
 Dinámica Semiclásica de los electrones
en un cristal
mv k

1 dE
v
 dk

v (-k )  - v (k )
v se puede interpretar como la velocidad de grupo de un paquete de
funciones de onda con k en el intervalo k  k
Dinámica Semiclásica de los electrones
en un cristal
mv k

1
v  k E


v (-k )  - v (k )

v también se puede interpretar

como el valor esperado del
operador i para una función de
m
Bloch
 Dinámica Semiclásica de los electrones
en un cristal
Segunda ley de Newton y momentum cristalino

dE(k )  
 e   v
dt

 
dE(k )  dk dE(k )  1
  k E(k ) 
dt dt  e    k E (k )
dt 


dk
  eε
dt
 Dinámica Semiclásica de los electrones
en un cristal
Segunda ley de Newton y momentum cristalino

dE(k )  
 e   v
dt

 
dE(k )  dk dE(k )  1
  k E(k )   e    k E (k )
dt dt dt 


dk 
  e
dt
 Dinámica Semiclásica de los electrones
en un cristal
Masa efectiva de un electrón en un
cristal
dv
a
dt
dv dk
a
dk dt
1 d 2E
a F
 dk 2

2

m  2
d E
dk 2
 Dinámica Semiclásica de los electrones en
un cristal
Masa efectiva de un electrón en un
cristal
dv
a
dt
dv dk
a
dk dt
1 d 2E
a F
 dk 2

2

m  2
d E
dk 2
 TENSOR MASA EFECTIVA

  1
v  v g  k E
  v x v x v x   dk 
 dv x    x
 dt   k x k y k z   dt 
    
 y     v y
dv v y v y   dk y 

dv x  v x dk x v x dk y v x dk z  dt    k k y k z   dt 
      x
  dk z 
dt  k x dt k y dt k z dt  dv z 

v z   
 v z v z
 dt   
 k x k y k z   dt 

  2E  2
E  2
E 
 2
 k X k x k y k x dk z 
 
  
 
2 2 2
E E E
m  12  
1

 k x k y k y 2
k z k y 
 
  E
2
 E  E 
2 2

 
 k x k z k z k y k z 
2
 CONCEPTO DE HUECO
Es una seudopartícula cuyas propiedades físicas son equivalentes a las de
todos los electrones de una banda con todos los estados ocupados
excepto uno.

El vector de onda del hueco es : kh   ke

La energía del hueco es: E h (k h )   E e (k e )

La velocidad del hueco es: v h (k h )  v e (k e )

m *
La masa efectiva del hueco es: h h(k )   m *
e (k e )

La carga eléctrica del hueco es: qh = e

 CONCEPTO DE HUECO

Ee Eh

kh
Eh (kh)
Ee(ke) Ke Kh
ke

kh   ke E h (k h )   E e (k e ) v h (k h )  v e (k e ) m*h (k h )   m*e  k e 

qh = e
 CONCEPTO DE HUECO

Ee Eh

ε kh

ke

dk h dk e
    eε  La carga eléctrica del hueco es positiva
dt dt
 Conductores, semicondutores y aislantes

BC
BC gap
gap

BV BV
BV

Conductores semicondutores aislantes

Conductores semicondutores y aislantes
Conductores semicondutores y aislantes

(a) aislantes, (b) metal o semimetal, (c) metal

Concentración de portadores en
semiconductores intrínsecos
 2k 2
E c (k )  E g 
E(electrones) 2m*e

 2k 2
E v (k ) 
2m*h

3/2
1  2m e 
g e (E)  2  2 
2π   
E  E  g
1/ 2

3/2
1  2m h 
g h (E)  2  2   E  1/ 2
2π   
E(huecos)
 Concentración de portadores en
semiconductores intrínsecos

E C2 EV 2
n  f (E) g(E) dE
E C1
p f h (E) g h (E) dE
EV 1
Concentración de portadores en
semiconductores intrínsecos
3/2 

 E  E 
1  2m e  E F /K BT
e  E/K BT dE
1/ 2
n  2 2  e g
2π    Eg


 1/2

0 e dx  2
1/ 2  x
x

3/2
 m e k BT  E F /K BT  E g /K BT
n  2 2 
e e
 2π 
 Concentración de portadores en
semiconductores intrínsecos

p   f h (E) g h (E) dE
0

1 1
f h  1  f (-E)  1 -   e  E F /K BT e -E/KBT
e  -E -E F  /K BT  1 e (E F  E)/KBT  1

3/2
1  2m h 
g h (E)  2  2 
 E  1/ 2
2π   

3/2
 m h k BT  -E /K T
p 2   e F B
 2π 
2
 3/2
 m e k BT  E F /K BT  E g /K BT
n  2 2 
e e
 2π 

3/2
 m h k BT  -E /K T
p 2   e F B
 2π 
2

3/2 3/2
 m e k BT  E /K T  E /K T  m h k BT  -E /K T
2  e F B e g B =2   e F B
 2π   2π 
2 2
 Concentración de portadores en
semiconductores intrínsecos
np
3/2 3/2
 m e k BT  E /K T  E /K T  m h k BT  -E /K T
2  e F B e g B =2   e F B
 2π   2π 
2 2

1 3  mh 
E F  E g  K BT Ln 
2 4  me 
1/T

3/2
 k BT 
ni  2  2 
 me mh  3/4
e
 E g /2K B T

2π 
 Estados de Impurezas

e2
V(r)  
4π ε r ε 0 r
1  m e   - m0e 4 
E d  2   
me ε r  m 0   2(4π ε 0  ) 2 
ε r  10 ,  0,2
m0
E d  0,01eV
 Estados de Impurezas

1  m h   - m 0e 4 
E a  2   
ε r  m 0   2(4π ε 0  ) 2 

Ordenes de magnitud:
E a  0,01 eV
Concentración de portadores en
semiconductores extrínsecos
 Evolución del nivel de Fermi con la temperatura

EF
B.C.
Ed

B.V. T
 Concentración de portadores en
semiconductores extrínsecos
3
 k BT 
np  4 2 
 me mh  e
3/2  E g /K B T

2π 

np  n 2
i
semiconductores extrínsecos tipo n

n i2
n  Nd  p
Nd

semiconductores extrínsecos tipo p

n i2
p  Na  n 
Na