Калин Физическое материаловедение Том 2 Основы материаловедения 2007
.pdfПодставляя (4.132) в (4.131), получаем уравнение, называемое
уравнением Гиббса– Дюгема: |
|
SdT – Vdp + Σnidμi = 0. |
(4.133) |
Оно определяет связи между количествами веществ ni и измене-
ниями химических потенциалов dμi |
в изобарно-изотермических |
процессах (в которых dT и dp равны нулю): |
|
Σnidμi = 0. |
(4.134) |
Поделив в (4.134) каждое слагаемое на Σni, получаем уравнение связи уже между концентрациями веществ и изменениями химиче-
ских потенциалов в изобарно-изотермических процессах: |
|
Σxidμi = 0. |
(4.135) |
Поделив теперь (4.135) на любое dxj и переходя к частным производным, имеем:
|
∂μ |
i |
|
|
|
Σx |
|
|
= 0 . |
(4.136) |
|
∂x |
|
||||
i |
|
|
|
|
|
|
|
j p,T |
|
|
|
Для бинарной смеси (4.136) запишется, например, так:
|
∂μ |
|
|
|
∂μ |
2 |
|
|
|
1 |
|
+ x2 |
|
|
|
= 0 , |
|
x1 |
|
|
|
∂x2 |
|
|||
|
∂x2 p,T |
|
|
p,T |
|
|||
или, поскольку в ней |
х1 = 1– х2 : |
|
|
|
|
|
||
|
∂μ |
|
|
∂μ |
2 |
|
|
|
|
1 |
|
|
|
|
= 0 . |
(4.137) |
|
(1− x2 ) |
|
|
+ x2 |
∂x2 |
|
|||
|
∂x2 p,T |
|
p,T |
|
|
|||
Дифференцируя (4.124) по концентрации второго компонента, для бинарной смеси имеем:
|
∂μ |
1 |
|
|
|
|
|
|
∂x2 |
|
|
|
p,T |
||
|
∂μ |
2 |
|
|
|
|
|
|
∂x2 |
|
|
|
p,T |
||
=kT
=kT
|
∂ln[γ1 (T , x)] |
|
||
|
|
|||
|
∂x |
2 |
|
|
|
|
p,T |
||
|
|
|
|
|
|
∂ln[γ2 (T , x)] |
|||
|
∂x |
2 |
|
|
|
|
|
p,T |
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
−1− x |
|
|
|||||
2 |
; |
|
|||||
|
|
|
|
|
|
(4.138) |
|
|
1 |
|
|
||||
|
|
||||||
|
|
|
|||||
+ x |
|
|
|
||||
2 |
. |
|
|
||||
|
|
|
|
|
|
||
Подставляя (4.138) в (4.137), получаем выражение, которое следует использовать в прикладных расчетах:
121
|
|
|
∂ ln[γ |
1 |
(T , x)] |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
∂ ln[γ |
2 |
(T , x)] |
|||||||||
(1− x2 ) |
|
|
|
|
|
|
|
|
|
|
− |
|
|
|
|
|
+ x2 |
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
∂x |
|
|
|
|
|
|
1 |
− x |
|
|
|
|
|
∂x |
|
|
|
|
|||||||||
|
|
|
|
2 |
|
p,T |
|
2 |
|
|
|
|
|
|
2 |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
∂ ln[γ |
1 |
(T , x |
, x |
2 |
)] |
|
|
|
|
|
|
|
∂ ln[γ |
2 |
(T , x |
|
, x |
2 |
)] |
||||||||
= (1− x2 ) |
|
|
|
1 |
|
|
|
|
|
+ x2 |
|
|
|
|
|
1 |
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
∂x2 |
|
|
|
|
|
|
|
|
|
|
|
∂x2 |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
p,T |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
1 |
|
||
|
= |
||||
|
+ |
||||
|
|||||
|
|
x |
|
|
|
p,T |
|
2 |
|||
|
|
|
|||
|
|
|
|
|
|
|
= 0, |
|
|||
|
|
||||
p,T |
|
|
|
|
|
или, оставляя как неизвестную концентрацию только второго компонента:
|
∂ln[γ1(T , x2 )] |
|
|
|
|
(1− x2 ) |
∂x2 |
|
|
|
|
|
∂ln[γ2 (T , x2 )] |
|
|
|
+ x2 |
|
|
= 0, |
(4.139) |
|
|
∂x2 |
|
|||
p,T |
|
p,T |
|
|
получаем уравнение связи между двумя коэффициентами активностей γ1 , γ2 через одну концентрацию x2. Формулу (4.139) часто так же называют уравнением Гиббса – Дюгема, хотя она менее общая,
чем (4.133).
На этом возможности феноменологической термодинамики по формальному математическому описанию равновесий в системах с реальными растворами исчерпываются – в ней нет и не может быть рецептов для теоретического определения зависимости активностей от температуры и концентраций.
4.6.3.Статистическая трактовка энергетики растворов
Для того, чтобы численно исследовать равновесие в системах с реальными растворами, необходимо иметь практическую возможность вычислять парциальные молярные потенциалы Гиббса по уравнению (4.119):
|
G |
(T, x) = |
Go(T ) + RT ln(a |
) = |
Go(T ) + RT ln(γ |
i |
x |
) . |
|
|
i |
i |
i |
|
i |
i |
|
||
Входящие в него коэффициенты активности γi или предварительно определяются из эксперимента, или рассчитываются по формуле
(4.120):
|
|
|
|
изб |
|
|
|
|
|
|
|
|
изб |
|
|
|
|
|
изб |
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
Gi |
|
(T , x) |
|
|
|
|
|
H i |
(T , x) − T S i |
(T , x) |
||||||||
γi(T, x) = exp |
|
|
RT |
|
= exp |
|
|
|
|
RT |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и где |
|
изб(T , x) |
≡ |
|
|
M (T , x); |
|
|
изб |
(T , x) |
= |
|
M |
(T , x) + R ln(x ) . |
|||||||
|
|
|
S |
S |
|||||||||||||||||
H |
H |
||||||||||||||||||||
|
i |
|
|
|
|
i |
|
|
i |
|
|
|
|
i |
|
|
|
i |
|||
|
|
|
|
|
|
|
|
|
122 |
|
|
|
|
|
|
|
|
|
|
||
Предложим теперь статистическую трактовку избыточных ве-
личин H iизб, S iизб , фактически определяющих коэффициенты
активности индивидуальных веществ. Для этого воспользуемся простейшим физическим представлением о конденсированных растворах на основе решеточной модели с учетом только парного взаимодействия и только ближайших частиц – так называемый ква-
зихимический подход.
Под решеточной моделью понимается, что движение частиц сводится к колебаниям вокруг некоторых фиксированных равновесных положений – узлов кристаллической решетки, упорядоченно расположенных в любой точке пространства (существует дальний порядок структуры вещества). Это допущение абсолютно верно для твердых растворов кристаллических тел, но является достаточно спорным для жидкостей. Однако ближний порядок приблизительно одинаков как в твердых телах, так и в жидкостях. Основное различие этих фаз в том, что в кристаллах число ближайших соседей – координационное число z – имеет четко определенную величину (см. п. 1.4.1), в то время как для жидкостей оно постоянно меняется во времени и можно говорить лишь о некоторых средних значениях (подробнее см. п. 6.2.2). Причем z всегда меньше (и уменьшается с температурой) для жидкости из-за ее более «рыхлой» структуры, чем фиксированное для того же твердого вещества. Считается, что решеточная модель применима к жидкому состоянию вдали от критической точки.
При рассмотрении взаимодействия двух частиц зависимость потенциальной энергии ε каждой частицы от расстояния r обычно представляют в виде функции с двумя масштабными коэффициентами: энергетическим ε° и геометрическим σ:
ε(r) = ε° ϕ(σ/r).
Исходя из такого представления, можно получить (что и будет сделано ниже, в разделах 4.6.4–4.6.7) основные термодинамические соотношения для бинарных конденсированных растворов, по свойствам не слишком сильно отличающихся от идеальных, но достаточно часто встречающихся на практике. Это – растворы с полной (неограниченной) взаимной растворимостью компонентов во всем
123
интервале концентраций. В данную категорию растворов входит подавляющее большинство смесей реальных газов и жидкостейнеэлектролитов, а так же многие твердые растворы замещения.
Обсудим физические аспекты принадлежности к растворам с неограниченной растворимостью, т.е. достаточно близким к идеальным. В этих реальных растворах не появляется значительного дополнительного силового воздействия на частицы веществ по сравнению с их чистым состоянием, что выполняется при условиях:
• геометрических:
а) диаметры атомов dA и dB компонентов А и В примерно одинаковы; есть эмпирическое правило Юм-Розери, что неограниченная растворимость наблюдается в твердых растворах при различии в диаметрах молекул не более 15%; в жидкостях это размерное ограничение значительно более слабое; б) для твердых веществ одинаковы сингонии кристаллических решеток;
• энергетических:
а) энергии парного взаимодействия εАВ , εАА и εВВ близки между собой по величине, и поэтому значение интегральной энтальпии смешения HM относительно невелико (| HM| < RT);
б) дополнительно предполагается близость электронных структур компонентов, чтобы вид силовых функций ϕ парного взаимодействия был примерно одинаков.
Выполнение приведенных условий означает, что во всем диапазоне концентраций образуется твердый раствор замещения, не имеющий существенных искажений кристаллической решетки и заметных внутренних механических напряжений.
В ряде твердых растворов с HM < 0, при достаточно низких температурах (| HM| > RT), устанавливается упорядочение распределения атомов каждого сорта в кристаллической решетке по узлам определенного типа. При этом между одинаковыми атомами начинает наблюдаться дальний порядок их расположения. Из-за возрастания симметрии в этом сплаве возрастает число связей типа А-В по сравнению с хаотическим распределением атомов. Сплав, находящийся в таком состоянии, называется упорядоченным. При по-
124
вышении температуры (| HM| > RT) наблюдается переход части атомов со «своих» узлов на «чужие», и концентрации атомов на узлах различных типов становятся одинаковыми, т.е. сплав становится неупорядоченным в смысле отсутствия в нем дальнего порядка в расположении разнородных атомов. Дальний порядок узлов решетки, естественно, сохраняется.
Сравнительно небольшие абсолютные величины относительных различий в энергиях взаимодействия однородных (εАА, εВВ) и разнородных (εАВ) частиц позволяют предполагать, что в таком растворе не существует предпочтительности одного из видов пар связей типа A-A, B-B или A-B, и атомы обоих компонентов случайным образом расположены в упорядоченном пространстве узлов решетки. Такие растворы с хаотическим, не имеющим дальнего порядка расположения атомов примеси называются неупорядоченными твердыми растворами. Неупорядоченность и практическое отсутствие дальнодействия дает основания учитывать силы только парного взаимодействия и только соседних атомов. Тогда молярная энтальпия раствора может рассматриваться как простейшая, линейная,
функция чисел пар связей типа A-A, B-B и A-B: |
|
Hm = PAA·εAA + PBB·εBB + PAB·εAB , |
(4.140) |
где PAA, PBB и PAB – числа пар атомов, дающих соответственно ком-
бинации A-A, B-B, A-B.
Несмотря на упрощенность модели, такой подход дает возможность удовлетворительного объяснения ряда явлений, происходящих в сплавах металлов. В частности, он позволяет наглядно пояснить появление основных видов диаграмм состояния. Более того, оказалось, что, несмотря на учет взаимодействия только ближайших соседей (сил ближнего порядка), он позволяет полуколичественно описать и явление упорядочения (установления дальнего порядка) в расположении атомов разнородных компонентов в узлах кристаллической решетки в некоторых твердых растворах. Полностью упорядоченное состояние соответствует определенным стехиометрическим составам сплавов с целочисленным соотношением компонентов и является своеобразным прообразом химического соединения с малой энергией образования. Поэтому рассмотрение только парного взаимодействия и только соседних атомов и назы-
125
вают квазихимическим подходом. Учет электронных сил только ближнего порядка и пренебрежение силами деформации из-за разности размеров атомов в основном определяют ограничения квазихимического подхода.
Энергии межатомной связи являются фундаментальными характеристиками веществ. Для простоты будем предполагать, что они не зависят от концентраций и температуры. Энергиям связей в (4.140) придадим в соответствии с термодинамическим правилом знаков отрицательные значения: знак «минус» указывает, что, при переходе от газообразного состояния с нулевой потенциальной энергией взаимодействия частиц к конденсированному энергия выделяется в окружающую среду. Заметим так же, что энергии парного взаимодействия εAA и εВВ можно грубо оценить именно из энергии испарения или сублимации как энергии разрыва всех связей. Пусть z – координационное число, NA – число Авогадро. Тогда:
εLAA = − |
|
|
H Aисп |
; |
εSAA = − |
|
|
H Aсубл |
; |
|
1/2 |
z NA |
1/2 |
z NA |
|||||||
|
|
|
|
|||||||
εBBL = − |
|
|
HBисп |
; |
εSBB = − |
|
|
HBсубл |
. |
|
1/2 |
z NA |
1/2 |
z NA |
|||||||
|
|
|
|
|||||||
Множитель «1/2» появляется в силу того, что связи одинаковых атомов типа А′-А′′ подсчитываются дважды: сначала из позиции А′, затем из позиции А′′. Это необходимо делать, поскольку положения А′-А′′ и А′′-А′, являясь энергетически эквивалентными, представляют собой разные позиции в пространстве и дают тем самым вклад в число возможных квантовых состояний системы – в пози-
ционную составляющую энтропии.
В табл. 4.4 приведены физические характеристики некоторых переходных металлов, сплавы которых широко используются в практике ядерного и ракетно-космического комплексов. Серым фоном в ней выделены металлы, которые при высоких температурах в принципе могут давать ( a/a < 15%), и во многих случаях
(Zr–U, Zr–Hf, Ta–Nb, Ta–W, Ta–Mo, Nb–W, Nb–Mo и др.) дают ме-
жду собой неограниченные твердые растворы. В то же время сплав Th–Cr ( a/a ≈ 35%) в твердом состоянии дает полное отсутствие взаимной растворимости.
126
Из табл. 4.4 и рис. 4.17 видно так же, что существует определенная корреляция между температурой плавления и энтальпиями плавления и сублимации (за исключением хрома).
|
|
|
|
|
Таблица 4.4 |
Физические свойства некоторых ОЦК фаз переходных металлов |
|||||
|
|
|
|
|
|
Металл |
Параметр |
Температура |
HTпл |
|
H298субл |
решетки |
плавления |
пл |
|
|
|
|
а, нм |
К |
|
кДж/моль |
|
β-Th |
0,4110 |
2023 |
19,2 |
|
598 |
β-Zr |
0,3616 |
2125 |
19,7 |
|
603 |
β-Hf |
0,3620 |
2503 |
25,5 |
|
619 |
γ-U |
0,3524 |
1405 |
15,5 |
|
480 |
Ta |
0,3303 |
3269 |
31,7 |
|
782 |
Nb |
0,3299 |
2741 |
27,2 |
|
729 |
W |
0,3165 |
3680 |
33,5 |
|
871 |
Mo |
0,3147 |
2890 |
27,6 |
|
659 |
Cr |
0,2885 |
2130 |
15,3 |
|
397 |
Вводя числа атомов NA и NB видов A и B в одном моле раствора (NA + NB = NA = 6,023·1023 – число Авогадро), имеем:
2PAA + PAB = z·NA ,
2PBB + PAB = z·NB .
Исключая из (4.140) величины PAA и PBB с помощью последних двух уравнений, получаем для интегральной молярной энтальпии раствора:
Hm = 1/2z·NA·εAA + 1/2z·NB ·εBB + PAB ·(εAB – 1/2 εAA – 1/2εBB). (4.141)
Очевидно, что энтальпия H Ao одного моля чистого вещества А
равна 1/2z·NA·εAA; для одного моля В имеем: H Bo = 1/2z·NA·εВВ. Так же очевидно, что хA= NA/NA, xB= NВ/NA. Тогда имеем:
для энтальпии NA атомов чистого A:
1/2z·NA·εAA = хA
для энтальпия NB атомов чистого B:
1/2z·NB·εBB = хВ
Отсюда для величины интегральной молярной теплоты смеше-
ния раствора
127
а
б
Рис. 4.17. Зависимость энтальпии плавления (а) и сублимации (б) от температуры плавления металлов
H M = H |
m |
−Σx |
H o |
|
i |
i |
(из энтальпии одного моля реального раствора вычитаются энтальпии чистых компонентов в соответствующих количествах), из
(4.141) получаем:
HM = PAB ·[εAB – 1/2 (εAA + εBB)]. |
(4.142) |
Следовательно, в квазихимическом приближении теплота смешения раствора заданного состава определяется числом пар связи PAB разнородных атомов.
В одном моле раствора имеется всего NA= 6,023·1023 позиций. Из них xA·NA позиций заняты атомами вида A. Это означает, что вероятность нахождения атома A в произвольно выбранном узле равна молярной доле хA. Если рассматривается конкретный атом вида A, то вероятность нахождения атома вида B в произвольной
128
позиции будет равна молярной доле хB. Вероятность такой комбинации, при которой один атом A находится в соседстве с одним атомом B (оба атома – уже в определенных позициях), равна поэтому хA xB. Тогда из 1/2 z NA связей в 1 моле полное число связей PAB между разнородными атомами будет
PAB = 2хA xB (1/2 z NA) = хA xB z NA |
(4.143) |
и |
|
HM = хA·xB·z·NA ·[εAB – 1/2 (εAA + εBB)]. |
(4.144) |
«Двойка» в (4.143) появляется в силу того, что позиции А–В и В–А статистически равноправны и должны подсчитываться раздельно.
Введем для удобства величину Ω, определяемую уравнением
Ω = z·NA·[εAB – 1/2(εAA + εBB)], |
(4.145) |
не зависящую от концентраций и называемую параметром взаимодействия (квазихимическим). Тогда интегральная молярная теплота смешения для одного моля раствора запишется как
HM = хA·xB·Ω = (1– хB)·xB·Ω . |
(4.146) |
Понятие «параметр взаимодействия» очень широко используется на практике; величины параметров взаимодействия определяются экспериментально и для бинарных систем табулируются в справочниках.
4.6.4.Совершенные растворы
Вконденсированном веществе часто предполагается, что потен-
циальная энергия εi частицы сорта i есть сумма ее парных взаимодействий εij с ближайшими z соседями различных сортов j, находящихся на расстояниях rij от нее:
Z
εi = ∑εij (rij ) . j=1
Полную потенциальную энергию частицы ε обычно представляют в виде функции с двумя масштабными коэффициентами ε° и r°:
ε(r) = ε° ϕ(r°/r),
и трактуют следующим образом: потенциальная энергия взаимодействия принимает минимальное (причем отрицательное) значение ε° при минимально возможном равновесном расстоянии r° ме-
129
жду центрами частиц. Символом «ϕ« обозначена некоторая силовая функция, в простейших случаях допускающая представление в виде двух степенных членов:
r o n |
r o m |
|
||||
ϕ = −α |
|
|
+ β |
|
|
, |
|
|
|||||
|
r |
|
|
r |
|
|
|
|
|
|
|
||
где α и β – константы, характеризующие соответственно силы притяжения и отталкивания; n, m – постоянные, характерные для взаимодействия данного типа.
Рассмотрим с этих простейших позиций модели специальных случаев растворов.
Идеальные растворы. Предположим, что законы притяжения и отталкивания частиц одинаковы для как для взаимодействия А-А, так и В-В и А-В:
ϕАА = ϕВВ = ϕАВ. |
(4.147) |
Для выполнения (4.147) из предидущих формул очевидно тре-
бование равенства размеров всех частиц: |
|
||
|
|
rºАА = rºВВ , |
(4.148) |
|
|
а для твердых тел неявно и требова- |
|
|
|
||
|
|
ние одинаковости сингоний. |
|
|
|
Допустим теперь, что и электрон- |
|
|
|
ные структуры компонентов, а с ни- |
|
|
|
ми и все парные энергии взаимодей- |
|
|
|
ствия частиц в растворе одинаковы: |
|
|
|
εºAA = εºBB = εºАB. Тогда энергии взаи- |
|
|
|
модействия всех частиц при любых |
|
|
Рис. 4.18. Зависимость |
расстояниях между частицами, т.е. |
|
|
энергий взаимодействия |
при всех температурах, будут одина- |
|
|
в идеальных растворах |
ковы: |
|
|
от межатомного расстояния: |
εAA = εBB = εАB , |
(4.149) |
|
εАА(r) = εBB(r) = εАВ(r) |
что иллюстрирует рис. 4.18. |
|
|
|
|
|
Условия (4.147)–(4.149) соответствуют формальному термодинамическому определению идеальных растворов (см. п. 4.5.3). По-
скольку из (4.149) видно, что εAB = ½(εAA + εBB), из (4.145), (4.146)
сразу формально математически следует: параметр взаимодействия
130
