CP2K计算吸附能

一、项目背景与科学问题

二氧化碳的捕集与活化是碳中和技术链条中的关键环节。碱性氧化物表面（如 MgO）对 CO₂ 的吸附是实现低温碳捕集的重要途径，而吸附能的准确计算对筛选吸附剂、设计捕集工艺具有直接指导意义。与平面波程序（VASP）不同，CP2K 采用高斯与平面波混合基组（GPW/GAPW）方法，能够在保持精度的同时实现大规模体系的线性标度计算，特别适合表面吸附与分子体系。

本项目的科学问题是：CO₂ 在 MgO(100) 表面的吸附能是多少？不同交换关联泛函（PBE、PBE-D3、revPBE-D3、HSE06-D3）对吸附能的计算有何系统差异？吸附能可以分解为哪些物理贡献（相互作用能、形变能）？CP2K 的混合基组方案能否给出与平面波方法一致的可靠结果？

二、计算方法与理论背景

2.1 CP2K 混合基组方法

CP2K 使用 GPW 表示法：波函数用高斯型原子轨道（GTO）展开，而电子密度与 Hartree 势用平面波辅助基组表示。平面波网格用多网格切分，通过快速傅里叶变换计算静电项，从而实现 O(N) 或近线性标度。GAPW 扩展则允许全电子势与杂化泛函计算。对表面体系，该方法可结合周期性边界与 PBE-D3 处理范德华作用。

2.2 吸附能分解

总吸附能 E_ads 可分解为 E_ads = ΔE_int + ΔE_def,CO₂ + ΔE_def,surf，其中 ΔE_int 为冻结几何下吸附质与表面的相互作用能，ΔE_def 为吸附质/表面从自由态弛豫到吸附态的形变能。分解有助于判断吸附的主导贡献是化学成键还是物理色散。

2.3 色散校正

CO₂ 是弱极性线性分子，其与氧化物表面的吸附包含显著的范德华贡献。DFT-D3 方法在泛函能量上叠加原子对色散修正，能显著改善弱吸附体系的吸附能精度；revPBE 作为改良的 GGA 泛函，对分子吸附能也有系统修正。

三、计算设置与模型构建

• 软件与方法：CP2K，GPW/GAPW 混合基组，分子优化基组 DZVP-MOLOPT-SR，辅助平面波网格截断 400 Ry。
• 表面模型：MgO(100) 三原子层板模型，2×3 表面超胞，真空层 15 Å，底层固定。
• 泛函：PBE、PBE-D3、revPBE-D3、HSE06-D3 四组对比计算。
• k 点：表面超胞采用 2×2×1 采样（Γ 点为中心）。
• 收敛标准：电子自洽 10⁻⁶ au，几何优化力收敛 0.001 a.u.。
• 分析：对吸附结构进行能量分解，比较各泛函的吸附能与平衡几何。

四、结果分析与案例图解读

图 3 的左面板（a）给出四种泛函计算的 CO₂/MgO(100) 吸附能柱状图。PBE 给出 −0.44 eV，PBE-D3 提高到 −0.62 eV，revPBE-D3 为 −0.57 eV，HSE06-D3 为 −0.66 eV。PBE 到 PBE-D3 的吸附能增加约 0.18 eV，清晰展示了色散校正在弱吸附体系中的关键作用；revPBE 的修正方向与 D3 叠加后介于两者之间；杂化泛函 HSE06 给出最强的吸附能 −0.66 eV，反映精确交换对电子局域化的改善。整体结果与文献中 CO₂ 在碱性氧化物上的化学吸附量级（约 −0.5 至 −0.7 eV）一致。

图 3 的右面板（b）给出 PBE-D3 构型下吸附能的物理分解：相互作用能 ΔE_int = −0.66 eV 为吸附的主体贡献，说明 CO₂ 与 MgO 表面存在显著的化学相互作用（CO₂ 的碳与表面氧、氧与表面镁成键）；CO₂ 形变能 ΔE_def,CO₂ = +0.03 eV，表面形变能 ΔE_def,surf = +0.01 eV，两者均为很小的正值，表明吸附引起的几何畸变很小，线性 CO₂ 分子仅略微弯曲。总吸附能 = −0.66 + 0.03 + 0.01 = −0.62 eV，与直接计算一致。

两面板结合，左面板比较了泛函方案对吸附能的系统影响，右面板揭示了吸附能的物理构成，共同说明 CP2K 混合基组方法能够可靠地描述 CO₂ 在氧化物表面的弱化学吸附，且色散校正是不可省略的关键环节。

五、结论与展望

本项目用 CP2K 的 GPW/GAPW 方法计算了 CO₂ 在 MgO(100) 表面的吸附能，PBE-D3 结果为 −0.62 eV，能量分解表明相互作用能主导、形变能微小。案例图以泛函对比柱状图与能量分解柱状图双面板呈现。

后续可拓展：(1) 考察 CO₂ 覆盖度与表面缺陷（氧空位）对吸附能的影响；(2) 结合 AIMD 研究 CO₂ 在表面上的动态吸附-脱附；(3) 用同一方法筛选取向不同的碱土金属氧化物吸附剂。