Модели молекулярной динамики: обзор EAM-потенциалов. Часть 1: потенциалы для однокомпонентных систем Текст научной статьи по специальности «Нанотехнологии»

Волегов П.С., Герасимов Р.М., Давлятшин Р.П. Модели молекулярной динамики: обзор EAM-потенциалов. Часть 1: Потенциалы для однокомпонентных систем // Вестник Пермского национального исследовательского политехнического университета. Механика. - 2017. - № 4. - С. 214-237. DOI: 10.15593/perm.mech/2017.4.14

Volegov P.S., Gerasimov R.M., Davlyatshin R.P. Models of molecular dynamics: a review of EAM potentials. Part 1: Potentials for single-component systems. PNRPU Mechanics Bulletin, 2017, no. 4, pp. 214-237. DOI: 10.15593/perm.mech/2017.4.14

ВЕСТНИК ПНИПУ. МЕХАНИКА

№ 4,2017 PNRPU MECHANICS BULLETIN

http://vestnik.pstu.ru/mechanics/about/ini7

001 10.15593/регш.шесЬ/2017.4.14 УДК 539.3

МОДЕЛИ МОЛЕКУЛЯРНОЙ ДИНАМИКИ: ОБЗОР EAM-ПОТЕНЦИАЛОВ. ЧАСТЬ 1: ПОТЕНЦИАЛЫ ДЛЯ ОДНОКОМПОНЕНТНЫХ СИСТЕМ

П.С. Волегов, Р.М. Герасимов, Р.П. Давлятшин

Пермский национальный исследовательский политехнический университет, Пермь, Россия

АННОТАЦИЯ

Статья является первой частью обзора современных подходов к построению потенциалов межатомного взаимодействия на базе метода погруженного атома (embedded atom method, EAM) для различных однокомпонентных материалов. Приведены основные положения и соотношения метода молекулярной (атомарной) динамики, отмечены недостатки подхода к моделированию атомных систем с использованием классических парных потенциалов взаимодействия. Приведены основные положения и соотношения метода погруженного атома, а также различные модификации метода. Рассмотрен алгоритм использования EAM-потенциалов в расчетах молекулярно-динамических систем, а также способы повышения эффективности вычислительных алгоритмов, в том числе, связанные с введением радиуса обрезания и обоснованием его выбора для тех или иных материалов. В рамках обзора особое внимание уделено способам построения и видам функции погружения для различных материалов, а также физико-механическим свойствам материалов, которые хорошо описываются тем или иным видом потенциала. Приведен краткий обзор способов получения параметров потенциалов: из квантово-механических расчетов (часто в иностранной литературе называемых «из первых принципов», «ab initio») либо с использованием экспериментально определяемых свойств (дифракционными, рентгеновскими методами или при помощи электронной микроскопии), а также перечень физико-механических свойств, которые возможно определять для материала с использованием методов молекулярной динамики с применением современных EAM-потенциалов. Сделана попытка подробно рассмотреть работы, в которых предложены потенциалы, в настоящий момент наиболее точно описывающие свойства таких сложных (для прямого моделирования) материалов, как бериллий, железо, вольфрам, ниобий, титан, уран, и некоторые другие. Актуальность исследований структуры и свойств именно этих материалов обусловлена во многом перспективностью их применения в деталях и конструкциях для атомной энергетики, аэрокосмической промышленности, в биомедицине. Отдельно рассмотрены работы, посвященные построению EAM-потенциалов, позволяющих описывать структуру и свойства ряда металлов (лития, никеля, меди, алюминия и других) как в состоянии расплава, так и после кристаллизации.

© ПНИПУ

О СТАТЬЕ

Получена: 30 августа 2017 г. Принята: 10 декабря 2017 г. Опубликована: 29 декабря 2017 г.

Ключевые слова:

молекулярная динамика, потенциал взаимодействия, метод погруженного атома, ЕАМ, МЕАМ.

© Волегов Павел Сергеевич - кандидат физико-математических наук, доцент, e-mail: [email protected] Герасимов Роман Михайлович - студент, e-mail: [email protected] Давлятшин Роман Позолович - студент, e-mail: [email protected]

Pavel S. Volegov - CSc in Physical and Mathematical Sciences, Associate Professor, e-mail: [email protected] Roman M. Gerasimov - Student, e-mail: [email protected] Roman P. Davlyatshin - Student, e-mail: [email protected]

MODELS OF MOLECULAR DYNAMICS: A REVIEW OF EAM POTENTIALS. PART 1: POTENTIALS FOR SINGLE-COMPONENT SYSTEMS

P.S. Volegov, R.M. Gerasimov, R.P. Davlyatshin

Perm National Research Polytechnic University, Perm, Russian Federation

ARTICLE INFO ABSTRACT

This article is the first part of the review devoted to the description of modern approaches to the construction of interatomic potentials based on the embedded atom method (EAM) for various one-component materials. We consider the general ideology of the molecular (atomic) dynamics method and outline the shortcomings of the approach to modeling atomic systems using the classical pair-interaction potentials. The basic ideas and relations of the embedded atom method are presented, as well as various modifications of the method, including MEAM. The algorithm for numerical implementation of EAM potentials in calculations of molecular dynamic systems is considered, as well as the ways to improve the efficiency of the computational algorithms, including those related to the introduction of the cut-off radius and the justification of its choice for various materials. The review focuses on the methods of construction and types of embedding functions for a variety of materials, as well as the physical and mechanical properties of materials, which are well described in some kind of potential. A brief review of the methods for obtaining potential parameters is given, e.g. by using the quantum mechanical calculations (often referred to as "first principles" or "ab initio") or using experimentally determined properties (diffraction, X-ray methods or electron microscopy), as well as a list of physical and mechanical properties that can be determined for a material using molecular dynamics methods with the application of modern potentials constructed by EAM. The authors have considered in detail the works which offer those potentials that (at present) most accurately describe the properties of such "difficult" (for forward modeling) materials, such as beryllium, iron, tungsten, niobium, titanium, uranium and others. The relevance of studying the structure and properties of these materials is largely due to the prospects of its application in materials for nuclear power engineering, aerospace industry, and biomedicine. Also, this review is devoted to the construction of EAM potentials capable for describing the structure and properties of a number of metals (lithium, nickel, copper, aluminum, etc.), which are considered both in the melt state and after crystallization.

© PNRPU

Введение

Состояние и эволюция внутренней структуры твердых тел играет ключевую роль с точки зрения их физико-механических характеристик. Для изучения закономерностей изменения внутренней структуры часто применяются различные дискретные модели материалов, к числу наиболее распространенных из которых относится метод молекулярной динамики (МД). В частности, в настоящее время метод молекулярной динамки широко применяется для изучения [63]:

• полимеров: цепочек, колец и разветвленных молекул, процессов релаксации и явлений переноса;

• биомолекул: структуры и динамики белков, стыковки молекул, мембран, мицеллы, сворачивания белков;

• динамики жидкости: ламинарного потока, пограничных слоев, реологии неньютоновских жидкостей, неустойчивых потоков;

• твердых тел, в том числе кристаллов: образования и миграции дефектов, разрушения, границ зерен, структурных превращений, влияния радиации, упругих и пластических свойств, трения, изучения поверхностных слоев, получения покрытий методом послойной осадки из газовой среды.

Received: 30 August 2017 Accepted: 10 December 2017 Published: 29 December 2017

Keywords:

molecular dynamics, interaction potential, embedded atom method, EAM, MEAM.

Кратко напомним основную идею и соотношения МД. Основная идея метода молекулярной (атомарной) динамики заключается в представлении тела в виде дискретных взаимодействующих частиц (атомов) и (численном) решении классических уравнений движения Ньютона, описывающих это взаимодействие. При этом уравнение движения частицы с радиусом-вектором г; имеет следующий вид (см., например, [8] или [3]):

N N

mr = Z ф(r)+ Z Т(r,r)+v(r,r),

j=1, i*i 1=1, i*i (1)

ф (ry ) = "VU (j),

где r. - радиус-вектор, соединяющий частицы с номерами i и j, r. = ri - rj; Ф (r) - консервативная составляющая взаимодействия частиц, описывающая силы взаимодействия, как правило, имеющие электромагнитный характер; U (r) - потенциал межатомного

взаимодействия; Т (r, r) описывает неконсервативную (диссипативную) составляющую взаимодействия частиц и часто вводится как некоторая фиктивная сила для описания «внутреннего трения», а у (r, r) - внешние силовые поля (которые иногда также разделяют на консервативную и неконсервативную составляющие). Так, например, для описания процесса осаждения частиц атомных размеров из газовой атмосферы часто необходимо учитывать силу тяжести, а для описания трения в области зоны контакта поверхностей - вводить «силу трения». Далее систему уравнений (1) необходимо снабдить совокупностью начальных (и в общем случае - граничных) условий и решить численно с применением той или иной схемы интегрирования.

Очевидно, что ключевым вопросом метода является выбор потенциала межатомного взаимодействия U (r), поскольку именно потенциал полностью описывает реальное физическое взаимодействие частиц друг с другом (в отличие от двух других слагаемых правой части (1), которые используются для учета внешнего воздействия или для описания объемной диссипации). В настоящее время существует множество потенциалов межатомного взаимодействия, которые описывают те или иные свойства исследуемых материалов, необходимые при исследовании различных процессов (например, упругие характеристики кристаллической решетки, коэффициенты теплового расширения, свойства точечных и линейных дефектов, свойства фазовых переходов, в том числе твердотельных, и многие другие).

Исторически сложилось, что идея и методологические вопросы метода молекулярной динамики: корректное определение граничных и начальных условий, решение вопросов о представительности расчетной области, применение устойчивых численных схем интегрирования уравнений движения и связанный с этим вопрос о выборе шага интегрирования, вопросы поддержания микроканонических ансамблей, например, изотермо-изобарических (NPT) или изотермо-изохорических (NVT) и других, - развивались с некоторым запозданием по сравнению с исследованиями, в которых рассматривались частные по отношению к МД вопросы описания характера взаимодействия атомов в различных системах (или, что то же самое, предлагались различные потенциалы межатомного взаимодействия). Так, если ключевыми работами, заложившими основы методов МД, являлись работы, опубликованные в середине XX века (а частности, статьи B.J. Alder и T.E. Waingwright [7] и J.B. Gibson, A. N. Goland, M.Milgram и G.H. Vineyard [41]), то частные вопросы описания взаимодействия частиц поднимались еще за 20-30 лет до этого.

Например, один из наиболее «популярных» среди исследователей потенциал Леннард-Джонса был предложен последним еще в 1924 г. [47] для описания парного взаимодействия сферических неполярных молекул, потенциалы типа Морзе, первоначально предложенные как приближение для энергии двухатомной молекулы - в конце 20-х гг. XX в. [60], а потенциал Штокмайера для взаимодействия молекул, обладающих постоянным дипольным моментом - в 1941 г. [68]. В последние же десятилетия ситуация существенным образом меняется: при уже достаточно «устойчивой» теоретической и вычислительной основе методов МД вновь встают вопросы, связанные с построением и идентификацией потенциалов межатомного взаимодействия. Уже давно известно, что использование классических парных потенциалов (Леннард-Джонса, Морзе, Ми и др.) дает не вполне корректные результаты при моделировании металлов и некоторых других материалов по причине учета лишь электромагнитных сил парного взаимодействия: такие потенциалы не позволяют получить устойчивую конфигурацию для кристаллических решеток, не относящихся к плотноупакованным; например, ОЦК или простая гексагональная решетки оказываются неустойчивыми при использовании большинства парных потенциалов. Начиная с 80-х гг. XX в. развиваются подходы к моделированию структуры и свойств металлов и сплавов с использованием потенциалов, основанных на так называемом методе погруженного атома (ЕАМ), в котором в отличие от парного потенциала учитываются не только взаимодействия между атомами, а еще и взаимодействия электронных облаков с атомами и, наоборот, атомов с облаками. Следует отметить, что такой подход к описанию энергии атомной системы оказался чрезвычайно удачным при моделировании внутренней структуры и физико-механических свойств металлов и полупроводников (в широком смысле - всех кристаллических твердых тел, а в некоторых случаях - и жидких и/или аморфных сред), что подтверждается огромным количеством исследований и публикаций по этому направлению. Цель настоящего обзора (не претендующего на всеобъемлющий охват всех существующих подходов к построению и анализу потенциалов, построенных на методологии ЕАМ) в том, чтобы дать читателю общее представление о текущем состоянии и тенденциях развития современной (последних 10-15 лет) молекулярной динамики в части разработки и верификации ЕАМ-потенциалов для различных материалов. Авторы обзора постарались, с одной стороны, обратить внимание на «классические» подходы к построению современных ЕАМ-потенциалов, а также выделить некоторые специфические работы, связанные с построением потенциалов и исследованиями свойств материалов, имеющих важное значение для современной техники (материалы для атомной энергетики, аэрокосмической промышленности, биомедицины и т.д.). Обзор состоит из двух логически связанных частей: вначале внимание обращено к работам, в которых разработаны методологические основы метода, в том числе решение важного вопроса о подходах к идентификации параметров потенциала, и в которых рассматриваются одноком-понентные атомные системы и материалы; вторая часть обзора посвящена современному состоянию МД в части описания структуры, свойств и поведения многокомпонентных (бинарных и тернарных) систем.

При этом важно понимать, что метод молекулярной динамики имеет свои, часто довольно существенные, ограничения. В силу того, что движение частиц описывается лишь с помощью второго закона Ньютона, МД неприменим для описания материалов, в которых квантовые эффекты играют основную роль во взаимодействии между частицами (что характерно для процессов, происходящих при экстремально высоких температурах, на сверхмалых расстояниях между атомами или при чрезвычайно малых характер-

ных временах процессов) [61, 72]. Второе существенное ограничение МД связано с невозможностью моделировать реальные технологические процессы из-за огромной вычислительной сложности метода, так как сегодняшний уровень развития вычислительной техники позволяет рассматривать системы размером до кубического микрометра [2]. Таким образом, в настоящее время с применением МД невозможно решать реальные технологические задачи, поскольку применимость метода даже с учетом возможностей суперкомпьютерной техники и технологий высокопроизводительных вычислений в настоящее время ограничена системами не более чем из 109-1010 частиц.

1. Обзор потенциалов для однокомпонентных систем

Первые работы, предвосхитившие появление модели погруженного атома (ЕАМ -embedded atom method), относятся к 1983 году и связаны с именами M.S. Daw и M.I. Baskes [30]. Основными предпосылками появления модели ЕАМ являлись различные теоретические работы в области фундаментальных исследований, направленных на изучение энергии и структуры примесей, поверхностей раздела (так называемых «интерфейсов») и других дефектов в кристаллических структурах многих материалов. Зачастую такие теории имели существенные практические ограничения из-за принятых при их построении «сильных» предположений. Так, например, методы теории групп [40, 66] и кластерные методы [56], как правило, ограничены размерами и требованием периодичности рассматриваемых множеств. Как в теоретико-групповом, так и в кластерном подходах полная энергия системы представляет собой сумму многих одноэлектронных энергий. Альтернативой таким моделям выступали парные потенциалы межатомного взаимодействия, построение («конструирование») которых проводилось из фундаментальных или эмпирических соображений. Однако использование парных потенциалов оказалось успешным лишь для относительно небольшого ряда задач, связанных в первую очередь с описанием инертных и химически неактивных элементов (например, He). Применение парных потенциалов для описания металлов с ОЦК или гексагональной структурой показало себя не вполне корректным: так, при термостатировании происходит самопроизвольное перестроение таких решеток в ГЦК-структуру. Кроме того, часто используемые в МД для описания однородного деформирования образца периодические граничные условия (ГУ) затрудняют моделирование процессов теплового расширения и накладывают ограничения на размеры области моделирования по сравнению с характерными масштабами дефектов, расположенных в этой области.

Основываясь на идеях квантово-механической теории функционала плотности (DFT -density functional theory), а точнее, теории квазиатома [69] и модели эффективной среды [62], M.S. Daw и M.I. Baskes в работе [29], опубликованной в 1984 г., предложили разделять полную потенциальную энергию системы взаимодействующих атомов на энергию парного взаимодействия и так называемую энергию «погружения». Предполагается, что каждый атом локально окружен электронным газом и, таким образом, испытывает влияние от соседних электронных облаков. Расчет электронных плотностей осуществляется согласно теории Хартри-Фока [5] (более подробные описания расчетов даны в работах [25, 53]). Энергия погружения определяется как избыточная энергия, которой обладает атом в электронном газе однородной плотности по сравнению с потенциальной энергией взаимодействия атома, отделенного от этого электронного газа. Полная энергия системы, вычисленная с помощью ЕАМ-потенциала, имеет вид

ЕШа1 I Е, , I

Е = Рг (р,)+ 2 (),

2 1,1

р| = !>( г. ),

1

где Е1 - потенциальная энергия отдельного атома во взаимодействии с окружающей его средой; ¥, - функция погружения; р, - электронная плотность вблизи атома I; ф - вклад

плотности заряда электронов от .-го атома в месте расположения 1-го атома; ф. - энергия

парного взаимодействия атомов I и. (включая их электронные оболочки, но без учета «потерянных» электронов внешних орбиталей, но не учитывая электронный газ); г. - расстояние между атомами I и .. Приведем описание имплементации подхода ЕАМ к уравнениям движения атома (1), используемых в молекулярной динамике [3]. Первое слагаемое в правой части (1) является важнейшим фактором движения частиц и в случае использования ЕАМ-подхода определяется следующим образом:

Ф (г ) = -^,ЕШа, = Е

-I

1

¥ (р,М¥ (р.)+!фу ()

1

1

Г ^ (Р) .¿¥ (р) > й(г) й ф.(г)

ф V ф р=й р=Р. ) йг йг г=ги г=ги к

где г.- расстояние между атомами , и..

Результаты, полученные при использовании такого подхода к построению потенциала, позволили с высокой точностью описать широкий ряд физико-механических свойств (например, упругие постоянные, свойства точечных, поверхностных дефектов, диффузионные свойства и др.) исследуемых материалов, что, в свою очередь, вызвало большой интерес со стороны исследователей и способствовало формированию нового направления в молекулярной динамике. Стоит отметить, что в литературе часто при использовании потенциала ЕАМ говорят о так называемом потенциале Етшв-БтсЫг (ББ) [36], который является частным случаем семейства потенциалов, построенных с применением методологии ЕАМ. Потенциалы типа ББ впервые дали «хорошее» согласование с экспериментальными данными для материалов с ОЦК кристаллической структурой (например, Бе, V, ЫЬ, Та, Мо, ') и в каком-то смысле явились предтечей ЕАМ-подхода.

Несмотря на то, что первые ЕАМ-потенциалы [37, 38] удовлетворительно описывали целый ряд свойств различных материалов (что не удавалось сделать с помощью парных потенциалов межатомного взаимодействия), было обнаружено, что они достаточно хорошо описывают лишь свойства материалов с ГЦК-структурой [10]. Упругие константы, полученные в ходе расчетов с использованием потенциалов ЕАМ для материалов с направленными связями (например, для кремния Б1), плохо согласовывались с экспериментальными значениями, что связано с тем, что формулировка модели ЕАМ предполагает сферическое усреднение плотности электронных облаков, вследствие чего теряется учет направленности связей. Поэтому практически сразу после опубликования пионерской ра-

боты в 1984 г. M.I. Baskes и др. продолжили попытки создания и внедрения более точной и универсальной модели, способной описать многие свойства исследуемых материалов, и уже в 1987 г. опубликовали работу [10]. Авторы предложили так называемую модифицированную модель погруженного атома (МЕАМ), которая основывалась на ранее разработанном ЕАМ-подходе. Основной новизной предложенной модели являлось включение различия между ковалентными и металлическими связями за счет явного учета угловых взаимодействий в функции электронной плотности. Функция погружения в рамках формализма MEAM имеет вид

F (р) = AEc ^ln Р

—0

где А - материальный параметр; Ес - энергия связи; р - электронная плотность структуры с минимальной потенциальной энергий (то есть атомной конфигурации, в которой атомы находятся в положениях устойчивого равновесия); р - суммарная электронная плотность, включающая парциальные вклады сферически-симметричной электронной плотности р(0) и плотности р(1), р(2), р(3), учитывающие угловые вклады 8-, р-, й-,

/-орбиталей соседних атомов [44]. Частичные электронные плотности объединяются с помощью вспомогательной функции

Г i =Z f

(k)

k=1

Р

(k)

Р;

(0)

где t(k) - материальные параметры, описывающие «несимметрию» (анизотропии) оболочек по трем направлениям. Полная электронная плотность вокруг частицы i в таком случае принимает вид

Pi =P(0)G (Гг),

где функция G (Г) выбирается исходя из особенностей исследуемого материала. Стоит отметить, что в формализме МЕАМ конкретная форма функции парного взаимодействия ф(г) не задается, а значение функции парного взаимодействия получается как разность

полной энергии и функции погружения (эта функция в «контрольном» соотношении будет зависеть только от расстояния, так как электронные плотности при расчете являются функциями от расстояния). Полная энергия, вычисленная в расчете на один атом, получается из уравнения, предложенного J.H. Rose и др. [64],

EW (r ) = -Ec (1 + a*) exp [-a* ],

a

= (r / r -1)

/ ч1/2

v Ec j

где Ес - энергия атомной связи; г0 - расстояние до первого соседа в эталонной структуре

(радиус первой координационной сферы); В - (изотропный) объемный модуль упругости; О - равновесный объем элементарной ячейки. Таким образом, парный потенциал имеет вид

ъ ( )_Г 2Е: )-Р (рг ))-р (р- ()

V

где 7. - число ближайших соседей атомов г и /.

Кроме учета угловых взаимодействий, при построении потенциала МЕАМ важным остается вопрос о так называемом радиусе обрезания гсШф (если расстояние между двумя атомами превышает гсиф, то считается, что взаимодействием между парой атомов можно пренебречь). Традиционно при использовании модели типа ЕАМ диапазон взаимодействия ограничивается за счет плавного усечения функций вблизи точки отсечения при помощи различного рода экспоненциальных или гиперболических зависимостей. Обоснованием такой процедуры является то, что силами взаимодействия, возникающими между атомами на больших расстояниях, можно пренебречь. Однако такой способ не является универсальным, так как на практике чаще всего спад на больших расстояниях не происходит достаточно быстро, поэтому потенциал взаимодействия необходимо искусственно продолжать. С.Ь. Би и др. [39] показали, что при описании релаксации атомной системы важны только взаимодействия с атомами, попадающими в первую координационную сферу. Это относится и к материалам с ОЦК-структурой, несмотря на то, что отношение расстояний между второй и первой координационными сферами составляет лишь 1,15. В работе [12] предложена схема, накладывающая дополнительные ограничения на функции взаимодействия. Обоснование такой схемы состоит в том, что атом, находясь между другими атомами, способен «экранировать» взаимодействие между этими «внешними» атомами. Таким образом, расчет взаимодействия происходит только между ближайшими соседями, а взаимодействия с более «отдаленными» атомами происходит путем переноса ближайших взаимодействий, отмасштабированных с помощью так называемого скринин-гового фактора, имеющего вид [59]

^ = П 3

Г

/ *г,к

/с

с — с

¡к/ тт,1к/'

с — с

у тах,¡к/ тт, ¡к/ J

2 2,2 2 4

г г + г г — г

с _ 1 I 2 _ Ь ¡к ¡1 ]к ¡1 -1

(( — г/к )2

г/ —(к — г

фс (х)_

0, х < О, [1 — (1 — х)4]2, О <х < 1,

1, х >

где г - номер рассматриваемого атома; / - один из ближайших соседей атома г; к - экранируемый атом; сгк- - скрининговые параметры; г- - расстояния от атома г до атома/ Учет скринингового фактора происходит путем его умножения на функцию парного взаимодействия, а также функцию электронной плотности.

В рамках рассмотрения подхода МЕАМ к построению потенциалов межатомного взаимодействия стоит отметить работу 2. Bangwei и др. [9], в которой предложен более простой аналитический МЕАМ-потенциал, являющийся модификацией потенциала, предложенного

M.I. Baskes [10]. В [10] электронная плотность каждого атома определяется сферическим усреднением и аппроксимируется суммой плотностей соседних атомов. Однако в действительности распределение электронов внешней оболочки отклоняется от сферической симметрии, следовательно, линейная аппроксимация электронной плотности является не вполне корректной. Предложенная в работе [9] модификация потенциала отличается от [10] тем, что электронная плотность рассматривается способом, аналогичным ЕАМ, однако к общему выражению энергии добавляется энергетический член M(P), который используется для выражения разности между фактической полной энергией системы атомов и рассчитанной из исходной ЕАМ-модели с использованием суперпозиции сферических электронных плотностей. При таком подходе полная энергия системы выражается следующим образом:

= £ F (pi ) + 2 ЕЕ Ф; ( ) + £ Mi (Pi ) ,

i 2 i i, i*i i

p.-=ЕФ( ri),

i

P=£Ф2 (r,.).

i

Функция M(P) не имеет универсальной формы, то есть для аппроксимации экспериментальных данных подбирается эмпирически. Результаты, полученные с использованием предложенной модели, хорошо согласуются с входными экспериментальными данными: постоянной решетки, энергией связи, модулем сдвига, коэффициентом анизотропии, а также позволяют достаточно точно определять энергию вакансии. Кроме того, упругие константы, рассчитанные с помощью модели для различных ОЦК-материалов, в точности совпадают с экспериментальными значениями.

Несмотря на стремление большинства авторов при построении ЕАМ-потенциалов описать максимальное количество физико-механических и геометрических свойств, которые не участвуют в верификации параметров, как правило, сделать это удается лишь для некоторой вполне определенной заранее группы свойств (например, для упругих свойств, характеристик фазового перехода, диффузионных свойств и др.). Достигается это путем создания набора «эталонных» параметров, полученных в ходе экспериментов или же в квантово-механических расчетах [5]. Следует отметить, что второй подход в литературе часто называют основанным на «первых принципах» (ab initio), однако использование этого термина представляется не вполне корректным, поскольку в любом случае при решении соответствующих уравнений исследователи вводят те или иные предположения о характере решения, поэтому далее при упоминании такого подхода более уместным представляется использование термина «квантово-механических расчеты». Чаще всего таким набором параметров выступают: парные корреляционные функции (PCF - pair correlation function) [22]; структурный фактор (TSF - total structure factor) [33] - характеристика рассеивания рентгеновского излучения; свойства переноса; характеристики фазового перехода; параметры решетки; силовые характеристики (полученные, как правило, с использованием квантово-механических расчетов) и др. В зависимости от выбранного набора параметров для верификации параметров потенциала исследователь получает возможность построить потенциал таким образом, чтобы при моделировании с помощью полученного потенциала (например, в расчетах молекулярной динамики) быть уверенным в точном воспроизведении нужных ему свойств.

В качестве примера использования такого подхода приведем работы, посвященные построению потенциалов для бериллия Ве. Основная особенность Ве заключается в особенностях его строения: межатомные связи являются направленными, поэтому работы, посвященные построению межатомных потенциалов для Ве, как правило, использовали потенциалы МЕАМ для включения эффекта угловых сил: например, в работе [13] авторы предложили МЕАМ-потенциал, который позволял рассчитать параметры решетки и энергии связи, хорошо согласующиеся с экспериментальными данными, однако при этом потенциал оказался непригодным для учета зависимости этих свойств от температуры, из чего был сделан вывод, что только лишь описания направленности связей недостаточно для корректного описания термодинамических свойств. В работе [6] предложен ЕАМ-потенциал для Ве, имеющего гексагональную структуру кристаллической решетки. Следует отметить, что предложенный потенциал является самым точным ЕАМ-потенциа-лом для Ве на настоящий момент, хотя в явном виде и не учитывает направления атомных связей. В качестве функции парного взаимодействия выбран потенциал, основанный на парном потенциале Морзе и имеющий вид [60]

Фм,. (r) = D (e-2a(r->- 2e-a(r-rc),

где Бс - энергия связи; гс - длина равновесной связи; а - вторая производная функции

потенциальной энергии в точке минимума. Результирующий парный потенциал имеет следующий вид:

Фе (Г) = Фмаг*е (Г)-ф(ГШа// )

cutoff

m

1 -

f \ r

r

V cutoff

d ф

dr

cutoff

где гС1йоф- радиус обрезания потенциала; т - материальная константа. Такая форма парного потенциала обеспечивает непрерывность функции парного взаимодействия и ее производных при различных значениях радиуса обрезания гсиф Функция электронной плотности выбрана в классическом экспоненциальном виде [46]:

ф(г)= Лв~В(г гс),

где Л, В- параметры потенциала, а в качестве функции погружения использовалась функция Джонсона [45]

F (р) = F

1 - ln

- F

1р0 у 1р0 у 1ро у

где и - константы; р0, Р, у - параметры потенциала.

Качество подгонки контролировалось с помощью взвешенной суммы квадратов разностей:

G ^ ' Wi fi(calc) fi(obs)

i=1

где wi - веса, определяемые исследователем; /(са1с)и ф(оЫ) тальные значения наблюдаемых параметров.

расчетные и эксперимен-

Подбор параметров потенциала проводился с помощью метода Ньютона-Рафсона. В качестве экспериментальных данных использовались значения упругих постоянных и энергия связи. Кроме того, в экспериментальные данные были внесены значения энергии образования вакансии, полученной из расчетов с применением моделей квантовой динамики. Полученный потенциал сравнивался с рядом предшествующих ему потенциалов [13, 48] и др. Потенциал правильно предсказывает стабильную конфигурацию кристаллической решетки и может быть использован для моделирования структуры материала, а также повреждения радиоактивным излучением (что является актуальной задачей в связи с частым использованием бериллия для производства радиационных окон рентгеновских трубок). Потенциал хорошо подходит для изучения нелинейного расширения, а также для описания свойств точечных и линейных дефектов, однако предсказываемая температура плавления для потенциала составляет 1200-1300 К, что не согласуется с экспериментальными данными (1560К). В силу того что большинство других потенциалов для Be строились методом МЕАМ, предложенный потенциал имеет преимущество в вычислительной эффективности.

В общем случае параметры ЕАМ-потенциала для конкретного материала идентифицируют один раз, а затем используют для исследования структуры и свойств материала как в твердом, так и в жидком состоянии. Однако такой подход не совсем корректен и вызывает ряд вопросов, связанных в первую очередь с адекватным описанием динамических свойств и свойств переноса, так как в большинстве случаев ЕАМ-потенциалы дают хорошее согласование с экспериментом в сравнительно небольшом диапазоне температур. В работе [27] авторы рассмотрели вопросы переносимости эффективных ЕАМ-потенциалов твердой фазы для изучения свойств жидкой фазы Cu. При построении потенциала использованы вид функции погружения, электронной плотности, а также парного взаимодействия, предложенные в работах Ю. Мишина и др. [58] и M. Doyama и др. [32], после их репараметризации (первоначально параметры этих потенциалов были предложены в работе [28]). Переносные свойства (самодиффузия) предложенных потенциалов сравнивались с помощью схемы, предложенной L.E. Gonzalez и др. [42]. Сравнивая потенциалы, авторы [27] пришли к выводу, что потенциал, предложенный M. Doyama и др., лучше согласуется с экспериментальными данными и его можно использовать для получения представления о понимании микроскопического поведения некоторых металлов и сплавов.

Работы Д.К. Белащенко и др. [15-18] представляют собой серию статей, в которых одним и тем же способом осуществляется построение ЕАМ-потенциалов для жидких (высокотемпературных) Li, Na, Cu, U и др. Для построения потенциала использовались данные, полученные с помощью метода дифракционной микроскопии при температуре, близкой к температуре плавления материалов. Согласно экспериментальным данным оценивались критические параметры материалов (температура, давление, плотность). Структура жидких материалов оценивалась с помощью парных корреляционных функций (PCF). Подбор параметров парной части потенциала выполнялся путем построения и идентификации прямой МД-модели по дифракционным данным с использованием алгоритма Шоммерса [65], который был продолжен на область малых расстояний (в случаях высоких давлений) с помощью экспоненциальной формулы Борна-Майера [31]. Подбор параметров осуществлялся с помощью оценки невязки между дифракционной и модельной PCF. Затем построенные потенциалы использовались для моделирования и оценки структуры, диффузионных свойств и температуры плавления исследуемых материалов. Сравнение полученных результатов с экспериментальными данными позволяет утверждать корректность предложенного потенциала.

М.И. Менделев и др. [43, 54] предложили несколько ЕАМ-потенциалов для ОЦК Fe, моделирование которого, как известно, сильно осложнено из-за наличия у него магнитных свойств. Ранее потенциалы для Fe были построены B.-J. Lee и др. [51] методом МЕАМ, однако температура плавления в расчетах получалась равной 2200 К (при эмпирическом значении 1812 К). В [54] на примерах работ предшественников было показано, что использование данных, относящихся к свойствам исключительно идеальных кристаллов, недостаточно для потенциалов типа ЕАМ ввиду их плохого согласования со значением структурного фактора (TSF) жидкости. Дополнить процедуру подбора параметров было предложено двумя способами: включением данных о межатомных силах, полученных из квантово-механических расчетов, или экспериментальных данных, полученных с помощью дифракционной микроскопии. Первый способ подразумевает следующее: генерируется атомная модель для жидкости с правильной (заранее известной) плотностью, затем выполняется квантово-механический расчет результирующих сил на всех атомах, и значения сил используются для определения параметров потенциала (вид потенциала выбирается заранее). Такой метод подбора параметров в литературе принято называть методом «согласования сил» (force-matched method). Процедура обеспечивает рассмотрение ряда атомных геометрий, включая случаи малых атомных расстояний. МД-модель состояла из 100 атомов, что связано с большими вычислительными затратами квантово-механических расчетов. Второй способ подразумевает использование для подбора параметров потенциала данных парной корреляционной функции (PCF), а также параметров ГЦК- и ОЦК-решетки. Преимущество первого подхода заключается в том, что данные дифракции доступны не всегда. Наоборот, преимущество потенциалов, построенных по дифракционным данным, обусловлено тем, что МД-модели лучше согласуются с экспериментальными данными, что связано с тем, что расчеты с использованием моделей квантовой динамики, как правило, выполняются для идеальных кристаллов. Исходя из предложенных в работе данных стоит отметить, что второй способ (то есть по данным дифракционной микроскопии) оказался несколько лучше, так как более точно воспроизводил значения энергии связи, а также параметры парной корреляционной функции. Однако нельзя утверждать, что использование данных, полученных из квантово-механических расчетов, позволяет строить ЕАМ-потенциалы для любых систем без использования дифракционных данных. В качестве примера стоит отметить работу [34], в которой процедура построения потенциала с помощью квантово-механических расчетов для чистого алюминия происходила похожим образом, однако согласие с экспериментальными данными по структурному фактору (TSF) было плохое (несмотря на хорошее согласование с температурой при заданном давлении).

В работе H. Chamati и др. [21] предложен ЕАМ-потенциал для Fe, разработка которого была нацелена на исследование свойств самодиффузии. Эффективная парная функция взаимодействия частиц строилась с помощью суммы двух функций взаимодействия Морзе. Функции погружения аппроксимируются полиномами 4-й степени по экспериментальным точкам. Подбор параметров происходил путем минимизации взвешенного квадрата отклонений значений выбранных свойств исследуемого материала от их целевых (экспериментально известных) значений. Экспериментальный набор параметров включал в себя параметры решетки, некоторые значения потенциальной энергии взаимодействия атомов, объемный модуль упругости, упругие константы, энергию связи ослабленной вакансии, энергию миграции точечных дефектов, а также среднюю поверхностную энергию. Полученный потенциал точно описывает свойства кристаллической решетки, тепловое расширение, а также свойства точечных дефектов для ОЦК и высокотемпературной ГЦК-фазы Fe. Построенный потенциал использовался для моделирования самодиффузии

одиночных атомов в температурном интервале 600-950 К. Температурные зависимости коэффициентов диффузии, полученные в численных экспериментах, согласуются с законом Аррениуса. Построенный потенциал также позволил обнаружить, что коэффициенты диффузии с диагональным механизмом обмена примерно на порядок выше, чем при механизме «прямого скачка» и механизме «недиагонального обмена».

В работе M.-C. Marinica и др. [52] предложены и сопоставлены 4 варианта потенциала для чистого W. Сначала вводится разделение наборов данных, по которым происходит построение потенциалов:

1. экспериментальные данные, под которыми понимаются экспериментально измеряемые свойства твердого тела: параметры решетки, энергия связи и упругие константы С11, С12, С44 (в нотации Фойгта [4]);

2. данные, полученные из квантово-механических моделей (или, точнее, при использовании теории функционала плотности): энергия образования основных точечных дефектов вольфрама (вакансии и междоузельного атома);

3. данные, полученные из квантово-механических моделей: характерное расположение атомов в жидких (случайных) состояниях.

Форма ЕАМ-потенциала соответствовала предложенной в работе М.И. Менделева [54]. Процедура подбора параметров заключалась в вычислении всех параметров трех функций (парной функции, функции электронной плотности и функции погружения) и основывалась на принципе вариационной ассимиляции [50]. Целевая функция, используемая для минимизации отклонения, выглядела следующим образом:

2

J(А) = Е[Yo,s - Yo,s (A)] + £W [Flnitio -Fi (A)]2,

obs i

где А - вектор набора параметров потенциала; Yobs и Yobs (A) - значения набора параметров, полученных экспериментально, и значения, вычисленные с использованием набора параметров А, F1cbinitio и F! (A) - силы, действующие на атом i, вычисленные с помощью квантово-механических расчетов и с помощью полученного потенциала ЕАМ; Wobs и Wi - подгоночные параметры. Варьирование параметрами Wobs и Wi в работе позволило получить четыре ЕАМ-потенциала одного типа, отличающихся значениями параметров. Все эти потенциалы хорошо описывают образование и характеристики винтовых дисклокаций вольфрама, а также возникновение радиационных дефектов, что не удавалось сделать в более ранних исследованиях. Данные, полученные с использованием первого потенциала, не согласовывались с некоторыми данными, полученными с использованием теории функционала плотности (DFT), поэтому в дальнейшем первый потенциал был исключен из использования.

В работе [24] проведено сравнение четырех потенциалов МЕАМ, предложенных Baskes [11], с двумя потенциалами ЕАМ, предложенными J.E. Angelo и др. [14] и J. Cai и др. [20], построенных для жидкого Ni, причем вид функций и методы подбора параметров потенциалов были различны. Подбор параметров потенциала J.E. Angelo происходил по экспериментальным данным, которые включали в себя упругие модули, параметр решетки, энергию образования вакансии, энергию дефекта упаковки и энергию устойчивости ГЦК-фазы к фазам ОЦК Ni, с помощью нелинейного метода наименьших квадратов, разработанного в Sandia National Laboratories. Построение потенциала J. Cai и др. проводилось по такому же набору экспериментальных данных, как и в случае с первым потенциалом, одна-

ко процедура подбора параметров осуществлялась путем минимизации отклонения экспериментальных и модельных данных с помощью симплекс-метода. Потенциалы, предложенные в работе [11], были построены по схеме, описанной в работе [12]. Для этих шести потенциалов (четыре потенциала M.I. Baskes и два потенциала J.E. Angelo и J. Cai) выполнялось МД-моделирование структуры и свойств Ni при низких давлениях и температурах, превышающих температуру плавления, в котором оценивались свойства переноса (коэффициенты диффузии, вязкость и др.), структурного фактора (TSF), значение температуры плавления. В результате было показано, что потенциал, предложенный J.E. Angelo и др., дает худшее согласование по температуре плавления, в то время как другие потенциалы описывали это значение достаточно точно; с другой стороны, потенциал J.E. Angelo и др. хорошо согласуется по экспериментальным значениям энергии связи. Проведенные численные эксперименты с использованием вышеперечисленных потенциалов позволяет судить о том, что коэффициент самодиффузии и вязкость зависят от особенностей структуры жидкой фазы в пределах радиусов второй-третьей координационных сфер.

В работе M.R.Fellinger и др. [35] предложен ЕАМ-потенциал для Nb. Актуальность исследования материала заключается в его свойствах: низкая плотность, высокая температура плавления, биосовместимость делают его привлекательным материалом для использования в сплавах различного рода (например, в медицине создаются прочные и нетоксичные соединения Ti-Nb). Парная функция, функции погружения и электронной плотности были построены в виде кубического сплайна. В работе использовался метод «согласования сил»: полученные в ходе квантово-механических расчетов значения межатомных сил и энергии, вычисленной на один атом, вносились в базу данных и использовались при подборе параметров ЕАМ-потенциала с помощью корректировки относительно отклонения параметров, полученных с помощью потенциала в квантово-механических расчетах и ЕАМ-потенциала [19]. Расширение наборов значений физических величин, полученных из квантово-механических расчетов, путем добавления значений этих же величин при различных температурных состояниях позволяет использовать потенциал для широкого диапазона температур. В целом подбор параметров ЕАМ-потенциалов с использованием метода «согласования сил» является задачей оптимизации в многомерном пространстве. Функции ЕАМ-потенциала параметризуются с помощью кубических сплайнов, а программа POT-FIT [19] оптимизирует сетку сплайн-узлов, используя комбинации имитируемых алгоритмов отжига [71] и сопряженных градиентов [1]. Стоит отметить, что полученный потенциал успешно прошел верификацию на свойствах, не включенных в набор данных для идентификации параметров: потенциал точно описывает структурные и упругие свойства, свойства точеных и линейных дефектов, однако не очень хорошо подходит для исследования ударной волны и радиационного воздействия. Кроме того, в цитируемой работе подчеркивается, что на основе построенного потенциала имеется возможность строить потенциалы и для сплавов ниобия с помощью метода ЕАМ.

В работе [55] предложены три ЕАМ-потенциала для чистого Ti, представляющего большой интерес в аэрокосмической и медицинской отрасли благодаря высокой прочности, низкой плотности и коррозийной стойкости. Несмотря на это, в литературе относительно мало работ, посвященных изучению внутренней структуры Ti, что связано с малым количеством информации (например, о критических значениях различных параметров), а также особенностями структуры, которые обусловлены наличием двух устойчивых фаз Ti: ГПУ (a-фаза) и ОЦК (Р-фаза). Исходные данные, по которым происходило построение и идентификация параметров потенциала, разделены на четыре группы:

1. базовые свойства кристалла (параметр решетки, энергия связи, упругие константы при 0 К);

2. энергия образования точечных дефектов;

3. энергия образования плоских (планарных) дефектов;

4. параметры фазового перехода.

Использование данных первой группы при подборе параметров гарантирует, что ГПУ-структура является более стабильной фазой и позволяет подобрать правильный энергетический «масштаб». Значения параметров второй группы позволяют воспроизводить свойства самодиффузии. Данные об энергии упаковки, относящиеся к третьей группе, обеспечивают корректные результаты при моделировании пластической деформации. Четвертая группа содержала информацию о температуре плавления ß-фазы, температуре фазового перехода a-ß и скрытые теплоты этих переходов. Стоит отметить, что построенные потенциалы были изначально предназначены для описания как жидкой, так и твердой фазы Ti. Для всех трех потенциалов вид и значения параметров функции электронной плотности совпадали, функции погружения и парного взаимодействия имели один и тот же вид в каждом из трех потенциалов, но значения подгоночных параметров были разными, что обусловлено выбором различных весовых коэффициентов при подгонке потенциала к экспериментальным данным. Основная задача при построении первого потенциала заключалась в воспроизведении данных четвертой группы, однако полученный потенциал плохо согласовывался с данными скрытой теплоты и теплового расширения. С использованием предложенного потенциала было показано, что ß-фаза титана более плотная, чем a-фаза. Основная цель разработки второго потенциала заключалась в возможности описания энергии образования точечных дефектов в a-фазе. В работе отмечено, что эта цель может быть достигнута только путем исключения данных о температуре плавления из набора параметров для подгонки параметров потенциала. Как результат второй потенциал плохо описывал температуры плавления ß-фазы, но все же корректно описывал температуру фазового перехода a-ß-фаз. Полученные энергии образования точечных дефектов довольно точно согласуются с данными, полученными из квантово-механических расчетов. Исходя из того, что второй потенциал не давал монотонно возрастающей зависимости параметра решетки от температуры (то есть, по сути, не описывал линейное тепловое расширение), авторами был предложен третий потенциал, который учитывал линейное тепловое расширение, однако хуже согласовывался с данными точечных дефектов в a-фазе. В цитируемой работе показано, что построенные потенциалы пригодны для описания пластической деформации, так как правильное описании энергии дефекта упаковки обеспечивает корректное описание поведения дислокаций.

В работе Д.Е. Смирновой и др. [67] предложен ЕАМ-потенциал для урана U. Сложности, связанные с моделированием урана, обусловлены его радиоактивной природой, а также большим набором фаз: орторомбической (a-фаза), тетрагональной (ß-фаза), высокотемпературной ОЦК (у-фаза) и жидкой фаз. Предшествующие работы по построению потенциалов для урана описывали, как правило, относительно небольшой спектр исследуемых свойств и были основаны на данных, вычисленных в квантово-механических расчетах. Стоит отметить, что работа Д.К. Белащенко [16] посвящена построению ЕАМ-потенциала для урана, который дает точное описание свойств жидкого урана, однако не позволяет описывать твердую фазу. Разработка потенциала осуществлялась с помощью метода «согласования сил» [34], предполагающего построение набора данных различных

параметров материала, полученных из первых принципов, а затем подбор параметров потенциала и его оптимизация. При построении потенциала авторы ориентировались на свойства а-, у-фаз, и жидкую фазу урана. Для построения справочного набора параметров использовались наборы эталонных структур, которые были получены в ходе прямого мо-лекулярно-динамического моделирования с использованием потенциала Д.К. Белащенко [16]. Затем, используя квантово-механические соображения, рассчитывались значения межатомных сил, а также значения энергии на атом для каждой конфигурации. Для воспроизведения внутренней структуры урана использовался псевдопотенциал, предложенный в работе [49]. Все функции ЕАМ-потенциала интерполированы кубическими сплайнами. Сплайн-узлы в ходе подбора параметров потенциала корректировались итерационным способом путем минимизации отклонений между значениями, полученными в ходе квантово-механических расчетов, и значениями, вычисленными с помощью ЕАМ-потенциала. После процедуры подбора параметров проводился ряд численных экспериментов для определения тестовых структур и свойств урана. Полученные результаты сравнивались с экспериментальными данными. Если соглашение являлось недостаточным, то производилась модификация базы данных: добавлялись или удалялись отдельные конфигурации, изменялись веса конфигураций на схеме минимизации. В цитируемой работе показано, что построенный потенциал дает хорошее количественное описание структуры и свойств твердого урана. Кроме того, показано, что построенный потенциал может быть использован для описания трансформации а-у урана.

В работе J.M. Winey [73] предложен потенциал для Al. Построение потенциала производилось по данным атомного объема [26] и упругих констант второго [70] (SOEC -second-order elastic constants) и третьего порядка (TOEC-third-order elastic constants), полученных из экспериментов [70] при комнатной температуре. Экспериментальные данные в ходе подбора параметров потенциала использовались после их экстраполяции до температуры T = 0 K с использованием классических термодинамических соотношений (квантовые эффекты не учитывались). Разработка предложенного потенциала была мотивирована тем фактом, что подавляющее большинство работ, связанных с построением потенциалов межатомного взаимодействия типа ЕАМ, не включают в себя оценки предсказания термоупругих свойств для различных термодинамических условий, следовательно, между результатами классического МД-моделирования и экспериментальными данными могут возникать существенные расхождения, ввиду того, что такие процессы, как упругопластическая деформация и фазовые переходы, определяются лежащей в их основе термоупругой реакцией материала. Разработанный потенциал сравнивался с экспериментальными данными [26], а также расчетными данными, полученными с использованием потенциалов предшественников [23, 34, 57]. Было показано, что экстраполяция экспериментальных данных до температуры T = 0 K позволяет точно описать атомный объем при более высоких температурах, что не удавалось сделать с помощью потенциалов, разработанных ранее: потенциалы [23, 34] разрабатывались на основе экспериментальных данных при низких температурах, поэтому классическое описание атомного объема при более высоких температурах является не вполне корректным, поскольку поведение измеренного низкотемпературного атомного объема отклоняется от классического высокотемпературного поведения из-за квантовых эффектов; потенциал [57] давал худшее согласование с атомным объемом ввиду того, что при построении потенциала использовались экспериментальные данные при комнатной температуре без экстраполяции.

В таблице приведена сравнительная характеристика некоторых потенциалов, предложенных в цитируемых выше работах, что позволяет использовать тот или иной потенциал в зависимости от особенностей исследования.

Сравнительная характеристика применяемых EAM-потенциалов Comparative characteristics of different EAM potentials

№ п/п Материал Авторы, ссылка на источник Физико-механические свойства, хорошо описываемые потенциалом Физико-механические свойства, плохо описываемые потенциалом, другие недостатки

1 Be M.I. Baskes [13] Параметры решетки и энергии связи Зависимость свойств от температуры

2 Be A. Agrawal [6] Структура кристаллической решетки, радиационные свойства, тепловое расширение Температура плавления

3 Li, Na, Cu, U Д.К. Белащенко [15-18] Критические температура, давление, плотность, фазовый переход первого рода, диффузионные свойства Требуется большое количество подгоночных параметров

4 Fe (ОЦК) М.И. Менделев [43, 54] Структура кристаллической решетки, фазовый переход первого рода, магнитные свойства Требуется большое количество подгоночных параметров

5 Fe (ОЦК) B.-J. Lee [51] Структура кристаллической решетки, магнитные свойства Температура плавления

6 Fe H. Chamati [21] Структура кристаллической решетки, свойства самодиффузии, тепловое расширение Требуется большое количество подгоночных параметров

7 W M.-C. Marinica [52] Образование и характеристики винтовых дислокаций, возникновение радиационных дефектов Структура кристаллической решетки

8 Ni M.I. Baskes [11] Фазовый переход первого рода, температура плавления Свойства переноса (коэффициенты диффузии, вязкость и др.)

9 Ni JE. Angelo [14] Фазовый переход первого рода, энергия связи Температура плавления

10 Ni J. Cai [20] Фазовый переход первого рода, температура плавления Энергия связи

11 Nb MR. Fellinger [35] Структурные и упругие свойства, свойства точеных и линейных дефектов Ударно-волновое воздействие, радиационное воздействие

12 Ti M.I. Mendelev [55] Линейное тепловое расширение, энергия дефекта упаковки Фазовый переход второго рода, свойства точечных дефектов в а-фазе

13 U Д.К. Белащенко [16] Свойства в состоянии расплава Свойства в твердом состоянии

14 U Д.Е. Смирнова [67] Структура и свойства твердого урана, фазовый переход а-у Свойства Р-фазы

15 Al J.M. Winey [73] Атомный объем при ненулевой температуре, тепловое расширение Фазовый переход первого рода

Заключение

В статье приведен краткий обзор работ последних 10-15 лет, в которых предлагаются различные подходы к построению, идентификации параметров и проверке (верификации) потенциалов межатомного взаимодействия для однокомпонентных атомных систем. В настоящее время подходы к построению потенциалов, основанные на методологии погруженного атома, все более широко используются в молекулярно-динамических моделей, а последние, в свою очередь, позволяют все более точно описывать как разнообразные свойства материала, так и процессы, происходящие на уровне атомарной структуры. Используя подходы молекулярной динамики, исследователи получают количественные оценки не только параметров кристаллической решетки, коэффициентов теплового расширения или упругих характеристик металлов, но и такие важные для понимания происходящих в материале процессов характеристики, как свойства точечных и линейных дефектов, энергию дефекта упаковки и межатомных границ, характеристики фазовых переходов первого и второго рода, магнитные свойства ферромагнетиков.

Благодарности

Работа выполнена при финансовой поддержке гранта Президента РФ № МК-1298.2017.1. Acknowledgments

This work was carried out with the financial support of the Grant of the President of the Russian Federation Nr. МК-1298.2017.1.

Библиографический список

1. Гилл Ф., Мюррей У., Райт М. Практическая оптимизация: пер. с англ. - М.: Мир, 1985. -509 с.

2. Кривцов А.М. Теоретическая механика. Упругие свойства одноатомных и двухатомных кристаллов: учеб. пособие. - СПб.: Изд-во Политехи. ун-та, 2009. - 126 с.

3. Кривцов А.М., Кривцова Н.В. Метод частиц и его использование в механике деформируемого твердого тела // Дальневосточный математический журнал. - 2002. - Т. 3, № 2. - С. 254-276.

4. Мейз Дж. Теория и задачи механики сплошных сред: учеб. - М.: Мир, 1974. - 319 с.

5. Хартри Д. Расчеты атомных структур. - М.: Изд-во Иностранной литературы, 1960. - 272 с.

6. An embedded atom method potential of beryllium / A. Agrawal, R. Mishra, L. Ward, K.M. Flores, W. Windl // Modelling and Simulation in Materials Science and Engineering. - 2013. - Vol. 21. -No 8. DOI: 10.1088/0965-0393/21/8/085001

7. Alder B.J., Waingwright T.E. Phase Transition for a Hard Sphere System // J. Chem. Phys. -1957. - Vol. 27. - No. 5. - P. 1208-1209.

8. Allen M.P., Tildesley A.K. Computer Simulation of Liquids. - Oxford: Clarendon Press, 1987. -385 p.

9. An analytic MEAM model for all BCC transition metals / Z. Bangwei, O. Yifang, L. Shuzhi, J. Zhanpeng // Physica B: Condensed Matter. - 1999. - Vol. 262. - No. 3-4. - P. 218-225. DOI: 10.1016/S0921-4526(98)01156-9

10. Baskes M.I. Application of the Embedded-Atom Method to Covalent Materials: A Semiempirical Potential for Silicon // Physical Review Letters. - 1987. - Vol. 59. - No. 23. - P. 2666-2669. DOI: 10.1103/PhysRevLett.59.2666

11. Baskes M.I. Determination of modified embedded atom method parameters for nickel // Materials Chemistry and Physics. - 1997. - Vol. 50. - No. 2. - P. 152-158. DOI: 10.1016/S0254-0584(97)80252-0

12. Baskes M.I. Modified embedded-atom potentials for cubic materials and impurities // Physical Review B. - 1992. - Vol. 46. - No. 5. - P. 2727-2742. DOI: 10.1103/PhysRevB.46.2727

13. Baskes M.I., Johnson R.A. Modified embedded atom potentials for HCP metals // Modelling and Simulation in Materials Science and Engineering. - 1994. - Vol. 2. - No. 1. - P. 147-163. DOI: 10.1088/0965-0393/2/1/011

14. Trapping of hydrogen to lattice defects in nickel / Baskes M.I., Sha X., Angelo J.E., Moody N.R. // Modelling and Simulation in Materials Science and Engineering. - 1995. - Vol. 5. - No. 6. -P. 651-652. DOI: 10.1088/0965-0393/3/3/001

15. Belashchenko D.K. Embedded Atom Method Potentials for Liquid Copper and Silver // Inorganic Materials. - 2012. - Vol. 48. - No. 9. - P. 940-947. DOI: 10.1134/S002016851209004X

16. Belashchenko D.K., Smirnova D.E., Ostrovski O.I. Molecular-Dynamic Simulation of the Thermophysical Properties of Liquid Uranium // High Temperature. - 2010. - Vol. 48. - No. 3. -P. 363-375. DOI: 10.1134/S0018151X10030107

17. Belashchenko D.K., Zhuravlev Yu.V. Application of the embedded-atom method to liquid copper // Inorganic Materials. - 2008. - Vol. 44. - No. 9. - P. 939-945. DOI: 10.1134/S0020168508090082

18. Belashchenko D.K., Ostrovskii O.I. Application of the embedded atom model to liquid metals: Liquid lithium // High Temperature. - 2009. - Vol. 47. - No. 2. - P. 231-237. DOI: 10.1134/S0018151X09020102

19. Brommer P., Gahler F. Effective potentials for quasicrystals from ab-initio data // Philosophical Magazine. - 2006. - Vol. 86. - No. 6-8. - P. 753-758. DOI: 10.1080/14786430500333349

20. Cai J., Ye. Y.Y. Simple analytical embedded-atom- potential model including a long-range force for fcc metals and their alloys // Physical Review B. - 1996. - Vol. 54. - No. 12. - P. 8398-8410. DOI: 10.1103/PhysRevB.54.8398

21. Embedded-atom potential for Fe and its application to self-diffusion on Fe(100) / H. Chamati, N.I. Papanicolaou, Y. Mishin, D.A. Papaconstantopoulos // Surface Science. - 2006. - Vol. 600. - No. 9. -P. 1793-1803. DOI: 10.1016/j.susc.2006.02.010

22. Chandler D. Introduction to Modern Statistical Mechanics. - New York, Oxford University Press, 1987. - 274 p.

23. Chantasiriwan S., Milstein F. Embedded-atom models of 12 cubic metals incorporating second-and third-order elastic-moduli data // Physical Review B. - 1998. - Vol. 58. - No. 10. - P. 5996-6005. DOI: 10.1103/PhysRevB.58.5996

24. Cherne F.J., Baskes M.I., Deymier P.A. Properties of liquid nickel: A critical comparison of EAM and MEAM calculations // Physical Review B. - 2002. - Vol. 65. - No. 2. - 024209. DOI:10.1103/PhysRevB.65.024209

25. Clementi E., Roetti C. Roothaan-Hartree-Fock atomic wavefonctions: Basis functions and their coefficients for ground and certain excited states of neutral and ionized atoms, Z<54 // Atomic Data and Nuclear Data Tables. - 1974. - Vol. 13. - No. 3-4. - P. 177-478. DOI: 10.1016/S0092-640X(74)80016-1

26. Copper A.S. Precise Lattice Constants of Germanium, Aluminum, Gallium Arsenide, Uranium, Sulphur, Quartz and Sapphire // Acta Crystallographica. - 1962. - Vol. 15. - P. 578-582. DOI: 10.1107/S0365110X62001474

27. Dalgic S. The EAM based effective pair potentials for the dynamic properties of liquid Cu // Materials Chemistry and Physics. - 2007. - Vol. 39. - No. 16. - P. 183-189. DOI: 10.1016/j.matchemphys.2007.02.013

28. Dalgic S., Colakorullari M. Self-Diffusion Coefficients in Liquid Ag Using the Embedded Atom Model Based Effective Pair Potentials // Turkish Journal of Physics. - 2006. - Vol. 30. - No. 4. -P.303-310.

29. Daw M.S., Baskes M.I. Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals // Physical Review B. - 1984. - Vol. 29. - No. 12. - P. 6443-6453.

30. Daw M.S., Baskes M.I. Semiempirical, Quantum Mechanical Calculation of Hydrogen Embrittlement in Metals // Physical Review Letters. - 1983. - Vol. 50. - No. 17. - P. 1285-1288. DOI: 10.1103/PhysRevLett.50.1285

31. Application of the Born-Mayer potential with a Hard-Sphere scattering kernel to rarefied hyperthermal gas flow modeling / W.L. Dimpfl, I.J. Wysong, S.F. Gimelshein, M. Braunstein, L.S. Bernstein // AIP Conference Proceedings. - 2009. - Vol. 1084. - P. 323-328. DOI: 10.1063/1.3076494

32. Doyama M., Kagure Y. Embedded atom potential in fcc and bcc metals // Computational Materials Science. - 1999. - Vol. 14. - No. 1-4. - P. 80-83. DOI: 10.1016/S0927-0256(98)00076-7

33. Calculation of the structure factors of harmonic and anharmonic Fibonacci chains with molecular dynamics simulations / M. Engel, S. Sonntag, H. Lipp, H.-R. Trebin // Physical Review B. - 2007. -Vol. 75. - No. 15. - 144203 p. DOI: 10.1103/PhysRevB.75.144203

34. Ercolessi F., Adams J.B. Interatomic Potentials from First-Principles Calculations: The Force-Matching Method // Europhysics Letters. - 1994. - Vol. 26. - No. 8. - P. 583-588.

35. Fellinger M.R., Park H., Wilkins J.W. Force-matched embedded-atom method potential for niobium // Physical Review B. - 2010. - Vol. 81. - No. 14. DOI: 10.1103/PhysRevB.81.144119

36. Finnis M.W., Sinclair J.E. A simple empirical N-body potential for transition metals // Philosophical Magazine A. - 1984. - Vol. 50. - No. 1. - P. 44-55. DOI: 10.1080/01418618408244210

37. Foiles S.M. Application of the embedded-atom method to liquid transition metals // Physical Review B. - 1985. - V. 32, No. 6. - pp. 3409-3415. doi: 10.1103/PhysRevB.32.3409

38. Foiles S.M., Baskes M.I., Daw M.S. Embedded-atom-method functions for the fcc metals Cu, Ag, Au, Ni, Pd, Pt, and their alloys // Physical Review B. - 1986. - Vol. 33. - No. 12. - P. 7983-7991. DOI: 10.1103/PhysRevB.33.7983

39. Energetics of Surface Multilayer Relaxation on W(QQ1): Evidence for Short-Range Screening / Fu C.L., Ohnishi S., Wimmer E., Freeman A.J. // Physical Review B. - 1984. - Vol. 53. - No. 7. - P. 675678. DOI: 10.1103/PhysRevLett.53.675

40. Gelatt C.D., Ehrenreich H., Weiss J.A. Transition-metal hydrides: Electronic structure and the heats offormation // Physical Review B. - 1978. - Vol. 17. - No. 4. - P. 1940-1957. DOI: 10.1103/PhysRevB.17.1940

41. Dynamics of Radiation Damage / J.B. Gibson, A.N. Goland, M. Milgram, G.H. Vineyard // Phys. Rev. - 1960. - Vol. 120. - No. 4. - P. 1229-1253. DOI: 10.1103/PhysRev.120.1229

42. Gonzalez L.E., Gonzalez D.J., Canales M. Atomic dynamics in liquid lithium // Zeitschrift fur Physik B Condensed Matter. - 1996. - Vol. 100. - No. 4. - P. 601-611. DOI: doi.org/10.1007/s002570050167

43. Gordon P.A., Neeraj T., Mendelev M.I. Screw dislocation mobility in BCC metals: a refined potential description for Fe // Philosophical Magazine. - 2011. - Vol. 91. - No. 30. - P. 3931-3945. DOI: 10.1080/14786435.2011.597947

44. Numerical tools for atomistic simulations / P.M. Gullett, G.J. Wagner, A. Slepoy, M.F. Horstemeyer, H. Fang, M. Li, M.I. Baskes. - Livermore, SANDIA National Laboratories, 2003. - 68 p.

45. Johnson R.A. Alloy metals with the embedded-atom method // Physical Review B. - 1989. -Vol. 39. - No. 17. - P. 12554-12559. DOI: 10.1103/PhysRevB.39.12554

46. Johnson R.A. Analytic nearest-neighbor model for fcc metals // Physical Review B. - 1988. -Vol. 37. - No. 8. - P. 3924-3931. DOI: 10.1103/PhysRevB.37.3924

47. Jones J.E. On the Determination of Molecular Fields. II. From the Equation of State of a Gas // Proc. R. Soc. Lond. A. - 1924. - Vol. 106. - No. 738. - P. 463-477. DOI: 10.1098/rspa.1924.0082

48. Application of the Embedded Atom Method to Pb and Be / M. Karimi, Z. Yang, P. Tibbits, D. Ila, I. Dalins, G. Vidali // MRS Online Proceedings Library. - 1990. - Vol. 193. - No. 1. - P. 83-88. DOI: 10.1557/PROC-193-83

49. Kresse G., Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method // Physical Review B. - 1999. - Vol. 59. - No. 3. - P. 1758-1775. DOI: 10.1103/PhysRevB.59.1758

50. Law K., Stuart A., Zygalakis K. Data assimilation: A mathematical introduction. - Springer International Publishing, 2015. - 242 p.

51. Second nearest-neighbor modified embedded atom method potentials for bcc transition metals / B.-J. Lee, M.I. Baskes, H. Kim, Y.K. Cho // Physical Review B. - 2001. - Vol. 64. - No. 18. - 184102 p. DOI: 10.1103/PhysRevB.64.184102

52. Interatomic potentials for modelling radiation defects and dislocations in tungsten / M.-C. Marinica, L. Ventelon, M.R. Gilbert, L. Proville, S.L. Dudarev, J. Marian, G. Bencteux, F. Willaime // J. Phys.: Condens. Matter. - 2013. - Vol. 25. - No. 39. - 395502 p. DOI: 10.1088/09538984/25/39/395502

53. McClean A.D., McClean R.S. Roothaan-Hartree-Fock atomic wave functions Slater basis-set expansions for Z = 55-92 // Atomic Data and Nuclear Data Tables. - 1981. - Vol. 26. - No. 3-4. -P. 197-381. DOI: 10.1016/0092-640X(81)90012-7

54. Development of new interatomic potentials appropriate for crystalline and liquid iron / M.I. Mendelev, D.J. Srolovitz, G.J. Ackland, D.Y. Sun, M. Asta // Philosophical Magazine. - 2003. -Vol. 83. - No. 35. - P. 3977-3994. DOI: 10.1080/14786430310001613264

55. Mendelev M.I., Underwood T.L., Ackland G.J. Development of an interatomic potential for the simulation of defects, plasticity, and phase transformations in titanium // The Journal of Chemical Physics. - 2016. - Vol. 146. - No. 15. - 154102 p. DOI: 10.1063/1.4964654

56. The interaction of atomic hydrogen with Ni, Pd, and Pt clusters / R.P. Messmer, D.R. Salahub, K.H. Johnson, C.Y. Yang // Chemical Physics Letters. - 1977. - Vol. 51. - No. 1. - P. 84-89. DOI: 10.1016/0009-2614(77)85360-8

57. Farkas D., Mehl M.J., Papaconstantopoulos D.A. Interatomic potentials for monoatomic metals from experimental data and ab initio calculations // Physical Review B. - 1999. - Vol. 59. - No. 5. -P. 3393-3407. DOI: 10.1103/PhysRevB.59.3393

58. Structural stability and lattice defects in copper: Ab initio, tight-binding and embedded atom method calculations / Y. Mishin, M.J. Mehl, D.A. Papaconstantopoulos, A.F. Voter, J.D. Kress // Physical Review B. - 2001. - Vol. 63. - No. 22. - 224106 p. DOI: 10.1103/PhysRevB.63.224106

59. Moitra A., Kim S.-G., Horstemeyer M.F. Structural and thermal properties of calcium using an MEAM potential // Calphad. - 2001. - Vol. 35. - No. 2. - P. 262-268. DOI: 10.1016/j.calphad.2011.01.002

60. Morse P.M. Diatomic Molecules According to the Wave Mechanics. II. Vibrational Levels // Physical Review. - 1929. - Vol. 34. - No. 1 - P. 57-64. DOI: https://doi.org/10.1103/PhysRev.34.57

61. Niklasson A.M.N., Cawkwell M.J. Fast method for quantum mechanical molecular dynamics // Physical Review B. - 2012. - Vol. 86. - No. 17. - 174308 p. DOI: 10.1103/PhysRevB.86.174308

62. Norskov J.K., Lang N.D. Effective-medium theory of chemical binding: Application to chemisorp-tion // Physical Review Letters. - 1980. - Vol. 21. - No. 6. - P. 2131-2136. DOI: 10.1103/PhysRevB.21.2131

63. Rapaport D.C. The Art of Molecular Dynamics Simulation. - Cambridge University Press, 2004. - 549 p.

64. Universal features of the equation of state of metals / J.H. Rose, J.R. Smith, F. Guinea, J. Ferrante // Physical Review B. - 1984. - Vol. 29. - No. 6. - P. 2963-2969. DOI: 10.1103/PhysRevB.29.2963

65. Schommers W. Pair potentials in disordered many-particle systems. A study for liquid gallium // Physical Review A. - 1983. - Vol. 28. - No. 6. - P. 3599-3605. DOI: 10.1103/PhysRevA.28.3599

66. Slater J.C. Quantum Theory of Matter, 2nd edition. - New York, McGraw-Hill, 1968. - 321 p.

67. Smirnova D.E., Starikov S.V., Stegailov V.V. Interatomic potential for uranium in a wide range of pressures and temperatures // J. Phys.: Condens. Matter. - 2012. - Vol. 24. - No. 14. - 015702. DOI: 10.1088/0953-8984/24/14/149501

68. Stockmayer W.H. Second Virial Coefficients of Polar Gases // J. Chem. Phys. - 1941. -Vol. 9. - No. 5. - P. 398-402. DOI: 10.1063/1.1750922

69. Stott M.J., Zaremba E. Quasiatoms: An approach to atoms in nonuniform electronic systems // Physical Review B. - 1980. - Vol. 22. - No. 4. - P. 1564-1583. DOI: 10.1103/PhysRevB.22.1564

70. Thomas J.F. Third-Order Elastic Constants of Aluminum // Physical Review. - 1968. -Vol. 175. - No. 3. - P. 955-962. DOI: 10.1103/PhysRev.175.955

71. Vicente J., Lanchares J., Hermida R. Placement by thermodynamic simulated annealing // Physics Letters A - 2003. - Vol. 317. - No. 5-6. - P. 415-423. DOI: 10.1016/j.physleta.2003.08.070

72. Warshel A, Levitt M.J. Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme // Journal of Molecular Biology. -1976. - Vol. 103. - No. 2. - P. 227-249. DOI: 10.1016/0022-2836(76)90311-9

73. Winey J.M., Kubota A., Gupta Y.M. A thermodynamic approach to determine accurate potentials for molecular dynamics simulations: thermoelastic response of aluminum // Modelling Simul. Mater. Sci. Eng. - 2010. - Vol. 18. - No. 5. - 029801. DOI: 10.1088/0965-0393/17/5/055004

References

1. Gill P.E., Murray W., Wright M.H. Practical Optimization. Michigan, Academic Press, 1981, 401 p.

2. Krivtsov A.M., Teoreticheskaya mekhanika. Uprugiye svoystva odnoatomnykh i dvukhatomnykh kristallov [Theoretical mechanics. Elastic properties of monatomic and diatomic crystals]. Saint Petersburg, Polytechnic University Publishing House, 2009, 126 p.

3. Krivtsov A.M., Krivtsova N.V. Metod chastits i yego ispol'zovaniye v mekhanike deformiruyemogo tverdogo tela [Particle method and its use in the mechanics of a deformable solid]. Dal'nevostochnyy matematicheskiy zhurnal [Far-Eastern Mathematical Journal]. 2002, vol. 3, no. 2, pp. 254-276.

4. Maize J. Teoriya i zadachi mekhaniki sploshnykh sred [Theory and Problems of Continuum Mechanics]. Moscow, Mir [World], 1974, 319 p.

5. Hartree D. Raschety atomnykh struktur [Calculations of atomic structures]. Moscow, Publishing House of Foreign Literature, 1960, 272 p.

6. Agrawal A., Mishra R., Ward L., Flores K. M, Windl W. An embedded atom method potential of beryllium. Modelling and Simulation in Materials Science and Engineering, 2013, 21 (8). DOI: 10.1088/09650393/21/8/085001

7. Alder B.J., Waingwright T.E. Phase Transition for a Hard Sphere System. J. Chem. Phys., 1957, 27 (5), pp.1208-1209.

8. Allen M. P., Tildesley A. K. Computer Simulation of Liquids. Oxford, Clarendon Press, 1987, 385 p.

9. Bangwei Z., Yifang O., Shuzhi L., Zhanpeng J. An analytic MEAM model for all BCC transition metals. Physica B: Condensed Matter, 1999, 262 (3-4), pp. 218-225. DOI: 10.1016/S0921-4526(98)01156-9

10. Baskes M.I. Application of the Embedded-Atom Method to Covalent Materials: A Semiempirical Potential for Silicon. Physical Review Letters, 1987, 59 (23), pp. 2666-2669. DOI: 10.1103/PhysRevLett.59.2666

11. Baskes M.I. Determination of modified embedded atom method parameters for nickel. Materials Chemistry and Physics, 1997, 50 (2), pp. 152-158. DOI:10.1016/S0254-0584(97)80252-0

12. Baskes M.I. Modified embedded-atom potentials for cubic materials and impurities. Physical Review B, 1992, 46 (5), pp. 2727-2742. DOI: 10.1103/PhysRevB.46.2727

13. Baskes M.I., Johnson R.A. Modified embedded atom potentials for HCP metals. Modelling and Simulation in Materials Science and Engineering, 1994, 2 (1), pp. 147-163. DOI: 10.1088/0965-0393/2/1/011

14. Baskes M.I., Sha X., Angelo J.E., Moody N.R. Trapping of hydrogen to lattice defects in nickel. Modelling and Simulation in Materials Science and Engineering, 1995, 5 (6), pp. 651-652. DOI: 10.1088/09650393/3/3/001

15. Belashchenko D. K. Embedded Atom Method Potentials for Liquid Copper and Silver. Inorganic Materials, 2012, 48 (9), pp. 940-947. DOI:10.1134/S002016851209004X

16. Belashchenko D. K., Smirnova D. E., Ostrovski O. I. Molecular-Dynamic Simulation of the Thermophysical Properties of Liquid Uranium. High Temperature, 2010, 48 (3), pp. 363-375. DOI: 10.1134/S0018151X10030107

17. Belashchenko D. K., Zhuravlev Yu. V. Application of the embedded-atom method to liquid copper. Inorganic Materials, 2008, 44 (9), pp. 939-945. DOI:10.1134/S0020168508090082

18. Belashchenko D.K., Ostrovskii O.I. Application of the embedded atom model to liquid metals: Liquid lithium. High Temperature, 2008, 47 (2), pp. 231-237. DOI: 10.1134/S0018151X09020102

19. Brommer P., Gahler F. Effective potentials for quasicrystals from ab-initio data. Philosophical Magazine, 2006, 86 (6-8), pp. 753-758. DOI:10.1080/14786430500333349

20. Cai J., Ye. Y.Y Simple analytical embedded-atom- potential model including a long-range force for fcc metals and their alloys. Physical Review B, 1996, 54 (12), pp. 8398-8410. DOI:10.1103/PhysRevB.54.8398

21. Chamati H., Papanicolaou N.I., Mishin Y., Papaconstantopoulos D.A. Embedded-atom potential for Fe and its application to self-diffusion on Fe(100). Surface Science, 2006, 600 (9), pp. 1793-1803. DOI: 10.1016/j.susc.2006.02.010

22. Chandler D. Introduction to Modern Statistical Mechanics. New York, Oxford University Press, 1987, 274 p.

23. Chantasiriwan S., Milstein F. Embedded-atom models of 12 cubic metals incorporating second- and third-order elastic-moduli data. Physical Review B, 1998, 58 (10), pp. 5996-6005. DOI: 10.1103/PhysRevB.58.5996

24. Cherne F.J., Baskes M.I., Deymier P.A. Properties of liquid nickel: A critical comparison of EAM and MEAM calculations. Physical Review B, 2002, 65 (2), 024209. DOI:10.1103/PhysRevB.65.024209

25. Clementi E., Roetti C. Roothaan-Hartree-Fock atomic wavefunctions: Basis functions and their coefficients for ground and certain excited states of neutral and ionized atoms, Z<54. Atomic Data and Nuclear Data Tables, 1974,13 (3-4), pp. 177-478. DOI: 10.1016/S0092-640X(74)80016-1

26. Copper A.S. Precise Lattice Constants of Germanium, Aluminum, Gallium Arsenide, Uranium, Sulphur, Quartz and Sapphire. Acta Crystallographica, 1962, 15 (6), pp. 578-582. DOI: 10.1107/S0365110X62001474

27. Dalgic S. The EAM based effective pair potentials for the dynamic properties of liquid Cu. Materials Chemistry and Physics, 2007, 39 (16), pp. 183-189. DOI: 10.1016/j.matchemphys.2007.02.013

28. Dalgic S., Colakorullari M. Self-Diffusion Coefficients in Liquid Ag Using the Embedded Atom Model Based Effective Pair Potentials. Turkish Journal of Physics, 2006, 30 (4), pp. 303-310.

29. Daw M.S., Baskes M.I. Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals. Physical Review B, 1984, 29 (12), pp. 6443-6453.

30. Daw M.S., Baskes M.I. Semiempirical, Quantum Mechanical Calculation of Hydrogen Embrittlement in Metals. Physical Review Letters, 1983, 50 (17), pp. 1285-1288. DOI: 10.1103/PhysRevLett.50.1285

31. Dimpfl W.L., Wysong I.J., Gimelshein S.F., Braunstein M., Bernstein L.S. Application of the BornMayer potential with a Hard-Sphere scattering kernel to rarefied hyperthermal gas flow modeling. AIP Conference Proceedings, 2009, 1084, pp. 323-328. DOI: 10.1063/1.3076494

32. Doyama M., Kagure Y. Embedded atom potential in fcc and bcc metals. Computational Materials Science, 1999, 14 (1-4), pp. 80-83. DOI: 10.1016/S0927-0256(98)00076-7

33. Engel M., Sonntag S., Lipp H., Trebin H.-R. Calculation of the structure factors of harmonic and anhar-monic Fibonacci chains with molecular dynamics simulations. Physical Review B, 2007, 75 (15), 144203. DOI: 10.1103/PhysRevB.75.144203

34. Ercolessi F., Adams J.B. Interatomic Potentials from First-Principles Calculations: The Force-Matching Method. Europhysics Letters, 1994, 26 (8), pp. 583- 588.

35. Fellinger M.R., Park H., Wilkins J.W. Force-matched embedded-atom method potential for niobium. Physical Review B, 2010, 81 (14), 134303. DOI: 10.1103/PhysRevB.81.144119

36. Finnis M.W., Sinclair J. E. A simple empirical N-body potential for transition metals. Philosophical Magazine A, 1984, 50 (1), pp. 44-55. DOI:10.1080/01418618408244210

37. Foiles S.M. Application of the embedded-atom method to liquid transition metals. Physical Review B, 1985, 32 (6), pp. 3409-3415. DOI: 10.1103/PhysRevB.32.3409

38. Foiles S.M., Baskes M.I., Daw M.S Embedded-atom-method functions for the fcc metals Cu, Ag, Au, Ni, Pd, Pt, and their alloys. Physical Review B, 1986, 33 (12), pp. 7983-7991. DOI: 10.1103/PhysRevB.33.7983

39. Fu C.L., Ohnishi S., Wimmer E., Freeman A.J. Energetics of Surface Multilayer Relaxation on W(QQ1): Evidence for Short-Range Screening. Physical Review B, 1984, 53 (7), pp. 675-678. DOI: 10.1103/PhysRevLett.53.675

40. Gelatt C.D., Ehrenreich H., Weiss J.A. Transition-metal hydrides: Electronic structure and the heats of formation. Physical Review B, 1978, 17 (4), pp.1940-1957. DOI: 10.1103/PhysRevB.17.1940

41. Gibson J.B., Goland A. N., Milgram M., Vineyard G.H. Dynamics of Radiation Damage. Physical Review, 1960, 120 (4), pp.1229-1253. DOI: 10.1103/PhysRev.120.1229

42. Gonzalez L.E., Gonzalez D.J., Canales M. Atomic dynamics in liquid lithium. Zeitschrift für Physik B Condensed Matter, 1996, 100 (4), pp. 601-611. DOI: doi.org/10.1007/s002570050167

43. Gordon P.A., Neeraj T., Mendelev M.I. Screw dislocation mobility in BCC metals: a refined potential description for Fe. Philosophical Magazine, 2011, 91 (30), pp. 3931-3945. DOI: 10.1080/14786435.2011.597947

44. Gullett P.M., Wagner G.J., Slepoy A., Horstemeyer M.F., Fang H., Li M., Baskes M.I. Numerical tools for atomistic simulations. Livermore, SANDIA National Laboratories, 2003, 68 p.

45. Johnson R.A. Alloy metals with the embedded-atom method. Physical Review B, 1989, 39 (17), pp. 12554-12559. DOI: 10.1103/PhysRevB.39.12554

46. Johnson R.A. Analytic nearest-neighbor model for fcc metals. Physical Review B, 1988, 37 (8), pp. 3924-3931. DOI: 10.1103/PhysRevB.37.3924

47. Jones J. E., On the Determination of Molecular Fields. II. From the Equation of State of a Gas. Proc. R. Soc. Lond. A, 1924, 106 (738), pp. 463-477. DOI: 10.1098/rspa.1924.0082

48. Karimi M., Yang Z., Tibbits P., Ila D., Dalins I., Vidali G. Application of the Embedded Atom Method to Pb and Be, MRS Online Proceedings Library, 1990, 193 (1), pp. 83-88. DOI: 10.1557/PROC-193-83

49. Kresse G., Joubert D. From ultrasoft pseudopotentials to the projector augmented-wave method. Physical Review B, 1999, 59 (3), pp. 1758-1775. DOI:10.1103/PhysRevB.59.1758

50. Law K., Stuart A., Zygalakis K. Data assimilation: A mathematical introduction. Springer International Publishing, 2015, 242 p.

51. Lee B.-J., Baskes M.I., Kim H., Cho Y.K. Second nearest-neighbor modified embedded atom method potentials for bcc transition metals. Physical Review B, 2001, 64 (18), 184102. DOI:10.1103/PhysRevB.64.184102

52. Marinica M.-C., Ventelon L., Gilbert M.R., Proville L., Dudarev S.L., Marian J., Bencteux G., Willaime F. Interatomic potentials for modelling radiation defects and dislocations in tungsten. J. Phys.: Condens. Matter, 2013, 25 (39), 395502. DOI: 10.1088/0953-8984/25/39/395502

53. McClean A.D., McClean R.S. Roothaan-Hartree-Fock atomic wave functions Slater basis-set expansions for Z = 55-92. Atomic Data and Nuclear Data Tables, 1981, 26 (3-4), pp. 197-381. DOI: 10.1016/0092-640X(81)90012-7

54. Mendelev M.I., Srolovitz D.J., Ackland G.J., Sun D.Y., Asta M. Development of new interatomic potentials appropriate for crystalline and liquid iron. Philosophical Magazine, 2003, 83 (35), pp. 3977-3994. DOI:10.1080/14786430310001613264

55. Mendelev M.I., Underwood T.L., Ackland G.J. Development of an interatomic potential for the simulation of defects, plasticity, and phase transformations in titanium. The Journal of Chemical Physics, 2016, 146 (15), 154102. DOI: 10.1063/1.4964654

56. Messmer R.P., Salahub D.R., Johnson K.H., Yang C.Y. The interaction of atomic hydrogen with Ni, Pd, and Pt clusters. Chemical Physics Letters, 1977, 51 (1), pp. 84-89. DOI: 10.1016/0009-2614(77)85360-8

57. Mishin Y., Farkas D., Mehl M.J., Papaconstantopoulos D.A. Interatomic potentials for monoatomic metals from experimental data and ab initio calculations. Physical Review B, 1999, 59 (5), pp. 3393-3407. DOI: 10.1103/PhysRevB.59.3393

58. Mishin Y., Mehl M.J., Papaconstantopoulos D.A., Voter A.F., Kress J.D. Structural stability and lattice defects in copper: Ab initio, tight-binding and embedded atom method calculations. Physical Review B, 2001, 63 (22), 224106. DOI: 10.1103/PhysRevB.63.224106

59. Moitra A., Kim S.-G., Horstemeyer M.F. Structural and thermal properties of calcium using an MEAM potential. Calphad, 2001, 35 (2), pp. 262-268. DOI: 10.1016/j.calphad.2011.01.002

60. Morse P.M. Diatomic Molecules According to the Wave Mechanics. II. Vibrational Levels. Physical Review, 1929, 34 (1), pp. 57-64. DOI:10.1103/PhysRev.34.57

61. Niklasson A.M.N., Cawkwell M.J. Fast method for quantum mechanical molecular dynamics. Physical Review B, 2012, 86 (17), 174308. DOI: 10.1103/PhysRevB.86.174308

62. Norskov J.K., Lang N.D. Effective-medium theory of chemical binding: Application to chemisorption. Physical Review Letters, 1980, 21 (6), pp. 2131-2136. DOI: 10.1103/PhysRevB.21.2131

63. Rapaport D.C. The Art of Molecular Dynamics Simulation. Cambridge, Cambridge University Press, 2004, 549 p.

64. Rose J.H., Smith J.R., Guinea F., Ferrante J. Universal features of the equation of state of metals. Physical Review B, 1984, 29 (6), pp. 2963-2969. DOI: 10.1103/PhysRevB.29.2963

65. Schommers W. Pair potentials in disordered many-particle systems. A study for liquid gallium. Physical Review A, 1984, 28 (6), pp. 3599-3605. DOI: 10.1103/PhysRevA.28.3599

66. Slater J.C. Quantum Theory of Matter, 2nd edition. New York, McGraw-Hill, 1968, 321 p.

67. Smirnova D.E., Starikov S.V., Stegailov V.V. Interatomic potential for uranium in a wide range of pressures and temperatures. J. Phys.: Condens. Matter, 2012, 24 (14), 015702. DOI: 10.1088/0953-8984/24/14/149501

68. Stockmayer W. H. Second Virial Coefficients of Polar Gases. J. Chem. Phys., 1941, 9 (5), pp. 398-402. DOI: 10.1063/1.1750922

69. Stott M.J., Zaremba E. Quasiatoms: An approach to atoms in nonuniform electronic systems. Physical Review B, 1980, 22 (4), pp. 1564-1583. DOI: 10.1103/PhysRevB.22.1564

70. Thomas J.F. Third-Order Elastic Constants of Aluminum. Physical Review, 1968, 175 (3), pp. 955-962. DOI: 10.1103/PhysRev.175.955

71. Vicente J., Lanchares J., Hermida R. Placement by thermodynamic simulated annealing. Physics Letters A, 2003, 317 (5-6), pp. 415-423. DOI:10.1016/j.physleta.2003.08.070

72. Warshel A, Levitt M. J. Theoretical studies of enzymic reactions: Dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme. Journal of Molecular Biology, 1976, 103 (2), pp. 227-249. DOI: 10.1016/0022-2836(76)90311-9

73. Winey J.M., Kubota A., Gupta Y.M. A thermodynamic approach to determine accurate potentials for molecular dynamics simulations: thermoelastic response of aluminum. Modelling Simul. Mater. Sci. Eng., 2010, 18 (5), 029801. DOI: 10.1088/0965-0393/17/5/055004