模型建错=三个月白干，催化计算的第一步，不是直接算吸附能！

单原子催化计算的核心挑战

单原子催化（SAC）以其100%原子利用率和独特的电子结构，成为近五年催化领域增长最快的研究方向之一。Nature Catalysis、JACS、Angew Chemie等顶刊上，SAC相关论文几乎每期可见。与此同时，审稿人对计算验证的要求也在持续提高——仅凭实验表征已不足以支撑一篇顶刊论文的构效关系论证。

与CO₂RR等以"路径选择"为核心的计算需求不同，单原子催化计算面临一个更基础的挑战：活性位点的真实配位环境是"看不见"的。XAFS仅能提供平均配位数和键长，无法区分M-N₄平面构型与M-N₃C₁轴向构型——而这一差异足以改变整个电子结构分析的结论。如果在第一步就建错了模型，后续所有计算都是无效的。

单原子催化计算不是"算一个吸附能"——它是从实验信号反推原子尺度真实结构、建立配位构型→电子结构→催化活性完整逻辑链的系统工程。每一步都直接影响审稿人是否接受你的构效关系论证。

测试狗自2019年起开展单原子催化计算，至今已完成150余个SAC课题，其中12篇支撑成果发表在Nature Catalysis、JACS、Angew Chemie、Advanced Materials等顶刊。基于这一积累，我们形成了从"结构确认"到"活性解释"再到"理性设计"的系统化解决方案。

解决方案：三步递进计算体系

本方案围绕单原子催化论文中最核心的三个科学问题设计，每一步对应一个明确的论文论证需求，既可独立使用，也可递进组合。

步骤一 配位结构解析

解决的问题：催化剂的真实活性位点结构是什么？

实验表征（XAFS、XPS）只能给出统计平均信息。本模块通过计算手段将实验数据转化为精确的原子尺度结构模型，为后续电子结构分析和反应路径计算建立可靠基础。

技术方案：

基于XAFS拟合数据（配位数、键长），构建2-3种候选配位构型，逐一进行DFT结构弛豫优化；

计算每种构型的形成能，给出热力学稳定性排序；

模拟XANES谱图并与实验谱直接对照，以谱图偏差最小者为最优构型；

输出最优构型的Bader电荷分布，量化金属-载体电荷转移。

交付成果：候选构型形成能对比表、XANES模拟-实验对照图（含偏差数据）、最优构型的完整结构文件与电荷分析。

典型应用：某课题组Co-N-C催化剂，原始假设为Co-N₄平面构型，XANES模拟偏差8.7%。经本模块重建为Co-N₃C₁后偏差降至2.3%，后续ORR活性计算与实验完全吻合。

步骤二 电子结构与活性起源

解决的问题：为什么这个催化剂好（或不好）？构效关系是什么？

这是审稿人最关注的部分。仅给出PDOS图和d带中心数值不足以建立令人信服的构效关系——需要从多个维度完整刻画电子结构，并将各维度分析串联成因果逻辑链。

技术方案：五层递进分析体系

PDOS分态密度——确定金属d轨道与吸附分子轨道间的杂化特征；

d带中心理论——解释吸附能趋势，建立金属-吸附物相互作用强度的定量标度；

自旋态分析——判断是否存在自旋转变驱动的额外活化机制（单原子催化区别于团簇催化的重要特征）；

COHP/ICOHP键合分析——定量计算金属-吸附物化学键的成键态与反键态填充，直接反映键的强弱；

差分电荷密度——可视化电荷在空间的重新分布，展示电子转移路径。

五层分析输出后，整合为"配位结构→电子结构→催化活性"的逻辑链——每一环有量化数据支撑，而非定性描述。

交付成果：PDOS组合图（含d轨道投影与吸附分子轨道）、d带中心对比表、自旋态分析报告、COHP/ICOHP键合分析、差分电荷密度图。所有图表按顶刊格式输出，可直接用于论文。

典型应用：Fe-N₄与Co-N₄ ORR活性对比研究。常规分析仅关注d带中心差异（-1.82 vs -2.31 eV），本方案进一步揭示自旋态转变机制——Fe-N₄在O₂吸附后发生低自旋→高自旋转变，释放~0.3 eV额外驱动力，而Co-N₄始终维持低自旋。这一发现与COHP键合定量共同构成三环证据链，发表于JACS。

步骤三 反应路径与高通量筛选

解决的问题：完整的催化机理是什么？如何从候选体系中筛选出最优催化剂？步骤三包含两个独立子模块，根据课题需求灵活组合。

子模块3A：反应路径完整计算

Climbing Image NEB过渡态搜索，覆盖多步质子-电子转移反应的全路径（ORR四步、NRR六步、OER四步等）；

构建自由能台阶图，标注各步ΔG与决速步，加恒电势修正（CHE模型）后可直接与实验LSV曲线对标；

双原子催化剂的协同效应量化——分别计算M₁-only、M₂-only和M₁M₂-together三种模式，定量分离协同增益。

子模块3B：数据驱动的高通量筛选

基于150+课题积累的结构数据库（3000+种配位环境），使用d带中心、配位数、Bader电荷和关键中间体吸附能作为描述符；

高斯过程回归快速筛选，典型效率：从200种候选构型缩减至5个推荐位点（1周内完成）；

筛选结果经DFT精算验证后输出，附带ML特征重要性排序，便于论文中讨论"为什么这些描述符起决定作用"。

交付成果：完整自由能台阶图（含电位修正）、NEB能垒柱状图（附虚频验证数据）、构型-吸附能热图、火山图预测、ML特征重要性分析。

典型应用：① 双原子Zn₂-N₆ CO₂RR——计算四种Zn-Zn间距下的C-C偶联能垒，间距3.2 Å时能垒从1.2 eV降至0.58 eV，COHP揭示协同机制为双金属d轨道通过桥连N的间接耦合，发表于Nature Catalysis（IF 41.8）。② 36种M-N₄ NRR高通量筛选——火山图定位Mo-N₄和Ru-N₄为最优，排除Pt/Rh（HER竞争过强），发表于ACS Catalysis（IF 13.5）。

方案优势

以下从六个维度说明测试狗单原子催化计算方案与常规计算服务的差异。

代表性案例

案例一：Fe-N₄酸性ORR活性起源——自旋态转变的电子结构证据

课题背景：客户合成了Fe-N-C和Co-N-C两种单原子ORR催化剂，Fe-N-C的酸性ORR活性显著高于Co-N-C。审稿人要求从电子结构层面解释"同一M-N₄配位下活性差异的根源"。

计算内容：

PDOS分析：Fe-N₄ d带中心-1.82 eV，Co-N₄ d带中心-2.31 eV，前者更接近费米能级，O₂吸附更强；

自旋态分析（关键发现）：Fe-N₄在O₂吸附后发生低自旋→高自旋转变，电子自旋翻转释放~0.3 eV驱动力，加速O-O键活化；Co-N₄始终维持低自旋，无此增益；

COHP键合定量：Fe-O键的成键态填充高于Co-O，O-O反键态填充更低；

三组证据链整合为"d带中心→自旋态转变→O-O活化→ORR活性"的完整构效关系。

发表结果：发表于JACS。

案例二：Zn₂-N₆双原子CO₂RR——协同机制的量化揭示

课题背景：客户合成了双原子Zn₂-N₆催化剂，实验显示其对CO₂→C₂₊产物具有选择性，但无法解释两个Zn原子的协同机制。

计算内容：

构建Zn₂-N₆四种候选间距构型（2.5~4.0 Å），计算C-C偶联能垒；

间距3.2 Å时能垒从单原子Zn-N₄的1.2 eV降至0.58 eV，为最优几何；

COHP揭示协同本质：两个Zn的d轨道通过桥连N原子形成间接耦合，向*OCCO中间体π*轨道反馈电子，稳定过渡态；

NEB+虚频验证+ICOHP键合定量，方法论完整。

发表结果：发表于Nature Catalysis。

案例三：M-N-C NRR催化剂高通量筛选——36种金属火山图

课题背景：客户拟开展单原子NRR研究，文献中Fe、Mo、Ru、W均有报道，不知如何选择金属中心，需理论预测支撑。

计算内容：

构建36种过渡金属M-N₄构型，高通量DFT计算N₂吸附能与*N₂→*N₂H自由能；

发现标度关系，Mo和Ru偏离标度线——即Sabatier原理的最优区间；

火山图预测：Mo-N₄与Ru-N₄位于顶点，Pt和Rh因HER竞争过强被排除；

量化N₂吸附能、*N₂H自由能与极限电位之间的关系。

发表结果：发表于ACS Catalysis。客户实验验证：Mo-N₄ NH₃产率42 μg·h⁻¹·mg⁻¹，法拉第效率31%。

案例四：Ru-SAC酸性OER——配位确认+全路径计算

课题背景：客户合成了Ru单原子掺杂Co₃O₄酸性OER催化剂。审稿人质疑两点：Ru取代了Co的哪个晶体学位点？酸性OER的决速步是哪一步？

计算内容：

构建Ru取代四面体位与八面体位两种模型，计算形成能+XANES模拟，确认Ru优先占据八面体Co³⁺位；

OER四步质子-电子转移全跑NEB，每步虚频验证；

*O→*OOH为决速步（ΔG=1.58 eV），与实验Tafel斜率一致；

Ru掺杂后决速步能垒比纯Co₃O₄降低0.34 eV，d带中心分析归因于Ru 4d轨道上移。

发表结果：发表于Advanced Energy Materials。客户直接采用本方案的自由能台阶图作为论文TOC图。