CC BY
56
УДК 669.112.228:669.788+538.91 -405
ВЗАИМОДЕЙСТВИЕ ВАКАНСИЙ И АТОМОВ, ОБРАЗУЮЩИХ В АЛЬФА-ЖЕЛЕЗЕ РАСТВОРЫ ЗАМЕЩЕНИЯ
А.А. Мирзоев, Г.Е. Рузанова, Д.А. Мирзаев, К.Ю. Окишев
Представлены результаты первопринципных расчётов в пакете WIEN2k энергии связи вакансии в а-железе с примесями Pd, Т^ V и Сг, располагающимися в первой - третьей координационных сферах вакансии. Энергия связи для Pd и Сг достигает 0,3 эВ, а для V и Сг составляет менее 0,1 эВ. При этом для атомов Pd и Сг она положительна и монотонно убывает по мере удаления от вакансии, а для И и V положительна в первой и третьей координационных сферах и отрицательна во второй. Разработана статистико-термодинамическая теория образования вакансий в разбавленных растворах замещения с учётом ближнего взаимодействия вакансий с растворёнными атомами. Расчёты показывают, что легирование железа палладием и титаном может приводить к многократному увеличению равновесной концентрации вакансий; эффект ванадия и хрома значительно слабее.
Ключевые слова: железо, вакансии, концентрация вакансий, примеси замещения, пер-вопринципные расчёты, термодинамика.
Вакансии играют важную роль во многих кинетических процессах в сплавах. Доля не заполненных атомами узлов кристаллической решётки,
- nV
или концентрация вакансий -J—, зависит прежде
N
всего от энергии их образования eV, которая представляет повышение полной энергии (точнее, энтальпии) кристалла при появлении одной вакансии. В выражении для свободной энергии чистого металла с вакансиями возрастание энергии компенсируется ростом конфигурационной и колебательной энтропии [1, 2], так что для определённой
( SV 1
(равновесной) концентрации порядка exp I-------I
I kT)
свободная энергия системы оказывается минимальной. Прекрасный обзор экспериментальных данных о величинах sv и методиках измерения приведён в [2].
В бинарных сплавах замещения, которые являются предметом рассмотрения в данной работе, вакансии, совершая случайные перескоки по решётке, могут оказаться вблизи атома растворённого элемента. Для некоторых элементов между их атомами и вакансиями возникает притяжение, а фактически захват, который как бы выводит такие вакансии из группы равновесных. Тем самым возрастает общее количество вакансий. Будем считать, что энергии образования вакансии вблизи и вдали от примесного атома X равны sV0 и svx соответственно. Тогда величина
Sb =SV0 ~SVX (1)
представляет энергию связи вакансии и растворённого атома. Обычно принимают, что sV0 совпадает с энергией образования вакансии в основе сплава, в данном случае в а-железе [2].
Экспериментальные измерения sb обычно проводят, сравнивая равновесные концентрации
вакансий в чистых металлах и сплавах на их основе. Погрешности таких измерений весьма значительны. Поэтому в настоящей работе было проведено первопринципное моделирование взаимодействия примесей замещения с вакансиями в а-железе, в основе которого лежат квантовомеханические методы расчёта энергии кристалла, содержащего дефекты, в рамках теории функционала электронной плотности.
Было изучено взаимодействие с вакансиями в а-железе атомов четырёх элементов: палладия, титана, ванадия и хрома. Расчёт электронной структуры проводился методом линейных присоединённых плоских волн (LAPW) с учётом обобщённого градиентного приближения (GGA) в программном пакете WIEN2k [3]. При моделировании были использованы следующие параметры: радиус тиЕйп-йп сферы атомов железа Rmt = 2,0 а.е.; параметр, контролирующий сходимость, RmtKmax = 10; количество к-точек, характеризующее дискретную сетку, используемую для численного интегрирования по зоне Бриллюэна, пк = 27. Был выбран параметр решётки а = 2,84 А, являющийся равновесным для чистого ОЦК-железа [4]. В качестве модели была выбрана суперячейка, состоящая из 54 атомов. Для оценки энергии образования моновакансии в ОЦК-железе примерно в центре суперячейки был удалён один из атомов. Было определено значение энергии образования вакансии в чистом железе еу = 2,15 эВ. Далее один из атомов Fe, находящийся в первом, втором либо третьем окружении относительно вакансии, поочерёдно замещался атомом примеси. Для получения более точного результата система, содержащая вакансию и примесь, была отрелаксирована.
Энергия образования точечного дефекта любого вида А определялась как разница между энергиями суперячейки, содержащей дефект, и суперячейки чистого железа:
в(А) = Е(А) - Ем - Е(А^):
(2)
где Е(А) - энергия системы, содержащей дефект, Еге{ - энергия чистого железа, а Е(Аге£) - энергия одиночного дефекта. В частности, при расчёте энергии образования вакансии в системе с примесью замещения X использовалась формула [5]:
ЄУХ ~ Е¥е52-Х | ЕРе53-Х 54 ЕРе54
(3)
где ЕРе52-х - энергия системы из 52 атомов Fe + вакансия + 1 атом X; ЕРе53-х - энергия системы из 53 атомов Fe + 1 атом X (без вакансии), а ЕРе54 -энергия системы из 54 атомов чистого железа.
На рис. 1 представлена зависимость энергии образования вакансии в ОЦК-железе от положения атома палладия. Добавление примеси Pd к ОЦК-же-лезу снижает энергию образования вакансии; следовательно, образуется связанный комплекс У-Рй. Взаимодействие между вакансией и растворённым
атомом наиболее сильно тогда, когда они расположены в ближайших соседних узлах решётки. Энергия связи в этом случае максимальна и равна
0,30 эВ. Если переместить атом палладия во вторую, третью и т. д. координационные сферы по отношению к вакансии, то энергия связи заметно уменьшается (см. таблицу и рис. 1).
Энергии связи с вакансиями в различных координационных сферах для атомов остальных трёх элементов (Т^ V, Сг) даны в таблице и на рис. 2. Отметим, что энергия связи и в этих случаях максимальна для позиций ближайших соседей. Удивительно, что во второй координационной сфере энергия связи У-Т и отрицательна, что отражает отталкивание, то есть перескок атомов в первую или третью сферу. У атома Сг энергия связи положительна и для второй сферы, но её значение мало. Полученные нами результаты согласуются с данными других расчётных работ [6, 7], а также
Рис. 1. Энергии образования вакансии при её расположении в различных координационных сферах относительно атома Pd в системе Fe52 + вакансия + Pd (пунктирная линия соответствует энергии образования вакансии в чистом ОЦК-железе)
Энергии связи комплекса вакансия - атом замещения M
Атом замещения Номер координационной сферы Энергия связи, эВ
Данная работа [6] [7] Эксперимент [8]
1 0,30 0,36
Pd 2 0,25 - 0,26 0,21
3 0,14
1 0,28 0,18
Ті 2 -0,15 0,22 -0,11 0,16
3 0,04
1 0,08 0,04
V 2 -0,07 0,04 -0,06 < 0,105
3 0,03
Сг 1 2 3 0,07 0,04 0,03 0,05 0,05 0,03 < 0,105
0,3 г
0,2
0,0
-0,1
-0,2
□
Ті
□
V
гь
Cr
Pd
I I Первое окружение
I I Второе окружение
I I Третье окружение
Рис. 2. Энергия связи вакансии и атома замещения, расположенного относительно неё в первой, второй и третьей координационной сфере
экспериментальными данными [1, 8]. Так, на основании данных об электросопротивлении сплавов Fe-Cг [1], энергия образования вакансии в железе увеличивается на 0,05 эВ на каждый атомный процент хрома; энергии связи примесей с вакансиями, найденные методом положительных мюонов [8], приведены в таблице.
Между вакансией и атомом замещения может возникать как упругое, обусловленное несоответствием размеров примеси и матричного атома, так и магнитное взаимодействие. Уменьшение энергии системы при добавлении к железу примесей, радиус которых превосходит радиус его атомов, связано с частичной компенсацией поля возмущений, создаваемого вакансией. Если бы взаимодействие между вакансией и атомом замещения было чисто упругим, то энергии связи вакансии с Т и V имели бы большие значения, чем с Рй, так как радиус Т (1,47 А) превосходит радиус Рй (1,37 А), а радиус V сопоставим с ним (1,35 А); напомним, что радиус железа составляет 1,24 А. При этом Т и V проявляют в такой системе антиферромагнитные свойства, тогда как атом Сг практически не взаимодействует с вакансией ввиду его близости к атомам Fe как по размеру, так и по свойствам. Таким образом, сложный характер взаимодействия примесей с вакансией вызван суперпозицией двух сопоставимых вкладов: упругого и химического.
Полученные данные об энергиях взаимодействия примесных атомов с вакансиями позволяют рассчитать эффект изменения равновесной концентрации вакансий при легировании железа исследованными элементами.
Предположим, что в твёрдом растворе замещения железа и элемента X, содержащем Аре и Ых атомов этих элементов, возникает пУ вакансий. Их
атомные доли тогда равны х^ =
N
Ее
N
N
X
N
М — МЕе + Мх + пу .
взаимодействием ближайших соседей и введём вероятности событий, что данный ближайший узел решётки, соседствующий с вакансией, занят атомом железа:
N
N
Чу
или атомом элемента X:
N
у - х
у - X
(4а)
(4 б)
где ЫУ_Ре и ЫУ_х представляют числа соответствующих атомных пар, а г - координационное число решётки. Поскольку концентрация вакансий мала, то соседством вакансия-вакансия можно пренебрегать, и тогда ЫУ_Ре + ЫУ_х = гЫхУ, а
РУ _Бе + ру _ х =1.
В таком случае изменение энергии кристалла при образовании вакансии должно отражать лишь эффекты её взаимодействия с атомами Fe и X. Обозначим через еь энергию связи вакансии и растворённого атома X (положительную при притяжении и отрицательную при отталкивании). Тогда изменение энергии кристалла при образовании пУвакансий
АЕ = ПуеУ _ znyPy_хеь , (5)
где еУ - энергия образования вакансии в чистом железе.
Происходящее при этом изменение свободной энергии твёрдого раствора, как обычно, равно
АО — АЕ - ИГ 1п Ж .
(6)
где Ж - термодинамическая вероятность рассматриваемой атомной конфигурации. Основываясь на статистически-термодинамическом подходе [9, 10], запишем Ж в следующем виде:
Ж —-
N!
Ж1 -Ж2
(Ы_Пу)!пу! Ж1неуп-Ж2неуп (7)
Здесь через Жх обозначено число способов разме-
гЫ-ре парам, которые
Ограничимся щения Nу-Ее пар У-Ее по
и
п
атомы железа образуют с соседними узлами решётки, а через Ж2 - число способов размещения Nу-X пар У-X по ІNX таким же парам, образуемым атомами X:
W =
w2 =
( zNFe ) I
(zNFe NFe-V ) ! NFe- V !
(zNx)!
(8а)
(8б)
(гЫх _Ых_у)!Ых_у !
При этом ЫУ_Ре и ЫУ_х определяются соотношениями (4), в которых значения вероятностей РУ_Ре и РУ_х отражают существующий в растворе ближний порядок в расположении атомов замещения и вакансий. В условиях же полностью неупорядоченного размещения вакансий и атомов X по решётке (Ж1неуп и Ж2неуп) вероятности равны атомным долям: РУ_Ре = хРе; РУ_х = хх .
Подставив эти значения в (8), (7) и далее в (6), получим с учётом формулы Стирлинга 1п п !« п 1п п _ п общее выражение для изменения свободной энергии раствора при образовании пу вакансий:
АG = пуеУ _ гпУРУ_хеь _ кТ [Ы 1п Ы _ (Ы _ пу ) х х 1п(Ы _ пу ) _ пу 1п пу _ (гЫре _ гпУ + гпУРУ_х) х х 1п(гЫ¥е _ гПу + гПуРу_х) _ (гПу _ гПуРу_х) х х 1п (ту _ гПуРу_х ) _ (гЫх _ гПуРу_х) х х 1п(гЫх _ гПуРу_х) _ гПуРу_х 1п(гПуРу_х) +
+(гЫРе _ гпух¥е) 1п( гЫ¥е _ гпу хРе) + гПу хРе х 1п(гпУхРе) + (гЫх _ гпУхх) 1п(гЫх _ гпУхх) +
+ гпухх 1п(гпухх)]. (9)
Равновесное значение РУ_х тогда можно най-
ти из соотношения
SAG
dP
= 0 . После дифферен-
V - X
цирования это соотношение приводит к квадратному уравнению относительно Ру-X, положительный корень которого равен
хЕе h - ху (h +1) +1
P
V-x (
-1 + 11 + -
2 xV h
4xx Xv h
(xFeh - xV (h +1) +1)
(10)
где обозначено h = exp |-----|-1. Поскольку
kT
xV ^ 1, можно использовать разложение \І1 + x да да 1 + x/2, и тогда выражение (10) приводится к виду
P
V-x
1 + xFeh
(11)
Это означает, что равновесные вероятности РУ_х и РУ_-ре практически не зависят от концентрации
вакансий, что значительно облегчает дальнейшее рассмотрение.
Теперь найдём равновесное значение концентрации самих вакансий, которое определяется также из условия минимума свободной энергии dAG Л ^
-----= 0. После дифференцирования с учётом
3%
соотношения (11) и малости xV получим выражение
XV = exp| -kTJI1 ~ Xx + Xx expI kT
в котором сомножитель
x° = exp
є
V
kT
(12)
(13)
представляет долю вакансий в чистом железе.
Результаты нашего первопринципного моделирования дали для энергии образования вакансии значение еУ = 2,15 эВ. Авторами [6] и [11] было получено значение 2,17 эВ. Следует отметить, что эти величины находятся на верхнем пределе диапазона значений, наблюдавшихся экспериментально (от 1,4 + 0,1 до 2,0 + 0,2 эВ [12, 13]). А. Зеегер в обзоре [14] на основании литературных данных и собственных экспериментов пришёл к значениям еУ = 1,61...1,75 эВ, а из результатов работы [15], где непосредственно измерялась концентрация вакансий после закалки от 1040 К, получил еУ = 1,59 эВ. М.А. Штремель [2] считает наиболее надёжным значение еУ = 1,6 эВ. Авторы [12] показали, что повышение концентрации углерода в а-железе с 2 10-5 до 110-3 мас. % приводит к снижению наблюдаемой энергии образования
вакансии с 1,79 + 0,10 до 1,53_1,60 эВ; согласно
первопринципным расчётам [6], энергия связи вакансии с атомом углерода составляет 0,60 эВ. Возможно, именно поэтому рассчитанное для чистого железа значение еУ выше, чем в большинстве экспериментов, выполненных на железе технической чистоты. Кроме того, точное выражение для доли вакантных узлов имеет вид
AG
~kT
xv = exp| | = exp| -— |exp| -k" I, (14)
AS
k
где ДО - изменение энтропии колебательной подсистемы, обусловленное тем, что частота колебаний атомов вблизи вакансии (у') и вдали от неё (у) значительно различается; при использовании мо-
у
дели независимых колебаний атомов А£ = гк 1п — .
у'
Следовательно, определяемая в эксперименте эффективная энергия образования вакансии, которая и входит в формулу (13), будет меньше истинной. Поэтому в дальнейших расчётах мы использовали экспериментальное значение еУ = 1,6 эВ.
Результаты расчётов температурной зависимости концентрации вакансий для двойных спла-
z
X
x
T, К T, К
Т, К Т, К
Рис. 3. Температурная зависимость доли вакантных узлов ху, рассчитанная по формуле (12) при еу = 1,6 эВ и энергиях связи еь из таблицы для различных концентраций примесей Pd (а), Т (б), V (в) и Сг (г)
вов Fe с Рй, Ti, V и Сг разной концентрации приведены на рис. 3. Из него видно, что притяжение вакансий к примесным атомам приводит к тому, что в системах Fe-Pd и Fe-Ti, где энергия связи примеси с вакансией относительно велика (см. таблицу), заметно возрастает их равновесная концентрация; более того, при понижении температуры снижение концентрации вакансий может даже смениться ростом (рис. 3, а, б).
В системах же Fe-V и Fe-Cг, где энергии связи примеси с вакансией в 3-4 раза меньше, даже значительно большие концентрации примеси оказывают относительно небольшое воздействие на равновесную концентрацию вакансий: она возрастает примерно на порядок при концентрации хх = 0,25 (впрочем, не очевидно, что значения энергии связи еь из таблицы, рассчитанные для концентрации хх = 1/54 и 0,019, сохранятся неиз-
менными вплоть до таких высоких значений концентрации примеси).
Работа выполнена в рамках Государственного задания Министерства образования и науки РФ и частично поддержана грантом РФФИ 13-02-00048_а.
Литература
1. Криштал, М.А. Механизм диффузии в железных сплавах /М.А. Криштал. - М.: Металлургия, 1972. - 400 с.
2. Штремель, М.А. Прочность сплавов: учеб. для вузов. Ч. I: Дефекты решётки / М.А. Штремель. - М.: МИСиС, 1999. - 384 с.
3. Schwarz, K. Solid State Calculations Using WIEN2k / K. Schwarz, P. Blaha // Computational Materials Science. - 2003. - Vol. 28, no. 2. - P. 259-273.
4. Урсаева, А.В. Выбор оптимальных параметров для построения максимально точной мо-
дели ОЦК-железа / А.В. Урсаева, Г.Е. Рузанова, А.А. Мирзоев // Вестник ЮУрГУ. Серия «Математика, физика, химия». - 2010. - Вып. 2, № 9. -С. 97-101.
5. Terentyev, D. / Diffusion of 3D-migrating selfinterstitial clusters in diluted and concentrated Fe-Cr alloys /D. Terentyev, P. Olsson, L. Malerba // Journal of Nuclear Materials. - 2009. - Vol. 386-388. -P. 140-152.
6. Ohnuma, T. First-principles calculations of vacancy-solute element interactions in body-centered cubic iron / T. Ohnuma, N. Soneda, M. Iwasawa. // Acta Materialia. - 2009. - Vol. 57, no. 20. - P. 59475955.
7. Vacancy-solute interactions in ferromagnetic and paramagnetic bcc iron: Ab initio calculations / O.I. Gorbatov, P.A. Korzhavyi, A.V. Ruban et al. // Journal of Nuclear Materials. - 2011. - Vol. 419, no. 1-3. - P. 248-255.
8. Doyama, M. Study of metallic properties by positive muons / M. Doyama // Bulletin of the Japan Institute of Metals. -1986. - Vol. 25, no. 10. - P. 808-815.
9. Fowler, R.H. Statistical Thermodynamics /
R.H. Fowler, E.A. Guggenheim. - Cambridge University Press, 1949. - 701 p.
10. Лесник, А.Г. Улучшенный вариант квази-химической теории упорядочения сплавов /А.Г. Лесник // Физика металлов и металловедение. - 1977. -Т. 44, вып. 5. - С. 945-955.
11. Counts, W.A. First-Principles Energetics of Hydrogen Traps in a-Fe: Point Defects / W.A. Counts, C. Wolverton, R. Gibala // Acta Materialia. - 2010. -Vol. 58, no. 14. - P. 4730-4741.
12. Positron Annihilation on Pure and Carbon-Doped a-Iron in Thermal Equilibrium. / L. De Schep-per, D. Segers, L. Dorikens-Vanpraet et al. // Physical Review B. - 1983. - Vol. 27, no. 9. - P. 5257-5269.
13. Atomic Defects in Metals / ed. by H. Ull-maier. - Landolt-Bornstein, New Series, vol. III/25. -Springer-Verlag, 1991. - P. 128.
14. Seeger, A. Lattice Vacancies in High-Purity a-Iron/A. Seeger //Physica Status Solidi (a). - 1998. -Vol. 167, no. 2. - P. 289-311.
15. Seydel, O. Quenching-in of Vacancies in Pure a-Iron / O. Seydel, G. Frohberg, H. Wever // Physica Status Solidi (a). - 1994. - Vol. 144, no. 1. - P. 69-79.
Мирзоев Александр Аминулаевич, доктор физико-математических наук, профессор кафедры общей и теоретической физики, Южно-Уральский государственный университет. 454080, г. Челябинск, пр. Ленина, 76. Тел.: (351)2654713. E-mail: [email protected].
Рузанова Галина Евгеньевна, аспирант кафедры общей и теоретической физики, Южно-Уральский государственный университет. 454080, г. Челябинск, пр. Ленина, 76. Тел.: (351)2654713. E-mail: [email protected]. ac.ru.
Мирзаев Джалал Аминулович, доктор физико-математических наук, профессор кафедры физического металловедения и физики твёрдого тела, Южно-Уральский государственный университет. 454080, г. Челябинск, пр. Ленина, 76. Тел.: (351)2679013. E-mail: [email protected].
Окишев Константин Юрьевич, кандидат физико-математических наук, доцент кафедры физического металловедения и физики твёрдого тела, Южно-Уральский государственный университет. 454080, г. Челябинск, пр. Ленина, 76. Тел.: (351)2679013. E-mail: [email protected].
Bulletin of the South Ural State University Series "Metallurgy”
___________________________________________2013, vol. 13, no. 1, pp. 108-114
INTERACTION OF VACANCIES WITH SUBSTITUTIONAL SOLUTE ATOMS IN ALPHA IRON
A.A. Mirzoev, G.E. Ruzanova, D.A. Mirzaev, K.Yu. Okishev
The paper presents the results of ab initio calculations in WIEN2k package of binding energy of a vacancy in a iron with Pd, Ti, V and Cr impurity atoms located in the first, second and third spheres of coordination with respect to the vacancy. The binding energy reaches 0.3 eV for Pd and Cr being less than 0.1 eV for V and Cr. It is positive and monotonically decreases with increasing distance from the vacancy to Pd and Cr atoms, while for Ti and V atoms it is positive in the first and the third spheres of coordination and negative in the second sphere. A statistical and thermodynamic theory of vacancy formation in dilute substitutional solutions is proposed that takes into account the short-range interaction of vacancies with impurity atoms. Calculations show that palladium and titanium additions to a iron may lead to multiple increase of equilibrium vacancy concentration, while the effect of vanadium and chromium is much smaller.
Keywords: iron, vacancies, vacancy concentration, substitutional impurities, ab initio calculation, thermodynamics.
Mirzoev Aleksandr Aminulaevich, doctor of physical and mathematical science, professor of the General and Theoretical Physics Department, South Ural State University. 76 Lenin avenue, Chelyabinsk, Russia 454080. Tel.: 7(351)2654713. E-mail: [email protected].
Ruzanova Galina Evgen’evna, post-graduate student of the General and Theoretical Physics Department, South Ural State University. 76 Lenin avenue, Chelyabinsk, Russia 454080. Tel.: 7(351)2654713. E-mail: [email protected].
Mirzaev Dzhalal Aminulovich, doctor of physical and mathematical science, professor of the Physical Metallurgy and Solid State Physics Department, South Ural State University. 76 Lenin avenue, Chelyabinsk, Russia 454080. Tel.: 7(351)2679013. E-mail: [email protected].
Okishev Konstantin Yur’evich, candidate of physical and mathematical science, associate professor of the Physical Metallurgy and Solid State Physics Department, South Ural State University. 76 Lenin avenue, Chelyabinsk, Russia 454080. Tel.: 7(351)2679013. E-mail: [email protected].
Поступила в редакцию 25 февраля 2013 г.
