二维异质结建模最容易卡在“失配率到底怎么算”。很多同学拿两个晶格常数一除，发现差 4% 或 7%，就开始纠结能不能拼。真正进入 VASP 前，问题要再细一点：匹配的不是单个 a 值，而是两层材料在平面内的两根基矢、夹角、旋转关系和最终超胞大小。异质结匹配的本质，是用可承受的应变和可计算的原子数，换一个周期性边界下能重复的界面。本文按这个顺序把建模逻辑讲清楚。
为什么不能只拿两个晶格常数相除？
A 层晶格常数是 aA，B 层晶格常数是 aB，两者差值除一下就是失配率。这个做法只适合很有限的情况，比如两层都是近似六方、夹角相同、打算做 1×1 对 1×1 的共格模型。只要遇到矩形晶胞、低对称材料或不同取向，单个 a 值就不够用了。
更稳的做法，是把每一层的平面基矢写成 a1、a2，同时看长度和夹角。超胞匹配时，真正比较的是m1a1 + m2a2 与 n1b1 + n2b2 这类组合矢量能不能接近。失配率不是材料本身唯一的属性，而是某一组超胞组合和旋转角下的结果。
图1展示的 MoS2-黑磷体系就很典型：两层材料对称性不同，匹配要同时考虑实空间超晶格和倒空间关系。若只盯原胞常数，很容易把一个可通过旋转和放大解决的问题误判为“失配太大”。
图1 MoS2-黑磷范德华异质结的超晶格、扭转角和实空间/倒空间匹配示意。
应变应该压给哪一层？
建好周期性异质结后，两层必须共享同一个平面晶胞。此时总会有人问：到底拉伸上层、压缩下层，还是两边各让一点？这不是审美问题，而是物理近似。刚性更强、实验上更接近基底的一层，通常不宜承担过大的应变；柔软层或吸附层可以适度让步，但也不能让键长偏离到不合理。
如果把 6% 的失配全部压给一层，优化后看到的能带、功函数或吸附能，可能已经不是原材料的界面效应，而是大应变材料的性质。应变分配决定了后续电子结构变化究竟来自界面耦合，还是来自人为拉伸/压缩。
式(1) εA = (Lmatch – LA) / LA × 100%
旋转角改变的是周期，还是界面堆垛？
旋转角θ 先改变的是两层基矢之间的相对方向。方向一变，原来不好凑的长度组合可能突然接近；原来 1×1 看似很美的组合，也可能因为角度不合适而产生很大的莫尔周期。对于二维异质结，θ 不只是一个几何装饰，它直接决定超胞面积、原子数和界面局域堆垛。
图2中的氧化物-TMD 异质结清楚展示了这一点：扭转角改变后，界面会出现不同尺度的莫尔图案。放到 VASP 里，这意味着同一对材料可能存在多个可行模型，0°、30° 或小角度模型回答的不是同一个界面问题。
旋转角的作用，是在应变、原子数和局域堆垛之间重新分配代价。因此建模时不要只问“哪个角度失配最低”，还要问这个角度下超胞是否太大、是否保留了你想研究的活性位或界面区域。
图2 氧化物-TMD 扭转异质结中莫尔图案形成和旋转角依赖的结构示意。
为什么超胞越大不一定越好？
从几何上说，超胞放大以后，总能找到更低的平均失配；从计算上说，原子数一旦涨到几百个，弛豫、静态、Bader、差分电荷和吸附构型枚举都会变得很重。更麻烦的是，太大的超胞往往让你只敢算一个构型，反而失去对不同堆垛和吸附位的比较。
先列出几组候选超胞，再比较失配率、应变分配、原子数、真空层、层间距和可保留的对称关系。一个能完成完整证据链的中等超胞，通常比一个算不动的大超胞更有价值。
式(2) [b’1 b’2] = R(θ)[b1 b2]
候选超胞应该先筛什么？
进入 VASP 前，建议先把候选方案做成一张表，而不是凭感觉选一个看起来顺眼的 POSCAR。至少要记录 A 层超胞、B 层超胞、θ、εA、εB、总原子数、层间初始距离和是否保留目标活性位。这样后面结果不理想时，才能回头判断是模型选择问题，还是计算设置问题。
图3展示了界面超晶格自动筛选时对候选结构、应变和弹性能的比较。这个思路很适合迁移到手动建模：先让几何方案排队，再让物理判断筛掉不合理项。超胞匹配不是为了找到唯一答案，而是为了留下几组代价可解释的候选模型。
图3 InterMatch 筛选二维/一维界面超晶格时的候选超胞、应变和弹性能比较。
弛豫后要重点看哪些异常？
第一看层内键长有没有被拉坏。若某一层的关键键长偏离原始材料太多，后续电荷转移、能带对齐和吸附能解释都要非常谨慎。第二看层间距是否合理，范德华异质结一般需要打开 DFT-D3、D4、optB86b-vdW 或相近的 vdW 描述，不能只靠普通 GGA 期待层间相互作用自己长出来。
第三看界面平移构型。即使晶格匹配相同，上层相对下层平移半个键长，活性位附近的配位环境也可能完全不同。图4中的扭转角和位移能量图说明，旋转只解决方向关系，平移还会继续改变局域能量谷。同一个失配率下，平移构型不同，界面结合能和局域电荷也可能不同。
最后再决定是否进入吸附计算。若空白异质结还没有稳定，直接放 *H、*OH 或 CO2 中间体，只会把结构误差放大。可靠的顺序是先确认裸界面，再做吸附位枚举，最后再做自由能和电荷分析。
图4 不同扭转角下界面位移映射和分离能分布，用于说明旋转与平移构型会共同改变界面能量。
结论怎样写？
写结果时，建议把“模型条件”放在结论前面。例如“在 3×3 MoS2 / 4×4 基底、εA 小于 2%、θ = 0° 的模型中”，比直接写“该异质结促进电荷转移”更稳。前者告诉别人你的判断边界，后者容易让人以为所有取向和所有堆垛都成立。
如果多个匹配方案给出相反趋势，也不要急着删掉不合意的模型。可以把差异写成界面构型敏感性：某些取向更利于电荷转移，某些堆垛更利于吸附中间体稳定。二维异质结的结论越依赖旋转角和堆垛，越要把模型边界写清楚。
最朴素的复查方法，是把 POSCAR 里的晶胞矢量、应变分配和层间构型截图保存下来。很多后期争议不是能量公式算错，而是最初那个“看起来差不多”的超胞，其实已经改变了你想讨论的界面。
表1 二维异质结匹配时必须记录的变量
|
变量 |
看什么 |
容易踩的坑 |
|
a1、a2 |
两层平面基矢长度和夹角 |
只比较一个晶格常数 |
|
θ |
两层相对旋转角 |
默认 0° 就是最合理 |
|
εA、εB |
应变分配到哪一层 |
把失配全部压给柔弱层 |
|
Natom |
总原子数和可算性 |
超胞过大导致只能算一个构型 |
|
dinter |
层间距和 vdW 描述 |
未开 vdW 就讨论界面结合 |
表2 从几何匹配到 VASP 计算的推荐顺序
|
步骤 |
要做的事 |
通过标准 |
|
1 |
列出 3-5 组候选超胞和旋转角 |
失配、原子数和目标位点都可解释 |
|
2 |
分配εA、εB 并生成 POSCAR |
关键键长没有异常拉伸 |
|
3 |
弛豫裸异质结 |
层间距稳定，界面没有塌陷或穿插 |
|
4 |
枚举吸附位或界面平移 |
能比较不同局域堆垛的影响 |
|
5 |
做差分电荷、结合能和吸附自由能 |
结构、电子和能量趋势相互托住 |
要点1：二维异质结匹配的不是单个晶格常数，而是两层平面基矢、夹角、旋转关系和超胞组合。
要点2：失配率依赖具体超胞方案；同一对材料在不同θ 和不同 m×n 组合下会得到不同结果。
要点3：应变分配会影响电子结构解释，不能把人为拉伸造成的变化直接写成界面效应。
要点4：VASP 建模应先确认裸界面稳定，再进入吸附位、差分电荷和自由能比较。
