如何处理二维异质结的晶格失配？

引言

二维异质结最容易被低估的地方，不是把两层材料叠起来，而是叠起来以后到底承受了多少应变、形成了多大的莫尔周期、层间距是否落在合理势能谷里。早期做几何匹配时，常把失配率小于某个数值当作结束；进入 VASP 计算前，还要把旋转角、超胞规模、应变分配和界面距离放到同一套判断中。晶格失配不是单个数字，而是超胞取舍、层内应变和层间相互作用共同决定的模型误差。

对已经能找到候选匹配结构的体系，更关键的问题会转向模型取舍：旋转角改变的是局部堆垛和莫尔周期，超胞尺寸决定计算成本与周期镜像距离，界面距离则需要通过能量扫描和弛豫确认。只有这三项同时可解释，后续吸附能、功函数、能带对齐和差分电荷密度才有可比性。

图1 低应变 MoTe₂/ZrS₂ 旋转异质层的结构、能带、能级对齐和差分电荷密度。图中展示了旋转匹配后异质层仍需同时检查几何、能带和界面电荷。DOI: 10.1021/acsami.2c13151。

晶格失配先看哪个量？

两个二维晶格能否放在同一个周期边界里，首先取决于面内基矢是否能通过整数倍复制和旋转近似重合。对六方、正交或低对称二维材料，单独比较晶格常数 a 往往不够，因为角度和 b 方向也会参与误差。更可靠的做法是把一层作为参考层，另一层经过旋转矩阵 R(θ) 后，在整数矩阵 m、n 下寻找最小形变。

δ = |A₁ − R(θ)A₂| / |A₁| × 100% 式(1)

式(1)里的 A₁ 和 A₂ 可理解为两个候选超胞的面内矢量集合。δ 不是只看一个方向的长度差，而应同时检查两个方向的长度、夹角和面积变化。若只让其中一层承担全部应变，后续吸附能和能带对齐会混入人为拉伸或压缩效应。

经验上，小于 1% 的面内应变更容易作为严谨比较的起点；1% 到 3% 之间要说明应变分配方式，并检查关键性质是否敏感；超过 3% 的模型仍可用于趋势探索，但不宜直接和未应变单层的能带、功函数或吸附自由能做定量比较。这些阈值不是通用定律，仍需结合材料刚度和目标能量差确定。

旋转角怎样影响超胞尺寸？

旋转角的作用不只是改变视觉上的堆叠取向。角度一变，两个晶格的可公度周期、局部堆垛区域和莫尔长度都会改变。小角度扭转常带来很大的超胞，原子数迅速增加；大角度候选虽然更便宜，却可能对应完全不同的局域堆叠环境。选择旋转角时，应同时记录角度、原子数、应变张量和界面堆垛类型。

实际筛选可先生成一组候选：0°、30° 或材料高对称角附近的低原子数方案，再加入若干低应变旋转方案。每个候选不要只看失配率，还要看是否保留主要研究对象。例如研究界面电荷转移时，层间接触面积和局部注册关系比超胞是否最小更重要；研究吸附时，吸附位点周围的局域堆叠要能代表真实界面。

图2 边际扭转 MoSe₂/WSe₂ 异质双层中的莫尔应变图和局域堆垛区域。图中显示旋转角会改变局域应变分布和电子限域环境。DOI: 10.1021/acs.nanolett.3c04427。

当旋转角导致上百到上千原子的超胞时，完整结构优化、HSE06 或 CI-NEB 通常不现实。这时可以采用两级策略：先用较小公度超胞建立基准趋势，再对关键角度做更大模型的静态或局部弛豫验证。超胞越大并不自动更可靠，只有当它降低的失配误差大于新增的采样和收敛误差时才值得使用。

超胞尺寸怎样取舍？

超胞尺寸的取舍要同时满足三件事：面内应变可解释、原子数可承受、目标性质不被周期镜像主导。若只是计算裸界面结合能，可以优先选择低应变、低原子数候选；若要放入吸附物、缺陷或单原子位点，还要确保这些局域对象之间的周期距离足够大。否则看似在研究异质界面，实际可能是在研究一个高浓度周期阵列。

E_str = (E_A^str − E_A⁰) + (E_B^str − E_B⁰) 式(2)

式(2)把两层被迫拉伸或压缩的能量代价单独列出来。若某个匹配方案的界面结合能看起来更负，但应变能也显著增大，就不能直接判定该堆叠的稳定性更高。二维异质结比较时，界面结合能、应变能和弛豫后结构应分开记录。

VASP 中可以先固定面内晶格，分别优化两个被应变后的单层，再优化异质层。这样得到的能量差能区分层内应变代价和层间结合贡献。若直接用自由单层能量作为参照，式(2)中的应变代价会被隐藏，界面稳定性容易被高估。

图3 MoX₂/WX₂ 异质双层中扭转角与静水应变的关系。图中说明角度变化会改变局部应变大小，不能只用单个失配率描述模型。DOI: 10.1021/acs.nanolett.3c04427。

界面距离怎样确定？

界面距离不能只靠手动放一个 3 Å 或 3.5 Å 的初值。不同二维材料的层间相互作用、端基、极性和 vdW 修正方式会改变平衡距离。比较可靠的流程是先做层间距扫描，固定面内晶格和层内构型，改变层间距离 d，得到 E(d) 曲线；再从能量谷附近选 2 到 3 个初始值做完整弛豫。

E_int(d) = E_AB(d) − E_A^str − E_B^str 式(3)

式(3)中的 A 和 B 应使用与异质层相同面内晶格的应变单层。d 的平衡值不是只看最小能量，还要检查弛豫后是否出现层内异常翘曲、原子穿插或局部键化。对 vdW 异质结，D3、D4、optB86b-vdW 或 optB88-vdW 等口径要在同一组比较中保持一致。改变 vdW 修正后重新扫描 d，比只继承旧距离更可靠。

图4 二维材料层间结合力和结合能的计算流程。图中以层间距离 d 为变量寻找结合能谷，为界面距离扫描提供了可执行思路。DOI: 10.34133/2022/9765121。

若异质结后续要放吸附物，界面距离还会影响吸附位点的局域电场、功函数差和电荷转移。此时建议至少保留裸异质层、吸附前应变单层和吸附后异质层三套结构，确认吸附能变化不是由界面距离未收敛造成的。

VASP 任务怎样安排？

一个可执行的顺序是：先用脚本枚举旋转角和整数超胞，输出每个候选的原子数、应变张量和面积差；再对 3 到 5 个候选做单层应变能评估；随后扫描界面距离，选择能量谷附近结构做完整优化；最后只对通过几何和能量筛选的结构计算吸附能、功函数、差分电荷密度或能带。

INCAR 参考示例 ENCUT = 500 EDIFF = 1E-5 EDIFFG = -0.02 IBRION = 2 NSW = 150 ISMEAR = 0 SIGMA = 0.05 IVDW = 12 LDIPOL = .TRUE. IDIPOL = 3 LREAL = Auto

上面的参数只是参考示例，ENCUT、k 点、vdW 修正和偶极修正都要按体系测试。二维异质层通常使用面内 k 点、真空方向 1 个点；大超胞可先用 Γ 点优化，再做更密 k 点静态对照。不同旋转角之间比较时，POTCAR、ENCUT、vdW 修正、偶极修正和能量参照必须保持同一口径。

图5 MoTe₂/ZrS₂ 异质层在外场和面外应变下的电荷转移与结合能变化。图中体现了界面距离和面外作用会改变层间电荷响应。DOI: 10.1021/acsami.2c13151。

候选结构的筛选应先排除应变过大和原子数不可承受的方案，再用 E_str 与 E_int(d) 区分层内代价和层间结合，并用目标性质检验不同候选是否给出相同趋势。若两种匹配方案的吸附能差小于 0.05 eV，但应变状态和界面距离不同，不宜把差异解释成材料本征活性变化。

本章要点总结

要点：晶格失配应同时检查长度、夹角、面积和应变分配，不能只看一个百分数。

要点：旋转角会改变超胞规模、局部堆垛和莫尔应变，候选结构要按目标性质筛选。

要点：界面距离应通过 E(d) 扫描和弛豫共同确认，vdW 修正口径必须一致。

要点：异质结能量比较应拆分应变能、界面结合能和吸附后结构变化。