2. 表面与界面结构建模基础

表面与界面第一性原理计算的第一步，不是点“开始计算”，而是把一块无限延伸的固体“切”成一个既能代表真实表面、又能在周期性边界条件下被正确求解的模型。这个模型就是 slab 模型。它决定了后续吸附能、功函数、界面结合强度等所有结果是否可信。本节围绕 slab 的构造逻辑，把真空层、偶极修正、对称性约束和周期性边界条件的选择原则串成一条可操作的建模思路。

一、slab 模型：把半无限表面变成可计算的周期板

真实表面在面内是无限延伸的，在垂直方向则终止于真空。第一性原理计算无法直接处理半无限体系，于是采用 slab 近似：沿表面法线方向取有限层原子，上下各留真空，面内仍保持周期性。这样得到的体系在面内周期重复，在法线方向被真空隔断，从而模拟出“表面”这一几何特征。

构造 slab 时，首先要确定三个关键量：层数、面内超胞大小和真空层厚度。层数决定表面效应能否被充分体现——通常需要让中间层原子接近体相环境，表面几层才允许弛豫。以常见的金属氧化物或半导体表面为例，教学示例中常取 5～9 个原子层，其中底部 1～2 层固定为体相位置，中间层弛豫，顶层完全自由。面内超胞则决定吸附质或界面失配的覆盖度，例如 2×2 或 √3×√3 超胞用于模拟不同覆盖度。

判断 slab 是否“够厚”的实用标准是：继续增加层数时，表面能、功函数或吸附能的变化小于设定阈值（例如 0.01 eV/Ų 或 0.05 eV）。这不是一次计算就能确定的，而需要做层数收敛测试。

二、真空层与偶极修正：解决周期性带来的虚假相互作用

真空层的作用是隔断相邻 slab 之间的电子云重叠。真空层太薄，上下两个周期镜像会互相“看见”，导致表面能偏低、功函数失真；真空层太厚，则计算量无谓增加。经验上，真空层厚度应使真空区域中电势趋于平坦，通常取 10～20 Å。对于带极性表面的 slab，真空层还需要配合偶极修正。

所谓极性表面，是指沿法线方向原子层由带电层交替堆叠构成，例如 ZnO(0001) 的 Zn 终止面和 O 终止面。这种 slab 上下表面不等价，会形成一个净偶极矩。在周期性边界条件下，这个偶极矩会与镜像产生虚假电场，使真空能级发生倾斜，功函数无法收敛。偶极修正的做法是在真空区域引入一个反向的补偿偶极层，抵消净偶极矩。常见实现方式是在 slab 两侧真空区加入一个锯齿形电势修正，或采用偶极层模型。判断是否需要偶极修正的简单标志是：不加修正时，真空能级在真空区出现明显斜率；加上修正后，真空能级恢复水平。

三、界面匹配：让两种材料在面内“对得上”

当计算对象从表面变成异质界面或吸附界面时，核心问题变成：两种材料的晶格常数不同，如何构造一个面内周期一致的超胞？这就是界面匹配。常用策略有三种。

第一种是最小公倍超胞法：寻找两种晶格的最小公倍数，使超胞在两个方向上同时接近整数倍。例如材料 A 的晶格常数为 3.0 Å，材料 B 为 4.2 Å，可取 7 倍 A 与 5 倍 B，得到约 21 Å 的匹配超胞。第二种是应变补偿法：固定一种材料，对另一种施加 1%～3% 的双轴应变，使其晶格常数匹配。第三种是近似匹配法：允许一定失配率（通常小于 5%），直接取相近的超胞。失配率定义为 |a_A − a_B| / a_A。教学示例中，若失配率超过 5%，界面应力会显著影响结合能，需要谨慎解释结果。

界面匹配完成后，还要确定界面间距和堆垛方式。常见做法是取若干不同间距和不同平移位置，计算界面结合能，找到能量最低的构型。这一步往往需要 5～10 个初始构型做筛选。

四、对称性约束与周期性边界条件的选择原则

周期性边界条件（PBC）是 slab 计算的前提：面内两个方向必须周期重复，法线方向通过真空隔断。选择 PBC 时要注意，k 点采样在面内方向需要足够密，法线方向只需 1 个 k 点（因为真空层隔断了色散）。例如面内 2×2 超胞，常用 4×4×1 或 6×6×1 的 Monkhorst-Pack 网格。

对称性约束则用于减少计算量并保持物理合理性。对于对称 slab（上下表面等价），可以启用镜面对称，使上下表面同时弛豫，节省一半计算量。但要注意：一旦吸附质只在一侧，对称性被破坏，就不能再强制镜面对称。对于极性 slab，通常不启用镜面对称，因为上下表面本身不等价。实际操作中，常见做法是：先判断 slab 是否中心对称，再决定是否开启对称性；若开启，需确认固定层和弛豫层的划分与对称性一致。

把以上几个量放在一起看，slab 建模其实是一个“用有限换无限”的近似过程：层数换体相深度，真空层换表面隔离，偶极修正换极性补偿，超胞换界面匹配，对称性换计算效率。每一个选择都有明确的物理判据，而不是凭经验拍脑袋。理解这些判据，才能在后续吸附能和界面能计算中，知道结果为什么可信、误差从哪里来。