Аб initio молекулярная динамика для жидких металлов

Представлены ab initio квантово-механические расчёты молекулярной динамики, основанные на вычислении электронного основного состояния и сил Хеллмана–Фейнмана в приближении локальной плотности на каждом шаге молекулярной динамики. Это становится возможным благодаря использованию методов сопряжённых градиентов для минимизации энергии и прогнозированию волновых функций для новых положений ионов с помощью выравнивания подпространств. Такой подход избегает нестабильностей, присущих квантово-механическим расчётам молекулярной динамики для металлов, основанным на использовании фиктивной ньютоновской динамики для электронных степеней свободы. Метод обеспечивает полный контроль адиабатичности и позволяет выполнять моделирование в течение нескольких пикосекунд.