说明：本文华算科技主要介绍第一性原理表面计算中为什么要按不同晶面建立 slab 模型，以及 Miller 指数、表面终止、配位环境、表面能和吸附能怎样共同影响计算判断。
Slab 建模需明确哪些参数？
晶体内部的原子在三维周期中重复排列，表面计算把这个周期体相沿某个 Miller 指数 (hkl) 方向截开，再在法向加入真空层，形成面内周期、面外隔开的 slab。切晶面定义的是暴露出来的原子排布，不是单纯给模型换一个名字；(100)、(110)、(111) 对应不同原子面、层间距和法向堆垛。
图1. 从实验体相结构生成表面取向、判断极性并进入形貌预测的计算流程。
Miller 指数决定晶面法向，表面终止决定最外层露出哪类原子。对于金属 fcc 结构，(111) 常形成密排平台，(100) 暴露方形网格，(110) 更开放；对于氧化物和层状材料，同一个 (hkl) 还可能出现阳离子终止、阴离子终止、混合终止或重构终止。表面终止和晶面指数必须成对记录，因为后续总能、偶极、吸附构型都绑定在这一个原子模型上。
slab 模型里还存在厚度、真空层、固定层和弛豫区域。厚度太薄时，两个表面之间会互相影响；真空层不足时，周期镜像会改变静电势和功函数；底层固定过多时，表层弛豫又会被限制。同一晶面之间的比较要共享这些模型条件，否则能量差里会混入设置带来的偏移。
实际建模时，体相结构、切面枚举、极性判断和形貌预测常常连在同一套计算任务中。一个表面取向可能对应多个终止，某些极性 slab 会积累面外偶极，某些重构又会降低表面能。多晶面计算的第一步是把候选表面模型登记完整，包括 (hkl)、终止、层数、真空层、两侧是否对称和表面是否重构。
对于同一个 (hkl)，如果上下两侧表面不等价，平面平均静电势会出现斜率，功函数和带边对齐还受面外偶极影响。对称 slab、非对称 slab 和偶极修正对应不同物理假设，它们不能混用在同一组表面能或吸附能排序中。
把图1中的每一步落实到输入模型时，切面选择不再只是一个 Miller 指数。模型表要写清 slab 两侧是否等价、终止层是否电荷中性、弛豫后是否发生重构；这些项目决定表面能、功函数和吸附能能不能放在同一张表里比较。
晶面差异如何影响表面性质？
切面一变，最外层原子的配位数、键长弛豫和局域电荷分布会随之改变。anatase TiO2 的 (001)、(010) 与 (101) 表面中，Ti5f、Ti4fR、O2f、O3f 等低配位原子的空间位置并不相同；这些低配位位点控制吸附物最先接触的局域环境，并改变表面附近的电子态分布。
表面原子的配位数减少后，未饱和轨道、表面偶极和局域态密度都会发生重排。同一个体相带隙或同一个 d 态中心，放到不同晶面上会对应不同表面态，吸附物与表面形成的键长、键角和电荷密度差分区域也随之改变。
图2. anatase TiO2 不同晶面的优化表面结构、平面平均静电势和局域态密度。
同一材料的不同晶面还会给出不同真空能级平台、费米能级参照和 PDOS 峰位。对功函数、带边对齐、表面偶极或吸附前后电子重排而言，晶面本身已经参与定义能级零点；把某个晶面的 EF、Evac 或 PDOS 直接套到另一个晶面，会把表面几何差异误写成材料本征差异。
在吸附计算里，同一吸附物面对不同表面原子时会形成不同初态和弛豫终态。CO2 放在 TiO2(001) 表面时，同一个晶面内仍有 O3f、Ti5f、hollow 和 O2f 等候选位置，结合能随构型改变；切出晶面后，应在该晶面上寻找非等价位点和分子取向。晶面筛选和构型筛选是两个相邻步骤，前者定义表面，后者定义吸附态和能量参考。
图3. CO2 在 anatase TiO2(001) 表面的多个吸附构型及对应结合能。
当 CO2 转到 TiO2(101) 表面，最稳定构型、分子弯曲程度和结合能范围都会变化。这里比较的不是孤立的 CO2 分子，也不是抽象的 TiO2 体相，而是“CO2 + 具体晶面 + 具体吸附位点”的整体模型。若目标性质依赖表面反应或表面能级，单一晶面模型只覆盖了一个局域表面状态，后续自由能图和电荷分析也只对应这个表面。
图4. CO2 在 anatase TiO2(101) 表面的多个吸附构型及对应结合能。
什么时候必须同时切多个晶面？
真实颗粒、薄膜和纳米晶通常会暴露多个表面，实验样品里写作同一种材料时，表面比例可能由晶体生长、退火气氛、吸附物覆盖和电化学电位共同改变。当研究对象是表面能、吸附能、功函数、界面反应或形貌演化，多个低指数晶面的并列计算比单个 slab 更接近真实样品的表面集合。
单晶实验、外延薄膜和择优取向纳米片可以把晶面范围收窄，多晶粉体和反应后重构表面则要保留更多候选面。实验表征给出的暴露面比例可以决定计算优先级：主暴露面用于基准能量，少量台阶面或高指数面用于检验低配位位点贡献。
材料存在极性表面时，多晶面筛选还承担模型排错功能。Li2S、TiO2、LiCoO2 和 LiFePO4 的不同终止会产生不同偶极矩，某些 slab 的表面能随厚度持续上升。厚度依赖和偶极矩是识别异常极性终止的关键量，它们能提示某个切面是否要采用对称 slab、表面补偿或重构模型。
图5. Li2S、TiO2、LiCoO2 和 LiFePO4 中不同 slab 终止的表面能、偶极矩和厚度依赖。
表面能排序和 Wulff 构型常用来估计哪些晶面更可能大面积暴露。低表面能晶面更容易进入平衡形貌，高表面能晶面可能只在台阶、边缘、缺陷或动力学生长条件下出现；但离子晶体中，断裂哪一类局域结构单元同样改变表面代价和局域应力分布，例如 PO4 四面体是否被切断，会影响 LiFePO4 不同表面的能量排序。
图6. 多种离子晶体的低表面能晶面、代表性终止结构和重构效应。
催化和电化学场景里的表面模型还要覆盖反应通道差异，平整低指数面适合建立基准吸附能，高指数面、台阶面和缺陷面更容易暴露低配位原子；吸附物覆盖后，表面能和可用位点又会更新。若实验活性来自多晶颗粒，计算中只保留最平整的低指数面会缩小可讨论的表面状态，尤其会漏掉台阶位、棱边位和重构位点上的反应能量。
多晶面计算该如何规范对比？
多晶面计算要把比较基准固定在同一组数值设置下。泛函、赝势、U 值、vdW 修正、k 点、截断能、slab 厚度、真空层、偶极修正和表面弛豫标准应保持一致；吸附计算还应固定覆盖度、吸附物参考态和自旋设置。只有模型条件一致，γ、Eads、ΔG 或 Φ 的差值才主要反映晶面差异，而不是数值设置造成的能量漂移。
比较顺序可以按“表面存在概率”和“目标反应能力”分开处理。表面能或表面自由能用于判断某晶面在给定环境下是否容易暴露；吸附能、反应自由能、NEB 能垒、PDOS、差分电荷密度和 Bader 电荷用于判断该晶面上的具体反应步骤。一个低表面能晶面可能稳定但反应位点少，一个高指数面可能面积小却提供低配位活性位。
高通量数据库把这种分工写成样本变量：材料组成、表面取向、吸附物、初始构型和弛豫任务分别登记，后续机器学习模型才有机会区分“成分变了”“晶面变了”还是“吸附构型变了”。Open Catalyst 2020 这类数据集覆盖大量 surfaces 与 adsorbates，其核心样本并非体相材料条目，而是表面-吸附物结构，每一个样本都对应明确的 slab 和弛豫路径。
图7. Open Catalyst 2020 数据集中不同材料表面、吸附物和弛豫任务的整体示意。
交付结果时，建议把每个晶面的角色写成可复查表：哪些晶面来自 Wulff 或实验暴露，哪些晶面只是低指数基准，哪些终止因极性或重构被剔除，哪些吸附构型进入反应路径。晶面角色表能限制后续 Eads、DOS、Δρ、ΔG 和 NEB 能垒的可比范围，也能让实验晶面比例、TEM 暴露面和电化学重构结果与计算模型建立一一对应。
若文章只讨论“为什么要切不同晶面”，结尾应停在模型对象和判断量上：切面决定表面原子，终止决定表层化学组成，弛豫决定局域几何，吸附态决定反应能量，这些对象分别对应表面能、吸附能、功函数、PDOS 和反应自由能。