Slab模型解析：从周期单元构造到表面能、吸附与电子态的计算边界

说明：本文华算科技主要介绍 slab模型在周期性第一性原理计算中代表的物理对象、上下表面与真空区的关系、有限厚度引入的数值影响，以及表面能、吸附和表面电子态的解释范围。

体相晶体可以沿三个晶格方向无限重复，表面却在法向上截断了平移对称。平面波 DFT 仍习惯处理周期单元，于是研究者取出有限厚度的晶体片层，在法向留出一段没有原子的空间，再让整个单元周期重复。这个由有限原子层与真空区共同组成的周期单元，就是表面计算中常说的 slab 模型。

slab：有限周期层替半无限表面

设表面法向为 z。模型在 x、y 方向保持晶格周期，因而描述一张横向延展的晶面；沿 z 方向只保留若干原子层。这种面内连续、法向有限的几何特征，把三维体相改写成二维周期表面。程序仍沿 z 方向复制单元，真空区分隔相邻片层的原子和电子密度。真空没有撤销周期性，计算对象依旧是“晶体片层—真空—晶体片层”的无限序列，只是相邻副本的耦合被压低。

图 1 的 LaFeO₃(001) 单元给出典型侧视结构：黑框圈定计算单元，原子层构成固体区，上方 11 Å 空白区构成真空。面内复制生成连续表面，法向复制后，上一个单元的真空紧接下一个单元的原子层。底部两层保持体相坐标，其余层松弛，用来近似表面附近逐层衰减的结构响应。

真实块体表面从外侧向内部延伸到宏观尺度，计算只能保留有限层数。slab 中部若已恢复接近体相的键长、电荷密度或层分辨态密度，中部原子层便可充当体相参照；外侧几层则承载配位降低、表面弛豫和电子重排。slab 是半无限表面的有限替身，模型厚度承担从外侧表面恢复至体相内部的数值任务；研究自由悬挂薄膜时，有限厚度本身才属于待研究的真实尺度。

一个slab同时规定了哪些物理对象？

暴露晶面的原子排布由切割体相所沿的密勒指数决定；切口位于不同原子层之间，又会产生不同终止。面内超胞规定缺陷、吸附物和重构单元的横向重复距离。对吸附体系而言，一个分子放进 1×1 与 4×4 表面单元对应不同覆盖度，分子间的横向相互作用也随之改变。晶面、终止与面内周期共同定义了所算的表面。

原子沿法向的位置还把“表面”和“内部”分开。图 2 中，H 原子位于 V(100) 的顶位、桥位、空位以及第一、第二次表层间隙时，局域配位各不相同。slab 为表面与次表层位置建立同一法向坐标，吸附能、层间迁移势垒和电荷转移才有明确的结构含义。模型底层固定时，它们模拟受体相约束的内部晶格；放松层数过少，则表面弛豫可能被固定区提前截断。

周期复制的有限片层同时包含上、下两个终止面，它们都会进入总能与电子态。上下终止相同、几何关于中面对称的模型称为对称 slab；两侧终止、吸附物或弛豫方式不同，则属于非对称 slab。对称结构可消去法向净偶极，并让两侧表面能相同。非对称结构适合保留单侧吸附或不同终止，但两个真空平台、两侧功函数和内建电场都可能不同。

图 3 的 AgCrSe₂ 对称超胞含有两个 Se1 终止面。上、下表面的反演对称破缺方向相反，表面态在能带图中分别以紫色和橙色投影标出。这里的“底面”仍是一个真实计算表面；对称构造借助第二个同类表面约束电势和电子态。

有限厚度/真空区为什么会改变计算量？

slab 的有限尺寸引出两条耦合路径。一条穿过真空：相邻周期副本之间仍有波函数尾部、库仑势、偶极场或吸附物之间的作用；另一条穿过固体：上、下表面的弛豫和电子态可在过薄片层中相遇。增大真空只削弱前一条路径，增加原子层才处理后一条路径，二者对应不同误差来源。

Al {111} 纳米片展示了厚度效应的直观形式。灰色表面层在三层片中占据大部分原子，随层数增加，绿色内部层比例上升，表面原子占比随厚度增加而下降，面内模量逐渐接近体相值。对 DFT 表面模型，同类有限厚度效应会出现在层间距、磁矩、带隙、表面态分裂与表面能中；某一总能差已经平稳，并不自动保证层分辨电子结构也已平稳。

真空厚度应由目标物理量检验。图 5 将 Mg 表面含水体系的法向平均静电势、电荷密度和总能变化并列。较短单元中，水分子附近的密度尾部与势能尚受相邻副本影响；z 轴长度增加后，总能变化趋近于零。这里得到的 25.4 Å 属于该结构和检验量，换成带电表面、强偶极吸附层或更长分子，所需间距会改变。

非对称 slab 的法向电荷分布带有净偶极，周期复制后会在真空区形成额外电场，平面平均势可能呈斜坡。偶极修正针对这项周期静电伪影，它不替代真空间距、片层厚度或终止方式的检验。极性表面还可能通过电荷补偿、重构或化学计量调整降低极化场，单纯增加真空无法消除固体内部持续存在的电场。

从slab结果能读出哪些表面物理？

对含两个等价表面的化学计量 slab，常见表面能写作 γ=(E_slab−NE_bulk)/(2A)，分母中的 2 来自上、下两个等价表面。若两面不等价，这个式子给出两侧贡献的平均值，一次总能无法分别求出各侧 γ。吸附能同样依赖清洁 slab 参照、面内面积与覆盖度；比较不同吸附构型时，这些基准应保持一致。

功函数取真空平台与费米能级之差 Φ=V_vac−E_F。非对称 slab 两侧的 V_vac 可不同，报告 Φ 时应注明对应表面。表面态可通过层投影能带或层分辨态密度识别：表面态权重集中在外侧若干层并向内部衰减，且能与体相投影区分。此时 slab 厚度还应让上下表面态的杂化分裂降到研究所容许的尺度。

图 6 中，六层 CsV₃Sb₅ 对称 slab 的外侧层弛豫后，表面局域势比内部低约 1 eV。把能带投影到最外一层或最外两层，可分辨表面量子阱态，并与 ARPES 谱和费米面比较。这里得到的是特定终止、层数与弛豫结构下的表面谱；若这些模型条件改变，局域势和表面带也可能随之移动。

把 slab 结果用于真实表面时，结论应连同模型条件一起陈述：晶面与终止对应哪种化学环境，面内单元对应多大覆盖度，上下表面是否等价，内部是否恢复体相特征，真空平台是否稳定。表面重构、台阶、溶剂、电极电势或有限温度没有写入模型时，计算量描述的是理想化周期假设下的表面，外推到实验体系还要核对缺失的结构与环境因素。