



看 COHP 时,真正要盯的是三个层次:费米能级以下成键态有没有被占据,反键态有没有被占据,积分后的 -ICOHP 是否和反应中间体脱附/转化相匹配。-ICOHP 描述局部键相互作用强弱,不直接等于反应活性。




不同软件和论文里,COHP 的正负号可能会以 COHP 或 -COHP 的形式出现。常见画法会把成键贡献放在正方向,把反键贡献放在负方向,方便直接看峰。读图前先看坐标轴写的是 COHP 还是 -COHP,不然成键/反键会被完全读反。
好的读法不是背正负号,而是看图注和方法说明,再结合能量区间。COHP 的物理意义来自轨道相互作用和占据情况,坐标正负只是展示约定。同一批图可以比较,跨论文的数值不能随手相减。
EF 以下通常是已占据态,对现有键强贡献最大;EF 以上是未占据态,更像潜在可占据的反键或成键通道。若反键峰已经落到 EF 以下,说明电子占据正在削弱该键。若反键峰靠近 EF,吸附后少量电荷变化就可能明显改变键强。
因此,读 COHP 不能只看峰高,还要看峰在 EF 附近的位置。靠近 EF 的反键态往往比深能级成键峰更影响吸附稳定性变化。

图1 COHP/键级相关方法用于分析局部化学键的示意。 DOI:10.1021/acs.jpcc.8b08934




ICOHP 是把某一对原子之间的 COHP 在能量上积分,通常积分到 EF。在同一体系、同一种键、同一设置下,-ICOHP 越大,常说明该键的净成键作用更强。这对比较 M-H、M-O、M-C 等局部键很有用。
-ICOHP = -∫EF COHP(E)dE
这个式子有两个隐藏前提:积分上限通常取 EF,而且比较对象最好是同类键。若一个是金属–氧键,另一个是金属–碳键,数值大小不一定能直接排活性。
催化常常需要中间体既能形成,又能离开。M-*CO 键过弱,*COOH 可能难以稳定;M-*CO 键过强,*CO 可能中毒。-ICOHP 大只能说明局部键更强,不能告诉你这个强度是否刚好落在火山图顶部。
催化活性通常要求适中吸附,而不是最大键强。所以 -ICOHP 应该和 Eads、ΔG、反应能垒一起读。若 -ICOHP 变大但脱附能垒也变高,性能未必改善。

图2 不同键分析指标与轨道相互作用之间的关系。 DOI:10.1021/acs.jpcc.8b08934




COHP 常用 LOBSTER 等工具从 VASP 波函数和结构中重构局域轨道相互作用。这里的关键不是只跑完 WAVECAR,而是保证结构、POTCAR、能带投影和基组选择能互相对应。投影质量不好时,COHP 峰形和积分值都会变得不可靠。
计算后要检查 charge spilling。若 spilling 太高,说明局域基组没有很好还原平面波结果。COHP 不是 VASP 总能的直接输出,它依赖后处理投影质量,因此方法部分不能只写“计算了 COHP”。
不要把所有键都丢进去平均。先根据吸附构型选关键键,例如 M-H、M-OOH、M-CCO、O-H 或相邻金属–金属键。若体系有多个等价位点,可以分别算再说明平均方式。
对反应中间体,建议比较同一键在初态、中间态和过渡态附近的变化。这样可以看出成键/反键占据如何跟反应坐标走,而不是只拿最终态一个 -ICOHP 数字讲故事。

图3 HER 催化体系中电子结构与吸附强度的对照。 DOI:10.1093/nsr/nwaa204




如果问题是“为什么这个吸附构型更稳”,COHP 可以是很强的主证据;如果问题是“为什么催化活性更高”,它最好是电子结构证据链中的一环。前者看局部键,后者还要看自由能路径、覆盖度、竞争吸附和动力学。
写法上可以说“COHP 表明某键成键作用增强,从而有助于稳定某中间体”。不要直接写“COHP 证明催化活性最高”。COHP 能解释局部成键,不足以单独证明完整反应性能。
论文里建议搭配三类图:吸附构型标出被分析的原子对;PDOS 或轨道图说明 EF 附近态的来源;自由能图说明这个键强变化对哪一步有利或不利。这样 COHP 不会变成孤零零的一张电子结构图。
如果 -ICOHP 变大但 ΔG 反而偏离最优值,就应把矛盾写出来。很多好结论不是“越大越好”,而是“键强被调到不再过强或过弱”。

图4 局域成键变化需要与反应自由能和实验性能共同判断。 DOI:10.1093/nsr/nwaa204




要点1:读 COHP 前必须确认图中显示的是 COHP 还是 -COHP,避免把成键/反键读反。
要点2:-ICOHP 可比较同类局部键强弱,但不能直接等同于催化活性。
要点3:EF 附近的反键态占据往往决定吸附强度是否容易被电荷变化调节。
要点4:COHP 最好和吸附能、自由能路径、PDOS 或构型图一起解释。
