Моделирование методом Монте-Карло распада кластеров хрома в FeCr

УДК 621.039.531.001.57

МОДЕЛИРОВАНИЕ МЕТОДОМ МОНТЕ-КАРЛО РАСПАДА КЛАСТЕРОВ ХРОМА В FeCr

© 2013 Д.Я. Вострецов, А.С. Кадочкин

Ульяновский государственный университет

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

В работе проведено моделирование методом Монте Карло процесса распада кластеров хрома в сплаве FeCr. При моделировании учитывался начальный размер кластера, температура, концентрация хрома в матрице. Для расчетов был использован многотельный потенциал межчастичного взаимодействия. В результате проведенных расчетов получены зависимости кинетики рапспада кластеров при разных температурах. Кроме того, было проведено моделирование распада кластера хрома путем решения диффузионного уравнения.

Ключевые слова: моделирование, метод Монте Карло, распад кластеров хрома, сплав FeCr.

ВВЕДЕНИЕ

Сплав FeCr является основой феррито-мар-тенситных высоко-хромистых сталей, которые, обладая хорошими термическими, механическими свойствами, высокой радиационной стойкостью, широко применяются в ядерной энергетике [1-8]. Данная работа посвящена моделированию распада кластеров хрома в сплаве FeCr. Моделирование проводится двумя способами: методом Монте-Карло и путем решения уравнения диффузии. Метод Монте-Карло является одним из общепринятых методов моделирования, в котором материал рассматривается как система, состоящая из отдельных атомов. Для описания взаимодействий атомов в металлах и сплавах обычно используют многотельные потенциалы межатомного взаимодействия, такие как метод погруженного атома и апроксимация второго момента [9, 10]. В диффузионном приближении материал рассматривается как непрерывная среда. Поэтому такое приближение корректно только при рассмотрении достаточно больших объектов, когда атомарным строением вещества можно принебречь.

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

В работе для моделирования распада кластеров хрома в альфа-железе был использован метод Монте-Карло. При моделировании были использованы многотельные потенциалы межа-

Вострецов Дмитрий Ярославович, кандидат физико-математических наук, старший научный сотрудник Научно-исследовательского технологического института. E-mail: Dimka97@mail.ru

Кадочкин Алексей Сергеевич, кандидат физико-математических наук, старший научный сотрудник Научно -исследовательского технологического института. E-mail: askadochkin@sv.ulsu.ru

томного взаимодеиствия.

Многотельный потенциал взаимодеИствия имеет вид [9, 10]:

Е = {г ) + Р [Ър, )] ,

j * i

V j *i

где ^ ] - позиции для взаимодействующих атомов A и B, V - парный потенциал взаимодеИствия, F- функция внедрения, р - электронная плотность.

Для железа был использован многотельный потенциал из работы Акланда, Менделева и др. [11], который получен из потенциала, ранее предложенного в работе [12], путем небольшой модификации его функции внедрения. Этот потенциал считается на сегодняшний день одним из самых удачных и широко используется различными исследователями. Для хрома мы использовали потенциал Валлениуса и др. из работы [13].

Парный потенциал для взаимодействия атомов железа и хрома вычисляется по формуле [14]:

УАВ { х, г ) = И { х ) 2 [УАА { г ) + Увв { г )] ,

где И (х) - полином четвертого порядка.

Для вычисления энергий взаимодействия был использован метод нормализации электронных плотностей, описанный в работе [14]. Используемые в данной работе потенциалы имеют радиус обрезания между 5-ой и 6-ой координационными сферами включительно. Ниже изложен алгоритм моделирования.

Первый шаг. Задается начальная конфигурация атомов. В ОЦК решетке создается сферический кластер с необходимыми параметрами.

Второй шаг. Рассчитывается общее количество скачков атомов для заданной конфигурации решетки при заданной температуре и времени. Время моделирования рассматриваемой системы

определяется следующим образом. Вероятность перемещения атома определяется следующим выражением:

С

Р = А • ехр

Е - ЕЬ

а Ь

кТ

V У

где Р - вероятность перескока, Еа ЕЬ - энергии соответствующих состояний атома, к - постоянная Больцмана, Т - температура. Коэффициент А в (3) вычисляется по формуле:

А = 80Г'

где Б - коэффициент диффузии, а - длина скачка, А £ - промежуток времени, для которого рассчитывается вероятность перемещения атома. При этом величина А £ задается таким образом, чтобы вероятность (3) была меньше единицы (оптимальное значение для максимальной вероятности в системе Р = 0.5 - 0.7 ). Коэффициент диффузии для хрома в железе был взят нами из работы [15]:

Э = 2.33 • ехрI -

57070 ЯТ

см 2/с.

^=2

Ж

- Э

да

Сг

Сг

дг

Третий шаг. Проводится моделирование скачков атомов. Для этого необходимое число раз перебираются координаты рассматриваемых атомов. Затем для рассматриваемых атомов разыгрывается направление скачка в пределах первой координационной сферы.

Четвёртый шаг. Разыгрывается вероятности перемещения атома. Для этого рассчитываются энергии рассматриваемого атома до и после перемещения. Вероятность определяется из выражения (3).

Третий и четвёртый шаги повторяются необходимое число раз до достижения заданного времени моделирования.

РЕШЕНИЕ ДИФФУЗИОННОГО УРАВНЕНИЯ

Рассмотрим задачу о вычислении скорости изменения размера сферических кластеров в двух-компонентном твердом растворе замещения БеСг. Предполагается, что кластер на 80% процентов состоит из атомов хрома т.е. СР = 0.8С0, где СР -концентрация атомов хрома в кластере, С0 - концентрация частиц в матрице железа без примеси, концентрация атомов хрома в матрице на большом расстоянии от кластера См полагается равной нулю. Скорость, с которой изменяется количество атомов хрома в кластере радиуса Я , может быть определено с помощью соотношения:

В литературе традиционно используются два приближения для анализа данного уравнения. Наиболее простым является рассмотрение стационарного случая, когда по мере роста кластера распределение концентрации вокруг кластера определяется исходя из решения стационарного уравнения диффузии с учетом сферической симметрии диффузионных потоков. Вторым более точным приближением является приближение стационарной концентрации на границе между кластером и матрицей [16], которое приводит к зависимости концентрации атомов примеси в матрице от расстояния до центра кластера г и времени:

Ссг М= См +

здесь С1 - равновесная концентрация атомов примеси в матрице железа. Точное решение уравнения диффузии в сферических координатах можно также найти в работе [16].

Приближение стационарных диффузионных потоков является одним из традиционных приближений, используемых в кинетике роста частиц второй фазы (см. например [17]). В соответствии с данным приближением концентрация атомов сорта В определяется решением стационарного уравнения диффузии с учетом сферической симметрии диффузионных потоков атомов вокруг сферического кластера:

Ссг

г

что соответствует пределу ^ —> го в формуле (7).

Таким образом, скорость изменения количества частиц сорта В с учетом (6) и (8) оказывается равной:

ж

Сг

= 4пт(См - С, (Я)).

Поскольку Ср > С,, то рост кластера будет наблюдаться при См > С1, растворение кластера См < С, . Поскольку в рассматриваемом случае См = 0 , то выполняется второе условие, что соответствует растворению кластеров.

Скорость изменения числа частиц примеси в кластере за счет их диффузионного переноса имеет очевидную связь с радиусом кластера и его составом:

ЖЫСг (Я ) = Ж

4кЯ3

,

Подставляя уравнение (9) в (8), получим соотношение, в котором могут быть разделены переменные:

г = Я

R2 д

7 'т дя ^ с ^+R(Cf с я с:^^

dR = Dct ,

где Я0 - радиус кластера в момент времени t = 0.

В случае нестационарного приближения необходимо рассматривать в уравнении градиент концентрации атомов примеси с учетом зависимости от времени. В этом случае вместо уравнения (10) получим

N

Cr

dt

=4 кЩСм - С,

1 1

— +-

r

или с учетом выражения (10) может быть получено нелинейное дифференциальное уравнение:

r dr [с

= D

R 4nD

t

dR

Данное нелинейное уравнение, в отличие от уравнения (11) не допускает разделения переменных и может быть решено численно.

РЕЗУЛЬТАТЫ МОДЕЛИРОВАНИЯ

На рис. 1 показаны стадии распада кластера хрома в матрице железа. Изображения получены при моделировании методом Монте- Карло.

На рис. 2 показана кинетическая кривая растворения кластера хрома для различных температур, полученная методом Монте-Карло и при помощи решения диффузионного уравнения. На рисунках R - радиус кластера, R - начальный радиус, t -время в секундах. Начальный радиус кластера 1нм, начальное содержание хрома в кластере 90%, начальное содержание хрома в матрице равно нулю. Следует отметить, что результаты, полученные методом Монте-Карло хорошо согласуются с результатами, полученными из диффузионного приближения. Время уменьшения числа частиц в кластере вдвое для Т=450°С составило больше 1000 часов, для Т=475°С - 260 часов, для Т=500°С - 55 часов.

ЗАКЛЮЧЕНИЕ

В работе проведено моделирование процесса распада кластеров хрома в сплаве FeCr методом Монте-Карло и в диффузионном приближении.

*4

л* >

'sr.^rvi-. t ; * ...  t

Рис. 1. Стадии распада кластера хрома. Температура - 475°С. а - начальный кластер, Ь - кластер в процессе распада

Рис. 2. Кинетика растворения кластера хрома. Точки - моделирование методом Монте-Карло, линии - диффузионная модель растворения кластера. (1) - 450°С, (2) - 475°С, (3) - 500°С

При моделировании учитывался начальный размер кластера, температура сплава, концентрация хрома в кластере и матрице. Получены зависимости, характеризующие кинетику распада кластеров при разных температурах. Показано хорошее согласие результатов, полученных методом Монте-Карло и диффузионным методом.

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

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

1. Kinetics and thermodynamics of Cr nanocluster formation in Fe-Cr system / V. Svetukhin, P. Lvov, E. Gaganidze, M. Tikhonchev, C. Dethloff // Journal of Nuclear Materials, 415, 2011, pp. 205 - 209

2. MD simulation of atomic displacement cascades in Fe-10 at.%Cr binary alloy / M. Tikhonchev, V. Svetukhin, A. Kadochkin, E. Gaganidze //Journal of Nuclear Materials 395 (2009) p. 50-57.

3. Моделирование процессов первичной радиационной повреждаемости сплава Fe-1.8ат.%Ni методом молекулярной динамики / М.Ю. Тихончев, В.В. Светухин, Д.В. Козлов, В.Н. Голованов // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. 2010. №3 (15). С. 143 - 155.

4 Светухин В.В., Тихончев М.Ю. Моделирование взаимодействия каскадов атомных смещений с обогащенными хромом преципитатами в сплаве FeCr // Известия ВУЗов: Поволжский регион. Физико-математические науки. 2012. №4. С. 162 - 173.

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

6. Modeling of chromium nanoclusters growth under neutron irradiation / V. Svetukhin, P. L'vov, M. Tikhonchev, E. Gaganidze, N. Krestina // Journal of Nuclear Materials, 442 (2013) S624-S627

7. Львов П.Е., Светухин В.В. Термодинамика фазового равновесия многокомпонентных твердых растворов // Физика твёрдого тела. 2013. Т. 55. С. 2256.

8. Моделирование образования кластеров в сплавах на основе Fe-Cr в процессе термического отжига и под облучением / В.В. Светухин, П.Е. Львов, E. Gaganidze, Н.С. Крестина // Вопросы материаловедения. 2013. № 1(73). 230-240.

9. Daw M.S., Baskes M.I. // Phys. Rev. B 29, 6443 (1984).

10. FinnisM.F., SinclairJ.E. // Philos. Mag., A 50, 45 (1984).

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

12. Development of new interatomic potentials appropriate for crystalline and liquid iron / M.I. Mendelev, S. Han, D.J. Srolovitz, G.J. Ackland, D.Y. Sun, M. Asta// Philos. Mag., V. 83, N. 35 (2003), pp. 3977-3994.

13. Modeling of chromium precipitation in Fe-Cr alloys /J. Wallenius, P. Olsson, C. Lagerstedt, N. Sandberg, R. Chakarova, and V. Pontikis// PHYSICAL REVIEW B 69, 094103, 2004.

14. Classical Many-Body Potential for Concentrated Alloys and the Inversion of Order in Iron-Chromium Alloys /A. Caro, D. A. Crowson, and M. Caro // PHYSICAL REVIEW LETTERS, 95, 075702, 2005.

15. Коэффициент диффузии хрома и ближний порядок в сплавах Fe-Cr / Д.А. Мирзаев, А.А. Мирзоев, К.Ю.

Окишев, А.С. Созыкина// Вестник ЮУрГУ. 2009. №14. 17. Ландау Л.Д., Лифшиц Е.М. Статистическая физика. 16. H.B. Aaron, D. Fainstein, G.R. Kotler // J. Appl. Phys. Ч.2. Теория конденсированного состояния . М.: Физ-

1970. V.41. P.4404. матлит, 2001. 493 с

MONTE CARLO SIMULATION DECAY OF CHROMIUM CLUSTERS IN FeCr

© 2013 D.Y. Vostretsov, A.S. Kadochkin

Ulyanovsk State University

In this paper, a simulation of decay process of chromium clusters in FeCr alloy by Monte Carlo Method is performed. Taken into account when modeling the initial cluster size , temperature, concentration of chromium in the matrix. For calculations a multibody potential interparticle interaction was used. As a result of these calculations the dependences of kinetic decay of clusters at different temperatures are obtained. In addition, the simulation cluster decay chromium was performed by solving the diffusion equation. Key words: simulation, Monte Carlo Method, decay process of chromium clusters, FeCr alloy.

This work was supported by the Ministry of Education in the framework of the public task for 2012-2014.

Dmitry Vostretsov, Candidate of Physics and Mathematics, Senior Research Fellow at the Research Institute o f Technology. E-mail: Dimka97@mail.ru

Alex Kadochkin, Candidate of Physics and Mathematics, Senior Research Fellow at the Research Institute of Technology. E-mail: askadochkin@sv.ulsu.ru