Влияние энергий взаимодействия между атомами углерода в ГЦК-железе на концентрационную зависимость активности углерода

В программном пакете WIEN2k проведён расчёт из первых принципов энергий взаимодействия между атомами углерода в ГЦК-железе. Они составили E1 = 0,06 эВ, E2 = 0,1 эВ и E3 = 0,005 эВ (энергии взаимодействия между атомами углерода, находящимися в первой, второй и третьей координационных сферах подрешётки ГЦК-железа). Для расчётов использовался полнопотенциальный метод линейных присоединённых плоских волн LAPW с учетом обобщенного градиентного приближения PBE-GGA в суперячейке из 32 атомов железа с периодическими граничными условиями. Это наиболее точный метод, используемый в рамках теории функционала плотности. Для проверки достоверности найденных значений энергий на их основе методом Монте-Карло построена концентрационная зависимость активности углерода в ГЦК-железе. Показано хорошее качественное согласие кривой активности, полученной при использовании найденных значений энергий взаимодействия C-C, с экспериментальной кривой, что говорит о достоверности полученных параметров взаимодействия и адекватности использованной модели парамагнитного состояния. Использование данных о взаимодействии в третьей координационной сфере не влияет существенно на результаты расчета. Это объясняется очень низким значением энергии взаимодействия между атомами углерода в ГЦК-железе в третьей координационной сфере.