CC BY
48
УДК 669.112.227.1:538.915
AB-INITIO МОДЕЛИРОВАНИЕ ВЛИЯНИЯ БЛИЖНЕГО ОКРУЖЕНИЯ ПРИМЕСЕЙ УГЛЕРОДА НА ЭНЕРГИЮ ИХ РАСТВОРЕНИЯ В ГЦК-ЖЕЛЕЗЕ1
Я.М. Ридный2, A.A. Мирзоев3, Д.А. Мирзаев4
В программном пакете WIEN2k проведено первопринципное моделирование равновесной структуры и свойств ГЦК-железа. Получены оптимальные параметры, позволяющие построить наиболее точную модель. Произведено вычисление энергий растворения атомов углерода, их взаимного расположения и вклада упругих эффектов в энергию системы для немагнитного (НМ) и антиферромагнитного двухслойного (АФМД) состояний ГЦК железа.
Ключевые слова: первопринципное моделирование, ГЦК-железо, примеси углерода, WIEN2k.
Введение
Сплавы на основе железа, несмотря на развитие производства пластмасс, керамики, композитов, остаются основными конструкционными материалами. Исключительным достоинством чёрных металлов является, прежде всего, их дешевизна, а также универсальность — при помощи небольших добавок можно резко изменить свойства их сплавов. Например, внедрение всего лишь 1 ат. % углерода радикально повышает прочность железа. Железо является полиморфным металлом [1, 2]: при нагревании выше 911 °С низкотемпературное а-железо с ОЦК решеткой превращается в у-железо с ГЦК решеткой. Углерод растворяется в обеих фазах железа, так что при температурах выше 723-911 °С существует аустенит - раствор углерода в y-Fe, ниже этого интервала возникает феррит -раствор углерода в a-Fe [2]. Растворимость углерода в y-Fe довольно высока - до 2 масс. %. Напротив, a-железо растворяет очень мало углерода. Наибольшая растворимость (0,025 масс. % ) наблюдается при 723 °С. В ходе медленного охлаждения аустенита он превращается в феррит, а избыток углерода сверх равновесной концентрации выделяется в виде карбида Fe3C. Поэтому превращение аустенита при охлаждении стали, помимо перестройки решетки, требует прохождения процесса диффузии углерода в аустените и феррите.
Углерод образует твердый раствор внедрения как с у-, так и с a-железом, что было установлено довольно давно [1]. По результатам совместного измерения плотности р и параметра кристаллической решетки а для серии сплавов удается рассчитать число атомов п, приходящихся на одну элементарную ячейку:
п = ра3 Na / A, (1)
где А - средняя атомная масса, а Na - число Авогадро. Наблюдавшееся возрастание п при увеличении содержания углерода в аустените и феррите однозначно свидетельствовало об образовании углеродом растворов внедрения в решетках железа. В структуре у- и a-железа существуют пустоты (поры), ограниченные атомами железа, расположенными в вершинах октаэдра или тетраэдра, в связи с чем их называют октаэдрическими и тетраэдрическими порами. Вопрос о том, в порах какого типа расположены атомы углерода в решетках y-Fe и a-Fe, вызвавший многолетнюю дискуссию, был решен на основе данных экспериментов по ядерному гамма-резонансу и нейтронографии [2]. Было показано, что атомы углерода в обеих решетках располагаются в октаэдрических порах (ОП).
1 Работа «Многомасштабное моделирование поведения водорода в сплавах железа с целью разработки энергоэффективных методов борьбы с водородными трещинами в сталях» выполнена в рамках госбюджетной НИР в Южно-Уральском государственном университете.
2 Ридный Ярослав Максимович - студент, Южно-Уральский государственный университет.
E-mail: yaroslav@physics.susu.ac.ru
3 Мирзоев Александр Аминулаевич - доктор физико-математических наук, профессор, кафедра общей и теоретической физики, Южно-Уральский государственный университет.
E-mail: mirzoev@physics.susu.ac.ru
4 Мирзаев Джалал Аминулович - профессор, доктор физико-математических наук, кафедра физического металловедения и материаловедения, Южно-Уральский государственный университет.
E-mail: mirzaev@physmet.susu.ac.ru
Ридный Я.М., Мирзоев A.A., Ab-initio моделирование влияния ближнего окружения
Мирзаев Д.А. примесей углерода на энергию их растворения в ГЦК-железе
В условиях быстрого охлаждения аустенита, когда диффузия углерода подавлена, происходит бездиффузионное мартенситное превращение, при котором все атомы углерода принудительно переходят в решетку мартенсита, а последняя, хотя и похожа на решетку a-железа, имеет слабую тетрагональность [1-3]. Возникновение тетрагональности и бездиффузионный характер мартенситного ГЦК-ОЦК превращения был объяснен Бейном, который обратил внимание, что ГЦК-решетку можно превратить в объемно-центрированный куб при небольшом сжатии вдоль оси OZ и растяжении в направлениях осей OX и OY (бейновская деформация). Октапоры ГЦК-решетки вместе с располагающимися в них атомами углерода переходят только в октапоры типа «Z» решетки ОЦК-железа. При этом все атомы углерода растягивают решетку в одном Z-направлении, и она из кубической превращается в тетрагональную. Хотя теория Бейна общепри-знана, вопрос теоретического обоснования ее справедливости остается открытым.
В связи с этим, вопрос о распределении атомов углерода в ГЦК-решетке железа и их упорядочении остаётся очень важным. Очевидно, что мартенсит наследует расположение атомов углерода в аустените, которое, в свою очередь, зависит от характера взаимодействия соседних атомов углерода. Среди многочисленных характеристик аустенита наиболее чувствительной к особенностям распределения атомов углерода по октапорам оказалась активность углерода в у-фазе [4]. Еще Вагнер с сотрудниками [5] показали, что при отсутствии взаимодействия между атомами углерода величина активности дается выражением:
а=тег, (2)
1 -РХс
где xc - доля атомов углерода в сплаве, а коэффициент ¡3 = 2. Однако результаты экспериментального изучения активности углерода в у-железе показали, что величина с увеличением содержания углерода возрастает сильнее, чем это предписывает формула (2). Было высказано предположение, что это обстоятельство связано с наличием отталкивания между ближайшими атомами углерода, что и приводит к уменьшению числа возможных мест их расположения. Даркен и Смит [6] математически оформили данную идею, введя величину s - энергию взаимодействия рассматриваемого атома углерода с другими атомами С, расположенными в любой из 12 октапор ближайшего окружения. Окончательный результат для активности оказался таким:
агс =-----Х-------, (3)
с 1 - (14 - Ш)хс
где в = exp(-e/kT). Сопоставление с экспериментальными данными для активности углерода позволило определить величину s = 0,065 эВ/атом, что для 1200 К дает /3 = 7,6. Видно, что учет отталкивания действительно увеличивает активность углерода в растворе.
Темкиным и Шварцманом [7], а затем Кауфманом [8] была развита модель блокировки, в которой выдвинуто предположение, что каждый атом С блокирует, т.е. делает недоступным для заполнения другими атомами углерода, некоторое число m октапор, ближайших к нему. В этом варианте:
x
lnacY = ln---с----. (4)
С 1 - (m + 1)хс V '
Параметру m придавали значение 5 в [9], и 4 в [7]. Последнее значение можно интерпретировать в модели, что атомы углерода занимают только октапоры, расположенные в центре кубов. Квазихимическая модель (модель мягкого отталкивания) была разработана В.А. Кожеуровым [10], и позднее Мак-Лелланом [9]. В ней принимается, что атом углерода отталкивает с энергией £i любой атом углерода, расположенный в первой координационной сфере исходного атома. Совпадение теории и экспериментальных данных достигается при е1 = 0,085 эВ/атом [11]. Позднее Большовым Л.А. [12] была развита более общая теория, учитывающая энергии взаимодействия пар атомов углерода не только в первой, но и последующих координационных сферах.
Поскольку все рассмотренные и сильно различающиеся модели позволяют довольно точно описать активность углерода в аустените, можно заключить, что данные об активности не позволяют провести детальную реконструкцию распределения углерода в ГЦК-железе.
Важные данные о распределении углерода получены в экспериментах по ядерному гамма-резонансу (ЯГР). ЯГР-спектр аустенита имеет особенности, заключающиеся в существовании
центрального максимума от атомов железа, не имеющих в своем окружении атомов углерода, и дублета от атомов железа, соседствующих с углеродом в ближайших междоузлиях. Обзор мес-сбауэровских исследований распределения углерода в аустените был сделан Гаврилюком В.Г. [13] и Надутовым В.М. [14]. Было показано, что, хотя в выводах исследователей имеются противоречия, в основном большинство авторов склоняются в пользу модели жесткого отталкивания, т.е. блокировки октапор не только в первой, но и во второй координационной сфере. Этот вывод согласуется с результатами расчета энергии взаимодействия пары атомов углерода при их расположении в первом и втором соседствах - +0,15 эВ и 0,4 эВ, полученными Гаврилюком В.Г. методом функций Грина [13]. В диссертации Надутова В.М. [14] методом Монте-Карло были получены области возможных значений энергий взаимодействия атомов углерода в первой и второй координационных сферах на основании мессбауровских экспериментов. Было показано, что энергия отталкивания атомов углерода во второй координационной сфере в разы превышает отталкивание в первой сфере и превосходит величину 0,08 эВ.
Таким образом, экспериментальными методами довольно легко определить концентрацию углерода, но достоверно установить позиции атомов углерода в решетке железа, а также структуру их ближнего порядка на данный момент времени представляет очень сложную и дорогостоящую задачу. Поэтому в последние годы все большее внимание уделяется методам перво-принципного компьютерного моделирования, главным достоинством которых является отсутствие каких-либо входных данных, влияющих на достоверность результатов при использовании минимального количества подгоночных параметров.
Джианг и др. исследовали энергию растворения углерода в ферромагнитном(ФМ) состоянии ГЦК железа в программном пакете VASP [15]. Энергия растворения вычислялась относительно графита и составила -0,17 эВ. Этот результат оказался существенно ниже экспериментального значения 0,36 эВ, полученного в работе [16] методом рассеяния электронов. В работе [17] в программном пакете SIESTA проводилось моделирование растворения углерода в АФМД состоянии ГЦК-железа, энергия составила 0,55 эВ. В ряде работ [18-20] рассматривали растворение второго атома углерода в ГЦК-решетке железа и показали, что атомы углерода предпочитают располагаться на максимально возможном расстоянии. Наиболее полный анализ растворения углерода проведён в работе Медведевой [18], но она ограничилась только немагнитным (НМ) состоянием железа.
Видно, что несмотря на значительное количество работ, посвященных изучению растворения углерода в железе, вопросы о расположении атомов углерода, энергиях их растворения, а также энергиях взаимодействия между атомами углерода остаются открытыми. В связи со всем вышесказанным, целью данной работы являлось ab-initio моделирование влияния расположения соседних примесей на энергию растворения атома углерода в различных магнитных фазах ГЦК-железа с помощью программного пакета WIEN2k [21].
Реальное ГЦК-железо существует в парамагнитном состоянии, при температурах порядка 1000 °С. Его магнитный момент равняется 0,6 цЬ [22]. Моделирование парамагнитного состояния методами зонной теории в пакете WIEN2k довольно затруднительно, поскольку они позволяют проводить моделирование лишь коллинеарных магнитных состояний. Поэтому мы будем описывать парамагнитное состояние с помощью двух приближений - немагнитного (НМ), в котором магнитные моменты равны нулю, и в антиферромагнитном двухслойном состоянии (АФМД). Согласно работе [23] оно является наиболее стабильным из коллинеарных магнитных состояний в ГЦК -железе.
Методы
В данной работе расчёты проводись первопринципным полнопотенциальным методом LAPW, с учетом обобщенного градиентного приближения PBE-GGA в программном пакете WIEN2k, обеспечивающим высокую точность расчета полной энергии при минимальном количестве подгоночных параметров.
В методе LAPW [24] волновые функции, зарядовая плотность и потенциал разлагаются по сферическим гармоникам в неперекрывающихся атомных сферах радиуса Rmt и по плоским волнам в остальной области элементарной ячейки. Предполагается, что внутри маффин-тин сферы зарядовая плотность сферически-симметрична. Волновые функции в междоузлии раскладываются по плоским волнам с максимальным волновым вектором Kmax. Эта величина является одним
Ab-initio моделирование влияния ближнего окружения примесей углерода на энергию их растворения в ГЦК-железе
из основных параметров, влияющих на точность, так как она определяет число базисных функций (размер матриц). Для интегрирования в обратном пространстве и вычисления электронной плотности использовалась схема Монхорста-Пака с сеткой из N точек в зоне Бриллюэна. Чем больше используется точек, тем выше точность интегрирования, однако при этом возрастает расчетное время, поэтому оптимальное значение определяется из условия достижения требуемой точности энергии системы.
В нашей работе использована ГЦК-суперячейка, состоящая из 32 атомов. Была проведена оптимизация параметров моделирования, которая показала, что для достижения требуемой точности расчета полной энергии в 1 мРб/систему, необходимо использование следующих параметров: Ят1(Бе) = 2,0 а.е., Ктах = 5 а.е.-1. Эти значения совпадают с параметрами, определёнными в работах [25, 26] для ОЦК железа. Для определения геометрических параметров необходимо использование N = 64 точки, а чтобы гарантировать точность вычисления полной энергии в 1 мРб, необходимо его увеличить до Ык =343.
Для определения энергии растворения углерода относительно фазы графита использовалась формула:
АН = Е (Бе32 С) - Е (Бе32) - Е (С), (5)
где Е(Бе32С1) - энергия отрелаксированной суперячейки, состоящей из 32 атомов железа и одного атома углерода, находящегося в октапоре, Е(Бе32) - энергия суперячейки состоящей из 32 атомов железа, а Е(С) - энергия одного атома углерода в решётке графита. Для нахождения энергии растворения второго атома углерода формула (1) преобразуется:
АН 2 = Е (Бе32С2) - Е (Бе32 С) - Е (С). (6)
В формуле (6) Е(Бе32С2) - энергия релаксированной суперячейки, состоящей из 32 атомов железа и двух атомов углерода находящихся в октапорах.
Для нахождения энергии взаимодействия между атомами углерода из энергии системы с двумя атомами углерода вычитается энергия двух суперячеек, содержащих одиночные атомы углерода, и суммируется с энергией беспримесной суперячейки ГЦК железа:
АУ = Е(Бе32С2) -2Е(Бе32С) + Е(Бе32) . (7)
Для нахождения энергии углерода в графите мы использовали значения параметров решётки, полученные из рентгеноструктурных данных (а = 2,462 А, с = 6,656 А [15]). В расчете использовались оптимизированные параметры Ктах = 5 а.е.-1, Ит1:(С) = 1,2 а.е. Количество к-точек было увеличено по сравнению с работой [15] для обеспечения точности в 1 мРб/систему от Ык = 256 до Ык = 300.
Растворение углерода и деформация решётки
Для получения равновесных параметров суперячейки ГЦК-железа с растворенными атомами углерода была проведена атомно-силовая релаксация системы и объемная минимизация энергии.
Была получена зависимость параметра решётки от концентрации углерода, приведенная на рис. 1.
При ее построении мы использовали результаты моделирования с растворенными одним и двумя атомами углерода.
Было получено, что растворение одиночного атома С в НМ ГЦК-железе увеличивает параметр решётки с 3,45 А до 3,475 А и расталкивает ближайшие атомы железа на 7,5 %. При растворении углерода в АФМД состоянии происходит увеличение параметра решётки с 3,54 А до 3,555 А, и образование в первом окружении слегка несимметричного октаэдра, у которого апикальный атом с противоположным направлением спина относительно других атомов железа
3,65
3,4 >—
0 2 4 6 8
Концентрация углерода в ат.%
Рис. 1. График зависимости параметра решётки от концентрации углерода. Штриховая линия - АФМД состояние, короткий пунктир - экспериментальные значения, сплошная линия - НМ состояние
смещается лишь на 4 %, другой апикальный атом - на 6,5 %, тогда как планарные атомы - на 7,1 %. Полученные результаты находятся в согласии с аналогичными результатами, полученными Н. Медведевой в работах [18, 23], где при растворении углерода в НМ состоянии было получено, что атом С увеличивает параметр решётки с 3,456 А до 3,472 А и расталкивает ближайшие атомы железа на 7,7 %, а при растворении углерода в АФМД состоянии с 3,53 А до 3,55 А.
Реальное ГЦК-железо существует при температурах выше 911 °С. При такой температуре атомы углерода имеют практически одинаковую вероятность нахождения во всех конфигурационных сферах. Поэтому для вычисления параметра решётки было взято среднее значение параметра решётки, которое составило для НМ <а> = 3,5 А и для АФМД <а> = 3,575 А. Экспериментальная зависимость параметра решётки от концентрации и температуры была получена в работе [27] нейтронной дифракцией.
Проведя сравнение наших результатов с экстраполяцией экспериментальных результатов (рис. 1), мы пришли к выводу, что АФМД-состояние намного лучше описывает железо, чем НМ. и для него наблюдается довольно хорошее согласие углового коэффициента зависимости параметра решётки от концентрации углерода.
Энергии растворения атомов углерода в ГЦК железе
При использовании формулы (1), нами была вычислена энергия растворения атома углерода. В табл. 1 приведено сравнение нашей энергии с данными других работ, среди которых и экспериментальные данные для парамагнитного железа [28, 29].
Таблица 1
Значения энергии растворения одного атома углерода в октапоре для немагнитного состояния
(Результаты представлены для концентрации 3 ат. % углерода)
Статья Состояние Энергия, мэВ
[28] эксп. (парамагнитное) 360
[29] эксп. (парамагнитное) 420
[28] НМ 145
Данная работа НМ 180
[17] АФМД 550
Данная работа АФМД 340
Видно, что энергия, полученная при растворении углерода в решётке ГЦК-железа, находящегося в АФМД-состоянии, практически совпадает с экспериментальным значением.
Наиболее полный анализ растворения второго атома углерода проведён в работе [18], однако он был проведен только для НМ-состояния матрицы. В суперячейке из 32 атомов железа существует 4 координационных сферы (рис. 2). Конфигурация с двумя ближайшими атомами углерода, образующими 180°
С-Бе-С позицию, является наименее энергетически выгодной.
Наши данные для НМ-состояния хорошо согласуются с данными Медведевой [18].
В табл. 2 приведено сравнение значений энергии растворения 2 атомов углерода, полученных в нашем расчете как для НМ, так и для АФМД-фазы ГЦК-железа. Для их получения использовалась формула (6).
Таблица 2
Значения энергии растворения второго атома углерода_______________________________
Номер конфигурации НМ, мэВ АФМД, мэВ
1 260 510
2 1400 840
3 150 470
4 310 250
Проведя анализ полученных результатов, был сделан вывод, что энергия растворения второго атома углерода близка к энергии растворения первого для всех конфигураций, кроме второй.
Рис. 2. Возможные расположения атомов углерода в ячейке Ре32
Ab-initio моделирование влияния ближнего окружения примесей углерода на энергию их растворения в ГЦК-железе
Для оценки упругих взаимодействий в рассматриваемой системе мы нашли разность полной энергии релаксированной системы и энергии системы, в которой не производилось релаксации. Назовем эту величину энергией упругой релаксации. Данные представлены в табл. 3.
Для атомов углерода, находящихся во второй координационной сфере, релаксация оказалась наиболее слабой, потому что атомы углерода находятся на одной линии и разделены одним атомом Бе, который по этой причине не двигается. Скорее всего, это обстоятельство и объясняет энергетическую невыгодность данной конфигурации.
Таблица 3
Значения энергий упругой релаксации (Erelax - Enereiax)__________________________
Номер конфигурации НМ, мэВ АФМД, мэВ
1 -4220 -2600
2 -3180 -2140
3 -4000 -2460
4 -3860 -2600
Для вычисления энергии взаимодействия между атомами была использована формула (7). В табл. 4 мы привели значения энергий растворения второго атома углерода для различных конфи-
гураций в суперячейки состоящей из 32 атомов железа.
Таблица 4
________________________Значения энергий взаимодействия между атомами углерода_________________________
Номер конфигурации ГаврилюкВ.Г. [13], мэВ Данная работа, НМ, мэВ Данная работа, АФМД, мэВ
1 150 80 170
2 400 1220 500
3 - -30 130
4 - 130 -90
Видно, что учет магнитных эффектов (АФМД-состояние матрицы) несколько уменьшает отталкивание атомов углерода во второй координационной сфере, но увеличивает его в первой и третьей сферах. Атомы углерода в ГЦК-железе отталкиваются и предпочитают занимать наиболее удалённые позиции. Полученные нами значения энергии взаимодействия между атомами углерода в АФМД-состоянии матрицы хорошо согласуются с работой Гаврилюка В.Г. [13], в которой он получил указанные энергии, используя метод функции Грина. Это говорит о том, что АФМД-приближение для описания магнитного состояния может быть использовано для расчета как структуры реального ГЦК-железа, так и энергии взаимодействия атомов углерода в нем.
Заключение
Таким образом, в результате проведенных исследований:
1. Были найдены: энергия растворения приходящаяся на один атом углерода (при концентрации 3 ат.%) и энергии растворения второго атома углерода (при концентрации 6 ат.%, в зависимости от их взаимного расположения) в суперячейке из 32 атомов ГЦК решётки железа.
2. Растворение одиночного атома С в НМ ГЦК-железе увеличивает параметр решётки с 3,45 А до 3,475 А и расталкивает ближайшие атомы железа на 7,5 %. При растворении углерода в АФМД-состоянии происходит увеличение параметра решётки с 3,54 А до 3,555 А и образование в первом окружении слегка несимметричного октаэдра.
3. Показано, что между соседними атомами углерода в ГЦК-железе наблюдается отталкивание при их расположении в первой, второй и третьей координационных сферах. При этом наиболее сильно отталкивание между вторыми ближайшими соседями.
4. Показано, что для моделирования свойств ГЦК-железа с примесями углерода лучше использовать АФМД-состояние матрицы, чем НМ-состояние.
Литература
1. Кристаллография, рентгенография и электронная микроскопия / Я.С. Уманский, Ю.А. Скаков, А.Н. Иванов, Л.Н. Расторгуев. - М.: Металлургия. - 1982. - 632 с.
2. Курдюмов, Г.В. Превращения в железе и стали / Г.В. Курдюмов, Л.М. Утевский, Р.И. Энтин. - М.: Наука. - 1977. - С. 236.
3. Хачатурян, А.Г. Углерод в мартенсите стали. Несовершенства кристаллического строения и мартенситные превращения / А.Г. Хачатурян. - М.: Наука, 1971. - C. 34-45.
4. Могутнов, Б.М. Термодинамика железо-углеродистых сплавов / Б.М. Могутнов, Н.А. Томилин, Л.А. Шварцман. - М.: Металлургия, 1972. - C. 328.
5. Dunwald, H. Thermodynamische Untersuchungen zum System Eisen-Kohlenstoff-Sauerstoff /
H. Dunwald, C. Wagner // Z. Anorg. Allgem Chem. - 1931. - Vol. 199. - P. 321-346.
6. Darken, L.S. Appendix to the paper by Smith R.P. Equilibrium of iron-carbon alloys / L.S. Darken, R.P. Smith // J.Amer.Chem. Soc. - 1946. - № 7. - P. 1163-1175.
7. Темкин, М.И. Активность углерода в аустените / М.И. Темкин, Л.А. Шварцман // Журнал физической химии. - 1949. - № 6. - С. 755-760.
8. Kaufman, L. Thermodynamics of bainite reaction / L. Kaufman, S.V. Radcliffe, M. Cohen / Decomposition of Austenite by Diffusional Processes: сб. науч. тр. - AIME, New York: Interscience Publishers, 1962. - P. 313-352.
9. McLellan, R.B. A quassi-chemical treatment of interstitial solid solutions: it application to carbon austenite / R.B. McLellan, W.W. Dunn // J.Phys.Chem.Solids. - 1969. - V. 30, № 11. - P. 2631-2637.
10. Кожеуров, В.А. К термодинамике растворов внедрения / В.А. Кожеуров // Известия ВУЗов. Черная металлургия. - 1965. - № 2. - С. 10-16.
11. Dunn, W.W. The Application of Quassi-chemical Solid Solution Model to Carbon Austenite / W.W. Dunn, R.B. McLellan // Metall Trans. - 1970. - Vol. 1, № 5. - P. 1263-1265.
12. Большов, Л.А. Изучение ближнего порядка в аустентите системы Fe-C методами статистической теории и Монте-Карло / Л.А. Большов, В.Н. Суслов. // Физика металлов и металловедение. - 2004. - Т. 98, № 6. - P. 3-7.
13. Гаврилюк, В.Г. Распределение углерода в стали / В.Г. Гаврилюк. - Киев: Наукова Думка, 1987.- C. 208.
14. Надутов, В.М. Межатомное взаимодействие и распределение атомов внедрения в железоазотистых и железо-углеродистых сплавах: автореферат дис. д-ра физ.-мат. наук / В.М. Надутов. - Киев: ИМФ НАН Украины им. Г.В.Курдюмова, 1997. - 55 с.
15. Jiang, D.E. Carbon dissolution and diffusion in ferrite and austenite from first principles /
D.E. Jiang, E.A. Carter // Physical Review B. - 2003. - Vol. 67. - P. 214103.
16. Slane, J.A. Experimental and Theoretical Evidence for Carbon-Vacancy Binding in Austenite / J.A. Slane, C. Wolverton, R. Gibala // Metallurgical and Materials Transactions A. - 2004. - Vol. 35, № 8. - P. 2239-2245.
17. Magnetism and Local Distortions near Carbon Impurity in g-Iron / D.W. Boukhvalov,
Y.N. Gornostyrev, M.I. Katsnelson, A.I. Lichtenstein. // Physical Review Letters. - 2007. - Vol. 99. -
P.247205.
18. Ивановский Л.И. Распределение примесей Mn, Al, Si, С в аустените и их влияние на энергии дефектов упаковки / Л.И. Ивановский, Н.И. Медведева. // Фазовые переходы, упорядоченные состояния и новые материалы. - 2012. - № 10. - P. 24-28.
19. Timoshevskii, A.N. Ab-initio modeling of the short range order in Fe-N and Fe-C austenitic alloys / A.N. Timoshevskii, S.O. Yablonovskii // Functional Materials. - 2011. - Vol. 18, № 4. - P. 517522.
20. Gavriljuk, V.G. Change in the electron structure caused by C, N and H atoms in iron and its effect on their interaction with dislocations / V.G. Gavriljuk, V.N. Shivanyuk, B.D. Shanina // Acta Mate-rialia. - 2005. - Vol. 53. - P. 5017-5024.
21. Schwarz, K. Electronic structure calculations of solids using the WIEN2k package for material science / K. Schwarz, P. Blaha, G.K.H. Madsen // Computer Physics Communications. - 2002. -Vol. 147. - P. 71-76.
22. The Role of the Nature of Magnetic Coupling on the Martensitic Transformation in Fe-Ni / M. Acet, E.F. Wassermann, K. Andersen et al. // Journal de Physique IV France. - 1997. - Vol. 7, № C5. - P. 401-404.
23. Medvedeva, N.I. Magnetism in bcc and fcc Fe with carbon and manganese / N.I. Medvedeva, D.V. Aken, J.E. Medvedeva // Journal of Physics: Condensed Matter. - 2010. - Vol. 22. - P. 316002.
24. Cottenier, S. Density Functional Theory and the family of (L)APW-methods: a step-by-step introduction/ S. Cottenier, 2004. [http://www.wien2k.at/reg_user/textbooks/DFT_and_LAPW-2_cottenier.pdf]
Ab-initio моделирование влияния ближнего окружения примесей углерода на энергию их растворения в ГЦК-железе
25. Урсаева, A.B. Выбор оптимальных параметров для построения максимально точной модели ОЦК-железа / A.B. Урсаева, Г.Е. Рузанова, A.A. Мирзоев. // Вестник ЮУрГУ. Серия «Математика. Механика. Физика». - 2010. - Вып. 2. - № 9(185). - С. 97-101.
26. Ракитин, M.C. Изменение электронной структуры a-железа, содержащего внедренные атомы водорода / M.C. Ракитин, A.A. Мирзоев, Д.А. Мирзаев // Вестник ЮУрГУ. Серия «Мета-лургия». - 2010. - Вып. 14. - № 13(189). - P. 67-71.
27. The lattice parameters of austenite and ferrite in Fe-C as functions of carbon concentration and temperature / M. Onink, C.M. Brakman, F.D. Tichelaar et al. // Scripta Metallurgica Et Materialia.-1993. - Vol. 29, № 8. - P. 1011-1016.
28. Slane, J.A. Experimental and Theoretical Evidence for Carbon-Vacancy Binding in Austenite / J.A. Slane, C. Wolverton, R. Gibala // Metallurgical and Materials Transactions A. - 2004. - Vol. 35, № 8. - P.2239-2245.
29. Lobo, J.A. Thermodynamics of carbon in austenite and Fe-Mo austenite / J.A. Lobo, G.H. Geiger // Metallurgical Transactions A. - 1976. - Vol. 7, № 8. - P. 1359-1364.
AB-INITIO SIMULATION OF INFLUENCE OF SHORT-RANGE ORDERING CARBON IMPURITIES ON THE ENERGY OF THEIR DISSOLUTION IN THE FCC-IRON
Ya.M. Ridnyi1, A.A. Mirzoev2, D.A. Mirzaevz
The first principle simulation of equilibrium structure and FCC-iron characteristics was carried out by the software package WIEN2k. The optimal parameters which allow building the most accurate model were generated. Energies of dissolution of carbon atoms, their relative positions and the contribution of the elastic effects to the energy of system were calculated for non-magnetic (NM) and doublelayer antiferromagnetic states (AFMD) of FCC-iron.
Keywords: first principle simulation, FCC-iron, carbon impurity, WIEN2k.
References
1. Umanskiy Ya.S., Skakov Yu.A., Ivanov A.N., Rastorguev L.N. Kristallografiya, rentgenografiya i elektronnaya mikroskopiya (Crystallography, radiographic imaging and electronic microscopy). Moscow: Metallurgiya. 1982. 632 p. (in Russ.).
2. Kurdyumov G.V., Utevskiy L.M., Entin R.I. Prevrashcheniya v zheleze i stali (Transformations in iron and steel). Moscow: Nauka. 1977. 236 p. (in Russ.).
3. Khachaturyan A.G. Uglerod v martensite stali. Nesovershenstva kristallicheskogo stroeniya i martensitnye prevrashcheniya (Carbon in martensite of steel. Imperfections of crystalline construction and martensitic transformations). Moscow: Nauka. 1971. pp. 34-45. (in Russ.).
4. Mogutnov B.M., Tomilin N.A., Shvartsman L.A. Termodinamika zhelezo-uglerodistykh splavov (Thermodynamics of iron-carbon alloys). Moscow: Metallurgiya. 1972. 328 p. (in Russ.).
5. Dünwald H., Wagner C. Thermodynamische Untersuchungen zum System Eisen-KohlenstoffSauerstoff. Z. Anorg. Allgem Chem. 1931. Vol. 199. pp. 321-346. DOI: 10.1002/zaac.19311990132
6. Darken L.S., Smith R.P. Appendix to the paper by Smith R.P. Equilibrium of iron - carbon alloys / // J.Amer.Chem. Soc. 1946. no. 7. pp. 1163-1175.
7. Temkin M.I., Shvartsman L.A. Zhurnal fizicheskoy khimii. 1949. no. 6. pp. 755-760. (in Russ.).
8. Kaufman L.R., Radcliffe S.V., Cohen M. Thermodynamics of the Bainite Reaction. Decomposition of Austenite by Diffusional Processes. AIME, New York: Interscience Publishers, 1962. pp. 313352.
1 Ridnyi Yaroslav Maksimovich is Student, South Ural State University.
E-mail: yaros lav@physics.susu.ac.ru
2 Mirzoev Aleksander Aminulaevich is Dr. Sc. (Physics and Mathematics), Professor, General and Theoretical Physics Department, South Ural State University.
E-mail: mirzoev@physics.susu.ac.ru
3 Mirzaev Dzhalal Aminulovich is Dr. Sc.(Physics and Mathematics), Professor, Department of Physical Metallurgy and Material Science, Southern Ural State University.
E-mail: mirzayev@physmet.susu.ac.ru
9. McLellan R.B., Dunn W.W. A quassi-chemical treatment of interstitial solid solutions: it application to carbon austenite. J.Phys.Chem.Solids. 1969. Vol. 30. no. 11. pp. 2631-2637.
10. Kozheurov V.A. Izvestiya VUZov. Chernaya metallurgiya. 1965. no. 2. pp. 10-16. (in Russ.).
11. Dunn W.W., McLellan R.B. The Application of Quassi-chemical Solid Solution Model to Carbon Austenite. Metall Trans. 1970. Vol. 1, no. 5. pp. 1263-1265.
12. Bol'shov L.A., Suslov V.N. Fizika metallov i metallovedenie. 2004. Vol. 98, no. 6. pp. 3-7. (in Russ).
13. Gavrilyuk V.G. Raspredelenie ugleroda v stali (Carbon distribution in steel). Kiev: Naukova Dumka. 1987. p. 208.
14. Nadutov V.M. Mezhatomnoe vzaimodeystvie i raspredelenie atomov vnedreniya v zhelezo-azotistykh i zhelezo-uglerodistykh splavakh. Avtoreferat na soiskanie uchyenoy stepeni doktora fiziko-matematicheskikh nauk (Interatomic interaction and distribution of interstitial atoms in iron-nitrogenuous and iron-carbon alloys. Dr. phys. and math. sci. synopsis of diss.). Kiev: IMF NAN Ukrainy im. G.V. Kurdyumova, 1997. 55 p.
15. Jiang D.E., Carter E.A. Carbon dissolution and diffusion in ferrite and austenite from first principles. Physical Review B. 2003. Vol. 67. pp. 214103.
16. Slane J.A., Wolverton C., Gibala R. Experimental and Theoretical Evidence for Carbon-Vacancy Binding in Austenite. Metallurgical and Materials Transactions A. 2004. Vol. 35, no. 8. pp. 2239-2245.
17. Boukhvalov D.W., Gornostyrev Y.N., Katsnelson M.I., Lichtenstein A.I. Magnetism and Local Distortions near Carbon Impurity in g-Iron. Physical Review Letters. 2007. Vol. 99. p. 247205.
18. Ivanovskiy L.I., Medvedeva N.I. Fazovyeperekhody, uporyadochennye sostoyaniya i novye ma-terialy. 2012. no. 10. pp. 24-28. (in Russ.).
19. Timoshevskii A.N., Yablonovskii S.O. Functional Materials. 2011. Vol. 18, no. 4. pp. 517-522.
20. Gavriljuk V.G., Shivanyuk V.N., Shanina B.D. ActaMaterialia. 2005. Vol. 53. pp. 5017-5024.
21. Schwarz K., Blaha P., Madsen G.K.H. Computer Physics Communications. 2002. Vol. 147. pp. 71-76.
22. Acet M., Wassermann E.F., Andersen K., Murani A., Scharpff O. Journal de Physique IV France. 1997. Vol. 7, no. C5. p. 401-404.
23. Medvedeva N.I., Aken D.V., Medvedeva J.E. Journal of Physics: Condensed Matter. 2010. Vol. 22. pp. 316002.
24. Cottenier S. Density Functional Theory and the family of (L)APW-methods: a step-by-step introduction. 2004. Available at: http://www.wien2k.at/reg_user/textbooks/DFT_and_LAPW-
2_cottenier.pdf
25. Ursaeva A.V., Ruzanova G.E., Mirzoev A.A. Vybor optimal'nykh parametrov dlya postroeniya maksimal'no tochnoy modeli OTsK-zheleza (Selection of optimal parameters for formation the most accurate model of BCC iron). Vestnik YuUrGU. Seriya «Matematika. Mekhanika. Fizika». 2010. Issue 2. no. 9(185). pp. 97-101. (in Russ.).
26. Rakitin M.C., Mirzoev A.A., Mirzaev D.A. Izmenenie elektronnoy struktury a-zheleza, soderz-hashchego vnedrennye atomy vodoroda (Change of electronic structure in iron containing interstitial atoms of hydrogen) // Vestnik YuUrGU. Seriya: Metalurgiya. 2010. Issue 14. no. 13(189). pp. 67-71. (in Russ).
27. Onink M., Brakman C.M., Tichelaar F.D., Mittemeijer E.J., Van der Zwaag S., Root J.H., Kon-yer N.B. ScriptaMetallurgica EtMaterialia. 1993. Vol. 29, no. 8. pp. 1011-1016.
28. Slane J.A., Wolverton C., Gibala R. Metallurgical and Materials Transactions A. 2004. Vol. 35, no. 8. pp. 2239-2245.
29. Lobo J.A., Geiger G.H. Metallurgical Transactions A. 1976. Vol. 7, no. 8. pp. 1359-1364.
Поступила в редакцию 11 июня 2013 г.
