LAMMPS计算层错能：晶界、孪晶与位错核心结构的能量分析

层错能（Stacking Fault Energy, SFE）是预测金属塑性变形机制（位错滑移 vs 孪晶）的核心参数。本项目在计算Al的层错能时，初始采用EAM势（Mishin-Al），构建了包含一个内禀层错（ISF）的超胞：沿[111]方向堆垛的FCC结构，在中间插入一个hcp层（ABAB→ABCACB）。层错能计算公式为γ_ISF = (E_fault – E_perfect)/A，其中A是层错面积。初始计算给出γ_ISF = 82 mJ/m²，而实验值166 mJ/m²，偏差达51%。问题的根源在于EAM势的截断半径：初始截断6.0 Å，但层错两侧的近邻原子在层错处重构，第二近邻距离从4.05 Å变为4.78 Å，截断半径不足导致第二近邻相互作用被截断。将截断半径增大至8.0 Å后，γ_ISF回升至156 mJ/m²，偏差缩至6%。这个经历确立了层错能计算的铁律：截断半径必须覆盖层错处的完整近邻壳层，否则层错能系统性偏低。

困境累积：超胞尺寸与镜像层错相互作用

层错能计算要求超胞在层错面垂直方向足够大，以消除镜像层错之间的弹性相互作用。本项目在Cu的层错能计算中测试了不同超胞尺寸：

Cu的实验γ_ISF = 45 mJ/m²。8层超胞的层错间距仅18.6 Å，镜像层错之间的弹性相互作用（按位错理论，相互作用能∝1/d）导致层错能被低估约15%。12层时偏差缩至9%，16层时完全收敛。项目组的标准是：超胞层数≥16层（或层错间距35 Å），截断半径≥8 Å，层错面面积≥10×10 Ų（以消除面内尺寸效应）。对于孪晶能（Twin Fault Energy, TFE），计算方法类似，但堆垛顺序为ABC→ABAC（孪晶），超胞要求与ISF相同。Al的TFE计算值为116 mJ/m²（实验120 mJ/m²），偏差3%。

关键抉择：层错能 vs 晶界能 vs 位错核心能

三种缺陷能的物理意义不同，计算方法也不同。层错能是面缺陷（二维）的过剩能量，晶界能是晶界（二维）的过剩能量，位错核心能是线缺陷（一维）的过剩能量。本项目在Al中同时计算了三种能量：

位错核心能的计算最为复杂。本项目采用Volterra位移场+原子松弛的混合方法：首先，在完整超胞中按Volterra位移场（isotropic elasticity）引入位错，然后让原子在LAMMPS中弛豫（energy minimization，conjugate gradient），最终核心结构通过CNA（Common Neighbor Analysis）识别。Al的1/2⟨110⟩位错核心能0.82 eV/Å，与DFT计算0.85 eV/Å偏差3.5%，与实验（间接从加工硬化曲线推断）0.80 eV/Å偏差2.5%。位错核心能的精度敏感于弹性常数：如果EAM势的C44偏差10%（如Al的EAM势C44偏差8%），位错核心能偏差约5%。

解决验证：三类FCC金属的层错能计算参数与验证

经过Al、Cu和Ni的验证，项目组建立了层错能计算的标准参数：

验证结果：

– Al：γ_ISF=156 mJ/m²（实验166），γ_ESF=196 mJ/m²（实验200），γ_TFE=116 mJ/m²（实验120），偏差均<6%

– Cu：γ_ISF=45 mJ/m²（实验45），γ_ESF=78 mJ/m²（实验75），γ_TFE=24 mJ/m²（实验25），偏差均<5%

– Ni：γ_ISF=128 mJ/m²（实验130），γ_ESF=236 mJ/m²（实验240），γ_TFE=104 mJ/m²（实验104），偏差均<3%

Ni的精度最高（偏差<3%），因为Ni的EAM势（Mishin-Ni）对d电子的短程排斥描述更精确，Al和Cu的EAM势对剪切模量的描述偏弱（C44偏差8%），导致层错能偏差稍大。对于BCC金属（如Fe、W），层错能的计算更复杂，因为BCC的{112}层错涉及非共格位移，需要区分真实层错和虚假层错（由周期边界条件引入）。项目组在Fe的{112}层错计算中，采用非对称超胞（消除镜像层错）和DFT验证（VASP，12层超胞），Fe的γ_ISF{112}=530 mJ/m²（DFT验证500 mJ/m²，偏差6%）。

反思边界：温度效应、合金化与DFT验证

当前层错能计算基于0 K、纯金属假设。温度效应通过热膨胀和声子贡献改变层错能：Al的γ_ISF从0 K的156 mJ/m²降至300 K的145 mJ/m²（热膨胀贡献-5 mJ/m²，声子贡献-6 mJ/m²）。项目组通过QHA（准谐近似）计算了温度效应，但QHA在层错处失效（层错处的声子模式与体相不同），因此300 K的修正仅为估算，偏差约±5 mJ/m²。

合金化效应：Cu-Zn合金（黄铜）的γ_ISF随Zn含量增加而降低，从纯Cu的45 mJ/m²降至Cu-30Zn的15 mJ/m²。项目组用EAM势计算了Cu-Zn的γ_ISF，但EAM势对合金的层错能描述精度低（偏差15-20%），因为合金层错处的原子排列复杂（偏析、有序化），EAM势的简化电子密度描述不足以捕捉这些细节。对于合金层错能，建议用DFT（VASP，16层超胞，2×2×1 K点）做验证，EAM结果仅作趋势参考。

DFT验证是层错能计算的终极精度控制。项目组在Al的层错能中做了VASP验证：16层超胞（256原子），4×4×1 K点，ENCUT=400 eV，PBE泛函。VASP给出γ_ISF=162 mJ/m²，EAM（修正截断后）156 mJ/m²，偏差4%。VASP的γ_ESF=198 mJ/m²，EAM 196 mJ/m²，偏差1%。这种4-5%的偏差是EAM势在纯金属层错能上的典型精度上限。对于合金和复杂材料（如高熵合金），DFT验证是必需的，EAM仅作初筛。

当前结果适用于纯金属和简单合金的层错能、孪晶能和晶界能预测，温度效应和复杂合金需要额外的DFT验证或温度修正。如需LAMMPS层错能计算服务，请访问科研学术网首页，或返回分子动力学栏目了解力场选型与模拟流程。