纳米压痕分子动力学模拟中势函数的选择：EAM与MEAM对位错形核与力学响应的差异化影响

纳米压痕模拟的原子尺度精度高度依赖于势函数对缺陷能量的描述能力。势函数选择在弹性变形阶段似乎无关紧要，但当位错形核成为主导机制时，EAM与MEAM之间的差异从隐蔽变成了决定性。

弹性段的虚假共识

单晶铜(111)面的纳米压痕模拟中，EAM势函数（Mishin等，2001）给出的力-位移曲线在弹性段与Hertz接触理论吻合良好。压深达到8Å时，载荷-压深关系的斜率与理论预测的偏差小于3%。这个阶段容易让人产生一种错觉：EAM已经足够精确。

MEAM势函数在弹性段表现出几乎相同的行为。两种势函数给出的弹性接触刚度差异不到2%，这个差异远小于模拟本身的统计涨落。但弹性段的一致性掩盖了塑性变形机制的系统性偏差——EAM对堆垛层错能（SFE）的描述偏高，导致其在预测初始塑性事件（pop-in）时发生延迟。

Pop-in载荷的系统性偏移

Pop-in现象对应着力-位移曲线上载荷的突然下降，这是位错首次从压痕区发射的宏观信号。EAM势下，单晶铜(111)面的pop-in载荷约为1.8μN（压头半径5nm，加载速率0.5m/s）。相同的模拟条件下，MEAM势给出的pop-in载荷为1.4μN，与实验观测值（1.2-1.5μN，Atomic Force Microscopy测定）更为接近。

这个0.4μN的差异并非偶然。EAM势的高堆垛层错能使得Shockley偏位错的解离宽度被低估，位错形核需要克服更高的能量势垒。MEAM势通过包含角依赖的第二近邻相互作用，将堆垛层错能的预测值从EAM的48mJ/m²降低到32mJ/m²，更接近实验值（约35-40mJ/m²）。

加载速率对pop-in载荷的影响在两种势函数下表现出不同的敏感度。当加载速率从0.5m/s提高到5m/s时，EAM的pop-in载荷增加了65%，而MEAM下仅增加38%。这种差异来源于两种势函数对位错形核激活体积的描述不同——EAM预测的值偏小，导致应变率敏感性被高估。

位错演化路径的分歧

DXA（Dislocation Extraction Algorithm）位错分析揭示了两种势函数给出的位错演化路径的根本差异。EAM势下，首次pop-in后发射的位错主要以完美的1/2<110>全位错形式存在，Shockley偏位错的占比不到20%。这表明EAM势低估了堆垛层错的形成倾向。

MEAM势下，首次pop-in后约60%的发射位错为Shockley偏位错，这与TEM原位观测到的位错结构一致。偏位错环在(111)面上的扩展宽度在MEAM下约为8.5Å，与实验测量的9-10Å较为接近。EAM下这个宽度仅为5Å，偏位错的扩展被过度约束。

位错增殖的动力学行为也呈现出系统性差异。EAM势下，压深达到15Å后位错密度从约10¹⁶ m⁻²跃升至5×10¹⁶ m⁻²，跃升幅度过大，反映了EAM对位错相互作用的描述过于简化。MEAM势下的位错密度增长更为平缓，从10¹⁶ m⁻²逐步增加到3×10¹⁶ m⁻²（压深25Å），更符合原位TEM观察到的渐进式位错增殖过程。

压头半径的放大效应

压头半径从5nm增加到10nm时，EAM与MEAM之间的差异被进一步放大。大压头下的应力场更接近连续介质解，理论上应该降低对原子尺度势函数的敏感度。但实际情况恰恰相反：压头半径10nm时，EAM与MEAM给出的压痕硬度差异从5nm时的12%扩大到21%。

这种放大效应来源于大压头下塑性区的空间扩展。塑性区直径从5nm压头下的约15nm扩展到10nm压头下的约30nm，更多的原子参与到位错相互作用中。EAM势对位错-位错相互作用的过度简化在大规模塑性变形中累积成系统性偏差。MEAM势通过更精确的描述，在塑性区扩展后仍能保持合理的位错结构预测。

Pile-up形貌的定量差异

压痕周边的pile-up（材料隆起）形貌是评估势函数性能的另一个定量指标。EAM势下，单晶铜(111)面纳米压痕的pile-up高度为压深的18%（压深20Å），隆起区域的径向扩展约为压头半径的1.5倍。MEAM势给出的pile-up高度为压深的24%，隆起区域的径向扩展约为压头半径的1.8倍。

实验观测的pile-up高度通常为压深的22-28%，EAM的预测明显偏低。这种偏差来源于EAM势对交叉滑移能量的高估——位错更容易从原始滑移面交叉滑移到其他滑移面，导致塑性变形在深度方向过度集中，表面隆起的幅度被压缩。

势函数选择的决策框架

纳米压痕模拟中，势函数选择不应仅基于弹性常数的匹配程度。弹性段的一致性具有欺骗性，真正的差异隐藏在塑性变形的初始阶段。EAM势函数在仅需弹性响应或定性趋势的研究中仍具价值，但当位错形核的具体细节成为核心关注点时，MEAM提供了不可替代的精度。

计算成本的权衡同样关键。MEAM势的计算量约为EAM的3-5倍（取决于截断半径和近邻列表设置）。在需要遍历大量压痕参数的筛选性研究中，EAM提供的定性趋势可能已足够；但在需要定量化位错结构的机理研究中，MEAM的额外开销无法省略。

压头几何形状的复杂性进一步影响势函数的选择。球形压头下EAM与MEAM的差异相对较小（约10-15%），但尖锐压头（如Berkovich压头）下差异可超过30%。尖锐压头产生的高应力梯度对势函数的高阶弹性常数更敏感，EAM在这些条件下的预测可靠性显著下降。

欢迎访问 keyanxueshu.com 了解分子动力学计算服务。

参考文献

1. Mishin, Y., et al. “Structural stability and lattice defects in copper: Ab initio calculations and embedded-atom theory.” Physical Review B 63, 224106 (2001).
2. Lee, B.-J., & Baskes, M. I. “Second nearest-neighbor modified embedded-atom-method potential for bcc transition metals.” Physical Review B 62, 8564 (2000).
3. Tsuru, T., & Shibutani, Y. “Anisotropic effects on pop-in and incipient plasticity in nanoindentation of cubic crystals.” Journal of Materials Research 23, 1251-1262 (2008).