Ближний порядок в сплавах Fe-Cr: моделирование методом решеточного Монте-Карло

Abstract

Проведено исследование ближнего порядка в сплавах Fe-Cr методом Монте-Карло. Моделирование проводилось в рамках алгоритма Метрополиса в программном пакете LAMMPS. Анализ данных моделирования осуществлялся с помощью программного обеспечения для визуализации и анализа данных Ovito. Модель сплава предполагает, что структура решетки фиксирована и что взаимодействия существуют между первыми соседями и вторыми соседями. Установлено взаимодействие Fe-Cr с применением потенциала межатомного взаимодействия Abell-Brenner-Tersoff (ABOP). Были изучены различные концентрации примесей внедрения хрома в железе, а именно 5-50 ат. %. Рассчитана энергия смешения системы Fe-Cr при различных концентрациях примесей внедрения. Расчеты показали, что выбранный потенциал взаимодействия верно воспроизводит изменения знака энергии смешения как функции концентрации Cr. При применении в кинетическом моделировании по методу Монте-Карло потенциал правильно предсказывает разложение первоначально случайных сплавов Fe-Cr в зависимости от концентрации Cr. Определен параметр ближнего порядка Каули, который используется для количественной оценки степени упорядочения. В соответствии с экспериментами наблюдается сильная тенденция упорядочения в распределении Cr при низких концентрациях, что проявляется в отрицательных значениях параметров ближнего порядка. Short-range order in Fe-Cr alloys was studied by the Monte Carlo method. The simulation was carried out in the framework of the Metropolis algorithm in LAMMPS software package. Simulation data analysis was carried out using Ovito data visualization and analysis software. The alloy model assumes that the lattice structure is fixed and that interactions exist between the first neighbors and the second neighbors. The FeCr interaction was determined using the Abell-Brenner-Tersoff interatomic interaction potential (ABOP). Various admixture concentrations of chromium interstitial in iron have been studied, namely, 5-50 at. %. The mixing energy of the Fe-Cr system was calculated for various concentrations of interstitial impurities. Calculations showed that the chosen interaction potential correctly simulates changes in the sign of the mixing energy as a function of Cr concentration. When used in Monte Carlo kinetic modeling, the potential correctly predicts decomposition of initially random Fe-Cr alloys, depending on Cr concentration. The Cowley short-range order parameter, which is used for quantitative assessment of the degree of order, has been determined. According to experiments, there is a strong tendency of ordering in distribution of Cr at low concentrations, which becomes evident in negative values of the short-range order parameters.