Новые корреляции поверхностного натяжения с объемными свойствами жидкостей Текст научной статьи по специальности «Физика»

КОЛЛОИДНАЯ ХИМИЯ

УДК 541.18

НОВЫЕ КОРРЕЛЯЦИИ ПОВЕРХНОСТНОГО НАТЯЖЕНИЯ С ОБЪЕМНЫМИ СВОЙСТВАМИ ЖИДКОСТЕЙ

Б. Д. Сумм

(кафедра коллоидной химии)

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

Систематические исследования поверхностного натяжения начались в середине Х1Х в. К этому времени была создана теория капиллярных явлений (Лаплас, Юнг, Гаусс, Пуассон, Кельвин) и были разработаны достаточно точные методы измерения поверхностного натяжения (капиллярный подъем, сталагмометрия, пластина Вильгельми). Вскоре был накоплен достаточно большой объем информации о поверхностном натяжении воды и органических веществ, в том числе о температурной зависимости (Менделеев). Эти сведения стимулировали поиск эмпирических соотношений для связи поверхностного натяжения жидкостей о с другими физическими свойствами.

Первую корреляцию поверхностного натяжения (1886) предложил Стефан:

о = О /^1/3 V 2/3),

- мольная теплота испарения, N - число Авогадро, Ут - атомный объем [1].

В 1893 г. Рамзай и Шилдс описали температурную зависимость о(Т) многих веществ уравнением

о (ц / р)2/3 = к (Тк - Т - 6),

где ц и р - молекулярная (атомная) масса и плотность жидкости, Т - критическая температура, к = (2.1-2.2). 10-Дж/К [1, 2].К

Метод эмпирических корреляций успешно применяли и позднее. Приведем несколько известных примеров.

Правило Антонова (1907 г.)

где о12 - межфазное натяжение на границе жидкостей 1 и 2 с собственными натяжениями (на границе с паром) о1 и о2 [2]. Исследования смачивания во многих системах показали, что правило Антонова в форме ок1 = |о8 - о1 I

хорошо выполняется и для межфазной поверхности жидкость - твердое тело, ок1 - межфазное натяжение, 08 и о1 - натяжения твердого тела и жидкости на границе с газовой фазой [3].

Уравнение Бачинского (1922 г.)

о = В (р - ру)4, где ру - плотность пара, В - константа [1].

Уравнение Сагдена (1924 г.)

Р = М В1'4 = М о1/4/ (р - ру),

где параметр Р с размерностью объема называется пара-хором. Для многих веществ он не зависит от температуры [1].

Правило уравнивания полярностей Ребиндера (1926 г.)

В соответствии с этим правилом поверхностное натяжение о12 на границе веществ 1 и 2 с симметричными молекулами является функцией их диэлектрической проница-емостей е1 и е2. Зависимость о12 = _/(е1, е2) спрямляется в координатах о = / (А П):

АП=-

1-1 1 + 2

2 -1 2 + 2

где ДП - разность полярностей [3].

Уравнение Кунина (1951 г.) для жидких металлов о = 444.5 ф /г 2,

где ф - работа выхода электрона, г - атомный радиус [1]. Приведенные корреляции носят ограниченный харак-

о12 = о1 - о2

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

По мере развития теории поверхностных явлений и более совершенных экспериментальных методов выдвигаются и новые корреляции поверхностного натяжения. В данной статье представлены корреляции, в основе которых -прямые подтверждения квазикристаллической структуры поверхностных слоев различных жидкостей [4, 5] и сеточная модель поверхностного слоя [6, 7].

Корреляция с теплотой плавления

В теории капиллярности Лапласа поверхностный слой (ПС) жидкости рассматривается как тонкая пленка, отличающаяся по своим механическим свойствам от объемной жидкости. Главная особенность ПС заключается в его способности выдерживать большие упругие (обратимые) деформации и в действии (в пределах ПС) сил поверхностного натяжения о. Модель упругого ПС была использована Лапласом при выводе закона капиллярного давления. Он применил для этого разработанный им ранее приближенный метод расчета напряжений в тонких оболочках под действием внутреннего или внешнего давления.

Сходную модель ПС использовал и Юнг при выводе закона смачивания: он моделировал ПС жидкости натянутой пленкой пренебрежимо малой толщины.

В термодинамике поверхностных явлений агрегатное состояние ПС жидкости на границе с собственным паром не обсуждается, хотя представление о ПС как о самостоятельной фазе представляет краеугольный камень теории Гиббса [1]. Термодинамика дает более общее описание поверхностных явлений, чем механическая теория капиллярности Лапласа-Юнга. По-видимому, по этой причине проблема агрегатного состояния ПС на длительное время оказалась в тени. Принципиальное значение вопросов и понятий, связанных с микроскопической структурой жидкости, однако, весьма велико [1].

До недавнего времени из-за отсутствия адекватных экспериментальных методов выяснение строения ПС жидкостей не проводилось. Только в последние 5-7 лет усовершенствованный метод малоуглового рассеяния позволил получить достоверную информацию о строении ПС различных жидкостей (ртуть, галлий, предельные углеводороды, начиная с С16Н34, некоторые другие). Главный результат: для всех изученных веществ при температурах на несколько градусов выше точки плавления Тт ПС представляет двумерную квазикристаллическую структуру. Ее толщина 5 оценивается в 1-2 молекулярных (атомных) диаметра; параметры двумерного кристалла отличаются от решетки «объемного» кристалла того же вещества [4, 5].

Таким образом, модель твердого ПС можно считать достаточно обоснованной как теоретически (представления Лапласа и Юнга), так и экспериментально.

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

о = 1 @ р 5 = , (1)

т ^т т ¿-'.у ' 4 '

где р. - плотность твердой фазы, 5т - толщина ПС при Т = Тт, безразмерный коэффициент 2 = 1 соответствует числу атомных слоев в ПС, величина Q/ имеет смысл удельной теплоты плавления ПС (в расчете на единицу площади).

Для проверки корреляции (1) и лежащей в ее основе модели ПС были проведены расчеты поверхностного натяжения от вблизи точки плавления для двух групп веществ - жидких металлов и сжиженных инертных газов. Эти вещества имеют преимущественно симметричную кристаллическую решетку - объемноцентрированную (ОЦК) и гранецентрированную (ГЦК). Это обстоятельство позволяет оценить толщину ПС 5 = 2ё, ё - параметр решетки. Необходимые для расчета данные были взяты из справочника [8].

Результаты расчетов представлены в координатах Q/ -от (рисунок). Точки соответствуют справочным значениям от; при отсутствии этих сведений натяжения от определялись линейной экстраполяцией температурной зависимости о(Т).

Для щелочных металлов (решетка ОЦК) линейная зависимость от = /(бт) выполняется вполне удовлетворительно (рисунок, а). По углу наклона расчетное значение числа атомных слоев в ПС 2 = 3.7, т.е. твердообраз-ный ПС включает возле точки плавления примерно 4 монослоя. Существенно, что коррелирующая прямая в соответствии с (1) выходит из начала координат. Отметим, что последовательность элементов на рисунке отвечает уменьшению их металлических свойств - от цезия к литию.

Для других металлов с ОЦК-решеткой (рисунок, б) линейная корреляция от = /йт) также имеет место, хотя отклонения от нее заметно больше, чем для щелочных металлов. Число атомных слоев в ПС 2 = 2.5.

Для металлов с решеткой ГЦК линейное соотношение от = / (бт) можно вполне считать удовлетворительным (рисунок, в); 2 = 2.

Для сжиженных инертных газов (решетка ГЦК) линейная корреляция от = / (2т) выполняется очень четко, элементы располагаются в той же последовательности, как и

Корреляции поверхностного натяжения от вблизи точки плавления с теплотой плавления Q ^ : а - щелочные металлы, б - металлы с решеткой ГЦК, в - металлы с решеткой ОЦК, г - сжиженные инертные газы

в соответствующей группе Периодической системы (рисунок, г). Число атомных слоев 2 = 0.89, т.е. весьма близко к монослою [9].

Таким образом, для перечисленных выше элементов корреляция (1) дает достаточно хорошее согласие с экспериментальными значениями поверхностных натяжений вблизи точки плавления.

Температурный коэффициент поверхностного натяжения

Для теоретического расчета температурных коэффициентов поверхностного натяжения может быть использовано термодинамическое уравнение -ёо/ёТ = п, П - удельная поверхностная энтропия. Как и другие термодинами-

ческие параметры ПС, величина п представляет алгебраический избыток; она представляет разность энтропий ПС и объемной жидкости, отнесенную к единице площади. Поверхностная энтропия определяется уравнением П = Ц// Т, где Ц/ - теплота образования ПС при температуре Т [1].

Какие-либо сведения о величине Ц/ (расчетные, экспериментальные) отсутствуют. Поэтому независимые расчеты поверхностной энтропии п без привлечения температурных коэффициентов ёо /ёТ до сих пор, по-видимому, не проводились.

Рассмотренная выше модель ПС позволяет рассчитать по уравнению (1) теплоту образования поверхностного слоя при температуре плавления Тт. Далее рассчитывается поверхностная энтропия п (в точке плавления)

Т а б л и ц а 2

Температурные коэффициенты -<&/ЛТ поверхностного натяжения сжиженных газов, мДж/(м2К)

п = О /Т ) р 5 .

ут т7 г^ т

(2)

Металл Расчет Эксперимент

Литий 0.135 0.14

Натрий 0.100 0.10

Калий 0.064 0.066

Рубидий 0.057 0.058

Цезий 0.049 0.048

Медь 0.290 0.24

Серебро 0.230 0.18

Золото 0.220 0.17

Ртуть 0.200 0.20

Железо 0.230 0.23

Никель 0.310 0.34

Алюминий 0.270 0.15

Олово 0.220 0.09

Германий 0.630 0.095

Газ Расчет Эксперимент

Азот 0.10 0.23

Аргон 0.11 0.24

Водород 0.09 0.16

Гелий 0.19 0.10

Криптон 0.14 0.23

Ксенон 0.12 0.17

Кислород 0.12 0.24

Неон 0.20 0.35

Фтор 0.13 0.27

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

Уравнение (2) было использовано для расчета поверхностной энтропии пт жидких металлов и сжиженных газов с последующим сопоставлением с экспериментальны-

Т а б л и ц а 1

Температурные коэффициенты -Ло/ЛТ поверхностного натяжения жидких металлов, мДж/(м2К)

ми значениями температурного коэффициента поверхностного натяжения -ёо/ёТ (табл. 1, 2). При обсуждении представленных данных нужно учесть следующие обстоятельства.

Уравнение (2) позволяет рассчитать поверхностную энтропию и соответственно температурный коэффициент поверхностного натяжения только в точке плавления. Экспериментальные значения -(ёо/ёТ) рассчитаны по температурной зависимости о(Т) для интервала, где эта зависимость линейна. Для большинства веществ минимальные температуры, при которых измерялись поверхностные натяжения, по крайней мере на несколько градусов превышают температуру плавления. Поэтому сопоставление расчетных и экспериментальных значений температурных коэффициентов требует допущения: пт = пТ, т.е. поверхностная энтропия, рассчитанная для температуры плавления Т , остается постоянной и при более высоких температурах Т > Тт.

В проводившихся расчетах принималось, что толщина ПС 5т = ё, т.е. затвердевший ПС является монослоем. Необходимые для расчета данные взяты из [8].

Данные, приведенные в табл. 1, 2, позволяют сделать следующие заключения. Для щелочных металлов расчет дает полное согласие с экспериментом. Расхождение не превышает 4%. Достаточно близкое соответствие имеет место и для большинства металлов с решетками ГЦК и ОЦК (исключение - алюминий). Сильное превышение расчетных величин характерно для элементов со сложными решетками (олово, германий).

Напомним, что расчеты по уравнению (2) проводились при допущении 5т = ё. Поэтому, строго говоря, можно сделать вывод только об удовлетворительном выполнении линейной корреляции температурных коэффициентов поверхностного натяжения и энтропии плавления, а не об абсолютной точности данной схемы. Действительно, даже в самой благоприятной группе (щелочные

Т а б л и ц а 3

Автоадсорбция ГА (кг/м2) на поверхности жидкий металл - насыщенный пар вблизи температуры плавления Тт. Столбцы А и В рассчитаны по уравнениям (4) и (5)

Металл Vs■ 10 3, м/с -ДУ/У-102 Га 109, кг/м2 Га/Гв

А В

№ 2.4-2.65 2.5 31 10.4 2.98

K 1.82-1.88 2.5 32 11.5 2.79

Rb 1.26 2.5 55 21.8 2.52

Cs 0.967 2.6 76.6 30.1 2.54

Од 3.46 4.1 114 105 1.09

Ag 2.71-2.77 4.4 120 117.6 1.02

А1 4.67-4.73 6.4 41 55 0.74

Оа 2,74-2,87 3.4 90 69.6 1.8

1п 2.21-2.31 2.5 110 69 1.58

Т1 1.625 3.2 175 146 1.2

гп 2.7-2.85 6.5 100 144.5 0.69

са 2.2-2.25 2.5 112 75.9 1.47

8п 1.45 2.7 105 64.3 1.63

РЬ 2.27-2.48 3.44 135 75.5 1.79

Ы 1.79-1.88 3.22 137 130.3 1.05

^ 1.65 3.5 221 174.2 1.27

металлы) толщина ПС 5т = ё отличается от приведенных выше оценок по уравнению (1), где 5т = 3.7ё. Скорее всего это расхождение можно объяснить различным характером температурной зависимости о(Т) вблизи точки плавления и вдалеке от нее. Тем не менее этот вопрос, несомненно, требует более глубокого изучения.

Для всех сжиженных газов (кроме гелия) расчеты занижают температурные коэффициенты примерно вдвое (табл. 2). Формально это означает, что толщина ПС 5т соответствует не одному монослою, а двум (5т = 2ё). Это снова расходится с расчетными 5т по поверхностному натяжению о (5 = 0.89ё). Таким образом, в отличие от

тт

жидких металлов, для совпадения расчета с экспериментом нужно увеличивать, а не уменьшать толщину ПС. Этот результат также показывает, что достаточно достоверно можно установить только линейность корреляции -

(ёо/ёТ) = /(0 / Т).

тт

Корреляция с автоадсорбцией

Равновесное состояние ПС однокомпонентных систем наряду с поверхностным натяжением о и поверхностной энтропией п характеризует автоадсорбция ГА. По определению, она представляет уплотнение массы вещества в ПС, т.е. избыток массы ПС по сравнению с объемом кон -денсированной фазы. Приближенно автоадсорбция на поверхности однокомпонентной жидкости

ГА = (р./- р)5, (3)

где р./ - плотность ПС, р - плотность жидкости, 5 - толщина ПС [1].

По определению автоадсорбция ГА представляет прирост массы ПС, а поверхностное натяжение о - прирост его энергии по сравнению с соответствующими объемными величинами (в расчете на единицу площади ПС).

Тогда для описания связи между величинами о и ГА можно использовать принцип Эйнштейна об эквивалент-

_ т^ 2

ности массы и энергии: о = Гаус , ус - предельная скорость переноса энергии.

Предположим, что vc = V., где V. - скорость распространения звука в жидкости (высказанное предположение можно аргументировать тем, что предельная скорость распространения капиллярных волн на поверхности жидкости близка к скорости звука, а максимальная скорость распространения упругих волн в конденсированной фазе - это скорость звука).

Таким образом, применение принципа эквивалентности к поверхностным явлениям дает корреляцию

^ = ГА V/

(4)

Уравнение (4) проверяли следующим образом. Вначале на основе справочных данных [8] рассчитывали ав-

тоадсорбцию ГА различных жидких металлов вблизи температуры плавления Тт (столбец А табл. 3). Эти значения сопоставляли с расчетами автоадсорбции по уравнению (3). В принятой модели плотность ПС р^ = р^, где р^ - плотность твердой фазы. Тогда уравнение (3) преобразуется к виду

ГА = р S (ДУ/У), (5)

где ДУ/У - изменение объема при плавлении. Значения ГА, вычисленные по уравнению (5), представлены в столбце В (табл. 3). Отношение автоадсорбций, рассчитанных по уравнениям (4) и (5) соответствует относительной толщине ПС Z = 5/d.

Данные табл. 3 показывают, что значения Z лежат примерно в тех же пределах, что и при расчетах по теплоте плавления. Это соответствие можно считать подтверждением уравнения (4).

СПИСОК ЛИТЕРАТУРЫ

1. Семенченко В.К. Поверхностные свойства металлов и

сплавов. М., 1957.

2. Адамсон А. Физическая химия поверхностей. М., 1979.

3. Абрамзон А.А. Поверхностно-активные вещества. Л., 1981.

4. Wu X.Z., Ocko B.M., Sirota E.B. et al. // Science. 1993. 261. P. 1018.

5. Regan M.J., Kawamoto E.H., Lee L. et al. // Phys. Rev. Let.

1995. 75. P. 2498.

6. Сумм Б.Д. // Вестн. Моск. ун-та. Сер. 2. Химия. 1993. 34. С. 59.

7. Сумм Б.Д. Рос. хим. журнал им. Менделеева. 1994. 38.

С. 19.

8. Физические величины. Справочник. Под ред. И.С. Григорьева, Е.З. Мейлихова. М., 1991.

9. Сумм Б.Д., Ильичев ЕЮ. // ЖФХ. 1996. 70. С. 291.

Поступила в редакцию 20.04.99