d带中心计算

我第一次做 d带中心计算，是帮客户解释”为什么 Pd 比 Ag 更容易吸附 CO”。当时直接报了一串吸附能，客户追问”背后的物理是什么”，我答不上来。那次之后我学会先把 d 带中心算出来——d带中心计算给的不是一个数，而是一个能讲清活性来源的物理量。

一、为什么 d带中心能当活性标度

按 Newns-Anderson 模型，金属 d 带和吸附物轨道杂化形成成键/反键态。d 带中心 ε_d 越靠近费米能级，反键态空出越多，吸附越强。我算过一排 Pt 族金属，ε_d 从 Ni 的 −1.3 eV 到 Pt 的 −2.5 eV，CO 吸附能跟着从强到弱单调变化，趋势和实验吸附热一致。这比逐个位点算吸附能高效得多，是催化剂初筛的利器。

二、核心原理：ε_d 怎么积分出来

d带中心计算 的公式是 ε_d = ∫ E·ρ_d(E) dE / ∫ ρ_d(E) dE，其中 ρ_d(E) 是 d 轨道分波态密度。实操上我先用 VASP 的 LORBIT=11 输出每个原子的 lm 分解，再用脚本只抽取 d 分量，在费米能级附近（通常 −10~0 eV）数值积分。关键是积分窗口要覆盖整个 d 带，不能只取主峰，否则会系统性偏移。我做d带中心计算时习惯把窗口画出来确认 d 带起止，再积分。

三、关键技术要点：表面和位点要对应

同一金属不同晶面的 ε_d 差得明显，Pt(111) 和 Pt(100) 能差 0.3 eV 以上，吸附能随之变。我算合金或单原子时，还要区分中心原子和配位原子的 d 带。曾经一个 Pt₃Ni 位点，Ni 的 ε_d 比纯 Ni 上移，解释了对 ORR 的协同增强。所以 d带中心计算 必须和具体表面、具体位点绑定，脱离结构谈数值没有意义。

四、实操流程：从 DOS 到描述符

我的步骤：结构弛豫 → 加大 k 网格自洽 → 输出 d-PDOS → 用脚本按原子、按轨道积分 ε_d → 与吸附能做线性拟合验证（一般斜率负，即 ε_d 越上移吸附越弱）。如果拟合 R² 低于 0.8，我会回去看是不是混入了磁矩异常或结构未收敛的位点。这套流程跑一遍，十几个候选的活性排序就清楚了。

实算一个 Pd/Au 合金：纯 Pd 的 ε_d 约 −1.6 eV，掺 25% Au 后上移到 −1.9 eV，CO 吸附能相应从 −1.6 eV 降到 −1.1 eV，趋势和实验吸附热完全一致。更有意思的是，合金表面 Au 位点的 ε_d 比体相 Au 更负，说明配位环境把 d 带拉了下来，单看体相值会误判。所以我现在报 ε_d 一定带位点信息，不报”这个金属”的笼统值——同一金属不同位点，能差出活性档位。

五、常见踩坑：忽略填充度与轨道对称性

d带中心计算 只看中心位置会漏掉填充度。两个 ε_d 相同的金属，d 带填充不同，反键态占据也不同，吸附能可以差很多。我吃过用 ε_d 单变量解释不了 Pd/Pt 差异的亏，补上填充度后才闭合。另一个是 p–d 杂化强的氧化物，纯 d 带中心失效，要结合晶体场和 eg 占据一起看。

六、复盘：d 带中心是”翻译层”

回过头看，d带中心计算 最值钱的是它把复杂的吸附能压成一维描述符，让催化剂筛选从”逐个算”变成”按规律排”。我接项目时默认先要确定是否过渡金属表面，再决定用不用这个描述符。把 d-PDOS 提取、积分窗口、位点对应三件事做对，ε_d 才真的能指导实验。被证明有用的，是愿意为一根 d 带反复核对积分边界的严谨。下次有人只给一个数就说”活性高”，我会先问——这个 ε_d 是哪个位点、哪个晶面的。