книги из ГПНТБ / Хачатурян, А. Г. Теория фазовых превращений и структура твердых растворов
.pdfподобно тому как они были поЯученМ выше, для случая ОЦК ре шетки Изинга.
Параметры дальнего порядка в распределениях (13.36) —
(13.38) описываются общим трансцендентным уравнением |
|
|||
in |
_ |
П к ) |
(13.42) |
|
( I - c+ t a H^ + VU |
яг Ts4s’ |
|||
|
||||
где s — 1,2,3, |
|
|
|
|
kj — 2л (a! -f а2 + а3), |
л (а* + а*), к3 — л (аі “Ь э2 -f- а3)I |
|||
|
|
|
(13.43) |
|
базисные векторы обратной решетки а*, |
а*, а* определены в § 11, |
|||
а V (ks) имеет вид (13.19). При обычной |
нормировке параметров |
|||
дальнего порядка (т]і = Цг = |
т]з — 1 в полностью упорядоченном |
|||
состоянии) Ys = 1/2. |
|
|
|
Уравнение (13.42), как и аналогичное уравнение (13.17), справедливо при наличии взаимодействия в произвольном числе координационных сфер.
Подставляя (13.43) в (13.19), получим, что для случая ОЦК решетки
V (kj) = |
2 |
V (R) е ~ і2л(аі +а2+ аз )К = |
_ |
8 u 7X + Öw2 + |
1 2 w 3 - |
Эи?* + . . . , |
|
R |
|
|
|
|
(13.44) |
|
|
|
|
|
|
|
V (k2) = |
2 |
V (R) er (V a3)R = _ |
2w2- 4u>3 + . . . , |
|
(13.45) |
|
|
R |
|
|
|
|
|
V (k3) = |
2 |
Г (R) е-іл(W a3)R = |
- |
6w2 + 12m;3 + |
. . . , |
(13.46) |
|
R |
|
|
|
|
|
где wi, w2, w3, ... — энергии смешения соответственно в первой, второй и т. д. координационных сферах.
Уравнение (13.42) при s — 1 описывает температурную и кон центрационную зависимость параметра tjj. в распределении (13.36), при s = 2 — параметра rj2 в распределении (13.37), при s = 3 — параметра ц3 в распределении (13.38).
Температурная и концентрационная зависимость параметров дальнего порядка и т)2 в распределении (13.39) описывается си стемой двух трансцендентных уравнений (13.31), в которых, одна ко, сверхструктурные векторы Ь* и к2 определены в (13.43). Зави симость параметров дальнего порядка г]і и г)2 в распределении (13.40) описывается системой двух уравнений с двумя неизвест ными:
1 (1 — с — гцті — бірта) (с |
+ |
ЛіТі |
- 2т]2Т2) __ |
87 (ki) |
|
|
||
(1 — с — TjiY i + |
2т]2Та(с) |
+ |
т ц Т і |
-4- б т р т з ) |
яТ |
І2 ‘ 2 ’ |
(13.47) |
|
1п (1 - с — т)іТі + |
2г|2Т2)(с - |
гцтх) |
_ 2Ѵ (ki) |
|
2V (k2) |
|||
11'1 |
І2‘2’ |
|||||||
(1 — с + rjiYi) (с + тцті — 2т]2Т2) |
v T |
я Т |
140
Обычной нормировке параметров дальнего порядка отвечает выбор коэффициентов Ѵі = Уг = 1/8. Векторы кх и к2 в (13.47) опреде лены в (13.43).
Наконец, зависимость параметров ■%и г|3 в (13.41) может быть найдена из системы уравнений:
1 |
(1 — с + |
тцті — т)зГз) (С - г)іТі — РзТз) |
_ 2Ѵ (k3) |
|
|
|
|
|||
1 |
(1 — с + |
rjiifi + |
г]зТз) (с — ЛіТі + |
W s) |
хЗ1 |
|
‘3 |3’ |
|
(13.48) |
|
|
|
|
|
|
|
|
|
|
||
. ■ |
(1 —с + |
тцті + |
груз) (с + rim) |
_ |
2.V (kt) |
|
|
V (k8) |
“ПвТа- |
|
1 |
(1 — с — тцті) (с — rjiTi — г}зТз) |
|
кТ |
11 |
1 |
хТ |
|
|||
|
|
|
||||||||
Нормировка |
параметров дальнего порядка, |
|
при |
которой |
= |
= t]8 -- 1, в полностью упорядоченном состоянии достигается при
Yi = — V*, Уз = Ѵ2. |
(13.49) |
Векторы кх и к3 в (13.48) определены в (13.43).
Распределение вероятностей (13.36) описывает сверхструкту-
ры типа В2: CuZn, GuBe, FeAl, CuPd, AuZn и т. д. Изоморфная
им фаза внедрения имеет структурную формулу Ме2Х. Распре деление (13.37) описывает сверхструктуру замещения, изобра женную на рис. 28,А ѵ Изоморфная ей фаза внедрения Ме2Х
141
изображена на рис. 28, В г. Соответствующая фаза внедрения Та20 была обнаружена в [6] (см. § 11).
Распределение (13.38) описывает упорядоченные фазы заме щения NaTl, LiAl и т. д. (рис. 28, И3). Изоморфная фаза внедре ния Ме2Х изображена на рис. 28, В3. Распределение (13.39) описывает сверхструктуру замещения A SB, изображенную на рис. 28,Лгу и изоморфную сверхструктуру МеАХ, изображенную на рис. 28, В 2. Фаза Та40 с такой структурой была описана в § 11. Распределение (13.40) относится к сверхструктуре замещения А 2В, в которой атомы сорта В образуют ОЦК решетку с удвоенным, по сравнению с решеткой Изинга, параметром. В изоморфной ей фазе внедрения МевХ атомы внедрения, располагаясь в октаэдри ческих междоузлиях одной подрешетки, также образуют ОЦК решетку удвоенного периода. Распределение (13:41) описывает сверхструктуру типа DOs (Fe3Al, Fe3Si), которая изображена на рис. 28, Ац, и изоморфную ей сверхструктуру внедрения Ме^Х
(рис. 28, В 4).
На рис. 28, как и на рис. 27, под каждой парой изоморфных сверхструктур замещения и внедрения изображена соответствую щая обратная решетка. Она одинакова для изоморфных сверхструк тур. Расположение сверхструктурных узлов обратной решетки определяется векторами kjs, входящими в соответствующие распре
деления вероятностей п. (R). Сравнивая наблюдаемые дифракци онные картины с расположениями структурных и сверхструктур ных узлов обратной решетки, приведенными на рис. 27, С и рис. 28, С, можно идентифицировать структуру исследуемых упорядоченных фаз.
§ 14. Метод статических концентрационных волн в уравнениях самосогласованного поля (сложные решетки Изинга)
Как уже упоминалось, сложные решетки Изинга представля ют собой несколько взаимно проникающих решеток Бравэ, сдви нутых относительно друг друга на расстояния hp (р — 1, 2,. . ., ѵ, где V — число решеток Бравэ). Иными словами, сложную решетку Изинга можно представить себе как решетку, в которой на каж дую элементарную ячейку Бравэ приходится ѵ узлов, сдвинутых относительно центра элементарной ячейки на те же векторы hlt h2, . . ., hv соответственно. Эти узлы образуют базис или мотив сложной решетки Изинга. Трансляция узлов базиса дает упомя нутые выше V подрешеток (взаимно проникающих простых ре шеток Бравэ). Всюду в дальнейшем мы будем рассматривать ча стный случай, когда все узлы решетки Изинга кристаллографи чески эквивалентны, т. е. могут быть совмещены друг с другом одним из преобразований симметрии неупорядоченного крис талла.
142
Таким образом, положение каждого узла г сложной решетки Изинга может быть задано двумя векторами (R, 1ір). Вектор R определяет положение центра элементарной ячейки, в которой находится данный узел, вектор h;, — положение данного узла от носительно центра элементарной ячейки. При этом
r = R + hp. |
(14.1) |
Так как парный потенциал взаимодействия двух атомов, на ходящихся в узлах г и г', не может измениться при смещении начала координат на вектор трансляции решетки Изинга, то его можно представить в видег
V (г, r') = Vpq (R - R'), |
где r' = R' |
+ h„ q = 1, 2, . . . , v. |
(14.2) |
Используя выражение (10.3) в (10.2), перепишем уравнение |
|||
самосогласованного поля в виде |
|
|
|
п (г) = [exp ( |
^ ft Н “jr 2 |
V (г, г') п (г')) + і ]_1 . |
(14.3) |
Оно, в частности, может быть получено, если приравнять нулю первую вариацию по п (г) от свободной энергии
*■ = 4 - 2 |
о * о-) * (О + |
|
|
г,г' |
|
|
|
+ Y.T 2 {п (г) In п (г) + [1 — п (г)] In [1 |
п (г)]} — ц 2 п (г)- |
(14-4) |
|
г |
|
г |
|
В (14.4) первое слагаемое |
|
|
|
|
= 4 - 2 F ( r , r ' ) n ( r ) n ( r ' ) |
(14.5) |
Г,І*'
есть выражение для внутренней энергии.
Используя представление (14.2) в уравнении (14.3), перепишем уравнение самосогласованного поля в форме
|
V |
|
п(р, R) = {ехр Г— JJr + ^ |
2 2 |
- К>(<7. R')l + l}"1 • |
|
|
(14.6) |
Уравнение (14.6) представляет собой нелинейное конечно разностное уравнение относительно п (р, R) — вероятности найти атом данного сорта в узле (р, R). Поэтому зависимость от коорди нат (р, R) решения уравнения (14.6), как и раньше (см. § 10), оп ределяет симметрию упорядоченной фазы. Кроме решения
п (р, R) = с,
где с — атомная доля атомов данного сорта, отнесенная к полному
143
числу узлов решетки Изинга (для растворов замещения с есть атомная доля данного компонента, для растворов внедрения с — отношение числа атомов внедрения к числу мест внедрения), урав нение (14.6), как правило, обладает другими решениями п (р , R), обнаруживающими зависимость от координат узлов (р , R). Каждое такое решение описывает свою сверхструктуру замещения или внедрения.
Уравнение (14.6) можно существенно упростить, если пред ставить функцию распределения атомов п {р, R) в виде линейной суперпозиции статических концентрационных волн, имеющих вид функций Блоха:
фок(Р, R) = Ѵак(р) ехр (— ikR), |
(14.7) |
где гок(р) представляет собой «вектор поляризации» волны с вол новым вектором к, находящимся в первой зоне Бриллюэна решет ки Изинга; а — номер «поляризации». Аналогично тому, как плоские волны ехр (—ikR) являются собственными функциями матрицы V (R — R'), образуемой потенциалами межатомного взаимодействия в простой решетке Изинга, волны (14.7) являются собственными функциями матрицы Vpq (R — R') потенциалов вза имодействия в сложной решетке Изинга:
2 |
2 Vvq (R - R') Ф«к (?, R') = К (к) фак ( Р , R), |
(14.8) |
9= 1 |
R ' |
|
где Ха (к) — спектр собственных значений матрицы Fpg(R — R'), индекс «поляризации» о выступает здесь в роли номера ветви спектра собственных значений Ä,0(k). Так как матрица Vpq(R —R') эрмитова, то ее собственные значения Я0(к) — действитель ные числа. Подставляя в уравнение (14.8) представление (14.7), получим, что «векторы поляризации» г0к (р) удовлетворяют урав нению для собственных значений:
2 Ѵ РЯ (к) |
(?) = ко (к) Увк (р), |
(14.9) |
Vpq (к) = 2 |
Vpq (R) ехр (— ikR). |
(14.10) |
R |
|
|
Из (14.9) следует, что каждому волновому вектору к отвечают ѵ ветвей XQ(к), а именно (к),Х2(к), . . ., Х„(к), и ѵ собственных «векторов» г0к (р), т. е. vl k (p), г2к (р), ■■., ѵѵк(р). Так как матрица Vpq (к), по определению, эрмитова, то все \ г ее собствен
144
ных «векторов» rok (р) ортогональны друг к другу:
V |
|
|
2 |
(Р) yo'k(Р) = бса-• |
(14.11) |
Р=1 |
|
|
Так как статические концентрационные волны фок (р, R), будучи собственными функциями эрмитовой матрицы VPQ(R — R'), образуют полную систему ортонормированных функций, то вероятности п (р , R) можно представить в виде линейной супер позиции этих концентрационных волн:
V |
|
|
п (р, R) = С + -і- 2 2 |
K?« (k) (р)e_ikR + |
Qa (k) ѵ'л (p) etkR], |
0=1 k |
|
(14.12) |
|
|
|
где Qa (k) — амплитуды |
концентрационных |
волн rak (p) e~ikR. |
Представление решения уравнения (14.6) в виде суперпозиции концентрационных волн (14.12) осуществляет переход от описа
ния распределения атомов с помощью |
vN вероятностей п (р, R) |
к описанию с помощью vN амплитуд |
(?e(k) (первая зона Брил |
люэна содержит N волновых векторов к, причем каждому векто |
|
ру к отвечают ѵ «поляризаций» ѵвк(р)). |
Этот переход аналогичен |
соответствующему переходу от смещений к амплитудам нормаль ных колебаний в теории колебаний сложных решеток.
Из симметрии решетки Изинга следует, что собственные значе ния %а(к) вырождены относительно всех собственных функций фок (р> R), относящихся к одной ветви а, но к разным векторам звезды вектора к. Поэтому разложение (14.12) можно перегруп
пировать таким образом, чтобы собрать вместе |
все концентра |
|||||||
ционные волны Vgk(p) exp (—ikR), |
относящиеся к одному и тому |
|||||||
же собственному значению Я0 (к). |
При этом получим: |
|
|
|||||
|
П {р, R) = |
С+ |
2 |
Ли Вза (р, R), |
|
(14.13) |
||
|
|
|
о , S |
|
|
|
|
|
где |
, . \ |
, . |
-ikj |
R , |
» .. . • ikj R. |
|
|
|
1 чгд г |
(14.14) |
|||||||
eie(p, R) = -2- 2 і[Т«(7.)У«кУі(р)в |
|
+ Та, (/,) Ѵак]е |
’• ], |
|||||
|
Q° (kj,) = ЛавТа« (Іг)- |
|
(14.15) |
|||||
Как и в § 10, индекс s нумерует звезды волновых векторов, ин |
||||||||
декс /, — волновые векторы в пределах звезды s. |
Суммирование |
|||||||
в (14.14) производится по всем векторам звезды s. Величины |
|
ті05 — |
||||||
параметры дальнего порядка. Из (14.15) |
следует, |
что параметры |
||||||
дальнего порядка определены неоднозначно с точностью до |
мно |
|||||||
жителей у оз (/в)- |
Д л я однозначного определения |
требуется |
еще |
145
одно условие. Им может быть либо требование нормировки
2 I Tos (/'s) |2— |
(14.16) |
либо же обычное условие нормировки, чтобы в полностью упоря доченном состоянии все параметры дальнего порядка сверхструк туры были бы равны единице.
*“ik*R
По определению все волны і;ак. (р) е 3s , входящие в функ- h
цию egs(p, R) (14.14), являются собственными функциями матри цы Vvq (R — R') и относятся к одному и тому же собственному значению %а(ks) — ^osПоэтому функция е„3(р , R) также явля ется собственной функцией матрицы FP9(R — R'), отвечающей собственному значению Âos:
V
2 2 Ѵѵя (R - R') e.. (P, R') = K#a. (p, R). |
(14.17) |
q=l R' |
|
Подставляя (14.13) в (14.6) и используя условие (14.17), получим уравнение самосогласованного поля:
С+ 2 Лавбв. (р, R ) =
0,8 |
|
= [ехр (— |
+ ~ w 2 tasW a, (Р, R ) j + l ]-1 . (14.18) |
Уравнение (14.18) является аналогом уравнения (10.15), спра ведливого для простой решетки Изинга. Уравнение (14.18), точно так же, как и (10.15), может быть сведено к системе трансцендент ных уравнений относительно параметров дальнего порядка т]03. Для этого необходимо придать координатам (р, R) узлов сложной решетки Изинга конкретные численные значения.
Обобщая рассуждения, приведенные в Приложении 2, на слу чай сложной решетки, можно показать, что коэффициенты ygs(js) в выражении (14.13) (как и коэффициенты ys(js) в выражении (10.9)) являются константами. Эти константы не зависят от темпе ратуры, состава и давления в однофазных областях диаграммы рав новесия. Воспользовавшись этим обстоятельством, можно пока зать, что функция п (р, R), имеющая вид (14.13), может быть ре шением уравнения (14.6), если она удовлетворяет тому же усло вию I, которое было сформулировано в § 10:
число значений, принимаемых функцией п (р, R) на множестве всех узлов сложной решетки Изинга, должно бытъ на единицу больше, чем число параметров дальнего порядка ц03.
Условие I позволяет найти численные значения коэффициентов уas (js) и, следовательно, определить атомное строение сверх структуры.
146
Температура фазового перехода порядок — беспорядок, как обычно, определяется точкой ветвления уравнения самосогласо ванного поля (14.6). Температура ветвления может быть получена путем линеаризации уравнения (14.6). Заменяя в (14.6) вероят ность п (р , R) на с + Ьп(р, R), где 6п(р, R) — малая вариация, и разлагая (14.6) относительно 6п(р, R), получим линеаризован ное уравнение:
Öп (Р, R) = - -с--^ с)- 2 Ѵрч (R - R') К (R')- |
(14.19) |
< 7 , R ' |
|
Умножая правую и левую части уравнения (14.19) на ^ok(p)exp(— ikR) и производя суммирование по (р, R), получим:
öQ* (к) = - |
К (к) бQa(к), |
(14.20) |
где
б(?а (к) = 2 Ö« (р, R) Уак (р) e~ikR.
T>,R
Уравнение (14.20) имеет нетривиальное решение 6@о ( к ) ^ 0
(оно отвечает случаю 8 n( p, R )^0 ), е с л и ----с^ ~ с^ |
А,д(к) = 1, |
||
или |
g (1 — <0 К (к) |
|
|
Ток = — |
(14.21) |
||
X |
|||
|
|
Уравнение (14.21) описывает континуум точек ветвления урав нения (14.6). Каждой такой температуре Таk отвечает потеря устойчивости неупорядоченного состояния относительно концен трационной волны с соответствующим волновым вектором к и «поляризацией» а. Температура фазового перехода порядок — беспорядок определяется максимальной температурой Гок:
Т0 = |
max |
c(l-c)XB(k) |
С (1 — с) |
min ко (к). (14.22) |
|
X |
|||||
|
|
X |
|
Структура упорядоченной фазы определяется поляризацией о = = о0 и звездой (kja), обеспечивающими минимум функции Я0(к):
(к,-§) = min ка(к).
В качестве примера приложения теории упорядочения в слож ных решетках Изинга, изложенной в настоящем параграфе, рас смотрим конкретный случай сверхструктуры типа Та20. При по лучении сверхструктуры Та20, изображенной на рис. 27, Blt мы полагали, что атомы кислорода располагаются только в одной (третьей) подрещетке октаэдрических междоузлий, и, следова тельно, пренебрегали эффектом перераспределения атомов внед рения между тремя подрешетками октаэдрических междоузлий. Для того чтобы учесть последнее обстоятельство (оно особенно су щественно вблизи температур фазового перехода порядок — бес
147
порядок), необходимо прибегнуть к описанию вероятностей рас пределения атомов кислорода с помощью выражения (14.13).
Сверхструктура типа Та20 связана с двумя векторами
к0= 0, кх = я (а^ + а3) |
(14.23) |
(см. § 11). Поэтому построение функции п (р , R) в виде (14.13) необходимо начать с определения собственных значений Ка(к) и собственных функций ѵ0]і (р) уравнения (14.9), отвечающих век торам (14.23). Вектору к0 = 0 отвечает матрица FPg(0), имеющая вид
(Ѵ пф ) |
Ѵѵш(0) |
Ѵи (0)\ |
где Vpq(0) = 2 v vq (R). (14.24) |
VPq(0) = va (0) |
Vn (0) |
vn (0) , |
|
\F i2(0) |
F12(0) |
Vu (0)/ |
R |
Свойства симметрии матрицы (14.24) связаны с кристаллографи ческой эквивалентностью трех подрешеток октаэдрических междо
узлий в ОЦК решетке. |
|
|
следующие |
собственные значения |
||||||
Матрице |
(14.24) отвечают |
|||||||||
%а(0) и соответствующие им собственные векторы |
0(р): |
|
||||||||
КФ)= Угг (0) + 2Ѵ12(0), |
|
КФ)= |
Я. (0) = |
Fu (0) - |
F12(0); |
|||||
|
|
|
|
|
|
|
|
|
(14.25) |
|
|
|
|
|
|
|
|
|
|
(14.26) |
|
Элементы столбцов в собственных векторах |
г?о>0(р) |
в (14.26) — |
||||||||
компоненты |
собственного |
вектора |
(vO)0(l), r0,o(2), нв)0(3)). |
Все |
||||||
собственные векторы г?о>0 (р) |
нормированы на единицу. |
|
||||||||
Вектору |
кх = я (а 2 + |
а3) |
отвечает матрица FPQ(n(aa + |
а3)), |
||||||
имеющая вид |
|
|
|
(Ѵ к |
0 |
0 \ |
|
|
||
|
|
|
|
|
|
|||||
|
Fpg (Я (а£ + |
аз)) = |
0 |
F22 |
Fas |
, |
(14.27) |
|||
где |
|
|
|
\ |
0 |
Fss |
Faa/ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Fpg = |
2 |
Vpg (R) еія(Ѵ аз)К. |
|
(14.28) |
|||||
|
|
R |
|
|
|
|
|
|
|
Легко показать, что общий вид матрицы (14.27) определяется только симметрией решетки Изинга и не связан с конкретным вы бором потенциалов Vpq (R — R'). Матрица (14.27) имеет следую щие собственные значения К (кх) и собственные векторы н„кі (р) '■
К(К)= У22 + Fa,, К(К)= F22 - |
F23, КІЮ = Fu ; (14.29) |
|
= |
^ a 0 >)= |
Ѵз^ р) = (о )' (14,30) |
148
Воспользовавшись собственными векторами (14.26), (14.30), можно записать выражение (14.13) в форме
„ (р, R) - (»(2; R)) = с ( і ) + % Л . - ± Г Q Q +
+ [лнТи |
( О + Л»Т» |
f l ) |
+ ЛиТп М ] |
«'■'*•***• |
||
Обозначая |
|
|
|
|
|
(14.31) |
|
|
|
|
|
|
|
^ з о Т г о - |
= Л ь Л зоТ зо у |
— = |
Ла» |
Л п Т и |
у ^ - = |
Лз» |
|
|
|
|
|
|
(14.32) |
|
Л 2 1 Т 2 1-у=г = |
1*>т |
’ Т и Т и = |
% |
|
и используя определение (11.2) радиуса-вектора R, получим:
/«(1, R)\ п(р, R) = I n(2, R) 1=
- c ( I ) + % ( ! ) + % ( ! ) + [ % 0 + ' 0 + ' Q J r * ” -
(14.33)
Легко убедиться в том, что функция (14.33) удовлетворяет ус ловию I : она принимает шесть значений на множестве всех узлов (р, R) решетки Изинга и зависит от пяти параметров дальнего по рядка. Из (14.33) следует, что упорядоченное распределение ато мов в каждой из трех решеток Бравэ, образующих решетку Изинга, описывается одним и тем же распределением (11.4):
п (1, R) |
= |
с (1) |
+ |
г\ (1) ехр я |
(у + |
z), |
|
п (2, R) |
= |
с (2) + |
т) (2) ехр я |
(у + |
z), |
(14.34) |
|
п (3, R) |
= |
с (3) |
+ |
ц (3) ехр я |
(у + |
z), |
|
где
с (1) = с — Г]! — ц2, с (2) = с — % + Т|2, с (3) = с + 2ци (14.35)
Л (1) = Jls. Л (2) = Лз — 4«. Л (3) = Лз + Л4-
Этот вывод следует из общих соображений симметрии и был сделан ранее в § 10 (см. выражение (10.45)). Температурные зави симости концентраций атомов внедрения с (р) и эффективных па раметров дальнего порядка rj(p) в трех подрешетках октаэдри ческих междоузлий могут быть определены с помощью соотноше ний (14.35), связывающих с{р) и т] (р) с истинными параметрами дальнего порядка. Последние же могут быть найдены с помощью уравнений самосогласованного поля (14.18).
149