这个概念并不抽象,它直接决定键长为何缩短、吸附为何增强、掺杂为何改变力学与电学响应。
从计算角度看,成键与反键不是凭图像直觉命名,而是由波函数线性组合后的能量变化决定。若两个原子轨道以同相方式叠加,核间电子密度增加,体系势能下降,形成成键轨道;若以反相方式叠加,核间出现节点,电子密度被排空,体系能量升高,形成反键轨道。
后续华算科技无论讨论 CO 分子、过渡金属配位,还是二维材料吸附,这一判断标准都不会改变。

一、轨道如何形成?
最基本的出发点是原子轨道线性组合。设两个相近能量的原子轨道分别为 ψA和 ψB,则分子轨道可写成两个组合:一个是相加,另一个是相减。相加并不只是“数值变大”,其物理本质是波函数相位一致,核间区域的概率幅增强;相减则对应相位相反,核间区域相互抵消,因此形成节点。

图1:过渡金属氧化物中电子结构描述符。(a) 能带结构示意图,展示了M 3d与O 2p轨道形成成键(σ, π)、非键(σ0)和反键(σ*, π*)轨道的过程。(b-f) 为各类电子描述符应用示例。DOI:10.1039/d2cs01068b。
ψb= cAψA+ cBψB
ψab= cAψA− cBψB
式中,ψb表示成键轨道,ψab表示反键轨道,ψA与 ψB是参与成键的两个原子轨道,cA和 cB是组合系数,反映各原子轨道对分子轨道的贡献大小。若两轨道能量接近、空间重叠显著,则 cA与 cB往往都较大,轨道分裂更明显;若能量差过大,即使几何距离合适,杂化也会受限。
因此,成键与反键的出现依赖两个条件:一是对称性匹配,二是能量匹配。例如两个 s 轨道沿核间方向重叠,可形成 σ 成键与 σ*反键;两个平行 p 轨道侧向重叠,则形成 π 成键与 π*反键。
若轨道对称性不相容,例如一个沿键轴的 pz与垂直键轴的 px,即使能量接近,重叠积分也接近零,几乎不会形成有效成键。

二、如何判断强弱?
仅知道“有成键、有反键”还不够,真正决定键是否稳定的是占据电子究竟落在哪些轨道上。最常用的量是键级,它把成键电子与反键电子的竞争直接写成一个可计算表达式。成键轨道被占据越多,键级越高;反键轨道被填充越多,键级越低,极端情况下甚至导致键断裂。

图2:电催化描述符体系。(a) 氮吸附自由能与交换电流密度的火山曲线关系。(b) d带中心理论示意图,展示了金属s/d轨道与吸附质轨道耦合形成成键与反键态的过程。DOI:10.1039/d2cs01068b。
BO = (Nb− Nab) / 2
这里 BO 是键级,Nb是占据在成键轨道中的电子数,Nab是占据在反键轨道中的电子数。以 H2为例,两个电子全部进入 σ 成键轨道,Nb=2、Nab=0,因此 BO=1,形成稳定单键;若是 He2,四个电子会同时填满 σ 与 σ*,得到 BO=0,所以常态下不形成稳定共价键。
在实际材料中,键级并不总是整数,因为周期体系中的电子态是离域的,且轨道杂化常跨越多个原子。此时更实用的判断方式是结合投影态密度、差分电荷密度和COHP/ICOHP。
如果某一能量区间内,核间区域电荷累积且对应态位于费米能级以下,通常说明成键贡献占优;若费米能级附近出现明显的反键态占据,则结构常处于可调控但不完全稳定的边缘。
成键强弱还与轨道分裂大小密切相关。在线性组合近似下,成键与反键能级分裂可写成 ΔE≈2|HAB|,其中 ΔE 是轨道分裂能,HAB是两原子轨道间的耦合矩阵元。HAB越大,说明重叠越强、相互作用越显著,成键态更低、反键态更高。
对于典型共价键,分裂常在 1–10 eV 范围内;若仅为弱吸附或范德华相互作用,相关态分裂往往不足 0.5 eV。

三、计算中怎么看?
在第一性原理计算里,单看总态密度往往无法识别成键来源,因为不同原子、不同角动量通道会混在一起。更有效的方法是先用PDOS判断哪些轨道在同一能区发生杂化,再用分波电荷密度或能带分解图观察核间电子密度是否增强。若两个峰在相近能量处同时出现,并且对应波函数在核间连续分布,这通常是成键态的直接证据。

图3:Pt(111)表面吸附H原子的多维度电子结构分析。(a) Bader电荷分析。(b) 差分电荷密度。(c) 电子局域化函数(ELF)。(d) 分波态密度(PDOS)。(e) 晶体轨道哈密尔顿布局(COHP),直接量化成键与反键贡献。DOI:10.1039/d2cs01068b。
进一步地,COHP方法比单纯看 DOS 更接近“键的能量分解”。其思想是把能带能量按原子对投影,区分成键和反键贡献。常见作图中,负的 COHP 或正的 −COHP 峰通常代表成键贡献,接近费米能级的反键峰则提示该键可能因电子注入而削弱。
积分后的 ICOHP 可用于比较不同键强,数值绝对值越大,通常表示该原子对的相互作用越强。
一个常见误区是把“电子多”直接等同于“成键强”。事实上,若新增电子进入的是反键轨道,结果恰好相反。以过渡金属表面吸附 CO 为例,C 端 5σ 轨道可向金属空态给电子,增强吸附;同时金属 d 电子会回馈到 CO 的 2π*反键轨道。
适度回馈有利于吸附活化,但若 2π*占据过多,C–O 键会显著拉长,振动频率从气相约 2143 cm−1下降到 1800–2100 cm−1区间,这正是反键占据增加的实验与计算对应关系。

四、如何用于材料分析?
在材料研究中,成键与反键分析最直接的用途,是解释“结构变化为什么发生”。例如在二维过渡金属硫化物中,引入电子掺杂后,若费米能级上移并开始占据金属–硫反键态,常会看到键长增加、声子软化甚至相变倾向增强。相反,若空穴掺杂优先抽空反键态,则局域键会收缩,某些畸变结构反而更加稳定。

图4:电催化祈氧反应(OER)的不同反应机理示意图。(a) 吸附质演化机制(AEM)。(b) 自由能台阶图。(c-d) 晶格氧氧化机制(LOM)。(e) 氧配对机制(OPM)。各机制均涉及表面成键状态的变化。DOI:10.1039/d2cs01068b。
一个具体计算案例是比较 O 原子在金属表面不同位点的吸附稳定性。
若在顶位吸附时,O 2p 与金属 d 态杂化较弱,反键峰贴近费米能级并部分占据,则吸附能往往偏弱;若在空位或桥位吸附时,O 2p 与多个金属 d 轨道形成更充分的成键分裂,反键态被推高到费米能级以上,则吸附更稳定。
此类分析通常结合吸附能、PDOS、差分电荷和 ICOHP 一起给出,结论会比只报一个能量更有说服力。
归根结底,成键轨道描述的是电子如何降低体系能量,反键轨道描述的是电子如何抬高体系能量。计算工作真正需要把握的,不是名词本身,而是三件事:哪些轨道发生了耦合、这些耦合是否满足对称与能量匹配、费米能级附近究竟占据了成键还是反键态。
只要这条主线清楚,分子轨道图、能带、DOS 和化学键分析就能彼此贯通,而不是各说各话。
