Научная статья на тему 'Учет свободной поверхностной энергии при расчете диаграммы состояния многокомпонентных твердых растворов А3В5, полученных методом эпитаксии из жидкой фазы на подложке фосфида Индия'

Учет свободной поверхностной энергии при расчете диаграммы состояния многокомпонентных твердых растворов А3В5, полученных методом эпитаксии из жидкой фазы на подложке фосфида Индия Текст научной статьи по специальности «Химические технологии»

CC BY
65
14
i Надоели баннеры? Вы всегда можете отключить рекламу.

Аннотация научной статьи по химическим технологиям, автор научной работы — Ратушный В. И., Шишков M. B., Уелин В. В.

Исследовано влияние кристаллографической ориентации подложки InP на состав твердого раствора GaxIn1-xPyAs1-y, полученного методом жидкостной гетероэпитаксии. Представлен метод расчета состава с учетом поверхностной энергии. Полученные результаты показывают величину вклада межфазной поверхностной энергии границы раздела кристалл-расплав в полную избыточную энергию смешения твердого раствора.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Учет свободной поверхностной энергии при расчете диаграммы состояния многокомпонентных твердых растворов А3В5, полученных методом эпитаксии из жидкой фазы на подложке фосфида Индия»

ПОЛУПРОВОДНИКОВЫЕ МАТЕРИАЛЫ И ЭЛЕМЕНТЫ

УДК 621.315.592:548.55

УЧЕТ СВОБОДНОЙ ПОВЕРХНОСТНОЙ ЭНЕРГИИ ПРИ РАСЧЕТЕ ДИАГРАММЫ СОСТОЯНИЯ МНОГОКОМПОНЕНТНЫХ ТВЕРДЫХ РАСТВОРОВ А3В5, ПОЛУЧЕННЫХ МЕТОДОМ ЭПИТАКСИИ ИЗ ЖИДКОЙ ФАЗЫ НА ПОДЛОЖКЕ ФОСФИДА ИНДИЯ (111)

© 2003 г. В.И. Ратушный, M.B. Шишков, В.В. Уелин

При росте эпитаксиальных пленок существует зависимость состава осаждаемого твердого раствора от кристаллографической ориентации подложки ((111) A

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

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

Gm _ Gm

+ Gml + G

ef

dn

i 1

Для четырехкомпонентной системы АхБ1-хСуД-у

уравнения, связывающие ликвидус и солидус, имеют вид

ЯТ Щу?, • х?,) = ЯТ 1п[4хгЬхЬ уЬ у Ь / у ?Ь у?Ь ] +

+ ASF (TF - T

нений в твердой фазе; х?, - мольная доля компонента

в твердом растворе; хЬ - концентрации компонентов

в жидкой фазе; у^ и - коэффициенты активности

компонентов соответственно в жидкой фазе и в переохлажденной жидкости стехиометрического состава; А?г^ - энтальпия плавления; АТгр, - температура

плавления.

Если рассматривать поверхностную энергию как меру недостатка энергии сцепления на мономолекулярной поверхности, то она будет равна [2]

Gef _

f _ a* AS (1 -ß)AH

Sa2 NA

(1)

где - избыточная энергия смешения ненапря-

женной твердой фазы, О^1 - упругая энергия деформированного эпитаксиального слоя, О^ - межфазная поверхностная энергия границы раздела «кристалл - расплав».

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

С практической точки зрения наиболее приемлемо квазихимическое приближение парного взаимодействия для расчета термодинамических характеристик твердых растворов А3В5. Согласно [1], коэффициенты активности бинарных компонентов у?, определяются выражением

ЯТ 1п у?, =3^ (Отп).

где а* = 1 / Ь , в свою очередь Ь следует рассматривать как количество уровней решетки в эпитаксиаль-ном слое (для данной работы а = 0,05); в - отношение количества ближайших соседей атома, находящегося на поверхности, к количеству соседей атома, находящегося внутри кристалла (в = 0,8 для (111)А и в = 0,6 для (111)Б [3]); А? - площадь поверхности, для ориентации (111)

AS _тт S 2

1 Sna

4

(2)

Подставляя (2) в (1), получим:

Gef _

f _ a(1 - ß)Aff

8

Выражение для энтальпии находится следующим образом:

lS £iiJ _ ABCD ABCD '

S

m ABCD

HABCD _ *HmSABCD + ^fc + (1 - х)УН'вС +

+ x(1 - y)hasd + (1 - x)(1 - y)hbD ;

mmSABCD _ ^ (x(1 - х)УаAC-BC + x(1 - x)(1 - У)аAD-BD +

где у- коэффициенты активности бинарных соеди-

+ xy(1 - y)aSAC-AD + (1 - x)y(1 - y)a

BC-BD

);

2

a

HAC - -4 а AC + "2 HA + "2 HC - TAC ASAC

hBc - 1 aBc + "2 HB + 1 Hc - TBC asbc ;

H AD - 1 aAD + "2 HAl + "2 HD - TAD ASAD

hbd - 1 aBD + "2 HBl + 2 HD' - TBD asBD ■

H'abcd - ^HWabcd + xAHAl + ^BHBl +

RT In yAD - [yaAc_BC + (1 - У)аAD_BD](1 - X)2 +

+ [xaSAC_AD + (1 - x)aSBC_BD ](1 - У)2 + + [aC + x(aAD BD - aAC _ BC ) —

- (1 - y)(aBc BD -a SAC_AD )](1 - x) У + ((1 - x)ya AC BC +

a(1 - ß) ( 1

8 v 2

+ (1 - x)asAD

BD + xya AC AD ) -

+ XlGH'Cl + xDHD1 ;

АНшЛВСЭ - ХЛХБ аЛв + ХЛХС аАС + ХБХС аВС +

_1_ _1_ VI VI п/1 _1_ VI п/1

+ Х ЛХ ОаАО + Х БХ Э аБЭ + Х СХ ОаСО •

Получаем конечную формулу для свободной поверхностной энергии

а" -^^СХЛХБаЛБ + ХЛХСаЛС +

VI VI /у/ 4- у' у' Г»' 4- у' у' (V 4- у' У"' (V 4-ХБХСавс + ХЛХЭаЛЭ + хбхэабэ + х сх эасэ +

+ xAHA1 + XBHB1 + XCHC + XDHD1 -

- (4 aAD - TFD ■ &SFD ) -1 (HA + HD)):

RT ln yBc - [yaSAC BC + (1 - y)aSAD BD]x2 +

+ [xaLC_LD + (1 - x)aSBC_BD](1 - У)2 +

+ [aC - (1 - x)(aAD BD - aAC BC ) +

+ y(aBC BD - aLC _ AD )]X(1 - У) +

+ a(18- ß) (- 2 (xaSAC BC + x(1 - y)aSAD BD +

+ (1 - x)(1 - y)aSbc bd ) -

- 2 (x(1 - x) ya Ac-bc + x(1 - x)(1 - y )a Ad - bd + + xy(1 - y)a Ac - ad + (1 - x) У(1 - y)a Bc - bd ) -

- xy(1 aAc - TFcASFc) - (1 - x)y(-r aBc - TgcASfc) -

- x(1 - y)(- alAD - TFdASFd ) - (1 - x)(1 - y) X

х (4 а Во - Т^э ) - - хНЛ1 - - уН'С -

- 1 (1 - х)НБ - 2 (1 - У)НЭ).

Таким образом, коэффициенты активности с учетом поверхностной энергии имеют вид:

Ш 1п уЛс - [УаЛс _ВС + (1 - У)аЛЭ _БЭ](1 - х)2 + + [хаЛс_лэ + (1 - х)аБс_БЭ](1 - У)2 -

- [аС + х(алэ БЭ - аЛс _ ВС ) +

+ У(абс_бэ -алс_лэ )](1 - х)(1 - У) +

, а(1 -в)( 1 ((1 х) Уа 5 +

+ -8-(- 2 ((1 - Х)уаЛС _ БС +

+ (1 - У)аЛс ЛЭ + (1 - х)(1 - У)абс БЭ ) -

- (4aac - tfc asac ) - j(hla + h^));

- (1 a bc - Tbfc ASfc ) - 2 (hB + hlc));

rt ln Ybd - [yaslc_bc + (1 - y)aLD_BD ]x2 +

+ [xaLC _ LD + (1 - x)abc _ BD ]У 2 -- [aC - (1 - x)(aLD _ BD - aLC _ BC ) -- (1 - y)(asbc _ BD - aLC _ LD )]xy + + a(1 ~ ß) (- .1 (x(1 - y)asld bd + xyaslc LD +

82 1

v4

+ ya bc bd ) - (1 a ld - tfd asfd ) - 2 (ha + hd)).

В работе были теоретически определены фазовые равновесия системы Оа!пРЛ&' на подложке 1пР при температуре эпитаксии Т - 923 К . Условие изопе-риодичности для данной системы имеет вид у - (0,1933 - 0,431 х) /(0,1933 + 0,0115х).

Термодинамические параметры, использованные при расчетах, приведены в таблице.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Результаты расчетов представлены на рис. 1-3.

Как видно из графиков, наиболее сильное влияние поверхностного фактора испытывает концентрация ва в твердом растворе и в гораздо меньшей степени по Р и лб. Различия в составах жидких фаз равновесных твердому раствору одного состава для ориентаций (111)А и (111)В могут достигать по ва 20 - 25 %. А для лб и Р этот же показатель лежит в пределах 4 - 5 %.

Термодинамические параметры

Т^ = 1740 (К) [1] АИ1ва = -270061,44 (Дж / моль) [2]

Т[пр = 1335 (К) [1] АН'Пп = -232023,44 (Дж / моль) [2]

Т= 1511 (К) [1] АНР = -207261,12 (Дж / моль) [2]

ТПА = 1215 (К) [1] АН=-305833,88 (Дж / моль) [2]

А?ППаР = 72,4 (Дж / моль ■ К) [1] аОаР = 8876 - 18,63 ■ Т (Дж / моль) [1]

А?Па = 69,7 (Дж / моль ■ К) [1] аьаА, = 21604 - 38,35 ■ Т (Дж / моль) [1]

А?ПпР = 63,64 (Дж / моль ■ К) [1] аЬпР = 14980 - 14,82 ■ Т (Дж / моль) [1]

А?^ = 60,80 (Дж / моль ■ К) [1] аЬА = 16161 - 41,87 ■ Т (Дж / моль) [1]

аОаР № = 14650 (Дж / моль) [1] аОаР _ ОА = 1670 (Дж 1 моль) [1]

аОА _ ЪАЪ = 12560 (Дж 1 моль) [1] а?пр _ ЬА = 1670 (Дж / моль) [1]

аЬа1п = 4380 (Дж / моль) [1] арА^ = 6280 (Дж / моль) [1]

0.36 мол.доли

Рис. 1. Зависимость содержания галлия в жидкой фазе от мольной доли х твердого раствора Оах1п1-хРуА^1-у для

Т = 923К;- равновесный случай;-----плоскость (111)А;

---- плоскость(111)В

0.36 Люл.Ьопм

Рис. 2. Зависимость равновесного содержания фосфора в жидкой фазе от мольной доли х твердого раствора Оах1п1-хРуАя1-у для Т = 923 К

Рис. 3. Зависимость содержания мышьяка в жидкой фазе от мольной доли х твердого раствора Оах1п1-хРуАз1-у

для Т = 923 К

Следует отметить, что в данной системе именно галлий является элементом с наиболее высоким коэффициентом сегрегации. То есть при синтезе твердого раствора заданного состава необходимо уделять повышенное внимание содержанию галлия в жидкой фазе. Следовательно, для корректного расчета параметров технологического процесса получения изопе-риодных эпитаксиальных слоев твердого раствора Оа1пРА^' учет поверхностной составляющей полной энергии системы является необходимым условием.

Выводы

Определено влияние межфазной поверхностной энергии на состав твердого раствора ОаЫРАз при жидкостной эпитаксии на подложке 1пР с кристаллографической ориентацией (111)А и (111) Б . Избыточная

составляющая межфазной поверхностной энергии возрастает в ряду (111)А - (111)Б , что адекватно

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

Волгодонский институт Южно-Российского государственного технического университета (НПИ)

Литература

1. Кузнецов В.В., Москвин П.П., Сорокин В.С. Неравновесные явления при жидкостной гетероэпитаксии полупроводниковых твердых растворов. М., 1991.

2. Nakajima K., Ujihara T. et al. //J. of Crystal Growth 220. 2000.

3. Ормонт Б.Ф. Введение в физическую химию и кристаллохимию полупроводников. М., 1982.

18 октября 2002 г.

УДК 621.315.592:548.55

СВЧ ДИОДЫ С БАРЬЕРОМ ШОТТКИ НА НАНОПОРИСТЫХ ПОДЛОЖКАХ

© 2003 г. И.Н. Арсентьев, Л.С. Вавилова, Л.С. Лунин, В.И. Ратушный, M.B. Шишков

Использование тонкопленочных гомо и гетеро-эпитаксиальных структур на основе соединений А3В5 в электронных приборах различных типов предъявляет повышенные требования к их структурному совершенству и стабильности. Такие требования особенно существенны для материалов, используемых при создании СВЧ-приборов, что выдвигает в число фундаментальных проблем разработку новых методов получения эпитаксиальных структур и приборов на их основе. Одним из таких методов является эпитакси-альное выращивание на «мягких» пористых подложках полупроводников А3В5. Наличие пористого подслоя позволяет существенно снизить внутренние механические напряжения в получаемых структурах и повысить качество получаемых слоев.

Для исследования использовалась двухслойная гомоэпитаксиальная структура 1пР (п++-п+-п), барьер Шоттки формировался титановой металлизацией с последующим золочением, омические контакты были сделаны на основе золотогерманиевой эвтектики.

Пористые слои в кристаллах 1пР (100) п-типа проводимости (п~1018 см-3) были получены в анодном электрохимическом процессе в режиме импульсной подачи напряжения на ячейку с частотой 2 Гц.

В качестве электролитов использовались кислые (рН=0,5^1,5) водные растворы фторидов (КН4Р+ИР), хлоридов (КН4С1+ИС1) и бромидов (КВг+НВг). Напряжение на ячейке в импульсе в несколько раз превышало пороговые напряжения начала порообразования и составляло 7^10 В.

Объем пор заполнен продуктами неполного окисления 1пР, возникшими в результате его взаимодействия с анионными компонентами растворов по механизму реакции нуклеофильного замещения (исходный состав этих продуктов - 1пРХ2+п; п>2;Х - Б, С1, Вг). В процессах порообразования продукты реакции выполняют функцию твердого электролита, обеспечивающего транспорт адсорбирующих анионов к фронту реакции на дне поры. Впоследствии не исключено разложение остающихся в порах продуктов с выделе-

нием элементарного фосфора и тригалогенидов, а возможно и более сложных смешанных галогенидов, сохраняющих связи 1п - Р. На рис. 1 представлен скол подложки со сформированной пористой структурой. Хорошо различимы нитевидные поры, заполненные продуктами взаимодействия с электролитом. Глубина пористости в зависимости от состава электролита и плотности тока имела толщину в пределах 5-20 мкм и явно выраженный градиент плотности пор, направленный из объема кристалла к поверхности. Для эпи-таксиальных процессов использовали подложки с глубиной пористости 10-12 мкм.

Рис. 1. Граница раздела: пористая структура - объем кристалла

На рис. 2 показана поверхность пористой подложки. Видно, что исходная поверхность сингулярной грани практически утрачена, и произошло замещение ее системой вицинальных холмиков.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.