CO2吸附热计算：从等量吸附热到DFT结合能的两条路径

做CO2吸附材料的人，大概率在报告里见过”等量吸附热（Isosteric Heat of Adsorption, Q_st）”这个量。它描述的是：在特定吸附量下，CO2分子从气相进入吸附相时释放的热量。

这个量有两种获取途径：实验上从吸附等温线算，计算上用DFT算结合能。两条路各有优劣，合起来用价值最大。

实验路线：从吸附等温线到Q_st

实验上测的是不同温度下的CO2吸附等温线（通常是273K、298K、323K三温）。用克劳修斯-克拉佩龙方程（Clausius-Clapeyron）或者更准确的奥斯瓦尔德-弗罗因德利希方程（Ostwald-de Waele），可以算出来等量吸附热：

Q_st = -R × [∂(ln P) / ∂(1/T)]_q

实际操作里，用杜比宁-拉德克维奇（Dubinin-Astakhov, D-A）方程拟合等温线，再从拟合参数里提取Q_st，比直接数值微分稳定得多。

实验路线的优点是直接反映真实材料的吸附行为，包括孔隙填充、协同吸附、湿度影响等计算里很难完全考虑的因素。缺点是低温/高压的数据难测，Q_st的误差在低热区（<20 kJ/mol）比较大。

计算路线：DFT结合能

DFT计算CO2吸附热的思路很直接：

E_ads = E(CO2@material) – E(material) – E(CO2)

这个差值就是CO2分子和材料之间的结合能。正值表示放热（稳定吸附），负值表示吸热（不稳定）。

VASP里的具体操作：

第一步：构建吸附体系。在材料的晶体结构里（比如MOF的cif文件），找一个合理的吸附位点。CO2通常吸附在：

• 开放金属位点（比如Mg-MOF-74里的Mg²⁺）
• 碱性官能团附近（─NH2, ─OH）
• 窄孔隙的重叠区域（重叠吸附，overlap adsorption）

第二步：结构弛豫。对吸附体系做几何优化，INCAR关键设置：

ISIF = 2        （固定晶格，只弛豫离子位置——适用于MOF）
EDIFF = 1E-5
EDIFFG = -0.03   （力收敛标准，eV/Å）

如果研究的是柔性MOF（gate-opening行为），需要ISIF=3（同时弛豫晶格），但计算量会大很多。

第三步：计算结合能。分别算E(CO2@material)、E(material)、E(CO2)，注意要用相同的k点和截断能，让误差抵消。

BSSE修正：不能忽略的事

DFT计算的结合能，如果不做基组重叠误差（Basis Set Superposition Error, BSSE）修正，通常会系统性偏大10-30%。

MOF材料里，做BSSE修正的计算量很大（需要算ghost atoms的基组贡献）。实用策略是：用较小的基组（比如PAW PBE，截断能400 eV）做初步筛选，对Top候选材料再用大基组+BSSE修正做精细计算。

VASP里没有内置的BSSE修正工具（因为它是平面波基组，不是原子中心基组），但可以用”counterpoise correction”的思路，把CO2的假原子加到材料里算能量修正。

等量吸附热 vs DFT结合能：怎么对应

实验的Q_st和DFT的E_ads，物理上不是同一个量。Q_st包含熵的贡献（吸附过程ΔS<0，所以Q_st > |E_ads|），而且反映的是体相吸附（多个CO2分子共同吸附），而DFT通常算的是孤立吸附位点（单一CO2）。

对应方法：用DFT算不同吸附覆盖度（coverage）下的E_ads，看覆盖度从0到0.5 ML（单层等效）的变化。覆盖度升高，E_ads通常下降（位点间相互作用、竞争吸附）。

把不同覆盖度下的E_ads外推到零覆盖度，可以和实验Q_st的低覆盖度极限比较。

常见吸附位点类型

常见坑点

坑一：忽略范德华修正。PBE泛函完全不考虑范德华相互作用，算出来的物理吸附能接近零。必须用vdW修正（DFT-D3, vdW-DF, optB88-vdW）。

坑二：用错误的k点密度。MOF的晶胞通常很大（>1000 atoms/unit cell），只能用Gamma点。但CO2分子的气相计算也要用Gamma点，保持一致性。

坑三：忽略温度效应。DFT给的是0K的结合能，实验Q_st是室温（或实验温度）下的值。如果有条件，用分子动力学（MD）或者热力学积分，把有限温度效应加进去。

实际应用建议

如果目标是筛选MOF材料用于CO2捕集，推荐的工作流是：

1. 用DFT-D3（低计算成本）筛选100-200个候选MOF，看结合能是否在40-60 kJ/mol的”甜点”区间（吸附强但再生不太耗能）。
2. 对Top10候选，用更精确的泛函（比如HSE06）做单点能修正。
3. 用GCMC（巨正则蒙特卡洛）模拟不同压力下的吸附等温线，和实验比对。
4. 实验上测这Top10的CO2吸附等温线，验证计算结果。

CO2吸附热计算的价值，不在于算出和实验完全一致的Q_st（这需要考虑熵、覆盖度、湿度等诸多因素），而在于让你快速排除明显不行的材料，把实验精力集中在最有希望的候选上。