Universidad Nacional Mayor de San Marcos

Universidad del Perú. Decana de América

Facultad de Ciencias Físicas
Escuela Profesional de Física

Estructura electrónica y ordenamiento de bandas en el

HgTe y CdTe

TRABAJO DE INVESTIGACIÓN
Para optar el Grado Académico de Bachiller en Física

AUTOR
Cristhian David HINOSTROZA VARGAS MACHUCA

ASESOR

Dr. Pablo Héctor RIVERA RIOFANO

Lima, Perú

2021
Reconocimiento - No Comercial - Compartir Igual - Sin restricciones adicionales

https://creativecommons.org/licenses/by-nc-sa/4.0/
Usted puede distribuir, remezclar, retocar, y crear a partir del documento original de modo no
comercial, siempre y cuando se dé crédito al autor del documento y se licencien las nuevas
creaciones bajo las mismas condiciones. No se permite aplicar términos legales o medidas
tecnológicas que restrinjan legalmente a otros a hacer cualquier cosa que permita esta licencia.
Referencia bibliográfica

Hinostroza, C. (2021). Estructura electrónica y ordenamiento de bandas en el HgTe

y CdTe. [Trabajo de investigación de bachiller, Universidad Nacional Mayor de San
Marcos, Facultad de Ciencias Físicas, Escuela Profesional de Física]. Repositorio
institucional Cybertesis UNMSM.
 HOJA DE METADATOS COMPLEMENTARIOS

Código ORCID del autor https://orcid.org/0000-0002-3290-8465

DNI o pasaporte del autor 73123370

Código ORCID del asesor https://orcid.org/0000-0003-3609-3889

DNI o pasaporte del asesor 08280449

Grupo de investigación Lab2d4

Perú.
Universidad Nacional Mayor de San
Marcos. Vicerrectorado de Investigación y
Agencia financiadora Posgrado. VRIP - UNMSM
Programa de Promoción de Trabajo de
Investigación para optar el grado
académico de Bachiller (Primera
Convocatoria).
Número de contrato: B20130360a
País: Perú
Departamento: Lima
Ubicación geográfica donde se
desarrolló la investigación Provincia: Lima

Latitud: -12.0453
Longitud: -77.0311
Año o rango de años en que
2020 - 2021
se realizó la investigación

Física de la Materia Condensada

Disciplinas OCDE
https://purl.org/pe-repo/ocde/ford#1.03.02
 Universidad Nacional Mayor de San Marcos
Universidad del Perú. Decana de América

FACULTAD DE CIENCIAS FÍSICAS

ESCUELA PROFESIONAL DE FÍSICA

ACTA DE SUSTENTACIÓN DE TRABAJO DE INVESTIGACIÓN PARA OPTAR EL

GRADO ACADÉMICO DE BACHILLER EN FÍSICA

Siendo las 18:00 horas del martes 22 de junio del 2021, en la Sala de Sesión Virtual de la
Facultad de Ciencias Físicas, bajo la Presidencia del Dr. Carlos Vladimir Landauro Sáenz
(Presidente), el Dr. Jesús Félix Sánchez Ortiz (Miembro), el Dr. Pablo Héctor Rivera Riofano
(Asesor) y la Lic. Rocío Amelia Montalvo Balarezo (suplente), se dio inicio a la Sesión
Pública de Sustentación de Trabajo de Investigación para optar el Grado Académico
de Bachiller en Física, del egresado:

Cristhian David Hinostroza Vargas Machuca

El Presidente del Jurado en primer lugar dio lectura al documento formal que designa al
Jurado y aprueba el acto de sustentación del Trabajo de Investigación, así como el
Resumen del Expediente e invitó al citado egresado a realizar la presentación y
exposición de su Trabajo de Investigación:

“Estructura electrónica y ordenamiento de bandas en el HgTe

y CdTe”
Concluida la exposición del candidato y luego de las preguntas de rigor por parte del
Jurado, el Presidente del Jurado con el apoyo de la Unidad de Informática, invitó al
egresado y al público a abandonar momentáneamente la Sala de Sesión Virtual, para
dar paso a la deliberación y calificación por parte del Jurado.

Al término de la deliberación del Jurado el Dr. Carlos Vladimir Landauro Sáenz invitó al
candidato y al público en general a pasar a la Sala de Sesión Virtual, para dar lectura
al Acta de Calificación, en cuyo contenido se ha registrado la nota obtenida por el
egresado, la misma que ha sido de:

19 diecinueve Aprobado con máximos honores

............................. ....................................... .............................................................
(NÚMERO) (LETRAS) (MENCIÓN)

Habiendo concluido la Sustentación del Trabajo de Investigación, el Presidente del

Jurado, remitirá el Acta al Vicedecano Académico de la Facultad de Ciencias Físicas,
quien la suscribirá dando fe de la realización del acto virtual de manera remota.

Dirección: Calle Germán Amézaga N.º 375 - Lima Escuela Profesional de Física
Facultad de Ciencias Físicas 6197000 Anexo 3805
Apartado Postal: 14-0149 – Lima 14 Perú E-mail [email protected]
 Universidad Nacional Mayor de San Marcos
Universidad del Perú. Decana de América

FACULTAD DE CIENCIAS FÍSICAS

ESCUELA PROFESIONAL DE FÍSICA

El Acta de Calificación será remitida al Director de la Escuela Profesional de Física quien

a su vez solicitará la aprobación del Consejo de Facultad para el otorgamiento del
Grado Académico de Bachiller en Física al egresado Cristhian David Hinostroza Vargas
Machuca.

Siendo las 19:25 horas, se dio por concluido el acto académico, expidiéndose tres (03)
Actas originales de la Sustentación del Trabajo de Investigación.

Ciudad Universitaria, martes 22 de junio de 2021.

Dr. Carlos Vladimir Landauro Sáenz Dr. Jesús Félix Sánchez Ortiz
PRESIDENTE MIEMBRO

Dr. Pablo Héctor Rivera Riofano Lic. Rocío Amelia Montalvo Balarezo
ASESOR SUPLENTE
Firmado digitalmente por ROJAS
ACUÑA Joel FAU 20148092282 soft
Motivo: Soy el autor del documento
Fecha: 30.06.2021 15:51:58 -05:00

Dr. Joel Rojas Acuña

VICEDECANO ACADÉMICO - FCF

Datos de la plataforma virtual institucional del acto de sustentación:

Enlace zoom:

https://zoom.us/j/91488222047?pwd=a2hRd3JEeUE1bFZPem5BbTFLQXZMdz09
ID de reunión: 914 8822 2047
Código de acceso: 909818

Grabación archivada en:

Dirección: Calle Germán Amézaga N.º 375 - Lima Escuela Profesional de Física
Facultad de Ciencias Físicas 6197000 Anexo 3805
Apartado Postal: 14-0149 – Lima 14 Perú E-mail [email protected]
 i

A mi abuela
Amalia y mi tío
Juan.
ii

Agradecimientos

Un profundo agradecimiento a mis padres y hermanos, quienes me

ayudaron a lo largo de mi formación, aún sin entender lo que hacía.
Al doctor Pablo Rivera que, como asesor, no solo me guió en el
desarrollo del presente trabajo de investigación, sino que amplió mi
visión de la realidad.
A la facultad de Ciencias Físicas, administrativos, docentes y alum-
nos, por formar una comunidad capaz de generar retos.
 iii

Don’t become a
scientist!
Jonathan I. Katz
iv

Resumen
En este trabajo de investigación, se presentan las nociones básicas para
entender la formación del pozo de potencial HgCdTe/HgTe/HgCdTe
estudiando la inversión de bandas del HgTe originada por el intenso
acoplamiento espín órbita del mercurio.
Se calcularon numéricamente las estructuras de bandas del HgTe y del
CdTe considerando el acoplamiento espín-órbita, y la interacción entre
bandas mediante la teoría de funcionales de densidad. Se observó que
el HgTe presenta una inversión de bandas en el punto de simetría Γ,
con la banda de conducción siendo conformada por estados orbitales
tipo p, mientras que el CdTe posee un ordenamiento normal, con la
banda de conducción conformada por estados orbitales tipo s; por lo
que ambos semiconductores pueden formar pozos cuánticos de tipo
III, cumpliendo, además, las condiciones básicas para la existencia del
estado de aislante topológico.
Palabras claves: Semiconductores, DFT, estructura de bandas, aco-
plamiento espín-órbita.
 v

Abstract
In this research work is explained the basic notions to understand how
HgCdTe/HTe/HgCdTe quantum well is obtained studying the band
invertion of HgTe, originated by the intense spin-orbit coupling from
mercury.
Band structures for bulk HgTe and CdTe were calculated numerically
considering spin-orbit coupling and band interaction using density fun-
ctional theory. It was observed that HgTe shows a band inversion in
the high symmetry point Γ, where the contribution to the conduction
band is asociated to the p-state orbital; while CdTe posses a regular
band ordering, having a contribution of the s-type orbital in the con-
duction band; thus it is possible to create a type III quantum well
using both semiconductors, moreover, satisfying the basic conditions
for a topological insulator state to exist.
Keywords: Semiconductors, DFT, band structure, spin-orbit
coupling.
Índice general

1. Introducción 1

2. Fundamento teórico 5
2.1. Estructuras cristalinas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
2.2. Espacio y red recíprocos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.3. Estructura de bandas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.4. Heteroestructuras . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.5. Pozos cuánticos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.6. Fermiones de Dirac, Weyl y Majorana . . . . . . . . . . . . . . . . . . . . . . . 13
2.6.1. Fermiones de Dirac . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.6.2. Fermiones de Weyl . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.6.3. Fermiones de Majorana . . . . . . . . . . . . . . . . . . . . . . . . . . . 15

3. Herramientas teórico-computacionales 17
3.1. Teoría de funcionales de densidad . . . . . . . . . . . . . . . . . . . . . . . . . 17
3.1.1. Aproximación local de densidad de espín (LSDA) . . . . . . . . . . . . 19
3.1.2. Expresión del estado y elección de base . . . . . . . . . . . . . . . . . . 20

4. Resultados y discusión 23

5. Conclusiones 33

Bibliografía 35

vii
Capítulo 1

Introducción

Durante los últimos años, los estados cuánticos de la materia fueron tomando mayor rele-
vancia debido a un notable avance en las técnicas de síntesis de heteroestructuras semiconduc-
toras, como el crecimiento de capas usando la técnica Molecular Beam Epitaxy, y los nuevos
y sofisticados métodos de espectrometría, como el Angle Resolved Photoemision Spectroscopy,
haciendo posible su realización experimental y estudio a detalle.

Entre estos estados cuánticos de la materia, los materiales topológicos poseen un prota-
gonismo especial por sus características únicas que, como su nombre lo indica, son descritas
usando la topología, una rama de la matemática que estudia las propiedades geométricas de
un sistema que se mantienen invariantes ante cambios infinitesimales, siendo la fase de aislante
topológico una de las más estudiadas en la actualidad.

La historia de los aislantes topológicos tiene su origen en el descubrimiento del efecto Hall
cuántico. El efecto Hall cuántico es una versión cuantizada del efecto Hall clásico. En este
efecto, y debido a la naturaleza cuántica de los portadores de carga, el número de los canales
conductores irá cambiando tomando valores enteros positivos, siendo este el motivo por el
que la conductancia Hall cuántica varía discretamente al ir elevando la magnitud del campo
eléctrico aplicado, formando estados de conductancia constante, o plateaus a muy bajas tem-
peraturas.

Este estado de la materia fue descubierto teóricamente en 1975, cuando Tsuneya Ando,
Yukio Matsumoto y Yasutada Uemura, de la Universidad de Tokio, predijeron la cuantización

1
2 Capítulo 1. Introducción

Figura 1.1: Resistencia Hall de un gas de electrones a bajas temperaturas en función del campo magnético
aplicado [2].

de la conductancia Hall a escala microscópica en un material bidimensional [3]. Sin embar-

go, el descubrimiento fue tan sopresivo que no lo consideraron como válido. Fue así que,
recién en 1980, Klaus von Klitzing descubre el efecto Hall cuántico mientras exploraba las
características y propiedades de un Transistor de Efecto de Campo de Metal-Óxido (MOS-
FET) producido a base de silicio, fabricado por Michael Pepper y Gerhard Dorda [4]. Gracias
a este descubrimiento, von Klitzing fue galardonado con el premio Nobel de Física el año 1985.

La primera formulación teórica para explicar la formación de los canales conductores fue
el trabajo de Laughlin [5], en el que usó un argumento sobre la invariancia Gauge del vec-
tor potencial del campo electromagnético. Inspirándose en este trabajo, Thoules, Qian Niu y
Yong-Shi Wu, describieron la conductancia Hall como un invariante topológico [6], al gene-
ralizar la invariancia Gauge a la evolución de la fase geométrica de los parámetros del sistema.

Usando este trabajo como punto de partida, se inició la búsqueda de materiales que pudie-
sen presentar el mismo comportamiento que el efecto Hall cuántico sin la necesidad de aplicar
un campo magnético externo, llegando de esta manera al efecto Hall cuántico de espín, en
el cual el campo magnético es reemplazado por un acoplamiento espín-órbita intenso, siendo
este estado, uno tipo aislante topológico. Fue así que Andrei Bernevig, Taylor L. Hughes y
Shou-Cheng Zhang, el año 2006, propusieron la existencia del estado de aislante topológico
 3

en pozos cuánticos del semiconductor CdTe/HgTe/CdTe [7]. Propuesta que se logró demos-
trar experimentalmente un año después, el 2007, por el equipo de Markus König [8], de la
universidad de Würztburg, quienes realizaron el experimento en un pozo cuántico de HgTe
con barreras de Hg0.3 Cd0.7 Te.

La existencia del aislante topológico bidimensional en el pozo cuántico de CdTe/HgTe/CdTe

se debe, principalmente, a la inversión de bandas existente en el HgTe y al intenso acopla-
miento espín-órbita en las aleaciones HgCdTe. Sin embargo, en la literatura en español no
existen muchas referencias que hayan estudiado estas aleaciones, o incluso, es cómun encontrar
aquellos que no consideran la interacción espín-órbita en los cálculos numéricos al emplear la
teoría de funcionales de densidad, DFT por sus siglas en inglés, siendo esta característica una
de las más importantes a considerar en estos materiales.

Por este motivo y teniendo en cuenta los antecedentes presentados, el presente traba-
jo se centrará en la obtención de la estructura de bandas electrónicas de la aleación triple
Hg1-x Cdx Te en los casos límites x = 0 y x = 1. Es decir, para los semiconductores HgTe y
CdTe. Para lograr este objetivo se hizo uso de la teoría de funcionales de densidad. Una vez
expuesta la relación de dispersión se procederá a analizar los gaps de energía, en función de
su localización en la zona de Brillouin, y su ancho, además de la obtención del aporte de los
orbitales atómicos para cada elemento cerca al nivel de la energía de Fermi.

Este trabajo está estructurado de la siguiente manera: en el capítulo 2 se hace una revi-
sión teórica de la estructura cristalina de los sólidos y la teoría de bandas. En el capítulo 3
se expone la metodología a usar, específicamente, la Teoría de Funcionales de Densidad. En
el capítulo 4 se exponen y discuten los resultados obtenidos. Finalmente, en el capitulo 5, se
dan a conocer las conclusiones de este trabajo.
Capítulo 2

Fundamento teórico

2.1. Estructuras cristalinas

En la naturaleza, los elementos químicos se unen formando diversas estructuras. Entre
estas estructuras, los cristales poseen una forma ordenada, siendo definidos por la Unión
Cristalográfica Internacional como "todo material que posee un patrón de difracción discreto"
[9].
Las estructuras periódicas que conforman los cristales poseen simetrías de traslación y
rotación discretas que influyen en el comportamiento del electrón en su interior [10]. Estas
estructuras se pueden representar como una red de puntos, también llamada red directa R:

R = n1 a 1 + n2 a 2 + n3 a 3 , (2.1)

donde n1 , n2 y n3 son números enteros, y los vectores a1 , a2 y a3 definen una celda en el

espacio de posiciones.
A las redes que pueden generarse en el espacio tridimensional, únicamente mediante la
traslación de una celda, se las conoce como las redes de Bravais [11]. Estas redes se dividen
en 7 sistemas: triclínico, monoclínico, ortorómbico, tetragonal, trigonal, hexagonal y cúbico.
En este último sistema existen 3 redes de Bravais: la cúbica simple (SC), de cara centrada
(FCC) y de cuerpo centrado (BCC), por sus siglas en inglés.
Las redes de Bravais son una abstracción matemática del ordenamiento geométrico exis-
tente en el sistema. Por lo que, para formar el cristal, es necesario que cada celda albergue
un número de átomos en posiciones definidas, lo que recibe el nombre de base atómica. De

5
6 Capítulo 2. Fundamento teórico

Figura 2.1: Redes de Bravais del sistema cúbico [11].

esta forma, un cristal se puede describir conociendo la celda que lo conforma, los elementos
existentes en la celda, y la posición de estos elementos dentro de la celda.
Teniendo en cuenta lo anterior, no existe una restricción para el tamaño máximo de la
celda. Sin embargo, es preferible trabajar con las celdas que poseen un volumen mínimo,
llamadas celdas primitivas [10], a partir de la cual se construye una celda que posee las
simetrías puntuales del cristal, llamada celda de Wigner-Seitz. Esta útima se obtiene a partir
de las mediatrices de los vectores que definen la celda primitiva.

2.2. Espacio y red recíprocos

Muchos de los problemas en física del estado sólido son difíciles, o hasta imposibles de
resolver, cuando solo se trabajan en el espacio de posiciones. Por esta razón, se traslada el
problema al espacio cuasimomentos k = p/~, llamado también espacio recíproco.
Como se vió anteriormente, las redes de Bravais forman una red de puntos en todo el
espacio de posiciones, como lo muestra la ecuación (2.1). Por este motivo, el cristal posee
una periodicidad que depende del módulo y dirección de los vectores que conforman la celda
primitiva. Esta periodicidad en el sistema permite el uso del teorema de Bloch. Este teorema
establece que la función de onda del electrón toma la forma de una onda plana multiplicada
por una función periódica u(r), con periodicidad R, en el espacio de posiciones [10]:

ψk (r) = eik·r uk (r). (2.2)

Por otro lado, se conoce que un estado cuántico puede trasladarse en el espacio al aplicar
el operador traslación T , el cual traslada a la función de onda en una dirección ~r:
2.2. Espacio y red recíprocos 7

T (r)φ(x) = φ(x + r), (2.3)

donde r puede tomar cualquier dirección y magnitud.
De esta manera, aplicando el operador de traslación sobre la función de Bloch, de modo
que la traslade una distancia R, que coincide con la periodicidad de la red directa, se tiene
que:

T (R)ψk (r) = eik·(r+R) uk (r + R). (2.4)

La periodicidad del sistema cristalino permite considerar las condiciones de frontera de
Born-Von Karman, en donde la función de onda es periódica en cada celda, o dicho de otro
modo, la función de onda es una función periódica con periodicidad R. De esta manera, la
ecuación (2.4) debe cumlplir con la condición:

eik·R = 1. (2.5)
Teniendo en cuenta que la ecuación (2.5) restringe lo valores de k, definimos como G a los
punto k, del espacio recíproco, que cumplan con la condición mostrada en la ecuación (2.5).
Así, se define matemáticamente a la red de puntos del espacio recíproco como los puntos que
cumplen [12]

eiR·G = 1. (2.6)
Suponiendo que el vector G puede expresarse como una combinación lineal de vectores
b1 , b2 y b3 , y considerando que se debe cumplir la relación

ai · bj = 2πδij , (2.7)

se obtiene una expresión para la base de vectores en el espacio recíproco:

a2 × a3 a3 × a1 a1 × a2
b1 = 2π , b2 = 2π , b3 , = 2π .
a1 · (a2 × a3 ) a1 · (a2 × a3 ) a1 · (a2 × a3 )
De la misma forma que la combinación lineal de los vectores a1 , a2 y a3 forma una red de
puntos R en el espacio de posiciones, la combinación lineal de los vectores b1 , b2 y b3 definen
una red de puntos G en el espacio recíproco denominada red recíproca.

G = m1 b 1 + m2 b 2 + m3 b 3 (2.8)
8 Capítulo 2. Fundamento teórico

Esta red define una celda equivalente a la de Wigner-Seitz, llamada zona de Brillouin, en el
que se encuentran los puntos de alta simetría de la red recíproca.

2.3. Estructura de bandas

Cuando se estudia al sólido es necesario considerar tanto el comportamiento individual de
los electrones e iones que conforman el cristal, como las interacciones existentes entre ellas
y entre sí mismas. Teniendo esto en consideración, el hamiltoniano H que describe al sólido
tiene la forma [13]

Hsol = He + Hion + He−ion , (2.9)

donde He , Hion y He−ion son los hamiltonianos que describen a los electrones, iones, e inter-
acciones electrón-ión, respectivamente.
Para poder resolver el hamiltoniano en la ecuación (2.9) es necesario usar la aproximación
de Born-Oppenheimer, la cual considera que la velocidad del electrón es mucho mayor que la
de los iones, y que la masa de los iones es mucho mayor que la masa de los electrones. De
esta forma, es posible desacoplar la ecuación (2.9) en dos ecuaciones semi-independientes; una
correspondiente solo a los iones, y otra solo a los electrones. Gracias a esta aproximación, la
ecuación que expresa el comportamiento de los electrones dentro del sólido es:

~2 2
− ∇ + Ve−e + Ve−ion = He , (2.10)
2me
en donde el primer término corresponde a la energía cinética de los electrones; Ve−e , a la in-
teracción electrón-electrón; y Ve−ion , a la interacción electrón-ión, la cual considera la posición
de los iones como parámetros.
Usando la función de onda en (2.2) es posible expresar los autovalores de energía como
función de los vectores de onda En (k), llamada relación de dispersión. El valor de k puede
variar en todo el espacio recíproco, sin embargo, debido a la simetría de traslación de la red
recíproca, el estudio puede reducirse a la primera zona de Brillouin, en especial a los puntos
y líneas de simetría que posee.
Para representar la relación de dispersión de una forma que pueda analizarse adecuada-
mente se grafican los valores de En (k) a lo largo de un camino que une los puntos de alta
simetría de la zona de Brillouin, obteniendo así una representación bidimensional en la que
2.4. Heteroestructuras 9

los autovalores de energía En (k) forman bandas de energía que, en su conjunto, forman la
llamada estructura electrónica, la cual brinda información de las propiedades del material.
Una de las propiedades más sencillas de observar es el comportamiento electrónico que
posee el sistema estudiado, que puede ser: metal, aislante o semiconductor. Esto dependerá
del comportamiento de la bandas de energía cerca del nivel de Fermi y del ancho que posee
el espacio entre el máximo de la última banda de valencia y el mínimo de la primera banda
de conducción, llamado gap global de energía.

2.4. Heteroestructuras
Las propiedades de un sistema dependen, entre muchos factores, de sus dimensiones y
composición química [14], las cuales afectan directamente las propiedades del material.
Una heteroestructura es aquel sistema que se obtiene al unir dos tipos de semiconductores
[12] al hacer crecer un material sobre otro a lo largo de un eje de crecimiento. Así, el cambio
abrupto en la composición del material, como el tamaño de la celda del cristal, generará una
discontinuidad en las bandas de energía [15]. Dicha discontinuidad puede clasificarse tomando
en cuenta el gap global de cada semiconductor, dividiendo las heteroestructuras en 3 tipos [16],
como muestra la figura (2.2).

Figura 2.2: Tipos de heteroestructuras de acuerdo a sus discontinuidades. Adaptado de [16].

Para que las heteroestructuras posean una buena calidad al unir dos semiconductores A y
B, es necesario que los parámetros de red de los materiales a unir sean similares. El parámetro
que caracteriza esta condición es la incompatibilidad de redes η, y se define como:
10 Capítulo 2. Fundamento teórico

|aA − aB |
η=2 , (2.11)
(aA + aB )
donde aA y aB son los parámetros de red de los cristales A y B. Si η es muy grande la
heteroestructura presentará dislocaciones en la unión de ambos cristales. Para evitar estas
dislocaciones se suele colocar entre los cristales A y B aleaciones conformadas con las mismas
especies atómicas que A y B, para que el cambio entre los tamaños de las celdas sea lo más
suave posible.

2.5. Pozos cuánticos

Como se vio anteriormente, el comportamiento de un electrón depende de la composición
del material en el que se encuentra, sin embargo, también depende de las dimensiones geo-
métricas en las que está confinado. Normalmente, estos portadores de carga se encuentran en
estructuras tridimensionales tipo bulk, pero al reducir el tamaño de una de las dimensiones
del sistema, los electrones manifiestan sus propiedades cuánticas con mayor intensidad, lo
cual altera las propiedades del sistema de manera significativa. En la figura (2.3) se observa
la densidad de estados de un electrón casi libre en una red cuadrada al ir reduciendo las
dimensiones.

3D
2D
1D
DOS

Energía

Figura 2.3: Densidad de estados para el electrón casi libre en una red cuadrada al reducir dimensiones.

Una forma de obtener estructuras bidimensonales es uniendo dos heteroestructuras con-

formadas por los mismos componentes a lo largo de un eje de crecimiento, de tal forma que el
2.5. Pozos cuánticos 11

semiconductor con gap de energía menor se encuentre entre dos capas del semiconductor con
gap de energía mayor. Se formará entonces un pozo de potencial [17] que seguirá la forma fun-
cional de la composición del semiconductor [15], pudiendo ser pozos cuadrados, parabólicos,
triangulares, etc.
Si el ancho del semiconductor que se encuentra entre las barreras es demasiado grande el
sistema se comportará como un sistema macroscópico. Si, por el contrario, se va reduciendo
el ancho hasta un tamaño igual o menor a la longitud de coherencia de fase, los electrones
manifestarán sus propiedades cuánticas [18]. Esto se logra al construir un pozo de potencial
del tamaño de algunas celdas unitarias del semiconductor.
En este régimen, al interior del pozo de potencial los portadores de carga tendrán una
función de onda que se reflejará en las paredes del pozo a lo largo del eje de crecimiento,
de modo que existen "modos normales", en su equivalente clásico, que corresponden a las
funciones de onda de Bloch. De esta manera, las partículas dentro del pozo de potencial
estarán bajo un confinamiento cuántico cuasi-bidimensional.
Las propiedades de un pozo cuántico dependerán de la altura del pozo, o de la composición
química de los materiales que lo conforman; y de su ancho, que es el grosor del semiconductor
que se encuentra en medio. Por otro lado, es posible clasificar los pozos cuánticos siguiendo
la clasificación de las heteroestructuras que fue mostrada en la figura (2.2).

Figura 2.4: Pozo cuántico formado por heteroestructuras tipo I.

La manera más sencilla de construir un pozo cuántico es uniendo heteroestructuras de

tipo I como se muestra en la figura (2.4), de modo que no exista un solapamiento entre el
confinamiento de electrones y huecos.
De forma ilustrativa, para una partícula libre con masa efectiva m∗ , dentro de un pozo de
potencial cuadrado infinito en la dirección del eje de crecimiento ẑ, se tendría [19]:
12 Capítulo 2. Fundamento teórico

kzn = 2π/λn = nπ/L, (2.12)

donde kzn son los vectors de onda de Bloch en la dirección ẑ, λn son las longitudes de onda
permitidas, y n son números enteros positivos mayores que cero. Mientras que la energía del
sistema en el eje ẑ, teniendo en cuenta (2.12), sería:

~2 kzn
2
~2 πn 2
 
Enz = = , (2.13)
2m∗ 2m∗ L
por lo que, al no existir restricciones en el plano perpendicular a ẑ, la energía total del sistema
es:

~2 2 ~2 πn 2
 
2
En = (k + k ) + . (2.14)
2m∗ x y
2m∗ L
Los estados energéticos que puede tomar la partícula son infinitos, lo que guarda relación
con el tamaño del pozo. Sin embargo, en el pozo de la figura (2.4), y en todos los pozos de
potencial no ideales, los autovalores de energía la partícula son finitos, y dependen del alto y
ancho del pozo, así como de la penetración de la autofunción asociada a la energía propia en
las barreras de potencial.

Figura 2.5: Pozo cuántico formado por heteroestructuras tipo III.

Un caso diferente es el pozo cuántico conformado por heteroestructuras tipo III. En este
tipo de pozo cuántico se tienen confinados a los electrones y huecos interactuando entre sí,
por lo que se necesita considerar ambos portadores de carga para describir el comportamiento
de las bandas de energía. Para una descripción adecuada del sistema, el hamiltoniano que
describe este tipo de pozo cuántico necesita ser construido tomando en cuenta la base de
orbitales atómicos que aportan a la estructura de bandas, cerca al gap de energía, de los
semiconductores que conforman el pozo.
2.6. Fermiones de Dirac, Weyl y Majorana 13

2.6. Fermiones de Dirac, Weyl y Majorana

La ecuación de Dirac expresa el comportamiento de las partículas cuánticas en el régimen

relativista. Surge de asumir la naturaleza cuadrática de la energía predicha en la relación de
Einstein, de modo que, para que pueda cumplirse esta relación, Dirac estableció que la función
de onda ψ(x) debe cumplir la ecuación [20]:

(iγ µ ∂µ − m) ψ(x) = 0, (2.15)

donde m es la masa de la partícula cuántica y γ µ cumplen las condiciones del álgebra de

Clifford:

{γ µ , γ ν } = 2g µν . (2.16)

Así, la elección de γ µ no es única, por lo que cada elección de valores corresponde a una
representación o base diferente que, a su vez, influye en las soluciones que se puedan obtener.
Dentro de las propiedades que deben cumplir γ µ es que no conmuten entre sí, lo que
requiere que sean matrices de mínimo 2 × 2. Por esta razón, las soluciones que se obtienen en
las representaciones usadas habitualmente son funciones de onda con cuatro componentes.
Dentro de las representaciones más usadas se encuentran las representaciones de Dirac,
Weyl y de Majorana [21].

2.6.1. Fermiones de Dirac

Los fermiones de Dirac son aquellas partículas que se presentan un comportamiento guiado
por la función de onda obtenida usando la representación de Dirac. Esa representación es una
de las más usadas al resolver la ecuación de Dirac y da como soluciones las funciones de onda
de una partícula y su antipartícula.
Las matrices γ µ que conforman la base de Dirac son:
14 Capítulo 2. Fundamento teórico

   
1 0 0 0 0 0 0 1
   
0 1 0 0 0 0 1 0
   
γ0 = γ1 =

   
0 0 −1 0 0 −1 0 0
   
 
   
0 0 0 −1 −1 0 0 0

   
0 0 0 −i 0 0 1 0
   
0 0 i 0  0 0 0 −1
   
γ2 = γ3 =

   
0 i 0 0 −1 0 0 0
   

   
−i 0 0 0 0 1 0 0

La relación de dispersión para estas partículas tiene la forma

q
E= p 2 c 2 + m 2 c4 , (2.17)

la cual se reduce a una ecuación lineal cuando la masa m = 0. Esto forma una relación de
dispersión lineal, o de forma cónica en tres dimensiones, llamado cono de Dirac:

E = ±|p|c. (2.18)

2.6.2. Fermiones de Weyl

Los fermiones de Weyl son soluciones de la ecuación de Dirac para partículas sin masa. A
su vez, la representación que usa para las matrices en la ecuación de Dirac es la representación
de Weyl, o también llamada representación de quiralidad.
Los biespinores que conforman la solución de la ecuación de Dirac son llamados también
espinores de Weyl.
Las matrices γ µ que conforman la base de Weyl son:
16 Capítulo 2. Fundamento teórico

ν̃(x) = ν̃ ∗ (x). (2.19)

Esta condición significa que se busca una solución para la función de onda, en la ecuación
de Dirac, que sea igual a su compleja conjugada. Físicamente, la función de onda describe
una partícula que es igual a su antipartícula.
Debido a esta condición, los dos biespinores que conforman la solución no serán indepen-
dientes, sino que estarán relacionadas entre sí.
Las matrices γ µ que conforman la base de Majorana son:

   
0 0 0 −i 0 i 0 0
   
0 0 i 0 i 0 0 0
   
γ0 =   γ1 =  
0 −i 0 0 0 0 0 i
   

   
i 0 0 0 0 0 i 0

   
0 0 0 −i i 0 0 0
   
0 0 i 0 0 −i 0 0
   
γ2 = γ3 =

.
   
0 i 0 0 0 0 i 0
  

   
−i 0 0 0 0 0 0 −i
Capítulo 3

Herramientas teórico-computacionales

3.1. Teoría de funcionales de densidad

La ecuación que describe a los electrones en un sólido, dada por (2.10), es imposible
de resolver analíticamente para más de 3 electrones, en especial debido a las interacciones
existentes entre ellos y el ión, y por la gran cantidad de variables que se necesitarían al
intentar resolver el problema para cada electrón de forma individual. Viendo este problema,
Hohenberg y Kohn [23] proponen dos teoremas que buscan deshacerse de estos problemas:

1. Los observables de un sistema son funcionales de la densidad electrónica del estado

fundamental.

2. Esta densidad electrónica es única para cada sistema.

Estas suposiciones simplifican considerablemente el problema de muchos cuerpos, ya que se

pasa de estudiar la función de onda dependiente de las posiciones de cada uno de los electrones
en el sistema, a requerir únicamente la densidad electrónica dependiente de la posición en el
cristal.
De esta manera es posible expresar el hamiltoniano (2.10) en función de la densidad
electrónica. Además de esto, se necesita reducir el hamiltoniano de muchas partículas, a uno
que exprese el comportamiento de partículas no interactuantes [24]. Primero expresamos la
energía como una funcional de la densidad electrónica n(r) = |ψ(r)|2 :

F [n(r)] = hΨ| H |Ψi . (3.1)

17
18 Capítulo 3. Herramientas teórico-computacionales

Así, la energía del sistema se convierte en un funcional de la densidad electrónica. Usando

el principio variacional, la energía que posee el sistema en su estado fundamental será aquella
que minimize la energía, además que esta densidad electrónica será única para el sistema [25].
Entonces, expresamos la energía en función de un conjunto de orbitales de partículas
independientes φi , llamados orbitales de Kohn-Sham. De esta manera se pude expresar el
funcional de la energía F [n(r)], como:

~2 e2 Z Z n(r)n(r ′ )
F [n(r)] = ∇r |φi i + dn(r)dn(r ′ ) + E XC [n(r)], (3.2)
X
hφi | −
i 2m 2 |n(r) − n(r ′ )|

donde la sumatoria se realiza sobre cada orbital i, mientras que el primer término representa
la energía cinética en la base de orbitales de Kohn-Sham; el segundo término, la interacción
coulumbiana semiclásica de electrón-electrón, expresada con la densidad electrónica n(r) =
|φi (r)|2 ; y el tercer término, la energía de intercambio y correlación, expresado también
P
i

como un funcional de la densidad electrónica, en donde están incluídos todos los demás efectos
cuánticos generados por las interacciones de muchos cuerpos en el sistema real.
Esta energía de intercambio y correlación es uno de los términos más importantes en la
teoría de funcionales de densidad. El intercambio se refiere al cambio del signo de la función
de onda, luego del intercambio de posición de los electrones en el sistema, debido al principio
de exclusión de Pauli; mientras que la correlación refiere al efecto del movimiento de los
electrones sobre las interacciones entre estas partículas.
Usando los mutiplicadores de Lagrange para minimizar la energía en función de la den-
sidad, se obtiene una ecuación similar a la ecuación de Schrödinger para partículas indepen-
dientes, llamada ecuación de Kohn-Sham:

~2 2
" #
− ∇ + Vef f (r) φi (r) = ǫi φi (r), (3.3)
2m r
donde φ(r) son las funciones de onda de Kohn-Sham, y Vef f (r) es el potencial efectivo. Es de
suma importancia notar que la ecuación (3.3) expresa el comportamiento de una cuasipartícula
independiente que imita el comportamiento de un electrón no independiente presente en un
cristal. Por lo tanto, la ecuación de Kohn-Sham no es una expresión para obtener la energía
de los electrones en el cristal de forma directa sino que lo hace a través de la simplificación
del problema de muchos cuerpos al de un problema de partículas independientes. Por otro
lado, esta no especifica la forma funcional que deben poseer los orbitales de Kohn-Sham.
3.1. Teoría de funcionales de densidad 19

Se tiene, entonces, que el potencial efectivo que actúa sobre la cuasipartícula independiente
está expresado en función de la densidad electrónica y la posición:

Z
n(r ′ ) δE XC [n(r)]
Vef f (r, n(r)) = Vext (r) + e2 dr ′
+ . (3.4)
|r − r ′ | δn(r)

El primer término de la ecuación (3.4) expresa el potencial iónico de la red, por lo que la
geometría del sistema se encuentra implícita en ese término.
En la ecuación (3.4) el término de intercambio y correlación E XC no es conocido. Además,
y debido a que el potencial efectivo depende de la densidad electrónica, el valor de la energía
en el estado fundamental se obtiene por autoconsistencia hasta encontrar una densidad que
lo minimize. De esta manera, es necesario escoger un término de intercambio y correlación
que logre acercarse al comportamiento real del sistema a estudiar.

3.1.1. Aproximación local de densidad de espín (LSDA)

Una de las energías de intercambio y correlación más usadas es la LSDA, como la formu-
laron Perdew y Wang [26]. Este modelo considera un gas de electrones homogéneo en el que
los efectos de intercambio y correlación poseen un comportamiento local.
Las energías de intercambio y correlación dependerán del parámetro de densidad del gas
de electrones rs , considerado como el radio del volumen que ocupa un electrón en el gas, y
de la polarización relativa de espín ζ [24]:

4π ↑ n↑ −, n↓
rs = (n + n↓ )−1/3 ζ= ↑ , (3.5)
3 n + n↓

donde n↑ y n↓ son las densidades de espín up y down, resectivamente.

Para esta aproximación, la energía de intercabio es conocida analíticamente, mientras que
la energía de correlación solo se conoce analíticamente en los casos límites rs → 0 y rs → ∞.
De esta manera, la energía de intercambio y correlación será la integral sobre el espacio
de la densidad de energía de intercambio y correlación suponiendo que es la misma que en el
gas de electrones homogéneos ǫhom :

Z
EXC [n , n ] =
↑ ↓
d3 rn(r)ǫhom
XC (n , n ).
↑ ↓
20 Capítulo 3. Herramientas teórico-computacionales

3.1.2. Expresión del estado y elección de base

Para poder resolverse, los orbitales de Kohn-Sham necesitan ser expresados en una base
que considere las características del hamiltoniano. En otras palabras, si el electrón está o no
localizado.
Una de las bases más usadas es la de ondas planas, considerando solo un orbital de Kohn-
Sham se tiene:

φ(r) = = Cq eiq·r . (3.6)

X

Por otro lado, el término del potencial efectivo en la ecuación de Kohn-Sham (3.3) puede
expresarse mediante su transformada de Fourier debido a la periodicidad que posee por la red
cristalina, teniendo la forma:

Vef f (r) = Vk′ Cq−k′ eiq·r . (3.7)

X

qk′

Reemplazando (3.6) y (3.7) en (3.3), y realizando el cambio de variable q = k − G, donde

G es el vector recíproco que cumple la misma función que R en el espacio de posiciones, se
obtiene:

1
 
|k − G|2 − ǫ Ck−G + VG′ −G Ck−G′ = 0. (3.8)
X
2 G ′

Esta ecuación se expande en forma de ecuación matricial, cuyo tamaño depende del tamaño
máximo del vector recíproco G que se elija, teniendo una expresión de la siguiente forma:
  
0
 ǫk−G1 VG2 −G1 VG3 −G1 · · · Ck−G1 
VG −G ǫ0k−G2 VG3 −G2 · · ·
 Ck−G2 
  
 
 1
= 0. (3.9)
2

ǫ0k−G3
  
V VG2 −G3 · · ·
 G1 −G3  Ck−G3 
 
 . .. .. .. ...   .. 
. . .
De esta forma, la precisión del cálculo, aplicando la teoría de funcional densidad, dependerá
primeramente de la cantidad de ondas planas que se utilicen para expandir el orbital del Kohn-
Sham, de cuántos vectores k decidamos utilizar, y del vector de onda G.
De una manera similar a la expuesta para una base de ondas planas, es posible construir
una matriz, similar a la de la ecuación (3.9), en una base de funciones localizadas, como las
de Wanier y los armónicos esféricos. Para este trabajo se usó el método de las ondas planas
3.1. Teoría de funcionales de densidad 21

lienarizadas y aumentadas LAPW, por sus siglas en inglés, el cual divide el espacio en dos
regiones: intersticial (I) y muffin-tin (MT). La región muffin-tin está conformada por esferas
centradas en los átomos, etiquetados con α, Rα , en las posiciones Rα , con radios RM T , que
contienen a los electrones más cercanos al núcleo atómico y cuyas funciones de onda puedan
ser descritas usando orbitales atómicos. Por otro lado, la región intersticial es todo el espacio
sobrante en la celda, en la que las funciones de onda de los electrones más externos son
descritos usando una combinación lineal de ondas planas.
 X


G+k
Almα ulα (rα )Ylα (r̂α ) , r α ≤ RM T
φG+k (r) =  lm (3.10)
 √1 e(G+k)·r , r ∈ I,
Ω

donde rα = r − Rα . La región muffin-tin es expandida usando esféricos armónicos Ylm ,

funciones radiales ulα , y los coeficientes AG+k
lmα , que suavizan la función de onda en el borde
de la región muffin-tin para que sea continua.
Capítulo 4

Resultados y discusión

En este capítulo se presentan los resultados obtenidos al realizar los cálculos numéricos
para la obtención de las bandas de energía en el teluro de mercurio (HgTe) y el teluro de
cadmio (CdTe) usando la teoría de funcionales de densidad como está implementada en el
software Exciting code [27].

Z Configuración electrónica Tipo

Hg 80 [Xe]6 s2 4 f14 5 d10 Metal
2 10
Cd 48 [Kr]5 s 4 d Metal
2 10 4
Te 52 [Kr]5 s 4 d 5 p No metal

Tabla 4.1: Números atómicos, configuración electrónica, y tipo de material, para el Hg, Cd y Te.

Los sistemas HgTe y CdTe son semiconductores tipo II-VI ya que están conformados por
un metal del grupo IIB y un no metal del grupo VIA que se unen formando un enlace mixto
en el que prevalece el enlace covalente. Estos semiconductores poseen una estructura tipo
blenda de zinc y pertenecen al grupo espacial F-43m en la notación de Hermann-Mauguin,
que corresponde a una red conformada por dos redes cúbicas de cara centrada monoatómicas
trasladadas 1/4 de la constante de red en la dirección [111]. Estas aleaciones poseen una
constante de red a temperatura ambiente de a = 6.492 Å y a = 6.467 Å para el CdTe y el
HgTe, respectivamente [28]. Como se puede apreciar en la tabla (4.1) el Hg y el Cd poseen dos
electrones de valencia, pertenecientes a los orbitales 6s y 5s, respectivamente, mientras que
el Te posee seis electrones de valencia, pertenecientes a los orbitales 5s y 5p. De esta forma,
cada celda posee 8 electrones de valencia.

23
24 Capítulo 4. Resultados y discusión

Una forma de saber si es necesario tomar en cuenta los efectos relativistas es considerando
la velocidad de los electrones tipo 1s en las especies atómicas. A medida que el núcleo atómico
es más pesado, los electrones más cercanos a este aumentarán su rapidez. Esta rapidez depende
del número atómico y es descrita usando la relación obtenida por Bohr [37]:

Z
 
hvs i ≈ c, (4.1)
137
donde c es la velocidad de la luz y Z el número atómico del elemento.
La velocidad que alcanza el electrón 1s en el mercurio es lo suficientemente grande como
para alterar el comportamiento de las demás bandas de energía en el átomo, al mismo tiempo
que modifica la masa efectiva del electrón mrel . Tomando en cuenta la fórmula de Bohr para
el radio atómico junto con la contracción lantánida, que se da por el llenado del orbital 4f, se
encuentra que el radio atómico del mercurio se reduce, lo que tendrá un efecto en sus bandas
de energía.
Además de la rapidez de los electrones es necesario tomar en cuenta la contribución por
elemento al gap generado por la interacción espín-órbita ∆. De la literatura se encuentra que
el teluro es el elemento que aporta más en este efecto. Estos valores se muestran en la tabla
(4.2) junto a las velocidades de los electrones en la capa s en razón de la velocidad de la luz,
y la masa relativista mrel de estos electrones en razón de su masa en reposo mrep .

vs /c mrel /mrep ∆
Hg 0.58 1.23 0.5
Cd 0.35 1.07 0.1
Te 0.38 1.08 1.1

Tabla 4.2: Velocidades y masas efectivas de los electrones tipo 1s junto con el aporte por elemento al gap de
energía, para el Hg, Cd y Te.

Conociendo los parámetros de red del HgTe y del CdTe, es posible obtener el parámetro
de incompatibilidad de ambos semiconductores y saber si es posible construir una heteroes-
tructura usando ambos materiales. En la tabla (4.3) se muestran los parámetros de red de los
semiconductores más conocidos, GaAs y AlAs, junto con su parámetro de incompatibilidad,
para compararlos con los semiconductores estudiados en este trabajo. Se observa que el pará-
metro de incompatibilidad entre el GaAs y el AlAs es de 0.12 %, siendo las heteroestructuras
formadas por estos materiales ampliamente estudiadas. Por otro lado, el parámtro de incom-
 25

patibilidad entre el HgTe y el CdTe es 0.39 %. Esto nos indica que ambos materiales poseen
parámetros de red suficientemente similares para poder construir heteroestructuras usando
ambos materiales.

Parámetro de Red
Parámetro de incompatibilidad
(Å)
GaAs 5.653
0.12 %
AlAs 5.660
HgTe 6.467
0.39 %
CdTe 6.492

Tabla 4.3: Parámetros de red junto con los parámetros de incompatibilidad de los materiales HgTe y CdTe,
estudiados en este trabajo, junto con los del GaAs y AlAs, colocados con fines comparativos.

Para construir la celda primitiva que nos permita estudiar ambos materiales es necesario
definir los vectores primitivos {a~j } de la red fcc en el espacio de posiciones. Escogiendo
expresar estos vectores como:

1 1
a~1 = aŷ + aẑ
2 2
1 1
a~2 = ax̂ + aẑ
2 2
1 1
a~3 = ax̂ + aŷ,
2 2
se puede expresar la posición de los átomos A, para el Hg(Cd); y B, para el Te, de modo que
la celda primitiva esté centrada en el Hg(Cd):

PA = 0a~1 + 0a~2 + 0a~3

1 1 1
PB = a~1 + a~2 + a~3 .
4 4 4
En la figura (4.1) se muestran las celdas convencionales del HgTe y del CdTe centradas
en los átomos de Hg o Cd. En ellas se puede apreciar la estructura tetragonal que posee la
estructura blenda de zinc.
La elección de los puntos k usados para los cálculos se realizó tomando en cuenta una red
de puntos equidistantes como lo propusieron Monkhorst y Pack [33], que debido a nuestra
26 Capítulo 4. Resultados y discusión

(a) (b)

Figura 4.1: Celdas convencionales del HgTe y CdTe.

elección de la celda primitiva se encuentra centrada en el punto Γ. Por esta razón, es necesario
dejar fuera del cálculo al punto Γ debido a que es un punto de muy alta simetría, por lo
que los estados obtenidos en cada iteración podrán presentar divergencias en los valores de
energía [34], dificultando el cálculo autoconsistente. Así, se eligió una red de puntos en la zona
de Brillouin de 8 × 8 × 8, que se redujeron a 43 puntos al realizar las operaciones de simetría
de la red de Bravais y del cristal. Se eligió además un radio muffin-tin de RM T = 2 Bohr para
todos los elementos Este valor se tómo considerando que los radios covalentes del Hg, Cd y
Te son muy similares entre sí. Por último se eligió un máximo para la base de ondas planas
al considerar RM T Gmax = 7.0.
La energía de intercambio y correlación escogida fue según la Aproximación local de densi-
dad de espín (LSDA) incluyendo la corrección de Hubbard para la interacción electrón-electrón
siguiendo los trabajos de Zhang [35], y Al-Rajoub y B. Hamad [36]. Debido a la presencia del
orbital 5d en el mercurio, como se puede apreciar en la tabla (4.1), se incluyó la corrección de
Hubbard para estos electrones tomando el valor de U = 6 Ha. Este valor se eligió de manera
empírica luego de realizar los cálculos y compararlos con la literatura, haciendo las pruebas
entre 2 y 9 Ha, teniendo que, el ordenamiento de las bandas en el HgTe no es el mismo que
se tiene en la literatura para valores U < 6, mientras que, para valores U > 6, las bandas de
energía se deforman. Además de la estructura de bandas, también se calculó la densidad de
estados para ambos casos.
Por último, se incluyó el acoplamieto espín-órbita al considerar los efectos relativistas
usando la aproximación regular de orden cero, para los electrones de valencia, y resolviendo
 27

la ecuación de Dirac para los electrones de coraza.

Una de las características más importantes a observar en las estructuras electrónicas de
ambos materiales es la ubicación y ancho del gap global de energía que, en ambos materiales,
es un gap directo ubicado en el punto Γ. Es importante resaltar que la densidad de estados
en ambos sistemas se calculó con cierto grado de exactitud para la banda de valencia, sin
embargo existe una discrepancia en la banda de conducción de CdTe y en el gap de energía
del HgTe. Esto sucede por la incapacidad de la teoría de funcionales de densidad para calcular
los estados excitados.

↓
4 4 ↑

2 2
E - EF (eV)

0 0

-2 -2

-4 -4

W L Γ X W DOS
Figura 4.2: Estructura de bandas y densidad de estados del CdTe. Las líneas rojas(negras) inidican la densidad
de estados de electrones con espín up(down).

En la figura (4.2) se presenta la estructura de bandas del CdTe junto a la densidad de

estados. Se aprecia el comportamiento de semiconductor con un gap global directo en el punto
Γ. Se aprecia además la naturaleza ligante de la banda de valencia tipo p, y la naturaleza
antiligante de la banda de conducción, tipo s. La relación de dispersión se puede aproximar a
una forma funcional cuadrática cerca del punto Γ. A medida que nos alejamos del centro de
la zona de Brillouin se rompe la degeneración de espín, en los ejes W - L y X - W. Este efecto
también se aprecia en la densidad de estados, donde se observa cómo las curvas que indican
la densidad por polarización de espín se separan en los mismos niveles de energía donde se
rompe la degeneración de bandas en la estructura electrónica.
Por otro lado, la figura (4.3) muestra la estructura de bandas y la densidad de estados
28 Capítulo 4. Resultados y discusión

↓
4 4 ↑

2 2
E - EF (eV)

0 0

-2 -2

-4 -4

W L Γ X W DOS
Figura 4.3: Estructura de bandas y densidad de estados del HgTe. Las líneas rojas(negras) inidican la densidad
de estados de electrones con espín up(down).

del HgTe. Al igual que en la estructura de bandas del CdTe se puede observar la naturaleza
ligante de la banda de valencia tipo p, y antiligante en la banda de conducción tipo s. A
simple vista el HgTe es un material sin gap de energía, por lo que se comporta como un
semimetal, en acuerdo con la literatura. Analizando la densidad de estados, el HgTe posee
una polarización de espín menos pronunciada que el CdTe, también visible en los ejes W - L
y X - W. Además se puede apreciar que la forma funcional de la banda de conducción en la
dirección L - Γ y Γ - X es casi lineal, lo que guarda relación con la velocidad relativista de los
electrones tipo s del Hg. Esta característica indica la posibilidad de encontrar portadores de
carga que se comporten como fermiones de Dirac o Weyl en la interfase de la heteroestructura
HgTe/CdTe, o al ir variando la composición química de la aleación ternaria Hg1-x Cdx Te.
Antes de calcular el gap global de energía es importante analizar el gap generado por la
interacción espin-órbita en ambas aleaciones. En las estructuras tipo blenda de zinc, este gap
de energía se puede apreciar en las bandas tipo p (l = 1). Cuando se toma en cuenta la
interacción espín-órbita es mejor usar el número cuántico momento angular total j = l + s,
por lo que, para este caso, tomará los valores j = 3/2, 1/2, siendo la interacción espín-
órbita el origen de la apertura de un gap entre ambos estados, que de otra forma estarían
degenerados. La tabla (4.4) muestra los valores del gap generado por la interacción espín-
órbita en la estructura de bandas del HgTe y del CdTe obtenidos en este trabajo. Ambos
 29

compuestos presentan valores cercanos a 1 eV. En comparación, el arsenuro de galio presenta

una apertura de bandas de energía por la interacción espín-órbita de 0.34 eV [38], casi la
tercera parte de los valores obtenidos.

ESO (Theo.) ESO (Exp.)

HgTe 1.012 1.15 [31]
CdTe 0.951 0.949 [32]

Tabla 4.4: Gaps de energía generados por el acoplamiento espin-orbita, entre los autoestados Γ8 y Γ7 , en el
HgTe y CdTe. Los valores están el eV.

Como se mencionó anteriormente, los semiconductores de los grupos III-V y II-VI poseen
una estructura blenda de zinc. Esta simetría permite conocer de antemano qué transforma-
ciones se puden aplicar a las funciones de onda en cada punto de la zona de Brillouin y, de
esta manera, clasificar a las bandas de energía usando la teoría de grupos. Así se obtiene que
el gap global de energía usualmente se calcula evaluando la expresión

Eg = E(Γ6 ) − E(Γ8 ), (4.2)

donde Γ6 y Γ8 representan los autoestados de energía, cercanos al nivel de energía de Fermi,

en el punto Γ con simetrías similares a la de los orbitales j = 1/2 y j = 3/2, respectivamente.
Es necesario entonces conocer el aporte de los orbitales sobre las bandas de energía para
comprobar lo dicho en la teoría de grupos y calcular el valor de gap de energía.
En las figuras mostradas a continuación el tamaño de los círculos indican la intensidad
del aporte de los orbitales atómicos, mientras que los colores azul y rojo indican la presencia
de orbitales tipo s o p, respectivamente. Se ve en las figuras (4.4) y (4.5) que el aporte de
los orbitales atómicos a la mínima banda de conducción se da por el elemento del grupo IIB,
mientras que el elemento del grupo VIA realiza su aporte a la máxima banda de valencia,
lo que nuevamente tiene que ver con el tipo de enlace químico entre las especies atómicas
presentes en las aleaciones.
En la figura (4.4a) se puede observar el aporte a la banda de conducción del orbital tipo s
del Cd, mientras que el aporte a la banda de valencia del orbital atómico tipo p del teluro se
muestra en la figura (4.4b). Usando el momento angular total j, la banda de conducción cerca
al punto Γ tendrá un aporte debido al orbital atómico j = 1/2, clasificada por la teoría de
grupos como Γ6 . Mientras tanto, la banda de valencia tendrá el aporte de los orbitales atómicos
30 Capítulo 4. Resultados y discusión

1
E - EF (eV)

Orbital s - Cd
0
Orbital p - Te

-1

-2

L Γ X L Γ X
(a) (b)

Figura 4.4: A la izquierda(derecha) el aporte a los autoestados de energía de los orbitales atómicos tipo s(p)
en las bandas de energía del CdTe. El color azul(rojo) representa al orbital tipo s(p), y el tamaño del círculo
la intensidad del aporte.

j = 3/2, correspondiente a las dos bandas más cercanas al nivel de Fermi degeneradas en el
punto Γ, y del orbital j = 1/2, correspondiente a la banda inferior de la figura (4.4b), ambas
clasificadas como Γ8 y Γ7 , respectivamente. La diferencia de energía E(Γ8 ) − E(Γ7 ) se mostró
en la tabla (4.4). Este ordenamiento de las bandas de energía es el que se encuentra usualmente
en las estructuras semiconductoras III-V y II-VI, lo que obedece al tipo de enlace existente
entre las especies atómicas presentes en las aleaciones, por lo que se esperaría encontrar el
mismo ordenamiento en el HgTe.
En la figura (4.5) se aprecia el aporte de los orbitales atómicos de ambas especies atómicas
cerca del gap de energía para el HgTe. A diferencia del CdTe, existe una superposición de los
orbitales tipo s (Γ6 ) del mercurio y p (Γ8 ) del teluro sobre la banda de conducción. El aporte
del orbital tipo s se hace más intenso a medida que nos alejamos del punto Γ, mientras que
el aporte del orbital tipo p se hace más intenso al acercarnos al mismo punto. Este aporte
se genera por la interacción espín-órbita existente, la cual incrementa la diferencia entre los
autoestados de energía j = 3/2 y j = 1/2, de los orbitales l = 1, llegando a colocar los
autovaloder de energía de los estados Γ8 por encima del de los estados Γ6 . De esta forma, el
HgTe es un material sin gap de energía debido a la degeneración de los estados j = 3/2 en el
punto Γ.
 31

Este fenómeno se conoce como una inversión de bandas de energía, ya que como se ve
en la figura (4.5), el autoestado con aporte del orbital tipo p (Γ8 ) se encuentra por encima
del autoestado con aporte del orbital tipo s (Γ6 ), lo que da lugar a un gap global de energía
"negativo.en el HgTe al aplicar la ecuación (4.2). Es imprtante notar que el autostado Γ7 no
participa directamente de esta inversión y mantiene el mismo orden en el CdTe y en el HgTe.

1
E - EF (eV)

Orbital p - Te
0
Orbital s - Hg

-1

-2

L Γ X L Γ X
(a) (b)

Figura 4.5: A la izquierda(derecha) el aporte a los autoestados de energía de los orbitales atómicos tipo s(p)
en las bandas de energía del HgTe. El color azul(rojo) representa al orbital tipo s(p), y el tamaño del círculo
la intensidad del aporte.

Eg (Teo.) Eg (Exp.)
HgTe -0.10404 -0.32 [31]
CdTe 1.05134 1.513 [32]

Tabla 4.5: Gaps de energía teóricos y experimentales, entre los autoestados Γ6 y Γ8 , del HgTe y CdTe. Los
valores están en eV.

En la tabla (4.5) se muestran los valores de los gaps globales de energía obtenidos en
este trabajo usando la ecuación (4.2) y se comparan con los valores experimentales que se
encuentran en la literatura. Como se puede apreciar, existe una diferencia notable entre el valor
del gap de energía experimental, dado por la literatura, y los obtenidos teóricamente. Esto
sucede, y es importante incidir en este punto, debido a que el DFT es una teoría que obtiene
32 Capítulo 4. Resultados y discusión

exactamente el estado fundamental de un sistema, mas no así los estados excitados. Los estados
excitados requieren de un tratamiento especial para su obtención, como la aproximación
GW o la teoría de funcionales de densidad dependiente del tiempo, sin embargo, el costo
computacional y el tiempo de cálculo se incrementan significativamente.
Capítulo 5

Conclusiones

Una de las formas para estudiar las propiedades de los sólidos de manera teórica es la
Teoría de Funcionales de Densidad. Sin embargo, esta técnica no otorga un valor exacto para
la energía de los electrones en las bandas de conducción por lo que necesita correcciones extra
que amuntan el tiempo de cálculo significativamente.
Es necesario conocer las diferentes características de los átomos con los que se trabajan,
en especial si poseen una masa atómica relativamente grande y presenten orbitales tipo d.
Además los efectos relativistas que surgen, como el acoplamiento espín-órbita, son de gran
importancia para entender el comportamiento de los electrones en el material.
El CdTe y el HgTe presentan un gap en las bandas tipo p, generado por el acoplamiento
espín-órbita, entre los autoestados Γ8 y Γ7 , de 0.95 y 1.01 eV, respectivamente. Estos valores
son relativamente altos comparándolos con otros semiconductores.
La intensa interacción espín-órbita presente en ambas estructuras es responsable de un or-
denamiento de bandas diferente en el HgTe, al compararlo con el CdTe. De esta forma el CdTe
y el HgTe presentan gaps de energía, E(Γ6 ) − E(Γ8 ), de 1.051 y -0.104 eV, respectivamente.
Por este motivo se considera al HgTe como un material con gap de energía invertido.
Si el HgTe y el CdTe se unen formando una heteroestructura, lo que es posible al tomar
en cuenta el parámetro de incompatibilidad η, el gap entre los estados Γ6 y Γ8 se tendrá que
cerrar en la interfase, formando estados de superficie metálicos. Este efecto también se puede
obtener al variar la fórmula estequiométrica Hg1-x Cdx Te. Además, si se forma un pozo de
potencial de CdTe/Hgte/CdTe, es posible obtener un estado de efecto Hall cuántico de espín.
Estos casos son el objetivo del trabajo de licenciatura.

33
Bibliografía

[1] M. Z. Hasan y C. L. Kane; Colloquium: topological insulators, Rev. Mod. Phys. 82(4),
3045 (2010).

[2] J. E. Avron, D. Osadchy y R. Seiler; A Topological Look at the Quantum Hall Effect,
Phys. Today 56 38 (2003)

[3] T. Ando, Y. Matsumoto y Y. Uemura; Theory of Hall effect in a two-dimensional

electron system, J. Phys. Soc. Japan, 39, 279-288 (1975).

[4] K. V. Klitzing, G. Dorda y M. Pepper; New method for high-accuracy determination of

the fine-structure constant based on quantized Hall resistance, Phys. Rev. Lett., 45, 494
(1980).

[5] R. B. Laughlin; Quantized Hall conductivity in two dimensions, Phys. Rev. B, 23, 5632
(1981).

[6] Q. Niu, D. J. Thouless, y Y. S. Wu; Quantized Hall conductance as a topological inva-

riant, Phys. Rev. B, 31, 3372 (1985).

[7] B. A. Bernevig, T. L. Hughes y S. C. Zhang; Quantum spin Hall effect and topological
phase transition in HgTe quantum wells, Science, 314, 1757-1761 (2006).

[8] M. Konig, S. Wiedmann, C. Brüne, A. Roth, H. Buhmann, L. W. Molenkamp, ... y

S. C. Zhang; Quantum spin Hall insulator state in HgTe quantum wells, Science, 318,
766-770 (2007).

[9] International Union of Crystallography; Report of the Executive Committee for 1991,
Acta Crystallogr. A., 48 (6): 922-946.

35
36 Bibliografía

[10] P. Hofmann; Solid state physics: an introduction, segunda edición, Wiley-VCH Verlag
GmbH & Co., Wenheim (2015).

[11] H. Ibach y H. Lüth; Solid-State Physics: An Introduction to Principles of Materials

Science, cuarta edición, Springer Science & Business Media, Berlin (2009).

[12] S. H. Simon; The Oxford solid state basics, primera edición, Oxford University Press,
Oxford (2013).

[13] W. Nolting; Fundamentals of many body physics, primera edición, Springer-Verlag, Ber-
lin (2009).

[14] C. Hamaguchi y C. Hamaguchi; Basic semiconductor physics, segunda edición, Springer-

Verlag, Berlin (2010).

[15] W. R Frensley; Heterostructure and quantum well physics, en VLSI Electronics Micros-
tructure Science, Vol. 24, pp. 1-24, Elsevier (1994).

[16] B. R. Nag; Physics of quantum well devices, Vol. 7, Springer Science & Business Media
(2001).

[17] T. Ando, A. B. Fowler y F. Stern; Electronic properties of two-dimensional systems,

Rev. Mod. Phys., 54, 437 (1982).

[18] P. Yu y M. Cardona; Fundamentals of Semiconductors: Physics and Materials Proper-

ties, cuarta edición, Springer-Verlag, Berlín (2010).

[19] R. J. Martín-Palma, F. Agulló-Rueda y J. Martínez-Duart; Nanotechnology for micro-

electronics and optoelectronics, primera edición, Elsevier, Oxford (2006).

[20] T. Lancaster, y S. J. Blundel; Quantum field theory for the gifted amateur, primera
edición, OUP Oxford, Oxford (2014).

[21] B. R. Desai; Quantum mechanics with basic field theory, primera edición, Cambridge
University Press, Cambridge (2010).

[22] F. R. Villatoro; Fermiones de Weyl y semimetales de Weyl, recuperado de: https://

francis.naukas.com/2015/07/23/fermiones-de-weyl-y-semimetales-de-weyl/
2015
Bibliografía 37

[23] P. Hohenberg y W. Kohn; Inhomogeneus electron gas, Phys. Rev., 136, 136 (1964).

[24] R. M. Martin; Electronic structure: Basic theory and practical metods; primera edición,
Cambridge University Press, Cambridge (2004).

[25] K. Burke; The abc of dft, Department of Chemistry, University of California (2007).

[26] J. P. Perdew y Y. Wang; Accurate and simple analytic representation of the electron-gas
correlation energy, Phys. Rev. B, 45, 23 (1992).

[27] A. Gulans, S. Kontur, C. Meisenbichler, D. Nabok, P. Pavone, S. Rigamonti, S. Sagmeis-

ter, U. Werner y C. Draxl; exciting: a full-potential all-electron package implementing
density-functional theory and many-body perturbation theory, J. Phys.: Condens. Matter,
26 (2014)

[28] P. Capper y J. Garland; Mercury Cadmium Telluride: Growth, Properties and Applica-
tions, primera edición, John Wiley & Sons, Wst Sussex (2011).

[29] E. Grilli, M. Guzzi, R. Zamboni y L. Pavesi; High-precision determination of the tem-

perature dependence of the fundamental energy gap in gallium arsenide, Phys. Rev. B,
45, 1638 (1992).

[30] S. Adachi; GaAs, AlAs, and Al x Ga1-x As: Material parameters for use in research and
device applications, J. Appl. Phys., 58, R1-R29 (1985).

[31] N. Orlowski, J. Augustin, Z. Gołacki, C. Janowitz y R. Manzke; Direct evidence for the
inverted band structure of HgTe, Phys. Rev. B, 61, R5058 (2000).

[32] G. Fonthal, L. Tirado-Mejıa, J. I. Marın-Hurtado, H. Ariza-Calderón y J. G. Mendoza-

Alvarez; Temperature dependence of the band gap energy of crystalline CdTe. J. Phys.
Chem. Solids, 61, 579-583 (2000).

[33] H. J. Monkhorst y J. D. Pack; Special points for Brillouin-zone integrations, Phys. rev.
B, 13, 5188.(1976).

[34] A. Baldereschi y E. Tosatti; Mean-value point and dielectric properties of semiconductors

and insulators, Phys. Rev. B, 17, 4710. (1978).
38 Bibliografía

[35] H. Zhang y S. C. Zhang; Topological insulators from the perspective of first-principles

calculations. Phys. Status Solidi RRL, 7, 72-81 (2013).

[36] S. Al-Rajoub, B. Hamad; Theoretical investigations of the structural, electronic and

optical properties of Hg1-x Cd x Te alloys, Philos. Mag, 95, 2466-2481 (2015).

[37] L. J. Norrby, Why is mercury liquid? Or, why do relativistic effects not get into chemistry
textbooks?. J. Chem. Educ., 68, 110 (1991).

[38] H. Mayer y U. Rössler; Spin splitting and anisotropy of cyclotron resonance in the
conduction band of GaAs, Phys. Rev., 44, 9048 (1991).

[39] O. Madelung, U. Rösslery M. Schulz; Landolt-Börnstein - Group III Conden-

sed Matter, volumen 41B, Springer-Verlag Berlin Heidelberg, Berlin (1999).