硫空位怎么建模？

引言

硫空位常见于 MoS2、WS2、CoSx 和金属硫化物表面。很多模型看起来只是少了一个 S，但它们回答的问题并不一样：面内空位看基面如何被激活，边缘空位看低配位边界如何成键，吸附物补位看反应物是否会占据缺陷。硫空位不是一个统一标签，位置、配位和是否被补位决定了它的物理含义。

所以，做 VASP 模型时不要把“去掉一个 S”直接当成一类数据。真正需要分清的是：删掉的是基面 S、边缘 S、桥连 S，还是界面 S；缺陷是否保持空着；吸附物是落在空位旁边，还是直接填进空位。三种情况的形成能、吸附能和电荷重排口径都不同。

面内空位先回答什么

基面 S-V 看的是惰性表面如何被打开

对二维 MoS2 来说，完整基面通常比较惰性，暴露的是饱和配位的 S 层；删掉一个面内 S 后，附近 Mo 原子变成低配位，d 轨道更容易参与吸附。这个模型适合回答“基面缺陷能不能提供新位点”，不适合直接代表边缘催化位点。

Eform(VS) = Edefect − Eperfect + μS 式(1)

式(1)里的 μS 可以来自体相硫、S8、H2S 或给定硫化学势边界，具体选择要和研究环境一致。若只想比较同一 slab 中不同 S 位点，化学势项会相互抵消；若要讨论硫贫或硫富条件下空位是否容易出现，就必须交代 μS 的来源。面内硫空位首先是一个缺陷形成问题，然后才是吸附活性问题。

图1 硫空位限制的 Co-Mo 位点用于 CO2 加氢的结构与性能示意，说明硫空位可改变 MoS2 基面位点。DOI：10.1038/s41467-026-69780-8

为什么一个面内空位和两个面内空位不能混比？

单空位主要改变一个局域 Mo 三角环境，双空位可能形成相邻低配位 Mo 对、空位团簇或局域相变。若两个 S 相距很远，它更像低浓度缺陷；若两个 S 相邻，它已经是空位聚集。两者的形成能要按每个空位归一化，也要看空位间相互作用。

因此，比较吸附能时应把单空位、邻近双空位、远离双空位分成不同组。若把它们放进同一张排序表，强吸附可能只是空位数量增加带来的配位变化，并不代表某个单空位更优。

边缘空位为什么是另一类模型

边缘位本来就低配位

边缘 S 和面内 S 的环境不同。边缘处金属原子本来就更低配位，可能存在 S 端、Mo 端、重构端或被氢/羟基覆盖的端位。删掉边缘 S 后，得到的是边界重排位点，而不是基面缺陷位点。它更适合讨论纳米片边缘、台阶边界或小尺寸硫化物的活性。

在同一个 MoS2 体系里，面内空位和边缘空位可以同时存在，但它们不能直接用一个 Eads 排名说明谁更好。因为边缘位的参考表面、局域应变、固定层和硫化学势边界都可能不同。边缘 S-V 回答边界位点问题，面内 S-V 回答基面激活问题，两者不应混成一类缺陷。

图2 MoS2 中 Mo-edge 和面内双硫空位的 DFT 模型及 O2 解离吸附能对比，展示边缘空位与面内空位的口径差异。DOI：10.1038/s41467-026-69780-8

边缘模型要先说清楚哪一种边缘

二维硫化物边缘并不是一条简单线。不同切边会暴露不同金属/硫比例；同一边缘还可能有 0%、50%、100% 的 S 覆盖。若只写“edge S vacancy”，没有说明边缘类型和覆盖度，后续吸附能很难复现。

更好的做法是先固定一套边缘模型，再在同一边缘上删不同位置的 S。若要比较面内与边缘，至少要分开讨论：面内用 slab 或二维超胞的缺陷形成能，边缘用纳米带或台阶模型的局域位点能量，最后再把反应中间体分别放上去。

图3 Co 和 Ru 单原子在缺陷 MoS2 上的吸附模式，红圈标出 S 空位位置，可作为空位位点吸附差异示例。DOI：10.3762/bjnano.12.56

吸附物补位为什么不能当普通吸附

补位意味着空位身份发生变化

如果 O、OH、H、CO、N 或金属原子直接进入 S 空位，它就不再是“空着的硫空位”，而是补位缺陷或替代位点。比如 O 填进 VS 后，更像 OS 替位；金属原子嵌入 VS 后，可能成为锚定位；H 或 OH 占据空位时，则可能形成反应中覆盖态。

ΔEfill = Efilled − Evac − μX 式(2)

式(2)回答的是 X 是否愿意填进空位，而不是普通表面吸附强度。若 X 只是停在空位旁边，可以用吸附能；若 X 与空位形成替位式结构，就应把它作为新缺陷状态单独命名。吸附物补位会改变缺陷本身，不能和旁位吸附放在同一列比较。

图4 MoS2−x 中三类空位缺陷模型示意，提醒不同空位类型对应不同局域结构和电子态。DOI：10.1039/d1ra04902j

同一批模型怎样保持可比

先分组，再比较

建议把硫空位模型分成四组：完整表面、面内 VS、边缘 VS、补位型 VS+X。每组内部可以比较不同位点、覆盖度或吸附构型；组与组之间只能做机制对照，不宜直接用一个最低吸附能下结论。

计算设置上，面内缺陷要检查超胞大小，防止空位周期性镜像相互作用；边缘模型要检查纳米带宽度、边缘覆盖和真空层；补位结构要做多个初始构型，防止吸附物滑出空位或形成另一种低能结构。若空位附近发生重构，要按优化后结构重新命名。

电子结构证据看什么？

硫空位常会引入缺陷态、改变 Mo d 态占据，并在空位附近形成电荷重排。可以结合 PDOS、差分电荷密度、Bader 电荷和键长来看。若只看到吸附能变强，但缺陷态、成键区电荷和吸附构型都没有明确变化，就要谨慎解释。

图5 MoS2−x 不同缺陷的差分电荷密度示例，显示 S-V 会造成局域电子重新分布。DOI：10.1039/d1ra04902j

结果怎样表达更清楚

不要只写“硫空位增强活性”

更好的表达是说明哪一种空位增强了哪一步。例如：面内 VS 暴露低配位 Mo，增强 CO2 或 H2O 的初始吸附；边缘 VS 改变边界金属配位，降低某一步断键能垒；O 填入 VS 后形成 OS 位点，改变后续 *OH 或 *OOH 稳定性。

最后要把形成难度和反应收益分开。一个空位很活泼，不代表它容易稳定存在；一个补位结构很低能，也可能说明原来的空位在反应环境中会被填掉。硫空位建模的核心，是用分组方式把“空位在哪里”“是否仍为空位”“吸附物有没有补位”讲清楚。

若需要把硫空位和反应路径连接起来，还应保留空位前后同一中间体的构型继承关系。比如先在完整表面、面内 VS、边缘 VS 上分别放 *H 或 *COOH，再检查优化后是否仍停在原位点；如果中间体滑入空位、补位或诱导 S 迁移，就要把它从普通吸附组移到补位/重构组。这样整理后，吸附能差异才不会被模型身份变化掩盖。

还可以把结构截图按同一视角保存下来，尤其是空位周围第一配位层。面内空位更关注局域配位缺失和电子态变化，边缘空位更关注边界配位数和边缘重排，补位模型则要关注新原子是否真正进入空位。图像、能量和电荷结果对应起来，后面解释活性差异会更清楚。

本章要点总结

要点一：面内硫空位、边缘硫空位和补位型缺陷回答不同问题，不能混成一类吸附位点。

要点二：单空位、双空位和空位团簇要分组比较，必要时按每个空位归一化形成能。

要点三：吸附物进入 VS 后，模型身份可能从普通吸附变成补位缺陷或替代位点。

要点四：硫空位结论最好同时给出形成能、优化后构型、吸附能和电子结构证据。