第一性原理计算层错能：超胞建多大、层错怎么引入、弛豫要不要开

层错能（SFE）是fcc合金变形机制的核心参数——SFE > 20 mJ/m²，位错交滑移容易，材料趋向波状滑移和位错胞结构；SFE < 20 mJ/m²，全位错倾向于分解成Shockley不全位错，堆垛层错变宽，TWIP和TRIP效应被激活。

一个做高锰TWIP钢的团队用DFT算出了Fe-22Mn-0.6C合金的层错能是42 mJ/m²，但实验从XRD层错概率反推的SFE只有约25 mJ/m²。差异超过了50%——不是DFT的精度问题，是计算方案决定了结果。

超胞法引入层错：周期镜像的鬼影作用

DFT算层错能的标准做法是超胞法（supercell approach）：建一个足够大的超胞，沿[111]方向（fcc的密排方向）引入一个本征层错——即抽掉一层(111)原子面，然后让上下两半晶体相对滑移a/6[11-2]。

引入层错后的超胞包含两个层错面——一个是真实引入的层错，另一个是超胞另一端的周期镜像——因为超胞在[111]方向具有周期性。如果超胞不够大，两个层错通过周期边界条件发生相互作用。

这个项目第一次用了一个6层密排面的超胞，层错间距只有3个密排面间距（约6.9Å）。算出来的层错能是42 mJ/m²。扩展超胞到9层、12层，层错能分别是35和29 mJ/m²。9到12层之间SFE降了6 mJ/m²——还在变！拉到15层后，SFE稳定在26 mJ/m²——从6层的42到15层的26，降了16 mJ/m²，接近40%。

层错与层错的周期镜像之间的弹性相互作用是以1/d²衰减的——d是层错间距。衰减慢——需要至少10层以上的密排面间距才能把镜像作用降到~2 mJ/m²以下。

结论：SFE < 50 mJ/m²的体系，超胞在层错方向至少需要12个密排面（即包含至少10个密排面间距的真空隔离）。6层超胞是完全不够的——但6层恰好是很多DFT计算中”看起来合理”的超胞选择。

是否弛豫原子位置

6层超胞不弛豫→42 mJ/m²；15层超胞不弛豫→26 mJ/m²。如果在这个基础上进一步弛豫原子位置（ISIF=2在层错面附近做离子弛豫），层错能从26降到了23 mJ/m²。

弛豫的贡献约3 mJ/m²——不算大，但方向一致：弛豫总是降低层错能。因为在VASP弛豫中，层错附近的原子会调整位置来释放积攒在错排面上的局部应力，这种弛豫能量增益是实验真实存在的——XRD层错概率法的SFE包含了弛豫效应。

如果你只算刚性层错（不弛豫），算出来的SFE会系统性偏高。

Fe-Mn-C体系的化学复杂性

层错能在Fe-Mn-C合金里不仅是几何问题——还是热力学问题。Mn和C在fcc Fe基体中的分布不是均匀的——层错面的局域化学组成和基体的名义组成不同（Suzuki效应：溶质原子偏聚到层错面降低层错能）。DFT的超胞通常假设溶质原子均匀分布——这是DFT层错能和实验值之间天然存在的系统性差距。

这个项目里，针对Fe-22Mn-0.6C（原子分数：Mn 20at%，C 2.5at%），用SQS（Special Quasirandom Structure）构建了化学无序的fcc超胞。16个Fe原子+4个Mn+1个C，穷尽了尽可能多的SQS排列，层错能的变化范围在21-32 mJ/m²之间。这11 mJ/m²的散布反映的是化学构型熵对层错能的统计影响——实验上，XRD层错概率测量的大量晶粒天然做了构型平均。DFT单次计算拿的只是一个微观构型的SFE。

均值26 mJ/m²和实验25 mJ/m²对上了，但散布的11 mJ/m²也告诉了你：这个成分区间的层错能对局域化学起伏极度敏感——TWIP效应在高Mn钢中有时效性变化，本质就是层错能的化学散布问题。