Файл: Лодиз, Р. Рост монокристаллов.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 14.10.2024

Просмотров: 142

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

472

Р. ПАРКЕР. МЕХАНИЗМЫ РОСТА КРИСТАЛЛОВ

параметров), очевидно: ступени, сильно удаленные друг от друга, будут развивать максимальную скорость; по мере сбли­ жения ступеней их диффузионные поля перекрываются и движе­ ние будет замедляться, как это и наблюдается при росте по спиральному механизму. Продифференцировав наклон (20.1) по t, а скорость роста (20.2) по х, получим с учетом условия непрерывности (в данном случае закона сохранения числа сту­ пеней), что

 

 

 

 

 

 

 

(20.4)

Определив

скорость

кинематической волны

как c(k)

=

dq/dk,

найдем из соотношения

(20.4)

следующее нелинейное

дифферен­

циальное уравнение:

 

 

 

 

 

у

 

 

 

 

^ ) | г

+ Ж = 0 -

 

 

(20-5)

Из уравнения (20.5)

следует, что в плоскости xt вдоль линии

с наклоном

dx/dt =

c(k)

= dq/dk плотность

ступеней

k

(и,

сле­

довательно, q) постоянна. Такая линия называется характери­ стикой. Кинематическая волна—это участок кристаллической поверхности с постоянной плотностью ступеней k (и постоянным потоком q), который перемещается со скоростью c(k) = dx/dt. Состав волнового пакета со временем меняется: одни ступени заменяются другими, поскольку волновая скорость с — dq/dk необязательно совпадает со скоростью ступеней v == q/k. Термин «кинематический» здесь используют потому, что фундаменталь­ ное уравнение (20.4) не учитывает никаких «динамических» фак­ торов (например, законов механики, движущих сил, диффузион­ ных полей).

Профили кристаллической поверхности. Франк [203] дока­ зал две теоремы, связывающие кинематическую теорию с про­ филями кристаллической поверхности при росте или растворе­

нии.

 

 

 

 

 

Первая теорема заключается в следующем: точки данной

ориентации на некоторой кристаллической поверхности

движутся

в

пространстве

по прямолинейной

траектории, так что dtj/dx =

=

ду/дх + (ду/dt)

(dt/dx); тогда из

соотношений (20.1)

и (20.2)

следует, что

 

 

 

 

 

 

dx

k

с

(20.6)

 

 

 

 

 

 

 

 

 

Эта величина вдоль характеристики dx/dt постоянна, поскольку постоянны k и q. Соотношение у — y(x,t) определяет некоторую поверхность, составленную из характеристик, которые имеют прямолинейные проекции как на плоскость xt, так и на плос-


V. СТРУКТУРА ПОВЕРХНОСТЕЙ РАЗДЕЛА И ПОВЕРХНОСТНАЯ КИНЕТИКА 473

кость ху; следовательно, характеристики, лежащие на поверхно­ сти y(x,t), прямолинейны.

Вторая теорема гласит следующее: если построить полярную диаграмму обратных скоростей роста, измеренных по нормали к поверхности кристалла, то прямолинейная траектория некото­ рой точки данной ориентации на указанной поверхности будет параллельна нормали к полярной диаграмме в точке этой ори­ ентации. На фиг. 29 [144] иллюстрируется развитие кристалличе­ ской формы во времени; видно, что некоторые ориентации посте­ пенно исчезают. Франк и Айвес [206] изучали растворение Ge, а

 

 

а

 

 

 

 

 

б

 

 

Ф и г .

29. Профиль кристаллической

поверхности по

данным кинематиче­

 

 

 

 

 

ской

теории

[144].

 

 

а — изменение

q~1 (17—поток

ступеней)

в

зависимости от ориентации для кубического

 

 

 

кристалла

(схема); б — профили поверхности.

 

Айвес [207]

растворение

LiF;

как

оказалось,

ход

растворения

согласуется

с обеими теоремами Франка.

 

 

Ударные волны. На фиг. 30 [203] приведена

типичная кривая

для

q(k)\

показано также, что

же именно происходит с эшело­

ном ступеней в случае такой зависимости q{k).

Из

соотношения

с=

(dx/dt)

 

== (dq/dk)

следует,

что

при малых

k величина dqjdk

у характеристик больше, благодаря чему их наклон в коорди­ натах t — х меньше (фиг. 30); следовательно, характеристики в левой (хвостовой) части эшелона сходятся и, как показано на фиг. 30, плотность k ступеней претерпевает разрыв. По Франку такой разрыв можно уподобить ударной волне в термодинамике;

как можно

показать, он

движется со скоростью с= (q2

q\)l

l{k% — k\),

где параметры

q и k характеризуют участки на

про­

филе

поверхности, непосредственно прилежащие к точке раз­

рыва

по обе стороны от нее (фиг. 30). Как отметил Франк, что­

бы ведущий

край эшелона

был резким (что желательно при со­

здании ямок

травления),

производная (d2q/dk2)

должна быть

положительной величиной

(это случай, обратный

изображенному


474

 

на фиг. 30). Франк предположил, что такой эффект вызывается

примесями

(см. гл. V I I ) .

Чернов

[17] продолжил исследования ударных волн; в част­

ности, он разработал теорию их трансформации во времени;

были получены

выражения для интервалов времени, необходи­

 

 

 

 

 

мых для развития и исчезновения

 

 

 

 

 

ударных волн.

 

 

 

 

 

 

 

 

 

 

Некоторые

новейшие

разра­

 

 

 

 

 

ботки

кинематической

теории.

 

 

 

 

 

1. Маллинз

и Хирс [208] разрабо­

 

 

 

 

 

тали

детальную

микроскопиче­

 

 

 

 

 

скую

кинематическую

 

теорию,

 

 

 

 

 

выписав в явном виде систему

 

 

 

 

 

дифференциальных

 

уравнений,

 

 

 

 

 

каждое из которых соответствует

 

 

 

 

 

движению одной ступени. Оказа­

 

 

 

 

 

лось,

что эта теория

согласуется

 

 

 

 

 

с франковской

«непрерывной» те­

 

 

 

 

 

орией, где q рассматривается как

 

 

 

 

 

непрерывная функция

k. Получе­

 

 

 

 

 

ны некоторые новые выводы от­

 

 

 

 

 

носительно

поведения

конечных

 

 

 

 

 

эшелонов

ступеней,

 

хвостовая

Ф и г. 30.

Трансформация эшело­

часть

которых,

как

 

оказалось,

на ступеней по данным кинема­

разбивается

на эшелон

с двойны­

тической

теории [203].

 

ми ступенями.

 

 

 

 

 

Показаны

кривая зависимости

q (k),

и др. [209] экспери­

характеристика

х—t,

распределение

2. Бартини

плотности

ступеней

k (х, t) и

разви­

ментально

исследовали

оптиче­

вающиеся профили поверхности с угловой

точкой на

хвостовой

части профиля.

скими

методами

профиль

поверх­

 

 

 

 

 

ности

кристаллов 8-метилнафта-

лина при их росте из спиртовых растворов;

утверждается, что

методика обладает чувствительностью 38 А

(около

двух

перио­

дов решетки данного вещества). Профили кинематических волн наблюдались непосредственно. Эти профили имели гладкий округленный ведущий край и острую угловую точку на хвосто­

вой части, как показано на фиг. 30; это качественно

согласуется

с теорией

Франка.

 

 

 

 

 

Кроме

того, наблюдался

новый

вид

кинематических волн,

прямо

противоположный только что рассмотренному:

скорость

отдельных

ступеней увеличивалась

по

мере их

сближения

(d2q/dk2

есть положительная

величина). Чернов и Будуров [210]

интерпретировали их как «диффузионно-кинематические

волны»,

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


V I

МОРФОЛОГИЧЕСКАЯ УСТОЙЧИВОСТЬ

21. В В Е Д Е Н И Е

Основная задача теории морфологической устойчивости со­

стоит

в том, чтобы

выяснить, устойчива

ли данная конкретная

форма

кристалла,

растущего из расплава, раствора или пара,

по отношению к

малым искажениям.

Теоретический анализ

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

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

В своей классической работе Маллинз и Секерка [211] ввели в теорию роста кристаллов понятие о количественной оценке устойчивости формы роста. До появления этой статьи считалось, что общий вид формы, сохраняющейся при росте, — это иссле­ дованный Хэмом [68, 69] эллипсоид, хотя устойчивость такого эллипсоида не была изучена. Необходимость исследования ус­ тойчивости вызвана тем, что кристаллы часто растут в виде дендритов. Дендриты, или древовидные кристаллы, состоящие из центрального ствола, первичных, вторичных и т. д. ветвей,— частое явление как в природе, так и в лабораторных исследо­ ваниях. Рост дендритов служит ярким примером неустойчивости эллипсоидальной и полиэдрической форм роста, проявляющейся при определенных условиях.

Аналогичный анализ устойчивости проводится в гидродина­ мике (см., например, монографии Ландау и Лифшица [212], Чандрасекара [213], Лина [214], Доннели и др. [215], Стюарта [216]). Чтобы установить момент перехода от ламинарного течения к турбулентному, к известной стационарной скорости потока добавляют малую нестационарную скорость и затем, рассма­ тривая линеаризованное дифференциальное уравнение, которому


476

Р. ПАРКЕР. МЕХАНИЗМЫ РОСТА КРИСТАЛЛОВ

подчиняется эта добавка, выясняют, увеличивается она со вре­ менем или же убывает. Примером задачи о неустойчивости ско­ рости течения может служить задача, решенная Боллом [217], который исследовал перенос вещества через границу газ — жид­ кость с учетом как потока жидкости, так и диффузионного переноса вещества. Оба эти процесса, а именно диффузия и кон­ векция, принимались во внимание и в классической задаче Бенара (см., например, монографию Чандрасекара [213]; см. также гл. V I I I ) , посвященной исследованию устойчивости течения на­ греваемой снизу жидкости.

22.ЗАДАЧА О МОРФОЛОГИЧЕСКОЙ УСТОЙЧИВОСТИ

ВЧИСТОМ ВИДЕ; Р А З Л И Ч Н Ы Е ФОРМЫ РОСТА

Шар. Ради удобства дальнейшего изложения повторим здесь приведенное в гл. I I I доказательство того, что вместо точного уравнения диффузии, учитывающего зависимость от времени, можно взять уравнение Лапласа. Точное уравнение диффузии для шара имеет вид

dC(r,t)

= D V 2 C >

( 2 2 Л )

где C(r,t) есть концентрация, a D — коэффициент диффузии. При этом должны выполняться следующие граничные условия:

С (со, t) =

Cx

 

(22.2)

и

 

 

 

С (R, t) =

С„.

'

(22.3)

Здесь предполагается, что концентрация на поверхности шара радиуса R равна равновесной концентрации С0 . Скорость кри­ сталлизации находится из условия непрерывности потока

D ^ L, = ^ - C o ) f

(22-4)

где Ст — концентрация растворенного вещества в выделяющейся фазе (плотность кристалла). Как было показано в гл. I I I и как можно проверить подстановкой, решение задачи (22.1) — (22.4) имеет вид

Г (г

Л-Г

-

( С ° ~ С о о ) {( R l r )

6 Х Р I ~

fls'V*)2] ~ h*'1 eric (V/l)}

К '

'

°°

 

exp(-Jl2

5 )-VVrfc* 5

где

 

 

 

 

(22.5)

 

 

 

 

 

и

 

 

R =

2X5(Dt)'12

(22.6)

 

 

 

 

 

2Щ1 -п'/2 Я5 ехр(Я2 )ег[с(Я5 )] = 5 ^ ( С ю - С 0 ) / ( С г - С 0 ) . (22.7)