slab真空层要设多厚？15 Å 的适用边界与收敛测试方法

引言

真空层设15 Å。这个经验值确实常见，但它不是万能答案。表面有强偶极、吸附物伸得很高、体系带电、界面含水层，或者要算功函数和能级对齐时，15 Å 可能不够；而对很简单的非极性表面，盲目加到 30 Å 又会浪费计算量。真空层厚度不是固定参数，而是需要随体系和目标性质做收敛测试的边界条件。本文就讲 slab 真空层到底该怎么判断。

15Å为什么只是起点？

真空层到底在隔开什么？

slab 模型沿 z 方向周期重复。真空层的作用，是把一个表面和它在下一周期里的镜像隔开，避免电子密度、静电势和吸附物尾巴互相“看见”。如果真空太薄，上下两个周期会产生假相互作用，表面能、吸附能、功函数和能级位置都可能被带偏。

图1的静电势曲线展示了为什么要看真空平台。对功函数或能级对齐来说，真空层不是只看几何距离，而是要看平面平均静电势在真空区有没有形成稳定平台。如果真空区没有平坦电势平台，功函数和真空能级就没有可靠读数。

 图1  MAPI(001) slab 的平面平均和宏观平均静电势、真空能级、VBM/CBM 能级对齐示意。DOI: 10.1021/acsomega.0c04420。

为什么同样 15 Å，有的体系够，有的不够？

简单、非极性、吸附物不高的金属表面，15 Å 往往可以作为不错的起点。但如果表面是极性终止、有强偶极、有外加电场、有带电吸附物，或者吸附分子本身伸出表面很远，镜像之间的长程静电作用就会增强。此时 15 Å 只是初猜，不是结论。

尤其当你关心Φ、E_vac、VBM/CBM 对齐或带电缺陷时，真空层厚度会直接影响参考能级。真空层是否足够，取决于目标性质对长程静电势的敏感程度，不能只按别人的默认参数照抄。

式(1)  L_vac = c – h_slab+ads

怎样判断真空层已经收敛？

只看总能差够不够？

最直接的办法是做 L_vac 收敛：比如 12、15、18、22 Å 依次测试。对吸附能来说，可以看 E_ads 的变化是否低于自己设定的容差；对功函数和能级对齐来说，还要看真空平台位置和 E_vac 是否稳定。

图2就是典型的真空厚度收敛例子。它说明真空层不是装饰性空白，而会影响相对于真空能级的 slab 态能级。判断真空层收敛时，至少要同时看能量差和静电势平台。

图2  MAPI(001) slab 态能级随真空层厚度变化的收敛测试。DOI: 10.1021/acsomega.0c04420。

偶极修正什么时候要一起开？

如果 slab 上下两面不对称，比如只在一面吸附分子、只改性一侧表面、或者界面本身有净偶极，周期重复会让相邻镜像之间产生人工电场。增大真空层可以减弱这个问题，但不一定彻底解决；这时通常需要考虑 LDIPOL、IDIPOL = 3 和 DIPOL 设置。

不要把偶极修正当成“开了就万事大吉”。更稳的做法是比较开关偶极修正后的功函数、吸附能和静电势曲线；如果差异很大，就说明体系对长程偶极敏感。对不对称 slab，真空层收敛和偶极修正应当一起检查。

式(2)  |E_ads(L_vac,2) – E_ads(L_vac,1)|

哪些体系最容易被15Å误导？

含水层或大分子吸附为什么更麻烦？

图3中的水合 RuO₂ 模型显式包含固体表面、水层和 15.0 Å 真空层。这里的 15 Å 不是随手填的空白，而是给液水界面、表面偶极和周期边界共同留下空间。若水层更厚、分子取向更强、界面电场更明显，就要重新测试。

大分子吸附也类似。吸附物越高，离下一周期表面越近；如果还有氢键网络、离子或外电场，真空层不足会把真实界面作用和镜像作用混在一起。真空层要按最高吸附物、界面偶极和目标性质一起判断，不是只量裸 slab 顶端到盒子边界的距离。

 图3  水合 RuO₂(110) DFT-MD 模型盒示意，包含固体、水层和 15.0 Å 真空层。DOI: 10.1002/chem.202102356。

二维材料、带电吸附物和外电场怎么处理？

图4展示单层材料与 OH^– 吸附相关模拟盒。带电或强极性吸附物会让长程库仑相互作用更明显，真空层、偶极修正、背景电荷处理和溶剂模型都会影响结果。此时只说“我用了 15 Å”并不够，最好报告收敛测试或至少说明边界处理。

实操时可以把 15 Å 当作第一版，随后按目标性质加测试：如果只比较几何稳定性，容差可以宽一些；如果要算功函数、电场、能级对齐或带电吸附能，建议更严格。越依赖真空参考和长程电势的结果，越不能省略真空层收敛。

目录整理也要跟着这个思路走。建议把 12A、15A、18A、22A 分成平行目录，每个目录只改 c 方向盒长，其余 POSCAR 内部坐标、INCAR、KPOINTS 和 POTCAR 保持一致；静态计算结束后统一提取 OUTCAR 能量、LOCPOT 平面平均和目标性质。这样看到趋势漂移时，能明确归因到 L_vac，而不是把弛豫结构、k 点或能量截断的变化混进来。

写方法部分时也别只留一句默认参数。更稳的表述是说明“真空层经 L_vac 收敛测试后选定为某个数值”，并给出 E_ads、Φ 或目标能级变化已经低于设定阈值；这样别人复算时知道你控制的是边界误差，而不是凭经验填了一个数字。

 图4  BN 单层吸附阴离子 OH^–时使用的不同模拟盒与吸附态设置。DOI: 10.1021/acs.jpclett.6b02248。

表1  不同体系的真空层判断重点

本章要点总结

要点1：15 Å 只是 slab 真空层的常用起点，不是所有体系的统一答案。

要点2：真空层是否足够，要看周期镜像相互作用、静电势平台和目标性质是否收敛。

要点3：功函数、能级对齐、极性 slab、含水界面、带电吸附物和外电场对真空层更敏感。

要点4：不对称 slab 应把真空层收敛和偶极修正一起检查，不能只靠加大真空层解决。