В термодинамике гетерогенных систем Гиббса все экстенсивные величины — энергия, энтропия, масса компонентов — определяются суммированием членов, пропорциональных объемам фазы и площадям поверхностей раздела. Это верно при условии, что ни в одном месте системы переходные слои, принадлежащие разным поверхностям раздела или разным участкам поверхности раздела (рис. 1) не перекрываются. Перекрытие будет отсутствовать тогда и только тогда, когда из одной поверхности раздела фаз нельзя достигнуть другой, не пересекая областей, гомогенных по интенсивным свойствам, Т. е. обладающих свойствами фазы. Именно поэтому Гиббс рассматривал только такие тонкие пленки, средние части которых обладают свойствами фазы.
Между тем в учении о коллоидах и дисперсных системах особое значение имеют ситуации, когда при сближении двух частиц возникает перекрытие их поверхностных зон.
Рассмотрим, какие эффекты, новые по сравнению с трактованными в труде Гиббса и его прямых последователей, появляются при таком перекрытии в случае тонкой плоскопараллельной прослойки жидкости (или газа). Прежде всего, при таком перекрытии, если оно наступает в изотермическом процессе
Рис. 1. Перекрытие поверхностных зон. |
Сближения объемов двух фаз, как правило, должна меняться свободная энергия Гиббса системы (свободная энтальпия), а следовательно, должны возникать сильные взаимодействия, зависящие от степени перекрытия. В том случае, когда наступает перекрытие в области плоскопараллельной прослойки, взаимодействие однозначно характеризуется изменением давления в ней по сравнению с давлением в той объемной фазе, частью которой является прослойка. Это избыточное давление было нами названо расклинивающим давлением [1]. Здесь, однако, требуется уточнение, поскольку в случае перекрытия исчезают части фазы, сохранявшие изотропность свойств, в том числе изотропный характер гидростатического давления. Давление в зоне перекрытия анизотропно и характеризуется тензором давления. В общем случае этот тензор меняет свои компоненты в зависимости от местоположения точки, от ее координаты, направленной поперек прослойки.
Существенное упрощение в трактовке равновесия’таких прослоек достигается, если исключить из рассмотрения силы, зависящие от расстояния, заменив их эквивалентными им силами близкодействия. Сведение сил электростатического взаимодействия, связанных с объемными зарядами ионного происхождения, к силам близкодействия производится с помощью тензора натяжений электростатического поля Максвелла. Заменить силами близкодействия силы молекулярного притяжения можно, пользуясь теорией, развитой Лифшйцём [2] и сводящей эти силы к флуктуациям электромагнитного поля. Включая и соответствующий этим последним тензор в тензор давления, мы исключаем силы дальнодействия, приложенные к элементам объема прослойки.
Из условия равновесия этих последних вытекает уравнение
Pzz — Const
Где Ргг — компонента тензора давлений, нормальная к поверхностям, ограничивающим прослойку (жидкости или газа) [3, 4]. Тогда для расклинивающего давления П (К) мы получаем определение
П (Ft)=PZ2—P0 (1)
Где Р0 — давление в объемной фазе, утоньшением которой получился слой толщины H.
Расклинивающее давление ЩИ) есть основная термодинамическая характеристика тонких прослоек, необходимая для развития термодинамики гетерогенных систем Гиббса в новом направлении.
Известно, что для конкретизации значения термодинамических функций, определяющих равновесие гетерогенных систем, необходимо привлечение статистических методов. Ими следует пользоваться, например, для вывода изотермы адсорбции или для вывода зависимости межфазного натяжения от концентрации компонентов. Аналогично возникает задача использования статистических методов для расчета расклинивающего давления в функции толщины межфазной прослойки. Подход к решению этой задачи зависит от природы тех поверхностных сил, которые определяют эффект перекрытия межфазных переходных слоев. В общем случае это весьма сложная задача. Однако возможность более простых решений ее по сравнению со случаем классических гиббсовских поверхностных явлений возникает тогда, когда толщина пленки существенно превышает размеры молекул и ее состояние при перекрытии может трактоваться в основном макроскопически. Примером могут служить расчеты дисперсионной слагающей расклинивающего давления на основе макроскопического подхода Лифшица. Другой пример — расчет электростатической слагающей, возникающей при перекрытии диффузных ионных атмосфер, или же при перекрытии диффузных адсорбционных слоев в бинарных растворах неэлектролитов [5].
(2) |
Остановимся на первом примере для того, чтобы проанализировать, какую роль в этих расчетах играют статистические и макроскопические методы. Для достижения максимальной ясности и простоты воспользуемся уравнением Гиббса—Дюгема, дополнив выражение для работы двумя членами
— Sn dh — f S (г|>! Dax — f г|>2 Da2)
Ip2
H
Ju
Рис. 2. К выводу электростатической слагающей раскли* нивающего давления.
Где 5 — площадь прослойки между поверхностями (прослойка предполагается плоскопараллельной, толщины И);
П — расклинивающее давление прослойки; г|)1>2 — потенциал одной из двух межфазных поверхностей (рис. 2) относительно бесконечно удаленного электрода;
А1)2 — поверхностная плотность заряда.
Первый член в (2) — обычное выражение работы, затрачиваемой против расклинивающего давления при уменьшении толщины прослойки. Второй член выражает работу переноса зарядов Dah2 из бесконечности на поверхности пластин. Вычитание из обоих частей дополненного членами (2) уравнения Гиббса—Дюгема суммы 5 (i^i + Дает нам выражение для дифференциала дополненного термодинамического потенциала Гиббса:
DGj Р = — Sll dh — S (Аг + Ъ (3)
Поскольку DG‘R, Р есть полный дифференциал, то можно написать:
/jn =/_ао
Ч> 1 |
V dh /ф,.1|>,
Интегрирование уравнения Пуассона—Больцмана для прослойки электролита дает:
Здесь потенциал |э есть функция координаты х, отсчитываемой от левой поверхности прослойки, по нормали к ней
8я 8я
А = Т" E°LZl = ~тЕс%г*
Где е — заряд электрона;
Е — диэлектрическая проницаемость; Сх и с2 — концентрация катионов и анионов; гх и z2 — их заряды;
B ~ E/KT, где K — постоянная Больцмана, а Т—температура по шкале Кельвина; С — постоянная интегрирования. Для х = О
Откуда следует:
А п — а Г EZlb^ — 1 . _ I -. , ^ 2
L——- Тх——— 1- —Тш— J ~ ~дГJ,=O (8)
При этом очевидно, что при H —* оо:
TOC o "1-3" h z С = Соо = 0 (8а)
Из уравнений (7)—(8) следует:
В2 а Г е2^1 , e~z*b
О»
П Г Р— zjoipi -1
—CI (9)
— ■ ( — B‘fl ( Дс nm
СТЧа H 32я*Ь V Dh Л,,^ ( ‘ С помощью общеизвестного тождества
И очевидного (см. рис. 2) соотношения
(Ж) = JSL (12)
V Dh }с,|>, е V ‘ Мы из уравнения (10) получим:
В A / Ас 8п ь |
■)» <13) |
= Dh /ъ |
Из выражения (13) и уравнений (4) и (8), интегрируя и учитывая, что при H —> оо справедливо (8а), получим хорошо известное выражение:
П = JL— с = W [Cl — 1 ) + с2 — 1)] — JL— Е (14)
Первый член выражает избыточное гидростатическое давление, второй — электростатическое натяжение Максвелла,