很多人画 DOS 时会顺手把横轴零点放到 EF，图一下子变得整齐，但问题也常从这里开始。金属、半导体、干净 slab 和吸附后 slab 里的 EF 并不是同一种解释场景；有时它反映电子化学势，有时它被带隙、表面态、缺陷态或表面偶极牵着走。EF 不是脱离模型就能横向比较的绝对能级，先说明体系和参照，后面的图才有讨论基础。本文把这个量放回 VASP 输出、DOS 横轴和真空能级里重新梳理。
金属里为什么常看 EF 附近的态密度？
在金属里，能带穿过 EF，EF 附近有连续电子态。这时讨论 EF 附近的 DOS、d 态分布或吸附物轨道耦合，通常和电子供受能力有关。但这个判断默认模型是金属性，并且计算已经收敛到合理电子结构。金属体系中 EF 附近有没有可用态，往往比 EF 本身的一个数值更重要。
OUTCAR 里的 E-fermi 可以看作电子化学势在零温附近的数值。它适合和同一计算中的 DOS、能带、占据情况放在一起读，不适合孤零零拿出来和另一个材料直接比高低。
式(1) EF = μe(T → 0 K)
半导体里为什么不能把 EF 当成带边？
半导体或绝缘体有带隙，EF 可能落在带隙中，也可能被掺杂、缺陷或表面态推向价带顶或导带底。这时它更像载流子浓度和缺陷环境的线索，而不是某个可以直接发生反应的“最高占据能级”。
带边位置来自能带或 DOS 的边缘，EF 来自电子占据和电荷中性条件。如果体系有强表面态，EF 可能被表面态钉扎；如果做了掺杂或缺陷，EF 的靠近方向也会改变。所以半导体催化表面不能只截图一条竖线就下结论，至少要同时说明带隙、带边和表面态来自哪里。
图1里的 DOS 和能带对比适合放在这里看：材料从氢氧化物到羟基氧化物，带隙、态密度和能带位置都在变化。半导体中的 EF 必须和 EVBM、ECBM、缺陷态或表面态一起判断。
图1 Ni(OH)2与 NiOOH 的 DOS 和能带结构对比。DOI: 10.1038/s41467-023-38978-5。
为什么峰位能比，绝对能级不能直接比？
画 DOS 时把 EF 移到 0 eV，是为了让读图更方便。这样可以很快看出费米能级附近有没有新态、吸附物轨道是否靠近 EF、d 态或 p 态有没有重排。但这一步同时把绝对能量参照改成了各自的 EF。
如果 clean slab 和 ads slab 都各自设成 EF = 0，看起来像在同一把尺子上，其实未必。吸附物可能改变表面偶极，真空平台也可能移动。DOS 图中的 0 eV 是画图零点，不自动等于不同模型共享同一个绝对参照。
式(2) Eplot = E – EF
吸附后 EF 漂移能说明什么？
吸附、缺陷或界面接触会改变电荷分布，EF 也可能跟着移动。这个现象可以提示电子结构确实发生了变化，但它还不能单独说明吸附键增强、反应台阶降低或催化活性提高。
图2把体相和表面 PDOS、自旋态放在一起，更接近真实判断方式。先看 EF 附近新态从哪里来，再看这些态是否属于活性位、吸附物或表面重构。EF 移动只是线索，机制判断还要回到 PDOS、结构、电荷和能量变化。
图2 NiOOH 体相与表面态 PDOS 及自旋态分析。DOI: 10.1038/s41467-023-38978-5。
只看 OUTCAR 的 E-fermi 够不够？
OUTCAR 给出 E-fermi，DOSCAR 或 vasprun.xml 负责画总 DOS 和 PDOS，LOCPOT 才能帮助读真空能级。如果只是看同一个体系内的态密度变化，OUTCAR 加 DOSCAR 通常够用；如果要谈功函数、能级对齐或界面电势，就不能少 LOCPOT。
还有一个容易被忽略的细节是展宽设置。金属体系常用展宽帮助收敛，半导体或绝缘体则要避免展宽把带隙附近的占据涂得过宽。如果 ISMEAR、SIGMA 或 NBANDS 改动很大，E-fermi 和 DOS 近费米区形状都可能跟着变。遇到趋势反常时，先复查这些输入项，比急着解释“电子转移增强”更靠谱。
OUTCAR 中的 E-fermi、DOS 图采用的零点、LOCPOT 平面平均后得到的 Evac。OUTCAR、DOSCAR 和 LOCPOT 回答的是三个不同层面的问题，不能互相替代。
式(3) Φ = Evac – EF
功函数为什么必须回到真空平台？
功函数比较的是电子从材料内部电子化学势逃到真空所需的能量，所以它必须有真空能级参照。对于 slab，Evac 通常来自 LOCPOT 沿表面法向的平面平均电势平台。如果真空区没有稳定平台，功函数的数字就很难可靠。
图3里的表面 DOS 和 EF 标记可以提醒一件事：DOS 横轴和真空参照不是同一个东西。讨论功函数或不同 slab 的能级对齐时，必须把 Evac 和 EF 放在同一张能量尺上。
图3 Cu2O 表面与吸附体系 DOS 中费米能级标记示例。DOI: 10.1021/acs.chemmater.7b04803。
clean slab 和 ads slab 怎么做到同参照比较？
slab 层数、真空层、固定层、k 点、偶极修正、电荷状态和自旋设置都要对齐。否则 EF 的变化可能来自边界条件，而不一定来自吸附本身。
如果吸附物本身带有强偶极，或者吸附只发生在 slab 的一侧，尤其要检查偶极修正前后 LOCPOT 是否稳定。带电 slab 更要谨慎，因为背景电荷会改变长程电势，EF 的读数和真空平台都可能变得敏感。这类体系最好不要把一个 E-fermi 数值写成决定性证据。
更稳的做法是同时给出 DOS 变化、吸附构型、关键键长、电荷重排和能量变化。如果这些证据指向同一趋势，再讨论 EF 附近态密度对电子供受的影响。吸附前后 EF 的比较，前提是结构边界和能量参照保持一致。
什么时候必须补带边、功函数或偶极修正？
只想说明某个吸附态在 EF 附近出现新态，DOS 对齐到 EF = 0 通常够用。但只要结论涉及电子逸出难度、半导体带边位置、不同材料界面对齐，或者吸附导致的表面偶极变化，就应补功函数、带边或电势平台。
图4展示层数依赖 slab 的能带变化，适合提醒 slab 收敛问题。薄 slab 的带隙、表面态和带边位置可能还没有稳定，过早解释 EF 会把模型尺寸效应写成材料机理。当结论跨越不同 slab、不同表面或半导体带边时，EF 必须和带边、Evac 或偶极修正一起复查。
先说 DOS 或 PDOS 在 EF 附近出现了什么变化，再说明这个变化是否经过同参照能级对齐验证。如果没有 Evac 或带边位置，就把结论限定在“相对 EF 的态密度变化”，不要扩展到电子逸出能力或界面能级匹配。
式(4) EVBM F CBM
图4 层数依赖 slab 模型的能带结构与带隙变化。DOI: 10.1021/acsomega.4c06465。
要点：费米能级要放回具体体系，金属、半导体、slab 和吸附态的解释边界不同。
要点： DOS 图把 EF 设为 0 eV 是画图处理，不代表不同模型自动共用绝对零点。
要点：功函数、带边和界面能级对齐需要 Evac、电势平台或带边位置支撑。
要点：吸附前后 EF 变化只能作为电子结构线索，不能单独写成活性机制。
#华算科技朱老师#VASP吸附与催化计算教程#费米能级#DOS#能带结构#功函数#能级对齐
