Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Арсенид галлия. Получение, свойства и применение

.pdf
Скачиваний:
39
Добавлен:
25.10.2023
Размер:
18.36 Mб
Скачать

80

З О Н Н А Я С Т Р У К Т У Р А

[ГЛ. 2

Рассматривая

антисимметричную часть

потенциала,

а также разность между симметричными частями потен­ циалов структур алмаза и сфалерита как добавочный воз­ мущающий потенциал, Херман сделал некоторые пред­ сказания о зонной структуре ряда полупроводниковых соединений со сфалерптпой структурой.

Успеху метода Хермана содействует близость потенциа­ лов, создаваемых в решетке элементами — соседями в пе­

риодической системе.

Общими выводами этого

метода

применительно к арсеииду галлия являются:

 

1. В первом порядке теории возмущений снимается вы­

рождение

состояний

\VX и \ ¥ г в

углах зоны Бриллюэна,

а также Zx

вдоль дпагонали квадратной грани и состояний

Хг и X j в центре квадратных граней *). Последнее

обстоя­

тельство приводит к

появлению

запрещенного

участка

внутри валентных зон.

2. Во втором порядке теории возмущений происходит сдвиг краев зон, причем верхипй край валентиой зоны сдвигается вниз, а нижний край зоны проводимости — вверх, в результате чего ширина запрещенной зоны растет при переходе от AIV к соответствующему изоэлектроииому аналогу AIUBVB горизонтальных рядах периодической си­ стемы. Коллауэй [6] на основе метода Хермана исследо­ вал зонную структуру арсенида галлия, выводя ее из из­ вестной зонной структуры гермапия.

Количественные расчеты Коллауэя в общем подтвер­ дили теоретические выводы Хермана с некоторыми уточ­ нениями и добавлениями. Оказалось, что изменение ширины запрещенной зоны при переходе от германия к арсениду галлия происходит в основном за счет сдвига валент­ ной зоны; кроме того, основной минимум зоны проводи­ мости перемещается из центра гексагональной грани зоны Бриллюэна в случае германия в центр зоны Бриллюэна в случае арсенида галлия.

Верхний край валентной зоны, по-видимому, остается в центре зоны Бриллюэна (без учета спин-орбитальных эффектов) или же может быть сдвинут к точке W. В основ­ ном все минимумы зоны проводимости имеют тенденцию подниматься, за исключением минимума {100}.

*) Обозначения представлений по Херрннгу [5].

2.3]

П Р И М Е Н Е Н И Е к р - М Е Т О Д А

81

 

Общие качественные выводы о зонной структуре

АШВУ',

в том числе и арсенида галлия, были сделаны также в ра­ боте Губанова и Нраньяна 17]. Они установили располо­ жение трехкратно вырожденного максимума валентной зоны в центре зоны Бриллюэиа и указали па наличие за­ прещенной зоны внутри валентных зон для кристаллов со структурой сфалерита.

2.3. Применение kp-метода для исследования спектра электропов вблизи симметричных точек

Из изложенного выше видно, что в общем зонная струк­ тура арсенида галлия должна быть похожа на зонную структуру германия. В частности, верхний край валентпой зоны в центре зоны Бриллюэиа трехкратно вырожден без учета спина и принадлежит к состояниям связи, по­ строенным на функциях, обладающих симметрией ( X , Y, Z); нижний край зоны проводимости в той же точке об­ ладает сферической симметрией «S-состояния. Кроме того, имеются дополнительные минимумы в направлениях {100}

и{111}.

Вработах Дрессельхауза [2] и Кэйна [8, 9] рассмот­ рен вид завнеимости энергии от волнового вектора вблизи симметричных точек зоны Бриллюэиа kp-методом. Идея метода состоит в том, что решение уравнения Шредингера для блоховского множителя Uk(r) ищется в виде разло­

жения по блоховским множителям симметричной точки к0 . При этом член/1 (к—к0 )р/т в гамильтониане рассматрива­ ется в качестве возмущения. Этот метод позволяет получить закон дисперсии для электрона вблизи рассматриваемой симметричной точки, kp-метод распространен п на случай учета спии-орбитального взаимодействия. Особая теория возмущений для состояний в центре зоны Бриллюэна была применена Кэйном, который точно учел взаимодей­ ствие верхней вырожденной валентной зоны с нижней зо­ ной проводимости п приближенно — вклад от остальных зон. Рассматривая совместно трехкратно вырояеденную ва­

лентную зону Г 1 5

и невырожденную зону проводимости Гх

в точке k = 0 , ои

получил секуляриый определитель 4-го

порядка для нахождения энергии. Коэффициенты, входя­ щие в этот секулярный определитель, можно найти из эксперимента. В этом смысле kp-метод является полуэм-

6 Арсишд галлля

82

З О Н Н А Я С Т Р У К Т У Р А

[ГЛ. 2

пирическим.

Впоследствии Поллак и Кардона

[10] обоб­

щили этот метод для нахождения зависимости энергии от волнового вектора по всей зоне Бриллюэна.

Учет спина удваивает порядок секулярпого опреде­ лителя и вводит дополнительные величины, подлежащие определению. Появляется линейный член в спектре, бла­ годаря которому экстремум валентной зоны смещается от центра. Количественно этот эффект незначителен *).

Аналитическое решепие в общем виде пз секулярного определителя 8-го порядка получить ие удается. Поэтому Кейи предложил приближенно диагонализпровать секулярный определитель, учитывая точно только взаимодей­ ствие трех зон — зоны проводимости, зоны легких дырок и зоны тяжелых дырок, отщепленной от верхней ва­ лентной зоны благодаря сшш-орбитальному взаимодей­ ствию. В этом приближении решение имеет довольно слож­ ный вид. Это решение до сих пор практически не исполь­

зовалось. Эффективные массы в

точке к = 0

для зоны про­

водимости—т с , зоны тяжелых

д ы р о к — т ^ н , зоны легких

дырок — mih и отщепленной благодаря

спин-орбнталь-

ному

взаимодействию зоны

тяжелых

дырок — п г в ^

имеют

вид

 

 

m h h

 

6

 

 

л

 

-±- = l + M + {L-M-N)

 

k*kl +k*f>+ klkl

f

_ L = 1 I L +2M

*

A

p

 

msh

i +

3

3

Bg+Ar'

 

Выражения

для

m h h и mih

незначительно

отличаются

от точных эффективных масс, полученных диагонализацией полного определителя 8-го порядка. В направлениях

{100} и {111} они точны. При малых значениях к

точный

закон дисперсии для тяжелых и

легких дырок имеет вид

е (к) = Ак2 ± Ув2+ С4 [к\к\

+ k2xk2z + к2ук\)

,

*) Иногда этому линейному члену приписывают даже наблю­ даемые эффекты, см., например, [11].

2.3]

П Р И М Е Н Е Н И Е к р - М Е Т О Д А

83

где

Для арсенида галлия величина снин-орбитальиого рас­ щепления зон Д0 =0,34 эв [12, 13]. Эта величина нахо­ дится по формуле [13].

д Ш - V _ 2 . Щ

1 . у

0

~з~ р

" Т

р '

где Д Р вычисленные

значения

спин-орбитального

расщепления р-соетояний валентных электронов свобод­

ного атома с поправкой на

возмущение

кристаллическим

полем.

Последняя поправка вводится в виде множителя:

 

Д Ge /Д Ge ~ 3/2.

 

 

 

 

крист/

атом

'

 

 

 

Параметры L , М,

N и F были оценены Кардоной [14],

исходя

из экспериментальных

значений

энергетических

расстояний в точкек=0 для большинства соединений

АШВУ.

Параметр F определяется

взаимодействием

валентной

зоны Г 1 5 0 и зоны проводимости Г 1 с

и может быть

выражен

через матричный элемент Р:

 

 

 

 

 

 

Р* [ 8 р 1 5 с ) - 6 р

1 5 д ) + в (Г1 5 ) - г 2 5 )]

 

 

 

2 ^ [ е р ( Г 1 5 с ) - е ( Г 1 5 1 ) ) ]

 

 

 

где е р

1 5 ) означает

положение

соответствующей

зоны

в арсениде галлия, е (Г1 5 )— в германии, индексы у и с относятся к валентной зоне и зоне проводимости соответ­ ственно.

Матричный

элемент, связывающий зоны

Г 1 5 В и

Г 1 с ,

равен

 

 

 

 

 

 

 

 

p i , =

J r l < r i o i ^ i r 1 5 u > | 2 .

 

 

и ГПараметр1 5 С

М

определяется

взаимодействием зон

Г 1 5 в

 

М =

 

С2

 

 

Параметр L

= F-\-2G,

е р ( Г 1 5 с ) - 6 р ( Г 1 5 « ) ] '

взаимодейст­

a G

определяется

вием зон Г1 5 „ и Г 1 2 С :

 

 

 

 

r _ _ ,

R

« р ( Г 1 5 с ) - е р ( Г 1 5 Р ) + е ( Г 1 5 ) - в ( Г 2 5 )

 

 

'

 

2 [ 6 р ( Г 1 5 с ) - 6 р ( Г 1 5 и ) ]

 

 

6*

 

 

 

 

 

 

S4

З О И И Л Я С Т Р У К Т У Р А

[ГЛ. 2

При этом для матричных элементов Р и Q брали зиачеппя Р2 = 23 эв и Q°- — 15,5 эв, используя соответствую­ щие данные для германия. Величина 1,6 также найдена пз значения G для германия.

Наконец, параметр N равен

N — F — G -\-

М,

 

 

если пренебречь вкладом высоколежащпх зон

Г 2 5 , кото­

рые отбрасывались также п при

вычислении

М. Вели­

чины F, G п М в таблицах Кардоиы не приведены, однако

их легко заново вычислить из данных табл.

1 работы [14]

и вышеприведенных формул. Для арсеиида галлия:

F = — 12,8;

£ = - 1 5 , 5 6 ;

N = -

14,11;

G = - l , 3 8 ;

Л / = - 3 , 6 9 .

 

 

 

Параметры L , М,

N были оценены также

Закировым

и Щтивельмаиом [15]. Их оценки

близки к

приведенным

в тексте. Они также обнаружили пепараболпчпость зоны тяжелых дырок.

Приведенные формулы дают зпачения эффективных

масс в точке к

= 0. Поскольку закон дисперсии для

всех

 

о /с2

 

зоп но имеет

простого вида: h — ^ а " 2 т а * '» т 0 с а м и

ЭФ~

фектпвныо массы оказываются зависящими от волнового вектора. В слаболегированном арсениде галлия этот эффект (непараболичность зоны) невелик, тем не менее приведем законы дисперсии для зоны проводимости и зоны легких дырок, учитывающие непараболичность:

 

Bc.ift — - T ± y

4" +

- g ^ - •

 

 

При

этом эффективная масса

выражается

формулой:

 

 

ае

 

 

 

 

При сверхнизких

температурах,

порядка

1° К,

сле­

дует

учитывать вклад

линейного

по к

члена в

закон

дис­

персии. Тогда энергетическими поверхностями для тя­ желых дырок являются эллипсоиды вращения с осью

2.3]

П Р И М Е Н Е Н И Е

к р - М Е Т О Д А

85

вдоль {111} с

массами:

 

 

± =

l + M +

±-(L-M-N),

 

J - =

4_ + i _ M - i - ( L - M - y V ) .

 

Максимум зоны тяжелых дырок расположен на расстоянии

к т я ь

Ю - 3

К от центра зоны Бриллюэиа, где К — макси­

мальный

k-вектор в первой зоне Бриллюэиа, на высоте

em ^

Ю -

5 эв от потолка валентной зоны в толке к = 0.

В законе дисперсии для электронов вблизи дополнитель­

ных

минимумов в направлениях {100} и {111} из-за

низ­

кой симметрии уравнения Шредиигера для функции

Z7k(r)

вклад от линейных членов в спектр гораздо более сущест­ вен, чем в центре зоны Бриллюэиа. Если без учета спинорбитального взаимодействия экстремумы находятся на краю зоны Бриллюэиа, то с учетом спин-орбитальных эф­ фектов экстремумы смещаются. В окрестности минимума поверхностями постоянной энергии являются

е (к) = е 0 +

ру*? + p,fc? + atkt,

где at пропорционально

еппп-орбптальному взаимодей­

ствию, е0 — начало отсчета

энергии в точке к0 , места рас-

положенпя экстремума без учета

спин-орбитального

взаи­

модействия.

 

 

 

 

 

 

Экстремум имеет место вдали от точки к0 на расстоянии

&j=]^i/2p\|

с

минимумом,

расположенным

ниже

е 0 на

Emin =

/ 4|5/|.

Экстремумы

при этом

сливаются,

располагаясь

на

окружности

&( = const.

Такие

торои­

дальные поверхиости энергии впервые были найдены тео­ ретико-групповыми методами Рашбой и Щекой [16] в кри­ сталлах со структурой вюрцита и Шекой [3] в кристал­ лах со структурой сфалерита.

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

дующую

оценку для величины

сдвига минимума зоны:

em in Ао/9(е0 —е; ), где А0 — спин-орбитальное расщепление

валентной зоны в точке к = 0, а

е;- — энергия края ва­

лентной

зоны

в рассматриваемой

точке к0 .

Для

GaAs

эффект чрезвычайно мал, порядка 10~3 эв,

и может сказаться при очень низких температурах.

86

З О Н Н А Я С Т Р У К Т У Р А

[ГЛ. 2

2.4. Количественные расчеты зоипон структуры арсенида галлия

Количественному расчету зонной структуры арсенида галлия посвящено несколько работ [10—261. За исклю­ чением работ [10, 25], они ие претендуют на точные чис­ ленные результаты, а скорее дают правильную качествен­ ную картину расположения зон.

В работе Поллака и Кардоны [10], представляющей собой развитие пдей, изложенных в предыдущем разделе, расчет ведется kp-методом прп заданном значении энергий 15 состояний в точке к = 0 для Ge. Из них 4 состоя­ ния принадлежат валентной зоне и 11 зоне проводи­ мости.

Определитель 1'5-го порядка решается для нахождения спектра Ge, а для GaAs добавляется в качестве возмуще­ ния ант исимметричный потенциал, учитывающий различие атомов в решетке. Параметры теории паходятся из экс­ периментов по оптическим измерениям и циклотронному резонансу.

Отсутствующие данные для межзонных расстояний брались из результатов расчетов по методу ортогоналпзованных плоских волн. Кроме межзониых расстояний, необходимо было зиать значение 10 матричных элементов оператора импульса. Из них два были известны из экспе­ риментов по циклотронному резонансу, один определялся пз условия вырождения состояния Х^зоны проводимости в Ge, а один вычислен методом псевдопотенциала. Осталь­ ные 6 независимых параметров определяли, задаваясь энергиями состояпий Z 1 ( L 3 , L z - , зазора Х4 Хг , энер­ гией состояния Д : и расположением Аг минимума, из­ вестным для Si.

Вычисления для GaAs в точке к = 0 проводились ме­ тодом псевдопотенциала. Сущность метода псевдопотен­ циала состоит в том, что к потенциалу кристаллической решетки, создаваемому ионными остовами, добавляется некоторый потенциал, имеющий характер отталкивания электронов от ионов. Этот суммарный потенциал получил название псевдопотенциала. Сам псевдопотенциал может быть либо найден из экспериментальных характерис­ тик интересующего нас полупроводника или же из спект­ роскопических данных об уровнях энергии атомов,

2.4] К О Л И Ч Е С Т В Е Н Н Ы Е Р А С Ч Е Т Ы 87

составляющих кристаллическую решетку рассматриваемо­ го полупроводника.

Симметричная часть потенциала бралась такая же, как для Ge. Антисимметричная часть определялась шестью не­ зависимыми параметрами, один из которых оказалось воз­ можным положить равным нулю, остальные параметры были определены нз экспериментальных значений зазоров

Lsv — £ з с | r i 5 t ) — Г ) 5 С , Г 1 5 , , — T i c , Хъ-о ХИ Азтз — A i c .

Полученная картина зон приведена на рис. 2.4. Там же указаны'экспериментальные значения межзонных рассто­ яний, использованные в расчетах.

Рпс. 2.4.

Зопиая структура арсенида галлия, рассчитанная ]

 

kp-методом [Ю].

£fg , -

В работе Коэна'и Бергстрессера [25]

зонная структура

арсенида

галлия рассчитана методом

псевдопотенциала

в направлениях {100},

{111} и

{ПО}. Симметричная часть

псевдопотенциала

для

GaAs

бралась такая

же, как

и

для

Ge. Эта часть потенциала определялась

тремя

не­

зависимыми параметрами. Антисимметричная

часть так­

же

определялась

тремя другими параметрами. Мерой

88

З О Н Н А Я С Т Р У К Т У Р А

[ГЛ.

велпчнпы антисимметричного потенциала служило расщеп­ ление зоны X j - * (XjjX^) при переходе от Ge к GaAs.В валент­ ной зоне это расщепление оказалось порядка 4 эв. При определенпп параметров метода псевдопотенцпала спип-орбп- тальное взаимодействие не учитывалось. Рассчптаппая методом псевдопотеицпала зонная структура GaAs показана на рис. 2.5. Точность расчета—0,3 эв вблизи запрещенной зоны, и — 1 эв для остальных эпергпй.

Рис. 2.5. Зонная структура арсенпда галлия, рассчитанная мето­ дом псевдопотенцпала [25].

Караваев и Поплавпой [2G], пользуясь псевдопотепцпалом, найденным в работе [25], продолжили расчет этих авторов для точки W на границе зоны Брнллюэна.

Таким образом, общая картина зонной структуры ар­ сенида галлия подобна зонной стрза<туре германия.

Край валентной зоны находится в центре зоны Брнл­ люэна. Имеется одна зона тяжелых дырок и одна зона легких дырок, смыкающиеся в центре зоны Брнллюэна, и еще одна зона, отделенная от них из-за спин-орбиталь­ ного взаимодействия. Зоны тяжелых и легких дырок анизотропны, хотя и анизотропия последней слабая. Спин­ орбитально отщепленная зона изотропна. Все три зоны непараболичны [15]. Край валентной зоны слегка сдви­ нут от центра в направлении {111}.

Ыаинпзший

минимум зоны проводимости расположен

в цептре зоны

Бриллюэна и обладает изотропной эффек-

2.4]

К О Л И Ч Е С Т В Е Н Н Ы Е Р А С Ч Е Т Ы

89

тнвиой массой. Следующие минимумы зоны проводи­ мости расположены в направлениях {111}, па краю зопы Брпллюэиа и {100}. Изоэпергетическимп поверхностями вблизи этих минимумов являются эллипсопды вращения с осями вдоль {111} и (100} направлений соответственно. Учет спин-орбитальных эффектов приводит к тороидаль­ ным поверхностям постоянной энергии.

Теоретически

рассчитанные

межзоиные

расстояния

и эффективные

массы

электрона

и дырок в

точке к = 0

прнведепы в табл. 2.2.

Т а б л и ц а 2.2

 

 

 

 

М е ж з о ш ш е р а с с т о я н и я , эо

Э ф ф е к т и в н ы е м а с с ы

Г 1 8 у - Г 1 С = 1 , 4

[25]

 

 

' T J B V — Г 1 5 с = 4 , 5

[25]

 

 

Tjsv—^ic—1,8

[25]

mc=0,067

[9Д4]

 

 

 

[25]

/ns/i=0,135

[9,14]

^154—Ajc—1 >8

^i 5

v —

, 7

[25]

wz/i/i=0,45

[9,14]

^sv—^-ic=2,6

[25]

Ш1Л=0,091

[9,14]

^-3v—^зс —G,0

[25]

m(=0,51

[9,14]

•^5v—Xlc=4,0

[25]

m<=0,23

[9,14]

Л 1 с

^ з с -: 0,3

[25]

 

 

 

Д„=0,34

[12,13]

 

 

 

Ab=0,21

[12,13]

 

 

Более полные и точные количественные расчеты зон­ ной структуры арсенида галлия позволят надежнее вы­ числить параметры, определяющие эффективные массы не только в центре зоны Бриллюэна, но и в других точ­ ках, а также учесть релятивистские поправки к спектру.

Соседние файлы в папке книги из ГПНТБ