Гибридный потенциал межчастичного взаимодействия и расчет линии плавления лития методом молекулярной динамики

Для металлических систем предложен многочастичный потенциал нового типа (ЕАМ-2), включающий в качестве параметра степень кристалличности атомов. Этот потенциал базируется на потенциале модели погруженного атома (ЕАМ). Степень кристалличности атома определяется путем разложения радиус-векторов соседних атомов по шаровым функциям и вычисления характеристики $q_6$. Потенциал ЕАМ-2 хорошо описывает свойства жидкого, ОЦК- и ГЦК-лития при обычных и высоких давлениях. Рассчитаны линии равновесия ОЦК-литий–жидкость и ГЦК-литий–жидкость при давлениях до $40$ ГПа и получено хорошее согласие с опытом, в том числе для максимума на кривой плавления. Приведены таблицы зависимости плотности и энергии моделей от давления и температуры. Обнаружено аномальное поведение плотности ОЦК- и ГЦК-лития в интервале давлений $25$–$30$ ГПа (рост плотности при изобарном нагреве).