金屬-有機骨架材料用于去除天然氣中H2S的計算研究

許 紅 童敏曼 吳 棟 肖 剛 陽慶元 劉大歡,* 仲崇立(北京化工大學化學工程學院,先進納微結構材料實驗室,北京0009;中國石油吉林石化公司,吉林吉林30)

金屬-有機骨架材料用于去除天然氣中H2S的計算研究

許 紅1童敏曼1吳 棟1肖 剛2,*陽慶元1劉大歡1,*仲崇立1
(1北京化工大學化學工程學院,先進納微結構材料實驗室,北京100029;2中國石油吉林石化公司,吉林吉林132022)

眾所周知,天然氣作為一種利用效率高的清潔能源,其需求量正與日俱增.但天然氣中包含的H2S等有害氣體會危害人類健康、腐蝕設備、污染生態環境等.為解決這一問題,尋找良好的H2S吸附劑,本文采用巨正則系綜蒙特卡羅(GCMC)模擬方法,針對天然氣中H2S/CH4混合氣分離,對33種具有代表性的穩定金屬-有機骨架(MOF)材料進行H2S選擇性和工作容量(變壓吸附(PSA)及真空變壓吸附(VSA)過程)的篩選.結果表明,ZIF-80,Zn2-bpydtc,CAU-1-(OH)2,CH3O-MOFa最適用于本體系VSA過程的氣體分離;而后兩者最適用于PSA過程的氣體分離.通過分析高選擇性和高工作容量材料的結構特征,發現改性官能基團以及小孔作用的出現是影響選擇性的關鍵因素,其中―Cl、―OH、―OCH3基團對H2S氣體的吸附作用力最強.具有高的工作容量材料的特點是選擇性高,對氣體吸附作用力大,吸附位置多.基于篩選出的高選擇性、高工作容量的穩定MOF材料總結出的強化H2S選擇性及工作容量的一般性規律,為MOF材料應用于天然氣脫硫提供了理論基礎.

H2S;金屬-有機骨架;天然氣;篩選;分離

1 引言

近年來,由于其成本低廉、利用效率高、CO2等污染物排放量少,1天然氣被認為是重要的替代能源之一,已成為世界上需求增長率最快的能源.預計到2020年,天然氣的需求量將超過煤炭.2基于其形成過程,天然氣中除含有大量CH4外,還包含一定量的CO2和H2S等酸性氣體,3,4其中H2S的危害最為嚴重.5濃度為0.0003 mol·mol-1的H2S氣體在30 min內會導致人員窒息,混合氣中H2S濃度處于10-5mol· mol-1以下時才能保證其基本使用安全.6在天然氣輸送過程中,H2S的存在易導致管道等設備腐蝕,為保證天然氣的管道安全運輸要求,需將H2S濃度降到4×10-6mol·mol-1以下.4同時,從環境保護的角度來講,含有高濃度H2S天然氣的燃燒將極大地污染大氣,產生酸雨,危害生態環境.然而,即使不同地域的天然氣形成環境不同,H2S含量基本也在10-3mol· mol-1以上.7因此,無論是基于使用安全還是環境保護的角度,都有必要將H2S從天然氣中分離出去,大幅降低天然氣中H2S濃度.

常用的氣體分離方法包括吸收和吸附,其中吸收方法耗能較高且容易導致設備腐蝕.8相比而言,多孔材料的吸附具有操作條件溫和、設備簡單、對環境影響低等優點,被認為是一種有前途的可替代方法.其中,變壓吸附(PSA)是一種高效、節能、操作簡單的工業分離凈化方法.8它是在溫度一定的條件下變換壓力,通過高壓吸附和常壓脫附來達到分離氣體的目的.真空變壓吸附(VSA)可以被認為是一種特殊的變壓吸附過程,它與PSA的區別在于吸附壓力為常壓,脫附壓力為真空.變壓(真空)吸附單元操作的核心和難點,是選擇同時具有高選擇性和高工作容量的吸附劑.目前常用的吸附劑有沸石、活性炭、介孔硅、金屬-有機骨架(MOF)材料等多孔材料.其中,MOF材料是由金屬離子與有機配體經配位鍵結合而成的納米多孔材料,具有比表面積高、孔大小可調節性、結構多樣性、較好的化學穩定性和熱穩定性等特點,是目前化工和材料領域研究的熱點之一.9-15MOF材料結構具有特殊性,在作為吸附劑時與常規吸附劑相比具備優勢,尤其是在氣體吸附分離方面:基于組合化學思想,很容易在分子層次上對MOF材料的孔道化學特性和孔尺寸進行設計和精確調控,進而調控其氣體吸附分離等性能.12,13因此,研究MOF材料用作天然氣中H2S的分離劑在實際工業操作過程中具有重要的意義.材料的穩定性在工業應用中至關重要,尤其是熱穩定性和化學穩定性,但某些MOF材料的性能并不夠穩定,比如MOF-5,14UiO-67,UiO-68-NH215等在水的存在下會導致結構坍塌,故此類材料不適用于含水的混合氣體吸附分離.而吸附分離過程中常常是在一定的濕度條件下進行的,所以有必要選取穩定的MOF材料來進行吸附分離.前人在研究混合氣體吸附分離時(尤其是在大規模材料篩選中),大都沒有考慮材料的穩定性,使材料在應用方面受到了限制.為克服這一限制,我們基于大量的MOF材料文獻,13,14選擇了33種性能穩定的MOF材料作為研究對象,確保了材料在實際操作中的可循環使用性.在研究方法的選擇上,考慮到采用實驗方法分離H2S/CH4混合氣體時,成本較高且耗時較長,而分子模擬方法可以高效且快速地對大量材料某一性能進行篩選,節省開支和時間,并且在新方法的開發上有很大優勢,已經成為材料科學領域中一個強大的研究工具.16因此,本研究采用分子模擬方法來對MOF材料的H2S吸附分離性能進行篩選.

目前,已經報道了很多關于天然氣中分離H2S的實驗工作,3,6,17,18為我們提供了很多有價值的參考信息.然而,這些研究只是針對某一種或者幾種沸石和金屬-有機骨架材料進行的,而本工作選擇較大的樣本數,篩選H2S/CH4混合氣體中H2S選擇性和工作容量高的材料,可以獲得更為準確的性能研究結果.

2 研究模型和方法

2.1 材料結構

作為研究對象的33種穩定MOF材料的結構均以其實驗XRD所得的晶體結構為基礎,在提供的數據基礎上,采用Materials Studio軟件19進行結構搭建并進行部分材料結構性質的計算.材料包括10種ZIFs:ZIF-3,20ZIF-8,20,21ZIF-68,22ZIF-69,22ZIF-70,22ZIF-78,23ZIF-79,23ZIF-80,23ZIF-81,23ZIF-9024;6種UiO-66-Zr及其改性系列材料:UiO-66-Zr,25UiO-66-2CH3,26,27UiO-66-2COOH,8UiO-66-Br,28UiO-66-NH2,28Hydrad-UiO-66-NO2;287種MILs:MIL-47(V),29MIL-68-Al,30MIL-120,31MIL-140A,B,C,D;32以及其它類型水穩定材料:Al(OH)(1,4-NDC),33Al-PMOF,34CAU-10-H,35CAU-1-(OH)2,36Zr-fumarate,37MOF-525,38CoNIm-RHO,39NOTT-300,40CH3O-MOFa,41Zn2-bpydtc.42在化學成分方面,這些材料涉及Zn(12種MOFs)、Al(8種MOFs)、Zr(12種MOFs)、Co(1種MOFs)四種金屬.從有機配體方面來講,包含羧酸類配體,咪唑類配體,及其各種官能團改性配體(―NH2、―NO2、―CHO、―CH3、―Br、―Cl、―COOH、―OCH3)和卟啉結構的配體.本文中選擇的這些材料包括多種常見的MOF材料,具有很好的代表性.材料結構基本參數表見Supporting Information中表S1.

2.2 材料電荷

在分子模擬過程中,作為主體的MOF材料的骨架與作為客體的流體分子之間的相互作用力由范德華力和庫侖靜電力來描述.準確描述MOF骨架原子的配分電荷才能夠準確地描述流體分子在材料內部的吸附、分離等行為.本實驗室開發的CBAC (connectivity-based atom contribution)計算方法,是一種快速估計MOF材料電荷的有效手段,43-47由此方法建立的CBAC電荷庫已經被成功應用于MOF材料吸附選擇性能的模擬計算中.48-50本文中材料的骨架電荷均取自CBAC數據庫.51-53對于CBAC庫中未包含的CAU-1-(OH)2中的O3和C5的電荷,采用量子化學計算得到,所用的金屬簇、原子標號圖和電荷表見Supporting Information中圖S1和表S2.電荷計算方法采用密度泛函理論(DFT)方法.所有的計算由Gaussian 03軟件54完成.在計算過程中,應用非限制性的B3LYP泛函理論,采用ChelpG的方法得到材料中骨架原子的電荷分布.51在此過程中,對金屬原子采用LANL2DZ的基組,而對其它原子采用6-31+G*的基組.

2.3 力場參數

在模擬過程中,對于流體分子,CH4分子被視為無電荷的單一作用點,使用的力場為TraPPE力場,由Lennard-Jones(L-J)勢能模型描述.55將H2S分子視為三點模型,用Kamath等56開發的Kamathpotoff力場描述.S原子是H2S分子唯一的L-J相互作用點,局部電荷分布在各個原子上.H―S鍵長為0.134 nm,H―S―H角度為92.5°.材料與氣體之間的相互作用由L-J勢能模型來描述,MOF骨架原子所采用的力場均為Dreiding力場,57力場參數見表1.對于Dreiding力場缺失的原子類型,如V、Cu、Zr等,其力場參數從UFF力場58中獲得.上述客體分子模型和Dreiding力場已經被廣泛用于MOF材料的包括H2S/CH4在內的氣體吸附分離.30

2.4 模擬方法

在模擬方法上,本文采用常用的巨正則蒙特卡羅(GCMC)模擬方法來計算混合氣體的吸附選擇性.模擬溫度為303 K.H2S/CH4混合氣中H2S濃度為10-3mol·mol-1.采用壓力為0.01-1 MPa.PSA過程的吸附壓力為1 MPa,脫附壓力為0.1 MPa;VSA過程的吸附壓力為0.1 MPa,脫附壓力為0.01 MPa.所有的MOF材料骨架均被視為剛性.L-J相互作用參數的截斷半徑為1.4 nm,長程靜電力作用均采用以此為邊界條件的Ewald加和法進行處理.在模擬過程中,模擬盒子所包含的晶胞的數目因MOF材料種類不同而不同,范圍變化從1×1×1到8×8×8,以保證計算效率和模擬的準確性.對每一個狀態點,GCMC模擬均運行1×107步來保證體系平衡,然后運行1× 107步來統計所需的熱力學性質.另外,某些MOF材料中具有客體分子(CH4和H2S)不能進入的GME和HPR類型59的孔道,如ZIF-68、ZIF-69、ZIF-70、ZIF-78、ZIF-79、ZIF-80、ZIF-81,因此,我們在分子模擬過程中采用了堵孔策略,即以GME和HPR的孔中心為球心,取適宜的球半徑將孔道堵塞,以保證模擬的準確性.這種堵孔策略已經廣泛應用于沸石和MOF材料領域的氣體吸附分離,44,60包括本組以前的研究工作.30

表1 本文中H2S、CH4和MOFs的L-J勢能參數Table 1 Lennard-Jones(L-J)potential parameters for H2S,CH4,and MOFs used in this work

組分A對組分B的選擇性(S)用公式S=(xA/xB)/(yB/yA)計算.xA和xB分別代表體相中組分A和組分B的摩爾百分數,本文中xA/xB即某壓力下被MOF材料所吸附的H2S與CH4的摩爾比;yA和yB分別代表組分A和組分B在吸附相中的摩爾百分數,本文中yB/yA即CH4與H2S的初始摩爾比.材料的工作容量是指:單位質量的MOF材料在吸附與脫附壓力下對目標氣體H2S的吸附量之差.

3 結果與討論

3.1 力場驗證

力場的可靠性在分子模擬中非常重要,因此,本文對H2S和CH4的模型和力場進行了驗證.從圖1的H2S和CH4的吸附等溫線中可以看出,H2S和CH4的吸附量模擬值均與已有的實驗值3,28,30,61,62吻合,表明本工作所用力場以及模擬方法是可靠的.

圖1(a1,a3)低壓區模擬值偏高的原因,可能是由于通過實驗得到的材料本身存在活化不夠好或者晶體缺陷等問題,并不是完美晶體,而模擬是針對完美晶體進行預測,所以產生一定的誤差.這個誤差是在可接受范圍內的.28,63,64

圖1 H2S和CH4的吸附等溫線Fig.1 Adsorption isotherms of H2S and CH4

3.2 材料的選擇性

材料對氣體的選擇性是PSA及VSA過程有效分離氣體的關鍵.一般情況下,選擇性會受到操作壓力的影響,因此,本文針對天然氣中H2S/CH4混合氣分離,計算了混合氣體在0.01-1 MPa時的選擇性,研究此時選擇性與操作壓力的關系.總體上,選擇性隨著壓力的變化趨勢基本可以分為兩種,一種是隨著壓力升高呈下降然后趨于平緩的趨勢,如 ZIF-80(圖2(a1));另外一種是隨著壓力升高呈先上升然后基本平緩的趨勢,如UiO-66-2COOH(圖2 (a2)).為了清晰起見,此處只給出典型材料的選擇性-壓力關系圖,其余材料數據見Supporting Information中圖S2.由圖2(a1)可以看出,低壓下ZIF-80具有很高的選擇性,但隨著壓力的升高,選擇性迅速下降,這種現象的產生是由于H2S吸附熱和CH4吸附熱之差的下降.對比圖1(a)與圖1(b)可知,材料選擇性隨壓力的變化趨勢與兩種氣體的吸附熱差值變化一致:吸附熱增加,材料的吸附分離能力隨之變大;反之,吸附分離能力下降.這種吸附熱和選擇性的關系與在其他氣體吸附分離中類似.64-66從吸附的角度來說,ZIF-80選擇性迅速下降是由于H2S吸附量的增加速度遠低于CH4吸附量的增加速度. MOF材料分離混合氣時的作用力包括靜電力和范德華力,H2S有很強的偶極矩,而CH4沒有偶極矩.因此,H2S既有范德華作用,也存在靜電作用;而CH4分子只有范德華作用比較明顯.低壓時,靜電力作用占主導地位,基團對H2S的靜電作用使材料吸附了大量的H2S氣體(0.1 MPa下的質心分布圖見圖3(a)).當壓力升高時,H2S吸附量趨于飽和,客體分子聚集在功能基團周圍,并使靜電力的影響逐漸下降,壓力的作用逐漸增強.由于CH4與H2S的分子體積相近,分子體積對吸附量的影響差異很小,而CH4在混合氣中所占比例很大(0.999),故在壓力的推動作用下,材料開始大量吸附CH4氣體,使得CH4吸附量的增長速度要比H2S氣體快很多,從而導致ZIF-80的選擇性在低壓階段下降非常劇烈.待壓力升高到一定程度后,CH4的吸附量也趨于飽和,于是隨著壓力的繼續升高,選擇性下降程度逐漸減小,曲線開始趨于平緩.圖4為ZIF-80材料CH4與H2S吸附量隨壓力變化曲線.由圖可以看到,低壓時H2S吸附量較高,壓力升高后CH4吸附量迅速增加.

圖2 ZIF-80、UiO-66-2COOH選擇性(a1,a2)和吸附熱差(b1,b2)圖Fig.2 Selectivity(a1,a2)and adsorption enthalpy(b1,b2)curves of ZIF-80 and UiO-66-2COOH

為驗證基團靜電力對極性分子的重要作用,我們將ZIF-80中的―Cl替換成―H,在其它模擬條件不變的情況下重新計算材料的選擇性,結果表明,替換之后材料的選擇性遠低于ZIF-80(前后兩種材料的選擇性對比、ZIF-80的結構圖和前后電荷對比表見Supporting Information中圖S3、圖S4和表S3).類似地,對比ZIF-68和ZIF-69的結構和選擇性關系時可以發現(兩材料選擇性見Supporting Information中圖S2(a)),ZIF-69的結構只是將ZIF-68結構中某位置的―H替換成為―Cl,此時ZIF-69的選擇性也明顯高于ZIF-68.由此可見,基團改性的確對于提高材料選擇性具有重要作用,同時說明―Cl基團對H2S這種極性分子的靜電作用可以提高材料對于H2S/CH4混合氣體中H2S的選擇性.UiO-66-2COOH的選擇性先上升再趨于平緩,低壓下其H2S吸附量隨著壓力的增加速度比CH4快,壓力升高到一定程度后,兩種氣體的吸附量均趨于飽和,于是選擇性隨壓力的升高基本不發生改變.

計算得到33種穩定材料的選擇性后,我們從中篩選出一部分選擇性高的材料.圖5中列出的是0.1和1 MPa下篩選出的選擇性最高的十種材料;并通過質心分布圖(圖3)對它們的結構特點進行分析.由圖5可知,選擇性較高的材料中包括三種ZIFs:ZIF-78,79,80.它們的吸附位置均環繞在官能團周圍,形成一個圓形.以ZIF-80為例,如圖3(a)所示,H2S基本全部被吸附在―Cl基團周圍.由CAU-1-(OH)2質心分布圖可知,被吸附的H2S分子排列于兩―OH的中間位置,形成一條直線,與金屬Al基本在同一平面上,如圖3(b)所示.CAU-1-(OH)2含有―OH和―CH3兩種基團,故分析是―OH和―CH3兩種基團的共同作用使CAU-1-(OH)2選擇性較高.為了確定哪種基團吸附作用更強,我們作出了材料的不同種類原子之間的徑向分布函數圖(RDF,見Supporting Information中圖S5),發現H2S距―OH基團的距離小于距―CH3基團距離.因此,CAU-1-(OH)2中的―OH基團是吸附H2S時起主要作用的基團,―OH對H2S的吸附作用比―CH3更強.CH3O-MOFa吸附位置均圍繞于甲氧基基團周圍,如圖3(c)所示.將材料中―OCH3基團替換為―H后,選擇性明顯降低(前后選擇性對比見Supporting Information中圖S6),故―OCH3基團改性可以明顯提高材料對H2S氣體的選擇性.官能團的加入使材料孔體積有所減小,但根據質心分布圖中的吸附位置,可以確定提高選擇性更重要的原因是―OCH3對H2S的強吸附作用.Al-PMOF吸附位置基本在兩卟啉平面中間(夾層吸附),如圖3(d)所示.圖5中其它材料的質心分布圖見Supporting Information中圖S7.

圖3 MOF材料在0.1 MPa下吸附H2S的質心分布圖Fig.3 Probability distributions of mass center for MOFs with adsorbed H2S at 0.1 MPa

圖4 ZIF-80材料的CH4與H2S吸附量Fig.4 Adsorption capacities of CH4and H2S for ZIF-80 material

綜上,除材料本身結構特征之外,官能團以及小孔效應是影響選擇性的關鍵因素.材料基本拓撲結構相同時,能夠更直觀地看出官能團和孔尺寸的影響.因此,我們進一步分析了UiO-66系列和具有相同拓撲結構的ZIFs的選擇性.

由UiO-66系列的選擇性可以看出官能團影響的重要性.官能團的加入使材料對客體分子的作用力增強,同時也縮小了孔體積.如圖6(a)所示,改性材料UiO-66-Zr-Xn(Xn=―2COOH,―2CH3,―NO2,―NH2,―Br)的選擇性均高于UiO-66-Zr的選擇性.其中UiO-66-2COOH選擇性居于首位.材料所含官能團極性順序為:―COOH＞―NO2＞―NH2＞―Br＞―CH3,由此可見,材料選擇性順序與官能團極性順序基本相同(―CH3除外).―CH3的特殊性可能是由于兩個―CH3基團占用了較多的材料空隙,使孔徑縮小較多,于是選擇性偏高.材料孔徑大小順序為: UiO-66-2COOH＜UiO-66-2CH3＜UiO-66-NO2≈UiO-66-NH2＜UiO-66-Br＜UiO-66-Zr,與選擇性順序恰好相反.本組以前的研究工作28中也研究了UiO-66-Zr改性材料對CO2/CH4混合氣體的吸附選擇性,其結果也是―COOH改性的選擇性最高.圖6(b)中ZIFs具有相同的拓撲結構,只是官能團有所差異.ZIF-78吸附孔中攜帶的是―NO2基團,ZIF-79為―CH3, ZIF-81為―Br,ZIF-69為―Cl,ZIF-68沒有特殊官能團.此時選擇性大小順序為:ZIF-78(―NO2)＞ZIF-79 (―CH3)＞ZIF-81(―Br)≈ZIF-69(―Cl)＞ZIF-68.基本與材料所帶官能團極性強弱順序相同:―NO2＞―Br≈―Cl＞―CH3＞―H(―CH3除外);但與材料孔徑66大小順序相反:ZIF-78＜ZIF-79＜ZIF-81≈ZIF-69＜ZIF-68.因此,在具有相似拓撲結構的系列材料中,官能團改性和孔徑是影響材料選擇性的主要因素:官能團極性越強,材料孔徑越小,選擇性越高.

圖5 (a)0.1 MPa和(b)1 MPa下選擇性(S)最高的十種MOF材料Fig.5 The top 10 MOFs screened based on selectivity(S) at(a)0.1 MPa and(b)1 MPa

圖6 0.1 MPa下拓撲結構相同MOFs的選擇性對比Fig.6 Selectivity comparison of MOFs with same topology structures at 0.1 MPa

表2給出了本文材料選擇性與文獻中研究分離H2S/CH4混合氣體中H2S材料選擇性的對比.其中, ZIF-80和CAU-1-(OH)2是本研究中選用的MOF材料,Si-CHA、NaY、NC100AC是文獻中的沸石材料(Si-CHA的選擇性是理想選擇性,沒有給出H2S的具體濃度).通過表2的對比表明,ZIF-80和CAU-1-(OH)2在分離H2S/CH4混合氣體中H2S時的選擇性較文獻中沸石材料有明顯的優勢.

表2 本文中MOFs與文獻中沸石選擇性對比Table 2 Selectivity comparison of the MOFs used in this work and zeolites in the literature

3.3 工作容量

在實際應用中,材料除了應具有穩定的性能和優良的選擇性,還需要具備較大的處理量和良好的再生性能.工作容量正是反映材料處理量大小和可再生性能優劣的性能參數.因此,本研究在基于選擇性對穩定材料進行篩選之后,又分別基于VSA、PSA兩種過程的工作容量對材料進行了篩選.工作容量由吸附和脫附兩個壓力下的吸附量相減得到,單一壓力下的吸附量大不能代表材料的工作容量高,只有吸附和脫附的吸附量差值大時,材料才具備高工作容量.吸附量在低壓時主要受作用力的影響,高壓時則主要受孔體積的制約,有特殊官能團的材料,對H2S有更強的吸附作用,故而低壓下吸附量更高.因此,作用力和孔體積作為影響吸附量的關鍵因素,間接地影響著材料的工作容量.

圖7中給出了VSA和PSA兩種過程下工作容量最高的十種材料.與圖4對比可知,0.1 MPa和1 MPa下選擇性最高的10種材料,工作容量基本也最高(除ZIF-790.1MPa/1MPa、MOF-5250.1MPa和UiO-66-Br1MPa),但相對大小有所變化.其中,最適用于VSA過程氣體分離的材料有ZIF-80、Zn2-bpydtc、CAU-1-(OH)2、CH3O-MOFa等,最適用于PSA過程氣體分離的材料有CH3O-MOFa和CAU-1-(OH)2等.VSA過程中,ZIF-80的選擇性非常高,工作容量也最高,但它在PSA過程的工作容量并不占優勢(見圖7).這是由于壓力升高時,ZIF-80的吸附位趨于飽和,H2S氣體的吸附量增加得很少,導致其1 MPa下工作容量比0.1 MPa時提高較少,于是PSA過程的工作容量不高.由此可見,像ZIF-80這種選擇性隨著壓力發生驟降,高壓下吸附位置趨于飽和的材料,更適用于VSA過程的H2S吸附.CAU-1-(OH)2、CH3OMOFa、UiO-66-2COOH、MIL-47(V)這些材料在壓力升高時,H2S的選擇性也升高,壓力升高時可供吸附的位置多,對H2S的吸附量增加較多,工作容量必然上升較快,它們在PSA過程的工作容量較VSA過程增加的就比較明顯.故當材料的選擇性隨著壓力升高而升高時,材料在PSA過程的工作容量比較高.總之,單純選擇性高的材料工作容量不一定高,還需要考慮選擇性隨壓力變化的趨勢.工作容量高的材料具有的特點是,材料對氣體的吸附作用力大、吸附位置多,壓力升高時,對目標氣體選擇性呈上升趨勢,吸附量增加明顯.

圖7 VSA(a)和PSA(b)過程中工作容量最高的十種MOF材料Fig.7 The top 10 MOFs based on the highest working capacity in VSA(a)and PSA(b)processes

4 結論

本文研究了33種具有代表性的穩定MOF材料對H2S/CH4混合體系的分離性質.結果表明,官能團改性和小孔作用的出現是提高選擇性的關鍵.材料各基團中―Cl、―OH、―OCH3基團對H2S氣體的吸附作用力最強,在今后實際吸附分離H2S氣體時,可以嘗試選擇帶有這幾種官能團且孔尺寸比較小的材料.選擇性高,對氣體吸附作用力大,吸附位置多的材料工作容量一般比較高.單純選擇性高的材料工作容量不一定高,如PSA過程的ZIF-80.此外,通過計算篩選出最適用于VSA過程氣體分離的材料有ZIF-80、Zn2-bpydtc、CAU-1-(OH)2、CH3O-MOFa等;最適用于PSA過程氣體分離的材料為CH3OMOFa和CAU-1-(OH)2等.這些結論對于今后研究從H2S/CH4混合氣體分離H2S有一定的指導意義,為MOF材料應用于天然氣脫硫提供了一定理論依據.

Supporting Information:Material structure parameters; the cluster of CAU-1-(OH)2and the structure of CAU-1-(OH)2; atomic partial charges of CAU-1-(OH)2;selectivity of MOFs; selectivity comparison of ZIF-80 and its modified material; structure of ZIF-80;atomic partial charge comparison of ZIF-80 and its modified material;the RDF of CAU-1-(OH)2;selectivity comparison of CH3O-MOFa and its modified material; mass center distribution map of H2S are included.This information is available free of charge via the internet at http://www. whxb.pku.edu.cn.

References

(1) Shahbaz,M.;Lean,H.H.;Farooq,A.Renew.Sust.Energ.Rev.2013,18,87.doi:10.1016/j.rser.2012.09.029

(2) Tagliabue,M.;Farrusseng,D.;Valencia,S.;Aguado,S.;Ravon, U.;Rizzo,C.;Corma,A.;Mirodatos,C.Chem.Eng.J.2009,155,553.doi:10.1016/j.cej.2009.09.010

(3) Vaesen,S.;Guillerm,V.;Yang,Q.Y.;Wiersum,A.D.; Marszalek,B.;Gil,B.;Vimont,A.;Daturi,T.;Llewellyn,P.L.; Serre,C.;Maurin,G.;Weireld,G.D.Chem.Commun.2013,49, 10082.doi:10.1039/c3cc45828h

(4) Belmabkhout,Y.;Heymans,N.;Weireld,G.D.;Sayari,A. Energy&Fuels2011,25,1310.doi:10.1021/ef1015704

(5) Guidotti,T.L.Occupational Medicine-Oxford1996,46, 367.doi:10.1093/occmed/46.5.367

(6) Cosoli,P.;Ferrone,M.;Pricl,S.;Fermeglia,M.Chem.Eng.J.2008,145,86.doi:10.1016/j.cej.2008.07.034

(7) Bhide,B.D.;Voskericyan,A.;Stern,S.A.J.Membr.Sci.1998,140,27.doi:10.1016/S0376-7388(97)00257-3

(8)Yang,Q.Y.;Vaesen,S.;Ragon,F.;Wiersum,A.D.;Wu,D.; Lago,A.;Devic,T.;Martineau,C.;Taulelle,F.;Llewellyn,P.L.; Jobic,H.;Zhong,C.L.;Serre,C.;Weireld,G.D.;Maurin,G.Angew.Chem.2013,125,10506.doi:10.1002/ange.201302682

(9) An,X.H.;Liu,D.H.;Zhong,C.L.Acta Phys.-Chim.Sin.2011,27,553.[安曉輝,劉大歡,仲崇立.物理化學學報,2011,27, 553.]doi:10.3866/PKU.WHXB20110319

(10) Mu,W.;Liu,D.H.;Yang,Q.Y.;Zhong,C.L.Acta Phys.-Chim. Sin.2010,26,1657. [穆 韡,劉大歡,陽慶元,仲崇立.物理化學學報,2010,26,1657.]doi:10.3866/PKU.WHXB20100616

(11) Zhong,C.L.;Liu,D.H.;Yang,Q.Y.The Structure-Activity Relationship and Design of Metal-Organic Frameworks;Science Press:Beijing,2013;pp 1-3. [仲崇立,劉大歡,陽慶元.金屬-有機骨架材料的構效關系及設計.北京:科學出版社,2013:1-3.]

(12) Lu,Z.;Xing,H.;Sun,R.;Bai,J.;Zheng,B.;Li,Y.Cryst. Growth Des.2012,12,1081.doi:10.1021/cg201594p

(13) Zhai,F.;Zheng,Q.;Chen,Z.CrystEngComm2013,15,2040. doi:10.1039/c2ce26701b

(14)Yang,J.;Grzech,A.;Mulder,F.M.;Dingemans,T.J.Chem. Commun.2011,47,5244.doi:10.1039/c1cc11054c

(15) Schaate,A.;Roy,P.;Godt,A.;Lippke,J.;Waltz,F.;Wiebcke, M.;Behrens,P.Chem.Eur.J.2011,17,6643.doi:10.1002/ chem.v17.24

(16) Ma,Q.T.;Xiao,Y.L.;Tong,M.M.;Li,Z.J.;Yang,Q.Y.;Liu, D.H.;Zhong,C.L.Scientia Sinica Chimica2014,44(6),812. [麻沁甜,肖遠龍,童敏曼,李正杰,陽慶元,劉大歡,仲崇立.中國科學:化學,2014,44(6),812.]

(17) Maghsoudi,H.;Soltanieh,M.;Bozorgzadeh,H.; Mohamadalizadeh,A.Adsorption2013,19,1045.doi:10.1007/ s10450-013-9528-1

(18) Hamon,L.;Frère,M.;Weireld,G.Adsorption2008,14,493. doi:10.1007/s10450-008-9113-1

(19)Accelrys,Materials studio,Version 3.0;Accelrys Software Inc: San Diego,2003.

(20) Park,K.S.;Ni,Z.;Cote,A.P.;Choi,J.Y.;Huang,R.;Uribe-Romo,F.J.;Chae,H.K.;O'Keeffe,M.;Yaghi,O.M.Proc. Natl.Acad.Sci.U.S.A.2006,103,10186.doi:10.1073/ pnas.0602439103

(21) Shi,Q.;Chen,Z.F.;Song,Z.W.;Li,J.P.;Dong,J.X.Angew. Chem.2011,123,698.doi:10.1002/ange.v123.3

(22) Banerjee,R.;Phan,A.;Wang,B.;Knobler,C.;Furukawa,H.; O'Keeffe,M.;Yaghi,O.M.Science2008,319,939.

(23) Banerjee,R.;Furukawa,H.;Britt,D.;Knobler,C.;O'Keeffe, M.;Yaghi,O.M.J.Am.Chem.Soc.2009,131,3875.doi: 10.1021/ja809459e

(24) Morris,W.;Doonan,C.J.;Furukawa,H.;Banerjee,R.;Yaghi, O.M.J.Am.Chem.Soc.2008,130,12626.doi:10.1021/ ja805222x

(25) Cavka,J.H.;Jakobsen,S.;Olsbye,U.;Guillou,N.;Lamberti, C.;Bordiga,S.;Lillerud,K.P.J.Am.Chem.Soc.2008,130, 13850.doi:10.1021/ja8057953

(26) Huang,Y.;Qin,W.;Li,Z.;Li,Y.Dalton Trans.2012,41,9283. doi:10.1039/c2dt30950e

(27) Jasuja,H.;Walton,K.S.J.Phys.Chem.C2013,117,7062.doi: 10.1021/jp311857e

(28)Yang,Q.Y.;Wiersum,A.D.;Llewellyn,P.L.;Guillerm,V; Serre,C.;Maurin,G.Chem.Commun.2011,47,9603.doi: 10.1039/c1cc13543k

(29) Barthelet,K.;Marrot,J.;Riou,D.;Ferey,G.Angew.Chem.2002,114,291.

(30) Yang,Q.Y.;Vaesen,S.;Vishnuvarthan,M.;Ragon,F.;Serre,C.; Vimont,A.;Daturi,M.;De Weireld,G.;Maurin,G.J.Mater. Chem.2012,22,10210.doi:10.1039/c2jm15609a

(31) Volkringer,C.;Loiseau,T.;Haouas,M.;Taulelle,F.;Popov,D.; Burghammer,M.;Riekel,C.;Zlotea,C.;Cuevas,F.;Latroche, M.;Phanon,D.;Knofelv,C.;Llewellyn,P.L.;Férey,G.Chem. Mater.2009,21,5783.doi:10.1021/cm9023106

(32) Guillerm,V.;Ragon,F.;Dan-Hardi,M.;Devic,T.; Vishnuvarthan,M.;Campo,B.;Vimont,A.;Clet,G.;Yang,Q.; Maurin,G.;Férey,G.;Vittadini,A.;Gross,S.;Serre,C.Angew. Chem.Int.Edit.2012,51,9267.doi:10.1002/anie.201204806

(33) Comotti,A.;Bracco,S.;Sozzani,P.;Horike,S.;Matsuda,R.; Chen,J.;Takata,M.;Kubota,Y.;Kitagawa,S.J.Am.Chem. Soc.2008,130,13664.doi:10.1021/ja802589u

(34) Fateeva,A.;Chater,P.A.;Ireland,C.P.;Tahir,A.A.;Khimyak, Y.Z.;Wiper,P.V.;Darwent,G.R.;Rosseinsky,M.G.Angew. Chem.2012,124,7558.doi:10.1002/ange.v124.30

(35) Reinsch,H.;van der Veen,M.A.;Gil,B.;Marszalek,B.; Verbiest,T.;de Vos,D.;Stock,N.Chem.Mater.2013,25, 17.doi:10.1021/cm3025445

(36) Ahnfeldt,T.;Guillou N.;Gunzelmann D.;Margiolaki,I.; Loiseau,T.;Ferey,G.;Senker,J.;Stock,N.Angew.Chem.2009,121,5265.doi:10.1002/ange.v121:28

(37) Wi?mann,G.;Schaate,A.;Lilienthal,S.;Bremer,I.;Schneider, A.M.;Behrens,P.Microporous Mesoporous Mat.2012,152, 64.doi:10.1016/j.micromeso.2011.12.010

(38) Morris,W.;Volosskiy,B.;Demir,S.;Gándara,F.;McGrier,P. L.;Furukawa,H.;Cascio,D.;Stoddart,J.F.;Yaghi,O.M.Inorg.Chem.2012,51,6443.

(39) Biswal,B.P.;Panda,T.;Banerjee,R.Chem.Commun.2012,11868.

(40) Yang,S.;Sun,J.;Ramirez-Cuesta,A.J.;Callear,S.K.;David, W.I.F.;Anderson,D.P.;Newby,R.;Blake,A.J.;Parker,J.E.; Tang,C.C.;Schr?der,M.Nat.Chem.2012,4,887.doi: 10.1038/nchem.1457

(41)Weston,M.H.;Delaquil,A.A.;Sarjeant,A.A.;Farha,O.K.; Hupp,J.T.;Nguyen,S.T.Cryst.Growth Des.2013,13,2938. doi:10.1021/cg400342m

(42) Lin,X.;Blake,A.J.;Wilson,C.;Sun,X.Z.;Champness,N.R.; George,M.W.;Hubberstey,P.;Mokaya,R.;Schr?der,M.J.Am.Chem.Soc.2006,128,10745.doi:10.1021/ja060946u

(43) Meek,S.T.;Greathouse,J.A.;Allendorf,M.D.Adv.Mater.2011,23,249.doi:10.1002/adma.201002854

(44) Krishna,R.;van Baten,J.M.Phys.Chem.Chem.Phys.2011,13,10593.doi:10.1039/c1cp20282k

(45) Li,J.;Ma,Y.;McCarthy,M.C.;Sculley,J.;Yu,J.;Jeong,H.; Balbuena,P.B.;Zhou,H.Coord.Chem.Rev.2011,255,1791. doi:10.1016/j.ccr.2011.02.012

(46)Wilmer,C.E.;Leaf,M.;Lee,C.Y.;Farha,O.K.;Hauser,B.G.; Hupp,J.T.;Snurr,R.Q.Nature Chemistry2011,4,83.

(47) Yang,Q.Y.;Liu,D.H.;Zhong,C.L.;Li,J.R.Chem.Rev.2013,113,8261.

(48) Keskin,S.J.Phys.Chem.C2011,115,800.doi:10.1021/ jp109743e

(49) Krishna,R.;van Baten,J.M.Langmuir2010,26,10854.doi: 10.1021/la100737c

(50) Lin,L.;Berger,A.H.;Martin,R.L.;Kim,J.;Swisher,J.A.; Jariwala,K.;Rycroft,C.H.;Bhown,A.S.;Deem,M.W.; Haranczyk,M.;Smit,B.Nature Materials2012,11,633.doi: 10.1038/nmat3336

(51) Xu,Q.;Zhong,C.L.J.Phys.Chem.C2010,114,5035.

(52) Zheng,C.C.;Zhong,C.J.Phys.Chem.C2010,114,9945.doi: 10.1021/jp102409a

(53) Xiao,Y.L.;Yang,Q.Y.;Liu,D.H.;Zhong,C.L.CrystEngComm2013,15,9588.doi:10.1039/c3ce41081a

(54) Frisch,M.J.;Trucks,G.W.;Schlegel,H.B.;et al.Gaussian 03,Revision D.01;Gaussian Inc.:Wallingford,CT,2004.

(55) Martin,M.G.;Siepmann,J.I.J.Phys.Chem.B1998,102,2569.doi:10.1021/jp972543+

(56) Kamath,G.;Lubna,N.;Potoff,J.J.J.Phys.Chem.2005,123,124505.doi:10.1063/1.2049278

(57) Mayo,S.L.;Olafson,B.D.;Goddard,W.A.,II,J.Phys.Chem.1990,94,8897.doi:10.1021/j100389a010

(58) Rappe,A.K.;Casewit,C.J.;Colwell,K.S.;Goddard,W.A.; Skiff,W.M.J.Am.Chem.Soc.1992,114,10024.doi:10.1021/ ja00051a040

(59) Babarao,R.;Dai,S.;Jiang,D.J.Phys.Chem.C2011,115,8126.

(60) Krishna,R.;van Baten,J.M.Langmuir2010,26,2975.doi: 10.1021/la9041875

(61) Hamon,L.;Serre,C.;Devic,T.;Loiseau,T.;Millange,F.;Férey, G.;Weireld,G.D.J.Am.Chem.Soc.2009,131,8775.doi: 10.1021/ja901587t

(62) Meng,Z.S.;Lu,R.F.;Rao,D.W.;Kan,E.J.;Xiao,C.Y.; Deng,K.M.Int.J.Hydrog.Energy2013,38,9811.doi: 10.1016/j.ijhydene.2013.05.140

(63) Nilton,R.,Jr.;Herve,J.;Aziz,G.;Fabrice,S.;Guillaume,M.; Sandrine,B.;Philip,L.L.;Thomas,D.;Christian,S.;Gerard,F.Angew.Chem.Int.Edit.2008,47,6611 doi:10.1002/anie.v47:35

(64) Zhang,W.J.;Huang,H.L.;Zhong,C.L.;Liu,D.H.Phys. Chem.Chem.Phys.2012,14,2317.doi:10.1039/c2cp23839j

(65) Wang,S.Y.;Yang,Q.Y.;Zhong,C.L.Sep.Purif.Technol.2008,60,30.doi:10.1016/j.seppur.2007.07.050

(66) Wu,D.;Wang,C.C.AIChE2011,58,2078.

Computational Study of Metal-Organic Frameworks for Removing H2S from Natural Gas

XU Hong1TONG Min-Man1WU Dong1XIAO Gang2,*YANG Qing-Yuan1LIU Da-Huan1,*ZHONG Chong-Li1
(1Laboratory of Advanced Nanostructure Materials,School of Chemical Engineering,Beijing University of Chemical Technology,Beijing 100029,P.R.China;2PetroChina Jilin Petrochemical Company,Jilin 132022,Jilin Province,P.R.China)

Natural gas is a highly efficient energy source subject to growing demand.Natural gas contains H2S,which can harm human health and cause equipment corrosion and environmental pollution.Effective H2S adsorbents are necessary to overcome these problems.Grand canonical Monte Carlo(GCMS)simulations were performed to study the selectivity and working capacity(pressure swing adsorption(PSA)and vacuum pressure swing adsorption(VSA)processes)of H2S in 33 kinds of stable metal-organic frameworks(MOFs),with the aim of separating H2S from H2S/CH4gas mixture.ZIF-80,Zn2-bpydtc,CAU-1-(OH)2,and CH3O-MOFa were suitable materials for the VSAprocess.CAU-1-(OH)2and CH3O-MOFa were suitable for the PSAprocess.The structures of materials exhibiting high selectivity and working capacity suggested that appropriate functionality and small pore sizes were important for high selectivity.MOFs with―Cl,―OH,and―OCH3functionality exhibited the strongest adsorption.Materials exhibiting high selectivity,strong interaction with H2S,and large numbers of adsorption sites may have high working capacities.High selectivity and high working capacity stable MOFs were screened and analyzed,to enhance the selectivity and working capacity toward H2S.This provides a theoreticalbasis for separating H2S from natural gas using MOFs.?Editorial office of Acta Physico-Chimica Sinica

Hydrogen sulfide;Metal-organic framework;Natural gas;Screen;Separation

O641

10.3866/PKU.WHXB201411132www.whxb.pku.edu.cn

Received:August 5,2014;Revised:November 13,2014;Published on Web:November 13,2014.

?Corresponding authors.LIU Da-Huan,Email:liudh@mail.buct.edu.cn;Tel:+86-10-64431705.XIAO Gang,Email:jh_xiaog@petrochina.com.cn; Tel:+86-10-64418935.

The project was supported by the National Key Basic Research Program of China(973)(2013CB733503),National Natural Science Foundation of China(21136001,21276008,21276009,21322603),and Doctoral Program of Higher Education of China(20110010130001).

國家重點基礎研究發展規劃項目(973)(2013CB733503),國家自然科學基金(21136001,21276008,21276009,21322603)和教育部博士點基金(20110010130001)資助