1. 课程导入与工程视角

金属材料在宏观上表现出的强度、塑性、疲劳寿命与抗蠕变能力，几乎都可以追溯到微观尺度上两类缺陷的行为：位错与晶界。位错是晶粒内部的线缺陷，是塑性变形的载体；晶界是不同取向晶粒之间的界面，是位错运动的障碍、吸收体和发射源。工程上常说的"细晶强化""加工硬化""晶界偏聚导致脆化"，本质上都是这两类缺陷在原子尺度上相互作用的结果。本课程的目标，是建立一套用分子动力学（MD）研究晶界与位错行为的工程化工作流，使读者能够从原子构型出发，定量地观察、计算并解释这些微观机制。

为什么必须用原子尺度方法

连续介质力学可以很好地描述位错弹性场和晶界能，但它无法回答几个关键问题：位错核在原子尺度如何重构？晶界在受载时是先发射位错还是先发生迁移？溶质原子偏聚到晶界后，界面结合强度如何变化？这些问题涉及原子重排、键断裂与局部能量起伏，必须借助能显式处理原子间相互作用的模拟方法。分子动力学通过求解牛顿运动方程，在飞秒时间步长上追踪每个原子的轨迹，配合经验势函数（如嵌入原子法EAM、MEAM或机器学习势），可以在百万原子规模、纳秒时间尺度上给出位错与晶界的动态演化图像。

需要明确的是，MD不是"万能显微镜"。它的时间尺度远小于实验蠕变，空间尺度远小于工程构件，势函数的精度直接决定结论可靠性。因此工程视角下的MD工作流，核心不是"跑一个模拟"，而是把模拟设计成一个可验证、可复现、可解释的数值实验：明确初始构型、边界条件、加载方式、观测量与判据，并始终与理论或实验做交叉校验。

本节核心关系图

下图给出本课程的整体逻辑：从工程性能需求出发，向下分解到位错与晶界两类缺陷，再进入MD工作流，最终回到性能解释与设计建议。图中同时标出了各环节的关键输入与输出，便于读者建立全局认知。

MD工程工作流的关键环节

上图把工作流拆成五个环节，每个环节都有明确的工程含义。构型构建阶段，需要根据研究目标选择双晶模型（含特定取向差与晶界面）或位错偶极模型（含特定Burgers矢量与滑移面）。以面心立方金属为例，构建Σ5(310)对称倾斜晶界时，取向差约36.87°，晶界面法向为[310]，这些参数直接决定后续晶界能与位错交互行为。势函数选择阶段，必须先用已知实验值校验：例如Cu的EAM势应能复现晶格常数3.615 Å、弹性常数C11≈170 GPa、层错能约45 mJ/m²；若层错能偏差过大，位错分解宽度和交滑移行为就会失真。

弛豫与加载阶段，通常先做能量最小化，再在目标温度下NPT弛豫使系统达到热平衡，随后施加剪切或单轴应变。应变率在MD中通常取10⁸～10¹⁰ s⁻¹，远高于实验值，这是尺度限制带来的必然妥协，解释结果时必须说明。轨迹分析阶段，常用共同邻域分析（CNA）、位错提取算法（DXA）和PTM识别局部结构，配合应力-应变曲线判断屈服、位错形核与晶界迁移的临界条件。最后是机制解释：把临界应力、激活能垒与理论模型（如Hall-Petch、Orowan、晶界滑动模型）对比，判断模拟结果是否物理合理。

学习目标与工具链概览

完成本课程后，读者应能独立完成以下任务：根据工程问题设计合理的双晶或位错构型；选择并验证势函数；设置弛豫、加载与边界条件；使用OVITO或Python脚本提取位错密度、晶界迁移距离、应力-应变响应等关键量；对结果做物理合理性判断并撰写分析报告。工具链以LAMMPS为核心计算引擎，OVITO负责可视化与结构分析，ASE用于构型生成与格式转换，Python（NumPy、Matplotlib、pymatgen）承担数据处理与绘图。这套组合在学术与工业界都足够成熟，且全部开源，便于复现与二次开发。