Исследование нанообработки на кончике зонда методом молекулярной динамики: обзор

Аннотация

Методы нанообработки на кончике зонда (Tip-Based Nanomachining, TBN) доказали свою эффективность для изготовления микроструктур. Моделирование методом молекулярной динамики (МД) широко применяется для изучения механизмов TBN, которые не удается полностью раскрыть экспериментально. В данном обзоре представлены недавние научные достижения в МД-моделировании TBN. Рассмотрены методы построения моделей для различных материалов, анализ механизма обработки (силы резания, удаление материала, дефекты подповерхностного слоя), а также текущие недостатки и перспективы развития метода.

Введение

Микро/нанотехнологии находят применение в экологике, энергетике, биологии, медицине, обороне и других сферах, играя ключевую роль в технологическом прогрессе. Для высокоточного нанофорирования метод TBN на базе зонда атомно-силовой микроскопии (АСМ) выделяется простотой оборудования, нанометровой точностью и низкими требованиями к среде. На сегодняшний день с его помощью успешно изготовлены наноточки, линии/канавки, 2D/3D-структуры и наноструктуры на изогнутых поверхностях — основные элементы наноприборов.

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

Ранние работы (Fang et al., Isono & Tanaka, Yan et al.) изучали влияние скорости сcriberов, угла конуса, температуры и геометрии кончика. Недавно появились новые подходы: Xiao et al. сравнили статическую и динамическую плуговую литографию, Geng et al. моделировали нагрузочно-контролируемое наноскребливание, Ren et al. предложили критерий «грубое-тонкое» для определения минимального шага подачи. Предыдущие обзоры не охватывали эти достижения, поэтому данная работа фокусируется на современном состоянии МД-моделирования TBN.

Методы моделирования

Построение МД-модели

Точность симуляции определяется атомной конфигурацией и потенциалами взаимодействия. Кристаллические материалы (Cu — ФЦК, Fe — ОЦК, Ti — УПГ, Si — алмазная структура) моделируют через периодическое дублирование элементарной ячейки (однокристаллики) или метод Вороного (поликристаллики). Наноблizkированные (NT) структуры, обладающие уникальной прочностью и пластичностью, строятся путем формирования многослойных пакетов с границами-близнецами и последующей генерации зерен по Вороного. Аморфные полимеры (например, полиэтилен) создают алгоритмом Монте-Карло самово избегающей случайной ходьбы.

Геометрия кончика (конус, треугольная пирамида, полусфера) критически влияет на процесс. Обычно кончик считают абсолютно твердым. Выбор потенциала — ключевой момент: EAM для металлов, Терсофф и Stillinger-Weber (SW) для ковалентных структур (Si), ABOP/REBO/AIREBO для SiC и углеродных наноструктур. Основные инструменты — LAMMPS, Atomsk (сложные структуры), Material Studio (полимеры), OVITO и VMD (визуализация). Ограничения: типичные размеры модели < 50 нм, отсутствие потенциалов для GaAs, Lu₂O₃ и др.

Конфигурация обработки

В экспериментах преобладает нагрузочно-контролируемый режим (постоянная нормальная нагрузка), обеспечивающий высокую точность на наклонных поверхностях. Многие МД-исследования используют смещенческий контроль, что вносит расхождения. Реализован МД-подход с прямым заданием нормальной нагрузки. Процесс делится на стадии: релаксация, проникновение, skreбливание. Учитывается наклон кончика из-за изгиба балки АСМ, наклона подложки и шероховатости. Liu et al. изучили влияние угла наклона. Износ кончика — важный фактор; смазка водяным слоем может его снизить (Ren et al.).

Для расширения канавок применяют многоходовое skreбливание. Статическая плуговая литография (контактный режим АСМ) приводит к износу при больших длинах. Динамическая плуговая литография (режим таппинга) снижает износ. Xiao et al. показали, что динамический метод позволяет изготавливать более мелкие наноструктуры. Sundaram et al. разработали метод VANILA (вибрационно-ассистированное наноимпульсное механообработка свободными абразивами), где кинетическая энергия абразивов создается вибрацией кончика АСМ.

Анализ результатов моделирования

Анализ сил резания

Сила резания отражает процесс удаления материала и эволюцию дефектов. МД позволяет отслеживать в реальном времени поперечные и касательные компоненты, а также средние и результирующие силы. Вариация сил раскрывает влияние структуры материала и параметров обработки. Li et al.: на поликристаллическом Cu силы выше из-за анизотропии зерен, на монокристалле — стабильное пластическое течение. Скорость резания увеличивает силы за счет большего объема стружки (металлы), но для SiC более высокая скорость делает материал более пластичным (аморфизация) и снижает силы. Yan et al.: билизерные Cu/Ni демонстрируют повышенные силы из-за барьера на границе раздела. Для полимеров скорость влияет через накопление материала перед кончиком, упрочнение при деформировании и термическое размягчение. Геометрия кончика: увеличение полуапексального угла растет площадь контакта и силы, коэффициент трения падает, твердость растет.

Специфические условия: водяной слой (Ren et al.) увеличивает сопротивление, а не смазывает (в отличие от макромасштаба). Наклон кончика (Liu et al.) сильнее влияет на нормальную твердость, чем на касательную силу. Динамическая плуговая литография дает колебательные силы с периодом таппинга, средняя сила в 2 раза ниже статического режима — меньше износ. Ориентация кончика влияет на эффективность: ориентация «боковая грань вперед» дает наибольшую глубину при сопоставимых усилиях.

Анализ состояния удаления материала

МД позволяет наблюдать процесс удаления в реальном времени, в отличие от пост-фактум SEM/AFM. Выделяют состояние плугования (материал выдавливается по бокам) и резания (формируется непрерывная стружка). Отношение этих состояний зависит от радиуса кончика и глубины. Анизотропия монокристалликов сильно влияет на удаление. Для поликристалликов ориентация зерен доминирует над их размером (Gao et al.). Геометрия кончика: конус — зависимость от полуугла (больше угол → больше стружка, гладче поверхность, но выше силы, трение, температура, подповерхностные дефекты). Пирамидальный кончик (реальная геометрия АСМ): сравнивают направления «ребро вперед», «грань вперед», «боковая грань вперед». Для SiC «ребро вперед» дает стабильный процесс. Полимеры чувствительны к температуре: локальный нагрев выше T_g обеспечивает вязкое удаление. Гибкость молекулярных цепей растет с температурой (Zhan et al.). Скорость обработки меняет механизм деформации полистирола: внутрицепные изменения при низких скоростях, межцепное скольжение при высоких (Du et al.).

Вода улучшает качество поверхности: уменьшает прилипание, снижает сдвиговые силы на границе стружка-инструмент. Новые методы: VANILA (Shockly et al.) — параметры удара определяют размер нанокавитаций; динамическая плуговая литография (Xiao et al.) — меньшие размеры канавок, управляемые амплитудой привода. Износ кончика: в большинстве работ кончик жесткий. Meng et al. использовали деформируемый кончик: износ начинается как адгезивный, переходит в атом-пo-атомный, слабо зависит от скорости. Перспективно для изучения реального износа.

Анализ эволюции дефектов

Экспериментально дефекты изучают ТЭМ после ФИБ-подготовки (сложно, дорого, материалозависимо). МД дает доступ к подповерхностным дефектам и их динамике. Алгоритмы: CNA, CSP, анализ вектора сдвига, Ackland-Jones. Для трассировки дисплокаций при больших их числах предложен slipping analysis (Chen et al.), фильтрующий атомы, скользившие относительно соседей. Xiao et al. применили его к статической/динамической плуговой литографии: динамический режим — материал течет вниз и вбок, статический — нет явного нисходящего потока. Распространение дисплокаций зависит от ориентации кончика. DXA (Dislocation Extraction Algorithm) идентифицирует тип и вектор Бургера. Gao et al. для железа: перестройка дисплокационной сети, линейный рост пластической зоны, затем падение плотности дисплокаций, концентрация активности на фронте skreбливания, вакансии от реакций дисплокаций и деформационное близнецование.

Поликристаллы: деформация через нуклеацию дисплокаций в зернах и на границах. Крупные зерна — нанопустоты и интерстициальные кластеры; мелкие — близнецование на границах. Li et al. сравнили Cu: монокристалл — взаимодействие дисплокаций; поликристалл — дисплокации + границы зерен; NT-поликристалл — дисплокации + границы зерен + границы-близнецы + близнецование/деблизнецование.

Полупроводники (Si, SiC, GaAs, AlN): кроме скольжения дисплокаций важна фазовая трансформация. Взаимодействие зависит от ориентации. Si (010) до «pop-in» проходит неупругую деформацию с переходом Si-I → Si-III/XII (для (110) и (111) нет). Доминирует переход Si-I → bct-5. Dai et al.: эволюция фаз коррелирует с гидростатическим напряжением и температурой. SiC: удаление через переход цинкбленда → аморфное + немного гексагонального алмаза. Высокая скорость → больше аморфных атомов, меньше гексагональных и дисплокаций. Meng et al.: коэффициент Шмида предсказывает скольжение 3C-SiC, если фазовый переход не доминирует. Термическое размягчение и скорость деформации управляют объемом удаления и SSD. Поликристаллический SiC: аморфный переход доминирует, на границах — гексагональный алмаз и дисплокации. Более мягкий, чем монокристалл (меньше силы, тоньше SSD). Недостаточно исследований сложных структур.

Перспективы и вызовы

1. Новые технологии TBN: нанофрезеровка на кончике АСМ (перспективна, механика не изучена), вибрационно-ассистированное наноскребливание с вибрацией образца (еще не моделировалось). МД может дать понимание силы, дефектов, удаления материала при вращении кончика.
2. Ограничения масштабов: МД не может количественно описать эксперимент из-за ограничений по длине и времени. Часто только качественное объяснение. Для полимеров (зависящих от времени) скорость зонда критична. Необходимы улучшение алгоритмов и вычислительной мощности.
3. Полимеры в TBN: много наноструктур на полимерах создано экспериментально, но МД-исследований мало. Термическое сканирующее зондовое литографирование развивается. Нужно раскрыть механизм удаления при термомеханическом воздействии.

Заключение

МД-моделирование играет все большую роль в нанообработке. Основные выводы обзора:

1. Обобщены методы построения МД-моделей различных материалов и потенциалы взаимодействия. Детально разобраны NT-поликристаллы и аморфные полимеры. Точность модели критична для результатов.
2. Рассмотрены новые технологии TBN: многоходовое skreбливание, динамическая плуговая литография, VANILA. МД раскрыло принципиальные отличия статического и динамического режимов по внутренним дефектам, морфологии и силам.
3. Суммирован анализ МД: силы резания, состояние удаления материала, дефекты. Введены новые методы анализа (slipping analysis). Показана материалозависимость механизмов TBN.

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

TBN

Tip-based nanomachining — нанообработка на кончике зонда

AFM

Atomic force microscopy — атомно-силовая микроскопия

MD

Molecular dynamics — молекулярная динамика

NT

Nanotwined — наноблizkированный

PE

Polyethylene — полиэтилен

EAM

Embedded atom method — метод вложенного атома

SW

Stillinger-Weber — потенциал Стиллингера-Вебера

ABOP

Analytical bond order potential — аналитический потенциал порядка связи

REBO

Reactive Empirical Bond Order — реактивный эмпирический порядок связи

AIREBO

Adaptive Intermolecular Reactive Empirical Bond Order — адаптивный интермолекулярный REBO

LAMMPS

Large-scale Atomic/Molecular Massively Parallel Simulator

MS

Material Studio

OVITO

Open Visualization Tool

VMD

Visual molecular dynamics

VANILA

Vibration-Assisted Nano Impact machining by Loose Abrasives

SEM

Scanning electron microscope — сканирующий электронный микроскоп

CNA

Common neighbor analysis — анализ общих соседей

CSP

Centro-symmetry parameter — параметр центрально-симметричности

DXA

Dislocation extraction algorithm — алгоритм извлечения дисплокаций

GB

Grain boundary — граница зерна

TB

Twin boundary — граница-близнеца

SSD

Subsurface damage — подповерхностное повреждение