金属-载体强相互作用怎么计算？

引言

金属-载体强相互作用常被用来解释单原子、团簇或纳米颗粒催化剂的稳定性和选择性，但在 VASP 计算里它不是一个单独输出量。一个可复查的 MSI/SMSI 判断，至少要同时回答三件事：金属与载体是否形成明确接触，界面是否发生电子重排，目标中间体的吸附自由能是否随接触方式改变。金属-载体强相互作用的计算核心，是让结构接触、界面电荷和吸附响应共享同一模型口径。这一步决定后续图表能否支撑“载体调控金属中心”的判断。

先把金属-载体接触面定义清楚

金属对象先分成哪几类？

金属对象可以是单原子 M₁、双原子 M₂、小团簇 Mₙ，也可以是连续金属薄层或纳米颗粒截面。载体可以是 Al₂O₃、CeO₂、TiO₂、ZrO₂、碳载体或含缺陷氧化物。图1给出 Cu 单原子在不同氧化物载体上的模型和性能背景，适合说明同一个金属中心放到不同支撑体上时，局域配位、金属价态和反应选择性会同步改变。

建模时不能只把金属放到一个表面位点后直接比较 E_ads。单原子模型应先枚举顶氧位、桥氧位、氧空位邻近位和金属暴露位；团簇模型还要检查团簇取向、接触原子数和是否发生铺展。MSI 计算的第一层变量是接触方式，而不是反应中间体。若 Cu₁/TiO₂ 与 Cu₁/CeO₂ 的 Cu 配位数不同，后续吸附差异就同时包含载体效应和几何配位效应。

图1 Cu 单原子负载在不同氧化物载体上的结构与反应性能背景。该图用于说明金属中心、载体种类和局域配位需要先在模型层面区分。DOI: 10.1038/s41467-025-57307-6。

界面结合能回答的是结构接触强度

结合能参照怎样写？

界面结合能用于比较金属与载体之间的结构稳定性。对单原子模型，可写作 M + support → M/support；对团簇模型，可写作 Mₙ + support → Mₙ/support。参照态要固定：金属孤立原子、金属团簇、金属体相化学势或前驱体化学势不能混用。若要讨论抗团聚，金属体相或团簇参照更有意义；若要讨论某一支撑位点是否锚定单原子，孤立金属原子参照更直接。

E_bind = E_M/support – E_support – E_M　　式(1)

式(1)采用越负越稳定的写法。若比较 Cu₁/Al₂O₃、Cu₁/CeO₂ 和 Cu₁/TiO₂，三者应使用相同泛函、PAW 赝势、ENCUT、k 点、自旋设置和色散修正；氧化物 slab 的层厚、真空层、固定层和偶极修正也要统一。结合能只能说明接触强度，不能单独证明吸附位点已经被电子结构调控。图2给出 Cu 电子状态的实验和计算证据，说明结合能之外还要确认金属电子结构是否随载体改变。

金属与载体接触面积增大时，E_bind 往往更负，但单个金属中心的反应可及性可能下降。团簇铺展构型可能更稳定，却会遮蔽部分金属位；单原子进入氧空位后锚定增强，也可能使金属中心过度配位。结构稳定和反应可及性需要分列记录，这样后续吸附能变化才有解释边界。

图2 不同载体调控 Cu 电子结构的谱学和计算证据。该图对应 MSI 中金属价态、局域配位和电子态变化的判读。DOI: 10.1038/s41467-025-57307-6。

Bader 电荷和差分电荷密度回答电子重排

电荷转移怎样绑定同一几何构型？

界面电荷分析应在相同几何构型下完成。操作时先优化 M/support，再做静态计算输出 CHGCAR、AECCAR0、AECCAR2 和 LAECHG 所需文件；Bader 分析得到 q_M 和载体邻近原子的电荷变化，差分电荷密度用同一超胞下的 M/support、support 片段和 M 片段相减。若拆分片段时原子位置或 FFT 网格改变，Δρ 图会引入数值差异。

Δρ = ρ_M/support – ρ_support – ρ_M　　式(2)

式(2)说明差分电荷密度不是普通电荷密度截图，而是界面成键前后的电子重排。图3包含 Bader 电荷、ELF 和电子结构分析，可用于解释 Cu 与载体之间的电荷偏移和局域键合。Bader 数值和 Δρ 图应共同指向同一片界面区域。若 Bader 显示 Cu 失电子，但差分图只在远离界面处出现颜色变化，通常要复查片段拆分和网格一致性；若差分图在 Cu-O 或 Cu-M 接触区域出现明显累积/耗尽，同时 PDOS 出现金属 d 态与载体 O 2p 态相互作用，MSI 的电子证据会更完整。

图3 Cu/氧化物载体中的 Bader 电荷、ELF 和电子结构分析。该图用于说明电荷转移、局域成键和载体调控不能分开判读。DOI: 10.1038/s41467-025-57307-6。

吸附能变化要回到同一反应中间体

为什么要选同一个中间体？

MSI 最终通常要解释某个反应趋势，例如 CO₂RR 中 *COOH、*CO、*CHO 的稳定性，HER 中 *H 和 H₂O 解离，或 CO 氧化中 *CO 与晶格氧反应。比较时应把吸附物放在同一类金属中心或同一类界面邻近位上，先枚举吸附构型，再用低能构型进入自由能路径。若一个模型使用 Cu 顶位 *COOH，另一个模型使用 Cu-O 桥位 *COOH，差异就包含吸附位点差异。

E_ads(*X) = E_M/support+X – E_M/support – E_X　　式(3)

图4给出 Cu 载体调控下的 CO₂RR 反应路径和中间体信息，适合把界面电子结构与吸附变化连接起来。若 Cu₁/CeO₂ 上 *COOH 更稳定，同时 Bader 与 PDOS 说明 Cu 周围电子态更利于 C-O 相关成键，这才构成“载体改变吸附”的证据链。相反，若吸附能改变但界面电荷没有对应变化，应检查几何位阻、覆盖度、溶剂模型或自旋态是否才是主因。吸附响应只有绑定同一中间体和同一位点口径，才能用于 MSI 归因。

把 *CO 和 *H 同时纳入吸附矩阵后，MSI 的反应结论必须经过竞争中间体复核。若载体让 *COOH 稳定，但也让 *H 更稳定，CO₂RR 选择性可能被 HER 竞争削弱；若载体降低 *CO 脱附自由能，同时保持 *H 不强吸附，则更容易解释产 CO 趋势。单个目标吸附物不足以决定材料排序。

图4 Cu/载体体系中 CO₂RR 中间体吸附和反应路径计算。该图用于说明 MSI 需要通过具体吸附物和自由能变化落到反应选择性。DOI: 10.1038/s41467-025-57307-6。

VASP 任务怎样组织成闭环

目录和数据矩阵怎样排？

一个可复查的任务树可以分成五层：第一层是 clean_support，包括不同终止面、氧空位或羟基化表面；第二层是 metal_support，枚举 M₁、M₂ 或 Mₙ 的接触方式；第三层是 static_charge，输出 CHGCAR、AECCAR 和 LOCPOT；第四层是 adsorbates，枚举 *COOH、*CO、*H、*OH 等目标中间体；第五层是 pathway，把自由能、Bader、PDOS 和 COHP/ICOHP 汇总。

参数可从 PBE 或 RPBE、D3 或 optB86b-vdW、PAW 赝势、自旋极化、ENCUT 至少不低于 POTCAR 推荐值的 1.3 倍、EDIFF = 10⁻⁵ 到 10⁻⁶ eV、力收敛 0.02 eV Å⁻¹ 量级开始测试，最终取值以体系收敛为准。图5把电子结构、反应选择性和产物结果联系起来，提醒 MSI 不是停在界面模型上的概念。结果矩阵至少保留 E_bind、q_M、Δρ 观察区域、E_ads 和关键 ΔG；只有这些量在同一模型口径下互相支持，金属-载体强相互作用才进入可计算、可复查的层面。

最后可增加一个结构稳定复算：把最低能吸附态去掉吸附物后重新优化 M/support，检查金属是否迁移、团簇是否重排、氧空位是否被填补。若反应前后金属中心身份已经改变，吸附能和界面结合能就不能继续放在原始 MSI 模型下解释。这个复算成本不高，却能把“电子调控”和“工作态结构变化”分开。

图5 Cu/氧化物载体体系中电子结构调控与 CO₂RR 产物选择性关联。该图用于说明 MSI 计算需要回到可观测反应结果和路径差异。DOI: 10.1038/s41467-025-57307-6。

本章要点总结

要点：MSI/SMSI 不是 VASP 的一个单独输出量，而是结构接触、界面电荷和吸附响应共同构成的判断。

要点：界面结合能适合比较金属与载体接触强度，但不能替代 Bader 电荷、差分电荷密度和 PDOS 证据。

要点：吸附能变化必须绑定同一个反应中间体、位点口径和覆盖度，才能用于金属-载体作用归因。