Jan T. Duda† 和 Janina Milewska-Duda*,‡ † 克拉科夫AGH科技大学 自动化学院 ‡ 克拉科夫AGH科技大学 燃料与能源学院 该研究于2002年8月20日在线发表于 Langmuir 期刊（第18卷，第20期，第7503-7514页）。
该研究属于物理化学与材料科学领域，主要关注微孔材料中的气体吸附现象。在表面科学和多孔材料表征中，低温气体吸附法是测量材料比表面积和孔结构的关键技术。经典的BET（Brunauer-Emmett-Teller）理论作为分析多层吸附的基础，被广泛用于通过低压段（相对压力通常小于0.3）的吸附等温线来计算材料的比表面积。然而，当吸附质为近似球形的小分子时，在较高的相对压力下，特别是在微孔材料中，BET方程的适用性变得很差。这是因为BET模型假设吸附质会形成不受限制的、堆叠状（stacklike）的分子团簇（clusters），而忽略了在真实微孔中存在的几何约束（geometrical constraints）以及可能形成的比BET堆叠更紧密的团簇结构，即非BET构型效应（non-BET configurational effects）。
为此，作者在过去的研究中进行了开创性工作，通过结合溶液热力学理论对BET方法进行了推广，提出了能够处理吸附质-吸附质相互作用和几何约束的通用多层吸附数学模型，并推导出了用于近似计算的LBET和LGBET方程。本文旨在进一步发展这一理论框架，其核心目标有四个：第一，将此前建立的复杂通用模型转化为一个更便于分析和计算的递推形式（recurrent form），并从形式化角度深入探究孔几何形状与吸附等温线形态之间的内在联系；第二，为更精确地描述受限多层吸附，提出一种全新的近似公式——LCBET方程；第三，尝试处理在最高相对压力处（即接近凝聚点）的吸附特性；第四，通过将模型应用于分析合成活性炭的氮气吸附实验数据，展示其在解释高比压吸附行为和定量表征孔结构方面的潜力。
研究工作的第一步是对此前建立的通用多层吸附模型进行重写。该模型的核心思想是将吸附空间中的分子视为在不同层上形成的团簇。一个团簇从一个初级吸附位点（primary site）开始（第0层），允许一个分子提供多个上层吸附位点（即允许非BET堆叠），从而可以描述更紧凑的团簇。模型定义了最大层数k来描述该位点的容量，并通过参数β_kn来表征第n-1层到第n层的平均孔宽，β_kn = 1表示狭窄的通道状或狭缝状孔，β_kn > 1表示更紧凑或树枝状孔。通过将能量项ΔH和构型熵项ΔS的基本热力学方程（公式1-3）与平衡条件结合，研究者成功推导出了一套描述第k类初级位点上第n层吸附行为的严格方程组（公式4-6）。这套新推导出的递归方程组（公式7-9）构成了整个研究的理论基础框架。它简洁地表达了局部吸附等温线，即在给定参数q_kn（摩尔吸附自由能）、β_kn和初级位点数量m_hAk后，可通过逐层求解覆盖率θ_kn，进而计算出整个团簇的吸附量。
在得到递推公式后，作者对模型进行了深入的形式化理论分析。他们引入了辅助序列s_i来研究在具有常数参数b_c和β（β≥1）的无限大孔中，上层吸附覆盖率的行为。分析揭示了一个关键的理论现象——存在一个“凝聚压力”（condensation pressure）π_c。当相对压力π小于或等于一个最大值π_m（由π_m = β^(-β)(β-1)^(β-1)定义）时，序列s_i会收敛到一个稳定值s_s，这意味着覆盖层会趋向一个稳态覆盖率θ∞ = π_∞*；而当π > π_m时，序列会发散，覆盖率趋向于1，导致孔隙被完全充满。这一发现从理论层面洞察了受限多孔结构中多层吸附的两种可能行为，即平稳增长和突跃式填充。
基于此分析，作者提出了模型在实际应用时应将吸附系统分为三个部分（图1）：可严格处理的小团簇（容量k有限）、大团簇中的底层（容量k很大，但底层1到l层行为不同），以及大团簇中的上层（采用统一平均值θ）。为描述大团簇容量的分布，他们采用了指数分布m_hAk = m_hA(1-α)α^(k-k-1)，这一处理方式极大地简化了模型在真实多孔材料上的应用。
这是本研究最关键的方法论创新。此前提出的LGBET模型将上层覆盖率θ直接近似为稳态值π_∞，虽然提供了一个有价值的参考点，但在高压区的预测偏高，且严格限制了其仅适用于π < π_m的范围。为克服这些局限，作者提出了一个新的加权平均公式来计算上层平均覆盖率θ：θ = π(1 + w_hπ) / (1 + π)，其中π是通过方程π(1-θ)^(β-1) = π/b_c定义的修正压力，w_h是依赖于分布参数α的权重因子。这个公式通过巧妙结合朗缪尔（Langmuir）等温式和广义亨利（Henry）定律，更真实地模拟了有限尺寸团簇上层覆盖率的中间状态。该模型（即LCBET模型）将适用范围扩展到了π_m以上，直至由新方程f(π_max*) = 0确定的新上限π_max，实现了对由递推模型生成的理论等温线在整个可行压力范围内的高精度近似（图5）。
为验证LCBET公式的优越性，作者使用递推严格模型生成了一系列具有代表性的理论等温线（参数化：m_hA=1，q_a/RT=-0.5，不同α和β），然后用“三点法”识别程序分别用LGBET和LCBET去拟合这些等温线。拟合区域设定为π≤π_c。对比结果（表1）令人信服地表明，尽管两种模型都能在拟合区内实现极好的拟合效果，但LCBET模型反演出的参数（尤其是α、β和b_c）远比LGBET模型更接近真实值，由此验证了其参数评价的可靠性。该识别程序的具体步骤是：先用低压数据（π<0.4）拟合朗缪尔成分，得到m_hA和b_a；然后在三个选定的高压点（π_r, π_b, π_β）上调整α, b_c, β以实现最优拟合。
最后，研究使用LCBET和LGBET模型对Choma和Jaroniec发表的两个合成活性炭（Ambersorb 563和572）在77.4K下的氮气吸附实验数据进行了分析。分析流程分为三步：首先，利用极低压数据（6×10^-3 < π < 8×10^-2）的朗缪尔线性图确定表面参数，这一操作将表面异质性效应排除在主分析之外；然后，分别用LBET、LGBET的两种模式、以及LCBET进行等温线全范围拟合，所比较的参数包括m_hA、q_a/RT、α、β、b_c以及累积参数αβ/b_c和α/b_c；最后，为了解释在高压段（π > 0.85）吸附量的急剧增长，引入了一种可变b_kn的修正机制，即当层覆盖率θ超过某个临界值θ_c时，假定由于分子间接触加剧，吸附能发生下降，导致参数b_kn从b_c线性下降至1。
形式化分析清晰地揭示了孔几何与等温线形状的关键联系。当β>1时，存在一个凝聚压力π_c，在此压力下，具有无限容量的理想孔会发生等温线的阶跃式垂直上升（图4），标志着由构型效应引起的孔隙快速填充。对于实际有限孔，等温线在π_c附近是连续的，但其斜率变化快慢由孔容大小决定。这部分理论阐明了从低压到高压吸附行为的完整图像。
与严格递推模型的对比研究（图5）证明了LCBET模型的卓越表现。在π≤π_c的范围内，LCBET与严格解几乎完全重合，其近似精度显著高于LGBET。参数识别结果（表1）定量地显示了LCBET的优越性：例如，对于一个β=1.9、α=0.75的严格模型，LGBET拟合得到的α_e值严重低估为0.607，而LCBET得到的0.685则接近真实；更重要的是，决定高压段等温线形状的累积参数αβ/b_c和α/b_c，LCBET的评估误差远小于LGBET。这表明LCBET能更好地区分几何约束（α和β）与能量约束（b_c）对吸附的独立影响。
对活性炭数据的分析（图7, 8, 表2）证明： - 模型的适用性：任何具有恒定b_c的模型都无法解释在π>0.85观察到的吸附量急剧上升现象，这确认了存在由高覆盖率引发的吸附能量变化过程。 - 吸附机理的揭示：通过引入随覆盖率变化（b_kn从b_c降至1）的修正机制（LCBET*），模型获得了与实验数据的出色吻合。计算出的临界覆盖率θ_c约为0.26，结合较低的β值（约1.06-1.07），提示该急剧上升主要源于大孔中形成了沿不规则孔壁铺展的扁平团簇，一旦覆盖率增加导致的分子间相互作用克服了初始的高能量壁垒（表现为b_c异常高），吸附量便迅速上升。 - 孔结构定量信息：LCBET模型给出了更为可靠的孔几何参数。Ambersorb 563的数据可以被指数分布孔模型接近完美地拟合，而Ambersorb 572则显示出对此简单分布的偏离，需要额外增加大孔（k=33）的孔容才能实现对高压点的精确拟合。这表明两个材料虽然都含有宽孔径谱系，但约束状态和孔容分布不尽相同。
该研究成功地将一个复杂的热力学多层吸附模型转化为便于分析的递推形式，并在此基础上深入阐释了从低压到高压的吸附机理。其科学价值主要在于：首次清晰地从理论上定义了“凝聚压力”π_c的概念，并将其与孔几何参数（β）联系起来；揭示了由于构型效应和非堆叠团簇形成导致的非BET行为的内在规律。
更重要的是其应用价值：新提出的LCBET公式是对现有LGBET模型的重大改进，它填补了LGBET在π接近和超过π_m时的应用空白，能够在更宽的压力范围内提供可靠的孔隙形状和尺寸信息。研究指出，对高压段等温线形状影响的关键是累积参数αβ/b_c，而LCBET模型由于其上层覆盖率公式中包含了更丰富的几何效应，使得对单个参数（α, β, b_c）的区分成为可能，从而让孔几何的识别更加可靠。这种综合分析策略意味着，通过将LBET、LGBET和LCBET模型交叉应用，并结合高压吸附数据和理论凝聚压力，研究人员能够获得比传统BET方法深刻得多的孔结构表征信息。
本研究的核心亮点包括： 1. 理论创新：全面形式化地分析了服从指数分布的多层吸附模型，首次清晰地阐明了吸附覆盖率的收敛与发散条件，并定义了“凝聚压力”π_c，为理解高压吸附行为提供了关键理论工具。 2. 方法突破：创造性地提出了LCBET近似公式，通过一个巧妙的加权平均方案，将朗缪尔吸附与受限BET吸附结合起来，在保持计算简便性的前提下，大幅提升了对真实多孔体系在整个可行高压范围的描述精度和参数反演可靠性。 3. 深刻的物理洞察：通过对两种活性炭的分析，不仅定量评估了其孔表面特征，还揭示了高压段吸附急剧上升的潜在物理机制——即由覆盖率驱动的高层吸附能量下降，并成功用参数b_c的动态化处理对这一现象进行了统一描述，区分了材料中占主导地位的孔形状类型。