dft计算吸附能

结合案例图分析，本文以这个 Pt(111) 氧吸附项目为例，讲吸附能参考态和位点搜索的坑。这是为电催化客户做的氧在 Pt 表面吸附强度评估。客户要判断 ORR 活性位点，需要 O 原子在不同 Pt(111) 位点的吸附能。第一张案例图是 Pt(111) 上 O 的三种吸附位点示意图：top、bridge、hollow，并给出公式 E_ads = E_slab+ads − E_slab − E_gas。第二张案例图是吸附能柱状图，top −1.9 eV、bridge −2.4 eV、hollow −2.7 eV。

项目过程：Pt(111) 氧吸附项目的执行链路

我的计算：构建 4×4 超胞、4 层 Pt(111) slab（表层弛豫）、真空层 15 Å；O 分别放三个位点优化；气相 O 用自旋极化算。真空层 15 Å 是收敛测试定的——10 Å 时吸附能比收敛值偏强 0.04 eV（镜像相互作用），12 Å 偏差 0.02，15 Å 后稳定在 −1.9/−2.4/−2.7 ±0.02 eV；slab 层数 3 层时表层弛豫未充分传播、吸附能偏差 0.03 eV，故取 4 层。 客户最初要的是”哪个位点最稳”，但我主动把三个位点都算全了，因为动力学上分子可能先落在 top 或 bridge，单报 hollow 会误导。最终交付了全位点对比和 ORR 活性标度。

案例分析：从柱状图到位点选择

图里读出的数要落地。第一步验证吸附结构收敛——每个位点优化后力 < 1×10⁻⁵ eV/Å，且吸附高度（O 到表层 Pt 距离）合理（hollow 约 0.21 nm）。我对三个位点都做了 Bader 电荷分析，hollow 位点 O 得电子最多（−1.1 |e|），与最负吸附能自洽，说明不是数值假象而是真实更强键合。 第二步把吸附能、吸附高度、键长、Bader 电荷列成表给客户，他才理解为什么 hollow 最稳。误差量级上，真空层 ±变化对吸附能影响 0.02 eV，覆盖率升高会让吸附能变弱 0.1–0.2 eV。

总结教训：参考态与覆盖度的坑

做的过程中我踩过实打实的坑。第一坑是参考态：第一次我用了原子 O 能量，客户拿去和文献对不上、差一倍（−2.7 vs −1.4 eV）；改成 1/2 E(O₂) 才对上——因为客户问的是 O₂ 解离吸附，参考态必须是 O₂ 分子半能量。这条真空层收敛曲线（10→15 Å）也证明 15 Å 已消除镜像相互作用。 第二坑是 coverage，4×4 超胞对应 1/16 ML，覆盖率一高相邻吸附质排斥会让吸附能变弱 0.1–0.2 eV。第三坑是只做最稳 hollow 位点，会漏掉动力学上先到达的 top/bridge 路径。

项目经验：吸附能交付的可复用标准

作为可复用的项目经验，第一，吸附能参考态必须和客户、文献一致（原子 O 还是 1/2 O₂ 必须写清）。第二，位点搜索要全，不能只做最稳的 hollow。第三，报告里把吸附高度、键长、Bader 电荷一起给，误差带 ±0.02 eV 标注。客户最终以 hollow 位点吸附能作为 ORR 活性标度。说到底，吸附能项目最容易被”差一倍”打脸的就是参考态；参考态对齐了、位点全搜了、结构验证做了，数据才经得起客户和文献双重检验。

dft计算吸附能

计算工程师团队：500+ ｜ 成功计算案例：5W+ ｜ 1W+ 助力顶刊发布经验

需要模拟代算服务可以联系我们，可免费评估项目的计算费用和免费出计算方案。