Компьютерное моделирование никеля и учет электронных вкладов в методе молекулярной динамики

Предложены два новых потенциала модели погруженного атома для никеля с учетом и без учета теплового электронного вклада в энергию. Параметры потенциалов найдены по свойствам никеля на изобаре $p = 0$ и в условиях ударного сжатия до давлений $\sim760$ ГПа. Наилучшее согласие с экспериментом получено при учете электронного вклада в энергию. Расчетная адиабата Гюгонио хорошо согласуется с фактической. Учет электронного вклада сильно уменьшает температуру на адиабате и увеличивает “холодное” давление. Расчетная линия плавления моделей никеля плавно поднимается до $4518$ К при $300$ ГПа и почти нечувствительна к наличию электронного вклада. Определены координаты участка плавления на ударной адиабате (начало – $275.8$ ГПа и $4422$ К, конец – $297.6$ ГПа и $4499$ К). Предложена фазовая диаграмма никеля с областью устойчивой ОЦК-фазы при давлениях выше $110$–$130$ ГПа. Приведены таблицы энергии и давления моделей при степенях сжатия до $1.8182$.