Калин Физическое материаловедение Том 2 Основы материаловедения 2007
.pdf
2) связанной с появляющейся именно в растворах возможностью обмена местами (позициями) двух разнородных атомов; это так называемая позиционная составляющая Sp энтропии.
С другой стороны, энтропию раствора можно разбить на две логические части:
1) абсолютные энтропии чистых веществ Σxi Sio , имеющие только тепловую природу;
2) интегральную энтропию смешения S M.
Для нашего идеализированного раствора взаимное энергетическое взаимодействие одинаково для частиц всех компонентов. Поэтому тепловая составляющая для раствора имеет ту же величину, что и для чистых веществ по отдельности:
Sheat = Σxi Sio .
Следовательно, интегральная молярная энтропия смешения S M идеального раствора в другом определении есть не что иное, как
«позиционная энтропия»: |
|
Sp = S М . |
(4.108) |
Вычислим позиционную энтропию для бинарного раствора, исходя из формулы Больцмана Sp = k·lnWp, где Wp – вероятность реализации данного энергетического макросостояния через различные способы пространственного размещения (позиции) атомов.
Поскольку в нашем случае простейшей физической модели идеального раствора выполняется энергетическое условие H M = 0, то пары А-В энергетически ничем не предпочтительнее парам А-А или В-В, и наоборот. Тогда понятно, что:
•атомы видов А и В распределены в различных узлах решетки случайным образом, равновероятно в любой позиции;
•возможен обмен любых атомов местами без изменения внутренней энергии и энтальпии, характеризующих макросостояние.
Вэтом случае величина Wp равна числу перестановок из общего числа атомов (NА + NВ), деленному на произведение числа перестановок из чисел атомов каждого вида:
Wp = (NA + NB )! .
NA!NB!
111
Используя формулу Стирлинга lnN! = N·lnN – N, получаем для позиционной энтропии
Sp=k·lnWp = k·[(NА+NВ)ln(NА+NВ) – NА·lnNА – NВ·lnNВ]
В 1 моле раствора полное число атомов NА+NВ равно числу Аво-
гадро NA. Тогда имеем |
NА = xА·NA; NВ = xВ·NA; k·NA = R.. Отсюда |
Sp = R·(–xА·lnxА – xВ·lnxВ) |
и из (4.108) S М = R·(–xА·lnxА – xВ·lnxВ). |
С помощью основных термодинамических соотношений, для парциальных молярных энтропий смешения, для каждого из компонентов идеального раствора имеем:
SiM = −Rxi ln xi .
Соответственно, парциальные молярные свободные энергии смешения компонентов идеального раствора будут равны:
GiM = RT ln xi ,
а парциальные молярные свободные энергии будут равны:
|
G |
(T ) = |
Go(T ) + RT ln x . |
(4.109) |
|
|
i |
i |
i |
|
|
С разбивкой энергии Гиббса на энтальпийную и энтропийную составляющие, имеем:
Gi (T ) = Hio (T ) −T [Sio(T ) − R ln(xi )].
Заметим, что энергия смешения при образовании идеального раствора, определяемая как
GM = −T S M = RT ∑xi ln xi ,
есть величина всегда отрицательная (ибо xi < 1 всегда).
Это означает, что образование идеального раствора происходит самопроизвольно. При этом на концентрации компонентов идеального раствора не накладывается никаких дополнительных ограничений – они могут принимать любые значения от нуля до единицы. Говорят, чтов системеобразуетсянепрерывный ряд растворовили что всистемеимеютсярастворыснеограниченнойрастворимостью.
Реальные растворы замещения с неограниченной растворимостью ближе всего подходят под определение идеальных. На практике модель идеальных растворов в подавляющем большинстве случаев (при обычных температурах и давлениях, когда расстояния между молекулами велико) с хорошей степенью точности применима к смесям газов. Смесь обычных, не очень тяжелых газов сле-
112
дует начинать рассматривать как реальную, не описываемую моделью идеальных растворов. Это происходит или при низких температурах (Т < 100–50 К), или при повышенных давлениях (р > 10– 100 атм). Смесь жидкостей тоже часто с достаточной степенью точности можно рассматривать как идеальный раствор.
Модель идеальных растворов часто используется при теоретическом исследовании не изученных систем просто из-за отсутствия каких-либо сведений о механизмах взаимодействия компонентов реального раствора. Поэтому и делается простейшее предположение.
4.5.4. Идеальные бинарные системы
Рассмотрим случай бинарной смеси, когда неограниченная растворимость наблюдается как в первой, так и во второй фазах. Изучать будем простейшую математическую модель, когда приближение идеальных растворов применимо к обеим фазам. Реальные растворы будут, в принципе, обладать похожими свойствами, но несколько искаженными энергетическим взаимодействием разнородных частиц между собой.
В бинарной смеси, при равновесии двух фаз – двух растворов – жидкого (Liquid) и твердого (Solid) – имеем, в общем случае шесть неизвестных:
четыре концентрации: |
xAL , xBL и xAS , xBS ; |
одну температуру |
Т ; |
одно давление |
р. |
Поскольку давление мало влияет на равновесие в конденсированных средах, то в дальнейшем исключим его из числа неизвестных, которых останется всего пять. Не нарушая общности рассуждений, будем считать первым и более легкоплавким компонентом вещество А, вторым – более тугоплавкое вещество В.
Закон сохранения материи, частным случаем записи которого является условие нормировки концентраций, дает два уравнения
связи равновесных концентраций: |
|
|
xAL + xBL = 1; |
xAS + xBS =1. |
(4.110) |
Одно из явных условий установления термодинамического равновесия, а именно равенство химических потенциалов компонен-
113
тов в разных фазах (см. раздел 3), дает еще два уравнения, записываемые через парциальные молярные потенциалы Гиббса при произвольной (неизвестной) температуре:
|
|
AS = |
|
AL ; |
|
BS = |
|
BL . |
G |
G |
G |
G |
Тогда, в случае идеальных растворов, имеем еще два уравнения связи, теперь уже равновесных концентраций и температур:
GAoS (T ) + RT ln(xAS ) GBoS (T ) + RT ln(xBS )
= |
GAoL (T ) + RT ln(xAL ); |
(4.111) |
|
= |
GBoL (T ) + RT ln(xBL ). |
||
|
Перегруппировав члены в уравнениях (4.111), используя соотношения (4.110), окончательно получим два уравнения относитель-
но трех неизвестных величин: xBS , xBL , Т:
RT ln(1− xBS ) − RT ln(1− xBL ) |
= |
GAoL (T ) − |
GAoS (T ); |
(4.112) |
|
RT ln(xBS )−RT ln(xBL )= |
GBoL (T ) − GBoS (T ). |
||||
|
|||||
Неполные математические системы уравнений имеют бесконечное множество решений. В данном случае, при допущении о независимости потенциалов от давления, геометрическим местом точек решения системы (4.112) будут две линии:
xBS (T ) – описывающая температурную зависимость концентрации в твердой фазе и называемая линией солидуса;
xBL (T ) – описывающая температурную зависимость концентра-
ции в жидкой фазе и называемая линией ликвидуса. Несоответствие числа уравнений и неизвестных – обычный слу-
чай в термодинамике. Выход из такого положения в том, что ка- ким-либо неизвестным придаются некоторые произвольные, но имеющие физический смысл, численные значения, а затем решают получающиеся полные системы уравнений. При этом следует стремиться выбирать задаваемые переменные так, чтобы образовывалась максимально простая система математических уравнений.
Заметим теперь, что если задаваться составом жидкой или твердой фазы и рассматривать систему (4.112) относительно температуры как неизвестной величины, то система уравнений будет трансцендентна и решается только численно. Обратная задача –
нахождение неизвестных концентраций xBS и xBL в идеальном рас-
114
творе при заданной температуре – имеет аналитическое решение. Поэтому ниже остановимся именно на этом варианте.
В уравнениях (4.112) разницу потенциалов Гиббса G°L – G°S для двух фаз можно назвать свободной энергией плавления (melting) компонентов G°melt(Т) при произвольной температуре.
Обращаем внимание на тот факт, что входящая в (4.112) температура в общем случае отличается от температур плавления TAmelt , TBmelt чистых компонентов. Поэтому G°melt(Т) здесь не равно нулю,
как это было в однокомпонентных системах. Более того, в двухкомпонентных (и более сложных) системах нет однозначной температуры плавления: растворы начинают плавиться при одной температуре, а заканчивается процесс плавления при другой. Для растворов разных составов температуры начала и окончания процесса плавления будут разными. Находить эти температуры теоретически как раз и следует из решения системы уравнений (4.112).
Проведем несложное математическое преобразование в (4.112):
|
|
|
omelt |
|
S |
L |
|
GB |
(T ) |
xB (T ) = xB |
(T ) exp |
RT |
; |
|
|
|
|
|
|
S |
(T ) |
|
omelt |
|
|
1− xB |
|
GA |
(T ) |
||
|
|
= exp |
|
|
, |
L |
(T ) |
RT |
|
||
1− xB |
|
|
|
||
(4.113)
(4.114)
и, подставляя (4.113) в (4.114), получаем явную функцию для расчета зависимости концентрации от температуры:
|
|
|
|
|
omelt |
(T ) |
|
|
|
|
|
|
|
|
|
|
exp |
GA |
|
|
−1 |
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
RT |
|
|
|
|
|
|
|
|
|
xBL (T ) = |
|
|
|
|
|
|
|
|
|
|
. |
(4.115) |
||
|
|
|
|
|
|
omelt |
|
|
||||||
|
omelt |
|
|
|
|
|
||||||||
|
|
GA |
|
(T ) |
|
GB |
(T ) |
|
|
|||||
|
exp |
|
|
|
|
−exp |
|
|
|
|
|
|
||
|
RT |
|
|
|
|
RT |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||
Вычисления остальных концентраций xBS , xAL , xAS тривиальны и
идут по формулам (4.113) и (4.110). Перепишем уравнение (4.113) в виде:
xL |
|
Gomelt (T ) |
|
|
|
B |
=exp |
B |
= kr , |
(4.116) |
|
xBS |
RT |
||||
|
|
|
115
где kr – отношение равновесных концентраций в разных фазах, на-
зываемое коэффициентом распределения. Уравнение (4.116) носит название закона распределения Нернста.
Поскольку GBomelt ≠ 0 при любой температуре, отличной от TBmelt , из (4.116) очевидно, что kr ≠1 в общем случае. Это означает,
что равновесные концентрации одного и того же вещества в разных фазах принципиально не равны между собой. Даже в идеальных растворах. Данное свойство очень широко используется на практике для глубокой очистки веществ путем многократного повторения процессов испарения – конденсации (переход жидкость – пар) или плавления – кристаллизации (переход жидкость – твердое). Процессы очистки, основанные на использовании неравномерного распределения веществ при фазовых переходах, без химических реак-
ций, получили название физических методов очистки.
Конечно, из-за многообразия веществ в природе иногда получается такое сочетание термохимических свойств (Н, S, сp), что вели-
чина GBmelt случайно может оказаться близкой к нулю. Обычно это
происходит для веществ, у которых близки температуры фазовых переходов, и коэффициент распределения kr будет здесь равен примерно единице. Вариант kr = 1 может (случайно!) реализоваться на практике и из-за реальности раствора. Эти варианты очень неприятны для технологов, так как здесь невозможно добиться разделения компонентов; для них физические методы очистки будут малоэффективны.
4.6.Физическая химия неидеальных растворов
4.6.1.Термодинамика неидеальных растворов
Исходя из общих положений термодинамики для парциального
молярного потенциала Гиббса |
|
|
|
|
|
|||
Gi (T , x) i-го компонента такого |
||||||||
раствора, согласно (4.103) имеем: |
|
|
|
|||||
|
|
(T , x) = |
Go(T ) + |
|
M (T , x) , |
(4.117) |
||
G |
G |
|||||||
|
i |
|
i |
i |
|
|||
|
|
|
116 |
|
|
|
||
где парциальная молярная энергия смешения Гиббса GiM (T , x) –
дополнительная энергия, приобретенная чистым веществом из-за перехода в раствор. Это – пока неизвестная функция температуры и вектора концентраций х всех компонентов раствора.
На практике, кроме парциальной энергии смешения компонента GiM (T , x) , используют так же понятие избыточной парциальной свободной энергии Giизб(T , x) – той части энергии смешения
GiM в реальном растворе, на которую она превышает RT·ln(xi) – ее значение для идеальных растворов:
|
|
|
M |
(T , x) = RT ln(x ) + |
|
изб(T , x) |
|
||
G |
G |
|
|||||||
|
i |
|
|
i |
i |
|
|||
или |
|
|
|
|
|
|
|||
|
|
|
изб(T , x) = |
|
M (T , x) − RT ln(x ) . |
(4.118) |
|||
|
G |
G |
|||||||
|
|
|
i |
|
i |
|
i |
|
|
Представив в (4.118) избыточную свободную энергию через ее энтальпийную и энтропийную составляющие:
Giизб(T , x) = Hiизб(T , x) −T Siизб(T , x) =
=HiM (T , x) −T SiM (T , x) − RT ln(xi ) =
=HiM (T , x) −T [ SiM (T , x) + R ln(xi )],
после перегруппировки членов получаем определения для избыточных парциальных энтальпии и энтропии:
|
|
|
|
iизб(T , x) ≡ |
|
|
|
HiM (T , x); |
(4.119) |
H |
|
|
|
||||||
|
|
изб(T , x) = |
|
|
M (T , x) + R ln(x ) . |
(4.120) |
|||
S |
S |
||||||||
|
i |
|
i |
i |
|
||||
Перепишем теперь парциальный молярный потенциал Гиббса через стандартную составляющую, энергию смешения (образования) идеального раствора и избыточную свободную энергию, появляющуюся в реальном растворе:
|
|
(T , x) |
= Go(T ) |
+ RT ln(x ) + |
|
|
изб(T , x) |
||
G |
G |
||||||||
|
i |
|
|
i |
i |
|
i |
||
|
|
|
стандартное |
добавка в |
избыточная часть |
||||
|
|
|
значение для |
идеальном |
для реального |
||||
|
|
|
чистого в-ва |
растворе |
раствора |
||||
и внесем величину |
|
|
изб под знак логарифма концентрации: |
||||||
G |
|||||||||
|
|
|
|
i |
|
|
|
|
|
117
|
|
|
|
|
|
|
изб(T, |
|
|
|
|
o |
G |
||||
G (T , x) = |
||||||||
G |
(T ) + RT ln exp |
i |
||||||
|
i |
i |
|
|
|
RT |
||
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
Введем обозначения:
x) |
|
|
xi . (4.121) |
|
|
|
|
|
|
|
|
|
изб(T , x) |
γi(T, x) = |
G |
||
exp |
i |
||
|
|
|
RT |
|
|
|
|
|
|
|
изб |
(T , x) −T |
|
|
|
||
|
=exp |
H i |
||
|
|
|
|
RT |
|
|
|
|
|
|
изб |
|
|
|
|
||
Si |
(T , x) |
; |
|
|
|
|
|
|
|
|
|
|
|
|
|
ai(T, x) = γixi |
(4.122) |
и назовем величину ai(T, x) термодинамической активностью,
или просто активностью компонента, а величину γi(T, x) – коэф-
фициентом активности. Тогда (4.121) запишется в виде:
|
|
(T, x) = |
Go(T ) + RT ln[γ |
(T , x) x ] = |
|
|
G |
|
|||||
|
i |
i |
i |
i |
(4.123) |
|
= Go(T ) + RT ln[a |
(T , x)], |
|
||||
|
|
|||||
|
|
i |
i |
|
|
|
совпадающем, по форме записи, с выражением (4.109) для парциального молярного потенциала Гиббса вещества в идеальном растворе, если в последнем мольные доли xi заменить на активности ai. С формальных позиций это означает, что за стандартное состояние вещества в реальном растворе выбрано его состояние в виде чистого вещества, а все описание неидеальности отнесено на счет активности а – неизвестной пока функции температуры Т и вектора концентраций х.
Поскольку по физическому смыслу парциальный молярный потенциал Гиббса есть не что иное, как химический потенциал, умноженный на число Авогадро, то можем записать:
μi (T , x) = μoi (T ) + kT ln(ai ) = μoi (T ) + kT ln(γi xi ) . (4.124)
Новые термодинамические функции – активности и коэффициенты активности, на которые в прикладных расчетах возлагается вся ответственность за описание неидеальности растворов, были введены Дж. Н. Льюисом. С их помощью формальное математическое изложение теории растворов общего типа становится столь же простым в записи, как и в теории идеальных растворов. Конечно, для осуществления практических расчетов, в рамках классической термодинамики требуется экспериментальное определение активностей и их зависимостей от температуры и концентраций.
118
Первый метод экспериментального определения активностей, предложенный самим Льюисом и широко используемый на практике до настоящего времени, основан на измерении давления насыщенных паров. Температурную зависимость давления насыщенного пара над идеальным раствором можно рассчитать, приравнивая значения парциальных молярных потенциалов Гиббса вещества в конденсированном и газообразном состояниях:
|
|
(s) (T , x) |
= |
Go(s) (T ) + RT ln(x ) |
= |
||||
G |
|||||||||
|
i |
|
|
i |
i |
|
|||
= |
|
( g ) (T , p |
) |
= |
Go(g) (T ) + RT ln( pid ). |
||||
G |
|||||||||
|
|
|
i |
i |
|
|
i |
i |
|
Откуда, для случая идеальных конденсированного и газового растворов:
ln( piid ) = |
Go(s) (T ) − |
Go(g) (T ) |
|
|
i |
i |
+ ln(xi ) |
[атм]. (4.125) |
|
RT |
|
|||
|
|
|
|
По Льюису, формальная температурная зависимость давления насыщенного пара над неидеальным раствором сразу получается из
(4.125):
ln( pireal ) = |
Go(s) (T ) − |
Go(g) (T ) |
|
i |
i |
+ ln(ai ) [атм], (4.126) |
|
RT |
|
||
|
|
|
куда вместо концентрации входит уже активность ai(T, x). Заметим теперь, что в (4.125) и (4.126) выражение
Go(s) (T ) − |
Go(g) (T ) |
o |
i |
i |
= ln[ pi (T )] |
RT |
|
|
|
|
определяет зависимость давления насыщенного пара pio компонента i над чистым веществом, для которого а = х ≡ 1. Поэтому для давления насыщенного пара pireal над реальным раствором из
(4.126) имеем:
ln[ pireal (T , x)] = |
ln[ pio(T )] + ln[ai (T , x)], |
|
pireal (T , x)] |
= ai (T , x) pio(T ). |
(4.127) |
Или, через коэффициенты активности:
pireal (T , x)] = γi (T , x) xi ·pio(T ). |
(4.128) |
Откуда для расчета температурной и концентрационной зависимости активности по сериям экспериментов с чистым веществом и с растворами из (4.127) получаем:
119
|
preal (T , x) |
|
|
||
ai (T , x) = |
i |
|
|
, |
(4.129) |
|
o |
(T ) |
|||
|
|
pi |
|
|
|
или из (4.128) для коэффициентов активности:
|
preal (T , x) |
|
|
||
γi (T , x) = |
i |
|
|
. |
(4.130) |
|
o |
(T ) |
|||
|
xi pi |
|
|
||
Сегодня для определения активностей используют и другие экспериментальные методы, например, очень широко – электрохимические. Для теоретического установления вида зависимости коэффициента активности от концентраций и температуры необходимо выходить за рамки термодинамического метода и использовать методы статистической и квантовой физики с привлечением понятий об атомной структуре материи.
4.6.2. Уравнение Гиббса–Дюгема
Для термодинамических активностей, вместо простого уравнения связи, имеющегося для концентраций (Σxi = 1), получаются более сложные соотношения, выводимые из определения потенциала Гиббса для смесей. Полный дифференциал потенциала Гиббса в
его естественных переменных T, p, n записывается как |
|
dG = – SdT + Vdp + Σμidni , |
(4.131) |
где символами химических потенциалов μi обозначены парциальные молярные потенциалы Гиббса – соответствующие частные производные:
∂
μ= .
i∂ni T , p,n j≠i
Сдругой стороны, поскольку потенциал Гиббса есть аддитивная функция по отношению к количеству веществ:G
|
∂G |
|
|
|
G = Σni· |
|
= Σniμi , |
|
|
|
|
|||
|
|
|
|
|
|
∂ni T , p,n j≠i |
|
||
то его полный дифференциал можно записать как |
|
|||
dG = Σnidμi |
+ Σμidni . |
(4.132) |
||
|
120 |
|
|
|
