Примеси кремния в ОЦК-железе: ab-initio моделирование свойств и энергетических параметров

В программном пакете WIEN2k проведено первопринципное моделирование равновесной структуры и свойств атомов кремния, растворенных в ОЦК-железе. Для расчетов использовался полнопотенциальный метод линейных присоединенных плоских волн LAPW, с учетом обобщенного градиентного приближения PBE-GGA, в суперячейке из 54 атомов железа с периодическими граничными условиями. Это наиболее точный метод, используемый в рамках теории функционала плотности. В работе определены оптимальные значения основных параметров моделирования, позволяющие рассчитывать энергетические характеристики системы с точностью не менее 0,01 эВ. Энергия растворения кремния в ферромагнитное ОЦК-железе составила -1,19 эВ. Атомы кремния испытывают сильное взаимное отталкивание в первых двух координационных сферах, которое далее является спадающим. После третей координационной сферы отталкивание атомов становится близким к нулю. При растворении кремний не изменяет параметр ОЦК-решетки и не влияет на решетку железа. Магнитный момент атомов железа в первом окружении с 2,23 µB, уменьшается до 2,045 µB. Растворение 1,85 ат. % кремния приводит к уменьшению среднего магнитного момента на атом железа на 0,02 μB и снижению объемного модуля сжатия.