Анализ фазовых превращений в системе Fe—Zr был проведён с помощью ab initio молекулярной динамики и термодинамического подхода. Представлено первопринципное моделирование методом молекулярной динамики процессов аморфизации и кристаллизации в системе Fe—Zr. Положения атомов в сверхъячейке Fe$_{29}$Zr$_{3}$ моделировались путём числового отжига методом функционала плотности в обобщённом градиентном приближении. Обсуждаются изменения плотности электронных состояний сверхъячейки Fe$_{29}$Zr$_{3}$ при жидкостно—аморфно—кристаллическом фазовом переходе. Наиболее заметным отличием между электронными спектрами жидкой и аморфной фаз является появление псевдощели на уровне Ферми, что коррелирует с электронным критерием термостабильности аморфных металлических сплавов Нагеля—Таука. Дальнейший отжиг в изотермично/изоэнтальпийном ансамбле при более высоких температурах приводит к резкому изменению электронного спектра и перегруппировке атомов, которую мы приписываем первой стадии кристаллизации аморфного сплава. В рамках термодинамического подхода была построена концентрационная зависимость относительной интегральной свободной энергии Гиббса для исходной аморфной фазы; эта зависимость имеет специфическую S-образную форму, что указывает на тенденцию сплавов к фазовому расслоению. Термодинамический метод подтвердил, что процессы формирования фаз в бинарных сплавах системы Fe—Zr проходят в 2 этапа: сначала происходят процессы фазового расслоения на две аморфные фазы по типу спинодального распада, затем – кристаллизации каждой аморфной фазы. Рассчитанная концентрационная область аморфизации бинарных сплавов системы Fe—Zr хорошо согласуется с экспериментальными данными.