CO2吸附热怎么计算：DFT框架下从slab建模到吸附能提取的完整路径

CO2吸附热怎么计算这个问题，在碳捕集领域几乎每个做第一性原理的人都反复遇到过。吸附热的大小决定了材料对CO2的选择性和再生能耗——吸得太弱，选择性不够；吸得太强，脱附需要的温度和能耗急剧升高。理想窗口在-0.4到-0.8 eV之间（对应约-39到-77 kJ/mol），这个范围既能保证一定选择性，又不过于耗能。但用DFT把吸附热算到这个精度，远不是”做一次弛豫减两个能量值”的事。

表面模型的收敛性

团队做过一个MgO(100)负载过渡金属的CO2吸附项目，计算吸附热的第一步就是确定slab模型的厚度和真空层。MgO(100)的slab层数从3层开始测试，总能量随层数增加的变化：3层→4层吸附能变化0.12 eV，4层→5层变化0.04 eV，5层→6层变化0.02 eV。这说明4层是能量收敛的边界，但考虑到表面重构和电荷密度分布，团队最终用了5层——底部2层固定、顶部3层弛豫。

真空层厚度的选择直接影响CO2分子和上方周期镜像之间的虚假静电作用。CO2分子吸附在表面后，分子高度约3.2 Å（C到表面的距离），加上分子本身~2.3 Å的O-C-O长度，真空层至少需要15 Å才能让相邻slab之间的CO2镜像间距足够远。团队测试了12→15→18 Å的真空层，吸附能在15 Å以上才收敛到0.02 eV以内，最终选用16 Å。

范德华修正的取舍

CO2是线性非极性分子，但表面吸附时分子和载体之间有显著的色散作用——纯PBE泛函给出的吸附热通常比实验值弱0.2-0.4 eV。CO2吸附热怎么计算的第一个关键抉择就在这里：加不加范德华修正、加哪种。

DFT-D3（Grimme的半经验色散修正）是最常用的选择，但零阻尼D3在短程有双计数色散的风险。团队比较了三种方案：PBE无修正、PBE-D3(BJ阻尼)、optB88-vdW泛函。PBE给出的MgO(100)上CO2吸附热是-0.31 eV，PBE-D3(BJ)给出-0.52 eV，optB88-vdW给出-0.49 eV。实验报道的MgO(100)上CO2吸附热在-0.45到-0.55 eV之间——D3(BJ)的结果最接近实验中值。

差距不会说谎：纯PBE对CO2这类弱吸附体系的吸附热系统性地低了约40%。团队后续所有CO2吸附计算都采用PBE-D3(BJ)方案，并在INCAR里明确设置IVDW=11（零阻尼D3用IVDW=10，BJ阻尼用IVDW=11，这是VASP 5.4.4以后的标准）。

吸附构型的采样策略

CO2在氧化镁表面可以有多种吸附方式：碳朝下（C-Mg）、氧朝下（O-Mg）、平行于表面（flat-lying）。以碳朝下为起点，C-Mg距离从2.5到4.0 Å扫描得出最优距离约2.85 Å，吸附时CO2分子从线性180°弯折到约167°——氧原子被表面静电场推开。

但真正麻烦的是，optB88-vdW泛函下的势能面非常平坦：在C-Mg距离2.8-3.2 Å范围内，能量变化不到0.05 eV。这种平坦的势能面意味着吸附构型对外界条件（温度、覆盖度）非常敏感。团队在不同覆盖度下（1×1、2×2、3×3超胞，分别对应1、1/4、1/9 ML覆盖度）做了测试，发现高覆盖度（1 ML）下CO2-CO2侧向排斥使吸附热从-0.52 eV减弱到-0.41 eV——这个差异在选择性预测中足以改变结论。

振动态修正

吸附热在0 K计算的是电子能之差，但室温下的吸附热需要考虑零点能和热焓修正。频率分析是必须的——CO2分子吸附后，平动和转动自由度转化为振动模式，零点能从~0.12 eV（自由CO₂）变化到~0.15 eV（吸附态）。按吸附能Eads = E(slab+CO2) – E(slab) – E(CO2) + ΔZPE计算，ΔZPE约0.03 eV的修正量级虽然不大，但和吸附能误差在同一数量级上，不能忽略。

这个项目最终给出的MgO(100)上CO2吸附热是-0.54 eV（含ZPE修正），与程序升温脱附实验值-0.52±0.05 eV高度一致。CO2吸附热怎么计算不是一个能量减法问题，而是从slab建模到构型采样到振动态修正的一条完整链条。对能源催化项目中吸附热力学精度有严格需求的同行，可访问站点了解更多实操方案。

更多内容请访问 https://www.keyanxueshu.com/