Методами теории функционала плотности исследована растворимость 3p-(Al, Si) и 4р-(Ga, Ge) элементов в ферромагнитном и парамагнитном состояниях ОЦК-железа. Для моделирования парамагнитного состояния использовали два различных подхода, реализованные в пакетах SIESTA и LSGF. Установлено, что энергия растворения всех рассмотренных примесей понижается при переходе в парамагнитное состояние. Энергии растворения, полученные усреднением по ансамблю неполяризованных магнитных конфигураций, хорошо согласуются со значениями, полученными в приближении когерентного потенциала. В то же время учет магнитной поляризации в окрестности примеси приводит к понижению энергии растворения, которое особенно выражено при температурах, близких к ТC. Обсуждается температурная зависимость энергий растворения в парамагнитном состоянии для рассматриваемых примесных элементов.