Термодинамическая мера неравновесности в процессе роста кристаллов Текст научной статьи по специальности «Физика»

Отклонение от термодинамического равновесия при росте кристаллов, выраженное в форме разности химических потенциалов, называют движущей силой кристаллизации, а как результат отклонения системы от термодинамического равновесия обычно воспринимается масса выросшего за определенное время кристалла. Однако только на основе термодинамических представлений нельзя функционально связать разность химических потенциалов и массовую скорость роста кристалла. Такая задача решается в рамках кинетических теорий роста кристаллов. Получается, что понятие «степень отклонения от равновесия», выраженное вполне точно и строго в рамках классической равновесной термодинамики к росту кристаллов как к кинетическому процессу, может быть применимо со многими оговорками. На практике специально подбираются условия, обеспечивающие технологически значимую скорость роста кристалла при определенной степени отклонения от равновесия, для обозначения которой принято пользоваться другими термодинамическими переменными — пересыщением Эи или переохлаждением ЭТ [1, 2]. При этом бывают случаи, когда некоторые грани кристалла в заведомо неравновесных условиях практически не растут [1]. Для описания такого феномена часто оперируют понятием «мертвая зона» или «метастабильная область». Таким образом, создается определенный парадокс: степень отклонения от термодинамического равновесия принимает не нулевое значение, а неравновесный процесс движения системы к равновесию не развивается. Поэтому в сознании многих исследователей понятие «степень отклонения от равновесия» носит устойчивый качественный оттенок.

Противоречие можно снять, если в качестве «степени неравновесно сти» для заведомо неравновесного процесса принять удельное производство эн-

ТЕРМОДИНАМИЧЕСКАЯ МЕРА НЕРАИНОВЕСНОСТИ В ПРОЦЕССЕ РОСТА КРИСТАЛЛОВ

Д. г.-м. н. В. И. Ракин

гакт @geo. кот 1зе. ги

тропии, включающее не только термодинамическую силу, но и термодинамический поток. Известно, что в неравновесной термодинамике, в так называемом линейном режиме существует единственный принцип экстремума — минимум производства энтропии открытой системой при выходе на стационарное состояние [3]. Значение производства энтропии в стационарном состоянии показывает истинную степень неравновесности такой системы. Даже в нестационарном состоянии системы с помощью величины производства энтропии можно оценивать ее «не-равновесность». В нелинейном режиме также возможны стационарные состояния — диссипативные структуры, в которых производство энтропии системой принимает вполне определенные фиксированные значения. Производство энтропии правильнее рассчитывать на единицу массы или объема системы. Эта величина является вполне определяемой для растущего кристалла и для многих других физико-химических систем. Их можно сравнивать и делать определенные заключения о неравновесности систем. В определенном смысле производство энтропии аналогично понятию «уровень метаболизма» в биологии.

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

1. Диффузия

в одномерной колонке

Рассмотрим гелевую колонку длиной I см, и площадью сечения 8П см2,

контактирующую с обоих концов с системами: слева раствор постоянной концентрации, 0.1 моль/л, справа — чистый растворитель. В результате изотермической диффузии в гелевой колонке со временем устанавливается постоянный градиент концентрации. Коэффициент диффузии принят равным 10-5 см2/с. Результаты моделирования представлены на рис 1.

Рис. 1. Распределение концентрации растворенного компонента в диффузионной колонке в разные моменты времени

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

Локальное производство энтропии для обсуждаемой открытой системы будет определяться выражением

2*2

т ах т ах

где термодинамические силы Б связаны друг с другом уравнением Гиббса— Дюгема

+П2с1|112 = 0.

Потоки растворенного вещества ^ и растворителя 12 также взаимозависимы через уравнения, выражающие неподвижность центра масс,

+ М2Т2 = 0

или +У212 =0, где Мк — молярная масса компонентов, Ук — парциальный молярный объем.

Градиенты химических потенциалов можно вычислить, пользуясь формулой Льюиса

ц = Цо + ИТ 1п а , где а — активность вещества. После дифференцирования и перехода к молярным концентрациям при условии постоянства коэффициентов активности, получим

Фі

дх

кг ап

щ дх

Диффузионный поток растворенного вещества определяется эмпирическим уравнением Фика

дп1

дх

Если количество растворителя в системе постоянно, окончательно удельное производство энтропии системой примет вид

1 а:Э

р =-

_ш'л

1 щ

Юо А М1п1

вался, поскольку на конечный результат он практически не влияет.

2. Рост кристалла из раствора

Рассмотрим открытую систему, включающую вертикально ориентированную поверхность кристалла, площадью 80 с бесконечно тонким слоем кристаллического вещества, и раствор в пограничном слое толщиной 8. Обычно для измерения кинетики роста кристалла добиваются выхода на стационарный режим, когда скорость роста стабилизируется. Одновременно пограничный слой раствора также устанавливается в стационарное состояние, близкое к изображенному на

М2п2 у

Физическая размерность производства энтропии в расчете на 1 литр объема системы — Дж-К-1-с-1-л-1. Результаты расчетов производства энтропии в разные моменты времени численного эксперимента (рис. 1) приведены на рис. 2.

Рис. 2. Зависимость производства энтропии при диффузии от времени процесса

Таким образом, производство энтропии действительно достигает низкого стационарного значения, равного 4.310-5 Дж/(К-с-л). Однако эта величина не оказалась минимальной в ходе установления стационарного состояния. Это несоответствие теореме При-гожина связано с выходом части растворителя из системы, поскольку объем системы был зафиксирован. Выход растворителя из системы достиг максимального значения к моменту времени 104 с и полностью прекратился в стационарном состоянии. Этот фактор в расчетах энтропии не учиты-

РрУ

м ’

-М-} =

[дх,

-Оп'я,

где Б — коэффициент диффузии растворенного вещества в растворе, М —

молярная масса, п§

п5

- гра-

ае8+М = 0.

(2)

Рис. 3. Распределение концентрации в пограничном слое раствора у поверхности растущего кристалла

рис. 3. [4]. При росте кристалла с постоянной скоростью будем перемещать границы выделенной системы вместе с поверхностью растущего кристалла и постоянным количеством растворителя со скоростью роста вглубь раствора. Считаем, что градиент концентрации в пределах выделенного пограничного слоя постоянный. Таким образом, с одной стороны в открытую систему поступает поток растворенного вещества 1Б, а с другой выходит равный ему поток кристаллического вещества 1с. Количество растворителя в пределах выделенной системы остается неизменным. Пусть описываемая система находится в изотермических условиях, и термический эффект от кристаллизации незначителен. Не будем учитывать также тангенциальный поток раствора вдоль грани кристалла, считая его ламинарным и его вклад в изменение энтропии незначительным.

Потоки вещества через единичную площадку, измеряемые в единицах (моль-с-1 см-2), составляют

диент концентрации в пограничном слое, р — плотность кристалла.

Учитывая, что стационарное состояние предполагает сохранение постоянной энтропии, производство энтропии внутри системы должно компенсироваться отрицательным потоком энтропии в систему [3]

<155 Л Л Л

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

в которой, первое слагаемое определяет вклад диффузии, второе — вклад процесса роста кристалла, и = —

молярная скорость роста кристалла, рассматриваемая как «скорость поверхностной реакции», N — число молей, А = (щс -\1С) — химическое сродство «поверхностной реакции», Ц — градиент химического потенциала в растворе.

С другой стороны, поток энтропии в систему определяется выражением

1 <ДЕ ц§ -цс

А

Т А

-^0, (4>

ёЕ

где = Дес1У — полный поток

энергии. Выражения (3) и (4) с учетом уравнения (2) дают два альтернативных способа измерения производства энтропии в стационарном состоянии. Однако проще выполнить расчет по формуле (3).

Очевидно, что в системе выполняется баланс всех потоков вещества

1^ = 1^1 = 1^!. Вклад диффузии в производство энтропии определяется выражением

Л

Б

Рв =-Ьб = т / П'|'.а^

т.

V * * V

Преобразуем его, учитывая одномерную диффузию, постоянный градиент концентрации в рассматриваемой системе и связь активности компонента в растворе с концентрацией

ак = У к-^Пь

где ук — коэффициент активности, пк — молярная концентрация компонента, N2 — общее число молей в объеме V Получим

п5 - п.

Дп

5

5 5

1 да

ая дх

:ЯТ

п5 - п. П.5

ЯТДп

5

П§5 •

Рб =

БЯ80Дп5 | dx =

52

о

П

1п а^ = 1п

ая0

(

1 +

Да

л

Да Дп

а%о пя0 :

где Ап — значение абсолютного пересыщения раствора непосредственно у поверхности кристалла.

Ро = (Ц8С/с) и8о =

ЯДп _ ЯДп БДп 5

=------^о =------------—5 Бо.

пя0 пя0 5

Для производства энтропии при стационарном росте участка грани площадью Бо, прилегающей к объему V = Бо8 пограничного слоя, получим окончательное, относительно простое выражение

= dt =

Р = ■

БИЯ

5

Дп д + ДпДц

п.

П.о

(5)

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

Удельное производство энтропии поверхностью кристалла в расчете на единичный объем системы будет определяться выражением

1 =

"аг=

Р* =-

5д

д

Дп5 + ДпДп5

П.о

БЯБоДп 5 п 5 БИБо . д

=------0—5 1п—- и-------0 Дп д.

5 Пяс Пяс5

Вклад роста кристаллической поверхности в производство энтропии можно выразить формулой

Ро = | Аиаэ И (|бс -|с) иБ0. о Б Т т 0

Б0

Разность химических потенциалов компонента в растворе у поверхности кристалла и в кристалле

1^с(Т) = |4о(Т) + РТ1п aSc,

1с(Т) = |.о (Т) + РТ1па8о .

При данной температуре раствора Т существует равновесный химический потенциал растворенного вещества в растворе (и соответствующая активность а80), который равен химическому потенциалу вещества в кристалле — цс.

Получим

а

|8с(Т) -|с(Т) = РТ1п-Зч

а§0

При низкотемпературном росте кристаллов из раствора обычно а,с > а.0

и азс - а.0 = Аазс << -о тогда

Р* =

ок Д д —-Дп д

БМ Дп

5

ра вдоль грани не рассматривается, результаты расчета производства энтропии по формуле (6) будут незначительно занижены.

На основе экспериментальных данных по росту граней октаэдра кристалла алюмокалиевых квасцов в стационарном режиме [4], рассчитано производство энтропии с ростом пересыщения (рис. 4). В частности, при температуре 20° С и абсолютном пересыщении (2-10-5 моль/л), обычном для условий выращивания кристаллов квасцов, удельное производство энтропии составляет 2-10-8 Дж/(К-с-л). Для сравнения, в режиме стационарной диффузии вещества в растворе при градиенте концентрации 0.1 моль/(л-см) удельное производство энтропии достигает гораздо больших значений — 410-5 Дж/(К-с-л).

(6)

В диффузионном режиме роста обычно полагают, что Ап << Ап§ и соответственно ц.с = п.0.

Тогда выражение (6) упрощается до следующего:

4

3.5 -3

2.5 -. 2

1.5 -1

0,5

0

Пя05"

В кинетическом режиме роста кристалла пересыщение раствора у поверхности кристалла и перепад концентрации в пограничном слое могут принимать разные и независимые друг от друга значения, но в предельном случае Ап << Ап§ и

Р* БР Д Д

Р =-------2 ДпДп 5.

Пя05

Из уравнения баланса вещества получим нормальную скорость роста кристалла

0 0,5 1 1,5 2 2,5 3

ДС, 10"5 моль/л

Рис. 4. Зависимость удельного производства энтропии при росте грани (111) алюмокалиевых квасцов от пересыщения

3. Рост кристалла из расплава

Рассмотрим аналогичную задачу для роста кристалла из расплава методом Чохральского. Метод Чохральско-го предполагает достижение стационарного режима роста кристаллической були плоским фронтом (рис. 5) [2]. Представим однокомпонентную двухфазную систему. Теплоотвод от растущей поверхности кристалла осуществляется через охлаждаемый кристалл [2]. Перед фронтом кристаллизации в

V = - .

Ро 5

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

Учитывая, что в рассматриваемую систему включен не весь пограничный слой, а только его линейный участок, а вклад тангенциального потока раство-

Рис. 5. Распределение температуры в расплаве у поверхности растущего кристалла

п

и

п

бс

а.0

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

(7)

где АН — разность парциальных молярных энтальпий кристалла и расплава, и = — молярная скорость

роста кристалла.

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

а = Г

дх

vTy

I (М'т Мх) ^ (g)

Т 8

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

При отсутствии диффузионного и конвективного переноса вещества в выделенной системе, поток внутренней энергии можно записать в виде

Разность парциальных молярных энергий расплава и кристалла связана с разностью химических потенциалов и разностью парциальных молярных энтропий

Расчет удельного производства энтропии системой производится с помощью нахождения интегралов

1

Р* =

S05 dt

6 5n

1 ?(Дц+ТДз) еТ

udx +

Т дх 1 ГДР- Л

+-г1^ги(к-5 о Т

Если учесть термическую зависимость химических потенциалов и произвести необходимые преобразования

интегралов, удельное производство энтропии сводится к виду

АН АТ

Л~ AHTSC As +----

in'*

Т

, (9)

где А 8 —разность парциальных молярных энтропий кристалла и расплава АТ = То - Т§с.

Если считать, что AT << T0, Tsc

то можно пренебречь одним слагаемым в квадратных скобках и выражение (9) упростить до вида

Г АН

Р =

Pcv

мт(

г

As +

о

(10)

где V — нормальная линейная скорость роста кристалла, Т' — градиент температуры в расплаве перед фронтом роста, рс — плотность кристалла, М — молярная масса.

Таким образом, удельное производство энтропии растущим кристаллом в расчете на единичный объем контактирующего с ним расплава определяется относительно простым выражением. Нормальная скорость роста в первом приближении линейно зависит от скорости теплоотвода и соответственно от градиента температуры, поэтому полученное выражение для удельного производства энтропии закономерно проявляет квадратичную зависимость от нормальной скорости роста, так же как выражение для роста кристалла из раствора в диффузионном режиме. Поскольку анализируемая расплавная термодинамическая система не включает сам кристалл, за исключением его поверхности, то производство энтропии при диффузии тепла в кристалле не вносит вклада в формулу (10).

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

Исследования проводились в рамках программы фундаментальных исследований Президиума РАН за 2006 г. «Разработка методов получения химических веществ и создание новых материалов», также поддержаны грантом Президента Российской Федерации НШ.2250.2003.5.

Литература

1. Выращивание кристаллов из растворов / Т. [Г. Шетров, СЕ. \Б. [Трейвус, Ю. [Q. Шунин, [А. Ш. [Касаткин. Л.: Недра, 1983. 200 с. 2. ВилъкеК. [Т. Выращивание кристаллов. Л.: Недра, 1977. 423 с. 3. При-гожинИ., КондепудиЩ. Современная термодинамика. От тепловых двигателей до диссипативных структур. М.: Мир, 2002. 461 с. 4. Ракин[В.Ш. Пространственные неоднородности в кристаллообразующей системе. Екатеринбург: УрО РАН, 2003. 370 с. 5. РакинВ. [И. Производство энтропии растущим кристаллом, как мера неравновесности в кристаллогенезисе // Электронный журнал «Исследовано в России», 89, 855—861, 2006. http:// zhurnal.ape.relarn.ru/articles/2006/089.pdf