多孔材料液体吸附仿真：GCMC模拟MOF-5对CO₂/N₂混合气选择性吸附

用GCMC（巨正则蒙特卡洛）模拟多孔材料吸附，工程上不算新东西，但要算准选择性这个指标，参数链上的每一环都值得较真。这个项目模拟MOF-5分离电厂烟气中的CO₂/N₂（比例15:85，温度298 K），目标是得到不同压力下的CO₂吸附量和CO₂/N₂选择性，给后续工艺设计提供理论吸附等温线。

力场：通用力场和专门力场差了30%

MOF-5是经典的Zn₄O(BDC)₃框架结构。初期用UFF通用力场描述框架原子，CO₂和N₂用TraPPE力场，Lennard-Jones交叉参数用Lorentz-Berthelot组合规则。跑完1 bar下的GCMC后，CO₂吸附量为2.8 mmol/g——但实验值是4.0 mmol/g，偏差30%，方向是偏低。

问题出在UFF对Zn-O簇的范德华参数过于”硬”。实际上MOF-5的Zn₄O节点中，O原子极化率被Zn的电荷效应增强，导致CO₂与框架节点之间的色散作用被UFF低估。切换到DREIDING力场+部分原子电荷修正（DDEC电荷），CO₂吸附量提升到3.9 mmol/g，与实验值的差距缩到3%。

这个偏差的启示是：MOF吸附模拟不是跑通了就完了。通用力场在描述刚性金属有机框架时，框架-吸附质作用能误差可能大到你没法用这个结果做工艺设计。

逸度-压力换算：高压下忽略这一条就错了

GCMC的输入是逸度不是压力。在1 bar以下，CO₂和N₂的逸度系数接近1，逸度≈压力，直接用压力做输入问题不大。但项目要算到20 bar——20 bar下CO₂的逸度系数约0.75（偏离理想气体25%），如果还直接用压力代替逸度，相当于人为抬高了气相化学势，导致GCMC计算的吸附量偏高15-20%。

项目用Peng-Robinson状态方程把每个压力点换算成对应逸度再输入RASPA，高压段的数据才和文献实验值对齐。这个细节在教程里很少被单独强调，但它对高压吸附模拟的结果影响量级是可观的。

超额吸附量和绝对吸附量的区别

另一个容易混淆的点是：GCMC输出的是绝对吸附量（孔内吸附质总量），而实验测量的是超额吸附量（减去纯气相贡献）。在高压下，孔内气相密度不可忽略——20 bar下CO₂本体密度已达~0.04 g/mL，如果不扣除这一部分，绝对吸附量会被高估。

项目按超额吸附量公式换算后，20 bar下CO₂的绝对吸附量为8.2 mmol/g，超额吸附量为6.9 mmol/g——差了16%。选择性（CO₂/N₂）在1 bar下达18.5，在20 bar下降到11.2，因为高压下竞争吸附效应减弱。这个趋势和文献中的理想吸附溶液理论（IAST）预测一致。

整套流程走完之后最深的体会是：GCMC模拟多孔材料吸附，算力需求不高（一个压力点几分钟），但参数链上的纠正步骤是决定性的。逸度换算、力场匹配和超额吸附修正——三道关卡跳过任何一个，结论都在另一条路上。

参考文献：Yazaydın et al., J. Am. Chem. Soc., 2009, 131, 18198-18199；Walton et al., J. Am. Chem. Soc., 2008, 130, 360-369.