d带中心位置为何决定吸附强弱？

说明：本文华算科技主要介绍 d 带中心位置的计算来源、它怎样改变金属—吸附物轨道耦合与反键态占据，以及带宽、轨道对称性和电子计数怎样改变吸附能相关性。

d带中心位置从哪一组电子态得到？

过渡金属表面的价电子态包含展宽明显的 d 带。计算中常把活性位金属原子的 d 轨道投影态密度记作 ρ_d(E)，再取能量的一阶矩：ε_d＝∫Eρ_d(E)dE/∫ρ_d(E)dE。这个数值是所选 d 态在能量轴上的加权平均，峰高、带宽和各能区谱重共同参与积分。最高峰所在能量与 ε_d 没有固定重合关系。

以 E_F 为界选择积分上限时，覆盖全部 d 态得到的是整条能带的一阶矩，只积分占据态得到的则是已填充谱重的平均能量。表面吸附引起电子转移后，占据态积分还会受到电子数变化影响。同组材料采用相同能窗和相同轨道投影，ε_d 的位移才对应同一类电子结构变化。

文献通常以费米能级 E_F＝0 eV 作为参照。ε_d 数值增大，表示 d 态的谱重平均位置向 E_F 靠近；数值减小，则表示平均位置进入更低能区。投影原子、积分窗口、自旋分量和轨道分量一旦改变，ε_d 代表的电子态集合随之改变。同一材料中顶位原子与次表层原子的 ε_d 可以相差数百毫电子伏，表面平均值会掩盖真实吸附位点。

图 1 把 3d、4d、5d 过渡金属吸附原子分组后，ε_d 与吸附能形成各自的线性趋势。同一周期内，轨道空间伸展程度和电子计数较为接近，中心位置对吸附能差异的解释更集中；把三个周期混在一条拟合线上，散点离散度明显增大。这里的 ε_d 描述金属吸附原子的局域 d 态，吸附能则对应它与 MoS₂ 基底的结合。

d 带中心与吸附能出现相关，并非两个计算结果的经验巧合。ε_d 改变了金属 d 态与吸附物前线轨道之间的能级间距，进而改变轨道混合后成键态、反键态的位置和占据。吸附强弱正是这些占据态对总能量贡献的结果。

能级间距怎样改变轨道耦合/能级分裂

孤立吸附物具有离散分子轨道。它接近金属表面时，先与展宽较大的 s、p 态作用，原有能级展宽成共振态；随后再与局域程度更高的 d 态杂化。具有相同对称性的两个轨道混合后分裂成低能成键组合和高能反键组合。耦合矩阵元 V_ad 越大、吸附物能级 ε_a 与金属 d 态能量间距越小，分裂幅度通常越大。

ε_d 上移时，晚过渡金属的 d 态平均位置更接近 CO 的 2π*、O 的 2p 或 H 的 1s 等参与成键的能级，轨道混合增强。成键态向低能移动带来稳定化，反键态则被推向较高能量。ε_d 下移会增大能级间距，分裂减弱，金属—吸附物共价作用的能量收益随之缩小。

金属 s、p 态带来的展宽与正交化排斥也计入吸附能。吸附物波函数和金属占据态必须保持正交，短距离接近时会产生排斥能；d 轨道杂化提供的稳定化必须超过这部分排斥，强吸附才会形成。ε_d 主要追踪随金属种类变化的 d 态贡献，表面极化和长程静电作用并未写进这个单一数值。

-COHP 把这种杂化作用分解到具体原子对和能量区间。图 2 中 Au、Co、Cr、Zr 吸附原子与 Mo、S 的曲线形状差异很大，费米能级以下既包含成键贡献，也可能保留反键贡献。图中谱峰来自特定原子对的哈密顿布居，和 PDOS 的“有多少电子态”不是同一个量。

合金表面中，不同金属原子的配位数和近邻元素改变局域晶场与轨道杂化。图 3 的 Pt₂Fe 表面给出了多组 Pt 位点的自旋分辨 d-PDOS，ε_d 从约 -2.95 eV 变化到 -2.48 eV。位点中心越接近 E_F，CO 吸附能整体越负；其中的散点偏离则记录了局域几何、轨道分量和自旋态留下的差异。

反键态占据怎样转化为吸附能差异？

轨道分裂只给出能级位置，电子占据才把分裂转化为净成键能。低能成键态通常位于 E_F 以下并获得电子，高能反键态是否跨过 E_F，决定其中有多少削弱键合的电子。反键态上移并减少占据时，成键与反键贡献的抵消变小，吸附结构相对洁净表面与孤立吸附物获得更大的能量降低。

在成分和位点类型相近的晚过渡金属表面，ε_d 靠近 E_F 往往伴随更强的能级排斥，反键组合进入未占据区的比例增加，吸附能变得更负。ε_d 处于较低能区时，反键组合也处于较低能量，费米能级以下的反键占据增多，净键强下降。这个电子占据过程构成 d 带中心预测吸附趋势的直接物理环节。

不同吸附物调用的供电子轨道和受电子轨道并不相同。CO 吸附同时包含 5σ 向金属供电子与金属 d 态向 2π* 回馈电子，O、N、H 的前线轨道能量与对称方式又各不相同。同一次 ε_d 位移对各吸附物产生的分裂幅度并不相等，一种中间体的拟合斜率无法直接移给另一种中间体。

图 4 对反键贡献进行了数值积分。吸附原子—S 相互作用中，已占据反键贡献的比例升高，吸附能沿弱吸附方向变化；纳入吸附原子与多个基底原子的总反键贡献后，相关仍然存在。相关系数没有达到 1，表明 Mo、S 两类近邻、轨道径向伸展和不同元素的电子计数仍在改变总能量。

Pt₂Fe 上 CO 吸附的 -COHP 进一步区分 C 2p 与 Pt 5d 各轨道的作用。1π、4σ 和 5σ 相关峰分布在不同能区，4% 压应变后峰位与强度一起变化；电荷密度差显示 C—Pt 键区的电子积累和耗尽重排。ε_d 给出平均能级移动，-COHP 标出发生变化的原子对与具体轨道，两组量分别承担趋势描述与成键归因。

带宽/轨道对称性/电子计数怎样改变相关斜率？

ε_d 把整条 d 带压缩成一个平均值，两条中心相同的 d 带仍可能具有不同宽度和填充。宽带在更大的能量范围内与吸附物轨道相互作用，窄带则把谱重集中在中心附近。d 电子填充数决定 E_F 穿过能带的位置，也控制杂化后反键态获得电子的数量。双原子位点中，相邻金属还会产生新的金属—金属耦合，使单原子表面的拟合斜率发生变化。

图 6 把 d 带宽 W_d、填充度 f_d 和吸附自由能放进同一组双原子位点数据。较高反键态与较低反键占据对应强结合，反键态进入占据区则对应弱结合。以 f_d/W_d^1/2 组合带宽和填充后，多个吸附物的自由能趋势比单独使用 ε_d 更集中，说明中心位置只保留了能带形状的一阶信息。

压应变缩短金属原子间距，增强 Pt—Pt 与 Pt—Fe 轨道重叠，d 带通常展宽，中心位置和轨道分量一同变化。图 7 中 Pt₂Fe 从无应变到 4% 压应变时，Pt d-PDOS 的主峰重排，ε_d 向低能移动，CO 吸附明显减弱。0.35%—2% 区间的点接近一条局部趋势，4% 压应变点则偏离原有外推线，强结构扰动已经改变带宽与吸附构型。

顶位、桥位和空位呈现不同吸附趋势，因为参与耦合的 d 态受到轨道对称性筛选：顶位吸附主要调用沿表面法向伸展的 d_z²，桥位和空位则增加 d_xz、d_yz 或 d_xy 的贡献。总 d 带中心移动很小而对称轨道中心大幅变化时，吸附能仍可能明显改变；自旋分裂表面还会让两个自旋分量形成不同的成键与反键占据。

同一吸附物、相同吸附位型、接近的吸附物数目与表面积比例和一致的能量参照构成 d 带中心相关分析的比较范围。在这个范围内，ε_d 对能级间距的变化敏感，常能预测吸附能的相对顺序；材料跨越不同周期、配位环境、磁态或成键对称性后，带宽、填充、轨道分辨 PDOS 与积分 COHP 决定新的相关斜率。原有线性关系失去相同斜率时，材料族内重新拟合的相关系数才对应真实的吸附能差异。