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

книги из ГПНТБ / Хачатурян, А. Г. Теория фазовых превращений и структура твердых растворов

.pdf
Скачиваний:
11
Добавлен:
22.10.2023
Размер:
14.34 Mб
Скачать

двухфазного равновесия (5.1). В неравновесном состоянии nepéохлажденного однородного раствора оба минимума не равны друг

ДРУГУ (Рис23. а)-

Рис. 22. Формирование равновес­ ного распределения концентра­ ции в двухфазном состоянии. а) Зависимость термодинамическогом потенциала от состава. б) Фазовая траектория, отвеча­ ющая равновесному распределе­ нию концентрации, в) Равновесное распределение концентрации в

двухфазном состоянии.

Рис. 23. Формирование одномерного зародыша H O B oâ фазы и периодических

распределений концентрации, а) Зави­ симость термодинамического потенци­ ала от состава, б) Фазовые траектории, отвечающие одномерному зародышу новой фазы (случай а) и периодиче­ скому распределению концентрации (случай ß). в) Кривая а описывает распределение концентрации в крити­ ческом зародыше, кривая ß — в перио­ дическом распределении концентрации.

Отмеченные свойства функции / (с) — цс позволяют проана­

лизировать решения уравнения

(8.3)

90

следующего из (8.2). Решения уравнения (8.3) имеют вид

(8.4)

где х0 — постоянная интегрирования, играющая роль начала отсчета.

Все решения уравнения (8.3) зависят от двух параметров — химического потенциала р, и интеграла движения Е. Параметр р можно выразить через среднюю концентрацию в сплаве с по­ мощью уравнения сохранения числа атомов (7.7). В одномерном случае оно имеет вид

+ L/2

І (8-5)

-L,2

где L — размер кристалла в направлении оси х, с0 (х, р, Е) — решение уравнения (8.3). Уравнение (8.5) позволяет определить функцию р = р (Е , с) и, следовательно, записать решение урав­ нения (8.3) в форме с0(х) = с0(х; с, Е). Каждому решению урав­ нения (8.3), таким образом, отвечает некоторое стационар­ ное состояние системы. Анализ уравнения (8.3), проводимый методом фазовых траекторий, показывает, что существуют три типа решения. В случае ß (см. рис. 23, б и в ) решение уравнения (8.3) — функция с0 (х) — является периодической функцией ко­ ординат X . Период распределения а определяется из уравнения

(8.4) :

 

ßic

+

а §

= 2

V m - E + f(p ) - ] u > )

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

Наконец, третий случай, изображенный на рис. 22, в, описы­ вает две равновесные фазы и разделяющий их пограничный слой. Такое распределение концентрации отвечает равновесному со­ стоянию системы. Для того чтобы убедиться в этом, докажем, что термодинамический потенциал (7.9) принимает минимальное зна­ чение Ф = Фтіп при‘распределении концентрации, изображенном на рис. 22, в.

Чтобы вычислить термодинамический потенциал Ф любого экстремального распределения концентрации, удовлетворяющего уравнению (8.1), необходимо с помощью (8.3) исключить / (с) в вы­ ражении (7.9) и учесть одномерный характер распределения с (я).

91

При этом получим:

Ь/2

ф = , , (£ + т S т ) -

(8-7)

—Ь/2

Из выражения (8.7) следует, что термодинамический потенциал Ф принимает минимальное значение, если одновременно выполняются два условия: 1) Е принимает наименьшее возможное значение, равное Е 0, и 2) второе слагаемое в (8.7), которое, по определению, является положительной величиной, стремится к нулю при Z.-> оо. Из (8.2) следует, что

Е > / (с) — рс.

(8.8)

Следовательно, минимальное значение Е есть Е 0 =

If (с) — pc]mjn,

т. е. равно минимальному значению плотности термодинамическо­ го потенциала / (с) — рс. При Е = Е 0 распределение концентра­ ции с (х) — с0(х\ с, Е 0) имеет вид, изображенный на рис. 22, в. Распределение с0 (х ; с, Е 0) зависит от координат только в преде­ лах переходного слоя между двумя фазами. Поэтому интеграл во втором слагаемом в (8.7) имеет порядок

Ь/2

где D — толщина переходного слоя между фазами, а само второе слагаемое, следовательно, стремится к нулю при L —>оо. Таким образом, оказываются выполненными оба условия, которые не­ обходимы для того, чтобы функция с(х) = с0(х; с, Е 0) обеспечи­ вала абсолютный минимум термодинамического потенциала Ф. Минимальное значение термодинамического потенциала равно

Фтіп = Е0Ѵ = min [/ (с) — pc] • V.

(8.9)

Второе слагаемое в (8.7),

 

Ь/2

 

Es = S jj

(8.10)

—L/2

 

пропорционально площади межфазной границы S = VfL и по­ этому описывает поверхностную энергию. Мы не учли его в (8.9), так как поверхностная энергия дает асимптотически малый вклад по сравнению с объемной энергией. Формулу (8.10) для поверх­ ностной энергии можно переписать в более удобной форме:

 

_

Ь/2

Ь/2

с0

с0

 

Т =

-J =

^

 

^p d c = { /2 ß (—£ „ + /(с) —ре) de,

 

 

 

- L2ß^

с"

с“

(8.11)

где

у — коэффициент поверхностного

натяжения.

 

92

Мы приходим к выводу, что ни распределение, изображенное на рис. 23, в (случай а), ни распределение, изображенное на рис. 22, е, не могут описывать метастабильные состояния сплава (первое описывает седловое состояние, второе — состояние абсо­ лютного минимума свободной энергии). Следовательно, метаста­ бильные состояния могут быть связаны только с периодическими распределениями с0 (х\ с, Е), изображенными на рис. 23, в (слу­ чай ß) и являющимися решениями уравнения (8.3).

Впредыдущем параграфе было показано, что распределение концентрации

с(X) обеспечивает условный минимум свободной энергии или термодинами­ ческого потенциала, если квадратичная форма (7.13) положительно опре­

делена. В одномерном случае уравнение (7.13) можно переписать в форме

 

Ь /2

 

&F = 62Ф = S

бс (х) НЬс (X) dx,

(8.12)

где

—Ь/2

 

 

 

Н ~ ~ Р dP" +

d*f (с) \

(8.13)

de2 Jс=с0(эс; с, Е ) "

 

Произвольную вариацию бс (х) можно всегда разложить в ряд по соб­

ственным функциям фп (х) оператора Н, так как собственные функции эрмитового оператора образуют полную систему ортонормированных функций:

ОО

 

6с(я) =

2

 

(8.14)

 

 

п = 0

 

 

где п — номер собственного значения,

Ап — коэффициент

разложения.

Подставляя (8.14) в (8.12) и используя

условие ортонормированности

Ь/2

 

 

 

 

S

К ( * )

(х ) d x = К т >

(8 .1 5 )

а также определение собственного

значения оператора Н

 

получим:

Н Ф г . = 8 пФп .

(8 .1 6 )

 

 

 

 

 

 

О О

 

 

w

= y 2 '8п М „ І 2-

(8 .1 7 )

 

п = 0

Так как уравнение (8.16) по форме совпадает с одномерным уравнением Шредингера для частицы во внешнем поле, то из общих теорем квантовой механики следует, что п = 0, 1, 2, . . . , оо и еп есть возрастающая функция аргумента п. Коэффициенты Ап в разложении (8.14) не могут быть произ­ вольными. Они должны удовлетворять соотношению

Ь /2

оо

 

j

6c(*)d* = 2 Лпв„ = 0,

(8.18)

—Ь/2

п=0

 

93

где

L/2

ап — 5

■'l'n (*)dx'

(8.19)

- m

 

 

следующему из условия сохранения

числа атомов в кристалле

(7.71).

Для того чтобы установить знак квадратичной формы (8.17), необходимо найти минимум б2F при дополнительном условии (8.18). Если min (б*F) > О, то квадратичная форма (8.17) всегда положительна и, следовательно, перио­ дическое распределение с (я) = со (х, с; Е) является метастабильным. Ми­ нимум квадратичной формы (8.17) при дополнительном условии (8.18) можно найти методом неопределенных множителей Лагранжа. Для этого необхо­ димо минимизировать по Ап выражение

г = ^ 2 8п4 - ^ У п.

(8.20)

В результате получим:

п=о

п=0

 

 

 

 

 

 

 

Ап - S е„

 

(8.21)

 

 

 

Подставляя (8.21) в (8.17), перепишем последнее

в виде

 

 

L/2

ііі

2

 

 

I

А'*

00

пѴ

]'

Ф„ (^х) dx

-U2

(8.22)

 

 

п=0

 

71=0

 

 

 

 

Таким образом, вопрос о

типе

экстремума,

отвечающего функции

со (х\ с, Е), сводится к вопросу о знаке суммы в (8.22).

Интересно отметить,

что сумма в (8.22) может быть выражена через функцию Грина, удовлет­

воряющую операторному

уравнению

 

 

 

HG = 1,

(8.23)

где

Н — оператор (8.13),

6 — функция

Грина в операторном представле­

нии,

1 — единичный оператор, матричные элементы которого есть дельта­

функция Дирака. Билинейное разложение функции Грина уравнения (8.23)

всегда может быть представлено в

виде

 

 

,

,8.24,

71— 0

П

 

где G (х, х') — матричные элементы оператора G в ^-представлении. Инте­ ресующую нас сумму в (8.22) можно выразить через функцию Грина (8.24) уравнения (8.23). Для этого необходимо проинтегрировать уравнение (8.24) по X и у?. В результате получим:

L /2

оо

I 5 % (х) dx I2

 

? = ^

dx'dx'G (ж, х') = 2

(8.25)

—Ь /2

Из предыдущего следует, что при g > 0 возникает метастабильное сос­ тояние системы, при g < 0 — нестабильное состояние, отвечающее седловой точке на гиперповерхности ДF = ДF ({ с (ж)}). Ниже мы рассмотрим кон­ кретный пример, когда вычисление величины g может быть проведено в явном виде.

94

Примем для свободной энергии системы выражение (7.23). При этом положим <*і = а2 = а. Тогда «потенциал», фигурирующий в «гамильтониа­ не» (8.13), будет иметь вид:

U= (4 ^)с= с0(х) = “ — “ (См — с°і)6 (Со (*) — с* )•

(8.26)

Периодическая функция с = с0 (ж) изображена на рис. 24, а. Она при­ нимает значения с0 (ж) = с* в точках, имеющих координаты

X™ = та, ж£> = та + А

(8.27)

(см. рис. 24, а), где а — период распределения концентрации, т — любое

Рис. 24. а) Периодическое

распределение концентрации, б) Зависимость

 

от координат эффективного «потенциала» U ) (8.26).

целое

число,

А — расстояние между двумя соседними точками,

в которых

с (ж) =

с*. Учитывая в (8.26) определения (8.27), можно, без потери общ­

ности,

представить потенциал (8.26) в виде:

 

 

 

д.Щ (с)

 

 

 

 

 

 

de*

)с=с0(*)

 

 

 

 

 

 

 

dcp (ж) \

О О

 

 

= а — а (cos — соі)

2

[б (ж — та) — 6 (ж — та — Л)].

dx

/Сд(х)=с*

 

 

 

 

 

 

 

 

 

 

 

 

 

(8.28)

На самом деле истинные пределы суммирования в (8.28) есть ^

Ыа, однако

в рассматриваемом нами асимптотическом предельном случае Ы а —*оо они

могут

быть

заменены

пределами

+

оо.

 

 

 

 

 

Из

(8.3)

следует,

что

 

 

 

 

 

 

 

 

 

/йсо(ж)\

 

I de \

 

р

 

1

--------------------------------

 

i " * r - J c o(x)=С* =

І1^)с=с* = Т

=

Т

Vm-E-pc'+Af (С-)] .

(8.29)

Подставляя

в (8.29)

выражение

(7.23) при

«і =

а2 = а,

получим:

 

( ^

г ) Со(*)=с* =

K l“

^ ~

ЦС*+4"

а (с*-

е> (с* +

с' - 2®«)] •

(8.30)

Перепишем выражение

(8.28)

для

«потенциала» U:

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

 

 

U (ж) =

а — Ъ 2

[Ö (ж— та) — б (ж — та — А)],

(8.31)

т = —оо

95

Где Константа Ь в

соответствии с (8.30)

равна

Ь = а(сог — coi)

£ — Цс*+ - | - а

(с*— с)(с*+с — 2coi)j| . (8.32)

Принимая во внимание выражение (8.31), можно представить «гамильто­ ниан» (8.13) в следующем виде:

где

Я = Я 0 +

V,

(8.33)

 

 

 

,

 

(8.34)

Яо==—

+ “•

 

ОО

 

 

 

(У )**■ = - & 2 [б (Xта) — Ö(г — та — А)] б (ж — х')1.

(8.35)

т = —оо

 

 

 

Интересно обратить внимание, что по форме он почти совпадает с гамильто­

нианом в модели Кронига — Пенни [47].

Умножая операторное уравнение

(Но'+ V) G =

1

(8.36)

для функции Грина G на невозмущенную функцию Грина G0, удовлетво-

ряющую равенству

 

 

GoH = 1,

 

(8.37)

получим:

 

 

G = Go — GoV G.

(8.38)

Уравнение (8.38) обычно называют уравнением Дайсона. В ^-представлении

уравнение Дайсона

(8.38)

имеет вид

 

 

 

 

 

 

 

 

 

 

 

 

L/2

 

 

 

 

 

 

 

 

G (х, х') =

Go (X — х')

^ dx" dxmGo (х — х") V ( х \

хт) G (хт, х').

(8.39)

 

 

 

 

 

- L / 2

 

 

 

 

 

 

 

 

Подставляя

в

(8.39)

выражение

для

«возмущения»

(8.35),

получим:

 

 

 

 

 

О О

 

 

 

 

 

 

 

 

G (а?, а;') = Со (ос— а:') +

Ь

^

— т а ) ^ (т а *х )

 

 

 

 

 

 

т

е — О О

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

 

 

 

 

Ь

2

 

Со та — А) С (та -}- Л, х').

(8.40)

 

 

 

 

 

т = —оо

 

 

 

 

 

 

Введем

следующие определения:

 

 

 

 

 

 

 

 

 

 

 

 

L/2

 

 

 

 

 

 

 

 

 

 

 

 

go =

\

Go(x) dx;

 

 

 

(8.41)

 

 

 

 

oo

 

-L/2

 

oo

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

got =

 

2

Go(ma),

got =

2

Go (та + A);

 

(8.42)

 

 

 

7TI=S—OO

 

 

 

—OO

 

 

 

L/2

 

oo

 

 

 

 

 

L/2

 

OO

 

 

 

■'*'

3

G (та, x'),

g i =

^

dx'

2

G (ma +

a:7).

(8.43)

—L/2

 

m=—1d o

 

 

 

 

—L/2

 

m*—OO

 

 

96

Интегрируя уравнение (8.40) по х и x' в интервале [— Л/2, L/2] и учи­

тывая определения (8.25),

(8.41)

— (8.43),

получим:

 

 

 

 

 

 

 

g =

Lgo +

bgoëi — bgoëf

 

 

(8.44)

Полагая' в уравнении Дайсона (8.40)

х =

па

и х = па + h, где п

произвольное

целое число,

перепишем его

в виде:

 

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

G (па, х') = Go (па — х') + Ь

^

G° (ne — та) G (,иа- х') ~

 

 

 

 

 

 

 

тп——оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ОО

 

 

 

 

 

 

 

 

 

 

 

 

 

— &

2

 

(ла ”

та — Л) G (та + Л, а;');

(8.45а)

 

 

 

 

 

 

т = —оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

оо

 

 

 

 

 

 

С (/!« + Л,

ж ') = Go (па +

h x') +

b

^

G0(na +

h — та) G (та, x')

 

 

 

 

 

 

 

 

 

т = —оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

00

 

 

 

 

 

 

 

 

 

 

 

 

 

Ь

2

Go (па та) G (та + h, x').

(8.456)

 

 

 

 

 

 

 

 

 

т = —оо

 

 

 

 

 

 

Произведем

интегрирование по х' и суммирование от— L/2a до + Ы2а

уравнений

(8.45а) и (8.456). При

этом получим два уравнения:

 

 

 

L

go + bgoigi — bgmgt, gi =

L

go + bgoigi — bgoigi,

(8.46)

gi = —

где g0, g0i,

go2i gi n Si

определены в (8.41)

— (8.43). Система уравнений (8.46)

имеет решения:

 

 

 

 

 

 

 

 

 

 

 

 

 

gi-

L

'

1 +

5

(goi gm)

 

 

 

L

 

1 b (goi gm)

(8.47)

a

* - h i

/„2

_

„2

 

 

 

 

 

 

 

 

80

 

 

a

^ L - b H g h - g l»

 

 

 

 

1 -

62 (*М -

gm)

 

 

Подставляя

(8.47)

в (8.44),

получим:

 

 

 

 

 

 

 

 

 

 

 

 

26*

 

goi gm

 

 

 

 

 

g = Lgo 1 + —

i-b*(/ü-e3L

]

(8.48)

 

 

 

 

 

 

 

 

 

 

 

 

501

8<н)

 

 

Из определения

невозмущенной

функции

Грина

 

 

 

 

 

 

[ — Р Hfl

+ “ ] Go (ж — x') = 6 — x')

 

 

следует, что

 

 

 

 

оо

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Go (к)

=

^

Go (х) е~,кх dx

= ^

^

а .

 

(8.49)

 

 

 

 

 

—ОО

 

 

 

 

 

 

 

 

Полагая

в (8.49)

к — 0,

получим:

 

 

 

 

 

 

 

 

 

 

 

 

 

go =

Go (0) =

1/а.

 

 

 

(8.50)

Воспользовавшись выражением (8.49) для фурье-образа невозмущенной функции Грина G0 (х), представим выражения (8.42) в виде ряда по «векторам

4 А. Г. Хачатурян

97

обратной

решетки»

kp — (2n/a) p

(p = 0, ± 1 ,

±

2,

. . . ,

± oo) периодиче-

CKoro

распределения

концентрации с = с „ ( ж ):

 

 

 

 

g01

а

2

n 1

\a

= T h abW j r ° '

 

 

 

P=—°° ß I

a

P J + a

 

 

 

 

 

 

 

 

 

 

 

 

(8.51)

 

 

 

 

 

2л:

 

 

 

 

 

 

 

1

 

 

cos — ph

 

i _

 

 

 

261

 

 

 

 

a r

 

 

 

 

g02 = —

2

 

 

2

 

2 /aß

1

, / 1 Г

 

 

 

 

 

! +ot

 

a

 

 

P = — - C O

 

 

 

sh

2

| /

ß

 

 

 

 

 

 

 

 

 

 

 

Подставляя (8.51) в (8.44), получим окончательное явное выражение для величины g:

 

——

ch

 

cth т / t *

 

? = — 1

Ь2

- Г ^ Т *

з_ і_

 

2 о2

aa2 ß 62 - , / t j {a~ 2h)

4?ß cth2 -

V j a ~

(8.52)

Из выражения (8.52) следует, что в зависимости от периода а исследуемого распределения концентрации с = с0 (х) величина g может принимать как положительные, так и отрицательные значения. Последнее показывает, что периодические распределения концентрации в некоторых случаях могут отвечать метастабильному состоянию системы.

Г Л А В А ІП

СТАТИСТИЧЕСКАЯ ТЕОРИЯ ФАЗОВЫХ ПЕРЕХОДОВ ПОРЯДОК - БЕСПОРЯДОК В ДВУХКОМПОНЕНТНЫХ РАСТВОРАХ ЗАМЕЩЕНИЯ И ВНЕДРЕНИЯ

§ 9. Модель твердого раствора в статистической теории упорядочивающихся сплавов

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

Чаще всего используется модель, в основе которой лежит предположение о том, что атомы сплава размещаются по узлам некоторой жесткой кристаллической решетки. При этом конфи­ гурационная энергия сплава представляется в виде суммы всех парных потенциалов межатомного взаимодействия. Предположе­ ние о парном характере межатомного взаимодействия может быть сравнительно строго обосновано для непереходных металлов и металлических сплавов. В последнее время появились работы, в которых методами теории псевдопотенциалов было показано, что полная энергия электрон-ионного, ион-ионного и электрон-элек- тронного взаимодействия в непереходных металлах и металличе­ ских сплавах непереходных элементов может быть довольно точно представлена в виде суммы всех парных межатомных взаимодей­ ствий [48, 49].

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

4* 99

Соседние файлы в папке книги из ГПНТБ