
The red dashed lines represent two thresholds for adsorption capacity (NH2S+CO2 > 5.0 mol/kg) and selectivity (SH2S+CO2/C1-C3 > 85). The figure contains the data of 6013 CoRE-MOFs
随着人们生活水平的提高及环保意识的增强, 大部分城市对天然气的需求日益增加, 由于其成本低, 利用效率高, 污染物排放量少[1].天然气在民用燃料:取暖、烹饪和工业燃料:发电、制造氮肥及代替汽车用油等方面有着不可磨灭的作用, 一方面可以减缓能源紧缺, 另一方面具有投资少、成本低、污染少等优点[2].但是大多数天然气中含有大量硫元素, 其中H2S含量大于1%的天然气储量占到总储量的1/4, 普光气田开采出的天然气中H2S含量甚至高达13%~18%, 有机硫化合物高达340.6 mg/m3, 此外还包括其他一些酸性气体[3, 4], 如: CO2, 这些酸性气体(H2S和CO2)的存在不仅会严重腐蚀管道设备[5], 燃烧后有害物质的排放还会对人体健康造成不可忽视的危害, 因此对天然气中H2S和CO2的分离和去除成为天然气净化领域的研究重点.
分离这类酸性气体, 主要有吸收和吸附两种方法, 吸收脱除的方法主要有湿法、干法和生物法[6], 但这种方法耗能高, 再生难且容易腐蚀设备.相比较而言, 使用多孔材料的吸附法具有节能、操作条件温和、设备简单、对环境影响低等优点, 被认为是一种有潜能的可替代方法[7].过去常用的吸附剂有沸石、活性炭、介孔硅.且普遍对多孔沸石进行研究, 例如:有采用不同的Si/Al比来调整对CH4中CO2的最佳吸附性能[8]. Shah等[9]用分子模拟的方法研究了386种沸石在CH4, C2H6, H2S混合气体中对H2S的吸附分离性能, 并在最接近实际环境状态下进一步筛选出了16种最优的沸石.在材料领域的快速发展中, 一种新型有机-无机杂化多孔材料——金属-有机框架(metal-organic frameworks, MOFs)深受人们关注, MOFs是通过含O或N等多齿有机配体与无机金属离子自组装而成的多孔材料, 因其与传统材料活性炭、沸石等相比具有比表面积大、孔隙率高、结构和性质可调控等特性[10, 11].这些优良性质使其广泛应用于气体吸附[12, 13]和存储[14, 15]、分离和捕获[16]、磁学、光学[17]、催化[18]等方面, 并展现了良好的应用前景.在利用MOFs的吸附和分离性能针对天然气净化的领域, 研究者们做了大量相关的工作.例如, 杨江峰等[19]采用一维菱形孔道的MOFs材料Cu(INA)2作为吸附剂, 实验和模拟相结合, 发现这种材料可以高选择性地吸附CH4/N2混合物中的CH4, 并且对CH4/N2的分离效果优于CO2/CH4.实验在三种不同温度(25 ℃, 15 ℃和0 ℃)下其选择性分别为: CH4/N2: 8, 9和10.6, CO2/CH4: 3.9, 4和5.范卫东等[20]对自组装棒状晶体Zr-MOFs的CH4/N2吸附分离性能进行研究, 结果表明在273 K和0.1 MPa时, 材料对CH4吸附效果最好, CH4的吸附热在20 kJ/mol左右, 易于吸附剂的再生. Wang等[21]采用酸功能化的UiO-66分离H2S/CH4, 研究表明在不同的H2S/CH4的浓度配比及不同压力下, 其吸附选择性范围约为50~130, 吸附选择性随压力变化趋势呈现S型. Joshi等[22]研究了六种不同的UiO-66和MILs系列材料,并探讨了H2S/CO2/CH4和H2S/CH4在这六种不同材料中对H2S的吸附量, 分别约为: 0.15~0.45 mol/kg和0.2~0.55 mol/kg. Belmabkhout等[23]研究发现具有Ga3+, Al(Ⅲ)骨架离子结构的MOFs对H2S有很高的耐腐蚀性, 并能显著增强H2S吸附和分离性能. Bhatt等[24]研究表明fcu-MOF对H2S和CO2有很高的吸附量和选择性, 三种不同fcu-MOF对H2S吸附量分别为: 0.9, 1.1和1.5 mol/kg.近二十年, 被合成的MOFs已有数千种, 仅通过实验来逐一测试CO2和H2S的分离能力几乎是不可实现的[25].因此, 除了实验工作以外, 通过分子模拟手段的高通量筛选技术已经被许多研究者用来深入探讨MOFs的气体吸附分离性能.例如, 仲崇立课题组[26]研究了105种MOFs分离CH4/H2, 研究中最优材料的选择性和吸附量都远高于传统活性炭、沸石等吸附剂. Smit和合作者[27]对ZIFs的CO2/CH4, CO2/N2和CH4/N2的吸附选择性进行了模拟研究.发现吸附选择性随压力增大而减少, 最后趋于稳定.卞磊等[28]采用理论计算和实验相结合的方式, 从2932种MOFs中筛选出了三种对甲醛有较高吸附量的材料, 进行深入地实验测试. Actintas等[29]采用模拟计算研究了4240种MOFs对H2/CH4的分离性能, 从中筛选出10种性能最佳的MOFs.许红等[30]对33种MOFs去除天然气中H2S的能力进行模拟计算, 研究结果表明改性官能团和小孔作用是影响选择性的关键因素, 其中-Cl、-OH、-OCH3基团对H2S有较强的吸附作用力.
如何高效选取最佳吸附剂是这一领域亟待解决的关键问题.本文采用分子模拟方法研究了MOFs对天然气中酸性气体的吸附净化能力, 以期寻找到拥有共吸附H2S和CO2能力的最佳吸附剂.首先, 我们对比了三种权衡吸附量和选择性的方法, 通过皮尔森相关系数(Pearson)择优选取一种方法用于之后的研究; 接着, 采用两种统计学方法(多元线性回归分析和决策树模型)对四种MOFs描述符[最大孔径(LCD), 孔隙率(ϕ), 体积比表面积(VSA)和吸附热(Qst0)]进行分析, 寻找一条吸附性能最优材料的路径.最后, 探讨了20种最优MOFs的组成对分离性能的影响.
此工作选择的研究对象为6013种实验已合成的MOFs (CoRE-MOFs), 其细节在实验部分被详细阐述.通过分子模拟, 我们计算了它们的三种几何描述符(LCD, ϕ和VSA), 一种能量描述符(等量吸附热, Qst0)和五元混合物中各组分的吸附量, 并分析这四种描述符对H2S和CO2的总吸附量(NH2S+CO2=NH2S+NCO2)和它们相对于CH4, C2H6和C3H8 (C1-C3)的选择性(SH2S+CO2/C1-C3 =NH2S+CO2×yC1-C3 /NC1-C3×yH2S+CO2)的对应关系. 图 1为两者的关系图, 选择性高的MOFs相对吸附量都较低, 随着吸附量的增大选择性逐渐降低.图中还可观察到, 大部分SH2S+CO2/C1-C3>1, 这是因为H2S是极性分子, 且CO2分子具有较强的四极矩, 因此相较于饱和烷烃, H2S/CO2与MOFs具有更强的作用力, 从而优先吸附.
众所周知, 优异的吸附剂需要同时具备高的吸附量以及足以分离目标气体的选择性.例如, 我们可以基于两个限制条件(NH2S+CO2>5.0 mol/kg和SH2S+CO2/C1-C3>85)来选取吸附分离性能较为优异的材料, 即图 1右上角的阴影部分包含的MOFs.但这样的方法太过绝对, 会遗漏一些性能优异的MOFs, 而且也无法比较被选择材料的优劣.为此, 针对吸附量和选择性, 我们希望寻找到一种与MOFs的结构/能量参数具有强相关性的权衡变量, 用以判断MOFs的综合性能.所以, 三种权衡吸附量和选择性的方法被采用和对比:权衡α法(Tradeoff between SH2S+CO2/C1-C3 and NH2S+CO2, TSN), 标准值法(Standard normalization method, SNM)和权衡β法(Tradeoff between selectivity and capacity, TSC), 计算公式如下:
$ {\rm{TSN = }}{N_{{{\rm{H}}_{\rm{2}}}{\rm{S + C}}{{\rm{O}}_{\rm{2}}}}} \cdot {\rm{ln}}\left( {{S_{_{{{\rm{H}}_{\rm{2}}}{\rm{S + C}}{{\rm{O}}_{\rm{2}}}{\rm{/C1-C3}}}}}} \right) $ |
(1) |
$ {\rm{SNM = }}{\log _{10}}\left( {{Z_{{N_{{{\rm{H}}_{\rm{2}}}{\rm{S + C}}{{\rm{O}}_{\rm{2}}}}}}} \cdot {Z_{{S_{{{\rm{H}}_{\rm{2}}}{\rm{S + C}}{{\rm{O}}_{\rm{2}}}{\rm{/C1 - C3}}}}}}} \right) $ |
(2) |
其中,
$ \overline {{X_j}} {\rm{ = }}\frac{{\sum\nolimits_{i{\rm{ = }}1}^n {{X_{ij}}} }}{n}, S_j^2{\rm{ = }}\frac{{\sum\nolimits_{i{\rm{ = }}1}^n {{{({X_{ij}} - \overline {{X_j}} )}^2}} }}{{n - 1}} \\ {\rm{TSC = ln}}\left( {{N_{{{\rm{H}}_{\rm{2}}}{\rm{S + C}}{{\rm{O}}_{\rm{2}}}}}} \right) \cdot {\rm{lo}}{{\rm{g}}_{10}}\left( {{\rm{lo}}{{\rm{g}}_{10}}\left( {{S_{{{\rm{H}}_{\rm{2}}}{\rm{S + C}}{{\rm{O}}_{\rm{2}}}{\rm{/C1 - C3}}}}} \right)} \right) $ |
(3) |
上述公式(2)中Z表示标准化之后的吸附量和选择性, Sj表示吸附量和选择性的标准偏差, Xij表示每一种MOF对应的初始吸附量或选择性.
本研究选取NH2S+CO2大于1.2 mol/kg的4368种CoRE-MOFs用于这三种方法的对比分析, 因为吸附量太低时, 即使选择性非常高, 这类材料也无法被实际应用.在这三种方法中, TSN法是通过对选择性进行对数计算来平衡选择性和吸附量的数值关系, 使两者相差不大, 能更好的体现哪种变量作用显著, 这种方法在Shah等[9]的工作中被用于分子筛吸附剂脱除H2S性能的评估. SNM法是一种对大数据进行预处理简便、准确的方法, 有简单缩放, 逐样本均值消减, 特征标准化, 加快最优解的求解速度及提高精度等特征, 并且采用SNM法分析, 当数据因操作条件略有变动或出现手动误差时, 对分析结果影响都很小, 目前在深度学习和机器学习方面有较多的应用.最后, TSC法是本工作首次提出的权衡方法.它是在TSN法的基础上进行修改, 对SH2S+CO2/C1-C3进行两次以10为底的对数计算, 主要目的还是进一步降低选择性的影响, 使吸附量和选择性对于权衡变量的影响程度保持平衡, 不会仅因为选择性过大而判断材料性能为优.值得注意的是, 当NH2S+CO2>1.2 mol/kg时, 所取的4368种CoRE-MOFs的SH2S+CO2/C1-C3范围是: (1.02~1.0)×107.由于选择性数值都大于1, 因此所有选择性均可进行两次以10为底的对数计算.仅当SH2S+CO2/C1-C3<10, TSC为负值, 且认为这些MOFs的性能不佳.三种方法进行比较可看出, 选择性与吸附量数量级相差不大时TSN和SNM方法适用; 当选择性的数量级远大于吸附量时, TSC更适用.
图 2展示了四种MOFs描述符(LCD, ϕ, VSA和Qst0)与TSN, SNM和TSC三种权衡变量的关系.四种描述符中Qst0采用H2S的吸附热进行研究, 采用CO2的吸附热也有类似的趋势, 如图S1所示.前三种描述符(LCD, ϕ和VSA)与所有权衡变量都呈现负相关关系, 而吸附热则呈正相关关系.相比TSN和SNM权衡方法, TSC对四种描述符的散点更加集中.进一步对图 2中数据做线性拟合并使用Pearson法计算其相关系数也证实了上述结论, 如图 2中红线和表 1. Pearson相关系数可以看出:前三种几何描述符和TSN, SNM和TSC的相关系数为负值, 只有能量描述符——吸附热的相关系数为正值.这主要是因为, 当孔径较小时, 即LCD较小, 小尺寸的H2S和CO2气体分子恰好能进入材料孔道内, 被材料孔壁所包围, 接触面积大导致强的相互作用力[31], 故而吸附分离性能较好, ϕ也有类似的关系.然而, 当孔径变大时, LCD和ϕ也逐渐增大, 气体分子与孔壁分子间作用力随着孔径的增大而减弱, 从而导致脱附以及无选择的吸附增多.对于VSA, 随气体分子与孔壁的可接触面积逐渐增大, 将会增强大分子的吸附, 例如C2H6或C3H8, 从而导致分离性能的下降.对比三种权衡法的Pearson相关系数, 可知TSC与四种描述符的相关程度更大.故TSC法更适合用于权衡本工作中H2S与CO2的吸附量之和以及它们相对于C1-C3的选择性, 其权衡后的结果置信度更高.因此, 在下一步的工作中我们选取TSC这一权衡变量进行深入研究.
(a), (d), (g), (j) represent the effects of four descriptors on TSN; (b), (e), (h), (k) represent the effects of the four descriptors on SNM; (c), (f), (i), (l) represent the effects of the four descriptors on TSC. Red lines represent the linear fitting line. The figures contain the data of 4368 CoRE-MOFs with NH2S+CO2 > 1.2 mol/kg
MOF descriptors |
Pearson correlation coefficient |
Standardized regression coefficient | |||
TSN | SNM | TSC | TSC | ||
LCD | -0.340 | -0.419 | -0.613 | -0.209 | |
ϕ | -0.393 | -0.569 | -0.717 | -0.327 | |
VSA | -0.379 | -0.530 | -0.673 | -0.113 | |
Qst0 | 0.480 | 0.434 | 0.536 | 0.368 |
在材料设计和选择时, 我们需要综合考虑每一种描述符对MOFs性能的影响程度, 从而选取影响最大的变量进行材料设计路线的制定.故进一步采用多元线性回归分析(multiple linear regression, MLR)的方法[32]研究了四种描述符共同作用下对TSC的影响程度, 公式如下:
$ Y = {\alpha _0} + {\alpha _1}{X_1} + {\alpha _2}{X_2} + \cdots + {\alpha _n}{X_n} + \varepsilon $ |
(4) |
其中Y是因变量, 表示TSC, Xn是自变量, 表示四种描述符, αn是Xn的标准回归系数, ε是Y值实际和预测值之间的残差.
接着采用决策树模型规划优异MOFs的设计路线, 这两种方法也在之前的工作中被有效地运用在硫醇分离材料[33]和MOF膜的设计[34]中.在MLR中, 把TSC作为一个因变量, 四种描述符作为自变量, 我们先对这五种变量进行标准化处理, 从而得出标准化时每种描述符对应的标准回归系数, 列于表 1.标准回归系数的绝对值越大说明这种自变量对因变量(TSC)的影响越大.从表 1知, 吸附热的绝对值最大, 且为正, 说明Qst0对其吸附性能的影响最大, 且为正相关; ϕ其次, 而VSA最小.因此, 相比其他三种描述符, VSA的影响可忽略.
从多元线性回归分析得到Qst0, LCD和ϕ三种描述符对MOFs性能影响较大.接着我们采用一个三层树状图来分析这三种描述符的最优范围以达到最大概率发现优异MOFs的目的, 并直观地建立一条优异材料的设计路线. 图 3展现了这三种描述符的决策树路径, 设置目标: TSC>0.46 (由图 1中阴影部分NH2S+CO2=5 mol/kg和SH2S+CO2/C1-C3=85, 计算得TSC≈0.46), 有531种MOFs (12.16%)被选择.决策树模型通过三种描述符Qst0, LCD和ϕ来筛分MOFs. 图 3中一条清晰的路径可以被发现:当Qst0>49.3 kJ/mol, ϕ<0.308, LCD<3.88 ,大约77.6%的材料具有TSC>0.46, 即在这个范围内的材料具有最大概率拥有良好的H2S和CO2的吸附能力.通过决策树模型可以清晰地观察在哪条路径下能最大概率地寻找到性能最佳的CoRE-MOFs.
基于以上多种统计学分析, 为筛选出性能最优的CoRE-MOFs, 以TSC≥1.35为限制条件, 20种材料脱颖而出, 它们的吸附数据被列于表 2.表中可得No. 15材料(QUQGAL)对H2S和CO2的总吸附量最大, 为13.82 mol/kg, 而选择性却相对较低.选择性最大超过107, 说明它几乎不吸附C1-C3, 仅吸附H2S和CO2.综合考虑选择性和吸附量, 这些MOFs均优于常见分子筛[9].分析这20种材料对H2S的吸附热, 发现除了No. 13材料(FASJAL), 其余都大于50 kJ/mol.这是由于FASJAL只吸附微量H2S.虽然此材料仍倾向于吸附CO2, 但是它不符合共吸附的原则, 因此这一材料不符合此研究的目的.对这些MOFs的组成——金属中心的分析发现, 它们当中75.0% (15/20) MOFs含有碱金属和碱土金属, 而整个数据库中仅有7.9% (476/6013)的MOFs包含这两族金属.因此, 证明金属性强的这两族金属在MOFs中对酸性气体H2S和CO2有更强的亲和力, 从而增强材料对H2S和CO2的吸附分离能力[35].进一步对这20种材料结构的分析发现: 7种MOFs拥有开放金属点位(open metal site, OMS), 1种MOFs拥有游离金属离子(non-framework ion, NFI).这些特殊的金属状态对于增强H2S和CO2的吸附和分离有进一步的促进作用, 正如常见的开放金属点位MOFs——Mg-MOF-74[36].
No. | CSD code | TSC | NH2S/ (mol•kg-1) |
NCO2/ (mol•kg-1) |
SH2S+CO2/C1-C3 (105) |
Qst, H2S0/ (kJ•mol-1) |
Qst, CO20/ (kJ•mol-1) |
ϕ | VSA/ (m2•cm-3) |
LCD/ | Metal | Metal statea |
1 | KAPHOZ | 1.83 | 3.68 | 5.08 | 101.62 | 83.65 | 55.76 | 0.42 | 971.54 | 3.65 | K/Zn | NFI |
2 | IQUNAJ01 | 1.67 | 2.64 | 4.68 | 84.86 | 67.79 | 62.41 | 0.25 | 0 | 2.83 | K/Y | OMS |
3 | VAYDOQ | 1.66 | 1.83 | 5.41 | 84.00 | 101.47 | 100.25 | 0.28 | 42.45 | 4.11 | Li | SS |
4 | IQUNAJ | 1.66 | 2.49 | 4.73 | 83.75 | 66.71 | 61.90 | 0.24 | 0 | 2.83 | K/Y | OMS |
5 | WENSIS | 1.65 | 1.59 | 5.57 | 83.08 | 65.95 | 61.73 | 0.25 | 0 | 2.75 | K/Ca/Zr | OMS |
6 | YUJVOQ01 | 1.63 | 2.10 | 5.56 | 22.18 | 74.57 | 91.51 | 0.28 | 11.33 | 3.88 | Li/Zr | SS |
7 | YUJVOQ | 1.62 | 1.73 | 5.52 | 35.79 | 73.05 | 81.12 | 0.26 | 3.15 | 3.80 | Li/Zr | SS |
8 | YUJWAD | 1.60 | 1.51 | 5.41 | 52.77 | 71.08 | 78.67 | 0.24 | 0 | 3.70 | Li/Zr | SS |
9 | TUMGOX | 1.55 | 1.29 | 5.09 | 74.01 | 63.91 | 50.10 | 0.27 | 0 | 3.47 | Ce | SS |
10 | PAMGIV | 1.55 | 2.80 | 5.29 | 3.22 | 80.26 | 88.17 | 0.36 | 399.52 | 4.35 | Mn | SS |
11 | BEVQUP | 1.52 | 2.51 | 3.67 | 71.67 | 76.35 | 53.91 | 0.14 | 0 | 2.88 | Rb/Y | SS |
12 | SOPPOD | 1.44 | 2.16 | 3.79 | 28.74 | 66.65 | 62.41 | 0.25 | 0 | 2.85 | K/Tb | OMS |
13 | FASJAL | 1.43 | 0.00080 | 5.78 | 35.07 | 28.33 | 48.57 | 0.22 | 0 | 3.28 | Ca/Zr | SS |
14 | VIXGEQ | 1.43 | 0.69 | 6.17 | 3.10 | 62.85 | 93.35 | 0.27 | 266.29 | 5.15 | Zr | SS |
15 | QUQGAL | 1.42 | 4.88 | 8.94 | 0.030 | 98.99 | 79.14 | 0.47 | 893.42 | 5.24 | Li | OMS |
16 | IDAZEU | 1.42 | 1.65 | 3.83 | 63.49 | 69.95 | 66.39 | 0.25 | 0 | 2.74 | K/Yb | OMS |
17 | ZADDAJ | 1.38 | 1.47 | 5.28 | 1.88 | 77.26 | 93.13 | 0.27 | 0 | 3.58 | Li/Fe | SS |
18 | VIXGAM | 1.37 | 0.63 | 6.53 | 0.98 | 60.80 | 86.40 | 0.26 | 284.90 | 5.20 | Zr | SS |
19 | HIYMAF | 1.37 | 1.62 | 3.56 | 60.00 | 83.97 | 88.44 | 0.29 | 442.95 | 4.38 | K/Dy | OMS |
20 | VIXGAM01 | 1.35 | 0.63 | 6.50 | 0.78 | 62.47 | 85.05 | 0.27 | 286.20 | 5.21 | Zr | SS |
a Metal state includes open metal site (OMS), non-framework ion (NFI) and saturated state (SS). |
本工作中, 我们模拟了6013种CoRE-MOFs在五元天然气混合物中捕获H2S和CO2的性能.对比了三种综合吸附量和选择性的权衡方法(TSN, SNM和TSC ), 通过Pearson相关系数法比较, 得到权衡变量TSC与四种描述符(LCD, ϕ, VSA和Qst0)的相关性最佳.接着采用多元线性回归分析法探讨这四种描述符对TSC的影响程度, 总体标准回归系数表明Qst0对TSC的影响程度最大, 而VSA最小.之后, 采用决策树模型对前三种最佳的描述符进行分析, 规划出了一条发现性能最优MOFs的路径: Qst0>49.3 kJ/mol, ϕ<0.308和LCD<3.88 ,根据此路径筛选出来的CoRE-MOFs有77.6%的概率拥有良好的吸附性能, 对新材料的合成有一定指导作用.最后我们列出20种TSC最大的MOFs并分析它们的组成发现, 其中15种MOFs含有碱金属和碱土金属元素.这一模拟工作, 从几何效应、能量效应以及MOFs的组成, 多方面系统地探讨了MOFs的结构-性能关系, 指明了最优的设计路线以及金属组成, 对MOFs在天然气同时脱除H2S和CO2领域的发展给予了清晰重要的理论指导.
采用6013种实验已经合成的MOFs为研究对象, 进行晶体模型搭建, 所有结构以及结构参数都来自Chung等[37, 38]的工作, 他们对超过350种拓扑结构的MOFs进行收集, 并移除了游离的溶剂分子.本工作使用的是最新的CoRE-MOFs版本2[38].接着, 采用不同的模拟计算方法对CoRE-MOFs的几何描述符(LCD, ϕ和VSA)进行计算.其中LCD用Zeo++[39]的方法计算, ϕ和VSA分别采用直径为2.58 的He和3.64 的N2作为探针在RASPA[40]软件包下进行计算, 此方法来自Wilmer等的工作[31, 41].最后, 吸附热是通过分别计算无限稀释状态下每个单一气体分子在不同MOFs中的吸附能, 计算方法为RASPA软件包[40]的NVT-MC方法.
采用巨正则系综蒙特卡洛(GCMC)模拟6013种CoRE-MOFs对天然气五元混合物的吸附分离性能.在GCMC模拟中, 整个材料是处于静态的, MOFs的原子固定不动[42], 故不需要考虑MOFs的键伸缩势能、键角弯曲势能及双面角扭曲势能等.但部分气体分子仍然包括这些势能.模拟时分子间的相互作用采用伦纳德-琼斯(Lennard-Jones, LJ)势能函数描述:
$ {u_{{\rm{LJ + elec}}}}(r) = \sum {4{\varepsilon _{ij}}} \left[{{{\left( {\frac{{{\sigma _{ij}}}}{{{r_{ij}}}}} \right)}^{12}}-{{\left( {\frac{{{\sigma _{ij}}}}{{{r_{ij}}}}} \right)}^6}} \right] + \frac{{{q_i}{q_j}}}{{4\pi {\varepsilon _{\rm{o}}}{r_{ij}}}} $ |
(5) |
其中uLJ+elec(r)表示原子i和j间非键相互作用能, 前一项是范德华非键结合势能, 后一项是库伦静电势能. σij表示原子间的平衡距离, rij表示原子间的距离, εij表示由能量最低点至能量为0的差, qi, qj是原子i和j所带的电荷, εo=8.8542×10-12 C2•N-1, 表示真空介电常数. CoRE-MOFs的LJ势能参数来自于UFF力场[43], 其力场参数列于表S1中. MOFs的原子电荷采用MEPO-Qeq的方法[44]能快速计算得到. H2S, CO2和C1-C3的力场参数来自于TraPPE力场[45], 列于表S2中.目前有很多研究已经验证了UFF力场, MEPO-Qeq电荷算法和TraPPE力场在不同MOFs材料中能精确预测气体的吸附和扩散[34, 46].同时, 我们也验证了此参数和计算方法的可靠性.从图S2可看出模拟的H2S吸附等温线均与已有的实验数据较吻合, 证明了本工作所用力场、电荷算法以及模拟方法的准确性.
对于CO2模型, C—O键长为1.16 ,其键角为180°. CO2和CH4是刚性的, 不具有键长、键角和二面角势能. CH4, C2H6和C3H8均使用联合原子模型, C3H8的键角弹性计算公式和参数列于表S3中.值得注意的是H2S的分子模型, 正如图 4所示, 它是一个四位模型, 其中S—H键长为1.34 ,在S和H原子上有L-J势能.另外还存在一个靠近S原子的虚原子(M), 且H原子和虚原子上都带有部分电荷, 而S原子未带电荷[47].这一模型已被证明能够准确模拟H2S在分子筛内的吸附分离性能[9].
这一工作中, 为了模拟真实状态, 天然气被假设为在298 K和10 bar状态下的五种混合气体(CH4, C2H6, C3H8, H2S和CO2), 它们的比例为70:10:5:5:10.采用GCMC方法分别模拟了6013个CoRE-MOFs对五元混合物的吸附性能, 每个GCMC模拟独立进行, 整个模拟过程中保持化学势、体积、温度恒定.通过Lorentz- Berthelot规则计算MOFs与吸附质分子间的相互作用.在三维体系中施加周期性边界, 模拟晶胞沿三维方向分别扩展到至少24 .为了计算LJ相互作用, 长程校正的球形截断半径设定为12 ,静电相互作用使用Ewald求和计算[48].每次模拟进行20000次循环, 前10000次用于平衡模拟体系, 后10000次用于对整体摆幅的平均.整个模拟在RASPA软件包[40]下操作.每个周期由n次试验运动组成(n:吸附分子的数量), 包括平移、旋转、再生长和交换, 其中交换运动即为插入和删除运动.为保证模拟准确性, 所使用的GCMC循环次数也被测试.图S3显示20种最佳CoRE-MOFs的气体吸附量在不同循环数的模拟中保持稳定, 因此证明进一步增加循环次数对模拟结果影响甚微.
Shahbaz, M.; Lean, H. H.; Farooq, A. Renew. Sust. Energy. Rev. 2013, 18, 87. doi: 10.1016/j.rser.2012.09.029
Schoots, K.; Rivera-Tinoco, R.; Verbong, G.; Zwaan, D. V. B. Int. J. Greenhouse Gas Control 2011, 5, 1614. doi: 10.1016/j.ijggc.2011.09.008
吴基荣, 毛红艳, 天然气工业, 2011, 31, 99.Wu, J. R.; Mao, H. Y. Nat. Gas. Ind. 2011, 31, 99.
杨婷婷, 熊运涛, 崔荣华, 肖俊, 韩淑怡, 天然气与石油, 2013, 31, 40. doi: 10.3969/j.issn.1006-5539.2013.02.011Yang, T. T.; Xiong, Y. T.; Cui, R. H.; Xiao, J.; Han, S. Y. Nat. Gas & Oil 2013, 31, 40. doi: 10.3969/j.issn.1006-5539.2013.02.011
樊辉, 陈陆建, 赵红, 甑建, 孙文成, 胡康宁, 天然气化工, 2011, 36, 34.Fan, H.; Chen, L. J.; Zhao, H.; Zeng, J.; Sun, W. C.; Hu, K. N. Nat. Gas. Ind. 2011, 36, 34.
王剑, 张晓萍, 李恩田, 马路, 王树立, 常州大学学报, 2013, 25, 88.Wang, J.; Zhang, X. P.; Li, E. T.; Ma, L.; Wang, S. L. J. Changzhou Univ. 2013, 25, 88.
张晓东, 李红欣, 侯扶林, 董寒, 朱正, 崔立峰, 功能材料, 2016, 47, 8178. doi: 10.3969/j.issn.1001-9731.2016.08.031Zhang, X. D.; Li, H. X.; Hou, F. L.; Dong, H.; Zhu, Z.; Cui, L. F. J. Funct. Mater. 2016, 47, 8178. doi: 10.3969/j.issn.1001-9731.2016.08.031
Palomino, M.; Corma, A.; Rey, F.; Valencia, S. Langmuir 2010, 26, 1910. doi: 10.1021/la9026656
Shah, M. S.; Tsapatsis, M.; Siepman, J. I. Angew. Chem. 2016, 128, 6042. doi: 10.1002/ange.201600612
Furukawa, H.; Cordova, K. E.; O'Keeffe, M.; Yaghi, O. M. Science 2013, 44, 974.
张秀芳, 安晓辉, 刘大欢, 阳庆元, 杨祝红, 仲崇立, 陆小华, 化学学报, 2011, 69, 84. http://manu19.magtech.com.cn/Jwk_hxxb/CN/abstract/abstract339696.shtmlZhang, X. F.; An, X. H.; Liu, D. H.; Yang, Q. Y.; Yang, Z. H.; Zhong, C. L.; Lu, X. H. Acta Chim. Sinica 2011, 69, 84. http://manu19.magtech.com.cn/Jwk_hxxb/CN/abstract/abstract339696.shtml
周建海, 赵会玲, 胡军, 刘洪来, 胡英, 化工学报, 2014, 65, 1680. doi: 10.3969/j.issn.0438-1157.2014.05.018Zhou, J. H.; Zhao, H. L.; Hu, J.; Liu, H. L.; Hu, Y. CIESC J. 2014, 65, 1680. doi: 10.3969/j.issn.0438-1157.2014.05.018
穆韡, 刘大欢, 阳庆元, 仲崇立, 物理化学学报, 2010, 26, 1657. doi: 10.3866/PKU.WHXB20100616Mu, W.; Liu, D. H.; Yang, Q. Y.; Zhong, C. L. Acta Phys.-Chim. Sin. 2010, 26, 1657. doi: 10.3866/PKU.WHXB20100616
吴选军, 赵鹏, 王杰, 刘保顺, 蔡卫权, 物理化学学报, 2014, 30, 2043. doi: 10.3866/PKU.WHXB201409222Wu, X. J.; Zhao, P.; Wang, J.; Liu, B. S.; Cai, W. Q. Acta Phys.-Chim. Sin. 2014, 30, 2043. doi: 10.3866/PKU.WHXB201409222
周子娥, 薛春瑜, 阳庆元, 仲崇立, 化学学报, 2009, 67, 477. http://manu19.magtech.com.cn/Jwk_hxxb/CN/abstract/abstract327080.shtmlZhou, Z. E.; Xue, C. Y.; Yang, Q. Y.; Zhong, C. L. Acta Chim. Sinica 2009, 67, 477. http://manu19.magtech.com.cn/Jwk_hxxb/CN/abstract/abstract327080.shtml
Qiao, Z. W.; Wang, N. Y.; Jiang, J. W.; Zhou, J. Chem. Commun. 2015, 52, 974.
Wu, P.; He, C.; Wang, J.; Peng, X. J. Am. Chem. Soc. 2012, 134, 14991. doi: 10.1021/ja305367j
Kong, G. Q.; Ou, S.; Zou, C.; Wu, C. D. J. Am. Chem. Soc. 2012, 134, 19851. doi: 10.1021/ja309158a
张倬铭, 杨江峰, 陈杨, 王勇, 李立博, 李晋平, 化工学报, 2015, 66, 3549Zhang, Z. M.; Yang, J. F.; Chen, Y.; Wang, Y.; Li, L. B.; Li, J. P. CIESC J. 2015, 66, 3549.
韩素英, 范卫东, 高荔, 曹运祥, 孙道峰, 石油与天然气化工, 2017, 46, 51.Han, S. Y.; Fan, W. D.; Gao, L.; Cao, Y. X.; Sun, D. F. Chem. Eng. Oil Gas. 2017, 46, 51.
Wang, S.; Wu, D.; Huang, H.; Yang, M.; Tong, M.; Liu, D.; Zhong, C. L. Chin. J. Chem. Eng. 2015, 23, 1291. doi: 10.1016/j.cjche.2015.04.017
Joshi, J.; Zhu, G.; Lee, J. J.; Carter, E. A.; Jones, C. W. Langmuir 2018, 34, 8443. doi: 10.1021/acs.langmuir.8b00889
Belmabkhout, Y.; Pillai, R. S.; Alezi, D.; Shekhah, O.; Bhatt, P. M.; Chen, Z.; Adil, K.; Vaesen, S.; Weireld, G. D.; Pang, M.; Suetin, M.; Cairns, A. J.; Solovyeva, V.; Shkurenko, A.; Tall, O. E.; Maurin, G.; Eddaoudi, M. J. Mater. Chem. A 2017, 5, 3293. doi: 10.1039/C6TA09406F
Bhatt, P. M.; Belmabkhout, Y.; Assen, A. H.; Weselinski, L. J.; Jiang, h.; Cadiau, A.; Xue, D. X.; Eddaoudi, M. Chem. Eng. J. 2017, 324, 392. doi: 10.1016/j.cej.2017.05.008
Li, J. R.; Sculley, J.; Zhou, H. C. Chem. Rev. 2012, 112, 869. doi: 10.1021/cr200190s
Wu, D.; Wang, C. C.; Liu, B.; Liu, D.; Yang, Q. Y.; Zhong, C. L. AIChE J. 2012, 58, 2078. doi: 10.1002/aic.v58.7
Liu, B.; Smit, B. J. Phys. Chem. C 2016, 114, 8515.
卞磊, 李炜, 魏振振, 刘晓威, 李松, 化学学报, 2018, 76, 303. doi: 10.3866/PKU.WHXB201708302Bian, L.; Li, W.; Wei, Z. Z.; Liu, X. W.; Li, S. Acta Chim. Sinica 2018, 76, 303. doi: 10.3866/PKU.WHXB201708302
Actintas, C.; Avci, G.; Daglar, H.; Gulcay, E.; Erucar, L.; Keskin, S. J. Mater. Chem. A 2018, 6, 5836. doi: 10.1039/C8TA01547C
许红, 童敏曼, 吴栋, 肖刚, 阳庆元, 刘大欢, 仲崇立, 物理化学学报, 2015, 31, 41. doi: 10.3866/PKU.WHXB201411132Xu, H.; Tong, M. M.; Wu, D.; Xiao, G.; Yang, Q. Y.; Liu, D. H.; Zhong, C. L. Acta Phys.-Chim. Sin. 2015, 31, 41. doi: 10.3866/PKU.WHXB201411132
Wilmer, C. E.; Farha, O. K.; Bae, Y. S.; Hupp, J. T.; Snurr, R. Q. Energy Environ. Sci. 2012, 5, 9849. doi: 10.1039/c2ee23201d
James, L. Technometrics 1983, 26, 415.
Qiao, Z. W.; Peng, C. W.; Zhou, J.; Jiang, J. W. J. Mater. Chem. A 2016, 4, 15904. doi: 10.1039/C6TA06262H
Qiao, Z. W.; Xu, Q.; Jiang, J. W. J. Membr. Sci. 2018, 551, 47. doi: 10.1016/j.memsci.2018.01.020
Glover, T. G.; Peterson, G. W.; Schindler, B. J.; Britt, D.; Yaghil, O. Chem. Eng. Sci. 2011, 66, 163. doi: 10.1016/j.ces.2010.10.002
Herm, Z. R.; Snisher, J. A.; Smit, B.; Krishna, R.; Long, J. R. J. Am. Chem. Soc. 2011, 133, 5664. doi: 10.1021/ja111411q
Chung, Y. G.; Camp, J.; Haranczy, M.; Sikora, B. J.; Bury, W.; Krungleviciute, V.; Yidirim, T.; Farha, O. K.; Sholl, D. S.; Snurr, R. Q. Chem. Mater. 2014, 26, 6185. doi: 10.1021/cm502594j
Willems, T. F.; Rycroft, C. H.; Kazi, M.; Meza, J. C.; Haranczyk, M. Microporous Mesoporous Mater. 2012, 149, 134. doi: 10.1016/j.micromeso.2011.08.020
Dubbeldam, D.; Ellis, D. E.; Snurr, R. Q. Mol. Simul. 2016, 42, 81. doi: 10.1080/08927022.2015.1010082
Wilmer, C. E.; Leaf, M.; Lee, C. Y.; Farha, O. K.; Hauser, B. G.; Hupp, J. T.; Snurr, R. Q. Nat. Chem. 2011, 4, 83.
吴选军, 郑佶, 李江, 蔡卫权, 物理化学学报, 2013, 29, 2207. doi: 10.3866/PKU.WHXB201307191Wu, X. J.; Zheng, J.; Li, J.; Cai, W. Q. Acta Phys.-Chim. Sin. 2013, 29, 2207. doi: 10.3866/PKU.WHXB201307191
Rappe, A. K.; Casewit, C. J.; Colwell, K. S. Ⅲ, W. A. G.; Skiff, W. M. J. Am. Chem. Soc. 1992, 114, 10024. doi: 10.1021/ja00051a040
Kadantsev, E. S.; Boyd, P. G.; Daff, T. D.; Woo, T. K. J. Phys. Chem. Lett. 2013, 4, 3056. doi: 10.1021/jz401479k
And, M. G.; Siepmann, J. I. J. Phys. Chem. B 1998, 102, 2569. doi: 10.1021/jp972543+
Qiao, Z. W.; Xu, Q.; Cheetham, A. K.; Jiang, J. W. J. Phys. Chem. C 2017, 121, 22208. doi: 10.1021/acs.jpcc.7b07758
Shah, M. S.; Tsapatsis, M.; Siepmann, J. I. J. Phys. Chem. B 2015, 119, 7041. doi: 10.1021/acs.jpcb.5b02536
Ewald, P. P. Ann. Phys. 1921, 369, 253. doi: 10.1002/(ISSN)1521-3889
图 2 四种描述符对TSN, SNM和TSC的关系
Figure 2 Relationship between the four descriptors and TSN, SNM and TSC
(a), (d), (g), (j) represent the effects of four descriptors on TSN; (b), (e), (h), (k) represent the effects of the four descriptors on SNM; (c), (f), (i), (l) represent the effects of the four descriptors on TSC. Red lines represent the linear fitting line. The figures contain the data of 4368 CoRE-MOFs with NH2S+CO2 > 1.2 mol/kg
表 1 Pearson相关系数和标准回归系数
Table 1. Pearson correlation coefficient and standardized regression coefficient
MOF descriptors |
Pearson correlation coefficient |
Standardized regression coefficient | |||
TSN | SNM | TSC | TSC | ||
LCD | -0.340 | -0.419 | -0.613 | -0.209 | |
ϕ | -0.393 | -0.569 | -0.717 | -0.327 | |
VSA | -0.379 | -0.530 | -0.673 | -0.113 | |
Qst0 | 0.480 | 0.434 | 0.536 | 0.368 |
表 2 20种吸附性能最优的MOFs
Table 2. 20 MOFs with best adsorption performance
No. | CSD code | TSC | NH2S/ (mol•kg-1) |
NCO2/ (mol•kg-1) |
SH2S+CO2/C1-C3 (105) |
Qst, H2S0/ (kJ•mol-1) |
Qst, CO20/ (kJ•mol-1) |
ϕ | VSA/ (m2•cm-3) |
LCD/ | Metal | Metal statea |
1 | KAPHOZ | 1.83 | 3.68 | 5.08 | 101.62 | 83.65 | 55.76 | 0.42 | 971.54 | 3.65 | K/Zn | NFI |
2 | IQUNAJ01 | 1.67 | 2.64 | 4.68 | 84.86 | 67.79 | 62.41 | 0.25 | 0 | 2.83 | K/Y | OMS |
3 | VAYDOQ | 1.66 | 1.83 | 5.41 | 84.00 | 101.47 | 100.25 | 0.28 | 42.45 | 4.11 | Li | SS |
4 | IQUNAJ | 1.66 | 2.49 | 4.73 | 83.75 | 66.71 | 61.90 | 0.24 | 0 | 2.83 | K/Y | OMS |
5 | WENSIS | 1.65 | 1.59 | 5.57 | 83.08 | 65.95 | 61.73 | 0.25 | 0 | 2.75 | K/Ca/Zr | OMS |
6 | YUJVOQ01 | 1.63 | 2.10 | 5.56 | 22.18 | 74.57 | 91.51 | 0.28 | 11.33 | 3.88 | Li/Zr | SS |
7 | YUJVOQ | 1.62 | 1.73 | 5.52 | 35.79 | 73.05 | 81.12 | 0.26 | 3.15 | 3.80 | Li/Zr | SS |
8 | YUJWAD | 1.60 | 1.51 | 5.41 | 52.77 | 71.08 | 78.67 | 0.24 | 0 | 3.70 | Li/Zr | SS |
9 | TUMGOX | 1.55 | 1.29 | 5.09 | 74.01 | 63.91 | 50.10 | 0.27 | 0 | 3.47 | Ce | SS |
10 | PAMGIV | 1.55 | 2.80 | 5.29 | 3.22 | 80.26 | 88.17 | 0.36 | 399.52 | 4.35 | Mn | SS |
11 | BEVQUP | 1.52 | 2.51 | 3.67 | 71.67 | 76.35 | 53.91 | 0.14 | 0 | 2.88 | Rb/Y | SS |
12 | SOPPOD | 1.44 | 2.16 | 3.79 | 28.74 | 66.65 | 62.41 | 0.25 | 0 | 2.85 | K/Tb | OMS |
13 | FASJAL | 1.43 | 0.00080 | 5.78 | 35.07 | 28.33 | 48.57 | 0.22 | 0 | 3.28 | Ca/Zr | SS |
14 | VIXGEQ | 1.43 | 0.69 | 6.17 | 3.10 | 62.85 | 93.35 | 0.27 | 266.29 | 5.15 | Zr | SS |
15 | QUQGAL | 1.42 | 4.88 | 8.94 | 0.030 | 98.99 | 79.14 | 0.47 | 893.42 | 5.24 | Li | OMS |
16 | IDAZEU | 1.42 | 1.65 | 3.83 | 63.49 | 69.95 | 66.39 | 0.25 | 0 | 2.74 | K/Yb | OMS |
17 | ZADDAJ | 1.38 | 1.47 | 5.28 | 1.88 | 77.26 | 93.13 | 0.27 | 0 | 3.58 | Li/Fe | SS |
18 | VIXGAM | 1.37 | 0.63 | 6.53 | 0.98 | 60.80 | 86.40 | 0.26 | 284.90 | 5.20 | Zr | SS |
19 | HIYMAF | 1.37 | 1.62 | 3.56 | 60.00 | 83.97 | 88.44 | 0.29 | 442.95 | 4.38 | K/Dy | OMS |
20 | VIXGAM01 | 1.35 | 0.63 | 6.50 | 0.78 | 62.47 | 85.05 | 0.27 | 286.20 | 5.21 | Zr | SS |
a Metal state includes open metal site (OMS), non-framework ion (NFI) and saturated state (SS). |