Физика и химия стекла, 2020, T. 46, № 2, стр. 148-163

Влияние переохлаждения и отклонений от стехиометрии на кинетику и механизм кристаллизации стекла As2Se3

Е. В. Школьников *

Санкт-Петербургский государственный лесотехнический университет им. С.М. Кирова
194021 Санкт-Петербург, Институтский пер., 5, Россия

* E-mail: eshkolnikov@yandex.ru

Поступила в редакцию 15.04.19
После доработки 13.11.19
Принята к публикации 05.12.19

Полный текст (PDF)

Аннотация

Проведено сравнительное исследование кинетики массовой объемной кристаллизации массивных стекол As2Se3, As1.9Se3, As2.1Se3 при 240°С с использованием модели Колмогорова–Аврами. Выполнен теоретический анализ экспериментальных зависимостей скорости линейного роста кристаллов As2Se3 от температуры и вязкости стекла As2Se3. Показано, что в стекле As2Se3 при переохлаждениях ∆T = 30–135°С вероятен дислокационный рост кристаллов в форме сферолитов. Для индивидуальных пластинчатых кристаллов As2Se3 с бездислокационными гранями внутри сферолитов при ∆T = 75–135°С вероятен механизм роста с поверхностной двумерной нуклеацией (2Dsg-model).

Ключевые слова: массовая объемная и поверхностная кристаллизация стекол, влияние переохлаждения и отклонений от стехиометрии на кинетику и механизм линейного роста кристаллов

ВВЕДЕНИЕ

В работах [1, 2] исследовано влияние малых добавок олова и свинца на кинетические параметры ступенчатых превращений при изотермической объемной кристаллизации стекол на основе As2Se3. Чистое стекло As2Se3, используемое в электронной технике и волоконной ИК оптике, при термообработке значительно выше температуры стеклования кристаллизуется преимущественно с поверхности в изотермических [2] и неизотермических условиях [3]. Кинетику изотермической кристаллизация стекла As2Se3 исследовали методами измерения плотности, электропроводности и рентгенофазового анализа [2], дифференциальной сканирующей калориметрии и микроскопии с учетом влияния условий синтеза и дисперсности стеклопорошков [46].

Влияние величины переохлаждения ∆Т на механизм и кинетические параметры кристаллизации стекла As2Se3 изучено недостаточно, а влияние отклонений от стехиометрии на кристаллизацию As2Se3 не исследовано.

Цель работы – теоретический анализ результатов экспериментальных исследований кинетики изотермического линейного роста кристаллов в стекле As2Se3, сравнительный анализ кинетики и механизма изотермической кристаллизации стекол As2Se3, As1.9Se3 и As2.1Se3 в интервале температур 240–330°С.

МЕТОДИКА ЭКСПЕРИМЕНТА

Стекла As2Se3, As1.9Se3 и As2.1Se3, называемые ниже стеклами соответственно I, II и III синтезировали методом вакуумной плавки из особо чистых элементных веществ с общей массой 7 г в интервале 500–700°С с последующей закалкой кварцевых ампул с расплавами в воздухе от 700°С. Контроль стеклообразования и однородности закаленных сплавов проводили аналогично [7]. При исследовании поэтапной объемной изотермической кристаллизации стекол в форме шлифованных дисков (4 образца в закрытых бюксах с воздухом) для исключения вклада поверхностной кристаллизации в измеряемые свойства стекол механически удаляли внешний закристаллизованный слой. Эффективную плотность d кристаллизующихся стекол определяли при комнатной температуре гидростатическим взвешиванием в толуоле, микротвердость H измеряли методом Виккерса на приборе ПМТ-3. Рентгенофаовый анализ (РФА) выполняли в основном на дифрактометре УРС-50ИМ аналогично [1, 2]. Дифференциальный термический анализ (ДТА) проводили в вакуумированных кварцевых ампулах на дериватографе системы F. Paulik, J. Paulik, L. Erdey. Удельную электропроводность σ полупроводниковых сплавов измеряли в изотермических условиях в интервале 20–100°С с помощью моста постоянного тока Р4060. Среднюю скорость линейного роста кристаллов определяли, как и ранее [8], с погрешностью ±5–10%, измеряя по ГОСТ 22290-76 с помощью микроскопа толщину закристаллизованного слоя в нескольких точках по периметру скола полированных пластин после термообработки в бюксе с силиконовым маслом.

РЕЗУЛЬТАТЫ И ИХ ОБСУЖДЕНИЕ

При выборе температурного интервала кристаллизации учитывали результаты ДТА при скорости нагревания 6°С/мин: температуры стеклования Тg (174, 170 и 172 ± 2°С), температуры максимума экзотермического эффекта кристаллизации Тmах (327, 322 и 338 ± 2°С) соответственно для стекол As2Se3, As1.9Se3 и As2.1Se3 и температуру плавления Тm(375 ± 2°С) выделяющейся фазы стабильного соединения α-As2Se3 [2].

Таблица 1 .  

Кинетические данные объемной кристаллизации стекол As2Se3, As1.9Se3 и As2.1Se3 при 240°С


этапа
Время,
ч
Н ± 3,
кг/мм2
–lg σ293 lg σ0 ± 0.1 Eσ ± 0.02,
эВ
d ± 0.002,
г/см3
Pd Px Фазовый состав по данным РФА
(Ом–1 см–1) %
Стекло As2Se3
1 0 140 12.1 3.2 1.80 4.607 0   Аморфная фаза
2 12 140 12.1 3.3 1.82 4.607 0 0
3 24 143 12.2 3.5 1.84 4.613 4.0   As2Se3 и аморфная фаза
4 33 141       4.616 6.5  
5 44 140 12.3 3.4 1.83 4.630 16  
6 55 112 12.4 3.6 1.85 4.642 25  
7 60 109       4.671 45  
8 76 106 12.5 3.6 1.86 4.695 60  
9 100 101 12.6 3.5 1.88 4.710 72 65
10 125 72 12.8 3.5 1.95 4.746 97 90
11 150 68       4.748 98.6  
12 190 63 13.4 3.6 2.01 4.749 99.3  
13 220 61 13.5 3.9 2.09 4.750 100 100 As2Se3
             Стекло As1/9Se3
1 0 136 11.9 3.2 1.78 4.606 0   Аморфная фаза
2 2.9 138 11.9 3.2 1.78 4.606 0 0
3 9 137 12.0 3.3 1.80 4.612 4.0   As2Se3 и аморфная фаза
4 12 137       4.622 11  
5 15 136 12.2 3.2 1.81 4.639 23  
6 22 130       4.662 39  
7 30 114 12.5 3.6 1.86 4.695 62  
8 35 105 12.6 3.6 1.87 4.708 71 65
Стекло As2.1Se3
1 0 137 12.0 3.1 1.78 4.609 0   Аморфная фаза
2 18 140 12.2 3.1 1.80 4.609 0 0
3 41 143 12.2 3.5 1.84 4.612 2.0   As2Se3 и аморфная фаза
4 50.5 142       4.614 3.2  
5 60 136 12.3 3.3 1.83 4.619 7.0  
6 69 134       4.625 11  
7 80 134 12.3 3.4 1.84 4.626 12 10

Влияние отклонений от стехиометрии на кинетику массовой изотермической кристаллизации стекла As2Se3.

Исследование поэтапной кристаллизации стекол I–III при 240°С показало, что отклонения от стехиометрии на 1.2 ат. % мышьяка (селена) увеличивали по сравнению с As2Se3 склонность стекол к поверхностной кристаллизации, протекавшей быстрее объемного превращения и затруднявшей его исследование. Удаление поверхностного закристаллизованного слоя у стекол As1.9Se3 и As2.1Se3 приводило к образованию трещин и разрушению образцов до завершения объемной кристаллизации. Полученные результаты, представленные в табл. 1 , проанализированы по методике, описанной в работе [1].

Из данных табл. 1 (см. этапы № 2) следует, что избыток 1.2 ат. % As в 1.5 раза увеличивает, а избыток 1.2 ат. % Se в 4 раза уменьшает индукционный период τинд объемной кристаллизации стекла As2Se3 при 240°С. В течение τинд исследуемые образцы не содержали кристаллических включений, видимых в инфракрасный микроскоп МИК-1, и были рентгеноаморфными; плотность, микротвердость и параметры электропроводности образцов практически не изменялись. Последующее выделение кристаллической фазы As2Se3 в стеклах I–III приводило к постепенному повышению эффективной плотности (на 3.1% при чувствительности метода 0.03%) и энергии активации собственной электропроводности Eσ полупроводниковых стеклокристаллических материалов, а также к понижению их микротвердости и удельной электроповодности. Определенные аналогично [1, 2] значения процента объемной кристаллизации по данным измерения плотности (Pd,) и данным РФА (Px) удовлетворительно согласуются, свидетельствуя о замедляющем действии избытка 1.2 ат. % As и ускоряющем действии избытка 1.2 ат. % Sе на объемную кристаллизацию стекла As2Se3 (табл. 1).

Анализ изотерм α (τ) на рис. 1 проведен на основе уравнения Колмогорова–Аврами для неполных [1] изотермических фазовых превращений

(1)
$\alpha ~ = 1--\exp ( - k\tau _{{{\text{кин}}}}^{n}),~$
где α – степень завершенности кристаллизации по отношению к анализируемой фазе, принятая равной ~0.01 Pd (табл. 1), так как для фазы As2Se3 предельное значение объемной доли равно 1.0 в стекле I и ~0.95 в стелах II и III; k, n – константа скорости и кинетический параметр массовой кристаллизации; τкин – кинетическое время (время превращения), определяемое разностью общего времени термообработки τ и индукционного периода τинд кристаллизации. При выводе уравнения Колмогорова–Аврами предполагались отсутствие индукционного периода, неизменность скорости нуклеации или концентрации готовых зародышей, а также скорости роста кристаллов в течение всего процесса превращения. Следует заметить, что при наличии индукционного периода на экспериментальных изотермах α (τ) использование общего времени термообработки τ вместо времени фиксируемого превращения τкин приводит к получению завышенных значений параметра n массовой кристаллизации стекол [2, 4, 5].

Рис. 1.

Изотермы степени объемной кристаллизации стекол As2Se3 (1), As1.9Se3 (2) и As2.1Se3 (3) при 240°С по данным измерения эффективной плотности закаленных образцов.

Найденные двукратным логарифмированием уравнения (1) с использованием для величины α данных Pd (табл. 1) и метода наименьших квадратов (рис. 2) значения кинетических параметров равны: n = 2.05 ± 0.11, 2.00 ± 0.13 и 2.1 ± 0.2, –lg k = 11.1 ± 0.6, 10.0 ± 0.6 и 12.1 ± 1.2 для кристаллизации стекол соответственно As2Se3, As1.9Se3 и As2.1Se3. Полученные значения кинетического параметра n указывают на преимущественно гетерогенное зарождение и двумерный рост пластинчатых кристаллов As2Se3 в стеклах I–III. Морфология роста в виде пластинок характерна для монокристаллов As2Se3 и установлена прямым методом электронной микроскопии для растущих индивидуальных кристаллов в стекле As2Se3 в интервале температур 240–350°С [4].

Рис. 2.

Графическое решение уравнения α = 1 – exp($ - k~\tau _{{{\text{кин}}}}^{n}$ ) для объемной изотермической кристаллизации стекол As2Se3 (1), As1.9Se3 (2) и As2.1Se3 (3) при 240°С.

При гетерогенном объемном зарождении и двумерном росте кристаллов (n = 2) константа скорости массовой изотермической кристаллизации определяется формулой [1, 2]

(2)
$k = \pi lN{{u}^{2}},~$
где l – толщина растущей пластинки кристалла(~10–7 cм для As2Se3), N – объемная концентрация готовых центров кристаллизации, u – cкорость линейного роста кристаллов. Оценка концентрации готовых центров N по формуле (2) с использованием приведенных выше значений k и экспериментальных значений u (см. табл. 2) дает для объемной изотермической кристаллизации стекла As2Se3 при 240°С 1.3 × 109 см–3, а по скорости гомогенного зарождения (I ~ 3 см –3 с–1) и индукционному периоду (~5 × 104 с) –2 × 105 см –3 [2]. Поэтому вкладом гомогенного зарождения (<0.02% при 240°С) кристаллов As2Se3 в концентрацию готовых центров N в формуле (2) можно пренебречь. Аналогичная оценка величины N для объемной кристаллизации стекол As1.9Se3 и As2.1Se3 при 240°С дает соответственно 1.8 × 109 и 5 × 108 см–3. Индукционный (скрытый) период изменяется антибатно концентрации готовых центров объемной кристаллизации в ряду стекол As1.9Se3 → As2Se3 → As2.1Se3 и равен, соответственно, 2.9 ± 0.3, 12 ± 1 и 18 ± 2 ч при 240°С. Что касается природы гетерогенных центров N, то ими могут быть наночастицы малорастворимых селенидов примесных тяжелых металлов (Sn, Pb, Bi, In по 10–4 мас. %) в использованном для синтеза селене квалификации “ос. ч. 17–3”. Согласно [1, 2], малые добавки селенидов этих металлов к стеклу As2Se3 способны инициировать гетерогенное зарождение кристаллов As2Se3 и значительно уменьшать индукционный период объемной кристаллизации стекла.

Влияние переохлаждения и отклонений от стехиометрии на кинетику и механизм линейного роста кристаллов в стекле As2Se3

Таблица 2.  

Сравнительное влияние температуры и отклонений от стехиометрии на индукционный период и скорость линейного изотермического роста индивидуальных кристаллов в сферолитах [4], сферолитов [6] и роста поверхностного кристаллического слоя в стекле ${\text{A}}{{{\text{s}}}_{2}}{\text{Se}}_{3}^{*}$

Стекло T ± 1 °C τинд, u uТ (ф.3 ) lg uТ (ф.3) lg u
ч с cм/c
As2Se3 240 27 9.7 × 104 2.0 × 10–7
3.4 × 10–7 [4]
1.3 × 10–7 [6]
9 × 10–8 –0.35
–0.58
–0.16
270 5.0 1.8 × 104
         ~3 × 103 [4]
1.4 × 10–6
2.4 × 10–6 [4]
1.2 × 10–6 [6]
8.2 × 10–7 –0.23
–0.46
–0.12
310 0.5 1.8 × 103 1.4 × 10–5
1.1 × 10–5 [6]
1.2 × 10–5 –0.06
+0.08
330 0.083 3.0 × 102 3.3 × 10–5
2.0 × 10–5 [6]
2.7 × 10–5 –0.11
+0.13
As1.9Se3 240 7 2.4 × 104 4.2 × 10–7 1.4 × 10–7 –0.46
270 2.0 7.2 × 103 2.8 × 10–6 1 × 10–6 –0.4
As2.1Se3 240 39 1.4 × 105 1.3 × 10–7 6.4 × 10–8 –0.30
270 7 2.5 × 104 9.8 × 10–7 6.3 × 10–7 –0.19

* u – экспериментальные данные Э.Ю. Бессоновой и автора в первых строках, в следующих строках – из работ [4] и [6], uТ – расчет автора по теоретической формуле (3).

Экспериментальные зависимости толщины l поверхностного закристаллизованного слоя стекол I–III от времени термообработки τ при 270°С представлены на рис. 3. Линейный характер этих и аналогичных зависимостей при других температурах, наблюдаемый по истечении индукционного периода τинд, указывает на кинетический контроль скорости роста поверхностного кристаллического слоя в стеклах I–III. Первичное зарождение кристаллов происходило вероятнее всего гетерогенно на дефектах полированной поверхности [9] и в нашей работе не исследовалось. При выборе температурного интервала кристаллизации стекол на основе As2Se3 учитывали возможность взаимосвязи (приближенного равенства) кинетических параметров роста кристаллов и вязкого течения стекол [8] при температурах выше температуры декаплинга диффузии и вязкости [10].

Рис. 3.

Зависимости толщины повехностного закристаллизованного слоя от времени термообработки при 270°C в стеклах: As2Se3 (1), As1.9Se3 (2), As2.1Se3 (3).

Найденные из зависимостей l(τ) значения τинд и средней скорости линейного изотермического роста u закристаллизованного слоя в стеклах I–III приведены в табл. 2 и на рис. 4. Как видно, экспериментальные значения u скорости поверхностной кристаллизации полированного стекла As2Se3 в настоящей работе являются промежуточными между результатами работ [6] и [4]. В работе [6] определяли линейную скорость роста сферолитов (сферических агрегатов пластинчатых кристаллов) в объеме стекла As2Se3 с использованием оптической микроскопии, а в работе [4] – максимальные значения скорости анизотропного роста индивидуальных кристаллов As2Se3 в сферолитах, полученные методом сканирующей электронной микроскопии с предварительным селективным травлением аморфной фазы. При наличии согласованности экспериментальных зависимостей u (T) наблюдается заметное различие скорости роста кристаллов u (рис. 4, кривые 1–3), вызванное в основном различием методик получения образцов As2Se3 и измерения скорости линейного роста кристаллов в них.

Рис. 4.

Температурные зависимости скорости роста u (1–3), приведенной скорости uc(см. текcт) роста кристаллов(4–6) и величины u/∆T2 (7–9) в стекле As2Se3 (анализ эксперименентальных данных u (см. табл. 2) из работ [4] (кривые 1, 4, 7), [6] (кривые 2, 5, 8), данных Э.Ю. Бессоновой и автора (кривые 3, 6, 9)).

В отличие от стекла стехиометрического состава As2Se3 с тригонально-пирамидальными структурными единицами (“с. е.”) AsSe3/2, избыточный мышьяк (1.2 ат. %) в стекле As2.1Se3 образует с. е. As2Se4/2, характерные для трудно кристаллизующегося стекла AsSe, а избыточный селен (1.2 ат. %) в стекле As1.9Se3 образует цепочечные с. е. SeSe2/2, характерные для тригонального и довольно легко кристаллизующегося стеклообразного селена. По сравнению с As2Se3 стекло As2.1Se3 имеет более высокие, а стекло As1.9Se3 – более низкие значения динамической вязкости η [11] и индукционного периода кристаллизации τинд (рис. 3, кривые 1–3, табл. 2). Скорость линейного роста поверхностного кристаллического слоя u при 270°С в стекле As2.1Se3 примерно в 1.5 раза меньше, а в стекле As1.9Se3 больше в 2 раза, чем в стекле стехиометрического состава As2Se3 (табл. 2).

Для определения механизма роста кристаллов в исследуемых стеклах использовали сравнение экспериментально найденных температурных зависимостей скорости роста кристаллов u (T) с теоретическими зависимостями, ожидаемыми для известных механизмов. Нормальный рост кристаллов с малой энтропией плавления (Sm < 2R) происходит равновероятно в любой точке атомно-шероховатой грани по нормали к ней и наблюдается в жидкостях с малой вязкостью [12]. Для кристаллов с ∆Sm > 4R, включая As2Se3 (63.5 Дж/моль К [13]) плотно упакованные грани согласно [12] должны быть атомно-гладкими, а анизотропия роста – достаточно выраженной и определяемой дислокациями реального кристалла.

Из сравнения результатов измерения [6] и расчета lgu (табл. 2, 3) следует, что скорость стационарного роста пластинчатых (plate-like) кристаллов As2Se3 в стекле того же состава при t = 240–345°С (∆T = 135–30°C) может быть удовлетворительно рассчитана, если в теоретическую формулу [14] ввести вероятностный фактор f [13] роста на ступенях винтовых дислокаций согласно

(3)
где R – газовая константа, Т – абсолютная температура, f (TТ) вероятностный фактор, или плотность точек роста (см. ниже и табл. 3); δ – элементарное перемещение фронта кристаллизации, принятое равным сумме эффективной толщины полимерно-сетчатого слоя и кратчайшего расстояния между слоями, параллельными плоскости спайности (010). Величина δ оценена по структурным данным [15] равной ~0.68 нм. Уменьшение свободной энергии при кристаллизации ∆G на моль As2Se3 рассчитано аналогично [8] и представлено в табл. 3. В качестве величины ν принята основная частота валентных колебаний мостиковых атомов cелена, равная согласно колебательным спектрам стекла As2Se3 (6.5–6.4) × 1012 с–1 в интервале температур 220–358°С [16]. Свободная энергия активации нереконструктивного роста кристаллов при расчете скорости дислокационного роста по формуле (3 ) принималась, как и ранее [8], равной свободной энергии активации вязкого течения стекол $\Delta G_{\eta }^{\# }\left( Т \right)$, найденной по формуле (4) и представленной в табл. 3. Величина $\Delta G_{\eta }^{\# }$ определена нами по экспериментальным данным вязкости η (Т) стекла As2Se3 [11] с использованием наиболее общего уравнения Эйринга
(4)
$\Delta G_{\eta }^{\# } = RT[{\text{ln}}\eta - {\text{ln}}({{{{N}_{{\text{A}}}}h} \mathord{\left/ {\vphantom {{{{N}_{{\text{A}}}}h} {{{V}_{\eta }}}}} \right. \kern-0em} {{{V}_{\eta }}}})],$
где Vη – мольный объем мостиковых атомов селена (~19 см3/ грамм-атом Se, найденный по структурным данным [15]; h – постоянная Планка, NA – число Авогадро.

Таблица 3.  

Результаты расчета и измерений скорости роста кристаллов в стеклах As2Se3, As2.1Se3 и As1.9Se3

T, °C lg η (η, П) [11] $G_{\eta }^{\# }$ ± 0.9 $G_{{\text{a}}}^{{''}}$ ± 0.8 –∆G ± 0.2 f × 102 –lg u (u, см/c)
кДж/моль расчет по формуле (3) данные измерений
[6] [4] *
Стекло As2Se3
240 7.44 109.5. 108.9 6.8 3.2 7.05 6.89 6.47 6.70
250 7.02 107.4. 108.1 6.4 3.0 6.71 6.66 6.17 6.42
260 6.65 105.6. 106.8 6.0 2.8 6.38 6.38 5.89  
270 6.30 104.0. 103.7 5.4 2.6 6.08 5.96 5.62 5.85
280 5.90 101.7 102.1 5.2 2.4 5.75 5.70 5.34  
290 5.50 99.2 100.3 4.7 2.2 5.41 5.44 5.10 5.27
300 5.19 97.6 98.6 4.2 2.0 5.28 5.22 4.90  
310 4.85 95.5 96.7. 3.7 1.7 4.91 4.99   4.85
317 4.60 93.8. 95.5 3.4 1.6 4.76 4.88    
330 4.23 91.6 93.1 2.7 1.3 4.56 4.69   4.48
335 4.08 90.1 92.6 2.4 1.1 4.47 4.67    
345 ~3.8 88.7 90.4 1.8 0.85 4.32 4.60    
350 ~3.7 88.1 89.4 1.6 0.74 4.43 4.59    
Стекло As2.1Se3
240 7.48 109.8   6.2 2.9 7.19     6.89
270 6.34 104.4   4.9 2.3 6.20     6.01
Стекло As1.9Se3
240 7.1 106.1   5.7 2.8 6.84     6.38
270 ~6.1 ~102   4.4 2.2 ~6.0     5.55

* Данные Э.Ю. Бессоновой и автора.

Вероятностный фактор f для дислокационного роста кристаллов возрастает с увеличением ТТ (табл. 3) согласно выражению [8]

(5)
$f = \frac{{\delta T\Delta T\Delta {{H}_{m}}}}{{4\pi \sigma T_{m}^{2}{{V}_{M}}}},$
где σ – поверхностное натяжение стекло–кристалл для торцевой грани пластинчатого кристалла As2Se3(1.5 × 10–5 Дж/cм2) определялось с учетом анизотропии роста кристаллов As2Se3 в интервале температур 240–350°С [4, 17] как утроенная величина σl для боковых граней, найденная по cоотношению Cкапского–Тернбалла (Skapski–Turnbull relation) σ = β∆Hm/($N_{A}^{{1/3}}V_{m}^{{2/3}}$) [14]. Здесь β – константа, равная 0.2; ∆Hm – молярная энтальпия плавления, Tm – температура плавления, Vm – мольный объем кристалла. Величины ∆Hm, Tm и Vm для кристалла As2Se3 приняты равными, соответственно, 41 кДж/моль, 648 K и 81.4 см3/моль. По данным микроскопии [17] рост кристаллов со слоистой структурой As2Se3 в стекле 0.9As2Se3 · 0.1Sb2Se3 в перпендикулярном направлении к плоскости спайности (010) происходил в три раза медленнее, чем вдоль плоскости спайности.

Для стекол As1.9Se3 и As2.1Se3 расчеты с использованием формул (3)–(5) являются приближенными (табл. 2, 3) из-за неизбежных допущений при переходе к ним от стеклообразующего соединения As2Se3 и осложнены недостатком количественной информации о параметрах , необходимых при выполнении подобных расчетов. Динамическую вязкость η указанных стекол при 240 и 270°С оценивали методом интерполяции экспериментальных данных [11] для стекол соответственно As1.84Se3 и As2Se3, As2.19Se3 и As2Se3 и использовали для расчета величины $\Delta G_{\eta }^{\# }$ по уравнению (4) и кинетического фактора exp() в формуле (3). Для определения величины переохлаждения ∆Т, необходимой для расчетов вероятностного фактора f по формуле (5) и термодинамического фактора 1 – exp(∆G/RT ) в формуле (3) проводили оценку температур ликвидуса (363 и 368°С) по известному логарифмическому уравнению Шредера в интегральной форме соответственно для составов As1.9Se3 и As2.1Se3 в области максимума при составе устойчивого соединения As2Se3 (375 ± 2°C) на диаграмме состояния системы As–Se [18]. Молярные доли основного компонента As23, равные 0.86 и 0.92, учитывали также в расчетах изменения химического потенциала µ при кристаллизации As2Se3 соответственно в стеклах As1.9Se3 и As2.1Se3 с заменой ∆G на ∆µ (As2Se3) в формуле (3).

Результаты измерения и теоретического расчета скорости u линейного роста кристаллов As2Se3 в стеклах As2.1Se3 и As1.9Se3 (табл. 2, 3) с использованием формул (3)–(5) указывают на возможность дислокационного механизма роста кристаллов в этих стеклах с экспериментальной энергией активации соответственно 154 и 146 ± 2 кДж/моль в исследованном интервале температур 240–270°С (∆Т = 93–123°). Последующий анализ кинетики роста кристаллов (табл. 4, рис. 4–7) выполнен для стекла стехиометрического состава As2Se3, изученного прямыми методами микроскопии в широком интервале температур.

Таблица 4.  

Экспериментальная энергия активации роста кристаллов в стекле ${\text{A}}{{{\text{s}}}_{2}}{\text{Se}}_{3}^{*}$

Энергия активации (кДж/моль) в указанном интервале переохлаждений
способ 1 способ 2 способ 3 cпособ 4
Еа Т 0, °С Т 0, °С Т 0, °С Т 0, °С
140–150 75–135 153 ± 4 45–135 176 ± 5 45–135 190 ± 5 45–135
142–158 75–135 156 ± 3 75–135 179 ± 4 75–135 192 ± 4 75–135
150–162 75–135 167 ± 5 30–135 193 ± 5 10–135 207 ± 5 10–135

* Анализ экспериментальных данных настоящей работы в первых строках, в следующих строках – анализ данных из работ [4] и [6].

Рис. 5.

Графическое решение уравнения uη = C(∆T)n по экспериментальным данным вязкости η [11] и скорости изотермического линейного роста кристаллов u в стекле As2Se3 (анализ данных u из работ [4] (1), [6] (2), Э.Ю. Бессоновой и автора (3)).

Рис. 6.

Графическое решение уравнения u/[1 – exp(∆G/RT)]TT = bm по экспериментальным данным вязкости η [11] и скорости изотермического линейного роста кристаллов u в стекле As2Se3 (обозначения кривых 1–3 см. рис. 5).

Рис. 7.

Графическое решение уравнения u = С ехр(${{--\Delta G_{\eta }^{\# }} \mathord{\left/ {\vphantom {{--\Delta G_{\eta }^{\# }} {RT}}} \right. \kern-0em} {RT}}$)exp(${{--B} \mathord{\left/ {\vphantom {{--B} {T\Delta T}}} \right. \kern-0em} {T\Delta T}}$) по экспериментальным данным вязкости [11] и скорости изотермического линейного роста кристаллов u в стекле As2Se3 (обозначения кривых 1–3 см. рис. 5).

Кинетика и механизм изотермического роста кристаллов в стекле As2Se3

Из сопоставления рассчитанных по формуле (3) и экcпериментальных значений скорости u линейного роста кристаллов в стекле As2Se3 (табл. 3) следует, что в интервале температур 300–240°С c увеличением переохлаждения от 75 до 135°C монотонно увеличивается различие результатов расчета и измерения в работе [4] максимальной скорости роста индивидуальных кристаллов в сферолитах As2Se3. Вместе с тем наблюдается удовлетворительное согласие результатов расчета и экспериментальных значений u, недавно полученных в работе [6], а также Э.Ю. Бессоновой и автором (табл. 2, 3).

Для определения кинетических параметров роста кристаллов в стекле As2Se3 выполнили, как и ранее [8], анализ различными способами экспериментальных зависимостей скорости роста кристаллов от температуры (рис. 4, кривые 1–3). Как видно из табл. 4, эффективная энергия активации роста кристаллов Еа, определенная из криволинейных зависимостей lg u(T1), возрастает при увеличении переохлаждения ∆Т (см. способ 1). Графические решения уравнения (3) в координатах lg uc = f(Т1), где uc = = u/[1 – exp(∆G/RT)] и в координатах lg uc/∆Т = f( 1) позволяют в соответствии с формулами (3) и (5) исключить влияние не только термодинамического фактора 1 – exp(∆G/RT) (способ 2), но и переохлаждения ∆Т (см. числитель в формуле (5), способ 3) на скорость роста кристаллов As2Se3 и получить прямолинейные зависимости в интервале температур 240–345°С (рис. 3, линии 4–6) и найти из них значения энтальпии активации линейного роста кристаллов (табл. 4, способы 2 и 3). Близкие к ним результаты (табл. 4, способ 4) получены графическим решением в координатах lg u/∆Т2 = φ(Т 1) (рис. 4, линии 7–9) уравнения

(6*)
соответствующего при небольших ∆G полученному ранее [13] комбинированием формул (3) и (5) приближенному теоретическому уравнению для дислокационного роста кристаллов
(6)
где Aэ = Aт exp (), Aт${{{{\delta }^{2}}\Delta S_{m}^{2}\nu } \mathord{\left/ {\vphantom {{{{\delta }^{2}}\Delta S_{m}^{2}\nu } {(4\pi \sigma {{V}_{M}}R{{T}_{m}})}}} \right. \kern-0em} {(4\pi \sigma {{V}_{M}}R{{T}_{m}})}}.$ Найденное из сопоставления экспериментального Aэ по данным [6] и теоретического множителя Aт среднее значение энтропии активации дислокационного роста кристаллов равно 195 ± 3 Дж/(моль К) в стекле As2Se3 (T = 240–345°C). Среднее значение энтропии активации вязкого течения $\Delta S_{\eta }^{\# }$ в этом интервале температур, определенное из температурной зависимости $\Delta G_{\eta }^{\# }$ (табл. 3), равно 198 ± 3 Дж/(моль К).

Используя найденные величины и известное соотношение = , можно рассчитать отсюда по значениям (табл. 4, способ 4) свободную энергию активации роста кристаллов . Полученные значения (107.0 и 86.5.0 кДж/моль, соответственно, при 240 и 345°С) удовлетворительно согласуются с значениями (табл. 3), вычисленными непосредственно по теоретическим значениям предэкспоненциального множителя и экспериментальным значениям скорости роста кристаллов u [6] в формуле (3) при f < 1 (см. табл. 3).

Для определения взаимосвязи кинетических параметров кристаллизации и вязкого течения провели анализ экспериментальных зависимостей u(η) для стекла As2Se3. С помощью известного соотношения Стокса–Эйнштейна между коэффициентом самодиффузии и вязкостью переохлажденной жидкости получили на основе теоретической модели роста кристаллов [14] выражение

(7)
$u = \frac{{fkT}}{{3\pi \eta {{\delta }^{2}}}}\left[ {1 - \exp ({{\Delta G} \mathord{\left/ {\vphantom {{\Delta G} {RT)}}} \right. \kern-0em} {RT)}}} \right],$
где f, δ, ∆G аналогичны величинам в формуле (3), k – константа Больцмана.

Выражение (6) при небольших ∆G переходит с учетом формулы (5) в приближенную зависимость

(8)
$u\eta ~\,\, = c{{(\Delta T)}^{n}},~$
где с – константа стекла, n ≈ 2 для дислокационного роста кристаллов. Графическое решение (рис. 5) уравнения (8) по экспериментальным данным вязкости η [11], скорости линейного роста индивидуальных кристаллов в сферолитах [4], скорости роста сферолитов [6] и поверхностного закристаллизованного слоя (данные Э.Ю. Бессоновой и автора, см. табл. 3) дает, соответственно, значения n = 2.6, 2.0 и 2.1 при ∆T = 75–135, 25–135 и 45–135°C для кристаллизации стекла As2Se3. Полученные значения n = 2.0 и 2.1 (линии 2 и 3, рис. 5) согласуются с дислокационным механизмом роста кристаллов в стеклах [8, 1214]. К такому же заключению о механизме роста кристаллов в стекле As2Se3 приводит анализ экспериментальных зависимостей u (η) графическим решением (рис. 6) уравнения
(9)
${u \mathord{\left/ {\vphantom {u {[1--\exp ({{\Delta G} \mathord{\left/ {\vphantom {{\Delta G} {RT}}} \right. \kern-0em} {RT}})]}}} \right. \kern-0em} {[1--\exp ({{\Delta G} \mathord{\left/ {\vphantom {{\Delta G} {RT}}} \right. \kern-0em} {RT}})]}}T\Delta T~\,\, = {b \mathord{\left/ {\vphantom {b {{{\eta }^{m}}}}} \right. \kern-0em} {{{\eta }^{m}}}},$
полученного комбинированием формул (7) и (5) с введением дополнительного параметра m (показателя степени при вязкости η) для проверки соотношения (7) между скоростью роста кристаллов и вязкостью стекла [13]. В уравнении (9) m и b – константы в нешироком температурном интервале, определяемые свойствами стекла и образующегося кристалла [19]. Найденные значения показателя m по экспериментальным данным u [6] Э.Ю. Бессоновой и автора равны, соответственно, 0.96 и 0.93 и близки к теоретическому значению m = 1 выше температуры декаплинга диффузии и вязкости. Показатель m, определенный по данным скорости роста индивидуальных кристаллов в сферолитах As2Se3 [4] и данным вязкости [11], заметно ниже (0.87). Согласно электронно-микроскопическим наблюдениям [4] в интервале 240–350°С изотермический рост индивидуальных кристаллов As2Se3 в сферолитах происходил двумерно в форме пластин постоянной толщины. Поэтому провели анализ экспериментальных данных u(T) в табл. 3 на возможность механизма роста кристаллов с поверхностной двумерной нуклеацией (two-dimensional surface nucleated growth), или 2Dsg-model [20, 21], экспериментально доказанного при изучении бездислокационных граней кристаллов серебра [22] и весьма вероятного для процесса поверхностной кристаллизации стекла TlAsSe2 при переохлаждениях больше 120°C [8], а также стекла PbO · SiO2, по крайней мере, при малых переохлаждениях [23].

2Dsg-механизм вероятен при ∆Sm > 4R для послойного роста анизотропных кристаллов, лимитируемого образованием двумерных зародышей (ступеней роста) на атомно-гладкой поверхности раздела кристалл–жидкость по уравнению

(10)
$u = C\exp ({{--E} \mathord{\left/ {\vphantom {{--E} {RT}}} \right. \kern-0em} {RT}})\exp ({{--B} \mathord{\left/ {\vphantom {{--B} {T\Delta T}}} \right. \kern-0em} {T\Delta T}}),$
где C – коэффициент, слабо зависящий от температуры, E – энергия активации кристаллизации (для As2Se3 принята равной $\Delta G_{\eta }^{\# }$, табл. 3); B = πδσ2VM/3kΔSm; δ – толщина двумерного зародыша. Графическое решение уравнения (10) в координатах ln u + ($\Delta G_{\eta }^{\# }$/RT) = = f(105/TT) при наличии 2Dsg-механизма роста кристаллов должно давать прямые линии [8]. Как видно из рис. 7 (линия 1), прямолинейная зависимость наблюдается лишь при использовании экспериментальных данных [4] для роста индивидуальных кристаллов As2Se3 в сферолитах. Аналогичные по форме зависимости получены при построении графика в координатах ln ucη/T = f (105/TT) по экспериментальным данным [4, 6, 11] и настоящей работы. Таким образом, при одинаковых переохлаждениях 75–135°C в стеклах As2Se3 послойный рост кристаллов может происходить по механизму поверхностной двумерной нуклеации [4] или по дислокационному механизму (анализ данных [6] и настоящей работы). По нашему мнению, отмеченное различие механизмов роста кристаллов связано с тем, что авторы работы [4] измеряли в стекле As2Se3, полученном медленным охлажденном расплава от 500°С, максимальную скорость (параллельно плоскости спайности) бокового двумерного роста совершенных и чистых (99.999%) индивидуальных кристаллов As2Se3 с плоскими треугольными гранями без готовых ступеней роста, образуемых дислокациями (ступени роста эти авторы наблюдали в электронный микроскоп только на некоторых кристаллах при ∆T < 25°C). В работе [6] и в настоящей работе измеряли усредненную скорость анизотропного роста в основном реальных кристаллов, определяющих соответственно линейный рост сферолитов и поверхностного закристаллизованного слоя.

Из прямолинейной зависимости на рис. 7 (линия 1) для роста кристаллов As2Se3 [4] при ∆Т = 75–135°C найдены величины В = 2.0 × 106 К2 и $\delta \sigma _{e}^{2}$ = 2.1 × 10–18 Дж2/cм3. Отсюда при δ = 6.8 × 10–8 см получили оценку поверхностного натяжения стекло–кристалл σ = 5.5 × 10–6 Дж/cм2 для поверхностной двумерной нуклеации и послойного роста кристаллов As2Se3. Графическое определение σ в координатах ln ucη/T = f(105/TT) по экспериментальным данным [4, 11] дает значение 4.8 × 10–6 Дж/cм2, близкое к значению, определенному с использованием рис. 7, а также к величине σl (5.0 × 10–6 Дж/cм2) для боковых граней пластин, найденной по cоотношению Cкапского–Тернбалла (см. выше).

Для изотермического линейного роста сферолитов [6] и поверхностного закристаллизованного слоя (данные Э.Ю. Бессоновой и автора) в стекле As2Se3 при переохлаждении 25–135°C (T < Tu max) поверхностная двумерная нуклеация мало вероятна (рис. 7, кривые 2 и 3) и более вероятен механизм роста кристаллов на винтовых дислокациях в соответствии с формулами (3) и (5). Наличие у выхода винтовой дислокации на грань кристалла незарастающей ступени способно обеспечить согласно [12] послойный рост кристаллов без образования двумерных зародышей.

При температурах выше Tu max(350°С) в массивных стеклах и стеклопорошках As2Se3 [4] преобладает одномерный рост индивидуальных кристаллов в форме прутков (rod-like) с ориентацией по радиальным направлениям в сферолитах.

ЗАКЛЮЧЕНИЕ

При отклонениях от стехиометрического состава As2Se3 избыток 1.2 ат. % Se уменьшает в 4 раза, а избыток 1.2 ат. % As в 1.5 раза увеличивает индукционный период объемной и поверхностной кристаллизации стекол при 240°С. С использованием модели Колмогорова–Аврами установлено преимущественно гетерогенное зарождение и двумерный рост пластинчатых кристаллов фазы As2Se3 в cтеклах.

Выполнен теоретический анализ экспериментальных зависимостей скорости линейного роста кристаллов As2Se3 от температуры и вязкости стекла As2Se3. В интервале температур 240–345°С (переохлаждение ∆T = 135–30°C) для линейного роста сферолитов и поверхностного закристаллизованного слоя в массивном стекле As2Se3 вероятен дислокационный механизм. Для индивидуальных кристаллов с бездислокационными гранями в сферолитах As2Se3 вероятен механизм роста в форме пластин (plate-like) с поверхностной двумерной нуклеацией (2Dsg-model).

Скорость линейного роста кристаллов в стекле As2Se3 при 240–345°C, а в стеклах As2.1Se3 и As1.9Se3, по крайней мере, при 240–270°С теоретически рассчитывается, если в формулу Тернбала–Коэна ввести вероятностный фактор f(ТT) роста на ступенях винтовых дислокаций, а свободную энергию активации роста кристаллов принять равной кинетическому барьеру при вязком течении стекла $\Delta G_{\eta }^{\# }(T)~$ в уравнении Эйринга.

Список литературы

  1. Школьников Е.В. Кинетика ступенчатой объемной кристаллизации стекол AsSe1.5Snx (x ≤ 0.28) // Физ. и хим. стекла. 2017. Т.43. №2. С. 172–183. [Shkol`nikov E.V. Kinetics of Stepwise Bulk Crystallization of AsSe1.5Snx (x ≤ 0.28) Glasses // Glass Phys. Chem. 2017. V.43. № 2. P. 137–145.].

  2. Школьников Е.В Влияние добавок Pb на кинетику и механизм кристаллизации стекла As2Se3 // Физ. и хим. стекла. 2019. Т. 45. № 1. С. 16–28. [Shkol`nikov E.V. Influence of Pb Additives on the Kinetics and Mechanism of As2Se3 Glass Crystallization // Glass Phys. Chem. 2019. V. 45. № 1. P. 10–18].

  3. Svoboda R.A., Malek J. Non-Isothermal crystallization kinetics of As2Se3 glass studied by DSC // Thermochem. Data. 2014. V. 579. № 1. P. 56–63.

  4. Henderson D.W., Ast D.G. Viscosity and crystallization kinetics of As2Se3 // J. Non-Crystalline Solids. 1984. V. 64. № 1. P. 43–70.

  5. Cernoskova E., Holubova J., Cernosek Z. Crystallization kinetics of glassy As2Se3 // J. Therm. Anal. Calorim. 1999. V. 56. № 2. P. 423–428.

  6. Malek J., Shanelova J., Martinkova S., Pilny P., Kostal P. Crystal growth velocity in As2Se3 supercooled liquid // Crystal Growth Design. 2017. V. 17. № 9. P. 4990–4999Š.

  7. Школьников Е.В. Структурно-химические особенности кристаллизующихся стекол систем Se–Sn и As2Se3–Sn // Физ. и хим. стекла. 2013. Т.39. №4. С. 529–537. [Shkol`nikov E.V. Snructural-Chemical Features of Crystallizing Glasses of Se–Sn and As2Se3–Sn Systems // Glass Phys. Chem. 2013. V. 39. № 4. Р. 364–371.].

  8. Школьников Е.В. Структурно-химические особенности и механизм кристаллизации стекол TlAsX2(X = S, Se) // Физ. и хим. стекла. 2016. Т.42. № 4. С. 500–511. [Shkol`nikov E.V. Snructural-Chemical Features and Mechanism of Crystallization of TlAsX2 Glasses (X = S, Se) // Glass Phys. Chem. 2016. V. 42. № 4. Р. 363–371.].

  9. Филипович В.Н., Калинина А.М., Фокин В.М., Юрицын Н.С., Сычева Г.А. Кинетические закономерности зарождения кристаллов в объеме и на поверхности силикатных стекол // Физ. и хим. стекла. 1999. Т. 25. № 2. С. 327–336.

  10. Schmelzer Ju. W.P., Abyzov A.S., Fokin V.M., Schick Ch., Zanotto E.D. Crystallization in glass-forming liquids: Effects of decoupling of diffusion and viscosity on crystal growth // J. Non-Crystalline Solids. 2015. V. 429. P. 45–53.

  11. Немилов С.В., Петровский Г.Т. Исследование вязкости стекол системы селен–мышьяк // Журн. прикл. химии. 1963. Т. 36. № 5. С. 977–981.

  12. Jackson K.A., Uhlmann D.R., Hunt J.D. On the nature of crystal growth from the melt // J. Crystal Growth. 1967. V. 1. № 1. P. 1–36.

  13. Школьников Е.В. О взаимосвязи структурно-химических особенностей и кинетических параметров кристаллизации стекла // Стеклообразное состояние. Тр. V11 Всес. Совещ. Л.: Наука, 1983. С. 131–135.

  14. Turnbull D., Cohen M.H. Crystallization kinetics and glass formation // In: Modern aspects of the vitreous state. London, 1960. P. 38–62.

  15. Zallen R., Slade M.L., Ward A.T. Lattice vibrations and interlayer interactions in crystalline As2S3 and As2Se3 // Phys. Rev. B. 1971. V. 3. № 12. P. 4257–4273.

  16. Arai T., Komiya S., Kudo K. Temperature dependence of vibrational spectra in crystalline, amorphous and liquid As2Se3 // J. Non-Crystalline Solids. 1975. V. 18. № 2. P. 289–294.

  17. Platakis N.S., Gatos R.C. Devitrification characteristics of the semiconductor system (1 –x)As2Se3xSb2Se3 // J. Electrochem. Soc. 1972. V. 119. № 7. P. 914–920.

  18. Виноградова Г.З. Стеклообразование и фазовые равновесия в халькогенидных системах. Двойцые и тройные системы. М.: Наука. 1984. 176 с.

  19. Школьников Е.В. К кинетике роста кристаллов в стеклах M2O ⋅ 2 SiO2 (M = Li, Na, K) // Физ. и хим. стекла. 1980. Т. 6. № 2. С. 153–163.

  20. Странский И.Н., Каишев Р. К теории роста кристаллов и образования кристаллических зародышей // Успехи физ. наук. 1939. Т. 21. Вып. 4. С. 408–465.

  21. Smith G.L., Weinberg M.C. Experimental rtest of surface nucleated crystal growth model in lithium diborate glass // Phys. Chem. Glass. 1994. V. 35. № 1. P. 6–9.

  22. Budewski E., Bostanoff B., Vitanoff T., Stoinoff Z., Kotzewa A., Kaischew R., Zweidimensionale Keimbildung und Ausbreitung von Monoatomaren Schichten versetzungsfreien (100) – Flächen von Silbereinkristallen // Phys. Status Solidi B. 1966. V. 13. № 2. P. 577–588.

  23. Neiman T.S., Yinnon H., Uhlmann D.R. Crystallization kinetics lead metasilicate // J. Non-Crystalline Solids. 1982. V. 48. № 2–3. P. 393–403.

Дополнительные материалы отсутствуют.