VASP 默认三维周期性边界条件 PBC(这是根本原因)
VASP 对晶胞三个方向全部自动周期性复制,你建的薄层 / 表面会在 z 轴无限堆叠,上下两个表面直接互相接触、会发生虚假层间相互作用,导致
范德华力、静电偶极、电子波函数会跨层重叠;
优化后两层表面吸附分子互相吸引,键长、吸附能、带隙全部失真;
而真空层就是在 z 方向隔开重复的表面,切断上下周期镜像的相互作用。
区分 “块体” 与 “表面” 体系
块体材料三维都是实原子,不需要真空;
表面模型是二维薄膜 + 真空,只保留 xy 面周期性,z 方向人为用真空隔绝,模拟真实孤立表面。
不加真空 = 等效多层堆叠块体,根本不是表面。
消除 z 方向周期性带来的虚假静电相互作用
表面存在悬挂键、电荷极化、偶极矩,无真空层时,上下镜像表面偶极同向叠加,电场无限延伸,能量剧烈震荡、电子难收敛。搭配真空层后才能开启LDIPOL偶极修正,抵消虚假偶极误差。
给吸附分子预留运动空间
做吸附、催化计算时,分子会在表面上下弛豫、翻转、振动;
真空厚度不足会出现:
分子撞到上层镜像表面,原子飞散、结构崩坏;
力收敛困难,EDIFFG 始终达不到标准。
保证波函数、电荷密度充分衰减
电子波函数在原子外会缓慢拖尾,真空层厚度足够才能让波函数在真空区域衰减至接近 0,避免相邻周期薄膜的电子云交叠:
交叠会导致基组重叠误差 BSSE 大幅升高;
能带、差分电荷、功函数计算完全不准。
正确计算功函数时,需真空能级做参考
功函数 = 真空能级 − 费米能级
真空区域是平直的真空能级平台,只有真空足够厚,才能在真空层中间取稳定的静电势做参考;
真空太薄,电势曲线起伏,无法读取准确真空能级,功函数完全错误。
真空层厚度选择经验
单纯金属 / 非金属洁净表面:≥15 Å
含吸附分子、多层异质结:≥20 Å
极性表面、偶极明显体系:建议 20–25 Å
判断标准:真空中间区域静电势平直无波动即可
新手常见误区
真空层只留 5–8 Å:波函数交叠,能量严重偏差;
二维材料(MoS₂、石墨烯)不加真空:层间范德华作用被周期重复放大,层间距弛豫错误;
表面计算误用 ISIF=3:优化时 z 轴晶格会压缩真空层,一定要用 ISIF=2 固定晶格,只弛豫原子。
207