CC BY
66
Теоретическое исследование распада в бинарных системах.
Часть I. Уравнения спинодального и бинодального распада бинарной
системы в изотропной фазе
12 1 Е.С. Ларин , О.В. Наскалова , Е.М. Кузнецова
1 Южный федеральный университет, г. Ростов-на-Дону 2Чеченский государственный университет, г. Грозный
Аннотация: В работе проведено теоретическое исследование спинодального и бинодального распада бинарных растворов. Так как при спинодальном распаде отсутствует термодинамический барьер, фазовое разложение определяется только диффузией. Таким образом, при распаде многие особенности фазового разложения могут быть описаны аналитическим решением. Получены уравнения спинодального и бинодального распада в растворах без структурных фазовых переходов. Ключевые слова: бинарные растворы, спинодальный распад, бинодальный распад.
Введение
При резком изменении внешних параметров на термостате система переходит первоначально в неравновесное состояние. Процесс релаксации к новому состоянию равновесия определяется различными факторами, например, если несколько различных фаз при данных условиях локально стабильны, или линии фазовых переходов первого рода пересекаются [1]. Исследование фазовых переходов и их кинетика интенсивно изучаются в течение длительного периода и являются в настоящее время классической темой термодинамической (неравновесной) физики.
Естественно, кинетика фазовых переходов является определяющей в понимании превращений в широком круге материалов, таких как жидкие смеси или металлические сплавы со структурными, электрическими или магнитными фазами в твердых телах, а также в жидких кристаллах. В настоящей работе мы рассмотрим теорию распада в изотропных бинарных системах, в которых не происходят фазовые переходы.
На диаграммах состояния бинарных жидких смесей возможно существование промежутка смешиваемости, в котором гомогенное состояние
становится метастабильным или нестабильным. Если начальная фаза находится в нестабильном равновесии, то новая равновесная фаза спонтанно возникает при динамической релаксации. Это процесс, известный как спинодальный распад (СР), при котором бесконечно малые флюктуации приводят к росту областей сосуществования конечных устойчивых фаз. Теоретическое описание СР изотропной смеси было развито J.W. Cahn и J.E. Hilliard [2, 3]. Спинодальный распад - это механизм, при котором раствор двух или более компонентов распадается на фазы с отчетливо различными химическими составами и физическими свойствами. Этот механизм отличается от классического образования зародышей новой фазы (nucleation), фазовое разделение за счет СР происходит намного более определенно и однородно во всем объеме материала, при котором отсутствуют дискретные области образования зародышей. СР интересен, прежде всего, своей относительной простотой, и представляет собой одно из возможных фазовых преобразований в твердых телах, жидких кристаллах, полимерах для которых возможна количественная теория.
С другой стороны, неравновесные метастабильные состояния, такие как гели, могут сформироваться через динамическую фиксацию СР [4-6]. Отметим, что СР и диффузионное образование зародышей фаз и их рост -это два различных механизма, благодаря которым гомогенная пересыщенная фаза может распадаться на состояние с более низкой свободной энергией, состоящей из гетерогенной смеси двух или больше фаз.
1. Идеальные растворы
Раствор - это однофазная система, состоящая из двух или большего числа различных веществ [1]. Возможность образования раствора при заданных значениях Р и Т определяется изменением термодинамического потенциала. Растворение возможно и протекает самопроизвольно, если
X > Ф, где Ф - термодинамический потенциал (свободная энергия Гиббса),
г
X Ф* - сумма термодинамических потенциалов Фг чистых веществ, при
г
смешении которых образуется раствор.
Если привести в соприкосновение два вещества, то можно получить, либо их механическую смесь, либо новое соединение, образующееся за счет их химического взаимодействия. Растворение представляет собой промежуточный процесс между приведенными случаями. Состав растворов в определенных интервалах изменения температуры и концентраций может меняться непрерывно. Понятие растворителя и растворенного вещества носит относительно условный характер. Выбор одного из компонентов раствора в качестве растворителя часто определяется его более высокой концентрацией. Иногда за растворитель принимают вещество, имеющее более высокую температуру кристаллизации.
Идеальным называется раствор, образование которого из компонентов не сопровождается тепловым эффектом и изменением объема: АН = 0, А V = 0. Необходимо отметить, что в действительности незначительные тепловые и объемные эффекты присутствуют всегда. В реальной ситуации можно говорить только о том, в какой степени раствор является идеальным. В случае, когда тепловые и объемные эффекты пренебрежимо малы, при описании свойств такого раствора можно его рассматривать как идеальный.
Идеальные растворы близки к идеальным газам, т. е. системам, получаемым простым механическим перемешиванием компонентов. В этом случае силами взаимодействия между частицами разных веществ (А и В) можно пренебречь по сравнению с силами взаимодействия между частицами каждого из компонентов (А и А или В и В). В качестве идеальных растворов можно рассматривать растворы веществ, близких по своему строению и свойствам.
Изменение энтропии бинарной системы идеального раствора определяется выражением: AS = -^[х1п х + (1 - х)1п(1 - х)].
Концепция идеального раствора оказывается полезной для изучения поведения реальных растворов в тех случаях, когда недостаток опытных данных вызывает необходимость использования приближенных значений и, в особенности, в случае разбавленных растворов, поведение которых может быть достаточно точно описано на основе гипотезы идеального раствора.
2. Модель регулярных бинарных растворов Реальные растворы не обладают свойствами идеальных и для них не соблюдаются законы Рауля. На рис. 1 и 2 представлены зависимости давления пара (р) от состава для двух реальных растворов. Прямые штриховые линии соответствуют давлению пара, которое наблюдалось, если бы растворы были идеальными. На рис. 1 показаны положительные отклонения, при которых давление пара выше, чем для идеального раствора такой же концентрации, а на рис. 2 - отрицательные отклонения, где давление пара меньше, чем для идеального раствора. Касательные к кривым указывают области концентраций, где выполняются законы Генри и Рауля.
Наличие положительных отклонений показывает, что образование раствора затруднено по сравнению с идеальным. При больших положительных отклонениях жидкости расслаиваются. Отрицательные отклонения указывают на более сильную тенденцию к образованию раствора, чем в случае идеального раствора.
Учет взаимодействия между частицами приводит к существенным и очень важным на практике изменениям свойств растворов. Простейшей моделью, учитывающей такое взаимодействие, является модель регулярных растворов, в которой предполагается, что конфигурационная часть термодинамического потенциала (свободной энергии) является квадратичной функцией состава раствора:
П 1 -1
ф = Ефс, + ЕЕФкС,ск + ЯТ£с 1пс . (1)
1=1 1=1 к=1 1=1
В (1) Ф1 - набор потенциалов фаз, все частицы которых однотипны, а Фк являются параметрами взаимодействия компонентов 1 и к в регулярном растворе, зависящими только от внешних условий на термостате в частности, от температуры Т.
р |
О Молярная доля хг 1
Рис. 1. - Положительные отклонения Рис. 2. - Отрицательные отклонения от закона Рауля от закона Рауля
Равновесное состояние твердого раствора определяется из условия минимума его свободной энергии. Гетерофазное состояние твердого раствора возникает в случае, если свободная энергия смеси двух фаз принимает более низкое значение, чем свободная энергия однофазного состояния. Концентрационная зависимость свободной энергии двухфазной смеси идеального раствора описывается линейным законом. Концентрационная зависимость термодинамического потенциала в случае регулярного раствора Ф(с) (1) носит немонотонный характер. В этом случае на зависимости Ф(с) возможны экстремумы. Тогда равновесные составы фаз определяются как абсциссы точек касания общей касательной к кривой Ф(с). Это соответствует равенству химических потенциалов двух фаз:
дФ
М1 = ^2; д =
дс
(2)
Из уравнений (2) получаем уравнение, определяющее условие бинодального
дФ дФ
распада в бинарной системе: — =-.
дс1 дс2
При этом на концентрационной зависимости Ф(с) существуют интервалы концентрации, при которых функция Ф(с) является выпуклой кривой (рис. 36). Равновесное состояние системы, описываемое выпуклой кривой на концентрационной зависимости термодинамического потенциала, всегда будет гетерофазным.
Рис. 3. - Фазовая диаграмма Т—с диаграмма (а) и зависимость свободной энергии от концентрации с при температуре Т2 < (б)
Процесс распада идет при непрерывном понижении свободной энергии и поэтому не требует активационного образования зародышей новой фазы. Условие СР определяется уравнением (точки 51 и Б2 на рис. 3а):
д 2Ф(Т, с) дс 2
= 0
(3)
:
Положение критической точки (Тк, ск) СР находится из уравнения (3) и
д 3Ф(Т, с)
условия: дс3
= 0
При температуре Тк для состава ск термодинамический потенциал обладает особенностью в результате того, что происходит исчезновение вогнутости кривой Ф(с), обуславливающее наличие расслоения. В этой точке кривые бинодали и спинодали касаются.
Для модели регулярного раствора бинарной системы имеем: Ф = Ф1С1 + Ф2С2 + Ф12С1С2 + Ш(С1\пС1 + С21пС2), С1 + С2 = 1. В этом случае химические потенциалы определяются уравнениями:
М1 = Ф1 + Ф12(1 - С1) + ЯТ(Ъс1 +1) М2 = Ф2 + Ф12С1 + ят[1п(1 - С1) +1)]. (4)
Как следует из уравнений (4), химические потенциалы смешиваемых веществ зависят от состава немонотонно. Рассмотрение распада в модели регулярных растворов носит сугубо качественный характер.
3. Уравнение состояния Ван-дер-Ваальса Первым уравнением состояния, описывающим поведение реальной жидкости, является уравнение состояния Ван-дер-Ваальса, предложенное в 1873 г. [7]. Ван-дер-Ваальс преобразовал простейшее уравнение состояния для идеального газа: р¥ = пЫкТ, где N - число Авогадро, к - постоянная
Больцмана (Я = kN, идеальная газовая постоянная), п - число молей. Он
2 2
заменил V на V - Ш и р на р + N а/V . Первая модификация связана с тем, что доступный объём для каждой молекулы должен быть уменьшен на объем, уже занятый N1 - 1 = N2 молекулами, существующими в объеме V (метод исключения объема). Вторая модификация связана с дальнодействующими силами притяжения между молекулами. Дополнительный вклад в давление определяется некомпенсированными
силами притяжения молекул, расположенных на поверхности жидкости. Тогда уравнение Ван-дер-Ваальса принимает следующий вид:
P = kT/(v -b) - a/v2, v = V/N. (5)
3 2
Уравнение (5) можно переписать в виде: pv - (pb + kT)v + av - ab = 0. Это уравнение может иметь при данном давлении р одно или три решения. Его можно преобразовать: pc(v - vc)3 = 0, где vc = 3b; pc = a/Hb1 kTc = 8a/27b.
Другим уравнением состояния «модельной» системы является уравнение состояния Бертло: P = kT/(v - b) - a/v2T.
Другое развитие теории Ван дер Ваальса для бинарной смеси было
предложено в [8]: P = kT/(v -b(x)) - a(x)/v , где x - концентрация одного из
22 веществ в смеси, a(x) = (1 - x) a11 + x(1- x)a12 + x a22.
4. Модель Флори - Хаггинса
Наиболее распространенной моделью свободной энергии смешивания изотропной бинарной смеси жидких кристаллов, суспензий анизотропных наночастиц (нанотрубки, графены) в жидких кристаллах, полимеров является теория Флори - Хаггинса (Flory - Huggins) [9]. Свободная энергия в этом случае имеет вид:
F = F(Ni,N2;P,T) = Nf (^ ,р2;P,T) = N
г \
<Рх л Фгл
— 1П0>1 +— ln^2
V ni n2 У
(6)
п. N. п2 N 2
где N = п.N. + п2N2, р. = N (2 = ^ ' ^ = 1' Пх - числ0 мест,
занимаемых одной жидкокристаллической молекулой; п2 - число мест, занимаемых одной полимерной цепью или асимметричной наночастицей; ф. и ф2 - объёмные фракции компонентов 1 и 2 соответственно [10]. Параметр взаимодействия в модели Флори - Хаггинса х является функцией только температуры: х = А + В/Т, где А и В - постоянные [10].
Химические потенциалы определяются выражениями:
:
У
У
ч ^ У
V 1 У N 2
ч ^ 2 У N.
п
п„
/ + N
д/ д?1 дф1 дN 1
п
/ + (1 -ф)
дфх
/ + ф
дф
(7)
2
В случае модели Флори - Хаггинса (6) имеем:
д/ = — - — +—1ПФ1 -—1П(1 -ф )-х( - 2ф);
- 2х;
дф1 п1 п2 п1 п2
д2/ 1 1 _— =__\--
дф2 пф п2(1 -ф)
д3 / =__1_
3 2
дф1 п1ф1 п
\
1
:(1 -ф )2
Тогда химические потенциалы компонентов имеют вид:
(1 -ф1)+ щх(1 -Ф1)2
У = 1пФ1 +
/ л
1 - п.
п2 У
/
У
1П(1 -ф)
+
п
1 —2
п
ф + п2ХФ\
1 У
На рис. 4 изображена зависимость свободной энергии в модели (6). Точки касания двойной общей касательной определяют равенство химических потенциалов двух фаз.
Как известно, условие равенства химических потенциалов определяет бинодаль:
УМф!) = ^1(^1)
Ла2(фа1) = УР2(фР1) Равенство химических потенциалов геометрически определяется
двойной общей касательной на зависимости ^(с) (рис. 4) [9]:
/(фа1)-/(фв1) = (фа1 -фв1) д/
дф1
д/
дФ1
д/
ф 1
Теперь найдем уравнение СР в модели Флори - Хаггинса, которое определяется равенством нулю полной производной химического потенциала
d^1
d(1
0
Рис. 4. - Концентрационная зависимость свободной энергии (6) и вид
двойной общей касательной
д¡лх д2 / Так как —— = п (1 - ()--, то уравнение СР имеет вид:
оЫх д(
д 2 / = + 1 _ 2 П 2 (1 _( )+ п( _ 2П1 П (1 _( ) = 0 д(2 п( п2 (1 _() п2 (1 _( )п(
п2 (1 _ () + п( _ 2пхпм (1 _ () = 0
2хпхп(( _ ((2хпп + п2 _ п ) + п2 = 0 ;(1 _()+ п(
Х =
2п1п2(1 (1 _( )
Система уравнений (7) определяет критическую точку на куполе СР:
-1
X
2пгп2 , Л/п2
,(„ = г-" V. (8)
21
'п1 + >/ п2
(( + 4^2 )2
Из (8) видно, что в общем случае в рамках модели Флори - Хаггинса спинодаль и, соответственно, бинодаль асимметричны на фазовой Т-ф диаграмме (рис. 5б). При п1 = п2 = 1 эта модель аналогична модели регулярных растворов (3).
Т
ф|
а)
Ъ
Рис. 5. - Зависимости и Т-ф1 фазовая диаграмма модели Флори -Хаггинса: а) прямые соответствуют двойным общим касательным при различных температурах Т. Значения Т и ф1 в точках касания определяют бинодальную кривую на плоскости Т-ф 1; б) Т-ф1 фазовая диаграмма модели Флори - Хаггинса (6) при п1 = 3, п2 = 2, В = 5, А = 0
Состояния системы выше кривой бинодали (1) соответствуют устойчивому однофазному состоянию (рис. 3а) при Т1, фазы между кривыми 1 и 2 также устойчивы (метастабильные), а ниже спинодади (2) - область несмешиваемости веществ. В этой области, по крайней мере, сосуществует две устойчивые фазы.
Выводы и заключение
Как было отмечено в [11], наблюдаемые в бинарных смесях (сплавах) процессы распада можно разбить на два различных типа. Особенностью распада первого типа является тот факт, что он происходит за счет взаимодействия атомов компонент смеси, приводящего при понижении температуры к стабильности механической смеси однородных фаз. Такой тип распада бинарной смеси называется собственным распадом. Второй тип распада определяется наличием в смесях (твердых растворах) структурных
фазовых переходов с образованием различных упорядоченных фаз. Такой тип распада, связанного с сосуществованием фаз разной симметрии называется несобственным.
Рассмотренный нами в данной работе термодинамический потенциал Ф описывает состояние системы в изотропной фазе, который описывает только собственный распад. Однако взаимодействие частиц бинарной системы может приводить к различным типам упорядочения. В части II нами будут исследованы различные модели распада методами теории катастроф. Как известно, в различных бинарных системах, таких как твердые растворы, смеси жидких кристаллов, полимеров, возможны и структурные фазовые переходы [12-15]. Такие переходы могут описываться одним или несколькими параметрами порядка [16]. Например, в случае жидких кристаллов, молекулы (стержневидные, дисковидные) за счет конкурирующих взаимодействий ассиметричного ван-дер-ваальсовского притяжения и сил стерического отталкивания (Онзагер) может приводить к нематическому упорядочению [17, 18]. При фазовых переходах первого рода возможны области сосуществования двух и более фаз. В этом случае также возможен распад смеси на фазы. При этом возможна реализация при изменении внешних условий (температура, давление) ситуации, при которой наблюдается конкуренция двух типов распада.
СР представляет большой интерес с практической точки зрения, например, при получении точной дисперсной микроструктуры, которая значительно может усиливать различные физические свойства материала. Существует много примеров материалов, в которых требуется контролированная и строго определенная фазовая сегрегация. Например, при производстве пластмассовых солнечных батареей фуллерены С60 и полимерная фаза должны быть хорошо разделены, создавая микро- и даже наномасштабную сегрегацию фаз с п- и ^-типом проводимости. В этом
случае удается существенно увеличить n-р контакт и таким образом обеспечить более эффективное для разделения заряда, вызванного действием света [19].
Ларин Е.С. благодарит за финансовую поддержку НИИ Физики ЮФУ (внутренний грант № БЧ0110-11/2017-42 (в рамках госзадания 3.6386.2017/БЧ).
Литература
1. Gibbs J.W. The Collected works. V. 1. New Haven: Yale University Press, 1948. p.105.
2. Cahn J.W., Hilliard J.E. Free Energy of a Nonuniform System. I. Interfacial Free Energy // J. Chem. Phys. 1958. V. 28. pp. 258-267.
3. Cahn J.W. Phase Separation by Spinodal Decomposition in Isotropic Systems // J. Chem. Phys. 1965. V. 42. pp. 93-99.
4. Zaccarelli E. Colloidal gels: equilibrium and non-equilibrium routes // J. Phys. Condens. Matter. 2007. V. 19. pp. 323101.
5. Wittkowski R. et al. Scalar ф field theory for active-particle phase separation // Nat. Commun. 2014. V. 5. pp. 4351.
6. Speck T., Bialke J., Menzel A.M. and Lowen H. Effective Cahn-Hilliard equation for the phase separation of active brownian particles // Phys. Rev. Lett. 2014. V. 112. pp. 218304.
7. Van der Waals J.D. Over de Continuiteit van den Gas- en Vloeistoftoestand. University of Leiden, 1873. р. 55.
8. Rigby M. and Alder B.J. Augmented van der Waals Theory for Mixtures // J. Chem. Phys. 1970. V. 52. pp. 3665-3670.
9. Flory P.J. Principles of Polymer Chemistry. Ithaca, NY: Cornell University Press, 1953. 688 p.
10. Olabisi O., Robeson L.M. and Shaw M.T. Polymer Polymer-Miscibility. New York: Academic, 1979. 370 p.
11. Гуфан А. Ю. Теория фазового перехода типа собственного распада бинарного твердого раствора // Кристаллография. 2004. том 49, № 3. с. 515523.
12. Shen C. and Kyu T. Spinodals in a polymer dispersed liquid crystal // J. Chem. Phys. 1995. V. 102(1). pp. 556-562.
13. Иванов В.В., Наскалова О.В., Ларин Е.С., Зубхаджиев М.-А.В. Исследование покрытий с антифрикционными свойствами на основе гидрооксидной пленки А1(0Н)3 на алюминии с нанесённым твердосмазочным покрытием MoS2 // Инженерный вестник Дона, 2017, №4. URL:ivdon.ru/ru/magazine/archive/n4y2017/4548.
14. Smith N. and Provatas N. Generalization of the binary structural phase field crystal model // Phys. Rev. Materials. 2017. V. 1. pp. 053407.
15. Gruhn T., Pogorelov E., Seiferling F. and Emmerich H. Analyzing spinodal decomposition of an anisotropic fluid mixture // J. Phys.: Condens. Matter. 2017. V. 29. pp. 055103.
16. Гуфан Ю.М. Структурные фазовые переходы. М.: Наука, 1982. 304 с.
17. Prostakov A.E., Larin E.S. and Stryukov M.B. Phenomenological Theory of Phase Transition from Isotropic to Uniaxial, Biaxial, and Reentrant Nematic Phases in Achiral Liquid Crystals // Crystallography Reports. 2002. V. 47(6). pp. 1041-1048.
18. Larin, E. S. Phase Diagram of Transitions from an Isotropic Phase to Nematic and Smectic (Uniaxial, Biaxial) Phases in Liquid Crystals with Achiral Molecules // Physics of the Solid State. 2004. V. 46(8). pp. 1560-1568.
19. Oono Y. and Puri S. Study of phase-separation dynamics by use of cell dynamical systems. I. Modeling // Phys. Rev. A. 1988. V. 38(1). pp. 434-453.
References
1. Gibbs J.W. The Collected works. V. 1. New Haven: Yale University Press, 1948. p. 105.
2. Cahn J.W., Hilliard J.E. J. Chem. Phys. 1958. V. 28. pp. 258-267.
3. Cahn J.W. J. Chem. Phys. 1965. V. 42. pp. 93-99.
4. Zaccarelli E. J. Phys. Condens. Matter. 2007. V. 19. pp. 323101.
5. Wittkowski R. et al. Nat. Commun. 2014. V. 5. pp. 4351.
6. Speck T., Bialke J., Menzel A.M. and Lowen H. Phys. Rev. Lett. 2014. V. 112.pp.218304.
7. Van der Waals J.D. Over de Continuiteit van den Gas- en Vloeistoftoestand. University of Leiden, 1873. p. 55.
8. Rigby M. and Alder B.J. J. Chem. Phys. 1970. V. 52. pp. 3665-3670.
9. Flory P.J. Principles of Polymer Chemistry. Ithaca, NY: Cornell University Press, 1953. 688 p.
10. Olabisi O., Robeson L.M. and Shaw M.T. Polymer Polymer-Miscibility. New York: Academic, 1979. 370 p.
11. Gufan A. Yu. Kristallografiya. 2004. t. 49, № 3. pp. 515-523.
12. Shen C. and Kyu T. J. Chem. Phys. 1995. V. 102(1) pp. 556-562.
13. Ivanov V.V., Naskalova O.V., Larin E.S., Zubkhadzhiev M.-A.V. Inzenernyj vestnik Dona (Rus), 2017, №4 URL:ivdon.ru/ru/magazine/archive/n4y2017/4548.
14. Smith N. and Provatas N. Phys. Rev. Materials. 2017. V. 1. pp. 053407.
15. Gruhn T., Pogorelov E., Seiferling F. and Emmerich H. J. Phys.: Condens. Matter. 2017. V. 29. pp. 055103.
16. Gufan Yu.M. Strukturnye fazovye perekhody [Structural phase transitions]. M.: Nauka, 1982. 304 p.
17. Prostakov A.E., Larin E.S. and Stryukov M.B. Crystallography Reports. 2002. V. 47(6). pp. 1041-1048.
18. Larin E. S. Physics of the Solid State. 2004. V. 46(8). pp. 1560-1568.
19. Oono Y. and Puri S. Phys. Rev. A. 1988. V. 38(1). pp. 434-453.
