Формостабильность металлических нанопластин: инсайты из молекулярно-динамических симуляций

Аннотация

Металлические нанопластины демонстрируют замечательную функциональную универсальность, неразрывно связанную с их морфологией. С помощью атомистической молекулярной динамики мы исследовали формостабильность нанопластин кубически-объёмно-центрированных (bcc) металлов с низкоиндексными поверхностями. Симуляции показали, что (110) поверхности интринзически наиболее стабильны, тогда как (111) и (001) нанопластины претерпевают температурно-обусловленные деформации. (001) пластина трансформируется в седлообразную структуру, состоящую из (110) граней, перед тем как коллапсировать в неправильную частицу. Эволюция управляется распределением сдвигающих напряжений, зависящих от ориентации граней. Модулируя ориентацию граней, размер пластины (диаметр и толщину) и материал, мы продемонстрировали настраиваемую формостабильность для нескольких bcc металлов (Fe, W, Nb, Mo, Cr). Результаты дают атомистическую основу для контроля морфологии нанопластин при синтезе и применении.

Введение

Металлические нанопластины — ключевые функциональные наноматериалы с применением от катализа до оптоэлектроники. Их производительность диктуется площадью поверхности, химией граней и общей формой. Хотя термодинамика предсказывает стремление к сферической минимизации поверхностной энергии, кинетические барьеры часто стабилизируют анизотропные пластины. Тонкие пластины (несколько атомных слоев) обладают высоким отношением площадь/объём, делая их особенно чувствительными к температуре и среде. Предыдущие in-situ TEM исследования наблюдали фрагментацию краев и реориентацию граней при нагреве, однако количественный атомистический понимание оставалось ограниченным. Здесь мы используем крупномасштабную молекулярную динамику (MD) для захвата полной температурно-зависимой морфологической эволюции bcc металлических нанопластин, фокусируясь на роли кристаллической ориентации, размера и состава.

Методы

Fe нанопластины построили с параметром решётки a = 2,8665 Å и диаметром d = 32a, каждая из трёх атомных слоев. Начальные поверхности задали ориентациями (111), (001) или (110). Аналогичные структуры сгенерировали для W, Nb, Mo и Cr. Модели собрали в LAMMPS и релаксировали сопряжёнными градиентами. Последующий нагрев проводили в NVT ансамбле, повышая температуру от 1 К до 300 К (или выше) с шагом 1 К, шаг по времени 2 фс и 200 пс уравновешивания на шаге. Потенциальная энергия и тензоры напряжений записывали; последние 8 пс каждого шага усредняли. Межатомные силы описывали EAM потенциалами, валидированными в предыдущих MD исследованиях. Локальное напряжение для атома i вычисляли стандартным вириальным выражением:

σ_{αβ}=1/(2Ω_i) Σ_{j≠i} F_{ij}^α R_{ij}^β

Результаты и обсуждение

Энергетическая иерархия и эволюция формы. Начальные потенциальные энергии на атом составляли −2,833, −3,457 и −3,668 эВ/атом для (111), (001) и (110) Fe пластин соответственно, зеркаля поверхностные энергии 2,58, 2,47 и 2,37 Дж/м². При повышении температуры (111) и (001) пластины быстро принимают искривлённые метастабильные формы, в то время как (110) пластина сохраняет плоскую геометрию до самых высоких смоделированных температур. (001) пластина развивает седлообразную форму, доминируемую (110) гранями, что подтверждается аппроксимацией среднего слоя квадрической поверхностью z = ax² + by² + cxy + d. Критические переходы происходят при 8, 129 и 205 К, с соответствующими изменениями коэффициентов аппроксимации и R² > 0,8, указывающими на надёжную седлообразную геометрию до 205 К.

Напряжённо-обусловленная деформация. Карты сдвигающих напряжений раскрывают анизотропное распределение, связанное с ориентацией граней: положительные напряжения на красных гранях и отрицательные на синих создают противоположные силы, вызывающие изгиб. При 9 К увеличенный изгиб устраняет напряжённые разрывы, снижая потенциальную энергию. Вторичный переход между 129–134 К соответствует выбухиванию центральной области, вызывая внезапное падение R² и последующий коллапс в неправильную частицу.

Модуляция граней и материальная зависимость. «Модулированная (110)» Fe пластина, состоящая из четырёх по-разному ориентированных граней, показывает похожее седлообразование при 1 К, но возвращается к диску при 179 К и окончательно коллапсирует при 277 К. В противовес, регулярная (110) пластина остаётся стабильной до 552 К, подчёркивая влияние дизайна граней. Расширение исследования на W, Nb, Mo и Cr пластины d = 32a демонстрирует, что стабильность седла варьируется от 62 К (Cr) до 582 К (W), коррелируя с начальной потенциальной энергией и структурной стабильностью каждого элемента.

Размерные и толщинные эффекты. Fe (001) пластины с диаметрами 12a, 40a и 50a демонстрируют различные окна стабильности: 12a пластина коллапсирует к 84 К, 40a удерживает седло до 190 К, а 50a — только до 134 К. Вариации толщины показывают, что 4-слойные пластины сохраняют форму до 98 К перед конвертацией в (110) грани; более тонкие пластины коллапсируют раньше. В целом, более крупные пластины улучшают базовую стабильность, но кинетические и энтропийные вклады могут снижать седлостабильность при очень больших диаметрах.

Заключение

Наше MD исследование раскрывает атомистические механизмы, управляющие формоэволюцией bcc металлических нанопластин. (110) поверхности обеспечивают наивысшую структурную стабильность, тогда как (111) и (001) пластины трансформируются в седлообразные структуры из (110) граней перед коллапсом. Сдвигающие напряжения, связанные с ориентацией граней, движитель этих деформаций. Инженеря ориентации граней, размеров пластин и выбирая состав материала, стабильность желаемых морфологий может быть тонко настроена. Эти инсайты закладывают теоретическую основу для контроля морфологии при синтезе металлических наноматериалов.

Доступность данных и материалов

Все симуляционные данные и скрипты анализа включены в данную публикацию.

Аббревиатуры

a

Параметр решётки

bcc

Кубически-объёмно-центрированная

d

Диаметр

EAM

Встроенный атомный метод

E_p

Потенциальная энергия

MD

Молекулярная динамика

R²

Коэффициент детерминации

TEM

Трансмиссионная электронная микроскопия