实验等温线如何对接计算？多孔材料气体/蒸气吸附的读图与模拟解读

说明：本文华算科技主要介绍多孔材料气体与蒸气吸附等温线的读图方法，以及图形特征与分子模拟、第一性原理计算之间的关系。

吸附等温线的定义与构造

吸附等温线把温度固定，让气相压力逐点改变，记录样品达到平衡后的吸附量。横轴可以是压力 p，也可以是相对压力 p/p₀；纵轴常用 mmol·g⁻¹ 或 cm³·g⁻¹。每个数据点都对应一个气体化学势和一个平均吸附量，横轴记录热力学状态，时间未出现在坐标中。对水蒸气，p/p₀ 与相对湿度对应；对超临界气体，通常使用绝对压力。

实验中的质量基准、样品活化状态和吸附量定义必须统一。干样质量与含水样质量不同；过量吸附量扣除了气相在孔体积内所占的量，与绝对吸附量采用不同定义。曲线高低受质量基准和吸附量换算影响，称量差异也会造成表观容量差异。吸附支由低压升到高压，脱附支沿反方向测量，两条路径在同一压力下未必重合。对强吸附气体，每级压力的平衡时间会影响慢扩散条件下的读数。

IRMOF-74-V 的氩气曲线把宏观吸附量与孔内占位连在一起：低压时分子靠近孔壁，压力增加后出现第二层，随后孔腔填充。模拟曲线来自巨正则蒙特卡洛对各个化学势下孔内占位的平均统计。图里的平台只表示所测压力范围内增量变小；若横轴继续延伸，材料外表面凝聚或其他过程仍可能改变读数。曲线的转折位置要和温度、气体种类一起报告；87 K 氩气与室温 CO₂ 的拐点压力不能用同一标尺评价孔壁亲和力。

等温线的吸附量通常按样品质量归一，但质量单位相同也不保证样品处于同一初始状态。残余溶剂可能占据孔位；活化温度过高则可能改变柔性框架。比较微孔材料的低压吸附量时，活化后孔道是否开放与吸附物的动力学直径都影响实际可进入空间。升压和降压过程的停留时间也应一致记录。

低压上升、拐点与平台分别怎样读

低压区常用 n≈K_Hp 描述，K_H 是亨利系数，反映稀疏分子进入材料时的初始响应。斜率大可来自较强相互作用，也受可进入位点数影响；它本身不是某一位点的吸附能。CO₂ 在 MOF 上的低压数据尤其需要检查测量压力是否仍在线性区，高亲和材料往往很早偏离亨利近似。

曲线出现膝部，可能对应高能位点逐渐占满，也可对应微孔填充开始转缓。单层 Langmuir 形式能描述某些简单表面，但多孔材料常有孔壁、电荷位点和不同孔径的贡献；拟合得到的“单层容量”与实测孔容采用不同定义。若低压骤升后很快趋平，测试气体的可进入孔径与独立孔径测量能约束填孔解释。

高压吸附量还取决于所用气体、测试温度和最大压力。两个材料在低压段交叉、高压段反向排序并不矛盾：前者可能有更强的初始结合，后者可能有更多可容纳空间。对于储气，工作压力范围内的可交付吸附量 Δn＝n(p_高)−n(p_低)比单一高压终点容量贴近循环工况。

Ⅰ型曲线的低压骤升与微孔填充相符，Ⅳ型曲线的高压阶跃常伴随介孔毛细凝聚。孔径仍需由气体分子尺寸、吸附温度和独立孔结构测量限定。氮气在 77 K 扩散受阻的超微孔可能呈现偏低的表观吸附量，换用氩气或较高温度的 CO₂，低压段的差异才可能显露。不同探针气体的曲线不宜直接按平台高度排序孔体积。

NU-1500-Cr 在较低湿度下缓慢吸水，约 33%–37% 相对湿度出现明显阶跃；起始水分子与框架位点结合，后续水分子形成氢键网络。该位置对应水分子从局部结合和孔道延伸转向主孔填充，阶跃压力表示集体填孔条件；第一批水分子已在更低湿度进入材料。

NU-1500-Cr 的水吸附焓在阶跃附近由较负值向较弱放热方向变化，水分子熵在这一湿度区间明显下降。主孔填充增加了水—水接触，也限制了分子的可访问构型；焓和熵曲线与吸附阶跃发生在相近湿度，支持集体填孔解释。单条等温线只有吸附量信息，分子占位还要由模拟密度或原位谱学辨认。

阶跃与滞后环的不同物理来源

介孔中的毛细凝聚、柔性框架的开孔转变及吸附物之间的协同作用都可能产生阶跃。NU-1500-Cr 的模拟密度显示低湿度水分布在连接主孔的孔道附近；水链增长后，主孔逐步填满。低湿度的位点占据和高湿度的分子间氢键对应不同的吸附阶段。

在介孔体系中，凝聚和蒸发的亚稳态、孔阻塞及空化会使吸附支与脱附支错开；柔性 MOF 在升压和降压时还可能经过不同的晶格状态。墨水瓶孔只是滞后环的一种来源，环面积受孔径分布和测量条件共同影响。原位衍射中的晶格变化与变温曲线可以区分孔内相变和框架形变。

狭窄孔颈会推迟介孔孔腔中的氮气排空，脱附支因此可能在较低压力突然下降；单个独立圆柱孔的计算模型无法表示这种网络效应。以脱附支反演孔径时，孔阻塞与空化会改变表观峰位，吸附支和脱附支给出的孔径分布需按各自模型假设解释。

Fe(bdp) 与 Co(bdp) 的 CH₄ 曲线出现吸脱附转变压力错位，框架响应可能参与其中。图中实心点和空心点分别对应吸附与脱附；曲线还随温度改变。测量若尚未达到充分平衡，动力学迟滞会增大两支曲线间距；压力点的平衡条件决定环形是否反映准平衡路径。相似环形可由孔网络效应产生，框架结构随气压变化的证据有助于区分来源。

温度变化与计算结果怎样相互约束

同一材料的多温度等温线能区分容量变化和吸附热效应。在固定吸附量 n 下，沿不同温度曲线找到相应平衡压力，可用 Q_st＝−R[∂lnp/∂(1/T)]_n 估计等量吸附热。斜率必须沿固定吸附量取得；阶跃附近若存在明显滞后或相变，局部斜率会混入相变与路径效应。

Fe(bdp) 和 Co(bdp) 的转变压力随温度移动，吸附与脱附比较要限定热力学状态。两支的转变宽度也在变化，工作压力窗口随温度改变。对于水蒸气，横轴若用相对湿度，饱和蒸气压 p₀ 也随温度变化；相同相对湿度对应的绝对压力并不相同。

巨正则蒙特卡洛模拟在给定温度和化学势下统计孔内分子数，适合与实验等温线逐点对照；第一性原理 DFT 常用于检验局部吸附构型、电荷分布与相互作用能。整条等温线还包含覆盖度增加后的分子间作用和构型数变化，单个 DFT 吸附能只描述某个给定状态。孔径分布反演中常见的“DFT 模型”属于经典密度泛函，名称相同，求解对象不同。

模拟与实验可比较低压斜率、阶跃压力和平台容量，但三处偏差对应的来源未必相同。低压斜率对框架电荷与气体—孔壁相互作用敏感；阶跃位置还受分子间作用和框架柔性影响；平台容量则受可进入孔体积与缺陷影响。固定框架 GCMC 的阶跃偏移可以提示势函数或结构模型需调整，不能直接用比例因子平移整条曲线。

IRMOF-74-V 的模拟吸附量出现跃迁，实验和模拟的晶胞参数也随压力改变。晶胞参数随氩气压力的变化为框架响应提供了独立信号。若模拟阶跃位置偏移，低压段偏差主要涉及气体—框架势函数，跃迁区还涉及框架柔性，平台区则需检查可进入孔体积、晶体缺陷和实验活化状态。