异质结界面建模中的晶格失配处理：supercell匹配、应变调节与缓冲层三种策略的选择

异质结界面建模中最让人头疼的问题不是能带排列的计算方法，而是两种材料拼在一起时那百分之几到十几的晶格失配。这个数字看起来不大，但它直接决定了后续所有计算的可靠性——忽略它，界面处的原子排布就是错的，基于错误结构算出的功函数差、电荷转移量和界面态密度全部失去意义。

项目组在搭建一个MoS₂/WS₂范德华异质结模型时首次深刻体会到了这个问题的重要性。

MoS₂/WS₂：小失配下的supercell方案

单层MoS₂的晶格常数约为3.160 Å，而单层WS₂约为3.153 Å。两者之间的晶格失配仅为0.22%，属于极小的范畴。对于这种程度的差异，最直接的策略是构建一个公度超胞：寻找两个晶格的最小公倍数关系，使得超胞内的两种材料都能保持接近其平衡晶格常数。

具体操作中，通过枚举MoS₂和WS₂的超胞尺寸组合（如√3×√3 R30°重构），可以找到一个(4×4)MoS₂/(4×4)WS₂的匹配超胞。在这个超胞中，两种材料的面内晶格参数被强制取为两者的加权平均值（约3.1565 Å），各自偏离其自由-standing值的幅度均小于0.12%。这种微小的应变对电子结构的影响在DFT计算精度范围内完全可以接受。

计算结果显示，采用该超胞匹配方案的异质结界面建模给出的层间结合能约为18 meV/Ų，与实验剥离能数据吻合良好。Type-II型的能带排列特征清晰可见：MoS₂的导带底低于WS₂的导带底，而WS₂的价带顶高于MoS₂的价带顶，光激发后电子和空穴将分别位于不同层——这正是这类异质结作为光催化材料的核心优势所在。

但supercell方案有一个明显的代价：超胞的原子数会迅速膨胀。上述(4×4)/(4×4)构型已经包含了96个原子（含真空层）。如果需要引入缺陷或掺杂进行更精细的研究，计算成本将急剧上升。这是每个采用该方法的研究者都必须面对的现实权衡。

石墨烯/hBN：中等失配下的应变策略

当异质结界面建模的对象换成石墨烯/hBN时，情况变得复杂得多。石墨烯的晶格常数为2.461 Å，hBN为2.501 Å，失配率达到了1.6%。这个数字仍然不算大，但已超出简单平均法的舒适区——如果强行取平均值施加于两者，石墨烯将被拉伸约1.6%，hBN被压缩约1.6%。对于石墨烯而言，1.6%的双轴拉伸会使其狄拉克点处产生约0.2 eV的带开口，这直接改变了材料的低能电子性质。

实际项目中采用了分层应变策略来处理这一困境：

方案A是将整个体系固定在石墨烯的晶格常数上。此时hBN承受约-1.6%的压缩应变。由于hBN的弹性模量较大（约0.77 TPa），这种压缩对应的应变能在可接受范围。计算表明，压缩后的hBN带隙从自由-standing时的5.97 eV略微增大至6.12 eV，变化幅度小于3%。

方案B则是将体系固定在hBN的晶格常数上。此时石墨烯承受约+1.6%的拉伸应变。如前所述，这将打开一个约0.2 eV的带隙并使Dirac锥移动至K点以外。对于研究石墨烯输运性质的工作来说，这种改变可能是不可接受的。

方案C是取中间值（如2.48 Å）作为共同晶格常数。此时双方各承担约一半的应变，形成一种”各退一步”的妥协格局。这种方法在实际应用中被广泛采纳，因为它避免了任何一方承受过大的畸变。

经过对比验证，方案C在该体系中给出了最合理的综合结果：石墨烯的线性色散保持完好（残余曲率可忽略），hBN的绝缘体特性未受显著影响，且界面处的莫尔条纹周期与扫描隧道显微镜观测结果一致[1]。

大失配体系的缓冲层策略

当晶格失配超过5%时，前述两种策略都会遇到瓶颈。supercell方案需要的超胞尺寸可能大到无法计算（比如需要几十乘几十的原胞才能找到合理的匹配），而均匀应变策略则会导致至少一方承受不可接受的畸变（>3%的应变通常会使键长偏离平衡值10%以上，严重破坏局域化学环境）。

以一个假想的Si/GaAs异质结为例。Si的金刚石结构晶格常数为5.431 Å，GaAs闪锌矿结构为5.653 Å，失配率约4.1%。虽然勉强还在可控范围内，但如果进一步考虑取向差异（如Si(001)/GaAs(001)），问题就变得更加棘手。

在这种场景下，引入缓冲层的策略开始展现其价值。基本思路是在两种主体材料之间插入一层或多层过渡性材料，这些过渡材料的晶格常数介于两端之间，起到”桥梁”的作用。对于Si/GaAs体系，常见的缓冲层选择包括渐变组分Si₁₋ₓGeₓ合金或超薄GaP插入层。

在DFT模拟层面，缓冲层策略的实现方式是在超胞模型中包含若干层缓冲层原子。这些原子的位置需要在初始阶段进行充分的几何优化（离子弛豫步数应不少于100步，力收敛标准设为0.01 eV/Å），以确保界面处的应力得到充分释放。计算完成后，应单独检查缓冲层区域的应力张量分布——如果某个方向的残余应力超过5 GPa，说明缓冲层厚度不够或材料选择不当。

值得注意的是，缓冲层策略在二维材料异质结中的应用相对较少。这是因为二维材料可以通过弯曲而非应变来适应失配，且其超胞匹配的灵活性远高于三维块材。但在涉及二维/三维混合界面（如二维材料生长在衬底上）的问题中，缓冲层的概念依然具有重要参考价值[2]。

决策流程图

总结以上分析，不同晶格失配程度下的最优策略选择可归纳如下：

失配率<1%：直接采用加权平均晶格常数，无需特殊处理。supercell方案可选但非必需。

失配率1%-3%：优先考虑中间值应变策略。若研究重点偏向其中某一方的电子性质，可选择将该方固定在其本征晶格常数上。

失配率3%-5%：进入临界区域。应同时测试supercell匹配和中间应变两种方案，对比关键物理量（如带隙、功函数差、界面偶极矩）的敏感性后再做决定。

失配率>5%：必须考虑缓冲层方案或重新评估体系选择的合理性。如果实验上确实存在如此大的失配界面，那么界面位错和弛豫效应将成为主导因素，简单的连续介质力学近似可能不再适用。

方法论反思

回顾整个过程，异质结界面建模中的晶格失配处理没有万能钥匙。每种策略都有其适用边界和隐含假设：supercell匹配假设了完美的公度性，应变调节忽略了非均匀分布的局部应力，缓冲层策略则需要额外的材料参数输入。真正有价值的不是找到”正确答案”，而是清楚地知道每种方法的局限在哪里，以及所选方法对最终结论的影响是否在可量化的范围内。

了解更多计算服务，请访问 keyanxueshu.com。