Архитектура Аудит Военная наука Иностранные языки Медицина Металлургия Метрология
Образование Политология Производство Психология Стандартизация Технологии


Матрица зонной структуры кристаллов со структурой алмаза



При составлении матрицы тензора взаимодействия необходимо учитывать следующие правила отбора по перестановочной и инверсионной симметрии состояний:

Перстановочная симметрия.

Поскольку любое взаимодействие в орбитальном пространстве является скаляром в пространстве спина, то эти взаимодействия отличны от нуля только для состояний с одинаковой спиновой перестановочной симметрией.

· Инверсионная симметрия.

Чётное взаимодействие имеет отличные от нуля матричные элементы между состояниями одинаковой чётности, не чётное взаимодействие – между состояниями противоположной чётности.

В Табл. 2.2 приведена матрица взаимодействия для основных состояний зонной структуры кристаллов типа алмаза. В качестве базисных функций были взяты базисные функции нормированных представлений группы куба.

 

Табл. 2.2. Матрица взаимодействия для основных состояний зонной структуры кристаллов типа алмаза.

E1 T1kx T1ky T1kz
E25 Qky Qkx Pkz Ö 2R kz P2kz
E25 Qkz Qky Pkx -1/Ö 2R kx -Ö 1.5R kx P2kx
E25 Qkz Qkx Pky -1/Ö 2R ky Ö 1.5R ky P2ky
T1kx Qky Qkz E15 Tkx Q1ky Q1kz
T1ky Qkx Qkz E15 Tky Q1kx Q1kz
T1kz Qky Qkx E15 Tkz Q1ky Q1kx
Pkz Pkx Pky E2 P1kz P1kx P1ky
Tkx Tky Tkz E1(2)
Ö 2R kz -1/Ö 2R kx -1/Ö 2R ky E12 Ö 2R1 kz -1/Ö 2R1 kx -1/Ö 2R1 ky
-Ö 1.5R kx Ö 1.5R ky E12 -Ö 1.5R1 kx Ö 1.5R1 ky
Q1ky Q1 kx P1kz Ö 2R1 kz E25(2) P3kz
Q1 kz Q1 ky P1kx -1/Ö 2R1 kx -Ö 1.5R1 kx E25(2) P3kx
Q1kz Q1 kx P1ky -1/Ö 2R1 ky Ö 1.5R1ky E25(2) P3ky
P2kz P2kx P2ky P3kz P3kx P3ky E2(2)

 

Введены следующие обозначения:

Е1=Е(Г1+)+kx2+ ky2+ kz2,

E25= Е(Г25+)+kx2+ ky2+ kz2,

E15=Е(Г15-)+kx2+ ky2+ kz2,

E2= Е(Г2-)+kx2+ ky2+ kz2,

k=2p/a*ab*ki,

ki=0.. 1 [001], ki=0.. 1/2 [111],

где а – параметр кристаллической решётки, аb-Боровский радиус.

Индексами в скобках отмечены повторяющиеся однотипные состояния.

Значения недиагональных матричных элементов даны в таблице 3.

 

Матричные элементы Кремний (Ry) Германий (Ry) Si0.25Ge0.75 (Ry)
P P1 P2 P3 Q Q1 R R1 T T1 1.2 -0.09 0.1 1.32 1.05 -0.807 0.83 1.21 1.08 0.206 1.36 0.1715 0.1 1.6231 1.07 -0.752 0.8049 1.4357 1.2003 0.5323 1.32 0.1061 0.1 1.5473 1.065 -0.7658 0.8112 1.3793 1.1702 0.4507

Таблица 3.Значения недиагональных матричных элементов.

2.8 Зонная структура кремния в направлениях [001], [111]

Для расчёта зонной структуры в направлении [001]вектора k надо положить kx= ky=0 и рассчитывать энергии относительно проекции kz: kz=2p/ab*a*ki, где ki=0.. 1. Для расчёта зонной структуры в направлении [111] нужно положить kx= ky= kz= k/Ö 3 и рассчитывать матрицу относительно величины k: k=2p/ab*a*ki, где ki=0.. ½.

Ниже приводятся данные расчёта зонной структуры кремния в направлении [001] и [111] в единицах электрон-вольт (отрицательные значения ki в направлении [111] приняты условно).

 

Табл. 2.3. Зонная структура кремния в направлении [001] и [111]

(единицы электрон-вольт).

 

k Г1+ Г25+ Г25+ Г25+ Г15- Г15- Г15- Г2- Г1+ (2) Г12- Г12-
-0, 5 -11, 71 -4, 84 -0, 90 -0, 90 2, 33 4, 45 4, 45 8, 18 8, 78 9, 27 9, 27
-0, 4 -12, 14 -3, 74 -0, 70 -0, 70 2, 55 4, 35 4, 35 7, 13 8, 87 8, 87 8, 93
-0, 3 -12, 48 -2, 59 -0, 48 -0, 48 2, 82 4, 11 4, 11 6, 02 8, 75 8, 75 8, 83
-0, 2 -12, 72 -1, 45 -0, 25 -0, 25 3, 11 3, 80 3, 80 4, 95 8, 01 9, 09 9, 09
-0, 1 -12, 87 -0, 46 -0, 07 -0, 07 3, 36 3, 53 3, 53 4, 01 7, 33 9, 46 9, 46
-12, 92 0, 00 0, 00 0, 00 3, 43 3, 43 3, 43 3, 60 7, 07 9, 66 9, 66
0, 1 -12, 87 -0, 30 -0, 16 -0, 16 3, 27 3, 64 3, 64 3, 90 7, 33 9, 21 9, 71
0, 2 -12, 72 -1, 00 -0, 55 -0, 55 2, 90 4, 19 4, 19 4, 54 8, 01 8, 30 9, 86
0, 3 -12, 48 -1, 85 -1, 03 -1, 03 2, 47 4, 93 4, 93 5, 03 7, 66 8, 96 10, 11
0, 4 -12, 14 -2, 77 -1, 50 -1, 50 2, 06 4, 75 5, 78 5, 78 7, 93 10, 10 10, 47
0, 5 -11, 70 -3, 69 -1, 94 -1, 94 1, 71 4, 02 6, 69 6, 69 8, 83 10, 93 11, 38
0, 6 -11, 16 -4, 61 -2, 32 -2, 32 1, 45 3, 25 7, 65 7, 65 9, 92 11, 49 12, 79
0, 7 -10, 53 -5, 51 -2, 62 -2, 62 1, 27 2, 58 8, 66 8, 66 11, 05 12, 15 14, 31
0, 8 -9, 80 -6, 39 -2, 85 -2, 85 1, 20 2, 04 9, 71 9, 71 12, 19 12, 91 15, 95
0, 9 -8, 98 -7, 22 -2, 98 -2, 98 1, 23 1, 64 10, 81 10, 81 13, 35 13, 78 17, 70
-8, 06 -8, 01 -3, 03 -3, 03 1, 36 1, 37 11, 97 11, 97 14, 54 14, 75 19, 55

 

 

k Г25+ (2) Г25+ (2) Г25+ (2) Г2- (2)
-0, 5 13, 12 18, 14 18, 14 20, 47
-0, 4 11, 69 16, 61 16, 61 18, 71
-0, 3 11, 01 15, 23 15, 23 17, 10
-0, 2 11, 32 14, 04 14, 04 15, 64
-0, 1 12, 07 13, 15 13, 15 14, 34
12, 78 12, 78 12, 78 13, 46
0, 1 12, 45 12, 88 12, 88 14, 50
0, 2 12, 45 13, 18 13, 18 16, 22
0, 3 12, 76 13, 68 13, 68 18, 21
0, 4 13, 28 14, 38 14, 38 20, 38
0, 5 14, 01 15, 28 15, 28 22, 70
0, 6 14, 96 16, 38 16, 38 25, 16
0, 7 16, 13 17, 66 17, 66 27, 74
0, 8 17, 51 19, 12 19, 12 30, 43
0, 9 19, 11 20, 76 20, 76 33, 25
20, 88 22, 54 22, 54 36, 18

 


Рис. 2.1. Зонная структура кремния в направлениях[001], [111]

Данные расчёта энергетической структуры кремния в области минимума энергии зоны проводимости для значений ki=0.8..0.85.

 

ki E(эв.) ki E(эв.) ki E(эв.)
0, 8 1, 1986 0, 82 1, 1958 0, 84 1, 1972
0, 805 1, 1975 0, 825 1, 1958 0, 845 1, 1982
0, 81 1, 1967 0, 83 1, 1960 0, 85 1, 1995
0, 815 1, 1961 0, 835 1, 1965    

 

Энергия минимума зоны проводимости в зависимости от ki приведена на Рис. 2.2. Энергия в минимуме равна 1.1957эв. при значении квазиимпульса ki=0.824.

 
 

Рис. 2.2. Энергия дна зоны проводимости кремния от ki [001].

 

 
 

Рис. 2.3. Энергия дна зоны проводимости кремния в направлении [100].

Формулы зависимости энергии зоны проводимости от величины квазиимпульса для двух направлений представлены ниже:

E(ki[001])=–1.2·105(ki–0.8235)4+33.33·(ki–0.8235)3+5.328·(ki–0.8235)2+ 1.195761

( 2.4) E(ki[100])= -128·ki4 +26.36·ki2 +1.19576.

По этим зависимостям можно определить эффективную массу электронов:

( 2.5) ,

где а, а0- параметр решётки полупроводникового материала и Боровский радиус, Е-энергия в атомных единицах.

Для кремния эти формулы будут:

( 2.6) .

Вторая производная энергии по квазиимпульсу определяется по формулам ( 2.4 ) для двух направлений [001], [100].

 

 
 

Откуда эффективные массы будут:

Эффективная масса определяется относительно массы свободного электрона.

Поскольку эффективная масса это вторая производная скалярной величины по полярному вектору k , тоона является тензором второго ранга с тремя главными значениями m[001], m[100]= m[010].

 

Данные расчёта энергетической структуры кремния в области максимума энергии валентной зоны для значений ki от 0 до 0.1 в двух направлениях зоны [001], [111] даны на Рис. 2.5.

 

Зависимости энергии этих состояний от величины квазиимпульса для направлений [001], [00-1]:

( 2.7)


Рис. 2.4. Энергия валентной зоны кремния в направлениях [111], [001].

 

Массы дырок двух типов рассчитываются аналогично массы электронов:

 
 

mp1=-0.374/1.08=-0.346, mp2=-0.374/2.58=-0.145.

Из рис.4 видно, что зависимости энергии валентной зоны в направлении [001], [111] имеют некоторую ассимметрию зависимости. Аналитическая зависимость энергии валентной зоны в этом направлении для двух типов электронов в валентной зоне будет:

( 2.8)

Вторая производная энергии по квазиимпульсу и массы двух типов будут:

 

 
 

Экспериментально по наблюдению циклотронного резонанса эффективная масса двух типов дырок равна 0.5 и 0.16, что находится в согласии со средними величинами масс (0.346, 0.7).

Следует отметить, что вторые производные энергии максимума валентной зоны отрицательны. Это приводит к тому, что эффективные массы электронов в валентной зоне отрицательны, поэтому переходят к формализму дырок. Дырка с положительной массой и положительным зарядом ведёт себя также как и электрон с отрицательной массой и отрицательным зарядом.

Зонная структура германия

Расчёт зонной структуры германия приведён в ( 2.7):

Табл. 2.4. Зонная структура германия (эв.)

 

k Г1 Г25 Г25 Г25 Г2 Г15 Г15 Г15 Г1 Г12 Г12 Г25 Г25 Г25 Г2
-0, 5 -11, 6 -7, 8 -1, 2 -1, 2 0, 84 4, 8 4, 8 10, 3 12, 7 12, 7 12, 8 21, 0 28, 3 28, 3 33, 2
-0, 4 -11, 8 -7, 2 -1, 1 -1, 1 0, 94 4, 9 4, 9 10, 1 10, 9 10, 9 11, 0 17, 3 25, 0 25, 0 29, 5
-0, 3 -12, 3 -5, 7 -0, 9 -0, 9 1, 22 4, 7 4, 7 7, 9 9, 8 9, 8 11, 3 14, 5 22, 0 22, 0 26, 0
-0, 2 -12, 7 -3, 8 -0, 6 -0, 6 1, 62 4, 1 4, 1 5, 7 9, 7 9, 7 10, 2 14, 1 19, 5 19, 5 22, 9
-0, 1 -13, 0 -1, 7 -0, 2 -0, 2 1, 84 3, 5 3, 5 3, 8 8, 5 10, 2 10, 2 15, 6 17, 7 17, 7 20, 2
-13, 1 0, 0 0, 0 0, 0 0, 99 3, 2 3, 2 3, 2 7, 8 10, 5 10, 5 17, 0 17, 0 17, 0 18, 4
0, 2 -13, 0 -1, 7 -0, 6 -0, 6 2, 69 2, 8 4, 0 4, 0 8, 6 9, 5 10, 7 16, 2 17, 3 17, 3 21, 0
0, 4 -12, 6 -3, 8 -1, 5 -1, 5 1, 89 4, 4 5, 8 5, 8 8, 1 10, 7 11, 2 16, 5 18, 2 18, 2 25, 3
0, 6 -11, 9 -5, 8 -2, 3 -2, 3 1, 30 3, 9 7, 9 7, 9 9, 1 12, 2 13, 4 17, 7 19, 7 19, 7 30, 3
0, 8 -10, 9 -7, 8 -2, 8 -2, 8 1, 10 2, 5 10, 2 10, 2 11, 6 13, 5 16, 6 19, 7 21, 8 21, 8 35, 8
-9, 8 -9, 7 -3, 0 -3, 0 1, 38 1, 4 12, 8 12, 8 14, 3 15, 2 20, 2 22, 4 24, 5 24, 5 41, 9

 

Графически эти зависимости представлены на Рис. 2.5.


Рис. 2.5. Энергия уровней зонной структуры германия в двух направлениях [001], [111].

 

Для германия минимум зоны проводимости наблюдается в направлении [111] почти на границе зоны. В Табл. 2.5 приведена энергия зоны проводимости в области минимума.

Табл. 2.5. Энергия зоны проводимости в области минимума энергии

 

k[111] Г2(эв.)
-0, 4 0, 94123
-0, 41 0, 92276
-0, 42 0, 90617
-0, 43 0, 89151
-0, 44 0, 87878
-0, 45 0, 86802
-0, 46 0, 85924
-0, 47 0, 85246
-0, 48 0, 84770
-0, 49 0, 84498
-0, 5 0, 84432

 

 

На Рис. 2.6 показана зависимости энергии минимума зоны проводимости в увеличенном масштабе.

 

 


Рис. 2.6. Энергия минимума зоны проводимости германия в направлении [111]Г2(эв.)

 

Аналитически зависимость энергии уровня Г2 от k в области минимума энергии представляется следующим степенным рядом:

(2.9) .

 

Эффективная масса электронов в зоне проводимости рассчитывается для параметра кристаллической решётки 0.566 нМ и второй производной энергии вблизи минимума зоны проводимости по формуле:

(2.10) .

 

 
 

Вторая производная определяется по уравнению (2.9):

 

Точность определения эффективной массы в этом случае недостаточна, так как при получении степенного ряда разложение определялось по одной половине минимума энергии.

 

Область энергии вблизи максимума валентной зоны более подробно показан на Рис. 2.7.

 

 
 

 

 

Рис. 2.7. Валентная зона германия в области максимума энергии для двух направлений зоны: [111], [001].

 


Поделиться:



Последнее изменение этой страницы: 2017-04-12; Просмотров: 376; Нарушение авторского права страницы


lektsia.com 2007 - 2024 год. Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав! (0.044 с.)
Главная | Случайная страница | Обратная связь