晶界不是普通台阶位的另一个名字。台阶位通常出现在同一晶面上的高度突变,而晶界来自两个晶粒取向拼接后的边界区域。边界附近会出现低配位原子、拉伸或压缩应变、空位富集、杂原子偏聚和局域电荷重排。晶界催化模型先要定义边界区域,而不是把所有低配位原子直接当成同一类活性位。
在 VASP 里,晶界常用周期超胞表示。模型里至少有一个晶界核心区和两侧近似晶内区域;如果超胞太窄,两个周期晶界会相互影响,吸附能里就会混入边界-边界耦合。晶界模型的可比性,取决于晶界类型、周期长度和晶内缓冲层是否说清楚。

图1 Au 晶界网络与局域结构示例。图中展示晶界区域和晶内区域具有不同原子排列。DOI: 10.1038/s41467-024-52919-w。
晶界区域到底怎样定义? 

晶界常按两侧晶粒取向差、旋转轴和边界面来描述,例如对称倾斜晶界、扭转晶界或高角度晶界。对催化计算来说,第一步不是马上放中间体,而是把原子分成三类:晶界核心原子、晶界附近过渡原子、远离晶界的晶内原子。吸附位点必须先绑定到这三类区域之一。

图2 二维材料晶界用于 HER 的结构示意。图中可见晶界处原子排列和配位环境不同于晶内区域。DOI: 10.1038/s41467-019-13631-2。
实际建模时,可以沿晶界方向保留足够长的周期,让边界上出现多个非等价位置。若只截取一个最小重复单元,可能把真实边界中的位点多样性压成一个位点。晶界不是一个点,而是一条或一片局域结构不均匀的区域。
EGB= (Ecell with GB− N·Ebulk) / AGB式(1)
式(1)给出晶界形成能的一种常见写法。AGB是晶界面积,N 是超胞中原子数。这个量不能直接替代吸附能,但它能告诉我们某一类晶界是否容易作为结构候选存在。若一个晶界形成代价很高,就需要解释它为什么能在样品或反应条件下出现。
低配位原子为什么不等于晶界活性位? 

晶界里常有低配位原子,但低配位只是一个结构特征,不自动等于活性位。真正需要比较的是中间体和哪些原子成键、键长是否合理、吸附后边界是否保持原有结构,以及同类晶内位点的参照吸附能。低配位原子只提供候选位置,吸附构型才决定它是不是反应位。
例如 *H、*CO、*OH 可能优先落在晶界的桥位,也可能被旁边晶内原子共同固定。若中间体跨过晶界两侧晶粒成键,这类位点就不同于单边低配位原子。晶界位点命名要看成键范围,而不是只看某一个原子的配位数。
Eads,GB= EGB+X− EGB− EX式(2)
式(2)适合在同一个晶界超胞里比较不同 X 吸附位。若要比较晶界和晶内吸附,最好在同一个超胞中选择远离晶界的晶内参照位,或另外构造晶内表面模型并说明应变、层数和覆盖度差异。晶界增强不能只靠一个最低吸附能,需要晶内参照。

图3 Au 晶界结构与反应性能关联示例。图中说明晶界密度和局域结构会影响催化表现。DOI: 10.1038/s41467-024-52919-w。
缺陷富集会怎样改变吸附? 

晶界常是空位、掺杂原子、杂质和应变集中的位置。缺陷富集后,吸附能变化可能来自两部分:一部分是晶界本身的低配位和应变,另一部分是缺陷或掺杂带来的局域电子态。如果不拆开比较,就很容易把缺陷效应写成晶界效应。
建议至少准备三组模型:干净晶界、有缺陷晶界、晶内缺陷参照。干净晶界告诉我们边界结构本身的影响;有缺陷晶界告诉我们边界富集缺陷后的变化;晶内缺陷参照则说明同一个缺陷离开晶界后效果是否还存在。晶界和缺陷要做交叉参照,而不是只算一个复杂模型。
Eseg= Edopant/GB− Edopant/bulk-like式(3)
式(3)可以用来描述掺杂或缺陷向晶界偏聚的倾向。若 Eseg为负,说明该缺陷更愿意停在晶界附近;若吸附增强只在 Eseg很有利时出现,就可以把结论写成晶界富集缺陷后形成新活性中心,而不是单纯晶界本征活性。

图4 Cu 晶界引导 CO2还原产物选择性的示例。图中展示晶界结构可改变 C-C 偶联相关过程。DOI: 10.1039/d3sc02647g。
晶界吸附计算怎样排? 

比较顺序建议从简单到复杂:先优化干净晶界,确认晶界没有异常闭合或大幅重排;再枚举晶界核心、晶界邻近和晶内参照位;随后加入缺陷、掺杂或覆盖物,观察中间体是否仍停留在原位。晶界计算最怕初始位点和优化后位点已经不是同一类位置。
对多原子中间体,还要检查吸附方向。CO2、*OOH、*CHO 或 *COCO 可能沿晶界方向排列,也可能横跨晶界两侧。若只放一种方向,会漏掉由晶界线性结构带来的特殊构型。晶界是一维延展位点,吸附方向本身就是变量。
如果晶界超胞里有两条周期晶界,还要确认它们是否等价。两条晶界不等价时,一个中间体放在上边界和下边界可能给出不同能量;若只报告一个位置,就会把周期构造带来的差异误写成晶界本征差异。晶界超胞里的每一条边界都要先确认身份。
还有一个细节是固定层。晶界模型往往比普通 slab 更容易发生整体剪切或局域滑移,固定过少会让晶界形貌改变,固定过多又会压住边界弛豫。更清楚的做法是报告固定区域离晶界多远,并比较关键吸附能是否随固定策略明显变化。

图5 高密度晶界 Mo2C 纳米片示例。图中展示晶界密度提升会带来更多低配位和缺陷相关区域。DOI: 10.1038/s41467-022-34976-1。
结果怎样判断? 

结果不要只写晶界吸附更强。更清楚的说法是:某个中间体在晶界核心桥位比晶内位更容易吸附,同时优化后边界结构没有明显破坏,且缺陷偏聚能说明该缺陷可能在晶界富集。晶界结论至少要同时回答位置、成键、结构保持和缺陷来源。
若晶界位吸附过强,可能导致中间体难脱附或边界重构;若吸附只强一点,但晶界密度高,则可能通过位点数量贡献活性。晶界优势既可能来自单个位点强度,也可能来自位点数量和缺陷富集。这几类作用不要混在一句话里。
跨材料比较时,应避免把一种材料的高角度晶界和另一种材料的低角度晶界直接排序。更可靠的做法是先在同一材料内完成晶界类型比较,再比较相同晶界口径下不同材料或不同缺陷的趋势。
如果晶界吸附能和晶内吸附能差别很小,不建议强行说晶界是主活性中心。可以写成晶界提供了额外候选位点,但其贡献还需要结合晶界密度、位点暴露比例和反应能垒来判断。晶界位点不是只靠一个能量差就能完成归因。
本章要点总结 

要点一:晶界模型要先定义晶界核心、近邻区域和晶内参照区域。
要点二:低配位原子只是候选位,真正的反应位要看优化后成键范围。
要点三:缺陷富集需要用干净晶界、缺陷晶界和晶内缺陷参照拆开判断。
要点四:多原子中间体在晶界上要考虑沿晶界和跨晶界两类方向。
要点五:晶界结论应同时说明位点身份、吸附强度、结构变化和缺陷来源。
