做VASP表面计算时，最容易把时间浪费在后期反复返工上，而根子往往不在INCAR参数，而在模型从一开始就没有搭稳。VASP官方教程和Wiki都强调，表面计算本质上是从体相切出slab，再沿表面法向拉出真空区，同时slab厚度和真空厚度都需要分别做收敛检查，不能只盯其中一个。
一、VASP表面模型怎么构建
表面模型构建这一步，真正关键的不是先把晶胞拉长，而是先把体相基准、切面方向和slab厚度按顺序理顺。官方表面教程给出的主线很明确，就是先从体相出发构造表面，再做弛豫和表面相关性质分析，所以前面的模型骨架必须先站稳。
1、先把体相结构收敛好
表面slab不是脱离体相单独存在的，切面以前最好先把体相晶格常数、磁性状态和基础能量收敛清楚。这样后面切出来的表面模型才有稳定参照，不会把体相误差直接带进表面能和吸附能里。
2、再按目标晶面切出slab
切表面时要同时定晶面方向和表面终止，而不是只写一个面指数就结束。VASP Wiki明确指出，干净表面由体相晶体、切割晶面和表面终止共同确定，所以上下表面到底是哪种终止，一开始就要看清。
3、slab厚度先从中等层数起步
模型不要一开始就切得太薄，也不要第一次就堆得很厚。更稳的做法，是先用中等层数建立测试模型，再看表面能、功函数或吸附能是否随层数继续明显变化，因为官方已经明确提示slab thickness需要单独收敛。
4、弛豫时保留部分底层作为体相支撑
VASP官方表面教程示例采用了五层镍表面、只放松其中上部层的做法，这种思路在实际表面建模里非常常见。这样做的目的，是让表面能弛豫，底部又保留接近体相的支撑骨架，后面结果通常更稳。
二、VASP表面真空层厚度怎么确定
真空层厚度不是越大越好，而是要大到足以让相邻周期镜像之间不再明显相互作用。VASP官方在功函数和表面条目里都强调，真空区必须形成足够干净的无电荷区域和近似平坦的真空势平台，否则表面相关性质就很难可靠。
1、先用可工作的初值起算
如果只是做第一轮测试，官方给出的经验范围是法向两侧各留大约8到12埃真空区。这个数不是最终答案，但足够让你先观察电势是否开始出现平台，也适合作为后续收敛测试的起点。
2、再用目标物性做真空层收敛
真正确定真空层厚度时，不要只看几何尺寸，而要看你关心的量是否稳定，例如表面能、功函数、吸附能或真空势平台。VASP Wiki已经明确指出，表面能需要同时对slab厚度和vacuum thickness做收敛检查。
3、非对称表面要同步考虑偶极修正
如果slab上下表面不对称，或者一侧带吸附分子，即使真空区看着够厚，也可能因为净偶极导致真空中仍然残留电场。VASP Wiki明确建议这类slab计算使用LDIPOL，并把IDIPOL设在表面法向方向，必要时再配合DIPOL。
4、真空区没有平台就继续加厚
如果你从LOCPOT或真空势结果里还看不到足够平坦的区域，不要急着解释功函数或表面势垒，先继续增大法向晶胞长度。官方关于功函数和静电修正的说明都把“先确保有足够真空”放在前面，这一步是前提，不是可选项。
三、VASP表面收敛还要检查什么
很多表面模型看起来已经能跑，但结果还是不稳，常见原因并不是参数太少，而是模型方向、偶极设置和分析口径没有一起核对。VASP官方关于表面、功函数和偶极修正的说明实际上已经把这些检查点给出来了，所以正式出结果前最好再把这几项连起来看一遍。
1、先确认表面法向和第三晶格矢量一致
表面计算里，真空区通常沿某一个晶格方向展开，而很多后续设置都默认你已经把表面法向定义清楚。若法向理解错了，真空层、偶极修正和电势分析都可能一起偏掉。
2、再确认模型是对称slab还是非对称slab
对称slab和非对称slab的后处理思路并不一样，尤其在功函数和偶极修正上差别很明显。非对称模型若不额外检查，往往更容易在真空区残留净场。
3、再看LDIPOL和IDIPOL是否与法向一致
如果已经决定做偶极修正，那就不要只开LDIPOL而不核对方向。VASP Wiki对IDIPOL的建议很直接，slab计算应把它设在表面法向对应方向上，否则修正本身就可能偏离目标。
4、最后同时复核层数和真空而不是只改单一参数
很多人会先把真空层加厚很多，却不再回头检查slab层数，或者反过来只加层不加真空。官方已经明确提醒，slab thickness和vacuum thickness都需要独立收敛，所以这两项最好一起做交叉检查，不要只改单边。
总结
VASP表面模型怎么构建，VASP表面真空层厚度怎么确定，真正稳的做法不是先把晶胞随手拉长，而是先从体相出发切出合理slab，再分别对层数和真空厚度做收敛检查。等模型骨架、真空平台和偶极方向都确认以后，再去看表面能、功函数和吸附结果，后面的解释才更容易站得住。
