Монте-Карло моделирование активности углерода в ОЦК-железе

Задача исследования заключалась в проверке предложенных в литературе значений энергии взаимодействия между атомами углерода в ОЦК-железе для различных координационных сфер. Цель работы: рассчитать активность углерода в ОЦК-железе, используя различные параметры взаимодействий между атомами углерода. Методы: для моделирования использовалась самописная программа по Монте-Карло, в которой была рассмотрена суперячейка 20×20×20 элементарных ячеек объёмно-центрированного (ОЦК) железа. Нами была проведена оптимизация, в ходе которой была получена сходимость результата усреднения при 500 шагах Монте-Карло. Расчеты проводили для широкого диапазона температур (Т = 955, 975, 1000, 1026, 1056, 1070, 1086, 1121 К). Результаты: наши расчеты показали, что почти все предложенные параметризации согласуются с экспериментальными данными, поскольку при низких концентрациях атомы углерода практически не взаимодействуют. Чтобы изучить влияние взаимодействия атомов углерода на кривую активности, нам пришлось расширить наши расчеты до области более высоких концентраций, чем наблюдаются на самом деле. Это позволило нам сравнить различные параметризации и получить ряд интересных выводов. Кроме того, мы продемонстрировали возможность воспроизведения экспериментальных кривых активности углерода для нескольких наборов параметров С-С взаимодействия. Новизна работы: мы продемонстрировали, что для расчета активности углерода в ОЦК-железе достаточно вычислить энергии С-С взаимодействия для первых четырех координационных сфер. Практическая значимость заключается в получении новой теоретической информации, необходимой для разработки новых марок сталей и режимов термической обработки их.