描述符混用，是催化计算里最常见的隐性错误

引言

d带中心到底适合解释什么？

它为什么常和吸附强弱放在一起？

d 带中心的直觉来自金属表面吸附理论：金属 d 态越靠近 E_F，与吸附物轨道相互作用的能量窗口越容易打开，吸附往往更强。这个说法在同类金属表面、相似晶面和相似覆盖度下比较有用，因为模型之间的差异主要来自 d 态位置。

但它不是万能开关。d 带中心靠近 E_F 只能提示吸附可能增强，不能直接说明反应一定更快。若 *H、*O 或 *CO 已经吸得太牢，d 态上移反而可能带来毒化或脱附困难。描述符只负责给趋势线索，真正判断活性还要回到反应自由能台阶。

式(1)  ε_d = ∫E·D_d(E)dE / ∫D_d(E)dE

图1可以放在这里提醒：描述符不能脱离结构使用，同一个 ε_d 变化要先回到具体位点、吸附构型和反应环境。

图1  HER 催化剂结构设计与描述符关系。DOI: 10.3390/nano16120769。

哪些体系里 d 带中心容易失效？

还有一种常见失效情况，是不同吸附物对应的关键轨道并不相同。比如 *H 更关心金属 d 态和 H 1s 的耦合，*O 或 *OH 还会牵涉 O 2p 态和金属 d 态的成键/反键占据。若把所有中间体都用同一个 ε_d 解释，很容易把不同化学键压成同一种趋势。

对单原子位点尤其要小心“平均值很好看，局部峰很关键”的情况。一个体系的 d 带中心可能整体下移，但 E_F 附近仍保留一个尖锐 d_z2 峰；另一个体系 d 带中心略上移，却没有与吸附物方向匹配的轨道。描述符必须回到具体轨道和吸附构型，否则趋势线很容易误导。

单原子位点、强配位 M-N-C、氧化物缺陷和强自旋体系都容易让 d 带中心变钝。原因很简单：这些体系的 d 态不再像金属表面那样形成连续能带，而是被配位环境、自旋分裂和局域畸变切成几个窄峰。一个平均 d 带中心可能把最关键的局域峰平均掉。

这时更可靠的做法，是把 d 带中心当成第一层描述，再补 PDOS 分轨道、COHP 或具体轨道占据。比如某个 d_z2 峰正好靠近 E_F，而 d 带中心整体变化并不大，吸附差异可能主要来自这个局域峰，而不是整个 d 带的平移。

p 带中心和 EF 分别怎么用？

p 带中心更适合看哪些对象？

p 带中心常用于 O 2p、N 2p、S 3p 等非金属态，也常用于吸附物端原子的 p 态。氧化物里讨论晶格氧活性，氮化物里讨论 N 位点调控，或者吸附物 O/N/C 端轨道位置变化时，p 带中心比 d 带中心更贴近问题本身。

需要注意的是，p 带中心不是 d 带中心的替代品。它关注的是 p 态位置和占据情况，适合解释非金属参与成键、吸附物端原子受电子调控、或晶格氧是否更容易参与反应。若把 p 带中心上移也直接写成“吸附更强”，结论会过于粗糙。

图2可以用来说明载体结构和金属位点会同时改变电子态；只报一个 d 带中心或 p 带中心，都容易把结构因素藏起来。

图2  金属–载体结构对吸附强度影响。DOI: 10.3390/nano16120769。

E_F 是描述符，还是能量参照？

如果研究异质结、半导体或光/电耦合体系，E_F 还会和功函数、带边位置、界面电势降联系在一起。此时只把 DOS 图的 E_F 对齐到 0 eV 不够，最好同时给出真空能级对齐或平面平均电势结果，说明电子从哪个区域更容易流向反应位点。

E_F 首先是电子化学势，在 DOS 图里又常被设成 0 eV 作为能量零点。金属体系里，E_F 附近态密度可以提示电子供受能力；半导体和带隙体系里，E_F 位置还会和带边、缺陷态、吸附诱导态联系在一起。

真正容易出错的是横向比较。不同体系都把 E_F 平移到 0 eV 后，绝对能级差异可能被隐藏。若要比较半导体、异质结或带电吸附体系，最好再说明真空能级、功函数或带边位置。E_F 能参与解释电子供受，但不能脱离能级对齐单独横比。

式(2)  ε_p = ∫E·D_p(E)dE / ∫D_p(E)dE

多个描述符怎样组合？

先用哪个，后补哪个？

比较正确的顺序是先问反应问题。若研究金属位点对 *H、*CO 或 *O 的吸附趋势，可以先看 d 带中心；若讨论氧化物、硫化物、氮化物或吸附物端原子的 p 态变化，再引入 p 带中心；若讨论电子供受、功函数、带边位置，再把 E_F 放进来。

不要把三个量堆到同一个结论里。d 带中心解释金属 d 态与吸附物相互作用，p 带中心解释非金属或吸附物 p 态位置，E_F 解释电子化学势和能级参照。描述符越多，越要把因果链拆开。

图3适合放在组合描述符之后理解：吸附窗口不是单个数值决定的，而是结构、电子态和反应台阶共同决定的。

图3  电子结构调控与 HER 性能关系。DOI: 10.3390/nano16120769。

怎样把描述符落回 VASP 输出？

在 VASP 后处理中，d 带中心和 p 带中心通常来自投影 DOS。投影对象要写清楚：是金属中心、邻近 N/O/S 原子，还是吸附物端原子。能量区间也要固定，同一组对比不能一个积分到 E_F，另一个积分到更宽区间。

如果体系自旋极化，还要分别检查 spin-up 和 spin-down。很多单原子位点的关键差异不是总 DOS，而是某个自旋通道里出现了靠近 E_F 的局域峰。此时一个总的 ε_d 数字会显得干净，却可能丢掉真正的物理信息。

论文里怎样写？

哪些句子不建议写？

“d 带中心上移，所以活性提高”这种句子太短，也太容易被质疑。更合适的写法是：d 态靠近 E_F 后，目标吸附物轨道与金属态的重叠窗口增加，同时吸附自由能向更合适范围移动，因此支持该位点对某一步的调控。

p 带中心和 E_F 也一样。p 带中心变化要说明是哪个 p 态、对应哪个原子；E_F 移动要说明是相对真空能级、带边，还是仅在 DOS 图中作为零点。指标的名字不能替代物理对象。

图4可以放在结果表达前提醒：描述符最好服务于具体反应步，而不是在结尾堆成一串看似完整的电子结构解释。

图4  描述符与催化性能关联。DOI: 10.3390/nano16120769。

什么时候该放弃某个描述符？

放弃某个描述符，并不等于放弃解释。更好的做法是把原因写清楚：体系存在强局域态、自旋分裂或配位畸变，导致单一 ε_d/ε_p 无法线性对应吸附能；因此本文转向 PDOS 分峰、键长和反应自由能共同判断。这样的写法反而更像真实的计算分析过程。

如果同一族模型里 ε_d 与吸附能完全不相关，或者不同构型下同一个描述符给出相反趋势，就不要硬解释。可以换成局域 PDOS、COHP、差分电荷、键长和吸附能联合判断。描述符失效不是坏事，它反而提醒模型已经超出简单经验关系。

最终写法可以更实在：某描述符给出初步趋势，但该体系存在离散态、强配位或自旋分裂，因此进一步结合 PDOS 和成键分析解释差异。这样的表述比硬把所有体系塞进同一个 ε_d 关系里更可信。

图5表明：多个描述符可以共同出现，但每个描述符都要对应一个明确问题。

图5  多描述符与催化性能综合判断。DOI: 10.3390/nano16120769。

本章要点总结

要点1：d 带中心、p 带中心和 E_F 不是同一种量。

要点2：d 带中心更适合同类金属位点趋势，不能无限外推到所有体系。

要点3：p 带中心适合非金属态、吸附物 p 态和氧化物环境。

要点4：E_F 常是能量参照，横向比较前要先处理能级对齐。