Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Калин Физическое материаловедение Том 2 Основы материаловедения 2007

.pdf
Скачиваний:
1641
Добавлен:
16.08.2013
Размер:
8 Мб
Скачать

Подставляя (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)

(1x2 )

 

 

+ 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

 

 

 

 

 

 

 

 

 

 

 

 

1x

 

 

2

;

 

 

 

 

 

 

 

(4.138)

 

1

 

 

 

 

 

 

 

+ x

 

 

 

2

.

 

 

 

 

 

 

 

 

Подставляя (4.138) в (4.137), получаем выражение, которое следует использовать в прикладных расчетах:

121

 

 

 

ln[γ

1

(T , x)]

 

 

 

 

 

 

1

 

 

 

 

 

 

ln[γ

2

(T , x)]

(1x2 )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

+ x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

1

x

 

 

 

 

 

x

 

 

 

 

 

 

 

 

2

 

p,T

 

2

 

 

 

 

 

 

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ln[γ

1

(T , x

, x

2

)]

 

 

 

 

 

 

 

ln[γ

2

(T , x

 

, x

2

)]

= (1x2 )

 

 

 

1

 

 

 

 

 

+ x2

 

 

 

 

 

1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

x2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

p,T

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

=

 

+

 

 

 

x

 

 

p,T

 

2

 

 

 

 

 

 

 

 

 

= 0,

 

 

 

p,T

 

 

 

 

или, оставляя как неизвестную концентрацию только второго компонента:

 

ln[γ1(T , x2 )]

 

 

 

(1x2 )

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

H Ao ,
H Bo .

Из табл. 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– хBxB·Ω .

(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