说明：本文华算科技主要介绍表面计算中真空层的模型来源、法向周期镜像、电势平台和收敛判断。

表面计算把体相晶体沿某个 Miller 指数截断，只保留有限厚度的 slab，再让模型在面内方向重复。平面波 DFT 的求解器仍按三维周期晶胞工作，真空层承担的是面外数值隔离，让相邻 slab 镜像在法向上保持足够距离。

真实表面外侧可能是真空、气体、溶剂或电解质；计算盒里的真空层只是在三维周期模型中给表面留出空区。表面能、吸附能、功函数、带边对齐和平面平均静电势对这段空区的敏感程度不同，真空层厚度应随目标读数核对。

表面slab模型的基本概念

表面 slab 由两类方向组成：a、b 方向保留晶体在表面内的周期性，c 方向把有限层数的原子片和空区放进同一个晶胞。面内周期是材料模型的一部分，面外重复则是计算方法带来的镜像排列，真空层正是为这组镜像留出的距离。

图1. In2O3 的 (111)、(110) 和 (100)-O 表面 slab 模型。

图1里的 In2O3 slab 对应有限原子层对半无限表面的近似。若 c 方向没有空区，相邻 slab 会首尾相接，模型会退回三维连续固体，表面态、吸附位点和外侧势场都失去原本的含义。

表面超胞常在 slab 的基础上继续放大面内尺寸，用来放置空位、吸附物或低覆盖度构型。UO2(111) 表面缺陷模型里，面内超胞控制缺陷之间的横向距离，真空层控制上下周期镜像之间的法向距离，两个距离分别影响不同误差来源。

图2. UO2(111) 表面 slab、2 × 3 × 1 超胞和氧空位位置。

真空层的数值大小通常从 10 Å、15 Å 或 20 Å 起步，但经验数值只给初始模型。金属表面电子云衰减较慢，极性吸附层会产生面外偶极，二维材料的介电响应还会随晶胞高度换算；同一个 20 Å 在不同读数里含义并不相同，应以能量、力和电势曲线的变化幅度来定。

二维 slab 计算把同一张原子层放进不同 c 轴长度时，晶胞高度会进入按体积归一化的电荷密度、极化、介电函数和光学响应。若直接使用三维体积，真空越厚，平均数值可能越小；材料本身没有改变，改变的是归一化体积。比较二维体系的电偶极矩、面电荷密度或光学吸收时，应换算到面积、层厚或统一的有效体积，再报告真空层厚度。

真空层怎样隔开法向周期镜像

三维周期边界会把 slab 沿 c 方向复制成无限串列，真空层过薄时，电子密度尾部、范德华作用和长程库仑势会跨过空区相互作用。吸附能、表面能和表面态能级都可能随镜像距离改变，尤其在单侧吸附、带电缺陷和极性表面中更敏感。

TiO2(110) 水合表面模型给出的是典型 slab 几何：底部原子层模拟体相约束，顶部暴露配位不饱和 Ti、桥氧和羟基位点。真空层把上表面和下一周期的底面隔开，吸附物才主要响应当前表面，而非响应上一镜像的原子层。

图3. 水合金红石 TiO2(110) 表面 slab 中的 Ti、O 与 H 位点。

表面外侧若存在显式水层或粗糙形貌，计算盒的法向尺寸还要容纳水密度衰减区和界面起伏。粗糙 Cu–水界面模型把金属表面、水相和密度剖面放进同一盒子，法向空间不足会改变水层排布和金属表面松弛，后续 DFT 子区截取也会继承这种几何条件。

图4. 粗糙 Cu–水界面的表面结构、模拟盒和水分布剖面。

真空层收敛测试的对象应直接对应目标量：吸附问题关注 Eads 随真空厚度的变化，表面缺陷关注 VO 形成能与缺陷态位置，二维材料关注带边和面外极化。当目标量变化小于预设能量容差，这段空区才算服务于当前科学问题。

带电 slab、极性终止和强吸附偶极在真空区产生缓慢衰减的电势，简单增加空区会降低镜像耦合强度，却可能留下斜率或常数项偏移。法向镜像距离与电势参考的处理应分开记录，尤其在带电缺陷、外加电场和非对称吸附构型中要写明补偿电荷或偶极修正设置。

表面吸附计算还受覆盖度和横向超胞控制。真空层只处理上下镜像距离，面内吸附物之间仍按超胞周期重复；横向镜像耦合和法向镜像耦合要分开收敛。H*、*OH、CO2 或有机链吸附在同一表面时，覆盖度变化会改写 Eads，厚真空不会消除面内吸附物相互作用。

真空区电势平台为何影响能级读取

功函数和绝对带边位置把晶体内部能级同真空能级相连，常用关系写作 Φ = Evac – EF。这里的 Evac 来自真空区平面平均静电势的平台，平台是否平坦直接限定能级参照，与总能是否已收敛属于两个问题。

图5. MoO3 和 MoO2.875 单层的真空界面平面平均静电势。

图5的 MoO3/vacuum 曲线在左右两侧标出 vacuum 区域，费米能级用 0 eV 作为参照。若真空层太薄，曲线两侧平台会相互干扰；若 slab 上下表面不同，左右平台可能给出不同 Evac。功函数读数必须说明取哪一侧平台，半导体带边对齐同样受这一步影响。金属表面、半导体 slab 和二维层状材料的真空平台读取口径应保持一致。

界面模型里，平面平均电势还承担势垒分解任务。CuS/MoS2 接触计算中，ΔV 来自界面势阶，Hads-MoS2 的表面不对称还引入 ΔVAsym；真空层给电势曲线提供可读的外侧参照，电势阶才有可比较的零点。

图6. CuS/Hint-MoS2 接触和平面平均静电势中的 ΔV 与 ΔVAsym。

非对称 slab 中，单侧吸附、Janus 结构或一侧钝化会让电荷重心偏向某一侧。真空区平面平均势出现斜率时，偶极修正处理的是周期排列产生的面外电场；真空层、slab 层数、表面终止和固定层仍按各自目标量核对。

电势平台的核对不只服务于功函数。半导体表面要把 VBM、CBM 与真空能级相连，异质结要比较接触前后的外侧参照，电化学近似中还常把电子化学势换算到 SHE 或 RHE 标尺。Evac 的读取位置若没有写明，带边位置、肖特基势垒和电极电势换算会混用不同零点。

真空层收敛怎样和slab厚度分开判断

真空层厚度改变 slab 与法向镜像的距离，slab 厚度改变有限原子片对半无限体相的近似程度。这两项收敛服务于不同物理误差：前者对应外侧空区，后者对应内部层间弛豫、表面态穿透深度和固定层约束。

图7. Al(111) 表面能随原子层数变化的收敛结果。

图7里的 Al(111) 表面能随原子层数变化，在 3 层和 4 层之间差异明显，后续层数趋于平缓。这个趋势来自 slab 厚度，与真空层无关；厚真空无法补偿过薄 slab 的内部弛豫误差，反过来，增加原子层也无法消除上下镜像靠得过近的问题。

一篇表面计算文章记录模型时，真空层可写成目标量收敛曲线，避免只列孤立参数。吸附能可列出 15、18、20 Å 的 Eads 差值，功函数可列出 Evac 平台和 EF 的差值，表面缺陷可列出缺陷态位置随 c 方向长度的变化。同一模型内部比较应保持真空层、覆盖度、k 点和固定层口径一致，跨论文比较则要核对这些条件。

真空层对电子结构图谱的影响也要按目标读数分开处理。DOS 或 PDOS 若只比较局域轨道杂化，主要关注能量零点和投影原子；绝对带边、功函数、界面势降和外加电场响应则依赖真空参照。平面平均势、总能差和结构弛豫分别给出不同证据，混用同一个收敛数字会掩盖误差来源。

真空层收敛表最好只列服务于正文读数的项目。若文章讨论吸附构型，表中可给出 Eads、关键键长和表面弛豫量；若文章讨论功函数，表中可给出 Evac 平台差、Φ 和平面平均势斜率。收敛指标跟随科学问题选择，参数表才不会变成没有物理指向的数值堆叠。极性或单侧吸附 slab 可另列有无偶极修正的 Φ、Eads 或 ΔV 差值，把法向电场处理和真空厚度收敛分开登记。

表面计算加真空层，实质上是在三维周期算法中给二维表面模型设置面外间隔。足够的真空层让 slab 镜像、外侧电势平台和目标物理量各自达到稳定读数；slab 层数、表面终止、偶极修正和吸附覆盖度仍按各自模型条件登记。