气体吸附焓计算这件事，在MOF、沸石、多孔碳材料的研究中几乎是绕不开的。但麻烦的地方在于，”吸附焓”这个词本身涵盖了不同层面的含义，不同来源的数据拿来直接比较，很容易得出误导性的结论。
等量吸附热(isosteric heat of adsorption, Q_st)是实验端最常用的表征量，通过不同温度下的吸附等温线用Clausius-Clapeyron方程导出。它描述的是覆盖度依赖的吸附强度变化，能反映材料吸附位点的分布情况。而理论计算端拿到的，通常是0 K下的结合能，两者之间差着有限温度修正和零点振动能贡献。
DFT计算气体吸附焓的基本框架
气体吸附焓计算的DFT路线，核心与吸附能计算类似：E_binding = E(宿主+气体) – E(宿主) – E(气体)。但气体分子的参考态需要格外注意——以CO2为例，孤立分子的能量需要在足够大的超胞里单独计算，并且保证与吸附体系一致的截断能和自旋设置。
对于MOF和多孔材料体系，van der Waals色散修正几乎是标配。这类体系的气体吸附以物理吸附为主，PBE泛函在不加色散修正的情况下会严重低估结合能，与实验值的偏差常常超过50%。DFT-D3或DFT-D3(BJ)是目前使用最广泛的修正方案，在CO2、CH4、H2O等小分子体系上表现稳定。
覆盖度效应与多位点吸附的处理
真实材料中的气体吸附焓并不是一个固定值，它随覆盖度变化。DFT单点计算给出的是理想孤立吸附位的结合能，对应实验低覆盖度极限下的Q_st初始值。
如果需要考察不同覆盖度下的吸附焓变化，可以在不同负载量下分别计算结合能，模拟逐步填充的过程。MOF体系中常见的情况是：强吸附位点先被占据，结合能随覆盖度单调下降;而某些具有协同效应的体系，低覆盖度时结合能反而偏弱，负载增加后因邻近分子间作用而有所增强。这两种行为在气体吸附焓计算结果中会留下截然不同的特征。
有限温度修正与零点能的贡献
DFT结合能对应0 K静态晶格的极限，要与实验等量吸附热比较，需要引入有限温度修正。标准做法是通过谐振近似计算振动贡献，将零点振动能(ZPE)和热容修正加入到能量中。
对于轻质分子(H2、CH4)，ZPE修正的量级在10-40 meV范围，不可忽略。CO2和水分子因质量较大，ZPE贡献相对小，但仍建议在高精度计算中纳入。气体吸附焓计算的精度最终取决于这些修正项是否被完整考虑，而不仅仅是DFT基态能量的高低。
结果验证与实验数据对比策略
气体吸附焓计算的结果拿去与实验对比时，需要确认对比的是同一物理量。实验给出的等量吸附热是焓变(包含PV功)，而DFT直接算出的是内能变化。在气相吸附场景下两者差一个RT项(约2.5 kJ/mol在室温)，通常可以接受，但在高压吸附体系中这个修正不能忽略。
多孔材料的气体吸附焓计算如果和实验偏差超过20%，先检查色散修正是否开启，再检查气体分子参考态的计算是否一致，这两个环节覆盖了大多数常见误差来源。
