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

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

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

Коэффициент активности в предельно разбавленных растворах для растворенного компонента при хi 0 стремится к некоторому конечному значению, отличному от единицы:

 

 

 

 

xi

)2

Ω

i

 

 

Ω

i

 

 

lim (γi ) =

lim exp (1

 

=

exp

 

 

1.

 

 

 

 

xi 0

xi

0

 

 

 

 

 

 

RT

 

RT

 

Подставляя (4.168) в (4.128), получаем:

 

 

 

 

 

pireal (T , xi )]

=

 

Ω

i

 

 

xi ·pio (T ) = K Г·xi ·pio (T ).

exp

 

 

 

RT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.169)

(4.170)

Как видим, и для растворенного компонента в предельно разбавленном растворе давление насыщенного пара подчиняется закону прямой пропорциональной зависимости от его мольной доли и давления пара над чистым компонентом, но уже с некоторым коэффициентом KГ, отличным от единицы. Эту закономерность феноменологически установил Генри; уравнение (4.170) называется

законом Генри, а коэффициент KГ =

 

Ω

exp

 

– константой Ген-

 

 

 

RT

ри. Мы здесь выяснили ее физическое содержание, из которого видно, что константа Генри вовсе не является константой, а зависит от температуры.

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

вается областью Генри, а раствор – предельно разбавленным.

Для умеренной положительной неидеальности (Ω/RT < 2 и компоненты растворяются неограниченно) граница применимости закона Генри с точностью лучше 1% оценивается как 0,25–0,50 мол.%. При более высокой положительной неидеальности, когда компоненты ограниченно растворимы (Ω/RT > 2), ширина области

141

Генри (предельно разбавленных растворов) уменьшается. И, наоборот, в растворах, более близких по свойствам веществ, область Генри расширяется. Например, для Ω/RT = 0,16 область предельно разбавленных растворов оценивается как диапазон концентраций до 3 мол.%. Заметим также, что выполнимость закона Рауля всегда наблюдается в большей области концентраций, чем закона Генри

(см. рис. 4.22).

4.6.7. Субрегулярные и квазирегулярные растворы

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

Например, энтропия колебаний твердого раствора не является в точности линейной функцией состава, но это второстепенный фактор; из всех допущений наиболее уязвимым является предположение о равенстве параметров решеток всех компонентов раствора. Разница размеров частиц не имеют существенного значения для непрерывных твердых растворов (растворов с неограниченной растворимостью), поскольку они могут существовать только при условии, что параметры решетки чистых металлов различаются незначительно (<15%).

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

142

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

Рассмотрим влияние размерного фактора на свойства жидких растворов. Прежде всего, заметим, что введение в жидкость частицы другого размера не вызовет упругих напряжений – частицы жидкости легко, без энергетического барьера, сдвигаются относительно друг друга. Так как в жидкости не все атомы находятся в минимальном энергетическом положении, которое соответствует максимальному числу ближайших соседей, энергия жидкости по абсолютной величине всегда больше, чем энергия кристалла; на плавление всегда надо подвести энергию извне. Но по модулю она меньше, чем у твердых тел, и поэтому меньше взаимное силовое влияние частиц жидкости друг на друга. Кроме того, плотность жидкостей обычно примерно на 10% меньше, а расстояния между частицами больше, чем в кристаллах, что еще уменьшает энергию взаимодействия разнородных частиц.

Все это приводит к тому, что интегральная энтальпия смешения HM, а с ней и парциальные молярные энтальпии HiM , и избыточные парциальные энтальпии Hiизб в жидкостях очень часто близки к нулю. Модель регулярных растворов в этом случае вырождается в модель совершенных (Ω = 0). Но размерный фактор существенно сказывается на энтропии смешения SM, и модель совер-

шенных растворов для жидкостей не всегда адекватна.

Пусть, например, атомы B значительно больше, чем атомы A. Тогда очевидно, что число ближайших соседей для разбавленного раствора металла B в металле A будет больше числа ближайших соседей для разбавленного раствора металла A в металле B. Поэто-

143

му число ближайших соседей (координационное число z) в таких растворах для каждой частицы изменяется на протяжении всего интервала концентраций. Эти позиционные изменения находят отражение в появлении у парциальной молярной энтропии смешения SiМ (T , x) избыточной составляющей Siизб(T , x) над позицион-

ной компонентой, равной R ln(xi ) для идеальных, совершенных и регулярных растворов:

 

 

М (T , x) = −R ln(x ) +

 

изб(T , x) ,

S

S

 

i

i

i

у которых координационное число не меняется с составом.

Следовательно, в растворе из частиц разных размеров коэффициенты активности должны отличаться от единицы, даже если интегральная молярная теплота смешения HM равна нулю. Действительно, из самого общего термодинамического соотношения (4.122), для бинарного раствора имеем:

 

 

изб

(T , x) T

 

 

изб

 

 

 

 

 

 

 

 

изб

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

H A

 

S A

(T , x)

 

 

S A

 

(T , x)

;

γA = exp

 

 

RT

 

 

 

 

 

 

 

= exp

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

изб

 

 

 

 

 

 

 

 

(4.171)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

γB = exp

 

S B

 

(T , x)

.

 

 

 

 

 

 

 

 

 

 

 

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

ния в степенной ряд по малому параметру:

 

ln(γA) = α1·(1– xA) + 1/2·α2·(1– xA)2 + 1/3·α3·(1– xA)3 +

 

+ 1/4·α4·(1– xA)4

+...

(4.172)

ln(γB) = β1·(1– xB) + 1/2·β2·(1– xB)2 + 1/3·β3·(1– xB)3 +

 

+ 1/4·β4·(1– xB)4

+…

(4.173)

где коэффициенты α, β – некоторые функции температуры. Поскольку существует уравнение Гиббса–Дюгема (4.139), то

между коэффициентами αi , βi имеется определенная связь. Под-

144

ставляя (4.172) и (4.173) в (4.139) и выполнив соответствующее дифференцирование, получим:

α1 = β1 = 0; β2 = α2 + α3 + α4 + α5 + ...; β3 = − α3 − 2α4 − 3α5 + ...

Если в (4.172), (4.173) ограничиться вторым, ненулевым, членом разложения, имеем: β2 = α2.

Положив, что α2 = β2 = 2Ω/RT, получим выражение, соответствующее формуле (4.163), выведенной выше для коэффициентов активности по модели регулярных растворов:

 

Ω (1x

 

 

)2

 

 

 

изб

 

 

ln(γA ) =

 

 

H

;

RT

A

 

=

 

 

 

A

 

 

 

 

 

RT

(4.174)

 

Ω (1x

 

 

)2

 

 

 

 

изб

 

 

 

 

 

H

 

 

ln(γB ) =

 

 

B

=

 

 

 

B

 

,

RT

 

 

RT

 

 

 

 

 

 

 

 

 

в которых энергетический параметр Ω считается не зависящим от температуры. Если теперь в (4.171) положить

 

 

 

избA (T , x)

 

1

 

 

 

 

избB (T , x)

 

1

 

 

S

=

α2 (1xA )2 ,

 

S

=

α2 (1xB )2

 

 

R

2

 

 

R

2

 

 

 

 

 

 

 

 

 

 

и допустить, что α2 не зависит от температуры, то, вводя обозначе-

 

 

 

2

 

 

изб

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ние s

= α2

=

S

, получим выражения,

 

сходные по форме для

R(1

x)2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

регулярных растворов:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

изб

 

 

 

s (1x

 

 

)2

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

 

ln(γA ) = −

 

 

 

A

 

=

 

 

 

A

 

;

 

 

 

 

 

 

 

R

 

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(4.175)

 

 

 

 

 

 

 

 

 

 

 

 

изб

 

 

s (1

x

 

 

)2

 

 

 

 

 

 

 

 

 

 

S

 

 

 

 

 

 

 

 

 

 

ln(γB ) =−

 

 

 

B

 

=

 

 

 

B

 

 

,

 

 

 

 

 

 

 

R

 

R

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

но в которые входят величины, связанные уже не с энергетическими, а энтропийными характеристиками процесса растворения. Растворы, поведение которых адекватно описывается формулами (4.175), получили название субрегулярных. Из (4.175) видно, что в субрегулярных растворах коэффициенты активности, а с ними и активности компонентов не зависят от температуры.

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

145

вместе. Если оставаться в пределах точности сделанных выше допущений, то для более полного описания температурной и концентрационной зависимости коэффициентов активности разумно объединить формулы (4.174) и (4.175):

ln(γA ) = (Ω − sT ) (1xA )2 ; ln(γB ) = (Ω − sT ) (1xB )2 . (4.176) RT RT

Данная модель растворов получила название квазирегулярных. В моделях квазирегулярных растворов сохраняют форму записи (4.163), но параметр взаимодействия считают линейно зависящим

от температуры:

 

Ω = h – sT.

(4.177)

Физическая природа такой зависимости была раскрыта выше: h

– это энтальпийный вклад в неидеальность раствора, s – энтропийный. Феноменологические коэффициенты h и s определяются экспериментально.

4.7.Термодинамическая теория диаграмм состояния

4.7.1.Постановка задачи

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

Главным из явных условий термодинамического равновесия в неидеальных системах остается равенство химических потенциалов компонентов в разных фазах. По аналогии с идеальными системами (см. выражение (4.111)), запишем условие равновесия твердой и жидкой фаз через парциальные молярные потенциалы Гиббса. Для этого, согласно рекомендации Льюиса, заменим в (4.111) концентрации на активности:

GAoS (T ) + RT ln(aAS )

=

GAoL (T ) + RT ln(aAL );

(4.178)

GBoS (T ) + RT ln(aBS )

=

GBoL (T ) + RT ln(aBL ).

 

 

146

 

 

Для вычисления активностей согласно теории регулярных растворов имеем уравнения:

S

S

 

(1

xS )2 Ω

S

 

 

 

 

 

A

 

 

 

;

aA = xA exp

 

 

RT

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

S

 

 

 

S

2

 

 

 

 

 

(1xB )

 

ΩS

 

aB

= xB exp

 

 

 

 

 

 

;

 

 

RT

 

 

 

 

 

 

 

 

 

 

 

 

 

L

L

 

 

L

2

 

 

 

 

(1xA )

 

 

ΩL

aA

= xA exp

 

 

 

 

 

 

;

 

RT

 

 

 

 

 

 

 

 

 

 

 

 

L

L

 

 

L

 

2

 

 

 

 

 

(1xB )

 

ΩL

aB

= xB exp

 

 

 

 

 

 

.

 

RT

 

 

 

 

 

 

 

 

 

 

 

(4.179)

(4.180)

Воспользовавшись справедливыми и для неидеальных систем условиями нормировки концентраций (4.110):

 

 

 

 

xAL + xBL = 1;

xAS + xBS = 1,

 

подставляя их в (4.179) вместе с (4.180) в (4.178), имеем:

 

RT ln

1xBS

+( GAoL (T) GAoS (T)) = (xBS )2 ΩS (xBL )2 ΩL ;

 

 

 

 

 

1xBliq

 

 

(4.181)

RT ln

xBS

+(

GBoL (T) GBoS (T))= (1xBS )2 ΩS (1xBL )2 ΩL.

 

 

 

xBL

 

 

 

Как и в случае идеальных систем, получено два уравнения относительно трех неизвестных величин: xBS , xBL , Т. Следовательно, и в

неидеальных системах равновесие жидкость–твердое будет на диаграмме состояния в координатах {х,Т} описываться двумя линиями: ликвидусом и солидусом. Перепишем (4.112) в том же виде, что и (4.181), более типичном для компьютерных расчетов:

RT ln

1xBS

+ ( GAoL (T ) GAoS (T ))= 0;

 

1xBL

 

(4.182)

RT ln

xBS

 

+ (

GBoL (T ) GBoS (T ))= 0.

xBL

 

 

 

 

Как видим, в отличие от системы (4.182) правые части у (4.181) являются не нулем, а алгебраическими функциями от концентраций и зависят так же от параметров взаимодействия в обеих фазах. Поскольку в левой части системы (4.181) стоят логарифмические функции, составляющие ее уравнения трансцендентны. Следовательно, сама система для описания равновесия неидеальных растворов не только не имеет аналитических решений, но и при произвольно заданной температуре Т вообще стоит вопрос о сущест-

147

вовании и числе ее решений относительно xBS , xBL в положительных действительных числах, меньших единицы.

4.7.2. Основы графической термодинамики

До появления компьютерной техники для расчета диаграмм состояния применялись методы графической термодинамики. Она создана еще Гиббсом и основана на взаимосвязи интегральных и молярных термодинамических характеристик растворов. С развитием компьютерной техники для этой цели стало проще и эффективней сразу использовать численные методы решения системы уравнений (4.181) или эквивалентных ей для других моделей растворов. Однако графическая термодинамика в силу ее наглядности не потеряла своего значения и сегодня. Мы будем использовать графические методы не для строгого расчета, а для иллюстрации теоретических принципов, отражающих взаимоотношения равновесных фаз.

Основой методов графической термодинамики является совместное рассмотрение кривых концентрационной зависимости молярных потенциалов Гиббса для двух растворов, находящихся в равновесии. Зависимость потенциала Гиббса Gm одного моля всего раствора от концентрации x второго компонента В, согласно (4.98), описывается формулой:

Gm(р, Т, x) = x· G В (р, Т, x) + xА· G А(р, Т, x) =

 

 

 

 

 

 

 

 

 

 

 

 

= x· G В (р, Т, x) + (1– xG А(р, Т, x).

(4.183)

Поскольку

 

 

 

 

 

 

 

 

 

А(р, Т, x) =

GAo

(р, Т) + RTln(aA),

 

G

 

 

 

В(р, Т, x) =

GBo

(р, Т) + RTln(aB),

 

 

G

 

а для любого чистого компонента его активность равна единице, понятно, что кривая, соответствующая уравнению (4.183), отсекает на ординатах x = 0 (aA = 1) и x = 1 (aB = 1) величины:

Gm(р,Т, x=0) =

G

 

А (р,Т, x=0) =

GAo (р,Т) + RTln(1) =

GAo

(р,Т);

Gm(р,Т, x=1) =

 

В (р,Т, x=1) =

GBo (р,Т) + RTln(1) =

GBo

(р,Т),

G

 

148

 

 

равные значениям стандартных потенциалов чистых компонентов

А и В.

Далее, частная производная от выражения (4.183) по х, равная:

G

m

(x)

 

 

 

 

 

= G В(р,Т, x) – G А(р,Т, x) = tg(α(x)), (4.184)

 

 

 

 

x

p,T

по смыслу есть тангенс угла наклона касательной к кривой концентрационной зависимости молярного потенциала Гиббса раствора. Записывая само уравнение этой касательной f(x) в произвольной

точке с концентрацией x*B как уравнение прямой в отрезках:

f(x) = G А(р,Т, x*B )+x ·G В(р,Т, x*B ) – G А(р,Т, x*B )] (4.185)

видим, что касательная отсекает на ординатах x = 0 и x = 1 значения парциальных потенциалов Гиббса компонентов в этом рас-

творе x*B :

f(0) =

G

А(р,Т, xB* ); f(1) =

G

В(р,Т, xB* ).

(4.186)

Поскольку равенство именно парциальных энергий Гиббса всех одноименных компонентов в разных фазах является условием равновесия этих фаз, то из (4.185) и (4.186) следует, что графической интерпретацией этого равенства является наличие общей касательной у кривых концентрационной зависимости молярных потенциалов Гиббса растворов, образующих разные фазы. Визуальный поиск таких общих касательных и составляет сущность графической термодинамики.

Поясним сказанное выше численным расчетом простейшей бинарной системы, в которой оба раствора являются идеальными. Для иллюстрации порядка численных значений реальных величин за термохимические данные для компонента А примем свойства тантала; данные для компонента В пусть соответствуют β- цирконию. Исходные данные для чистых реальных компонентов при разных температурах приведены в табл. 4.5.

На рис. 4.24 для выбранной условной (виртуальной) системы представлены расчеты интегральной молярной энергии Гиббса Gm по модели идеального твердого раствора при температуре 2500 К. Заметим, что в действительности в системе Ta–Zr наблюдается су-

149

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

 

 

 

 

 

 

 

 

Таблица 4.5

 

Некоторые свойства тантала и циркония

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Пара-

Тплавл,

 

Энергия Гиббса

Go,

кДж/моль

 

Ме-

метр

 

 

 

 

 

 

 

2000 К

2500 К

 

3500 К

талл

решетки,

К

тв

ж

тв

ж

 

тв

 

ж

 

нм

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Та

0,3303

3269

–140

–128

– 189

– 181

 

–297

 

–299

 

 

 

 

 

 

 

 

 

 

 

Zr

0,3616

2125

–142

–141

– 193

– 196

 

–309

 

–319

 

 

 

 

 

 

 

 

 

 

 

Рис. 4.24. Энергетические характеристики:

1 – Gm виртуального твердого раствора Ta-Zr, принятого за идеальный; 2 – касательная к кривой 1

при х = 0,5

Рис. 4.24, прежде всего, иллюстрирует тот факт, что обсуждаемая кривая несимметрична – касательная в точке x = 0,5 не является горизонталью. Численные значения рассчитанных парциальных мольных величин (для х = 0,5) составляют:

GAs = –203,5 кДж/моль; GBs = –207,5 кДж/моль.

Далее, из рисунка 4.27 видно, что концентрационная зависимость молярных потенциалов Гиббса идеальных растворов описывается вогнутыми кривыми с единственным минимумом. Данный

150