半导体表面吸附遇到带边弯曲怎么办？

引言

半导体表面吸附后，吸附物会改变表面电荷和局域偶极，费米能级、真空能级、价带顶和导带底不再像体相那样直观。带边弯曲不是单独看 DOS 峰位就能处理的问题，它需要把 LOCPOT 平面平均、真空平台、层分辨态密度和吸附态投影放在同一能量尺上。

能级对齐先选什么参照？

真空能级为什么常用？

slab 模型中，真空区平面平均静电势可作为能级对齐参照。若真空层足够厚，势能曲线会出现平台；价带顶、导带底和费米能级相对于该平台的位置，可以用于比较不同表面终止、缺陷或吸附前后的能级变化。

真空能级尤其适合比较同一类 slab 的吸附前后变化，因为吸附物带来的表面偶极会直接改变真空平台。若两个 slab 的厚度、终止面或偶极修正不同，平台差异就可能包含模型差异。此时应先统一 slab 构造，再讨论带边相对于真空的移动。

图1 把 ZnO 和 TiO₂ 的 Hartree 势与带边放在一起，说明带边位置必须依附于明确的能量零点。体相带边不能直接贴到任意 slab 图上，需要通过平均势或核心能级等参照转移。否则，一个看似上移的 VBM 可能只是零点选择改变。

图1 ZnO 与 TiO₂ 体相结构、平面平均 Hartree 势和价带/导带位置。该图用于说明带边对齐需要共同能量参照。 DOI: 10.3390/nano11061581。

费米能级能直接说明弯曲吗？

费米能级是电子占据的化学势，对金属和半导体的含义不同。半导体 slab 中，E_F 可能落在带隙、表面态或缺陷态附近；吸附物引入局域态后，E_F 移动不等同于整体带边弯曲。

图2 显示 slab 面向真空时的 Hartree 势剖面。判断带边弯曲应先看空间方向上的势能变化，再看层分辨 DOS 中 VBM/CBM 在不同原子层的投影。如果只看总 DOS，表面层和体相内部层的贡献会叠在一起，容易把表面态当作体相带边。

图2 ZnO 与 anatase TiO₂ slab 面向真空时的 Hartree 势剖面及带边位置。该图对应真空能级对齐的 slab 处理。 DOI: 10.3390/nano11061581。

吸附后带边为什么会变？

表面偶极怎样进入？

吸附物带有偶极矩或发生电荷转移时，会改变表面外侧静电势。对极性分子、卤素空位、氧空位和有机受体分子，局域电荷重排可使真空平台移动，也可使表面附近能带发生弯曲。此时功函数变化和带边变化需要分开记录。

图3 的界面 slab 把结构、Hartree 势和带边位置联合展示。吸附后的能级变化同时包含局域势变化和电子态重排。若吸附物只改变真空平台而内部带边相对平均势几乎不变，主要表现为表面偶极变化；若表面层 DOS 的 VBM/CBM 相对内部层改变，则可讨论带边弯曲。

图3 ZnO|anatase TiO₂ 混合 slab 结构、Hartree 势和带边位置。该图用于说明界面或吸附后局域势变化会影响带边读取。 DOI: 10.3390/nano11061581。

电荷重排怎样辅助判断？

差分电荷密度能显示吸附前后电子密度在哪些区域增加或减少。对半导体表面，它可以帮助判断电荷转移发生在吸附物-表面键、表面缺陷附近，还是延伸到多层 slab 内部。空间范围不同，带边弯曲和表面态归因也不同。

界面电子密度和差分密度剖面可用于核对电荷重新分布的空间范围。电荷重排的方向和范围应与 LOCPOT 势能变化对应。若差分电荷只局域在吸附键附近，而层分辨势能没有明显倾斜，则更适合讨论局域成键和表面偶极；若多层内部出现连续势能梯度，才更接近带弯曲描述。

图4 ZnO|anatase TiO₂ 界面的电子密度及差分密度剖面。该图说明电荷重排会改变界面偶极和局域势。 DOI: 10.3390/nano11061581。

表面态和吸附态怎样区分？

投影态密度看哪几层？

表面态通常局域在最外层原子或缺陷附近，吸附态则主要投影到吸附分子或吸附物原子。区分二者时，总 DOS 不够，应做原子层分辨 PDOS、吸附物分辨 PDOS 和关键能带的 PARCHG。对有机分子吸附半导体，HOMO/LUMO 与半导体 VBM/CBM 的相对位置尤其重要。

层分辨 PDOS 可以把最外层、次表层和内部类体相层拆开。若某个带隙态只在最外层出现，并随着 slab 厚度增加仍局域在表面，它更接近表面态；若内部层也出现连续贡献，则可能是带边整体移动。吸附物投影能进一步判断该态是否来自分子轨道。

若一个带隙中峰主要来自吸附物的 LUMO，它说明吸附态引入了新能级；若同一能区主要来自表面未饱和键，则更适合归为表面态。态的空间来源决定了能级图的物理含义。同一个能量峰不应同时被描述为体相导带、表面缺陷态和吸附物轨道。

投影能带怎样使用？

图5 中，分子、石墨烯和 SiC 的投影能带及轨道空间分布被放在同一图组里，适合说明吸附态和衬底态的区分。投影能带能看出某条能带主要来自吸附物、二维覆盖层还是半导体表面；PARCHG 能进一步确认该态在实空间的位置。

当吸附态与半导体带边强杂化时，单纯给 HOMO/LUMO 标签已经不够。需要比较吸附前后能带、PDOS、PARCHG 和差分电荷密度。强杂化态应被描述为吸附物-表面耦合态，而不是简单说成某个孤立分子轨道进入带隙。

图5 F6TCNNQ/graphene/SiC 表面体系的投影能带与 HOMO/LUMO 空间分布。该图用于区分吸附态、石墨烯态和半导体带边。 DOI: 10.1021/acs.jpcc.4c06409。

VASP计算怎样安排？

slab 模型怎样减少假弯曲？

半导体 slab 要先检查层厚、真空层和上下表面对称性。若 slab 太薄，两个表面会相互影响；若上下终止不同，内建偶极会让真空平台和内部势能都出现斜率。非对称吸附通常需要 LDIPOL 与 IDIPOL 设置，并检查偶极修正后真空区是否出现平台。

固定层也会影响势能剖面。底层固定过多时，内部原子层可能无法恢复接近体相的局域环境；全部放松时，薄 slab 又可能整体弯曲。层厚和固定层对照用于确认内部参考层是否稳定，然后再读取吸附造成的额外变化。

结构优化后执行静态计算并输出 LOCPOT、CHGCAR、DOSCAR 和 PROCAR。沿表面法向提取平面平均势，再把 VBM/CBM、E_F 和吸附态峰位对齐到同一个参照。所有参数，如 ENCUT、k 点、真空层和偶极修正，都应按体系收敛测试确定。

结果矩阵保留哪些量？

结果矩阵可列出清洁 slab、吸附态、缺陷 slab 和缺陷吸附态的真空平台、E_F、VBM、CBM、表面态峰位、吸附态峰位、Bader 电荷和偶极矩变化。能级移动的来源需要在这些记录项之间核对：真空零点变化、表面态引入、带边弯曲和吸附物轨道耦合对应不同的可观测特征。

若同一吸附物在完整表面和缺陷表面上都能稳定存在，应把两套 slab 的内部参考层分开检查。缺陷表面可能已经存在带隙态，吸附后新峰不一定来自吸附物本身。此时 PARCHG 能帮助定位该峰是局域在缺陷邻近原子，还是分布在吸附分子和表面成键区域。

对极性半导体表面，还要比较上下表面终止和 slab 厚度。若随着层数增加，内部势能斜率逐渐减小，说明原先的弯曲可能包含有限尺寸影响；若层数增加后斜率仍保留，并与差分电荷方向一致，带边弯曲解释才更有依据。

吸附态能级也不能脱离覆盖度。低覆盖度下单个偶极的局域影响可能较弱，高覆盖度下吸附物偶极会形成更连续的表面电场。若要比较覆盖度效应，应在同一超胞系列中改变吸附物数量，并同步记录真空平台和层分辨 DOS。

半导体表面吸附的能级解释，应从空间参照开始，再回到电子态来源。只有当 LOCPOT 显示势能沿 slab 法向变化、层分辨 PDOS 显示带边随层移动，并且差分电荷密度支持电荷重排方向时，带边弯曲的判断才更充分。若只看到 E_F 改变，应把结论限定为费米能级或表面态变化。

本章要点总结

要点：半导体表面能级对齐要先确定真空能级、平均势或核心能级参照。

要点：费米能级移动不能单独说明带边弯曲，需结合 LOCPOT 和层分辨 PDOS。

要点：吸附态、表面态和体相带边要通过投影能带、PDOS 与 PARCHG 区分。

要点：非对称 slab 应检查偶极修正、真空平台、层厚和上下终止面影响。