表面切面怎么切：晶面指数定方向、终止方式定原子、极性/台阶表面怎么建模型

说明：本文华算科技主要介绍晶面指数如何确定切割方向、同一取向下切割位置与表面终止的差异，以及极性和实际台阶表面对计算模型的影响。

晶面指数怎样定义切割方向？

晶体表面计算从周期性体相结构出发。晶胞的三个晶格矢量给出周期，原子基元给出每个周期中的具体位置。选定晶面后，沿它的法向截取有限厚度的晶体，面内保留周期性，法向两侧留真空，才得到用于电子结构计算的表面模型。切面是穿过三维原子排列的一组平行平面；同一个晶面指数只固定这组平面的取向，尚未指定刀口穿过哪一层原子。

Miller 指数 (hkl) 来自晶面与三个晶轴的截距：把截距换成晶格常数的倍数，取倒数，再约成最小整数。某轴平行于晶面时，该轴截距无穷远，对应指数为 0。负指数常写成数字上方的横线。以立方晶体为例，(001) 平行于前两个晶轴，法向沿 [001]；(110) 斜穿前两个晶轴，形成另一种面内原子网格。(hkl) 写晶面，[uvw] 写方向；非立方晶系的表面法向需要从倒易晶格求出。同一晶面族写作 {hkl}，只有晶体对称操作能互相映射的成员才具有相同的理想晶格环境。六方晶体常用 (hkil) 四指数记法，第三个面内指数满足 i = −h − k；ZnO(0001) 中的四个数正是这种记法。

更一般地，法向由倒易矢量 G_hkl = h b₁ + k b₂ + l b₃ 确定，其中 b_i 是倒易晶格基矢。面内则需找出两条处于该平面的周期矢量；它们决定表面晶胞面积、原子重复方式，也决定异质结两侧能否形成低应变匹配。对于 β-Ga₂O₃，解理面穿过的键数和键种影响裸露表面的初始结构，需用晶体原子坐标逐层确定。单斜晶体的晶轴倾斜，倒易矢量给出真实的表面法向。

SrTiO₃ 的不同取向把同一块体结构投影成不同面内网格：(001) 和 (110) 的排列与 (111) 的三角网格明显不同。图中的 WS₂ 与 SrTiO₃(111) 叠加后出现莫尔周期，扭转角变化还会改变周期尺寸。选晶面同时选中了界面的二维几何条件；这里的莫尔图案对应 WS₂ 与 SrTiO₃ 的晶格参数，换用其他氧化物时须重新比较面内周期。

同一晶面沿法向移动，会切出何差别？

选定法向后，平面还可以沿法向平移。切割位置 δ 决定哪些键被截断、哪些原子留在最外层。硅的同一 (111) 取向内，候选切割平面彼此平行；在一个层间周期内移动切口，局部断键构型仍可能变化。下面的 Si 图将晶向分布、对应形貌与 [111] 方向的切割位置并列；右侧阶梯状能量曲线把 δ 的影响单独显出。图中能量由通用原子间势估计，描述初始切割的几何与断键差别；经过 DFT 弛豫的平衡表面能另需计算。

沿 SrTiO₃(001) 法向，SrO 层与 TiO₂ 层交替堆叠；把切口移过半个层序，就可能从 SrO 终止转成 TiO₂ 终止。图中的表面并非只有一种理想终止：SrO 可以局部覆盖在另一种表面区域上，不同发射角的光电子路径也使 XPS 对表面层更加敏感。相同 (001) 指数不保证相同的配位不饱和原子，吸附位点、功函数和表面电荷分布都可能随终止改变。建立吸附模型时，Sr 顶位、Ti 顶位和桥位的化学环境须分别对应所选终止。

比较两种终止时，还得核对模型两侧。上、下表面同为一种终止的对称 slab，原子数常偏离整数倍体相化学式；一个 SrO 面配一个 TiO₂ 面虽可保持体相计量，却把两种表面贡献与可能的法向电势差放在同一个模型里。常见的 γ = (E_slab − N E_bulk)/(2A) 只适合两面等价、计量关系明确的情形。非化学计量终止的表面能必须写明元素化学势，比较条件随生长气氛或化学环境而变。

体相坐标截出的最外层通常会发生表面弛豫，部分原子也可能重构或吸附环境中的分子；未弛豫结构用于定义初始终止，弛豫后的键长和层间距才进入能量比较。同一切割位置若允许不同重构，模型之间还需对齐面内超胞与原子计量，才有可比的表面能或吸附能。若表层原子发生横向位移，原始 (1×1) 周期还可能容不下实际重构单元。

极性切面为什么容易偏离理想截面？

在 ZnO 纤锌矿结构中，沿 c 轴的 (0001) 与 (0001) 对应相反的表面方向。理想化离子层可视作带相反形式电荷的 Zn、O 层依次堆叠；直接截取并周期复制，可能留下很大的法向偶极和内部电势变化。增厚真空仅拉大周期镜像距离，原有带电层堆叠仍在。表面可以通过空位、吸附物、电子转移或重构降低相应能量，实际路径取决于氧分压、温度和化学势。

Zn 终止 ZnO(0001) 的高分辨透射电镜图给出了一个具体情形：靠近表面的 Zn 原子列强度低于内部层，研究中结合成像模拟分析了厚度变化，支持表面 Zn 缺失的判断。原子列强度指向该实验条件下的局部占位变化；单个暗列还受样品厚度和成像条件影响，Zn 空位的分布需依靠多列统计与成像模拟约束。室温空气中的 ZnO 表面也可能带有羟基和吸附水，其终止条件另需测定。

原位图追踪 ZnO 在 300 ℃ 下的极性、半极性和非极性边缘：Zn 终止 (0001) 边缘出现准液态表层后，位置变化加快；其余边缘的时间轨迹不同。静态 DFT 裸表面模型描述的是给定原子排布的能量，这类随时间变化的蒸发还涉及缺陷生成、扩散和温度。图中的位移随时间变化是显微镜观测，蒸发快慢需结合 Zn 脱附、表面扩散及缺陷生成的速率分析；与计算比较时应对应终止、空位浓度和温度。

含偶极的非对称 slab 还会使真空区电势出现斜率，影响功函数以及不同覆盖度下的能量比较。看平面平均静电势时，应区分表面本身的电势阶跃与周期镜像间的虚假电场；检查 slab 加厚和真空加宽后的收敛，并明确是否使用偶极修正。极性问题来自终止、计量和电荷补偿方式；加厚原子层后仍需检查内部电场与总能量是否收敛。

计算中的平整晶面怎样对应实验表面？

XRD 给出的 (hkl) 通常是参与衍射的晶格面族信息；晶体取向或外形测量才能确定样品外表面的朝向。晶片相对理想晶面有小偏切角时，会形成由台阶和平台交替组成的表面，局部台阶边缘的配位不饱和原子与平台原子不同。偏切角越大，在相同步高下平均平台宽度越小；台阶还可能因退火、抛光或吸附而并合。对步高为 H、平均平台宽度为 L 的简单等间距几何模型，有 tanθ≈H/L；实际台阶聚集会使单一 L 难以代表整个表面。

蓝宝石 c 面抛光过程的 AFM 图保留了这种差别：左侧台阶较模糊，右侧出现更有序的平行台阶。AFM 量到的是一定扫描范围内的高度起伏，条纹间距可以对应平台宽度；它本身还不能判定台阶最外沿究竟是 Al 还是 O。若研究反应发生在台阶边缘，平整 Al₂O₃(0001) slab 会漏掉边缘配位环境，应另外考虑具有明确步高和台阶方向的模型。

晶体取向测量给出表面法向，AFM 或 STM 给出台阶与平台尺寸，LEED 可检查面内周期和重构，表面敏感的 XPS 可约束元素覆盖与终止。衍射面指数确定几何取向，显微图定位台阶，谱学约束表层组成。若目标是平台吸附，平整 slab 的面积与覆盖度要对应平台尺度；若目标是台阶活性，吸附位点应放在台阶边缘，并检查台阶间距是否足以隔开周期镜像。台阶模型沿边缘方向仍需保持合理的原子周期，否则重复排列的吸附物会产生并不存在的横向相互作用。