手机版

孔道内吸附量计算

发布时间：2026-08-31 来源：科研学术网

小中大

字号：

一、项目背景与科学问题

多孔材料（金属有机框架 MOF、沸石、共价有机框架 COF）具有高比表面积与可调孔道结构，在气体吸附与分离、碳捕集、储氢与催化领域具有广阔应用。吸附量（uptake）是评价多孔材料吸附性能的核心指标，它取决于孔道尺寸、比表面积、孔壁化学环境与吸附质性质。准确预测不同温度、压力下的吸附等温线，对吸附工艺设计与材料筛选至关重要。

巨正则蒙特卡洛（GCMC）模拟是计算多孔材料中气体吸附量的标准方法：在固定温度、压力（化学势）与孔道体积下，通过随机插入、删除、平移吸附分子来采样吸附平衡态，获得平衡吸附量与吸附热等热力学量。本项目以 CO₂ 在 MOF-5 中的吸附为对象，用 GCMC 计算不同温度下的吸附等温线与等量吸附热。

二、计算方法与理论背景

2.1 GCMC 模拟原理

GCMC 在巨正则系综（μVT）中采样：体系与虚拟的无限大气体储库处于化学势平衡，通过三种随机尝试——插入、删除与平移吸附分子——满足细致平衡，达到吸附平衡。压力通过理想气体化学势 μ = μ⁰ + RT·ln(P/P⁰) 关联。模拟输出的平均粒子数即平衡吸附量，可换算为质量吸附量（mmol/g）与体积吸附量。

2.2 吸附等温线类型

物理吸附等温线通常为 I 型（Langmuir 型，微孔饱和）或 I/II 复合型。对 CO₂ 在 MOF-5 的吸附，低压区吸附量随压力近似线性上升（Henry 区），随后趋于饱和平台。温度升高使吸附量下降（吸附为放热过程），因此低温等温线位于高温等温线上方。

2.3 等量吸附热

等量吸附热 Q_st 描述吸附的放热强度，由不同温度等温线的克劳修斯-克拉佩龙关系求得：Q_st = −R[∂lnP/∂(1/T)]_n。Q_st 反映吸附质-孔壁相互作用的强弱，对区分物理吸附（20–50 kJ/mol）与化学吸附（>60 kJ/mol）具有指示意义。对 CO₂ 在 MOF-5 中，Q_st 约 17–25 kJ/mol，属于典型物理吸附。

三、计算设置与模型构建

软件：RASPA / 自编 GCMC 程序；力场：UFF/通用力场描述 MOF 骨架，CO₂ 采用 EPM2 或 TraPPE 模型。
结构：MOF-5 实验单胞（含客体水去除），框架原子固定。
模拟：温度 273/298/323 K；压力 0–5 bar；每个状态点 5×10⁶ 步平衡 + 5×10⁶ 步采样。
截断：12 Å 以上，长程静电用 Ewald 求和。
输出：绝对/过量吸附量、等量吸附热、吸附位点分布。

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

图 7 的左面板给出 CO₂ 在 MOF-5 中 273 K（蓝色）、298 K（红色）与 323 K（绿色）三条过量吸附等温线。以 298 K 为例，吸附量随压力单调上升：0.1 bar 时约 5 mmol/g，1 bar 时约 19 mmol/g，5 bar 时约 38 mmol/g。等温线在低压区斜率较大（Henry 区），高压区趋于平缓，呈现典型的 I 型微孔吸附特征。温度效应显著：同一压力下，273 K 的吸附量明显高于 323 K，例如 1 bar 时 273 K 约 24 mmol/g、298 K 约 19 mmol/g、323 K 约 15 mmol/g，符合物理吸附放热的普遍规律。该等温线族为 MOF-5 的 CO₂ 捕集工艺提供了直接设计数据——低温高压有利于提高吸附容量，而变温吸附（TSA）再生正是利用这一温度依赖性。

图 7 的右面板给出等量吸附热 Q_st 随吸附量的变化（橙色曲线）。Q_st 在低覆盖度时约为 −25 kJ/mol，随吸附量增加先小幅波动后趋缓，整体在 −20 至 −26 kJ/mol 区间。该量级远低于化学吸附阈值（60 kJ/mol），确证 CO₂ 在 MOF-5 中为物理吸附，脱附容易、再生能耗低。Q_st 随覆盖度变化平缓，说明孔道内吸附位点相对均匀，无明显强吸附位优先占据现象。这一信息对吸附工艺的经济性评估至关重要——物理吸附的再生只需小幅升温或降压即可实现。

两面板结合，左面板给出吸附容量与温度压力的定量关系，右面板揭示吸附强度与吸附机制，共同完成 MOF-5 对 CO₂ 吸附性能的系统评估，为碳捕集吸附剂筛选与工艺参数优化提供完整数据。