Теория электролитов Дебая и Гюккеля: вывод первого приближения теории.
Основные положения современной теории растворов электролитов были сформулированы в 1923 г. Дебаем и Гюккелем. Для статистической теории электролитов исходным является следующее положение: ионы распределены в объеме раствора не хаотически, а в соответствии с законом кулоновского взаимодействия. Вокруг каждого отдельного иона существует ионная атмосфера (ионное облако) - сфера, состоящая из ионов противоположного знака. Ионы, входящие в состав сферы, непрерывно обмениваются местами с другими ионами. Все ионы раствора равноценны, каждый из них окружен ионной атмосферой, и в то же время каждый центральный ион входит в состав ионной атмосферы какого-либо другого иона. Существование ионных атмосфер и есть тот характерный признак, который, по Дебаю и Гюккелю, отличает реальные растворы электролитов от идеальных.
Энергия, связанная с ионной атмосферой, имеет электрическое происхождение, поэтому она должна быть функцией плотности электрического заряда и потенциала, создаваемого ионной атмосферой. Тепловое движение ионов в ионной атмосфере приводит к тому, что дискретные заряды этих ионов как бы размазываются. В результате ионную атмосферу, состоящую из отдельных ионов, в среднем за некоторый промежуток времени можно моделировать облаком размазанного заряда, плотность которого r уменьшается по мере удаления от центрального иона. Для связи между средней плотностью заряда r и отвечающим ему средним значением потенциала y можно использовать уравнение Пуассона
где Ñ2 – оператор Лапласа Ñ2 = + + .
Поскольку ионная атмосфера обладает шаровой симметрией, то уравнение Пуассона, записанное в сферической системе координат, принимает вид
и связывает объемную плотность заряда r с потенциалом y на расстоянии r от центрального иона. В уравнение Пуассона входят две неизвестные величины: r и y. Для их определения необходимо иметь второе уравнение, связывающее обе эти переменные. Дебай и Гюккель получили нужное уравнение следующим образом.
Для раствора в целом и для любого заданного его объема, достаточно большого по сравнению с размерами иона, справедливо условие электронейтральности
å niezi = 0 ,
где ni– число ионовi-го сорта в единице объема (если выбран объем раствора, равный единице). Однако заряд некоторого элемента объема dV, находящегося вблизи какого-либо иона, будет отличаться от нуля вследствие существования ионной атмосферы. Если центральный ион, расположенный в начале координат, заряжен положительно, то элемент объема dV будет обладать избыточным отрицательным зарядом. Предполагая, что к распределению ионов в растворе применим принцип Больцмана и что силы, действующие между ионами, по своей природе электростатические, число отрицательных и число положительных ионов в элементарном объеме dV можно выразить как
они ограничились только двумя его первыми членами, получив для r следующее уравнение:
r = e å zi ni – å zi2 ni y = – å zi2 ni y
(так как первый член правой части уравнения по условию электронейтральности равен нулю). Если ввести обозначение
то после подстановки выражения для r в уравнение Пуассона получаем дифференциальное уравнение вида
Общее решение этого уравнения второго порядка имеет вид
Константы интегрирования А1 и А2 находят из пограничных условий. Как следует из основных законов электростатики, при r ®¥ y® 0; это условие выполняется в том случае, если А2 = 0, так как иначе с удалением от центрального иона y будет стремиться к бесконечности. Для определения А1 в первом приближении теории Дебая – Гюккеля предполагается, что ионы можно отождествить с материальными точками, обладающими определенным зарядом. В этом случае при r ® 0 y должен стремиться к потенциалу самого иона yi , то есть
что является вторым пограничным условием. Если е–cr разложить в ряд, то получим
y = А1 = А1 (1 – cr + (cr)2 – …) .
При r ® 0 всеми членами ряда можно пренебречь по сравнению с 1. Тогда
yr® 0 ® yi = А1/r , A1 = , y = exp (–cr) .
Величина y в последнем уравнении представляет собой среднее значение потенциала в точке r, создаваемое ионной атмосферой и центральным ионом. Для вычисления ион-ионного взаимодействия представляет интерес не общий потенциал y, а та его часть yа , которая создается ионной атмосферой в месте расположения центрального иона. yа называется потенциалом ионной атмосферы. Его находят по правилу суперпозиции (то есть наложения) потенциалов
yа = y – yi = (exp (–cr) – 1) .
Для вычисления энергии взаимодействия необходимо определить потенциал ионной атмосферы в точке нахождения центрального иона, то есть найти предел yа при r ® 0. Это можно сделать, используя тот же прием – разлагая показательную функцию в ряд и пренебрегая высшими членами разложения:
Величину y¢ можно рассматривать как потенциал, создаваемый в точке нахождения центрального иона другим ионом с противоположным знаком, находящимся от центрального иона на расстоянии 1/c. Величина 1/c называется характеристической длиной. Так как потенциал создается не единичным ионом, а всей ионной атмосферой, то 1/c можно отождествить с радиусом ионной атмосферы. Величину c можно рассчитать по вышеприведенному уравнению.
В первом приближении теории Дебая и Гюккеля предполагается, что энергия межионного взаимодействия имеет электростатическое происхождение. Энергия взаимодействия ионов может быть определена поэтому как энергия заряжения центрального иона в электрическом поле ионной атмосферы. Известно, что энергия заряжения частицы от 0 до заряда q в поле y равна
(вместо q подставили заряд центрального иона, а вместо y – величину y¢).
Найдем выражение для коэффициента активности электролита. Коэффициенты активности первоначально были введены как эмпирические величины; необходимость их введения в уравнения для идеальных систем обусловлена различием между реальными и идеальными растворами. Коэффициенты активности не зависят от метода их определения и при заданных условиях и составе будут иметь одно и то же значение для любых типов равновесия. Таблицы, в которых приведены эмпирические коэффициенты активности, представляют большую ценность. В то же время в пределах теории Льюиса коэффициенты активности оказываются лишь формальными поправочными множителями, не связанными непосредственно с природой растворов и не поддающимися теоретическим расчетам.
Реальные растворы отличаются от идеальных энергией взаимодействия образующими их частиц. Сопоставление выражений для химических потенциалов частиц в реальных растворах
mi = mio + RT ln Ni + RT ln gi
и в идеальных растворах = mio + RT ln Ni
показывает, что RT ln gi = mi – = Dmi .
Таким образом, коэффициент активности отражает изменение энергетического состояния ионов, обусловленное межионным (межчастичным) взаимодействием. Расчет взаимодействия, а следовательно, и коэффициента активности возможен, если известны структура раствора и природа сил взаимодействия между образующими его частицами. Ни то, ни другое не известно с достаточной степенью точности, поэтому для решения проблемы приходится прибегать к некоторым предположениям и создавать модели растворов электролитов.
По модели Дебая и Гюккеля, энергия взаимодействия иона с ионной атмосферой, равная энергии заряжения, определяется выражением
Заряд элемента объема будет равен
rdV = ez– n– exp dV – ez+ n+ exp dV ,
а плотность заряда
r = ez– n– exp – ez+ n+ exp = e å ziniexp
(учитывая, что в растворе находятся различные ионы, и приписывая zi знак, отвечающий заряду иона). Подстановка r из этого выражения в уравнение Пуассона дает уравнение, которое в общем виде нельзя проинтегрировать. Чтобы упростить его, Дебай и Гюккель предположили, что eziy<< kT , то есть что электростатическое взаимодействие мало по сравнению с термической энергией. Разложив показательную функцию в ряд
Для коэффициента активности
Подстановка значения c приводит к
где I – ионная сила раствора, определяемая соотношением I = å сizi2 .
[ni] = ион /м 3 ; [ni/ NA] = моль / м 3 ; [ni/(1000×NA) = ci] = моль / л ).
После перехода к десятичному логарифму и подстановки численных значений постоянных величин
lggi = – 1,8246×106 (eT)– 3/2 zi2 .
Для данного растворителя и определенной Т можно написать
Для среднего ионного коэффициента активности, который в отличие от коэффициента активности отдельного вида ионов определяется экспериментально, справедливо уравнение
Для водных растворов при 25оС, допуская равенство диэлектрических проницаемостей раствора и растворителя (78,54), h = 0,51. Для 1-1 зарядных электролитов