本文以UiO-66为代表性材料,介绍采用氮气孔隙分析法对金属有机框架材料进行表征的方法。
本文以UiO-66为代表性材料,介绍采用氮气孔隙分析法对金属有机框架材料进行表征的方法。
金属有机框架(MOF)的比表面积和孔体积可为其结构特征及潜在应用提供重要信息。这两个参数通常通过氮气吸附实验数据测定,用于此类测量的商用仪器也广泛可用。这些仪器能够自动计算结构参数,但理解如何选择输入数据以及何时适用特定计算方法对目标MOF样品至关重要。本文概述了采用Brunauer-Emmett-Teller(BET)方法和Barrett-Joyner-Halenda(BJH)方法分别计算比表面积和孔体积的具体步骤,并以典型的MOF材料UiO-66为例进行计算演示。尽管这些方法普遍适用于MOF材料,但要获得准确的计算结果,除了需进行恰当的样品前处理外,样品本身及吸附数据还必须满足一定的条件。本文还讨论了这些方法的基本假设与局限性,并介绍了可用于MOF孔道结构表征的其他替代或互补技术。
表面积和孔体积的相关性
准确表征多孔材料对于理解其潜在应用至关重要。表面积和孔体积是重要的定量指标,可为金属有机框架材料(MOF)在气体吸附、分离、催化和传感等多种应用中的性能提供深入见解1。
金属有机框架(MOF)的比表面积是一个参数,用于量化其可用于与客体分子相互作用的表面大小,并可影响其在各种应用中的性能2,3。在气体吸附应用中,MOF的比表面积反映了结合位点的可及性与亲和力,这与其分离性能直接相关4。在催化应用中,MOF的比表面积可影响活性位点的数量及其对反应物分子的可及性,从而影响其催化活性5。在传感应用中,活性位点的数量和可及性同样重要,因为客体分子与活性位点的相互作用越多,灵敏度越高(选择性也可能随之提高)6。比表面积还可能影响MOF在极端条件下的稳定性,因为较大的比表面积可能意味着表面缺陷数量更多7。
金属有机框架(MOF)的孔体积是一个量化其多孔结构内空隙空间多少的参数。该参数定义为MOF中所有孔隙的总体积,包括开放(可接触)孔隙和封闭(不可接触)孔隙。MOF的孔体积会影响其在气体吸附、分离和催化等多种应用中的性能。与比表面积类似,MOF的孔体积直接关系到其气体吸附与储存能力,以及允许客体分子到达吸附或催化位点的能力8。
利用氮气吸附测定比表面积和孔体积
比表面积和孔体积通常采用气体吸附技术测定,最常用的是氮气吸附法。在Brunauer-Emmett-Teller(BET)分析中选择氮气作为吸附质,是因为氮分子具有四极矩,其分子取向依赖于吸附剂的表面化学性质,从而有利于单分子层的形成。通过绘制氮气吸附量随压力变化的等温曲线,可以获得关于金属有机框架材料(MOF)表面和孔径尺寸的信息。利用吸附数据可计算材料的比表面积和总孔体积9。本方法的总体目标是获取氮气吸附数据,并利用这些数据计算MOF的比表面积和孔体积。
BET 方法10 是一种广泛用于测定多孔材料比表面积的技术,其原理基于气体在固体表面的吸附量与表面积、气体分子特性以及系统性质相关。在一定压力范围内,向样品材料引入已知量的吸附质气体(如氮气),并在每个压力增量下测量吸附到表面的气体量。通过关联吸附质的吸附量、压力和单层吸附容量,利用所得数据计算比表面积,该关系由 BET 方程9 表示:
(公式 1;eq. 1)
其中:
p = 吸附质的平衡压力(Pa)
p0 = 吸附质的饱和压力(Pa)
n = 吸附质吸附量(m3/g)
nm = 单层吸附容量(m3/g)
C = BET 常数(无量纲)
单层容量与总表面积之间的关系由以下公式给出:
(方程 2;eq. 2)
其中:
St = 金属有机框架(MOF)总比表面积(m2)
nm = 单层吸附容量(m3/g)
NAv = 阿伏伽德罗常数(molecule/mol)
scs = 吸附质分子的横截面积(m2/molecule)
Vmolar = 吸附质的摩尔体积(m3/mol)
Barrett-Joyner-Halenda(BJH)法11 是一种利用脱附数据计算总孔体积的常用方法。与 BET 法类似,将已知量的吸附气体(通常为氮气)引入样品中。随后逐步降低吸附质的分压,并测量每一步脱附的气体体积。该方法假设每个孔道中的脱附首先发生在毛细孔体积内,随后是吸附层厚度的减小。BJH 方程将脱附体积与吸附层厚度、孔半径和孔体积关联起来。这种关系可通过 BJH 孔径分布图表示,该图以孔半径为横坐标,孔体积为纵坐标。通过对孔径分布曲线在孔尺寸上进行积分,可确定总孔体积。BJH 方程12 可表示为:
(方程3;eq. 3)
其中:
n = 脱附步骤(无量纲)
vn = 毛细管凝聚液被排空的孔隙体积(m3)
ΔVn = 从孔隙中脱附的吸附质体积(m3)
Δtn = 吸附层厚度的变化(m)
A = 参与脱附过程的孔隙表面积(m2)
Rn = 依赖于平均孔径的 BJH 常数(无量纲)
c = 依赖于平均吸附层厚度的 BJH 常数(无量纲)
1. 样品制备
2. 实验文件设置
3. 进行吸附量测量
4. 氮气吸附测量
5. 数据分析
按照本实验方案操作后,可对获得的吸附等温线进行分析,并由此推导出材料的关键特性。氮气吸附实验的结果能够提供有关特定吸附剂的比表面积、孔体积和孔结构的关键信息。本实验旨在研究利用氮气吸附法测定纳米多孔金属有机框架材料UiO-66的比表面积和孔体积。UiO-66是一种典型的锆基MOF,具有高比表面积和优异的稳定性。尽管许多MOF材料在热、机械和化学稳定性方面较弱,但UiO-66由于其锆氧立方八面体金属节点的存在而表现出极高的稳定性,该结构可在BDC连接配体上提供12个配位延伸位点。该材料的结构包含7.5 Å的四面体孔笼和12 Å的八面体孔笼14,15。
无缺陷的UiO-66表现出I型吸附等温线特征16。I型等温线表明材料为微孔固体,其外表面相对较小。在I型等温线中,吸附量迅速趋近于一个极限值,说明氮气吸附量主要由吸附质可进入的微孔体积决定,而非由内表面积决定。在低P/P0下出现急剧的吸附上升,表明吸附剂与吸附质在狭窄的微孔中存在较强的相互作用17。I型等温线中通常不出现滞后环,因为滞后环常见于物理吸附的多层吸附区域,与孔道中的毛细凝聚现象相关。在低P/P0范围内氮气在吸附剂表面形成单层吸附,与吸附剂的比表面积有关;而在P/P0接近1时的孔道填充则与材料的总孔体积相关17。
BET 方法的应用通常在吸附仪器的软件中完成。然而,分析和计算也可以手动进行,或使用其他可调整的计算程序和方法,以获得关键结果。将 BET 模型应用于获得的氮气吸附等温线时,有两个关键步骤。首先,必须将氮气等温线转换为 BET 图,从中推导出 BET 单层容量。接下来,通过单层容量并选择适当的分子横截面积值来计算 BET 比表面积17。这一步骤通常在氮气吸附仪器的软件中完成。图 2展示了 UiO-66 的氮气吸附等温线。该等温线为 I 型,表明具有微孔结构以及氮气单层的形成。在较高相对压力下出现的陡峭台阶,导致呈现轻微的 II 型等温线特征,这表明存在多层吸附,同时也反映出由于 UiO-66 中的缺陷工程而形成了较大的介孔或大孔。在较高相对压力下观察到的滞后现象,表明存在较大的介孔和大孔结构。表 1列出了 BET 分析所得的数值。
使用 BET 方法时,必须满足 Rouquerol 判据。Rouquerol 判据指出:对变换后的 BET 数据进行线性拟合必须得到一条直线;若方法处于合适的分析范围内,则 C 值应始终为正值;Rouquerol 变换值必须随相对压力的增加而增加;单层吸附容量必须在用于拟合 BET 参数的数据范围之内18。为将 BET 方法应用于特定吸附等温线,需选取拐点处的线性区间。对于介孔材料,该区间通常位于 P/P0 为 0.05–0.30 的范围内,而对于许多微孔材料,则通常取自 P/P0 为 0.005–0.03 的范围。然而,实际的线性区间往往更为受限,因其取决于材料本身及分析温度。因此,线性区间的选取需要进行定性评估,类似于表 1中所示参数(C 值为正且相关系数接近 1,表明分析范围合适)。同样,线性区间内必须有足够的实验数据点(至少 10 个),以确保分析的可靠性。这些考虑也揭示了 BET 方法的固有局限性。C 是一个与单层形成时相对压力相关的常数,是用于定义 BET 方法中统计学单层形成假设下表面未被单层覆盖部分比例的指标。因此,较大的 C 值对应更高的表面覆盖率和更均匀的单层形成。当 C 值小于 2 时,等温线属于第 3 型或第 5 型,此时 BET 方法不适用。当 C 值小于 50 时,单层与多层吸附存在显著重叠。C 值至少为 80 表明等温线拐点较陡峭,单层吸附已基本完成,多层吸附开始发生。C 值大于 150 通常与窄微孔的填充或高能表面位点上的吸附相关17。
UiO-66 是一种微孔金属有机框架材料(MOF),通常存在缺陷,这些缺陷可增加比表面积并改善某些理想的吸附性能,但可能导致稳定性与结晶度降低15。通过不同程度的缺陷调控,缺陷型 UiO-66 框架的 BET 比表面积可达 1000–1800 m2/g,孔体积为 0.40–0.90 cm3/g15,16。
对于所测得的UiO-66,当采用0.01–0.05的线性P/P0范围时,其BET比表面积为1211 m2/g,C值为457。模拟的无缺陷UiO-66的理论比表面积为1200 m2/g14。在UiO-66所呈现的I型吸附等温线中,由于BET模型无法验证BET单层容量的有效性,因此BET比表面积应被视为表观比表面积17。所测得的比表面积处于UiO-66的预期范围内,结合C值,表明该材料具有微孔结构,并实现了均匀的单层吸附及孔道填充。
材料的孔体积通常在 P/P0 为 0.80–0.95 的条件下进行分析。如果材料中存在大孔,则在 P/P0 接近 1 时,氮气吸附等温线不会接近水平,因此无法评估总孔体积17。在此情况下测得的孔体积仅反映微孔和介孔的孔体积。
在 P/P0 为 0.80 时测得的 UiO-66 孔体积为 0.86 cm3/g。UiO-66 的理论孔体积为 0.77 cm3/g15。所测得的 UiO-66 样品具有较高的孔体积,最可能的原因是 UiO-66 结构中存在缺陷。该材料并非仅含有微孔,而是因存在缺陷而形成了更大的介孔或大孔,从而导致孔体积增大。这一结论得到了氮气吸附等温线形状的支持,即在较高相对压力下出现急剧上升和滞后现象,并呈现 I-2 型等温线特征。
对于特定材料,其测得的 BET 比表面积和孔体积通常处于某一特定范围内。已有研究表明,氮气吸附等温线及比表面积测量的重复性在不同文献中存在显著差异19。这主要是由于所选 BET 范围不同、材料缺陷、未进行重复实验以及模型本身的固有特性所致。BET 表面积识别(BETSI)程序等工具可根据扩展的选择标准自动确定线性范围,从而对 BET 比表面积进行明确评估。尽管 BET 模型已成为材料表征的标准方法,但它并非为微孔材料中的吸附过程而开发,原因在于其基于单层吸附和理想化吸附行为的概念13。BET 模型假设吸附过程均匀且表面均一。然而,对于表面不均一或孔径极小的材料,这些假设可能不成立,因此必须针对每种具体材料评估 BET 模型的适用性。
氮气吸附实验与分析结果表明,成功形成了具有轻微缺陷的UiO-66微孔晶体结构。计算得到的比表面积和孔体积处于文献报道的范围内15,19,因此可以得出结论:BET模型适用于金属有机框架材料UiO-66,并且在满足相应假设和条件的情况下,该模型也可推广应用于其他纳米多孔材料。

图 1:吸附仪示意图。 密封的样品管连接至压力传感器、真空系统以及自由空间/分析气体源。 请点击此处查看此图的放大版本。

图 2:UiO-66 在 77 K 下的氮气吸附与脱附等温线。 金属有机框架材料 UiO-66 在 77 K 下的氮气等温吸附曲线,测得其 BET 比表面积为 1211 m2/g,孔体积为 0.86 cm3/g。 请点击此处查看该图的放大版本。
| BET 比表面积 | 1211 m2/g |
| 斜率 | 0.0035 g/cm3 STP |
| Y 轴截距 | 0.000008 g/cm3 STP |
| C 值 | 457 |
| 单层吸附容量 | 278 cm3/g STP |
| 分子横截面积 | 0.1620 nm2 |
| 相关系数 | 0.9999 |
表1:UiO-66在77 K下BET分析所得结果的汇总表。 该表格总结了金属有机框架材料UiO-66在P/P0为0.01–0.05范围内进行BET分析所获得的关键参数。正值的C常数和y轴截距,以及0.9999的相关系数,表明为BET分析所选的线性区域是可接受的。
适用性与局限性
BET 方法需要满足几个关键假设:(1)表面为平面且均匀;(2)表面是均质的,所有吸附位点具有相同的能量;(3)吸附质形成单分子层。因此,BET 方法可能不适用于非多孔材料、具有复杂表面结构的材料(如不同类型的表面位点、不规则的表面形貌、能量差异较大的位点),或不表现出单层吸附行为的材料。若实际情况与上述假设存在较大偏差,可能会影响比表面积计算的准确性。与 BET 方法类似,BJH 方法也假设吸附均匀且表面均质,并进一步假设孔道为刚性圆柱形孔。因此,该方法同样可能不适用于具有复杂表面或可渗透结构的材料20。此外,由于压汞法或多孔分析需要进入孔隙空间,其计算结果不会包含封闭孔隙的体积。
对于微孔材料,应谨慎使用 BET 和 BJH 方法。BJH 方法未考虑流体-表面相互作用以及孔道内吸附质分子之间的相互作用,而这两种作用在较小的孔道中更为显著。因此,BJH 方法仅适用于介孔和小的宏孔。由于微孔通常表现出孔道填充行为,难以在等温吸附曲线上确定进行 BET 计算所需的线性区域21。
两种方法的另一个局限性在于其对样品制备方法的敏感性。样品必须处于分散状态,例如粉末或薄膜,而均匀制备此类样品可能具有挑战性。这会引入测量误差,并导致重复性困难。比表面积和孔体积也可能受到样品制备方法及条件的影响,例如材料合成技术、活化方法/条件,或干燥温度/时间22。
相对于替代方法的意义
由于氮气具有四极矩,其分子取向依赖于吸附剂的表面化学性质,从而有利于单分子层的形成,加之成本较低,因此氮气是 BET 和 BJH 分析中的标准吸附质17。然而,氩气和二氧化碳23也可用于此类分析,尤其适用于微孔结构。 氩气化学性质惰性,且为对称的单原子分子;但其三相点温度高于 77 K,因此体相参考状态存在疑问,且氩气单分子层的结构在很大程度上依赖于吸附剂的表面化学性质17。
由于 BET 和 BJH 方法并非普遍适用,因此应考虑其他测量比表面积和孔体积的方法。可分别采用 Langmuir 作图法、t-作图法或 Horvath-Kawazoe 方法来测定微孔比表面积、孔体积和孔径分布。非局域密度泛函理论(NLDFT)建模也可用于分析孔径分布,尤其适用于微孔,因其考虑了流体密度随孔径变化的情况。压汞法可用于测定孔隙率和孔体积,但需注意该技术的可测范围,因其无法进入微孔。计算方法可用于计算理论表征参数,并为实验结果提供对比依据,这对具有封闭孔隙的材料尤为有用。尽管 BJH 法可得到孔径分布,但其未考虑孔径的非均匀分布,也无法全面表征孔之间的连通性。可结合使用其他表征手段,如扫描电镜(SEM)、透射电镜(TEM24)或 X 射线衍射(XRD),以更全面地理解多孔材料的结构。即使某些材料无法通过 BET 或 BJH 方法完全表征,这两种方法仍可用于不同材料之间的定性比较。氮气孔隙率分析法若与其他技术结合使用,可成为一种非常有效的工具。12
作者无竞争性财务利益。
本工作由美国能源部科学办公室基础能源科学项目资助的能源前沿研究中心——用于能源的材料酸性气体诱导演化机制理解与控制中心(UNCAGE-ME)支持,项目编号为 DE-SC0012577。J.S. 认明,本材料基于美国国家科学基金会研究生研究奖学金资助项目(资助号 DGE-2039655)支持的工作。本文所表达的任何观点、发现、结论或建议均为作者(们)的观点,不一定反映美国国家科学基金会的观点。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 吸附仪 | Micromeritics | TriStar II Plus | |
| 吸附分析软件 | Micromeritics | TriStar II Plus 版本 3.03 | |
| 天平 | |||
| 杜瓦瓶 | 液氮 N2 杜瓦瓶 | ||
| 二甲基甲酰胺(DMF) | 费舍尔科技公司 | D119-1 | |
| 氦气 | Airgas | HE UHP300 | 超高纯度 |
| 氮气 | Airgas | NI 230LT22 | 工业级液氮 N2 |
| 氮气 | Airgas | NI UHP300 | 超高纯度气态氮 N2 |
| 样品管 | Micromeritics | 302-61001-02 | 玻璃样品管 |
| 样品预处理系统 | Micromeritics | 061-00030-00 | VacPrep 061 |
| 对苯二甲酸(H2BDC) | 西格玛奥德里奇 | 185361 | |
| ZrCl4 | 西格玛奥德里奇 | 221880 | 四氯化锆,≥99.5% 痕量金属含量 |