

引言


HOMO/LUMO 投影电荷密度常出现在半导体、钙钛矿、分子晶体和缺陷体系的 DFT 结果中。它展示的不是总电子密度,也不是差分电荷密度,而是某个能级或能量窗口对应电子态在实空间中的分布。HOMO/LUMO 图回答的是带边态主要分布在哪里,而不是体系总共有多少电荷。
在周期性半导体里,更常用的说法是 VBM 和 CBM;在分子或有限团簇里,则常写 HOMO 和 LUMO。两者背后都指向前线电子态:价带顶对应最容易失去电子的占据态,导带底对应最容易接受电子的未占据态。图像能帮助判断带边来自无机骨架、有机配体、缺陷附近还是吸附物附近。

图1 低维极性杂化钙钛矿 IMP 的晶体结构。结构中的有机层和无机 Pb-I 骨架为后续 HOMO/LUMO 空间分布提供了不同区域。DOI: 10.1039/d0sc06112c。


HOMO/LUMO 图和普通电荷密度差在哪里?


普通电荷密度看什么?
普通电荷密度来自所有占据电子态的总和,主要反映电子云整体在哪里。它适合看成键区域、离子性/共价性差异和空间电子富集,但不容易告诉我们带隙边缘由哪些态控制。对于半导体光吸收、载流子分离、缺陷复合和界面电荷转移,带边态比总电荷更直接。
HOMO/LUMO 或 VBM/CBM 投影电荷密度则只选取特定能级或能量窗口。若 LUMO 主要分布在 Pb 原子附近,说明导带底主要由 Pb 相关轨道贡献;若 HOMO 主要分布在 I 原子附近,说明价带顶主要由 I 相关轨道贡献。投影电荷密度把能量信息投射到实空间。
它和差分电荷有什么不同?
差分电荷密度比较的是两个或多个体系相减后的电子重排,例如吸附前后哪里电子累积、哪里电子耗尽。HOMO/LUMO 图不做体系相减,而是从同一个体系中抽取某个电子态。二者可以互补,但不能互相替代。
图2 直接展示了 IMP 的能带、HOMO 和 LUMO 电荷密度等值面。HOMO 主要围绕 I 原子,LUMO 主要分散在 Pb 原子附近,这说明带边由无机骨架主导,而不是简单由有机层决定。带边态空间来源可以帮助解释带隙、光吸收和载流子迁移方向。

图2 IMP 的吸收谱、能带结构、HOMO/LUMO 电荷密度和电输运结果。图中 HOMO 主要位于 I 附近,LUMO 主要位于 Pb 附近。DOI: 10.1039/d0sc06112c。


带边态空间分布能说明什么?


离域带边意味着什么?
若 VBM 或 CBM 在晶体骨架中较为离域,说明该带边态具有较强的周期性延展。对于光电半导体,这通常有利于载流子在对应方向上传输;对于催化半导体,离域态可能更容易参与界面电子传递。但离域不必然代表迁移率高,需要结合能带曲率、有效质量和散射机制判断。
如果 VBM 和 CBM 分布在不同空间区域,例如一个偏有机层、一个偏无机层,可能形成空间分离,有利于降低电子-空穴复合;但过度分离也可能降低跃迁强度或界面反应耦合。空间分离本身不是好坏判断,必须结合目标功能。
局域带边意味着什么?
若带边态集中在某个缺陷、掺杂原子、表面吸附物或晶界附近,说明该局域结构对带隙边缘有强影响。对于光电材料,这可能意味着缺陷诱导复合中心;对于催化体系,这也可能意味着反应物诱导的新活性态。判断方向取决于能级位置、占据情况和反应对象。
ρn,k(r) = |ψn,k(r)|2 式(1)
式(1)说明投影电荷密度本质上来自某一能带 n、某一 k 点波函数的模平方。若选取的是一个能量窗口,则结果是多个态的加和。图像只说明该态的空间分布,不能单独给出电荷转移量。

图3 钴基硝酸盐转化催化剂的结构设计与机制示意。该图中的金属-配位结构可作为理解局域电子态与反应位点对应关系的例子。DOI: 10.1007/s40820-025-01877-z。


局域化是否等于缺陷复合中心?


为什么不能直接下结论?
带边态局域化经常被解读为缺陷态,但这一步需要小心。真正的缺陷复合中心通常要求缺陷态位于带隙内,且能够捕获电子或空穴;如果某个局域态只是价带顶的一部分,或者与主体带强烈杂化,就不能直接说明它会导致复合。
判断缺陷态时,应同时看缺陷前后 DOS、能带、PARCHG、形成能和电荷态。若缺陷引入孤立带隙态,并且 PARCHG 显示该态局域在缺陷邻近原子上,才可以说它具有缺陷诱导特征。局域化是缺陷态判断的证据之一,不是完整结论。
吸附态怎样区分?
催化或传感体系中,吸附物可能在带边附近引入新态。若 LUMO 或 CBM 由吸附物轨道贡献,说明未占据态可能参与受电子过程;若 HOMO 或 VBM 被吸附物占据态抬高,说明吸附物会改变空穴相关过程。此时需要把吸附前后能带、PDOS 和 PARCHG 放在同一能量参照下比较。
图4 和图5 展示了钴基催化体系中结构、电子态和反应路径之间的联系。虽然它们不是 HOMO/LUMO 专图,但能说明一个普遍原则:电子态空间分布必须回到具体结构和反应位点,才有机制意义。只看一张等值面图,无法判断该态是主体带边、缺陷态还是吸附诱导态。

图4 钴基硝酸盐转化催化剂的结构调控与电子结构表征。图中体现局域配位环境会改变反应相关电子态。DOI: 10.1007/s40820-025-01877-z。

图5 钴基催化剂上硝酸盐转化路径与机制分析。图中结构、电子态和中间体稳定性共同决定反应判断。DOI: 10.1007/s40820-025-01877-z。


怎样和 DOS、能带、Bader 联合判断?


DOS 负责定位能级来源吗?
DOS 和 PDOS 可以告诉我们带边由哪些元素和轨道贡献,例如 Pb-6p、I-5p、O-2p、过渡金属 d 态或有机 π 态。PARCHG 则把这些态放到实空间中。若 PDOS 说 CBM 来自金属 d 态,而投影电荷密度却主要分布在配体上,就需要检查能带编号、k 点和能量窗口是否选错。
能带负责判断态的位置、带隙和离散程度。平坦能带通常对应更强局域性,陡峭能带对应更强离域性;但最终仍要用 PARCHG 看空间分布。DOS、能带和投影电荷密度分别回答来源、能量位置和空间分布。
Bader 电荷能替代它吗?
Bader 电荷统计的是原子盆中的总电荷变化,适合看吸附或缺陷引起的净电荷转移。它不能直接告诉某个带边态长在哪里。相反,HOMO/LUMO 图也不能直接量化净电荷转移。联合判断时,Bader 给出电荷变化,PARCHG 给出态分布,PDOS 给出轨道来源。
PA = ∫ΩA ρedge(r) dr / ∫ ρedge(r) dr 式(2)
式(2)表示带边态在某个原子区域 ΩA 内的分布比例。实际工作中不一定都做积分,但它给出了定量比较 HOMO/LUMO 局域位置的口径:等值面图应与投影权重或原子分区积分结合。


VASP输出中怎样避免误读?


VASP 绘制带边投影电荷密度时,首先要确认体系是周期半导体还是有限分子。周期体系应从能带和 DOS 中定位 VBM/CBM 的能带编号、k 点和自旋通道;分子体系则可围绕 HOMO/LUMO 编号输出。若存在简并态,应把相关态一起纳入窗口,而不是只选其中一个能带。
其次要统一等值面阈值和视角。同一个体系比较 VBM 与 CBM,可以用相近阈值;跨体系比较时,等值面体积受归一化、超胞大小和态数影响,不能直接用“图更大”判断态更多。等值面大小受显示阈值控制,不等同于电子数量。
ρwindow(r) = Σn,k∈W fn,k |ψn,k(r)|2 式(3)
式(3)表示能量窗口 W 内多个态的投影电荷密度求和。若窗口太宽,主体带、缺陷态和吸附态可能混在一起;若窗口太窄,又可能漏掉简并态或同一带边的多个 k 点贡献。因此窗口选择要和能带、DOS、PDOS 对齐。
结果解释应限定在带边态空间来源和局域程度上,例如 HOMO 主要分布于 I 相关骨架,LUMO 主要分布于 Pb 相关骨架;或 VBM 局域于缺陷邻近 O 原子,CBM 仍保持主体骨架离域。这类图不能直接等同于价态变化、吸附增强或载流子寿命提高。HOMO/LUMO 投影电荷密度的结论边界,是空间来源和局域程度。
若讨论光吸收或光催化,还应把 HOMO/LUMO 分布与跃迁矩阵元和带隙大小分开。HOMO 与 LUMO 空间重叠较强,可能对应较强跃迁概率;空间分离较明显,则可能有利于电荷分离,但也可能降低直接跃迁强度。投影电荷密度只能提供空间重叠线索,不能单独替代吸收谱、介电函数或非绝热复合分析。


本章要点总结


要点:HOMO/LUMO 或 VBM/CBM 投影电荷密度展示的是特定电子态的空间分布。
要点:普通电荷密度、差分电荷密度和投影电荷密度回答的问题不同。
要点:局域化不直接等于缺陷复合中心,需要结合 DOS、能带、形成能和电荷态。
要点:等值面图不能单独推出净电荷转移、价态变化或载流子寿命。
夜雨聆风