где 'У'і — поток і-компонента |
в направлении х через единичную |
|
|
|
dot |
— градиент концентрации |
площадку, перпендикулярную х, —L |
і-компонента в направлении х, |
дх |
■— |
коэффициент диффузии. |
а А |
О о о |
о о |
о о о о |
о |
о о о о о |
о о о |
о о |
о о о о о |
о о о о |
о |
о о о |
о о |
о о |
|
о |
о о |
о |
о о о |
о о |
о о о о а |
о о о о о |
а |
|
|
б |
|
в |
|
Рис. 4.1. Различные механизмы диффузии: а) диффузия по междо узлиям; б) диссоциативная диффузия; в) диффузия по вакансионному -механизму
Поскольку диффузия совершается путем перескока атомов из одних кристаллографических позиций в другие, то коэффициент диффузии можно выразить соотношением
где Г — число перескоков атомов в единицу времени, a s — длина отдельного скачка. Величина s может быть выражена через посто янные решетки, а величина Г является функцией нескольких пере менных. При вакансионном механизме диффузии Г пропорцио нально концентрации вакансий [V] и частоте перескока атомов из регулярных в соседние вакантные узлы (со).
Для кристалла с кубической решеткой
где а0 — постоянная решетки, а а — коэффициент, зависящий от геометрии кристалла. Очевидно, что
/ AGv \ |
ASv |
/ |
АНѵ \ |
(4.4) |
[У] = ехр ( — — |
) - ехр— |
ехр (— — ) , |
где А Gy, А5у и АНу — изменение свободной |
энергии, энтропии |
и энтальпии процесса образования вакансий соответственно. |
|
Согласно теории активных комплексов [8] |
|
|
со = £ѵехр(— ^ L |
^ |
= £Ve x p ^ e x p ( — |
(4.5) |
где k — трансмиссионный коэффициент, характеризующий вероят ность того, что атом с достаточной для скачка энергией действи-
тельно совершит перескок, ѵ — частота, а АGm, ASm и АНт — из менение свободной энергии, энтропии и энтальпии активации про цесса соответственно. Величина k, по-видимому, близка к единице [8], а частота ѵ в первом приближении равна дебаевской час тоте [4, 8].
Комбинируя уравнения (4.3) — (4.5), находим
ASV -р AS„ |
А/К |
АН„ |
D — аа\ kv exp |
exp |
(4.6) |
|
|
RT |
Соотношение (4.6) показывает, |
что интенсивность диффузион |
ного массопереноса через кристалл зависит от легкости образова ния в нем вакансий и становится тем больше, чем выше концент рация последних. Однако выражаемая уравнением (4.3) простей шая взаимосвязь между коэффициентом диффузии и концентра цией вакансий справедлива лишь при условии, что концентрация вакансий сравнительно невелика и -перескоки атомов в вакантные узлы решетки взаимонезависимы.
В противном случае следует учитывать, что коэффициент диффузии пропорционален величине, выражающей вероятность заполнения соседнего с вакансией узла. Поэтому для нестехиомет рического кристалла МХі_ѵ, в котором дефицит компонента X связан с образованием вакансий в Х-подрешетке, коэффициент диффузии атомов X
где Dy — коэффициент диффузии вакансий.
Допустим теперь, что между вакансиями действуют отталки вающие силы, так что при некоторой концентрации у' вакансии полностью упорядочиваются, образуя новое соединение МХі_г . В этом случае
|
Dx = Dvy ( 1 - ^ t ) . |
(4.8) |
При у = у ' -Db — 0> |
т- е- вакансии становятся неподвижными и теря |
ют способность |
участвовать в |
массопереносе |
компонента В1 |
(во всяком случае по вакансионному механизму). |
распределения |
Рассмотрим, |
как эффект |
небеспорядочного |
вакансий влияет на соотношение между коэффициентами диффу зии, измеренными различными методами. Коэффициент диффузии можно определить, измеряя частоту перескоков атомов или ионов (методами внутреннего трения, ЯМР или нейтронной спектроско пии [6, 7]). С другой стороны, используют микроскопические мето ды, основанные на измерении потока материала в концентрацион ном поле. Очевидно, что в последнем случае имеет место взаимная
1 Напомним, что полное упорядочение.вакансий приводит к образованию новой структуры, относительно которой упорядоченные пустые узлы, строго го
воря, |
являются не вакансиями, а междоузлиями. |
|
1 9 1/ 2 |
ю . Д . Т р е т і . я к о з |
289 |
диффузия обоих компонентов, совершаемая с различной скоростью и сопровождаемая смещением первоначальной поверхности между составляющими диффузионной пары (эффект Киркендаля [4, 7]).
Будем рассматривать бинарное соединение, в котором подвиж ность одного компонента намного превышает подвижность второго. Возможны два рода экспериментов по измерению коэффициентов самодиффузии:
а) компонент і диффундирует навстречу незанятым узлам (вакансиям или межузельным позициям), приводя к изменению нестехиометрии кристалла;
б) стабильные атомы компонента і диффундируют навстречу атомам радиоактивного изотопа того же компонента. Получаемые в этих опытах коэффициенты диффузии обозначим соответственно
символами D{ и D f. Чтобы связать значения Di и D f , напом ним, что движущей силой изотермической диффузии является гра диент электрохимического потенциала, так что скорость диффун дирующих частиц
|
ѵі ----- Jі = |
— ut grad Pi, |
(4.9) |
где иI — подвижность атомов і, не зависящая от концентрации. |
Учитывая, что ці = |
ці0+ £ ’7'1паі, из уравнений (4.1) и (4.9) находим |
|
А |
/ _СЧП£І_\ |
(4.10) |
|
V\ д\lnпсCj\ J |
|
|
|
Подвижность |
радиоизотопа |
(и*) в общем |
случае отличается |
от подвижности стабильного изотопа (щ), т. е. |
|
|
|
fu„ |
(4.11) |
где f — некоторый корреляционный множитель.
Учитывая, что смесь радиоактивного и стабильного изотопов
|
одного и того же элемента является идеальной, находим |
|
|
D f = |
ux RT = fu{RT. |
(4.12) |
|
Из сопоставления уравнений (4.10) и (4.12) следует |
|
|
ъ , _ ! _ д х / _ 3 Іпщ_\ |
(4.13) |
|
1 |
/ |
V а)in cj |
|
|
Для нестехиометрического кристалла, являющегося достаточно разбавленным раствором нейтральных вакансий, справедливо со отношение
|
д Іпаі |
и, следовательно, |
|
|
|
д In о |
|
|
|
Yi |
|
(4.14) |
|
|
Dx = /AYi. |
где Yi |
мольная доля вакансий. |
г. |
|
Чтобы оценить, насколько коэффициент / отличается от I,
|
|
|
|
|
|
|
|
вернемся к уравнению |
(4.11). В случае беспорядочно перемещаю |
щихся вакансий, подвижность |
радиоактивного изотопа меньше |
ожидаемой |
(/< 1), так как после обмена |
местами вакансии и ра |
диоизотопа |
возникает |
ситуация, |
в которой |
велика |
вероятность |
того, что следующий |
скачок |
радиоизотопа |
будет |
совершаться |
в исходное положение, т. е. имеет |
место |
известная |
корреляция |
скачков. Доказательства кооперативного (взаимозависимого) дви
жения атомов (вакансий) были получены |
Мапингом |
[9] путем |
сравнения значений коэффициента самодиффузии, |
измеренного |
с использованием радиоизотопа и рассчитанного по |
уравнению |
Нернста—Эйнштейна. |
(4.13) на примере типич |
Рассмотрим применимость уравнения |
ного нестехиометрического окисла, каким является вюстит FeOy. Объемная диффузия кислорода в вюстите пренебрежимо мала [10], тогда как коэффициент диффузии Fe, установленный методом
|
изотопного обмена |
[11—13], для FeOi;09 |
выражается |
уравнением |
|
Dpe ==4,1- ІСМехр |
27 800 |
(4.15) |
|
RT |
|
|
|
|
|
Примечательно, |
что величина Dpe |
значительно |
возрастает |
с увеличением нестехиометрии вюстита, особенно при температу рах выше 900° С [14]. Известные в литературе данные [15—17] об активности Fe в нестехиометрическом вюстите можно выразить уравнением
31ngpe |
= 28,5 у2 |
(при 1050°С). |
(4.16) |
д InС; |
|
|
|
|
Используя значения |
Dpe И |
д 1ПЯре |
|
|
выражаемые уравнениями |
(4.15) и (4.16) соответственно, |
|
д ln |
(значение для |
и полагая / = 0,78 |
диффузии по вакансионному механизму в кристаллах с объемно центрированной решеткой), легко рассчитать по уравнению (4.13) химический коэффициент самодиффузии. При температуре 1050° С
Dpe = 4,8 • ІО-7 см2сек~1. Это более чем на порядок превышает все известные значения химического коэффициента самодиффузии
|
|
|
|
в вюстите. По |
данным различных авторов [19—21] для FeOi,i<> |
при 1050°С Dpe |
составляет от 0,9-ІО-7 до 5-10^7 см2-сек і. Более |
того, величина Dpe, измеренная |
экспериментально, |
уменьшается |
по мере отклонения вюстита от |
стехиометрического |
состава, что |
противоречит уравнениям (4.13) |
и (4.14), полученным для случая |
невзаимодействующих дефектов. Это может означать, что в окис лах с сильно взаимодействующими дефектами (к их числу при надлежит вюстит, имеющий кластерную структуру), диффузия имеет специфический механизм. Так, массоперенос при окислении вюстита требует или движения кластеров, или их разрушения
с последующим переносом индивидуальных атомов — и то и дру гое энергетически затруднено. Диффузию же радиоизотопа можно описать в рамках моделей, не связанных с разрушением кластеров.
Другой эффект, вносящий существенный вклад в диффузию составных частей нестехиометрических окислов, связан с возникно вением напряжений кристаллической решетки, сопутствующих процессу окисления (или восстановления), во всяком случае доста точно быстрому. В таком кристалле химический потенциал под вижного компонента в решетке предложено [18] выражать урав нением
Ві = Вт + RT 1па° 6[х, |
(4.17) |
где а0 — активность і компонента в кристалле, не имеющем напря жений, а бр — изменение химического потенциала, обусловленное
напряжениями. С учетом величины бр уравнение (4.17) |
принимает |
"■ид |
|
|
|
|
|
|
D, - |
д In а? |
J |
d(8pi) |
(4.18) |
|
DC |
RT |
d In с. |
|
|
д In с. |
|
|
1іо оценке О’Кифа [18], значения слагаемых |
и ш а г |
|
---------- и |
|
|
|
|
|
д Inc, |
—1— |
d(fyj)- сравнимы друг с |
другом, |
составляя 20 |
и 30 соот- |
RT |
д In с\ |
|
|
|
|
ветственно. Вместе с тем для окислов с узкой областью гомоген ности вклад напряжений в коэффициент диффузии, выражаемый
„ |
1 |
д (брі) |
величиной |
----- • ———, незначителен. |
|
RT |
д ln Cj |
В нестехиометрических фазах с сильно упорядоченными дефек тами (структуры сдвига) диффузия, несомненно, носит кооператив ный характер и не всегда требует участия точечных дефектов. Этим, в частности, объясняют необычайно высокую скорость взаимодействия Мо03 и Nb2Os, наблюдавшуюся Андерсоном [22], V
Представления о комплексном механизме реакционной диффу зии нашли развитие в работах Архарова с сотрудниками [23—27]. Они. обратили внимание на то, что значения кажущейся энергии активации, вычисленные из экспериментальных данных о темпе ратурной зависимости скорости процесса, в ряде случаев сущест венно меньше значений, которые получаются для того же процесса в другой области температур.
Пусть энергия образования MX достаточно велика и обозна чается символом Е в расчете на моль или &в расчете на молекулу MX, а энергия образования точечного дефекта (им может быть вакансия в металлической подрешетке или внедренный атом неме талла) равна U. При образовании в решетке точечных дефектов последняя искажается, причем из-за взаимного отталкивания меж ду одноименными дефектами суммарная энергия искажения тем выше, чем ближе друг к другу расположены дефекты
^ с б л > £ ^ .д ., 1
где f/сбл—энергия искажения решетки при образовании п дефектов
п
в достаточно малом объеме кристалла, а ^ Н т.д. — суммарная
1
энергия искажения решетки для случая, когда возникающие де фекты настолько удалены друг от друга, что их взаимодействием можно пренебречь.
При достаточно высоких температурах, когда энергия актива ции диффузии (qM и qx в расчете на 1 атом металла и неметалла соответственно) соизмерима с величиной kT, точечные дефекты, возникающие на границе окалина — газ, путем классических эле ментарных перескоков уходят в глубь окалины. Напротив, при более низких температурах, когда kT<^qM и kT<g.qx, диффузия путем элементарных перескоков практически невозможна, и точеч ные дефекты, образующиеся за счет химической реакции на гра нице окалина — газ, должны накапливаться в некоторой погра ничной зоне, называемой зоной химической реакции.
Вхождение атомов в зону стимулируется достаточно большим значением энергии АЕ (или Ае), выделяющейся при образовании новых связей между атомами М и X в зоне реакции. Можно полагать, что атомы X поступают в эту зону посредством абсорб ции, а атомы М выходят из наружных слоев кристаллической решетки. Накопление катионных вакансий (или внедрение избы точных атомов X в междоузлия решетки) создает в последней сильные искажения, благодаря чему зона химической реакции рас
ширяется. Условием |
существования таких искажений |
является |
Ае>Нт.д и А £> Нсбл- |
|
|
Очевидно, что расположение атомов в зоне реакции, сущест |
венно отличающееся |
от эквивалентного участка внутри |
решетки |
и характеризующееся энергией упругого искажения Ниск. не влия ет на более глубокие слои кристаллической решетки лишь до неко торого критического значения Ниск, равного Накт. При Ниск> [ /аІ(т начинается передача искажений вглубь окалины. Этот момент наступает, когда благодаря накоплению нестехиометрических де фектов в зоне реакции на границе кристаллическая решетка — зона реакции атомы оказываются настолько смещенными относи тельно регулярных узлов, что становятся возможными переходы этих атомов через энергетические барьеры, разделяющие положе ние атомов в регулярной решетке и в зоне реакции.
В результате атомных переходов произойдет разрядка упругих перенапряжений в области, расположенной вблизи наружной стороны окалины, и усиление напряжений в соседних участках, расположенных в глубине окалины. Иначе говоря, область иска жения как целое (своего рода многоатомная квазичастица или
19 ю . д . 1 реи-яков |
293 |
активированный комплекс) переходит из зоны реакции вглубь окалины путем последовательных перемещений атомов.
Этот процесс описывается Архаровым как передача упругого импульса искажений через кристаллическую решетку окалины. На границе окалина — металл происходит разрядка этого импульса путем вхождения в искаженную область атомов М, недостающих в ней до стехиометрического состава, или за счет удаления избы точных атомов X, внедренных в междоузлия. Тем самым устраня ются дефекты, образовавшиеся первоначально на границе окали н а— газ и эстафетно переданные к границе окалина — металл. Нет сомнений, что перемещение активированного комплекса приводит к массопереносу внутри окалины, так как оно связано с перекач кой определенного числа атомов М в направлении от металла к наружнему слою окалины или атомов X в противоположном на правлении.
Архаров [27] считает, что перемещение области искажений активированного комплекса через окалину, будучи по своей при роде бездиффузионным упругим процессом, в котором коллективно участвует большая группа атомов, может происходить со ско ростью, значительно превышающей скорость обычных диффузион ных процессов, слагающихся из независимых одноатомных элемен тарных актов. Вместе с тем время формирования активированного комплекса в зоне реакции и время его разрядки на противополож ной стороне окалины достаточно велико, так как эти стадии про цесса определяются классическими диффузионными актами. Это позволило Архарову [27], используя принцип локальных равно весий, дать термодинамическое ойосниьание участия активирован ных комплексов в массопереносе.
В заключение следует отметить, что упругоэстафетный меха низм реакционной диффузии может проявляться лишь в опреде ленном интервале температур (средние температуры). При высо ких температурах активированные комплексы не образуются, так как точечные дефекты могут легко перемещаться вглубь окалины путем элементарных независимых перескоков. При пониженных температурах упругоэстафетный массоперенос затруднен из-за медлительности элементарных актов, посредством которых проис ходит формирование активированного комплекса в исходной зоне химической реакции. Таким образом, участие активированных комплексов в процессе реакционной диффузии заметно сказывает ся на скорости реакции лишь в определенном интервале темпера тур, специфическом для каждой химической системы.
Итак, дефекты нестехиометрии играют первостепенную роль в диффузионных процессах, будь то элементарные перескоки ато мов или коллективное движение атомов в окислах с упорядочен
ными дефектами (кластеры, структура сдвига) |
или массоперенос |
в процессах реакционной диффузии с участием |
активированных |
комплексов. Доступные экспериментальные данные о зависимости
диффузии |
от |
нестехиометрии окислов |
подробно |
анализируются |
в главе II. |
|
|
|
|
|
Реакции между кристаллами и газовой фазой |
Наиболее |
важную |
в техническом |
отношении |
группу реакций |
«твердое |
тело — газ» с |
участием нестехиометрических окислов со |
ставляют процессы окисления металлов и сплавов. Значительный вклад в теорию высокотемпературного окисления металлов и спла вов был сделан Вагнером [28—30], работы которого стали фунда ментом для контролируемо го изменения скорости окис ления металлов и повыше ния их сопротивления к кор розии.
Теория Вагнера приме нима к плотной окалине, образованной продуктами реакции. В основе теории лежит предположение о том, что скорость роста окалины лимитируется диффузией в ней ионов, тогда как реак ции на границе фаз проис ходят достаточно быстро и
|
на поверхностях |
раздела |
|
|
|
|
|
|
«окисел — кислород» и «ме |
Р и с . 4.2. |
С х ем а т и ч еск о е |
и зо б р а ж е н и е |
|
талл — окисел» |
устанавли |
|
гр а д и е н т а |
к о н ц ен т р ац и й |
точ еч ны х д е |
|
вается локальное термоди |
ф ек тов и |
п р о ц есс о в |
п е р е н о с а , п р о и с |
|
намическое равновесие. |
х о д я щ и х |
в о к ал и н е |
M X , в |
к отор ой |
|
Движущей силой реак |
д о м и н и р у ю т к ати он н ы е |
в ак ан си и (а) |
|
или |
ан и он н ы е в ак ан си и |
(б ) |
|
ции является изменение сво |
|
|
|
|
|
|
бодной энергии при образо вании окисла МО из металла М и газообразного кислорода, в ре
зультате чего в окисле возникают концентрационные градиенты М и О. Скорость роста окалины определяется этими градиента ми, а также диффузионными характеристиками составных частей решетки, которые в свою очередь зависят от характера дефектно
сти структуры.
Чтобы выяснить роль дефектов нестехиометрии в процессе образования окалины, воспользуемся трактовкой Крёгера [31]. Рассмотрим в качестве примера образование бинарного соедине ния MX, характеризующегося разупорядочением типа Шоттки с преобладающей «-проводимостью. Очевидно, что термодинами ческое состояние соединения в пограничных слоях «MX — газ» и «MX—М» неодинаково (рис. 4.2). На границе «MX—М», обозна чаемой далее как граница I, устанавливается локальное равнове сие, которое можно выразить уравнением
19* 295
м , ет+ м й + ѵ£,
или, если Ѵх действует как донор, уравнением
Ммех ІИ Мм + Ѵх -\- e'.
В состоянии равновесия
Р ( М ) мет = Р (М)соед. I = І (М м )і Г S (Ѵх)і + I (<?')b |
(4.19) |
где |(і) — эффективный термодинамический потенциал і-струк- турного элемента.
На границе MX — газ (обозначаемой в дальнейшем как гра ница II) локальное равновесие можно выразить уравнением
Хгаз^5Хх + Ѵм,
или, если Ѵм действует как акцептор, уравнением
Хгаз^ Х х + Ѵм — ß'-
В состоянии равновесия
Р (Х )газ — р (Х)соед. II = ъ ( Х х ) ц + І ( Ѵ х ) іІ — I (е ' ) п ■ |
( 4 . 2 0 ) |
Учитывая, что |(і) = £°(і) + kT \u [і], легко убедиться, что
Р (Х)і — р (Х)ц = kT {ln [Ѵм]і •—• ln [VmJii + ln «и — ln m). (4.21)
Согласно этому уравнению движущей силой реакции является градиент концентраций точечных дефектов решетки, а не нормаль ных структурных элементов. Так как доминирующие атомные
дефекты решетки связаны соотношением [Ѵм] - [Ѵх] = As, то урав нение (4.21) можно преобразовать
Р (Х)і — р (Х)п = kT {ln [Vx]ii — ln [Vx]i + lnrtn — lntii}. (4.22)
Итак, процессы диффузии, ведущие к образованию соединения, описываются как результат перемещения дефектов.
В соответствии с обобщенным уравнением Фика поток заря женных частиц (или дефектов) типа і через единичную площадку, перпендикулярную направлению движения
|
/і |
А |
ЦТ |
&Чі |
(4.23) |
|
kT |
' 1 |
дх |
|
|
|
|
где Di — коэффициент диффузии |
частиц (или дефектов) |
типа і, |
|
а Гц — электрохимический потенциал. |
|
валент |
|
Учитывая, что, по определению, г]і= рі+ 2Гіі7ф, где Zi — |
ность перемещающихся частиц, q — заряд электрона, а ф — по тенциал поля, возникающего при движении заряженных частиц или дефектов, и используя уравнение (4.23), находим потоки дефектов, участвующих в диффузии
j(e') = - D e ^ |
+ neüe - f - , |
(4.24) |
|
OX |
o x |
|
;(V m) = - D |
m ^ ^ + [ Ѵ м ] о м - ^ - , |
(4.25) |
|
O X |
o x |
|
y'(Vx) =-- - D |
x ^ |
- [Vx] Ox - I е-. |
(4.26) |
|
ax |
ax |
|
Очевидно, что в стационарном состоянии потоки разнозарядных дефектов уравновешивают друг друга, т. е.
|
|
|
|
/ (е>) + |
/ (Ѵм) |
= І (V x). |
|
|
(4-27) |
|
Уравнения |
(4.24) — (4.27) |
в совокупности |
с уравнением Нерн- |
|
|
П. |
|
|
|
|
|
|
|
|
u |
ста—Эйнштейна —- = — !— позволяют |
рассчитать диффузионный |
потенциал |
кТ |
г\ц |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ом-а [ѵм1 |
|
|
|
|
дф |
|
|
КГ_ |
д[Ѵх] |
|
д п |
'і |
(4.28) |
|
|
|
|
|
|
|
|
|
|
|
|
д х |
|
|
a |
|
д х |
|
|
д х |
д х |
f ’ |
|
где о = q{vx [Ѵх] |
г г'м [Ѵм] |
+ |
ѵсп) |
(полная проводимость |
соедине |
ния). |
|
|
реакции |
(т. е. скорость |
образования |
молекул |
MX) |
Общая скорость |
определяется |
разностью |
потоков разноименных вакансий, т. е. |
|
|
|
|
/ р = |
/ (Ѵм) |
— / (Ѵ х ), |
|
|
(4.29) |
или, |
с учетом |
уравнения (4.27) |
|
|
|
|
|
|
|
|
|
|
|
/ |
р = |
І (е')- |
|
|
|
(4.30) |
|
Комбинируя уравнения (4.30), (4,24) и (4.28), находим |
|
|
|
|
|
/р |
кТ |
|
|
д I n |
[ Ѵ £ ] |
|
|
(4.31) |
|
|
|
|
Я2 |
|
|
д х |
|
|
|
|
|
|
|
|
|
|
|
|
где te — число переноса электронов, 6 = Ум-Их — суммарное число переноса катионов и анионов (или, что то же катионных и анион ных вакансий), а a — электропроводность кристалла. В соответ ствии с уравнением (4.31) скорость реакции определяется гради ентом концентрации анионных вакансий в неионизированной форме.
Сделанный вывод можно распространить на более общий случай, когда взаимодействие с парообразным неметаллом ведет к образованию соединения МаХь, характеризующегося дефектами
Ѵх , Ѵ м , e' и /г, причем атфЬг (т. е. степень ионизации не равна валентности составных частей кристалла). Легко показать, что в последнем случае