Энергетические характеристики взаимодействия точечных дефектов с симметрично-наклонной границей ∑5(210)[001] в сплаве FeCr Текст научной статьи по специальности «Физика»

УДК 544.022.342, 544.022.344.2

ЭНЕРГЕТИЧЕСКИЕ ХАРАКТЕРИСТИКИ ВЗАИМОДЕЙСТВИЯ ТОЧЕЧНЫХ ДЕФЕКТОВ С СИММЕТРИЧНО-НАКЛОННОЙ ГРАНИЦЕЙ 25(210)[001] В СПЛАВЕ FeCr

© 2013 М.Ю. Тихончев, А.Б. Муралев, В.В. Светухин

Научно-исследовательский технологический институт им. С.П.Капицы Ульяновского государственного университета

Поступила в редакцию 26.11.2013

Работа посвящена изучению энергетических свойства симметрично-наклонной межзеренной границы 2 5(210)[001] в а -Ре, содержащем хром. Исследование проводится путем моделирование методом молекулярной статики на основе многочастичного межатомного взаимодействия. Получены оценки удельной энергии границы и размеров соответствующих им зернограничных областей в чистом железе при нулевой температуре и в бинарном сплаве Ре-9ат.%Сг при температурах 0 и 300 К. Рассчитаны энергии формирования и связи с зернограничной областью вакансии, собственных меж-доузленных атомов и замещающего атома Сг в железе. Рассмотрено два типа собственных междоу-зельных конфигураций: гантель с ориентацией <110> из двух атомов железа и гантель смешанного Ре-Сг типа той же ориентации. Установлено, что энергия связи междоузельных конфигураций с зернограничными областями положительна. Энергия связи вакансии с большинством позиций в зернограничных областях также положительна. Энергия связи замещающего атома меняет знак в зависимости от позиции атома Сг в межзеренной области, максимальная положительная энергия связи составляет 0.05 эВ.

Ключевые слова: энергетические свойства, симметрично-наклонная межзеренная граница, сплав Ре-Сг.

1. ВВЕДЕНИЕ

Межзеренные границы (ГЗ) традиционно рассматриваются как стоки для точечных дефектов и примесных атомов. Однако, детальные исследования механизмов взаимодействия ГЗ с указанными дефектами на атомарном уровне начали проводиться сравнительно недавно. Знание структуры и энергетических характеристик межзерен-ных границ (ГЗ - граница зерна) важно для понимания таких явлений, происходящих в материалах, как рост зерен, зернограничная диффузия, сегрегация примесей, деформация и разрушение (см. работу [1] и дальнейшие ссылки в ней).

В работе [1] Zhang et al., используя многочастичный потенциал межатомного взаимодействия на основе модифицированного аналитического метода погруженного атома [2], рассчитали удельные энергии для 27 различных симметрично наклонных ГЗ в a -Fe. Однако их оценки оказались несколько завышенными по сравнению с более поздними исследованиями (см., например,

Тихончев Михаил Юрьевич, кандидат физико-математических наук, начальник лаборатории моделирования поведения неорганических материалов. E-mail: [email protected]

Муралев Артем Борисович, младший научный сотрудник лаборатории тмпьютерного моделирования поведения неорганических материалов. E-mail: [email protected] Светухин Вячеслав Викторович, доктор физико-математических наук, профессор, директор научно-исследовательского технологического института. E-mail: [email protected]

работы [6,8,9]). В 1999-2001 годах Perez and Smith [3, 7] на основе потенциала Финниса-Синклера одними из первых провели молекулярно-дина-мическое (МД) моделирование каскадов атомных смещений возле СНГЗ в чистом железе. Была установлена склонность произведенных в каскаде точечных дефектов (вакансий и СМА - собственных междоузедльных атомов) аккумулироваться в тонком слое, прилегающем к ГЗ. Аналогичный вывод был сделан и в недавней работе [21]. Martino et al. [4] провели МД-моделирование каскадов атомных смещений возле симметрично наклонных ГЗ в железе, содержащем примесный фосфор. Моделирование проведено на основе многочастичного межатомного взаимодействия для энергий первично-выбитого атома (ПВА) 0.5 - 6 кэВ в интервале температур 400 - 700 °С. Было установлено, что ГЗ является барьером на пути движения и рекомбинации точечных дефектов, которые скапливаются и сохраняют стабильность вблизи ГЗ. Отмечено тенденция к миграции междоузельного фосфора в направлении к ГЗ. В работе [5] методом МД исследовалась диффузия одиночной вакансии в приграничной зоне ГЗ в a -Fe. В работе [11] Terentyev et al. применили метод Монте-Карло для исследования сегрегации Cr в сплавах Fe-xCr, x=5, 10, 14 ат. % в температурном интервале 300 - 900 К. В работе [12] проведено расчетное исследование особенностей прохождения каскадов через границу раздела фаз ГПУ Zr и ОЦК Nb и пороговых энергий смещения в межфазной области.

Настоящая работа посвящена исследованию ряда энергетических свойств симметрично наклонной ГЗ Е 5(210)[001] в а-Бе и сплаве Ре-9ат.%Сг. Исследование проводится путем моделирования методом молекулярной динамики. Выбор сплава обусловлен тем, что высокохромистые феррито-мартенситные сплавы рассматриваются в качестве кандидатных материалов первой стенки и элементов бланкета будущих реакторов синтеза [13-15].

2. МЕТОД МОДЕЛИРОВАНИЯ

Симметрично наклонная или двойниковая ГЗ представляет собой плоскость зеркальной симметрии между двумя монокристаллами и считается наиболее простым типом ГЗ. Такие ГЗ нарушают трансляционную инвариантность в направлении, перпендикулярном ее плоскости. В направлениях параллельных границе трансляционная симметрия сохраняется (рис. 1). В настоящем исследовании мы ограничились рассмотрением одного типа симметрично наклонных ГЗ для ОЦК структуры: Е 5(210) с осью вращения [001]. Модельный бикристаллит строили так, чтобы ось 2 была коллинеарна оси вращения монокристаллов, ось X - нормаль к СНГЗ, ось У - третья ось декартовой системы координат. На граничных плоскостях перпендикулярных ГЗ использовали периодические граничные условия. Координаты атомов, попадающих в слои шириной ~ 0.6 нм, прилегающие к границам параллельным ГЗ, фиксировали после построения бикристаллита и не меняли в ходе дальнейших расчетов. Т.о. на граничных плоскостях параллельных ГЗ использовались жесткие граничные условия (см. рис. 2).

Рассматривали два типа кристаллита: чистое железо и случайный бинарный сплав Бе-9ат.%Сг. Для железа использовали многочастичный потенциал из работы [16]. Для описания межатомного взаимодействия в системе БеСг ис-

пользовали модифицированную версию многочастичного потенциала, предложенного А. Каро и др. [17]. Описание модифицированной версии этого потенциала можно найти в нашей работе [18]. Параметры модельных бикристаллитов представлены ниже в табл. 1.

После построения составного бикристалла методом молекулярной статики проводили релаксацию системы при нулевой температуре. При необходимости изменения температуры или состава размер релаксированного бикристалли-та изменялся в соответствие с параметром решетки (а0), который обеспечивает используемый потенциал в идеальном монокристалле такого же состава и температуры. Измененный бикристал-лит релаксировали методом молекулярной динамики с использованием термостата Берендсена в течении 10 пс.

Удельную энергию границы рассчитывали по формуле

Уоб = (еоб - ее ■ п)/Я , (1)

где EGB - энергия модельного бикристаллита, Ec - энергия связи идеального кристалла на один атом, п - число атомов в системе, 5 - площадь ГЗ в модельном бикристаллите.

В ходе релаксации бикристаллита атомы в нескольких атомных слоях, прилегающих с обеих сторон к ГЗ, заметно смещаются из узлов кристаллических решеток. Эти атомные слои образуют зернограничную область. Размер зерногранич-ной области определялся методом, основанном на анализе потенциальной энергии кристаллита в прилегающих к ГЗ областях (см., например, работу [10]). В этом методе исследуемый бикристаллит разбивается, начиная от ГЗ, на тонкие слои параллельные плоскости распространения ГЗ. В каждом таком слое вычисляется значение потенциальной энергии Е, приходящееся на один атом. Слои не содержащие атомов игнорируются. Граница зернограничной области определяется как плоскости посередине между двумя слоями, в ко-

нокристаллит I Монокристаллит II

Г1 = гп

X

Рис. 1. Схематичное изображение симметрично наклонной границы зерна

ГЗ

Фиксированные атомы

Рис. 2. Схема задания расчетной области с симметрично наклонной границы зерна

Таблица 1. Характеристики ГЗ. T- температура, у - удельная энергия границы, wGB - ширина зернограничной области. В скобках представлена оценка других исследователей.

Размер бикристалла, aoxaoxao Число атомов Чистое Fe, T = 0K Fe - 9ат.%Сг

T = 0K T = 300K

jGB, Дж/м2 Wgb, Á jGB, Дж/м2 Wgb, Á jGB, Дж/м2 Wgb, Á

80.5x58.1x59.0 550706 112 138 (1.39[19]) ' 1.18±0.02* 11.9±1.3 1.11±0.04 11.7±1.2

* Здесь и далее погрешности носят статистический характер и определяются как одно стандартное отклонение (1СТ) для выборочного среднего

торых отклонение Ep от Ec соответственно превышает и не превышает заданной величины А Е. В данной работе мы рассматривали слои шириной 0.2 А и АE = 0.01 эВ.

3. РЕЗУЛЬТАТЫ

Оценки удельных энергий рассматриваемой ГЗ и размера зернограничной области представлены в табл. 1. В табл.1 также приводятся оценки удельной энергии ГЗ в чистом железе, полученные ранее другими исследователями. Несмотря на имеющиеся отличия, в целом наблюдается удовлетворительное согласие наших расчетов с результатами других исследователей. Согласно полученным результатам происходит небольшое (на 0.06 Дж/м2) повышение энергии ГЗ при переходе от чистого железа к сплаву Ре-9ат.%Сг. С повышением температуры до 300 К удельная энергия ГЗ снижается на 0.07 Дж/м2. Ширина зернограничной области составляет ~ 1 нм, слабо меняется с изменением состава и практически не меняется с ростом температуры системы от 0 до 300 К.

Для практически важных приложений связанных с изучением радиационной повреждаемости представляют интерес энергии формирования точечных дефектов (вакансии и СМА) и энергии замещения атомов матрицы примесными атомами, а также энергии связи дефектов с зер-нограничной областью.

В настоящем исследовании рассчитаны энергии связи с зернограничной областью вакансии, СМА, и Сг как атома замещения. Все оценки проведены в матрице чистого Fe при нулевой температуре. Наиболее энергетическая выгодная собственная междоузельная конфигурация в железе и хроме имеет вид гантели с ориентацией <110> (см., например, работы [16, 20] и дальнейшие ссылки в них). Поэтому здесь рассматриваются только гантельные конфигурации с ориентацией <110> двух типов: Fe-Fe и Fe-Cr. Энергия связи рассчитывалась как разность

EGB - Ebulk EGB

(2)

энергий формирования вдали от ГЗ (Efuk) и

вблизи ГЗ (Е^). Энергии формирования вакансии и СМА рассчитывали по формуле

Ef

EGB

Ed

EG

+ k ■ Ec,

(3)

'вв вв

где Е'вв - энергия реласксированного бикрис-таллита с точечным дефектом (вакансия или СМА), Евв - энергия реласксированного бик-ристаллита без точечных дефектов,

k-

+1, для вакансии -1, для СМА ,

Ес - энергия связи на один атом идеального кристалла Бе ( ЕРг ) для случая вакансии или СМА

типа Fe-Fe либо Сг ( Ec ) для СМА смешанного

типа Fe-Cr. Вблизи ГЗ энергия формирования СМА зависит от его ориентации относительно границы. Мы рассматривали все неэквивалентные ориентации и в качестве искомой энергии выбирали наименьшее из полученных значений. Энергию замещения вычисляли по формуле

T^subst ECr

E

Fe+1Cr

EFe + EF

E

Cr

(4)

где EPe+1Cr - энергия кристаллита железа с одним атомом замещения хрома, EFe - энергия кристаллита железа содержащего такое же количество атомов.

Для используемого потенциала значения Ec составляют -4.013 эВ/атом для Fe [16] и -3,837 эВ/атом для Сг [20]. В монокристалле Fe релак-сированная энергия формирования вакансии равна 1,72 эВ [16], энергии формирования СМА с ориентацией <110> составляют 3.51 и 3.31 эВ для гантелей Fe-Fe и Fe-Cr соответственно [18], энергия замещения одного атома железа атомом хрома равна -0.108 эВ [18].

Полученные энергии связи с зернограничной областью вакансии, СМА двух типов и атома замещения Сг представлены на рис. 3-6 соответственно.

Для полученных оценок наблюдаются скачкообразные изменения по мере удаления от ГЗ. Причем в случае вакансии и атома замещения эти скачки происходят со сменой знака. Из получен-

1.О п

m

о

s"

«

к m

0

к s

£1 Ф

1

О

О. О

-1.О -I

О

1О

2 4 6 8

Расстояние от ГЗ, A

Рис. 3. Энергия связи вакансии с зернограничной областью в бикристаллите железа

З.О

m

о

3Ï 2. О м

CÛ

о

к s

il 1.О о

О

О. О -I

О 2 4 б В 1О

Расстояние от ГЗ, A

Рис. 4. Энергия связи СМА Fe-Fe<110> с зернограничной областью в бикристаллите железа

З.О

m

о

JS 2. О

M

s

о

к s

Л1.О ф

о

О.О 4

2 4 б В

Расстояние от ГЗ, A

1О

Рис. 5. Энергия связи СМА Fe-Cr<110> с зернограничной областью в бикристаллите железа

ных результатов следует, что вакансия имеет положительную энергию связи с ГЗ. Максимальное значение энергии связи с межзеренной областью составляет 0.5 эВ. Наименьшие значения не опускаются ниже -0.08 эВ. СМА обоих типов имеют положительную энергию связи с межзерен-ной областью достигающую 2.9 эВ. При этом максимальное значение Eqb для гантелей Fe-Fe

превосходит соответствующее значение для смешанных Fe-Cr СМА на 0.1 эВ.

4. ЗАКЛЮЧЕНИЕ

В работе рассмотрена модель симметрично наклонной границы Z 5(210) с осью вращения [001] в железе и случайном сплаве Fe-9ат.%Cr. Моделирование проведено методом молекуляр-

О.1 п

О. О

m

о

s

M

к m

о -О.1 к

s |_

а а)

X -0.2

О

-О.З

2 4 б В

Расстояние от ГЗ, A

1О

Рис. 6. Энергия связи атома замещения Сг с зернограничной областью в бикристаллите железа

ной статики с использованием многочастичного межатомного взаимодействия. Получены оценки удельной энергии ГЗ и размера соответствующей зернограничной области. Для чистого железа при нулевой температуре отмечено удовлетворительное согласие результатов расчетов с результатами других исследователей. Установлено, что при переходе от чистого железа к сплаву Бе-9ат.%Сг происходит незначительное повышение удельной энергии ГЗ (на ~0.1 Дж/м2). С повышением температуры сплава до 300 К удельная энергия ГЗ незначительно снижается на ~0.1 Дж/м2. Ширина зернограничной области составляет ~ 1 нм, слабо меняется с изменением состава и практически не меняется с ростом температуры системы.

Для бикристалла железа рассчитаны энергии связи с зернограничной областью вакансии, двух типов междоузельных конфигураций (гантель с

ориентацией <110> из двух атомов железа и гантель смешанного Fe-Cr типа той же ориентации) и замещающего атома Cr. Отмечены скачкообразные изменения полученных оценок по мере удаления от ГЗ, вплоть до периодической смены знака в случае вакансии и атома замещения. Обнаружена высокая положительная связь (~3 эВ) СМА обеих конфигураций с зернограничными областями всех рассматриваемых типов. Присоединение к зернограничной области вакансии так же следует рассматривать как энергетически выгодное (энергия связи 0.5 эВ). Эти результаты обосновывают свойство зернограничной области аккумулировать точечные дефекты в процессе радиационного повреждения материала, отмеченное, в частности, в работах [3, 4, 7, 21]. При этом накопление междоузельных атомов происходит более интенсивно как по причине более высокой энергии связи, так и вследствие их высокой подвижности.

Согласно полученным результатам присоединение с зернограничной области замещающего атома хрома нельзя однозначно трактовать как энергетически выгодное или невыгодное. Это объясняется сменой знака энергии связи при рассмотрении различных позиций в зерногранич-ной области. Чаще эта энергия все же отрицательна. Этот результат следует сопоставить с работой [11], где получены соответствующие оценки энергии связи для границ Z19(331), Z 9(221), Z 3(111), Z 3(112), Z 11(113) и Z 9(114) с осью вращения [110], а расчет проведен с двух-полосным (two band) многочастичным эмпирическим потенциалом. Авторами этой работы установлена положительная энергия связи замещающего атома Cr со всеми рассмотренными границами. Однако и в этом случае наблюдаются отрицательные значения энергии связи для отдельных позиций вблизи ГЗ. Поэтому, учитывая различие в типах рассмотренных границ, можно признать качественное согласие наших результатов с результатами из ранее опубликованных работ.

Работа выполнена при поддержке Минобрна-уки в рамках государственного задания на 20122014 гг.

СПИСОК ЛИТЕРАТУРЫ

1. J.-M. Zhang et al. Energy calculation for symmetrical tilt grain boundaries in iron // Applied Surface Science 252 (2006) p. 4936-4942

2. W.Y. Hu, B.W. Zhang, X.L. Shu, B.Y. Huang, Calculation of formation enthalpies and phase stability for Ru-Al alloys using an analytic embedded atom model // J. Alloys Comp. 287 (1999) p. 159-162

3. F.Javier Pérez Pérez, Roger Smith, Modelling radiation

effects at grain boundaries in bcc iron / / Nuclear Instruments and Methods in Physics Research B 153 (1999) p. 136-141

4. S.F. Di Martino, R.G. Faulkner, R. Smith, Modelling radiation damage effects on a bcc iron lattice containing phosphorous impurity atoms near symmetrical tilt boundaries // Journal of Nuclear Materials 417 (2011), p. 1058-1062

5. Yan-Ni Wena, Yan Zhang, Jian-Min Zhang, Ke-Wei Xu, Atomic diffusion in the Fe [0 0 1] S= 5 (3 1 0) and (2 1 0) symmetric tilt grain boundary / / Computational Materials Science 50 (2011) p. 2087-2095

6. Y. Shibuta, Sh. Takamoto, T. Suzuki, A Molecular Dynamics Study of the Energy and Structure of the Symmetric Tilt Boundary of Iron // ISIJ International, Vol. 48 (2008), No. 11, p. 1582-1591

7. F.Javier Pérez Pérez, Roger Smith, Preferential damage at symmetrical tilt grain boundaries in bcc iron // Nuclear Instruments and Methods in Physics Research B 180 (2001) p. 322-328

8. Ning Gao, Chu-Chun Fu, M. Samaras, R. Schaublin, M. Victoria, W. Hoffelner, Multiscale modelling of bi-crystal grain boundaries in bcc iron // Journal of Nuclear Materials 385 (2009) 262-267

9. E. Wachowicz, T. Ossowski, and A. Kiejna, Cohesive and magnetic properties of grain boundaries in bcc Fe with Cr additions // // Phys. Rev., B 81, 094104, 2010

10. С. Г. Псахье, К. П. Зольников, Д. С. Крыжевич, А. В. Железняков, В. М. Чернов, Развитие каскадов атомных соударений в кристалле ванадия с внутренней структурой // Кристаллография, 2009, том 54, № 6, с. 1053-1062

11. D. Terentyev, X. Hea, E. Zhurkin, A. Bakaev, Segregation of Cr at tilt grain boundaries in Fe-Cr alloys: A Metropolis Monte Carlo study // Journal of Nuclear Materials 408 (2011) p. 161-170

12. Тихончев М.Ю., Светухин В.В. Расчетное определение пороговых энергий смещения и исследование особенностей развития каскадов атомных смещений вблизи протяженной границы раздела фаз циркония и ниобия: молекулярно-динамическое моделирование // Вопросы материаловедения. 2011. № 4(68). C. 140-152

13. Bloom E.E., Zinkle S.J., Wiffen F.W., Nucl J. // Mater. 329-333 (2004) 12.

14. L.K. Mansur, A.F. Rowcliffe, R.K. Nanstad, S.J. Zinkle, W.R. Corwin, R.E. Stoller // J. Nucl. Mater. 329-333 (2004) 166.

15. N. Baluc, Plasma Phys. Control. Fusion 48 (2006) B165.

16. G.J. Ackland, M.I. Mendelev, D.J. Srolovitz, S.W. Han, A.V. Barashev. Development of an interatomic potential for phosphorus impurities in a-iron. J. Phys.: Condens. Matter 16 (2004) pp. S2629-S2642

17. A. Caro, D. A. Crowson, and M. Caro, Phys. Rev. Lett. 95 (2005) 075702, 4 pp.

18. M. Tikhonchev, V. Svetukhin, E. Gaganidze, MD simulation of atomic displacement cascades near chromium-rich clusters in FeCr alloy // Journal of Nuclear Materials, 442 (2013) S618-S623.

19. D. Terentyev, Xinfu He. Properties of grain boundaries in BCC iron and iron-based alloys: An atomistic study. Open report of the Belgian nuclear research centre SCK^CEN-BLG-1072, Belgium (2010). 70 p.

20. J. Wallenius, P. Olsson, C. Lagerstedt, N. Sandberg, R.

Chakarova, V. Pontikis, Phys. Rev. B 69 (2004) 094103. 21. Муралев А.Б., Тихончев М.Ю, Светухин В.В. Моделирование каскадов атомных смещений в альфа-железе,

содержащем симметрично-наклонную межзеренную границу // Известия ВУЗов: Поволжский регион. Физико-математические науки, №1, 2013 с. 144 - 158.

ENERGY PROPERTIES OF INTERACTION BETWEEN POINT DEFECTS AND 25 (210) [001] TILT GRAIN BOUNDARY IN FECR ALLOY

© 2013 M.Yu. Tihonchev, A.B. Muralev, V.V. Svetukhin

Research Technological Institute named after S.P. Kapitsa of Ulyanovsk State University

This paper studies the energy properties of symmetrically inclined grain boundary 2 5(210)[ 001] in a-Fe, containing chromium. The study is carried out by the method of molecular statics simulations based on a many-interatomic interaction . The estimates of the specific energy boundaries and sizes of the corresponding grain boundary regions are obtained in pure iron at zero temperature and in the binary alloy Fe- 9at .% Cr at temperatures between 0 and 300 K. The energy of the formation and binding of vacancy, self-interstitial atoms, and substitute Cr atoms with grain boundary region in iron are evaluated. Two types of self-interstitial dumbbell configurations with orientation of < 110 >: Fe-Fe and mixed Fe-Cr dumbbells. It is revealed that the binding energy of interstitial configurations with grain boundary regions is positive. The binding energy of a vacancy with a majority position in the grain boundary region is positive too. The binding energy of the substitution Cr atom changes sign depending on the position in the boundary region, the maximum positive binding energy is of 0.05 eV.

Keywords: energy properties, symmetrically inclined grain boundary, Fe-Cr alloy.

This work has been partly supported by the Russian Ministry of Education and Science in frame of government assignment in 2012-2014.

Mikhail Tikhonchev, Candidate of Physics and Mathematics,

Head at the Non-Organic Materials Behaviour Modelling

Laboratory. E-mail: [email protected]

Artem Muralev, Associate Research Fellow at the Computer

Modeling Laboratory of Inorganic Materials.

E-mail: [email protected]

Viacheslav Svetukhin, Doctor Physics and Mathematics, Professor, Director of Research Institute of technology. E-mail: [email protected]