Метод погружённого атома

Метод погружённого атома (англ. Embedded Atom Method, EAM) — это полуэмпирический метод атомистического моделирования, используемый для расчёта полной энергии и межатомных взаимодействий в металлических системах. Метод основан на представлении, что энергия атома в кристаллической решётке складывается из энергии парных взаимодействий с соседями и энергии «погружения» атома в электронное облако, создаваемое всеми остальными атомами. EAM широко применяется в материаловедении, физике твёрдого тела и вычислительной химии для изучения механических свойств, дефектов, фазовых переходов и поверхностных явлений в металлах и сплавах.

История

Метод погружённого атома был разработан в 1983 году американскими физиками Мюрреем Доусом и Майком Баскесом в Ливерморской национальной лаборатории имени Лоуренса (США). Исходной мотивацией стало преодоление ограничений классических парных потенциалов (например, потенциала Леннарда-Джонса), которые неадекватно описывали упругие свойства и энергию дефектов в металлах. Доус и Баскес предложили добавить в функционал энергии нелинейный член, зависящий от локальной электронной плотности, что позволило более точно воспроизводить экспериментальные данные для переходных металлов.

В 1984 году метод был обобщён на сплавы и получил название «метод погружённого атома для сплавов» (EAM-alloy). В последующие десятилетия EAM был расширен на системы с неметаллическими элементами (например, водород в металлах), а также дополнен учётом угловых зависимостей (модификации MEAM, Modified Embedded Atom Method). В России метод активно применяется в Институте физики прочности и материаловедения СО РАН и других научных центрах.

Основные принципы

Физическая модель

В основе EAM лежит предположение, что полная энергия системы атомов \(E_{tot}\) может быть записана как сумма двух вкладов:

\[ E_{tot} = \sum_i F_i(\rho_i) + \frac{1}{2} \sum_{i \neq j} \phi_{ij}(r_{ij}) \]

где:

• \(F_i(\rho_i)\) — функция погружения, зависящая от локальной электронной плотности \(\rho_i\) в точке расположения атома \(i\);
• \(\phi_{ij}(r_{ij})\) — парный потенциал, описывающий отталкивание и притяжение между атомами \(i\) и \(j\) на расстоянии \(r_{ij}\).

Локальная электронная плотность \(\rho_i\) вычисляется как суперпозиция вкладов от всех соседних атомов:

\[ \rho_i = \sum_{j \neq i} \rho_j^{atom}(r_{ij}) \]

где \(\rho_j^{atom}(r)\) — электронная плотность, создаваемая атомом \(j\) на расстоянии \(r\). Обычно эта функция аппроксимируется экспоненциальной или степенной зависимостью.

Функция погружения

Функция погружения \(F(\rho)\) является ключевым элементом EAM. Она отражает энергетический выигрыш или проигрыш при помещении атома в электронное облако заданной плотности. Для большинства металлов \(F(\rho)\) имеет минимум при некотором равновесном значении плотности, соответствующем равновесному межатомному расстоянию. Вблизи минимума функция часто аппроксимируется квадратичной зависимостью, что приводит к линейной зависимости силы от смещения (упругое поведение).

Парный потенциал

Парный потенциал \(\phi_{ij}(r)\) обычно берётся в виде короткодействующего отталкивательного члена (например, потенциал Борна-Майера) и дальнодействующего притяжения (например, потенциал Морзе). В EAM парный потенциал не является единственным вкладом в энергию, что отличает метод от классических парных потенциалов.

Параметризация

Для практического применения EAM необходимо задать три функции: \(\rho^{atom}(r)\), \(F(\rho)\) и \(\phi(r)\). Эти функции не берутся из первых принципов, а подгоняются под экспериментальные данные или результаты квантово-механических расчётов. Типичный набор данных для подгонки включает:

• равновесную постоянную решётки;
• энергию когезии;
• упругие постоянные (\(C_{11}, C_{12}, C_{44}\));
• энергию образования вакансии;
• поверхностную энергию.

Для сплавов дополнительно подгоняются параметры смешения, такие как энергия растворения и энергия образования упорядоченных фаз. Существуют стандартные библиотеки потенциалов EAM для многих металлов (например, для алюминия, меди, никеля, железа, вольфрама), разработанные разными исследовательскими группами.

Применение

Механические свойства

EAM широко используется для моделирования пластической деформации, распространения дислокаций, зарождения трещин и деформационного упрочнения. Например, с помощью EAM изучают поведение границ зёрен в поликристаллах, что важно для понимания прочности и пластичности металлов.

Дефекты кристаллической решётки

Метод позволяет рассчитывать энергию и структуру точечных дефектов (вакансий, межузельных атомов), линейных дефектов (дислокаций) и объёмных дефектов (пор, включений). Результаты EAM хорошо согласуются с экспериментальными данными для многих металлов.

Фазовые переходы

EAM применяется для моделирования плавления, кристаллизации, мартенситных превращений и полиморфных переходов в металлах. Например, для железа метод воспроизводит переход из объёмно-центрированной кубической (ОЦК) в гранецентрированную кубическую (ГЦК) структуру при высоких температурах.

Поверхности и тонкие плёнки

Метод позволяет изучать реконструкцию поверхностей, адсорбцию атомов, рост тонких плёнок и эпитаксию. EAM использовался для моделирования роста наночастиц и нанопроволок.

Радиационные повреждения

В атомной энергетике EAM применяется для моделирования каскадов смещений, вызванных облучением, и образования дефектов в конструкционных материалах (например, в циркониевых сплавах и аустенитных сталях).

Ограничения и критика

Несмотря на широкое распространение, метод погружённого атома имеет ряд ограничений:

• Эмпирический характер: параметры потенциалов подгоняются под ограниченный набор данных, что может приводить к ошибкам при экстраполяции на нетипичные условия.
• Отсутствие угловых зависимостей: классический EAM не учитывает направленность химических связей, что делает его малопригодным для полупроводников и керамик.
• Неучёт электронных корреляций: метод не описывает эффекты, связанные с перераспределением электронной плотности при деформации (например, изменение магнитного момента в ферромагнетиках).
• Вычислительная сложность: хотя EAM быстрее методов квантовой механики (DFT), он всё же требует значительных ресурсов при моделировании систем с миллионами атомов.

Для преодоления этих ограничений были разработаны модификации, такие как MEAM (Modified Embedded Atom Method), который включает угловые зависимости, и COMB (Charge-Optimized Many-Body) потенциалы, учитывающие переменный заряд атомов.

Интересные факты

• Метод погружённого атома стал основой для создания многих коммерческих и открытых программных пакетов молекулярной динамики, включая LAMMPS, NAMD и DL_POLY.
• В 2010 году группа российских учёных из Института физики прочности и материаловедения СО РАН (Томск) разработала потенциал EAM для сплава титан-алюминий, который используется для моделирования жаропрочных материалов.
• EAM позволил впервые предсказать существование «сверхпластичности» в нанокристаллических металлах, что было подтверждено экспериментально.

Источники

• Daw, M. S., & Baskes, M. I. (1984). Embedded-atom method: Derivation and application to impurities, surfaces, and other defects in metals. Physical Review B, 29(12), 6443–6453.
• Foiles, S. M., Baskes, M. I., & Daw, M. S. (1986). Embedded-atom-method functions for the fcc metals Cu, Ag, Au, Ni, Pd, Pt, and their alloys. Physical Review B, 33(12), 7983–7991.
• Baskes, M. I. (1992). Modified embedded-atom potentials for cubic materials and impurities. Physical Review B, 46(5), 2727–2742.
• Lee, B. J., & Shim, J. H. (2004). A modified embedded-atom method interatomic potential for the Fe–Al system. Calphad, 28(2), 125–132.
• Харченко, В. В., & Старцев, С. В. (2015). Моделирование методом погружённого атома дефектов в металлах. Физика твёрдого тела, 57(10), 1951–1958.