CC BY
53
УДК 669.111.31
АНАЛИЗ ВОЗМОЖНОСТИ ПЕРЕРАСПРЕДЕЛЕНИЯ АТОМОВ УГЛЕРОДА В РЕШЁТКЕ ЦЕМЕНТИТА
К.Ю. Окишев
ANALYSIS OF POSSIBILITY OF CARBON REDISTRIBUTION IN CEMENTITE LATTICE
K.Yu. Okishev
На основании результатов существующего первопринципного расчёта энергии атома углерода в «нормальной» октаэдрической поре решётки цементита (по отношению к «нормальной» призматической) оценены его энергии в «искажённых» призматической и октаэдрической поре. Построена теория, позволяющая находить концентрацию атомов углерода в порах различных типов. Показано, что предположение о переходе части атомов углерода в «искажённые» поры не может объяснить наблюдаемых изменений EELFS-спектров цементита при длительном отжиге перлита ниже Ac1.
Ключевые слова: цементит, углерод, кристаллическая решётка, поры решётки, перераспределение атомов.
On the basis of existing ab initio calculation of energy of a carbon atom in the “normal” octahedral pore of cementite lattice (related to “normal” prismatic pore) its energy in the “distorted” prismatic and octahedral pores is estimated. A theory is developed that permits to find concentration of carbon atoms in various types of pores. It is shown that supposition that some of carbon atoms are transferred to “distorted” pores cannot explain the observed changes in EELFS spectra of cementite after long annealing of pearlite below Ac1.
Keywords: cementite, carbon, crystal lattice, lattice pores, atomic redistribution.
Цементит Бе3С является одной из основных фаз углеродистых сталей и чугунов. Он имеет ор-торомбическую решётку с параметрами a = 4,524; Ь = 5,088; c = 6,742 А [1] и принадлежит к пространственной группе Pnma. В элементарную ячейку входят 12 атомов железа и 4 атома углерода. Положение атомов железа в решётке установлено надёжно, однако для атомов углерода рассматриваются два варианта размещения: в октаэдрических либо призматических порах железной подрешётки [2].
В последние годы появился ряд работ, в частности, Ю.В. Раца с сотрудниками [3, 4], в которых обнаружено заметное изменение ближайшего ок-
ружения атомов в цементите при термообработке. На рис. 1 представлена функция радиального распределения атомов железа в ближайшем окружении атомов углерода (парная корреляционная функция - ПКФ С-Бе) в цементите стали У15, полученная в [4] на основе анализа спектров энергетических потерь электронов (метод ББЬБ8) после двух вариантов термической обработки: 1) аусте-нитизация при 1050 °С и последующий перлитный распад при 500 °С в течение 1 мин (режим 1); 2) то же плюс дополнительный отжиг при 700 °С в течение 20 ч (режим 2). В обоих случаях наблюдаются два максимума, соответствующие ~2,0 и
ос
1,50 1,75 2,00 2,25 2,50 2,75
R, А
Рис. 1. Парная корреляционная функция О-Ре, полученная из анализа ЕЕЬРБ-спектров цементита в стали У15 после термообработки по режимам 1 и 2 по данным [4]
2,3...2,4 А, однако после дополнительного отжига высота второго максимума сильно увеличивается.
В работах [5, 6] нами было показано, что в железной подрешётке цементита можно выделить поры четырёх типов: помимо упомянутых «нормальных» октаэдрических и призматических пор имеются также «искажённые» октаэдрические и призматические поры, и было высказано предположение, что наблюдаемые изменения ближайшего окружения атомов обусловлены переходом атомов углерода из «нормальных» призматических пор в поры иных типов. Целью данной работы является анализ возможности такого перераспределения атомов углерода в цементите.
В принципе, любые варианты расположения атомов углерода вне «нормальных» призматических пор должны повышать энергию цементита, так как радиус последних составляет 1,99.2,04 А, что соответствует сумме атомных радиусов углерода и железа (0,77 + 1,24 = 2,01 А [7]), а радиус пор прочих типов заметно меньше (см. таблицу). Однако перераспределение атомов углерода должно повышать и энтропию системы, поэтому состоянию равновесия (минимуму свободной энергии) должно соответствовать наличие какого-то количества атомов углерода в порах всех трёх типов. Поскольку длительный отжиг при 700 °С должен приводить систему в более равновесное состояние, то можно предположить, что именно такими переходами и объясняется наблюдаемое изменение ПКФ С-Бе.
Итак, в структуре цементита имеются четыре типа пор, причём количество пор каждого типа одинаково и втрое меньше числа атомов железа. Предположим, что в кристалле цементита, состоящем из 3N атомов железа и N атомов углерода, п атомов углерода перешло из «нормальных» призматических пор, где они первоначально располагались, в «нормальные» октаэдрические; п2 атомов углерода - в «искажённые» призматические и п3 атомов - в «искажённые» октаэдрические, тогда количество атомов углерода, оставшихся в «нормальных» призматических порах, равно
п0 = N - п1 - п2 - п3. (1)
Примем энергию исходного состояния за ноль (е0 = 0) и обозначим энергию остальных трёх состояний атомов углерода по отношению к исходному £1, е2 и е3; будем приближённо считать, что взаимодействия между атомами углерода нет, то есть эти энергии остаются постоянными для любой конфигурации углеродной подрешётки. Тогда изменение энергии при переходе атомов составит ЛE = п1£1 + п2£2 + п3£3 . (2)
Происходящее при этом изменение энтропии найдём с помощью формулы Больцмана
Л£ = k 1п W , (3)
где № - термодинамическая вероятность данного состояния системы, а к - постоянная Больцмана. Если ограничиться рассмотрением конфигурационной составляющей энтропии, то № равно числу способов размещения п0 атомов углерода по N «нормальным» призматическим порам, п атомов углерода - по N «нормальным» октаэдрическим порам и т. д., а именно:
N! N!
W = -
(N-n0)!n0! (N-n1)!n1! N! N!
. (4)
(N-n2)!n2! (N - n3)!n3!
Изменение свободной энергии системы при
перераспределении атомов углерода составит
AF = AE - TAS = AE - kT ln W . (5)
Равновесию отвечает равенство нулю всех трёх
частных производных свободной энергии:
dAF dAF dAF Л
-----=------=------= 0 . (6)
dn1 dn2 dn3
Эти производные нетрудно найти, если учесть соотношение (1) и избавиться от факториалов в (4) с помощью формулы Стирлинга ln M! ~ M ln M - M. Так, например,
dAF
dn1
= е1 - kT ln
'0
N - n.
= 0 .
(7)
0 y
Если обозначить доли атомов углерода, занимающих соответствующие позиции, X- = Пг^ (г = 0, 1, 2, 3) и ввести новые переменные
Атомы железа в ближайшем окружении (до 3 А) атомов углерода, находящихся в порах различного типа [5, 6]
Тип поры
Количество атомов и расстояния до них в Á
Тип поры
Количество атомов
Á
2 1,804
2 1,858
2 1,870
2 1,247
2 2,270
2 2,279
2 2,981
«Нормальная» призматическая (z = 6)
1,989
1,992
2,018
2,039
2,390
2,836
«Нормальная» октаэдрическая (z = 6)
«Искажённая» призматическая (Z = 4)
1,622
1,624
1,626
2,516
«Искажённая» октаэдрическая (z = 2)
Уг = (1 - Хг )/хг , то окончательно получаем следующую систему уравнений, описывающих равновесное распределение атомов углерода по порам всех четырёх типов:
' У1 = У0/ 4;
' у2 = у0 /В2 ; (8)
„ у3 = у0 / В3 ,
где Вг = ехр(-£г /кТ). В таком случае условие х0 + х1 + х2 + х = 1 можно переписать в виде
у0 - у0 (В1 + В2 + В3 + В1В2 + В2В3 + В1В3 ) -
-2у0 (В1В2 + В2В3 + В1В3 + В1В2В3 ) -
-3В1В2В3 = 0 . (9)
Решение этого уравнения позволяет найти величину у0, а затем при помощи соотношений (8) и остальные уг и концентрации хг.
Расчёт концентраций можно упростить, если учесть, что энергии £г вряд ли могут быть слишком малы, а рассматриваемые температуры не превышают 1000 К. По этой причине коэффициенты Вг должны быть существенно меньше единицы и их двойными и тройными произведениями в (9) можно пренебречь. Тогда равновесные значения
уг “V В1 + В2 + В3 / Вг , то есть
В (10)
Вг + ^ В1 + В2 + В3
Поскольку сам вопрос о том, в каких позициях железной подрешётки цементита размещаются атомы углерода до сих пор не имеет окончательного ответа, то прямых экспериментальных измерений энергии соответствующих состояний, разумеется, нет. Однако в последние годы появились работы, где энергии кристаллической решётки цементита были рассчитаны «из первых принципов». В первой работе такого рода [8] было найдено, что при размещении всех атомов углерода в «нормальных» призматических позициях энергия когезии (энергия, выделяющаяся при образовании кристаллической решётки из изолированных атомов, взятая с обратным знаком) цементита Бе3С составляет 7,33 эВ/атом, а при их размеще-
нии в «нормальных» октаэдрических позициях -7,16 эВ/атом. Следовательно, если принять за нулевой уровень энергию атома углерода в «нормальной» призматической поре (е0 = 0), то его энергия в «нормальной» октаэдрической будет составлять е1 =-(7,16 - 7,33)/хСе3С = 0,68 эВ; здесь
х,
£е3С = 0,25 - атомная доля углерода в цементите.
В более поздней работе [9] была получена разность энергий решётки цементита с размещением атомов углерода в «нормальных» октаэдрических и призматических порах 0,068 эВ/атом, откуда
81 =- 0,068/ хСе3С = 0,272 эВ.
Расчёты энергии атома углерода в «искажённых» призматических и октаэдрических порах в [8, 9] не делались, поэтому попытаемся получить приближённую оценку на основании имеющихся результатов. Будем считать, что энергия атома углерода пропорциональна числу атомов железа в первой координационной сфере, а зависимость этой энергии от радиуса первой координационной сферы в неискажённой решётке линейна. Согласно таблице для «нормальных» призматических и октаэдрических пор число таких атомов железа 2 = 6, а расстояния до них составляют Яг = 1,99.2,04 и 1,80.1,87 А
соответственно. Среднее расстояние составляет для «нормальных» призматических позиций — 1 2 1
Я0 =_. У Я1 = --12,066 = 2,011 А, (11)
2 г=1 6
а для «нормальных» октаэдрических -
— 1 2 1
Я1 =--У Я =— 11,064 = 1,844 А. (12)
2 г =1 6
«Искажённую» призматическую пору окружают г = 4 атома железа на расстояниях Яг = 1,62 .1,63 А, поэтому для неё
Я2 =
= _. 6,494 = 1,624
(13)
г=1
Наконец, в первой координационной сфере «искажённой» октаэдрической поры располагаются всего 2 = 2 атома железа на очень малом расстоянии Яг = 1,247 А, так что
со
СО N
СО
R = Ж / г, А
Рис. 2. Зависимость энергии атома углерода в решётке цементита £, отнесённой к числу атомов железа в первой координационной сфере z, от радиуса первой координационной сферы Я
— 1 2 1
Я3 = - • У Я = - - 2,494 = 1,247 А.
(14)
г=1
На рис. 2 представлена зависимость энергии этих состояний в расчёте на один атом железа в первой координационной сфере от среднего радиуса этой сферы Я . Прямая, проходящая через значения для «нормальных» призматических и октаэдрических позиций, описывается уравнением
е/2 = 0,5459 - 0,2715 - Я , эВ. (15)
Это даёт для энергии атома углерода в «искажённой» призматической поре оценку е2/ 2 = 0,105 эВ, откуда
е2 = 0,421 эВ, (16)
а в «искажённой» октаэдрической е3 /2 = 0,207 эВ и е3 = 0,415 эВ. (17)
Отметим, что значения энергий для «искажённых» пор обоих типов оказались очень близки и не слишком значительно превышают значение для «нормальных» октаэдрических пор.
При указанных значениях энергий коэффициенты Вг при 700 °С будут равны В1 = 0,0390 ; В2 = 0,0066; В3 = 0,0071, а концентрации
х1 = 0,145; х2 = 0,028; х3 = 0,030. (18)
Очевидно, что столь малые количества атомов углерода в порах разных типов не могут вызвать
наблюдавшиеся изменения ББЬБ8-спектров (см. рис. 1). Конечно, приведённые оценки являются очень приблизительными, так как действительные значения энергий е2 и е3 могут быть иными.
Поэтому будем считать эти энергии свободными параметрами и попытаемся оценить, как их значения будут влиять на степень заполнения пор разных типов. На рис. 3 приведены значения концентраций хг, рассчитанных по формулам (9) и (8) при условии, что е1 = 0,272 эВ и е2 =е3. Как видно, доля атомов углерода, находящихся в «искажённых» порах, становится заметной при значениях е2 и е3 менее ~0,3 эВ.
Вернёмся к результатам [4] (см. рис. 1). Фактически после обоих режимов обработки у ПКФ наблюдались два пика - при ~2,0 и 2,3.2,4 А. Как видно из таблицы, атомы железа могут находиться на таком расстоянии от атома углерода при различных вариантах расположения последнего. Поскольку пиков на расстояниях меньших чем 2 А не наблюдается, то можно предположить, что даже при размещении атома углерода в «тесных» порах он смещает окружающие атомы железа, так что они оказываются от него на расстоянии, близком к 2 А (которое, как было отмечено выше, равно сумме атомных радиусов углерода и железа). В работе [10] показано, что как бы ни было мало расстояние между атомами в
1,0
0,8
0,6
0,4
0,2
. • "1 * 1 1 - : "Искажённые" ; октаэдрические : поры 1 1 ' 1 1 1 • . ’ .....^"Нормальные" окт. поры;
/ "Нормальные" : / призматические поры " *
Г "Искажённые" ; призматические / ~ поры / : . л .
х го ч о
<и
ЕЕ
ш о 2 О Iго
к
о 0,0 с! -0,6 -0,4
-0,2 0,0 0,2 0,4 0,6 0,8 1,0
е2 = £3- 36
Рис. 3. Зависимость доли атомов углерода в порах каждого из четырёх типов от избыточной энергии атомов углерода в «искажённых» призматических и октаэдрических порах е2 = е3
ез, эВ
Рис. 4. Зависимость отношения высот первого и второго максимума ПКФ О-Ре от избыточной энергии атомов углерода в «искажённых» призматических и октаэдрических порах £2 = £3
е
2
решётке железа при помещении между ними атома углерода расстояние С-Бе всегда оказывается не менее 1,75.1,90 А.
Если считать, что смещаются атомы только первой координационной сферы, то на основании данных таблицы можно записать, что отношение интенсивностей первого и второго пиков ПКФ
6 х0 + 6 х, + 4 х2 + 2 х3 3 - х2 - 2 х3
а = —0-----1----2----3 =-----2----3. (19)
2 х0 + 2 х2 + 4 х 1 - х1 + х3
На рис. 4 приведены значения этого отношения, рассчитанные на основании данных рис. 3. Из рис. 4 следует, что при любых значениях
е2 = е3 > 0 первый пик должен быть минимум в
полтора раза выше второго.
В то же время можно заметить, что некоторые первопринципные расчёты [11] указывают на весьма малые значения энергии образования вакансии в углеродной подрешётке цементита (~0,2 эВ). Это указывает, во-первых, на то, что и другие искажения углеродной подрешётки могут приводить лишь к небольшим изменениям энергии кристалла Бе3С, а во-вторых, что при приближении к равновесию в процессе отжига цементит может терять часть атомов углерода из пор разных типов, что скажется на ближайшем окружении оставшихся атомов. Поэтому для ответа на вопрос о возможности возникновения разных конфигураций углеродной подрешётки требуются прямые расчёты энергий соответствующих состояний цементита.
Работа выполнена в рамках ФЦП «Научные и научно-педагогические кадры инновационной России».
Литература
1. Эндрюс, К. Электронограммы и их интерпретация / К. Эндрюс, Д. Дайсон, С. Киоун. - М.: Мир, 1971. - 256 с.
2. Гольдшмидт, Х.Дж. Сплавы внедрения / Х. Дж. Гольдшмидт. - М.: Мир, 1971. - Т. 1. - 424 с.
3. Влияние термической обработки на локальную атомную структуру цементита в стали / А.Н. Маратканова, Ю.В. Рац, Д,В. Сурнин и др. // ФММ. - 2000. - Т. 89, № 6. - С. 76-81.
4. Маратканова, А.Н. Исследование локальной атомной структуры цементита / А.Н. Маратканова, И.Л. Яковлева, Ю.В. Рац // ФММ. - 2004. -Т. 98, № 3. - С. 72-79.
5. О возможных позициях атомов углерода в решетке цементита / В.М. Счастливцев, И.Л. Яковлева, Д.А. Мирзаев, К.Ю. Окишев // ФММ. - 2003. -Т. 96, № 3. - С. 75-82.
6. Поры в кристаллической решётке цементита и положение атомов углерода /Д.А. Мирзаев, К.Ю. Окишев, В.М. Счастливцев, И.Л. Яковлева // Вестник ЮУрГУ. Серия «Математика, физика, химия». - 2003. - Вып. 3, № 6. - С. 79-85.
7. Григорович, В.К. Дисперсионное упрочнение тугоплавких металлов / В.К. Григорович, Е.Н. Шеф-тель. - М.: Наука, 1980. - 304 с.
8. Влияние эффектов атомного разупорядо-чения и нестехиометрии по углеродной подре-шётке на зонную структуру цементита Fe3C / Н.И. Медведева, Л.Е. Карькина, А.Л. Ивановский // ФММ. - 2003. - Т. 96, № 5. - С. 16-20.
9. Arzhnikov, A.K. Structural peculiarities of cementite and their influence on magnetic characteristics /A.K. Arzhnikov, L. V. Dobysheva, C. Demangeat // J. Phys.: Cond. Matter. - 2007. - Vol. 19. - 196214.
10. Jiang, D.E. Carbon Dissolution and Diffusion in Ferrite and Austenite from First Principles / D.E. Jiang, E.A. Carter // Phys. Rev. B. - 2003. -Vol. 67. - 214103.
11. Состав цементита в зависимости от температуры. In-situ нейтронография и результаты ab-initio расчётов / В.И. Воронин, И.Ф. Бергер, Ю.Н. Горностырёв и др. // Письма в ЖЭТФ. -2010. - Т. 91. - Вып. 3. - С. 154-157.
Поступила в редакцию 18 февраля 2011 г.
