книги из ГПНТБ / Трусов Л.И. Островковые металлические пленки
.pdfзначениям в массивном металле. При этих допущениях
С\ ~ ІО16-г ІО1 8 дин-1сек~\
Энергетический барьер, соответствующий образова нию критического зародыша, определяется по формуле (3), конкретизированной для случая роста из пара и ку полообразной формы зародыша (см. например, [1 0 ]), причем для оценки пересыщения используются экспери ментальные зависимости упругости паров р от температу
ры для пленок висмута и серебра. |
Подстановка числен |
||||
ных значений параметров в формулу |
(54) приводит к |
||||
выражению, в котором зависимость скорости |
зарожде |
||||
ния от температуры |
при |
данной |
плотности |
атомного |
|
пучка определяется |
единственным |
параметром L = |
|||
= 2Едес—Едиф. Типичные |
зависимости |
log/c—Т, |
вычис |
ленные на ЭВМ при различных значениях L, показаны на рис. 4. Сравнение с экспериментальными кривыми по
Рис. 4. Сравнение температурных зависимостей скоростей зарождения (экс периментальных и рассчитанных на основании термодинамической теории) для различных значений С =2Ядес—£ диф, ккал/моль:
а — 40 (/); 30 (2); 20 (3); б — 30 |
(/); 24 |
(2); 20 (3); 16 |
(4); 10 |
(5). Сплошные |
линии — теоретические кривые; |
темные |
кружки — экспериментальные данные |
||
для висмута на подложке С; светлые |
кружки — для |
висмута |
на подложке |
|
SiO [30] |
|
|
|
|
называет, что согласие с опытом неудовлетворительное. Как отмечает Поппа [30], наблюдаемое несоответствие можно уменьшить, если учесть поправку к поверхност ному натяжению, связанную с высокой кривизной по верхности малых зародышей. Кривые log/c—Т, раесчи-
40
тайные по формуле |
(54), в которой значение а заменено |
|||||||
на а/4, |
показаны |
на рис. 4, |
б. |
Если |
Ет$ < £ Дее, |
то |
||
Едеc~L/2, и из кривых на рис. 4 видно, что для |
пленок |
|||||||
висмута |
на |
подложке С |
величина |
£дес = |
12-г |
14 |
||
ккал/моль, что |
приблизительно |
в 1,7 раза выше зна |
чения, полученного из сравнения со статистической тео рией. В работе [30] показано, что для пленок серебра на
подложке С расхождение |
несколько меньше, однако |
для пленок висмута на |
подложке SiO значения Едес |
различаются почти в два раза. Расхождение результатов термодинамического расчета е экспериментом, а также с зависимостями, построенными на основе статистической
теории, существенно |
увеличивается в области |
низких |
||
температур подложки. |
Это |
неудивительно, |
поскольку |
|
оценка по формуле типа (2 ) |
показывает, что |
іс = |
2 + 3 |
|
(для висмута) и іс— 1/2-Ь2 |
(для серебра). |
Очевидно, |
что для критических зародышей такой величины термо динамический анализ не является корректным.
Другой характеристикой, для которой возможно сра внение с экспериментом, является плотность зародышей при насыщении. Анализируя зависимости, подобные за висимостям, показанным на рис. 1 и 2, Поппа заключил,
что имеет место соотношение |
|
|
N ' ^ N i |
|
(55) |
Используя формулу (45), можно получить: |
|
|
log Ns — A (ic, а...) 4 іе log 1k |
*£ £дес A E[ |
1000 |
+ |
.(56) |
|
|
4Д 8 |
|
Экспериментальные зависимости для висмута на под ложке С показаны на рис. 5. По наклону низко- (2) и
Рис. 5. Температурная зависимость логарифма плотности зародышей при насыщении при конденсации висмута на графите [30] (/ =1,6*
• 10м см— 2 сек— 1 )
41
высокотемпертурной (1) ветви графика определяли величину ісЕЛы-\-Еіс, которая оказалась равной 10,5
и 52 ккал/моль, |
соответственно. |
Отсюда |
получаются |
|||
разумные |
значения £ дес и Е3, если, как и прежде |
счи |
||||
тать, что |
для |
высокотемпературной ветви |
(1 ) |
іс= 3, а |
||
для низкотемпературной ветви (2 ) іс= 1 . |
|
что |
фор |
|||
В работе Льюиса [96, с. 30] |
указывается, |
мула (82), широко используемая в теоретических рабо тах, противоречит экспериментальным результатам Поппа. Действительно, согласно формуле (82), плот ность стабильных зародышей на поверхности при насы щении iVs не зависит от плотности потока Д и возрастает с повышением температуры. Практически [30] Ns уменьшается с повышением температуры и уменьшени
ем Д. Льюис отмечает также, |
что |
хотя наблюдалась |
||||
пропорциональность между іѴ5 |
и |
N ic , |
но по абсолютной |
|||
величине они различаются |
в 1 |
0 |
5 |
раз |
и объяснить |
этот |
факт до -сих пор не удалось. |
|
|
|
|
|
|
Более строгое рассмотрение процесса захвата адато |
||||||
мов зародышами, следствием которого является |
насы |
|||||
щение зависимости ns(x), |
провел |
Гальперн [61]. Он |
рассмотрел случай, когда длина поверхностной диффу зии адатомов ds сравнима с диаметром островка Do.
Вэтом случае*
*(D, + 1,2 ds)2
Однако оказалось, что незначительный рост D0 с уве личением температуры не может компенсировать умень шение dg, связанное с реиспарением.
Ротлидж и Стоуэл [91], используя метод ячеек (см. гл. II, п. 7), рассмотрели вопрос о распределении оди ночных адатомов на подложке в промежутках между островками. В отличие от Гальперна, они учли, что об ласти захвата для различных островков могут перекры ваться, а также то, что часть адатомов может реиспариться до того, как они будут захвачены островками. Соответствующее уравнение для П\(г) имеет вид
d2n1 |
1 |
drii |
ч tii |
Ik |
(58) |
|
dr2 + |
7 |
Tr |
ТГ + “ГГ = |
|||
|
Граничные условия отражают то обстоятельство, что островки являются идеальными стоками для адатомов
42
[п1 ( ^ ) = 0 ] |
и то, что границы ячеек |
(площадь которых |
|||||||||
лL2= n s~l) |
непроницаемы |
d г = |
О, |
при r = L ) . |
|
||||||
Решением уравнения (58) является функция |
|
|
|||||||||
|
|
|
/ о |
(kr) |
К 1 (Щ + |
h |
(Щ |
Ко (kr) |
' |
(59) |
|
|
|
|
Io (kR) Ki (Щ + h |
(Щ |
Ko |
(kR) |
’ |
||||
|
|
|
|
||||||||
где Іі |
и Ki — модифицированные |
функции |
Бесселя |
по |
|||||||
рядка |
і и k2=v/D. При L -* оо |
из |
выражения |
(59) |
по |
||||||
лучается формула Гальперна [61]. |
-* ® |
) из выраже |
|||||||||
В случае полной конденсации |
(т8 |
||||||||||
ния |
(59) следует, что |
|
|
|
|
|
|
|
|
||
" . м |
- |
Ч т г 111 ( т ) - |
-Й > |
|
|
|
|
|
<6 °) |
||
Как указал Марков [97; 98, с. 119], это |
решение не |
||||||||||
учитывает тот факт, что в результате поступления |
ато |
мов из пара плотность одиночных адатомов п\ линейно зависит от времени и, следовательно, необходимо рас сматривать нестационарное уравнение диффузии. Реше ние, предложенное Ротлиджем и Стоузлом, справедливо в случае когда характерное время поступления атомов из
пучка N0/Ik |
больше |
характерного |
времени |
диффу |
зии L2jD. |
анализ, |
проведенный |
в работе |
[91], и |
Численный |
сравнение с результатами Гальперна и Льюиса, показали, что учет перекрытия областей захвата адатомов различ ными растущими островками приводит к существенному уменьшению плотности адатомов по сравнению со ста ционарным значением, при котором nl= Ikxs.
Приняв во внимание эти факторы, Логан [99] и Стоуэл [100] вывели более строгую зависимость плот
ности |
островков |
при |
насыщении от параметров |
роста, |
|
в частности, от размера и энергии образования |
крити |
||||
ческого зародыша, а также от плотности потока: |
|
||||
|
|
і+ 1 |
|
|
|
N, |
к і+ъ |
|
Е, |
(61) |
|
D |
exp |
(I + |
3) k Т |
||
|
|
|
Применив для численных расчетов ЭВМ, Ротлидж и Стоуэл [91] получили зависимости плотности стабиль ных зародышей на поверхности при насыщении Ns от интенсивности потока частиц при испарении Іи, темпера
43
туры Т, энергии активации адсорбции £ адс, диффузион ной энергии £диф и энергии диссоциации E t в виде, удоб
ном для сравнения с экспериментом. Типичные |
зависи- |
|||||||
|
Е |
|
|
|
|
|
|
|
мости logyVjs—logІи и Ns ----- показаны на рис. 6 |
|
|||||||
и 7. В частности, на рис. 6 видно, что |
плотность стабиль |
|||||||
ных зародышей Ns возрастает при увеличении |
плотно- |
|||||||
Рис. 6. Логарифмическая зависимость плотно |
||||||||
сти стабильных |
зародышей |
при |
насыщении |
|||||
от |
интенсивности потока |
частиц |
1 ^ для |
|||||
различных параметров роста: |
|
|
|
|||||
|
|
аде |
■- |
15,5; |
диф |
|
|
|
|
|
kT |
|
|
kT |
|
|
|
|
Е |
аде |
- |
15,5; |
Е |
4. |
|
3,9; |
= 2- |
|
Диф - - |
||||||
|
k Т |
|
|
k T |
|
|
||
Ег _ |
15,5; |
3 —/ = 3 ; |
Еаде =16,9; |
|||||
k Т |
|
|
|
|
k Т |
|
|
|
£Диф |
■—4,2; |
Ез —16,9 |
|
|
|
|
||
k т |
|
k T |
|
|
|
|
|
|
сти потока. Это качественно согласуется |
с эксперимен |
|||||||
тальными наблюдениями Поппа [30]. |
|
|
|
|
|
Анализируя зависимости Ns от параметров роста, ав торы работы [91] пришли к заключению, что следует проявлять осторожность при определении величины іс на основании измерений плотности островков при насы-
Рис. |
7. Зависимость N |
от |
Диф |
|
А Т |
||||
при |
S |
|
||
различных размерах |
|
критиче |
ского зародыша по данным работы
[911: |
= 1; |
2- |
іс =2; |
3 - |
= |
|
1 - |
і с |
|||||
= 3 |
(/к |
=10" |
„ -2, |
- 1 |
|
|
см |
* сек |
|
|
|||
Ег=4Е Диф’ Е3= |
10 Е ДИф' |
|
щении. Например, при типичных значениях параметров
(4 = Ю 1 3 см~2-сек-\ — 5, =20) неста
ционарная плотность при насыщении достигает мак симума за время порядка ІО2 —ІО3 сек (рис. 8 ). В тече-
44
ние этого времени условия Hâ поверхности существенно изменяются и определяемый по положению максималь ного значения величины Ns размер критического зароды ша не обязательно должен совпадать с величиной іс на начальной стадии зарождения. Таким образом, если кри тический зародыш содержит более, чем один атом, опре деляемое этим методом значение іс является некоторым эффективным параметром.
Рис. 8. Зависимость плотности стабильных зародышей/! от времени конден
аде
сации по данным даботы [91], при различных отношениях
к Т
1—16,9; 2—18,6; 3—20,7; 4—23,0; 5—25,3; 5—27,6
/ , = 1013 |
сек −1. |
=3; |
диф = 5; Ез = 20 |
к |
|
|
k Т |
Льюис считает, что более удовлетворительное объяс нение зависимостей плотностей стабильных зародышей от времени и температуры п3(т) и NS(T) может быть дано на основе гипотезы о миграции многоатомных островков. Рассуждая так же, как и при выводе фор мулы (82), он получает выражение для скорости взаим ного захвата островков, содержащих по s атомов и миг рирующих с относительной скоростью щ:
/ = |
«s üd |
(62) |
Nо
Если воспользоваться экспериментальными значе ниями параметров для пленок висмута на подложке С, то получится, что
45
|
« 1 ,5 9 -1 0 9 |
CM~2 сек~\ |
|
|
N0 |
|
|
Это |
означает, |
что |
миграционная подвижность ост- |
ровков |
должна |
|
О |
составлять приблизительно 30 А/сек. |
Льюис считает это значение разумным, поскольку Ет ф для отдельного адатома висмута на подложке С равна всего лишь 3 ккал/моль. Если эта гипотеза верна, то уменьшение Ns с увеличением температуры обусловлено тем, что активизируется миграция и последующее объе динение малых островков на стадии зарождения [96, с. 30].
Помимо этого, Льюис показал, что согласие резуль татов Поппа со статистической теорией становится превосходным, если предположить, что частота колеба ний адатома на поверхности равна 5 -10й сект1 (вместо традиционного значения ІО1 3 сек-1). В случае атомов висмута, связанных с нейтральной подложкой С слабы
ми силами физической адсорбции, |
это |
допущение |
не |
|
кажется неправдоподобным [1 0 1 ]. |
|
|
|
|
В настоящее время имеется ряд работ, в которых |
||||
приведены экспериментальные доказательства |
справед |
|||
ливости как термодинамической [10; 18, |
с. 15], |
так |
и |
|
статистической теории [85; 8 6 ; 87, |
с. 3], |
причем работы |
первой группы выполнены в условиях низких, а второй — относительно более высоких пересыщений. Однако очень мало число работ, в которых определенная система ис
следуется в широком диапазоне значений |
пересыщения |
|||
и экспериментально |
подтверждаются |
преимущества |
||
каждой теории. |
Здесь, |
помимо работ Поппа, |
следует |
|
упомянуть работу |
Корбетта и Боусвелла |
[102, |
с. 2663], |
в которой исследовался процесс зарождения пленок се
ребра на молибдените MoS2. |
Плотность атомного пучка |
|||||
варьировали |
в интервале 8 |
-101 1 -У- 7 -ІО1 3 |
см~2 сект1, а |
|||
температуру |
подложки изменяли |
от 30° до 200°С. |
При |
|||
высоком пересыщении |
(/& = 7- ІО1 3 |
смг2 сек-1, 7’ = 30°С) |
||||
радиус критического |
зародыша |
должен |
быть |
равен |
~ 1 0 А. Если воспользоваться численными значениями параметров этого эксперимента и провести оценки на основе термодинамической теории, то получатся ненеправдоподобно низкие значения скорости зарождения. Помимо этого, как и в работе [30], не удается описать
46
зависимость /ДГ), используя какие-либо определенные значения поверхностного натяжения а и контактного угла Ѳ. Однако, по мнению Корбетта и Боусвелла, в этом случае отсутствует согласие и с теорией Уолтона. Дей ствительно, наблюдавшуюся линейную зависимость /с от Іи, можно объяснить, если критический зародыш со держит всего один атом. Другими словами, необходимо допустить, что в этом эксперименте наблюдалось заро ждение на точечных дефектах, причем плотность послед них должна быть порядка ІО1 1 см~2. Но Корбетт и Боусвелл считают, что подложки молибденита были зна чительно более совершенны. Более определенно о раз мере критического зародыша можно судить по зависи мости log Іи—1IT, однако в работе [102, с. 2663] эти зависимости не были получены.
Корбетт и Боусвелл полагают, что количественное со поставление с экспериментом возможно лишь после, то го, как будут строго учтены объемные и поверхностные статистические суммы, а также члены в выражении для свободной энергии, учитывающие энтропийные ограниче ния. Роль последних, как известно, особенно велика при
гомогенном зарождении |
[18,. с. 15; |
103], а также в не |
||
которых случаях при гетерогенном |
зарождении (ртуть |
|||
на кварцевом стекле [104], кадмий |
на |
меди |
[ГО, 105], |
|
натрий на подложке CsCl [106]). |
Однако, как показал |
|||
Поппа [30] при высоких |
пересыщениях |
эти |
поправки |
|
мало существенны. |
|
|
|
|
Специальные исследования влияния точечных дефек тов на процесс образования зародышей провел В. С. Пост ников с сотрудниками [107, е. 87; 108; 109]. Для этого скол в вакууме кристалла NaCl, на который конденси ровали золото, предварительно облучали у-излучением в течение 20 мин. Одновременно конденсацию проводили и на необлученный скол. По электронномикроскопичес ким снимкам были построены функции распределения островков по размерам, которые имели характерный вид кривых с максимумом. Сравнение функций распределе ния показало, что плотность островков на облученной подложке приблизительно в 4 раза выше, чем на необлу ченной. Увеличение плотности сопровождается относи тельным уменьшением дисперсии функции распределе ния и смещением ее максимума в область меньших размеров. Наблюдавшееся увеличение плотности Пост
47
ников и другие исследователи связывают е возрастанием поверхностной концентрации точечных дефектов под влиянием радиации. Уменьшение среднего размера мо жет быть обусловлено заряженностью образующихся дефектов. Действительно, с вакансиями на поверхности щелочногаллоидного диэлектрика могут быть связаны ло кальные заряженные микрообласти [НО]. При этом ло кальное пересыщение при зарождении в такой области может быть существенно выше среднего пересыщения над поверхностью подложки [1 1 1 , 1 1 2 ].
Очевидно, что при прочих равных условиях увеличе ние плотности островков должно быть более существен ным для подложек с меньшей энергией образования дефектов. Этот эффект действительно наблюдался [109], когда в качестве подложек использовали ряд щелочногаллоидных кристаллов (NaCl, KCl, KBr).
Оказалось, что условия зарождения на ступенях скола после облучения изменялись незначительно. Это может быть связано с относительно малым вкладом зарядов вследствие высокой поверхностной энергии ступени. Кроме того, мала вероятность образования де фектов точно на ступени из-за малого сечения рассея ния у-квантов.
В работах В. С. Постникова и др. было также пока зано, что плотность островков зависит не только от со стояния поверхности подложки, но и от состояния кон денсируемых атомов металла. Так, ионизация атомов золота в пучке, осаждаемом на кристалле NaCl элект ронной бомбардировкой, приводит к такому же увеличе нию плотности островков, как и повышение дефектности подложки. Эти эксперименты показали, что ионизация атомов потока способствует совершенству ориентации, уменьшает плотность дефектов в пленке и приводит к об разованию сплошной пленки на более ранних стадиях. Кроме того, они (качественно согласуются с предполо жением, (высказанным в работе [ 1 0 2 . с. 2663], относитель но роли точечных дефектов. Однако, поскольку в (работе [109] не приведены кинетические зависимости скорости образования зародышей, то невозможно количественно оценить размер критического зародыша при росте на де фектной поверхности.
Таким образом, в настоящее время корректность рас смотренных теорий, имеющих и внутренние ограничения,
48
трудно оценить. Вместе с тем,. требуется дальнейшее совершенствование современной техники эксперимента, чтобы получить однозначные результаты о механизме зарождения пленок.
Как известно, одной из важнейших задач статисти ческой теории кристаллизации является «расшифровка» смысла параметров, используемых в феноменологичес ких теориях роста. Это положение можно проиллюстри ровать на примере работ Д. Е. Темкина [45 с., 63, т. 5, с. 89; 44, с. 15], Джексона и других авторов [31, 36], в которых методами статистической механики рассматри валась теория движения фазовой границы при росте пленки (ср. с п. 2). Джексон [31] в приближении моле кулярного поля рассчитал изменение свободной энергии первоначальной плоской фазовой границы при присоеди нении к ней новых молекул. Учитывая взаимодействие ближайших соседей и предполагая, что новые частицы располагаются на поверхности подложки неупорядочен но, он получил формулу для изменения свободной энер гии при увеличении доли ц занятых поверхностных узлов:
- A L - = |
ССі 7) |
(1 — 7)) + 7J ln Y| + (1 — 7]) ln (1 — 7)), |
(63) |
N k T e |
|
|
|
где |
|
|
|
LT |
|
NA |
|
0l ~ ~RT~e |
У |
И ^ — ~No~ ’ |
|
LT — скрытая теплота превращения; Te — равновесная температура;
NQ— общее число поверхностных узлов на первоначаль
но плоской грани подложки при равновесной |
тем |
пературе; |
по |
NA — число частиц, присоединившихся к плоской |
|
верхности; |
|
q— часть общего числа ближайших соседей, располо женных во вновь образующемся слое.
Зависимости, полученные из формулы (63), показа ны на рис. 9. Оказывается, что при а і< 2 минимум сво бодной энергии соответствует поверхностям, у которых сорбированными частицами занята в среднем полови на общего числа узлов: такие поверхности шерохо ваты. При а і> 2 минимумы свободной энергии соответ
49