Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,441

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

Мураева О.А. 1 Панаетова Т.Д. 1
1 Харьковский национальный университет городского хозяйства имени А.Н. Бекетова
В этой работе в рамках одной модели была сделана попытка дать полуколичественное описание физических явлений, связанных с сольватацией электролитов в разбавленных водных и неводных растворах. Предложено полуэмпирическое уравнение для описания чисел сольватации электролитов Z. Построение уравнения потребовало введения обобщенного радиуса Онзагера для случаев конкуренции трех вкладов в энергию молекулы растворителя: энергии взаимодействия молекул растворителя друг с другом, энергии взаимодействия растворителя с ионами и кинетической энергии теплового движения. Проведено сравнение выводов, полученных при анализе полуэмпирического уравнения и экспериментальных данных о числах сольватации электролитов Z, которые были рассчитаны из данных о динамической вязкости растворов электролитов. Полученное уравнение адекватно описывает положительную и отрицательную сольватации, а также позволяет предсказать зависимость чисел сольватации от ряда факторов (заряд и размер ионов, диэлектрическая постоянная растворителя). Подтверждено, что граница между положительной и отрицательной гидратацией находится между катионами натрия и калия и зависит от природы аниона. Вычислены критические значения радиусов катионов, при которых положительная сольватация невозможна. Полученные результаты полностью согласуются с положениями теории сольватации О.Я. Самойлова.
растворы электролитов
положительная и отрицательная сольватация
числа сольватации (гидратации)
заряд иона
радиус иона
диэлектрическая постоянная растворителя
1. Самойлов О.Я. Структура водных растворов электролитов и гидратация ионов / О.Я. Самойлов. – М.: Наука, 1976. – 256 с.
2. Вода: структура, состояние, сольватация. Достижения последних лет / под ред. А.М. Кутепова. – М.: Наука, 2003. – 402 с.
3. Смирнов П.Р. Структурные параметры ближнего окружения ионов в водных растворах неорганических электролитов / П.Р. Смирнов, В.Н. Тростин. – Иваново: ОАО «Издательство «Иваново», 2011. – 400 с.
4. Бакеев М.И. Теория гидратации и физико-химические свойства растворов электролитов / М.И. Бакеев. – Караганда: Изд-во КарГУ им. Е. А. Букетова, 2005. – 180 с.
5. Афанасьев В.Н. Верификация сольватации электролитов от разбавленных до концентрированных в водных растворах / В.Н. Афанасьев, А.Н. Устинов // Вестн. Моск. ун-та. – Сер. 2. Химия. – 2011. – Т. 52(5). – С. 323–340.
6. Абросимов В.К. Развитие физической химии жидких растворов в трудах отечественных ученых (Научные школы К.П. Мищенко, Н.А. Измайлова, О.Я. Самойлова и Г.А. Крестова) / В.К. Абросимов // Изв. вузов. Химия и хим. технология. – 2008. – Т. 51, № 9. – С. 5–10.
7. Кошкин В.М. Относительная вязкость и числа сольватации в растворах / В.М. Кошкин, О.А. Мураева // Теорет. и эксперим. химия. – 1985. – Т. 21, № 5.– С. 627–631.
8. Кошкин В.М. Примеси в полупроводниках со стехиометрическими вакансиями / В.М. Кошкин, Л.В. Атрощенко, Ю.А. Фрейман // Докл. АН СССР. 1968.– Т. 183, № 1. – С. 83–86.
9. Родникова М.Н. Отрицательная гидратация ионов / М.Н. Родникова // Электрохимия. – 2008. – Т. 39, № 2. – С. 214–219.
10. Бутырская Е.В. Квантовохимический расчет гидратации солей щелочных металлов / Е.В. Бутырская, В.А. Шапошник, А.М. Бутырский, А.Г. Рожкова // Ж. структ. химии. – 2006. – № 47. – С. 91.

Сольватация (для водных растворов гидратация) – одно из определяющих явлений в растворах электролитов, ответственное практически за все их свойства. На сегодняшний день различными методами анализа открыты многие закономерности, связывающие изменение какого-либо физического свойства с определенным параметром сольватации [1–3]. Большое количество работ посвящены разработке методов определения или расчета важного параметра сольватации – числа сольватации [4, 5]. Известны несколько теоретических моделей строения растворов, обзор и анализ которых приведен в [6]. Однако следует признать, что все эти модели основаны на различных приближениях, содержащих ряд эмпирических параметров, физический смысл которых часто носит формальный характер. Они не способны охватить и описать все известные проявления сольватации, зачастую применимы либо в относительно узком интервале концентраций и температур, либо для описания отдельных систем. Поэтому поиск сколько-нибудь общих закономерностей, даже полуколичественных, связывающих явление сольватации с микроскопическими параметрами растворенного вещества или растворителя, по-прежнему остается актуальной задачей.

Цель данной работы – попытаться в рамках молекулярно-динамической теории растворов О.Я. Самойлова получить полуэмпирическое уравнение для описания чисел сольватации электролитов Z в разбавленных водных и неводных растворах, провести анализ этого уравнения и сопоставить полученные выводы с экспериментальными данными о числах сольватации электролитов Z, которые были рассчитаны по данным о динамической вязкости растворов методом, предложенным авторами в работе [7].

Теория. Основные формулы

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

Полная энергия молекулы растворителя в ионном растворе W может быть представлена как

mur01.wmf (1)

где Uсв. – энергия межмолекулярного взаимодействия молекул растворителя; Екин. – кинетическая энергия теплового движения молекул растворителя, равная Екин. mur02.wmf kT; k – константа Больцмана; Т – абсолютная температура; mur03.wmf – энергия ион-дипольного взаимодействия; q – заряд иона; ε – диэлектрическая постоянная растворителя; r – расстояние между ионом и молекулой растворителя; μ – дипольный момент молекулы растворителя.

В соответствии с тем, какой из трех энергетических вкладов преобладает в уравнении (1) при данном значении r, строение раствора определяется именно этим фактором. Если преобладает вклад Uсв., то расположение молекул растворителя подобно строению кристалла (у воды это льдоподобная структура). Если преобладающей оказывается энергия ион-дипольного взаимодействия, то молекулы растворителя выстраиваются так, чтобы направление их дипольных моментов совпадало с вектором напряженности электрического поля иона. Будем называть такое упорядоченное построение молекул растворителя вблизи иона сольватной структурой. Наконец, если энергия теплового движения существенно превышает все другие члены в уравнении (1), то любое упорядочение отсутствует, молекулы растворителя распределены в растворе хаотично.

Определим радиус сольватации r0 как такое максимальное расстояние, при котором r < r0, и энергия ион-дипольного взаимодействия превышает абсолютное значение суммы двух других вкладов в общую энергию молекулы в уравнении (1):

mur04.wmf, (2)

тогда

mur05.wmf (3)

Очевидно, что при Uсв. = 0 радиус сольватации есть радиус Онзагера

mur06.wmf (4)

Выражение (3), таким образом, представляет собою обобщенный радиус Онзагера. Разумеется, уравнение (3) неприменимо вблизи температуры mur07.wmf. При mur08.wmf предлагаемая в данной работе качественная модель должна быть существенно уточнена учетом флуктуации параметров порядка.

Как в уравнении (3), так и в уравнении (4), вообще говоря, нужно ввести полное значение дипольного момента

mur09.wmf, (5)

где α – поляризуемость молекулы растворителя; μ0 – жесткий дипольный момент молекулы растворителя; Е – напряженность электрического поля, создаваемого ионом в точке r0, равная

mur10.wmf (6)

Кроме того, следует учесть также то обстоятельство, что значение ε выражается через α и μ0 по уравнению Клаузиуса – Мосотти – Дебая:

mur11.wmf (7)

где d – плотность растворителя; М – молярная масса растворителя.

С учетом уравнений (5) и (6), уравнение (3) приобретает следующий вид:

mur12.wmf (8)

Тогда для r0 получаем уравнение

mur13.wmf

в окончательном виде

mur14.wmf, (9)

если Uсв. << kT, то

mur15.wmf (10)

Если поляризуемость α мала, то уравнения (9) и (10) сводятся, естественно, к уравнениям (3) и (4). Если собственный дипольный момент молекулы растворителя μ0 мал и в выражении (8) им можно пренебречь, а поляризуемость α значительна, то уравнение (9) будет иметь следующий вид:

mur16.wmf (11)

и радиус Онзагера (при Uсв. = 0)

mur17.wmf (12)

Исходя из сделанных выше определений, перейдем к математическому выражению для чисел сольватации Z. Представим объем сольватированной оболочки иона как

mur18.wmf, (13)

а количество молекул растворителя, входящих в сольватную оболочку иона (число сольватации) Z с учетом уравнения (3) в виде

mur19.wmf (14)

где V0 – объем, приходящийся на одну молекулу растворителя.

В случае структурированной жидкости, например воды при температурах, близких к комнатной температуре, ионы, как известно, могут локализоваться либо в позиции «замещения» молекулы растворителя, либо в позиции «внедрения» в межмолекулярные пустоты [1, 2]. При любом характере локализации размеры соответствующего иона нужно соотнести с характерным эффективным радиусом ρ того пространства в структуре растворителя, которое занимает ион (размерами пустот в случае замещения или размерами молекул растворителя в случае внедрения). Для случая воды крупные ионы (размер иона больше ρ) разрушают структуру растворителя – межмолекулярные связи вблизи внедрившегося иона разрываются или, по крайней мере, ослабляются. Подобное происходит в кристалле вблизи крупного внедренного атома: связи в кристалле-матрице оказываются растянутыми, и в ряде моделей это хорошо описывается упругими напряжениями [8]. Если радиус иона меньше, чем ρ, то структура и межмолекулярные связи растворителя также нарушаются, молекулы «стягиваются» к маленькому иону, нарушая связи между собой. Только если радиус иона в точности совпадает с ρ, то напряжения и разрыва связей не происходит. Логично предположить, что размеры области, где межмолекулярные связи растворителя нарушены за счет объемных эффектов при внедрении ионов в структуру растворителя увеличиваются с ростом различий в объеме mur20.wmf и объеме внедренных ионов. Поскольку обычно корректно не удается разделить вклад в сольватацию от каждого из ионов электролита, и реально в эксперименте определяется брутто-значение Z для электролита KmAn в целом, то объем растворителя Vd, нарушенный введением m катионов с радиусом rк и n анионов с радиусом ra равен

mur21.wmf (15)

где B – константа.

Ясно, что если размеры введенных ионов совпадают с ρ, то Vd = 0. Введение модуля разности в уравнение (15) учитывает то обстоятельство, что как очень малые, так и очень крупные ионы разрушают структуру растворителя. Если область сольватного упорядочения VS простирается дальше от иона, чем область «геометрически» разрушенной структуры жидкости Vd (т.е. VS > Vd), то имеет место положительная сольватация. При обратном неравенстве VS < Vd, когда структура растворителя разрушена на расстоянии большем чем то, на котором проявляется сольватное упорядочение, имеет место отрицательная сольватация по Самойлову.

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

mur22.wmf (16)

С учетом уравнения (3)

mur23.wmf (17)

Очевидно, что при VS > Vd число гидратации Z будет больше нуля, а в случае VS < Vd – меньше нуля.

Особо следует остановиться на случае очень больших значений r0 в структурированных жидкостях. В проведенных выше оценках не учитывался тот факт, что характерные размеры D льдоподобных структурных областей в жидкой воде, например, ограничены. За пределами такой области жидкость неупорядочена. Если r0 < D, то граница гидратного упорядочения проходит в пределах структурного кластера воды. Имеет место конкуренция двух типов упорядочения, что и отражается в обобщенном радиусе Онзагера. Если же r0 > D, то граница гидратного упорядочения проходит в области, где льдоподобная структура уже отсутствует. Поскольку с ростом температуры значение D(T) уменьшается, размер гидратного образования при больших r0 и Z определяется конкуренцией гидратного упорядочения уже не с межмолекулярными связями, а с энергией тепловых колебаний. При больших Z это соответствует случаю Uсв. = 0.

Сопоставление значений mur27.wmf с величиной mur28.wmf (водные растворы)

Электролит (1): электролит (2)

mur27.wmf

mur30.wmf

Электролит (1): электролит (2)

mur27.wmf

mur31.wmf

mur32.wmf:mur33.wmf

2,6

2,8

mur34.wmf:mur35.wmf

1,3

1

mur36.wmf:mur37.wmf

2,7

2,8

mur38.wmf:mur39.wmf

1,7

1

mur40.wmf:mur41.wmf

2,6

2,8

mur42.wmf:mur43.wmf

0,97

1

mur44.wmf:mur45.wmf

2,5

2,8

mur46.wmf:mur47.wmf

1,1

1

mur48.wmf:mur49.wmf

2,8

2,8

mur50.wmf:mur51.wmf

1,4

1

mur52.wmf:LiBr

2,9

2,8

mur53.wmf:mur54.wmf

1,3

1

mur55.wmf:mur56.wmf

3,0

2,8

mur57.wmf:mur58.wmf

1,7

1

mur59.wmfmur60.wmf

2,9

2,8

mur61.wmf:mur62.wmf

1,1

1

mur63.wmfmur64.wmf

2,4

2,8

mur65.wmf:mur66.wmf

1,4

1

mur67.wmf:LiCl

4,0

5,2

mur68.wmf:mur69.wmf

1,4

1

mur70.wmf:mur71.wmf

1,5

1,8

mur72.wmf:mur73.wmf

1,7

1

mur74.wmf:mur75.wmf

1,4

1,8

mur76.wmf:mur77.wmf

0,99

1

mur78.wmf:mur79.wmf

1,4

1,8

mur80.wmf:mur81.wmf

1

1

mur82.wmf:mur83.wmf

1,5

1,8

mur84.wmf:mur85.wmf

1,4

1

mur86.wmf:mur87.wmf

2,0

1,8

mur88.wmf:mur89.wmf

1,2

1

mur90.wmfmur91.wmf

mur92.wmf:mur93.wmf

2,0

1,8

mur94.wmf:mur95.wmf

1,2

1

mur96.wmf:mur97.wmf

2,3

1,8

mur98.wmf:mur99.wmf

1,3

1

mur100.wmf:mur101.wmf

0,99

1

mur102.wmf:mur103.wmf

1,2

1

mur104.wmf:mur105.wmf

0,96

1

mur106.wmf:mur107.wmf

1,2

1

mur108.wmf:mur109.wmf

I,1

1

mur110.wmf:mur111.wmf

1,2

1

 

muraev1.wmf

Рис. 1. Зависимость чисел сольватации Z от диэлектрической постоянной растворителя ε при 273 К: а – LiCl в изопентаноле (1), метаноле (2), N,N-диметилформамиде (3), диметилсульфоксиде (4), воде (5); б – LiNO3 в изопентаноле (1), бутаноле (2), метаноле (3), воде (4)

Результаты исследования и их обсуждение

Проведем анализ уравнения (17) и сопоставим полученные выводы с экспериментальными данными о числах сольватации электролитов Z, которые были рассчитаны по данным о динамической вязкости растворов методом, предложенным авторами в работе [7]. Определяемые этим методом числа сольватации представляют собой среднестатистическую величину, суммарно отражающую сольватацию катиона и аниона.

Зависимость чисел сольватации Z от заряда ионов. Из уравнения (17) следует, что Z должно расти при увеличении заряда ионов. При больших значениях Z, когда вторым членом в уравнении (17) можно пренебречь, появляется возможность количественной проверки. Действительно, в этом случае для двух солей

mur24.wmf

Эта зависимость неизменна, как в случае полярных растворителей, так и в случае неполярных, но поляризуемых. В таблице приведены сопоставления, полученных значений mur25.wmf с величиной mur26.wmf для водных растворов.

Соответствие следует признать хорошим. Для неводных растворов имеет место аналогичная закономерность.

Зависимость Z от диэлектрической постоянной растворителя. В соответствии с уравнением (17) с ростом диэлектрической постоянной растворителя ε величина Z должна уменьшаться, что действительно, имеет место (рис. 1).

Зависимость чисел сольватации от радиуса иона. Существование отрицательной гидратации и то, что граница между положительной и отрицательной гидратацией находится между ионами Na+ и K+, давно общепризнано. Впервые это было показано в работе [1], а затем подтверждено различными методами [2, 4, 5, 9, 10]. Было установлено критическое значение радиуса катиона, при котором отрицательная гидратация переходит в положительную, при 293 К оно равно mur112.wmf [1], при 298 К – mur113.wmf [1, 6, 9], для однозарядных анионов mur114.wmf [6].

muraev2.wmf

Рис. 2. Зависимость Z от радиуса ионов при 293 К: а – LiI(1), NaI(2), KI(3); б – LiBr(1), NaBr(2), KBr(3), RbBr(4); в – LiCl(1), NaCl(2), KCl(3), RbCl(4), CsCl(5); г – KF(1), KCl(2), KBr(3), KI(4); д – NaCl(1), NaBr(2), NaI(3); е – LiCl(1), LiBr(2), LiI(3)

Рассмотрим, насколько адекватно уравнение (17) позволяет оценить зависимость чисел сольватации от размера ионов и определить значение критических радиусов катиона mur115.wmf и аниона mur116.wmf, при которых положительная гидратация переходит в отрицательную. Из уравнения (17) следует, что значения Z могут принимать как положительные, так и отрицательные значения. Из условия Z = 0 можно определить mur117.wmf и mur118.wmf.

Для катиона

mur119.wmf, (18)

аналогично для аниона

mur120.wmf (19)

Из уравнений (18) и (19) следует, что величина mur121.wmf должна уменьшаться с ростом размера аниона mur122.wmf mur123.wmf, а mur124.wmf – с ростом размера катиона mur125.wmf

График зависимости величин Z различных электролитов от кристаллографических радиусов катионов и анионов (по Полингу) действительно демонстрирует переход от положительной к отрицательной гидратации при увеличении как радиуса катиона, так и радиуса аниона (рис. 2).

На рис. 2, а–в, представлены зависимости Z от радиуса катионов, которые подтверждают, что граница между положительной и отрицательной гидратациями действительно находится между катионами Na+ и K+, зависит от природы аниона (чем больше радиус аниона, тем меньше mur126.wmf) и приходится на критический радиус катиона mur127.wmf (mur128.wmf), равный для хлоридов 1,27, нитратов 1,21, бромидов 1,1, иодидов 1,07, что совпадает с данными работ [1, 6, 9].

Для нахождения критического радиуса аниона, к сожалению, недостаточно данных, есть только данные зависимости Z от радиуса анионов для хлоридов, бромидов, иодидов лития, натрия и калия (рис. 2, д–е). Видно, что переход от положительной гидратации к отрицательной имеет место только для галогенидов калия, лежит между ионами фтора и хлора и происходит при mur129.wmf, что согласуется с данными работ Г.А. Крестова [6]. Следует полагать, для галогенидов рубидия и цезия критическое значение радиуса аниона будет меньше.

Выводы

Предложено полуэмпирическое уравнение (17) для описания чисел сольватации электролитов в водных и неводных растворителях. Проведена апробация этого уравнения, показавшая, что уравнение адекватно описывает положительную и отрицательную сольватацию в разбавленных растворах, подтверждает наличие критического радиуса катиона и аниона, при которых происходит переход от положительной сольватации к отрицательной, позволяет прогнозировать зависимость чисел сольватации от заряда и размера ионов, диэлектрической постоянной растворителя. Полученные результаты полностью согласуются с положениями теории сольватации О.Я. Самойлова, демонстрируют не только валидность предложенного полуэмпирического уравнения, но и надежность вискозиметрического метода определения Z.


Библиографическая ссылка

Мураева О.А., Панаетова Т.Д. ПОЛУЭМПИРИЧЕСКОЕ УРАВНЕНИЕ ДЛЯ ОПИСАНИЯ ЧИСЕЛ СОЛЬВАТАЦИИ ЭЛЕКТРОЛИТОВ // Фундаментальные исследования. – 2018. – № 4. – С. 21-27;
URL: http://www.fundamental-research.ru/ru/article/view?id=42112 (дата обращения: 23.04.2021).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074