晶体分子动力学模拟 — 晶界建模与位错-缺陷相互作用的原子级分析

晶体分子动力学模拟和生物系统的MD有本质的区别：生物体系关心的是氢键网络和疏水接触——作用力在kcal/mol量级；晶体体系关心的是位错运动和晶界迁移——作用力在eV/原子量级，和键能同一个量级。这意味着力场精度要求更高，势函数必须能够捕捉键的断裂和重建、非简谐晶格振动、以及不同晶格缺陷之间的长程弹性相互作用。本文围绕晶体MD的核心——晶界建模、位错分析、和缺陷演化展开。

一、双晶模型与晶界构建

晶体分子动力学模拟中最常用的配置是双晶（bicrystal）模型——两个不同取向的完美晶体接触在一起形成一个晶界。晶界建模的基本流程：先用LAMMPS或Atomsk生成两个各自取向的完美晶体（各自有自己的一套晶格取向基矢），将它们接合——界面上原子位置重叠的删除，原子间距过近的调整——形成初始晶界构型。

晶界的重合位置点阵（CSL，Coincidence Site Lattice）理论描述了晶界的几何特征。Σ3（60°绕<111>旋转）是FCC中的低能孪晶界——密排面上有一半原子重合。Σ5（36.87°绕<100>）是非孪晶低能晶界，但晶界能高于Σ3。

晶界能γ_GB的计算：γ_GB = (E_GB – N×E_bulk)/A，其中E_GB是含晶界体系的总能量，N是原子数，E_bulk是完美晶体的原子内聚能，A是晶界面积×2（体系有两个晶界面——周期性边界）。对于FCC铝的Σ3孪晶界，γ_GB约0.08-0.12 J/m²；对于Σ5非孪晶界，约0.3-0.5 J/m²。晶界能可以用MD做能量最小化（0 K弛豫）获得——高温的振动熵贡献需要通过准简谐近似或热力学积分来补充。

二、位错分析的核心工具

晶体分子动力学模拟中位错的识别和表征需要专门的算法——不能用肉眼在原子快照中找。几种主流分析工具：

位错提取算法（DXA，Dislocation Extraction Algorithm）：从原子构型中通过Burgers回路分析自动检测位错线位置、Burgers矢量和位错类型（刃型、螺型、混合型）。DXA输出的是位错网络的拓扑结构——每条位错线的长度、滑移面、和交点。OVITO中的DXA插件是晶体MD的最核心可视化/分析工具。

共近邻分析（CNA，Common Neighbor Analysis）：基于原子对的近邻拓扑关系来判断原子局部环境——FCC、BCC、HCP或其他。位错核心处的原子通常不属于任何理想晶格类型——CNA值标记为”Other”——使得位错核心位置在可视化中一目了然。

中心对称参数（Centrosymmetry Parameter，CSP）：专门为FCC设计的缺陷识别参数，利用FCC晶格的12个最近邻原子的反演对称性来计算每个原子的CSP值。完美FCC原子的CSP≈0，位错核心>5，自由表面>20。CSP的计算速度远快于CNA——对于百万原子级的大体系尤其明显。

三、纳米压痕模拟

晶体分子动力学模拟中纳米压痕是研究材料塑性变形微观机制的标准MD实验。用LAMMPS的`fix indent`命令施加一个虚拟的球形压头，以恒定速度（v_indent）压入晶体表面——载荷-位移曲线和位错-晶界相互作用的原子级现场。

压痕的载荷P和压痕深度h的关系对于FCC金：P∝h²——这是Hertz弹性接触理论+弹塑性过渡的特征。载荷-位移曲线上出现的”pop-in”（位移突然跳变）对应着位错的首次形核——从完美晶格中形核第一个位错环所需的应力接近理论强度（约G/10，G为剪切模量）。在MD中看到的pop-in载荷通常是300-500 nN量级（对于10 nm半径压头）。

纳米压痕MD的精度陷阱：压头半径和压痕深度必须远大于晶格常数才有连续体力学意义。r_indent>10a（a=晶格常数），h_max<r_indent/2——对于FCC铜（a=0.361 nm），压头半径至少3-5 nm。如果r_indent只有2 nm——压头只接触了不到100个表面原子——不可能得到连续的弹塑性响应，pop-in是离散原子重排的噪声而非真正的位错形核。

四、辐射损伤MD模拟

晶体分子动力学模拟在核材料研究中的代表应用是级联碰撞（Cascade Collision）——一个高能粒子（如中子或重离子）撞击晶体，产生一连串的原子位移级联，最终形成一个包含空位、间隙原子、位错环的损伤区。

级联碰撞MD的设置：从晶格中随机选一个原子作为初级撞击原子（PKA，Primary Knock-on Atom），给它一个动能E_PKA（对于裂变中子，典型平均E_PKA约10-30 keV，可通过SRIM/TRIM程序计算）。PKA以这个速度飞出去——撞击沿途原子形成碰撞级联。级联持续时间<10 ps（弹道相）→10-100 ps（热峰相，局部温度可达5000 K以上）→>100 ps（冷却和缺陷稳定相）。

级联碰撞MD的核心输出：存活缺陷数（Frenkel pairs）与PKA能量的关系。在能量<1 keV时缺陷数线性增长（每个级联产生几个Frenkel pair）；在>10 keV时缺陷数与PKA能量的非线性次幂增长（约N_FP∝E_PKA^0.7-0.8），因为位错环和缺陷团簇的形成改变了缺陷复合的概率。

五、势函数选择

晶体分子动力学模拟的势函数是精度的第一决定因素。经典选择速查：FCC金属→EAM（Embedded Atom Method）势，由Foiles、Mishin等开发；BCC金属→EAM也可用，但注意BCC铁中螺位错的非Schmid行为——不是所有EAM势都预测准确；共价晶体（Si/SiC）→Tersoff势或Stillinger-Weber势；金属-非金属混合体系→MEAM（Modified EAM）或COMB（Charge-Optimized Many-Body）势。

第一原理精度势（如SNAP、GAP、DeepMD、NequIP等机器学习势）正在逐步替代经典势，精度接近DFT、速度接近EAM。SNAP势用一个展开在4D超球面上的双光谱分量来描述每个原子的局部环境——拟合到DFT数据的精度在力误差约50 meV/Å量级。对于需要描述键断裂/形成或强非简谐效应的场景——机器学习势的精度优势显著。

六、专业晶体分子动力学模拟服务

 需要晶体分子动力学模拟服务？

科研学术网提供专业的晶体MD模拟分析服务：

– ✅ 晶界与界面：双晶建模、晶界能计算、晶界迁移/滑动动力学

– ✅ 力学行为：纳米压痕、单轴拉伸/压缩、循环加载的位错演化

– ✅ 辐射与缺陷：级联碰撞、缺陷团簇演化、氦泡/空洞生长

– ✅ 高精度势函数：EAM/MEAM经典势+SNAP/DeepMD机器学习势

立即咨询报价 →