成键态与反键态全解析：能量来源、键强变化及缺陷判断方法

说明：本文华算科技主要介绍成键态和反键态的能量来源、电子占据如何改变键强，以及带边缺陷和晶格振动中的判断方法。

计算得到的一组 Kohn–Sham 本征态没有自带“成键”或“反键”标签。这两个词描述的是：某个电子态被占据后，它对指定原子对的相互作用能产生怎样的贡献。同一个能带可以对 A–B 原子对呈成键特征，对另一组原子对却接近非键态。讨论前应先说清原子对、轨道成分、能量区间和占据数。

成键态与反键态的能量来源是什么？

成键或反键来自轨道重叠后的能级分裂。在最简单的双原子近似中，两个可相互作用的原子轨道形成一低一高两组组合。低能组合在核间区域具有有利的相位关系，电子占据后可降低体系能量；高能组合在核间出现节点或相消，电子占据后抵消部分成键收益。这里比较的是同一组轨道相互作用前后的能量，轨道正负号本身没有可观测意义。

能级分裂幅度取决于能量匹配、空间重叠和对称性。两个轨道能量接近、指向合适且重叠充分，分裂通常更明显；能量相差很大或对称性不允许时，相互作用很弱。能量高低无法单独判定成键性质，判断对象仍是具体原子对在该本征态中的相位与哈密顿量贡献。

双原子模型常把键级写成 (N_b − N_ab)/2，其中 N_b 和 N_ab 分别是成键与反键轨道的电子数。这个计数揭示了占据差的作用，也隐含了轨道可清楚分组、电子局域在该原子对等条件。周期晶体的能带彼此交叠并含多中心杂化，此时积分 −COHP、Wannier 轨道或密度矩阵键级更适合量化指定原子对；不同方法的绝对数值不可混成同一标尺。

晶体中的本征态遍布布里渊区，成键特征还会随波矢 k 改变。一条能带在 Γ 点和布里渊区边缘可具有不同相位排列，元素投影态密度只能回答“由哪些原子轨道贡献”，无法独立给出“这两个原子之间成键还是反键”。实空间波函数、带分解电荷密度和成对的 COHP/COOP 分析各自补充相位、空间分布与能量分辨信息。

图 1b 中，GeSe 价带顶同时含有 Se 4p、Ge 4p 和可辨认的 Ge 4s 权重。图 1c 把这些成分还原为轨道相互作用：Ge 4s 与 Se 4p 耦合产生的高能组合已被电子占据，构成价带顶的一部分；导带底还有 Ge 4p–Se 4p 反键成分。带边名称给出能量位置，轨道耦合才给出成键性质，两类描述不能混写。

加电子和抽电子为何会让键长朝相反方向变化？

占据数把电子态与结构力联系起来。向成键态加电子通常增加恢复力，向反键态加电子通常削弱恢复力；从成键态抽走电子会损失稳定化能量，从已占据的反键态抽走电子则可能解除抵消。分子里的键级规则正是这四种变化的简化表达，晶体中还需计入电荷屏蔽、其他能带的重新杂化和全结构弛豫。

力的方向可由本征能对原子位移的导数理解。某个态的能量随键长缩短而降低，它被更多占据时会偏向短键；若该态的能量随原子靠近而升高，增加占据会偏向长键。结构优化读到的是所有占据态、离子静电作用和交换相关能共同形成的总力，单条轨道只解释其中一项。

HfO₂ 给出了一个多原子晶体中的具体例子。图 2a–c 显示空穴主要分布在 O 2p 价带，三配位 O^III 与四配位 O^IV 在不同晶相中的分配不同。图 2h 的 −COHP 表明，费米能级附近的 Hf–O^III 和 Hf–O^IV 成键贡献也不相同。同样抽走 0.2 个电子，空穴位于哪类氧及哪种 Hf–O 键上，会改变结构响应。

图 3 把电子态变化转成可计算的结构量。空穴浓度增加时，多个晶相的总能斜率和晶格参数发生分化，约 0.16–0.26 h⁺/f.u. 范围内极性正交相进入最低能区间；电子掺杂则使部分长 Hf–O^IV 键不稳定，并在较高浓度触发新的反极性或二维结构。键长突变和空间群改变说明体系已换到另一组局部配位，此后不宜继续沿用掺杂前的单键解释。

图 4b 中，空穴增加使极化翻转势垒持续下降；图 4c 中，电子掺杂先抬高势垒，接近结构转变后又快速降低，载流子浓度与势垒没有统一的线性关系。占据成键态或反键态可以说明局部力的来源，却不能代替完整能量面。比较不同电荷浓度时应同时检查晶相、磁态、体积、补偿电荷处理和费米能级位置。

价带顶是反键态，为何缺陷却可能较浅

“反键”这个标签只描述指定电子态对某组原子相互作用的贡献。若反键组合已占据，它会抬高价带顶；形成阳离子空位或抽走价带电子时，局部反键占据随之减少，缺陷能级可能留在价带顶附近。缺陷深浅取决于缺陷态相对带边的位置，宏观稳定还受其余化学键、缺陷形成能和竞争相控制。

GeSe 的价带顶含已占据的 Ge 4s–Se 4p 反键成分。图 1e、f 中，形成能最低的 V_Ge 转变能级距价带顶约 0.05 和 0.15 eV；形成深能级的几类缺陷往往具有更高形成能。高价带顶把受主态压在价带边附近，载流子较容易热激发回到延展态，深能级复合受到抑制。

带电缺陷的转变能级来自不同电荷态形成能的交点，数值还受超胞尺寸、静电有限尺寸修正、带隙误差和化学势范围影响。仅在价带投影中看到反键峰，尚不足以判定缺陷能级浅；应计算各电荷态形成能，并检查缺陷波函数是否延展、校正后的能级是否贴近价带顶或导带底。

GeSe 的体相浅能级与表面悬挂键属于两类局部环境。图 5c 的表面出现局域电荷，钝化后图 5d 的分布明显展开；C–V 与 DLCP 的差值对应的界面缺陷密度由约 2 × 10¹² cm⁻² 降至 2 × 10¹¹ cm⁻²。带边反键特征描述体相成键网络，表面配位不饱和需要单独建模。判断“缺陷容忍”时，体缺陷形成能、热力学转变能级和表面局域态应分开核对。

怎样把轨道占据与声子/力常数对应起来

静态键长只反映能量极小点，声子频率和原子间力常数还测量极小点附近的曲率。光激发能在晶格尚未来得及充分弛豫时改变价带与导带占据，适合检验电子态如何改写原子受力。若载流子从成键价带跃迁到反键导带，原有恢复力通常减小；具体声子变硬或变软仍取决于振动波矢和位移模式与这些能带的耦合。

静电掺杂模型加入或移除电子后，通常允许晶胞、原子坐标和电荷密度达到新的基态；瞬态光激发则可在固定电子—空穴数下求受限占据势能面。两类计算应分别报告载流子浓度、占据约束和弛豫自由度；光激发力常数对应非平衡电子分布，直接套到平衡掺杂结构会混入载流子寿命和晶格热化带来的差异。

PbTe 的价带和导带边缘具有离域的 p 轨道成分。图 6 用颜色标出光激发后的空穴和电子占据，并比较平衡结构与横光学声子位移后的能带。Γ 点位移使 L 点附近价带与导带分开，电子从价带移到导带后，这项畸变获得的电子能降低量变小，Γ 点软模随之硬化；X 点位移对带边的作用不同，对应的模式反而软化。相同载流子占据变化可对不同声子分支产生相反频移。

图 7 给出 216 原子超胞中前 11 近邻的谐性力常数。0.5% 光激发后，沿 方向关联的第一、第四和第八近邻大项明显减小，说明原有长程相互作用被削弱。此处的证据次序很清楚：受限占据 DFT 指定电子与空穴分布，能带位移识别电子—声子耦合，力常数和声子频率记录晶格响应。成键态与反键态的解释成立于指定占据、结构和振动模式，换泛函、载流子浓度或晶相后应重新计算，不能把标签当作固定材料常数。