氣體在金屬-有機(jī)骨架材料中的吸附分離：經(jīng)典密度泛函理論的應(yīng)用

劉宇，趙雙良，胡軍，劉洪來(lái)

(1華東理工大學(xué)化工學(xué)院，國(guó)家化學(xué)工程重點(diǎn)實(shí)驗(yàn)室，上海 200237；2華東理工大學(xué)化學(xué)院，國(guó)家化學(xué)工程重點(diǎn)實(shí)驗(yàn)室，上海 200237)

?

氣體在金屬-有機(jī)骨架材料中的吸附分離：經(jīng)典密度泛函理論的應(yīng)用

劉宇1，趙雙良1，胡軍2，劉洪來(lái)2

(1華東理工大學(xué)化工學(xué)院，國(guó)家化學(xué)工程重點(diǎn)實(shí)驗(yàn)室，上海 200237；2華東理工大學(xué)化學(xué)院，國(guó)家化學(xué)工程重點(diǎn)實(shí)驗(yàn)室，上海 200237)

摘要：回顧了近年來(lái)經(jīng)典密度泛函理論在預(yù)測(cè)金屬-有機(jī)骨架材料吸附特性方面的研究進(jìn)展，重點(diǎn)介紹了經(jīng)過(guò)快速傅里葉變換加速后的經(jīng)典密度泛函理論在MOF吸附材料的大規(guī)模篩選方面的應(yīng)用。相較傳統(tǒng)的計(jì)算機(jī)分子模擬，加速后的經(jīng)典密度泛函理論的優(yōu)勢(shì)在于計(jì)算效率，對(duì)于簡(jiǎn)單的小分子氣體系統(tǒng)尤其具有優(yōu)勢(shì)，對(duì)MOF吸附材料進(jìn)行大規(guī)模篩選是可行的；但對(duì)于對(duì)復(fù)雜分子的處理尚沒有特別有效的方法，如何合理構(gòu)建復(fù)雜流體自由能泛函是它面臨的主要挑戰(zhàn)。

關(guān)鍵詞：金屬-有機(jī)骨架；經(jīng)典密度泛函理論；吸附（作用）；統(tǒng)計(jì)熱力學(xué)；氣體

2015-09-06收到初稿，2015-11-19收到修改稿。

聯(lián)系人：劉洪來(lái)。第一作者：劉宇(1984—)，男，博士。

Received date: 2015-09-06.

引 言

金屬-有機(jī)骨架（metal organic framework，MOF）是以金屬離子或金屬氧化物為節(jié)點(diǎn)，依靠有機(jī)配位體構(gòu)架出的具有三維多孔結(jié)構(gòu)的配位化合物。自從20世紀(jì)90年代被合成以來(lái)[1-3]，因其優(yōu)良的特性被認(rèn)為在吸附分離、催化以及傳感檢測(cè)等領(lǐng)域具有很好的潛在應(yīng)用前景。迄今為止，已經(jīng)合成出數(shù)千種金屬-有機(jī)骨架材料，并提出了數(shù)十億種假想的MOF材料結(jié)構(gòu)[4-5]。MOF材料之所以被認(rèn)為具有廣泛的應(yīng)用前景，關(guān)鍵在于其優(yōu)良的吸附性能，因此對(duì)MOF材料吸附特性的預(yù)測(cè)在MOF材料合成和應(yīng)用中有至關(guān)重要的意義。

對(duì)MOF材料吸附特性的預(yù)測(cè)，傳統(tǒng)的方法是通過(guò)計(jì)算機(jī)模擬。例如，Snurr等[4-5]采用Monte Carlo (MC)模擬預(yù)測(cè)了甲烷在137000種模型MOF材料中的吸附特性，并根據(jù)預(yù)測(cè)結(jié)果合成了著名的儲(chǔ)能MOF材料NU-100。Zhong等[6]采用分子動(dòng)力學(xué)(MD)模擬預(yù)測(cè)出MIL-47材料具有呼吸相變現(xiàn)象，隨后被實(shí)驗(yàn)所證實(shí)。Zhang等[7-8]采用MD模擬預(yù)測(cè)認(rèn)為吸附過(guò)程中ZIF-8材料的結(jié)構(gòu)可以發(fā)生轉(zhuǎn)變等。計(jì)算機(jī)模擬是目前預(yù)測(cè)MOF吸附最為準(zhǔn)確的方法，但由于其計(jì)算工作量巨大，效率不高，在實(shí)際應(yīng)用中受到一定的限制。例如，Snurr等對(duì)137000種MOF材料的大規(guī)模預(yù)測(cè)即建立在消耗大量計(jì)算資源的基礎(chǔ)上。開發(fā)更高效的計(jì)算方法是預(yù)測(cè)MOF材料吸附特性的當(dāng)務(wù)之急。經(jīng)典密度泛函理論(CDFT)似乎是一種較好的選擇。CDFT[9-15]建立于20世紀(jì)70年代，是一種基于統(tǒng)計(jì)力學(xué)的理論方法，被廣泛應(yīng)用于預(yù)測(cè)非均勻流體的微介觀性質(zhì)，如表界面效應(yīng)、浸濕、溶解等。CDFT在保證相同計(jì)算精度的前提下，往往能展現(xiàn)出比計(jì)算機(jī)模擬更高的計(jì)算效率，在一些特殊場(chǎng)合有望成為一種取代計(jì)算機(jī)模擬的理論方法。本文重點(diǎn)介紹應(yīng)用CDFT預(yù)測(cè)MOF材料的吸附特性以及目前的進(jìn)展。

1 氣體在MOF材料中吸附的經(jīng)典密度泛函理論

CDFT的基本思想是根據(jù)熱力學(xué)平衡條件——自由能(封閉系統(tǒng))或巨勢(shì)(開放系統(tǒng))極小原理，通過(guò)變分方法求解流體密度分布，進(jìn)而得到系統(tǒng)的其他物理量。在這個(gè)過(guò)程中，如何構(gòu)建自由能F或巨勢(shì)W關(guān)于流體密度分布的泛函形式是CDFT的核心問(wèn)題。具體的構(gòu)建方法與研究的體系有關(guān)，但一般說(shuō)來(lái)，均是將自由能分為理想氣體、外場(chǎng)作用、流體分子作用的剩余貢獻(xiàn)3項(xiàng)之和，即

式中，F(xiàn)id為被吸附氣體的理想自由能，即流體的動(dòng)能項(xiàng)貢獻(xiàn)；Fext為外場(chǎng)作用對(duì)自由能的貢獻(xiàn)，即流體與外場(chǎng)相互作用的貢獻(xiàn)；Fex為剩余自由能，即流體分子間相互作用的貢獻(xiàn)。

根據(jù)統(tǒng)計(jì)力學(xué)的相關(guān)理論，理想自由能可以精確地表示為流體密度分布r(r)的泛函形式

式中，L為de Broglie波長(zhǎng)；kB為Boltzmann常數(shù)；T為溫度。

外場(chǎng)作用對(duì)自由能的貢獻(xiàn)考慮的是外部條件，如外部電場(chǎng)、限制空間等的影響。如果外部條件不受流體的反饋影響，這部分貢獻(xiàn)可簡(jiǎn)單表示為

式中，Vext(r)為外勢(shì)場(chǎng)。以本文討論的MOF材料為例，如果將MOF結(jié)構(gòu)看成剛性的，則Vext(r)為定值，并可簡(jiǎn)單地表示為MOF中各原子與流體(吸附質(zhì))相互作用ui-fluid(r)的加和

而分子間作用勢(shì)，一般由短程的van der Waals相互作用和長(zhǎng)程的靜電相互作用構(gòu)成，最常用的勢(shì)能函數(shù)為L(zhǎng)J12-6 + Coulombic勢(shì)

式中，eij和sij為L(zhǎng)J參數(shù)，可由標(biāo)準(zhǔn)力場(chǎng)(Dreiding、UFF等)以及相應(yīng)的結(jié)合規(guī)則得到；z為電荷數(shù)；e為元電荷電量；e 為介電常數(shù)。

剩余自由能Fex反映了分子間的復(fù)雜相互作用，一般沒有精確的表達(dá)式，只能采用近似方法。常用的近似方法有基本度量理論(FMT)[9-11]、加權(quán)密度近似(WDA)[12-15]、泛函展開近似(FEA)[16-17]、平均場(chǎng)近似(MFA)[15,18]等。一般是將剩余自由能分為硬球排斥和吸引兩部分貢獻(xiàn)之和

Yu采用FMT、MFA以及WDA分別對(duì)式(7)中的硬球排斥項(xiàng)、吸引項(xiàng)和相關(guān)項(xiàng)進(jìn)行近似，并將其應(yīng)用于研究狹縫中LJ流體的相行為，得到了與計(jì)算機(jī)模擬高度吻合的結(jié)果。最近，該方法又被Liu等[14]進(jìn)一步推廣到研究MOF材料中的吸附現(xiàn)象，并得到了很好的結(jié)果。

圖1 H2在MOF-5中的吸附等溫線及吸附熱曲線[21]Fig.1 Adsorption isotherm of H2in MOF-5[21]solid line: double WDA; dashed line: WDA + MFA;：MC simulation results；：experimental results

2 氣體在MOF材料中吸附的CDFT

CDFT自從在20世紀(jì)70年代被提出以來(lái)，長(zhǎng)期被應(yīng)用于具有高對(duì)稱幾何結(jié)構(gòu)的體系，如狹縫型孔道、圓柱型孔道等。而被應(yīng)用于具有三維復(fù)雜空間結(jié)構(gòu)的體系，如本文所討論的MOF，則是最近五六年才發(fā)展起來(lái)的。

2009年，Siderius等[20]將CDFT應(yīng)用于MOF材料的吸附研究。他們對(duì)硬球貢獻(xiàn)項(xiàng)采用WDA+ MFA方法研究了高溫下H2在MOF-5中的吸附現(xiàn)象，并得到了與實(shí)驗(yàn)和模擬相吻合的結(jié)果。但Liu等[21]的計(jì)算結(jié)果表明該方法對(duì)于低溫系統(tǒng)的計(jì)算效果并不令人滿意。為此，Liu等[21]采用雙WDA近似改進(jìn)了該方法，使得CDFT能夠在全溫度區(qū)間內(nèi)準(zhǔn)確預(yù)測(cè)MOF材料的吸附特性(圖1)。求解CDFT可以得到的是流體密度分布，因而可以很方便地得到MOF材料的吸附位點(diǎn)等微觀結(jié)構(gòu)信息。Siderius和Liu等的工作均預(yù)測(cè)出H2在MOF-5中的最佳吸附位點(diǎn)位于金屬離子附近(圖2)。Liu等[14]進(jìn)一步將CDFT推廣到混合氣體在MOF材料中吸附的研究，考察了CO2/CH4和CO2/N2兩類混合氣體在MOF材料中的吸附及分離特征，并得到了與MC模擬一致的結(jié)果。

3 快速傅里葉變換與吸附材料的大規(guī)模篩選

圖2 H2在MOF-5中的密度分布（等密度面）[20]Fig. 2 Density profiles of H2in MOF-5 (isosurface)[20]gray: 1.33 mol·L-1; green: 19.9 mol·L-1

3.1 快速傅里葉變換在CDFT中的應(yīng)用

在采用CDFT預(yù)測(cè)MOF材料對(duì)氣體及其混合物的吸附現(xiàn)象時(shí)，需要對(duì)密度分布方程進(jìn)行三維數(shù)值求解，由于MOF材料具有復(fù)雜的三維孔道結(jié)構(gòu)，需要根據(jù)不同MOF材料的具體結(jié)構(gòu)將孔道空間劃分成不同的計(jì)算格點(diǎn)，并采用合適的數(shù)值求解方法得到這些格點(diǎn)上的氣體密度，高精度的計(jì)算需要高密度的格點(diǎn)劃分，而高密度的格點(diǎn)劃分則導(dǎo)致計(jì)算量的猛增。實(shí)際計(jì)算中發(fā)現(xiàn)，CDFT需要的平均時(shí)間較MC模擬高達(dá)兩個(gè)數(shù)量級(jí)以上。為解決這個(gè)問(wèn)題，Wu等[15]引入了快速傅里葉變換方法。

采用CDFT處理具有復(fù)雜三維結(jié)構(gòu)的密度分布時(shí)，最為耗時(shí)的部分在于剩余自由能Fex的計(jì)算。以WDA近似為例，剩余自由能可表示為

如果直接用數(shù)值積分的方法計(jì)算式(9)，其計(jì)算復(fù)雜性為N2D（N為每個(gè)維度上的區(qū)間數(shù)，D為維數(shù)）。對(duì)于D=3的MOF體系來(lái)說(shuō)，其復(fù)雜性相當(dāng)高。而通過(guò)傅里葉變換，式(9)可在傅里葉空間中簡(jiǎn)化為

其中r%(k )表示函數(shù)r(r )的傅里葉變換。由于快速傅里葉變換的計(jì)算復(fù)雜度僅為NDlog2ND，采用式(10)取代式(9)可以大幅度提升計(jì)算效率。這是一個(gè)通用的方法，CDFT的其他近似如FMT、泛函展開等均可被快速傅里葉變換加速。而實(shí)踐表明，采用快速傅里葉變換加速后的CDFT計(jì)算所需時(shí)間比MC模擬少兩個(gè)數(shù)量級(jí)以上。

快速傅里葉變換的引入，使得CDFT得以超越MC成為預(yù)測(cè)MOF材料吸附性能更為高效的工具，也極大拓廣了CDFT的應(yīng)用范圍。其中一個(gè)最重要的應(yīng)用即是對(duì)MOF吸附材料的大規(guī)模篩選。

3.2 大規(guī)模篩選

MOF材料種類千差萬(wàn)別，通過(guò)不同點(diǎn)群、不同金屬離子節(jié)點(diǎn)和有機(jī)連接體的排列組合可以生成的結(jié)構(gòu)數(shù)以億計(jì)，現(xiàn)今已提出數(shù)十億種假想MOF。如何才能從數(shù)量如此龐大的MOF材料中篩選出具有所需要吸附性能的結(jié)構(gòu)呢？采用實(shí)驗(yàn)手段顯然是不現(xiàn)實(shí)的，必須采用“先篩選，后合成”的理念：即首先對(duì)這些可選的MOF材料進(jìn)行一輪理論篩選，找出少數(shù)幾種最具潛力的結(jié)構(gòu)，然后僅對(duì)這少數(shù)幾種結(jié)構(gòu)進(jìn)行合成和測(cè)試。在這個(gè)過(guò)程中間，篩選方法至關(guān)重要，必須同時(shí)具有準(zhǔn)確性和高效率的特點(diǎn)。傳統(tǒng)的做法是通過(guò)計(jì)算機(jī)模擬，如Snurr等[4-5,22-28]在MOF材料儲(chǔ)氫和儲(chǔ)甲烷等方面的工作以及Zhong等[29-31]在MOF材料吸附分離二氧化碳等方面的工作。但計(jì)算機(jī)模擬耗時(shí)量大，效率較低，對(duì)于大規(guī)模篩選來(lái)說(shuō)難以承受。傅里葉變換加速后的CDFT則有著明顯的優(yōu)勢(shì)，由于在計(jì)算速度上有著兩個(gè)數(shù)量級(jí)的優(yōu)勢(shì)，CDFT能更好地勝任大規(guī)模篩選。

將CDFT應(yīng)用于MOF吸附材料的大規(guī)模篩選，最初是由Wu等[15]完成的。他們考察了氫氣和甲烷分別在1200種MOF材料中的吸附特性，得到了與計(jì)算機(jī)模擬高度一致的結(jié)果(圖3)。

Wu等采用MFMT對(duì)硬球項(xiàng)進(jìn)行近似，而分別采用了4種不同的近似方法計(jì)算吸引項(xiàng)貢獻(xiàn)，包括兩種不同版本的WDA近似(Yu版[19]和Liu版[15,21])、MFA以及泛函展開近似。結(jié)果表明，4種方法的計(jì)算效果在高溫下難分伯仲，而在低溫下，Liu版WDA最為準(zhǔn)確。計(jì)算結(jié)果還表明，在WDA中，狀態(tài)方程的優(yōu)劣對(duì)計(jì)算結(jié)果有至關(guān)重要的影響。如圖4所示，采用更為精確的MBWR (modified benedict-Webb-Rubin) 狀態(tài)方程的預(yù)測(cè)結(jié)果明顯優(yōu)于FMSA(first-order mean spherical approximation)。

圖3 CDFT和MC預(yù)測(cè)氫氣在1200種MOF材料中的吸附量的比較[15]Fig. 3 Comparison between CDFT and MC prediction of H2adsorption in 1200 types MOF materials[15]

圖4 選用不同狀態(tài)方程對(duì)預(yù)測(cè)結(jié)果的影響[15]Fig. 4 Effect of different equation of state on predicting results[15]

Wu等以及Snurr等和Zhong等的計(jì)算機(jī)模擬都主要集中在少數(shù)幾個(gè)特定溫度和壓力條件下預(yù)測(cè)MOF材料的吸附特征。最近，Liu等[32]將預(yù)測(cè)范圍進(jìn)一步拓廣，預(yù)測(cè)了712種MOF材料分別在441個(gè)不同溫度和壓力條件下的吸附特征，篩選出不同外部條件下的最佳儲(chǔ)氫MOF材料。根據(jù)Liu等的計(jì)算結(jié)果(圖5)，在77 K、10 MPa條件下，MOF 對(duì)H2的吸附量可達(dá)37 mol·L-1，甚至超出液氫的的密度(35 mol·L-1)；但在高溫條件下，MOF材料的儲(chǔ)氫性能卻并未有質(zhì)的提升，依然在5 mol·L-1以下。

圖5 不同溫度和壓力條件下712種MOF材料中的最大H2吸附量[32]Fig. 5 Highest adsorption amount for 712 types of MOF material on different temperatures and pressures[32]

Liu等發(fā)現(xiàn)，H2在MOF材料中的吸附等溫線可以用Langmuir吸附方程進(jìn)行擬合，即

發(fā)現(xiàn)不同MOF材料的兩個(gè)Langmuir參數(shù)p0和N0呈現(xiàn)出明顯的相關(guān)性(圖6)，尤其在高溫下，這種相關(guān)性十分明顯而且是線性的。

最近，Liu等[33]又將CDFT預(yù)測(cè)MOF吸附的方法進(jìn)一步應(yīng)用于吸附脫硫領(lǐng)域。雖然該工作中涉及的噻吩類硫化物存在于液相，但其行為在低密度下與氣體分子基本一致。該工作表明CDFT不僅能預(yù)測(cè)小分子氣體在MOF中的吸附特性，而且能預(yù)測(cè)諸如噻吩類硫化物等較為復(fù)雜的分子的吸附特征。與預(yù)測(cè)小分子氣體類似，該工作也通過(guò)大規(guī)模CDFT計(jì)算篩選出了針對(duì)不同硫濃度下的最佳脫硫吸附劑，可以為實(shí)驗(yàn)合成新型吸附脫硫材料提供指導(dǎo)。

在計(jì)算速度方面，通過(guò)Wu等和Liu等的工作，可以明顯看出CDFT的優(yōu)越性。例如，Liu等對(duì)MOF儲(chǔ)氫方面的工作包含了712 × 441 = 313992 次CDFT計(jì)算而僅耗1500 CPU小時(shí)，這是目前所有分子模擬方法所望塵莫及的。采用CDFT對(duì)MOF材料進(jìn)行大規(guī)模篩選是一種極具潛力的方法。

圖6 712種MOF材料的Langmuir參數(shù)的相關(guān)性[32]Fig. 6 Correlation between Langmuir parameters[32]

4 對(duì)擴(kuò)散系數(shù)的預(yù)測(cè)

在MOF材料應(yīng)用于氣體吸附分離時(shí)，除了吸附量和吸附選擇性，吸脫附速率也是評(píng)價(jià)材料性能的重要指標(biāo)。例如，以MOF材料儲(chǔ)氫取代汽油作為發(fā)動(dòng)機(jī)燃料時(shí)，雖然總的儲(chǔ)能量取決于MOF材料的吸附量，但其輸出功率則主要依賴于MOF材料的吸脫附速率，即MOF材料的吸附動(dòng)力學(xué)屬性。

擴(kuò)散系數(shù)是考察MOF材料吸附動(dòng)力學(xué)屬性的一個(gè)極為關(guān)鍵的指標(biāo)。一般來(lái)說(shuō)，擴(kuò)散系數(shù)越大，則吸脫附速率越高。對(duì)氣體在多孔材料中擴(kuò)散系數(shù)的預(yù)測(cè)，目前主要是通過(guò)分子動(dòng)力學(xué)模擬完成[34-36]。然而，與預(yù)測(cè)吸附量類似，速度依然是計(jì)算機(jī)模擬的短板。尤其是對(duì)于擴(kuò)散系數(shù)的預(yù)測(cè)，模擬所消耗的計(jì)算資源還將遠(yuǎn)大于預(yù)測(cè)吸附量。如何快速預(yù)測(cè)擴(kuò)散系數(shù)，依然是目前亟待解決的問(wèn)題。

雖然CDFT不能直接預(yù)測(cè)擴(kuò)散系數(shù)這種動(dòng)力學(xué)物理量，但通過(guò)剩余熵標(biāo)度理論[37-38]，可以將CDFT和擴(kuò)散系數(shù)聯(lián)系起來(lái)。剩余熵標(biāo)度理論給出了擴(kuò)散系數(shù)和剩余熵的普適關(guān)系式

式中，D*為量綱1化后的擴(kuò)散系數(shù)；sex*為量綱1化后的剩余熵；A和B為普適常數(shù)。該關(guān)系式最早由Rosenfeld[32]于1977年提出，并通過(guò)擬合得到A = 0.585，B = 0.788。Dzugutov[38-39]指出，式(12)的指數(shù)關(guān)系源于等概率原理和Boltzmann關(guān)系

如果各態(tài)歷經(jīng)假設(shè)能夠成立的話，普適常數(shù)B應(yīng)為1，否則B< 1。

圖7 298 K下H2在4種多孔材料中的擴(kuò)散系數(shù)[47]Fig. 7 Diffusion coefficient of H2in four types of porous material[47]

剩余熵標(biāo)度理論建立了熱力學(xué)性質(zhì)和動(dòng)力學(xué)參數(shù)之間的橋梁，使得人們可以通過(guò)計(jì)算熱力學(xué)的剩余熵來(lái)預(yù)測(cè)動(dòng)力學(xué)的擴(kuò)散系數(shù)。該方法最初被應(yīng)用于均勻流體[37,40-42]：首先通過(guò)狀態(tài)方程計(jì)算剩余熵，然后根據(jù)剩余熵標(biāo)度理論預(yù)測(cè)擴(kuò)散系數(shù)。后來(lái)，剩余熵標(biāo)度理論又被進(jìn)一步推廣到非均勻流體[41,43-46]，即通過(guò)計(jì)算機(jī)模擬計(jì)算非均勻流體的剩余熵，然后依靠剩余熵標(biāo)度理論來(lái)預(yù)測(cè)非均勻流體的擴(kuò)散系數(shù)。由于在計(jì)算機(jī)模擬過(guò)程中，熵的計(jì)算亦需要較大的計(jì)算開銷，因此該方法并未得到廣泛的應(yīng)用。最近，Liu等[47]采用CDFT計(jì)算剩余熵，并結(jié)合Knudsen理論、平均自由體積理論成功預(yù)測(cè)了MOF材料中的擴(kuò)散系數(shù)(圖7)。

不同于計(jì)算機(jī)模擬，CDFT能很方便地計(jì)算非均勻流體的熵和剩余熵等熱力學(xué)信息

因此，采用CDFT替代模擬計(jì)算剩余熵，然后通過(guò)剩余熵標(biāo)度理論計(jì)算擴(kuò)散系數(shù)能避免大量耗時(shí)的模擬抽樣，大幅度提高計(jì)算效率。從目前來(lái)看，該方法亦是將剩余熵標(biāo)度理論推廣到非均勻流體的最高效的辦法。

5 結(jié)論與展望

MOF作為一種高性能納米孔道材料有著良好的應(yīng)用前景。對(duì)其吸附性能的預(yù)測(cè)是材料篩選的關(guān)鍵。傳統(tǒng)的計(jì)算機(jī)模擬由于計(jì)算速度的限制，計(jì)算效率不高，最近發(fā)展起來(lái)的快速傅里葉變換CDFT則能很好地解決這一問(wèn)題，它可以在保持精算精度不變的前提下比分子模擬方法大幅提高計(jì)算效率。但CDFT也存在一些不足，其中最主要的問(wèn)題是如何處理具有復(fù)雜空間結(jié)構(gòu)分子的吸附問(wèn)題，尤其是當(dāng)分子偶極矩很大的時(shí)候。以水分子為例，如果簡(jiǎn)單地將水分子考慮成具有介電常數(shù)e 的背景介質(zhì)，則無(wú)法反映實(shí)際MOF材料吸附過(guò)程中因水分子占據(jù)最佳吸附位點(diǎn)而阻礙其他分子(如CO2等)吸附的現(xiàn)象；如果將水分子粗?；癁楹?jiǎn)單的球形分子，則無(wú)法正確反映出水分子偶極矩和靜電屏蔽效應(yīng)對(duì)吸附的影響；如果采用全原子模型(如SPC、TIP3P等)，雖然最為接近實(shí)際情況，但如何構(gòu)建相應(yīng)的剩余自由能泛函Fex又是目前CDFT中尚未解決的問(wèn)題。因此，對(duì)CDFT理論本身的改進(jìn)是將CDFT推廣到更為復(fù)雜體系的關(guān)鍵。在當(dāng)前CDFT的理論基礎(chǔ)上，結(jié)合分子模擬來(lái)預(yù)測(cè)流體在MOF材料中的行為可能是一種較有潛力的方法。例如，通過(guò)分子模擬計(jì)算復(fù)雜分子的構(gòu)象而通過(guò)CDFT計(jì)算簡(jiǎn)單流體的密度分布；通過(guò)分子模擬計(jì)算流體的密度分布而通過(guò)CDFT計(jì)算體系的自由能等宏觀熱力學(xué)量等。總而言之，CDFT理論本身的完善以及CDFT與分子模擬的結(jié)合是CDFT預(yù)測(cè)氣體在MOF中吸附現(xiàn)象的兩條主要途徑。

References

[1] YOUNG K M, MILLER T M, WRIGHTON M S. Comparison of the thermal and photochemical reactions of (.eta.1-cyclopentadienyl) rhenium pentacarbonyl and (.eta.1-9-fluorenyl)rhenium pentacarbonyl: nonthermal chemical reactions from the lowest excited state [J]. Journal of the American Chemical Society, 1990, 112(4): 1529-1537.

[2] FUJITA M, KWON Y J, WASHIZU S, et al. Preparation, clathration ability, and catalysis of a two-dimensional square network material composed of cadmium (Ⅱ) and 4,4'-bipyridine [J]. Journal of the American Chemical Society, 1994, 116(3): 1151-1152.

[3] VENKATARAMAN D, GARDNER G B, LEE S, et al. Zeolite-like behavior of a coordination network [J]. Journal of the American Chemical Society, 1995, 117(46): 11600-11601.

[4] FARHA O K, YAZAYD?N A ?, ERYAZICI I, et al. De novo synthesis of a metal-organic framework material featuring ultrahigh surface area and gas storage capacities [J]. Nature Chemistry, 2010, 2(11): 944-948.

[5] FERNANDEZ M, WOO T K, WILMER C E, et al. Large-scale quantitative structure-property relationship (QSPR) analysis of methane storage in metal-organic frameworks [J]. Journal of Physical Chemistry C, 2013, 117(15): 7681-7689.

[6] MA Q T, YANG Q Y, GHOUFI A, et al. Guest-modulation of the mechanical properties of flexible porous metal-organic frameworks [J]. Journal of Materials Chemistry A, 2014, 2(25): 9691-9698.

[7] ZHANG L L, HU Z Q, JIANG J W. Sorption-induced structural transition of zeolitic imidazolate framework-8: a hybrid molecular simulation study [J]. Journal of the American Chemical Society, 2013, 135(9): 3722-3728.

[8] ZHANG L L, WU G, JIANG J W. Adsorption and diffusion of CO2and CH4in zeolitic imidazolate framework-8: effect of structural flexibility [J]. Journal of Physical Chemistry C, 2014, 118(17): 8788-8794.

[9] ROSENFELD Y. Free-energy model for the inhomogeneous hard-sphere fluid mixture and density-functional theory of freezing [J]. Phys. Rev. Lett., 1989, 63(9): 980-983.

[10] ROSENFELD Y. Phase separation of asymmetric binary hard-sphere fluids: self-consistent density functional theory [J]. Physical Review Letters, 1994, 72(24): 3831-3834.

[11] YU Y X, WU J Z. Structures of hard-sphere fluids from a modified fundamental-measure theory [J]. Journal of Chemical Physics, 2002, 117(22): 10156-10164.

[12] PARR R G. Density-functional Theory of Atoms and Molecules [M]. New York, Oxford: Oxford University Press, 1989.

[13] KIERLIK E, ROSINBERG M L. Density-functional theory for inhomogeneous fluids — adsorption of binary mixtures [J]. Physical Review A, 1991, 44(8): 5025-5037.

[14] LIU Y, LIU H L, HU Y, et al. Density functional theory for adsorption of gas mixtures in metal-organic frameworks [J]. Journal of Physical Chemistry B, 2010, 114(8): 2820-2827.

[15] FU J, LIU Y, TIAN Y, et al. Density functional methods for fast screening of metal-organic frameworks for hydrogen storage [J]. Journal of Physical Chemistry C, 2015, 119: 5374-5385.

[16] YU Y X, YOU F Q, TANG Y P, et al. Structure and adsorption of a hard-core multi-Yukawa fluid confined in a slitlike pore: grand canonical Monte Carlo simulation and density functional study [J]. Journal of Physical Chemistry B, 2006, 110(1): 334-341.

[17] TANG Y P, LU B C Y. Analytical description of the Lennard-Jones fluid and its application [J]. AIChE Journal, 1997, 43(9): 2215-2226.

[18] NEIMARK A V, RAVIKOVITCH P I. Capillary condensation in MMS and pore structure characterization [J]. Microporous and Mesoporous Materials, 2001, 44: 697-707.

[19] YU Y X. A novel weighted density functional theory for adsorption, fluid-solid interfacial tension, and disjoining properties of simple liquid films on planar solid surfaces [J]. Journal of Chemical Physics, 2009, 131(2): 024704.

[20] SIDERIUS D W, GELB L D. Predicting gas adsorption in complex microporous and mesoporous materials using a new density functional theory of finely discretized lattice fluids [J]. Langmuir, 2009, 25(3): 1296-1299.

[21] LIU Y, LIU H L, HU Y, et al. Development of a density functional theory in three-dimensional nanoconfined space: H2storage in metal organic frameworks [J]. Journal of Physical Chemistry B, 2009, 113(36): 12326-12331.

[22] BORAH B, ZHANG H, SNURR R Q. Diffusion of methane and other alkanes in metal-organic frameworks for natural gas storage [J]. Chemical Engineering Science, 2015, 124: 135-143.

[23] KIM K C, MOGHADAM P Z, FAIREN-JIMENEZ D, et al. Computational screening of metal catecholates for ammonia capture in metal-organic frameworks [J]. Industrial & Engineering Chemistry Research, 2015, 54(13): 3257-3267.

[24] SIMON C M, KIM J, GOMEZ-GUALDRON D A, et al. The materials genome in action: identifying the performance limits for methane storage [J]. Energy & Environmental Science, 2015, 8(4): 1190-1199.

[25] ZHANG H D, DERIA P, FARHA O K, et al. A thermodynamic tank model for studying the effect of higher hydrocarbons on natural gas storage in metal-organic frameworks [J]. Energy & Environmental Science, 2015, 8(5): 1501-1510.

[26] CHUNG Y G, CAMP J, HARANCZYK M, et al. Computation-ready, experimental metal-organic frameworks: a tool to enable high-throughput screening of nanoporous crystals [J]. Chemistry of Materials, 2014, 26(21): 6185-6192.

[27] COLON Y J, SNURR R Q. High-throughput computational screening of metal-organic frameworks [J]. Chemical Society Reviews, 2014, 43(16): 5735-5749.

[28] GOMEZ-GUALDRON D A, GUTOV O V, KRUNGLEVICIUTE V, et al. Computational design of metal-organic frameworks based on stable zirconium building units for storage and delivery of methane [J]. Chemistry of Materials, 2014, 26(19): 5632-5639.

[29] LIU D H, ZHONG C L. Understanding gas separation in metal-organic frameworks using computer modeling [J]. Journal of Materials Chemistry, 2010, 20(46): 10308-10318.

[30] WU D, WANG C C, LIU B, et al. Large-scale computational screening of metal-organic frameworks for CH4/H2separation [J]. AIChE Journal, 2012, 58(7): 2078-2084.

[31] LI Z J, XIAO G, YANG Q Y, et al. Computational exploration of metal-organic frameworks for CO2/CH4separation via temperature swing adsorption [J]. Chemical Engineering Science, 2014, 120: 59-66.

[32] LIU Y, ZHAO S L, LIU H L, et al. High-throughput and comprehensive prediction of H2adsorption in metal-organic frameworks under various conditions [J]. AIChE Journal, 2015, 61(9): 2951-2957. DOI: 10.1002/aic.14842.

[33] LIU Y, GUO F Y, HU J, et al. Screening of desulfurization adsorbent in metal-organic frameworks: a classical density functional approach [J]. Chemical Engineering Science, 2015, 137(1): 170-177.

[34] SKOULIDAS A I, SHOLL D S. Transport diffusivities of CH4, CF4, He, Ne, Ar, Xe, and SF6in silicalite from atomistic simulations [J]. Journal of Physical Chemistry B, 2002, 106(19): 5058-5067.

[35] BABARAO R, JIANG J W. Diffusion and separation of CO2and CH4in silicalite, C-168 Schwarzite,and IRMOF-1: a comparative study from molecular dynamics simulation [J]. Langmuir, 2008, 24(10): 5474-5484.

[36] SKOULIDAS A I, SHOLL D S. Self-diffusion and transport diffusion of light gases in metal-organic framework materials assessed using molecular dynamics simulations [J]. Journal of Physical Chemistry B, 2005, 109(33): 15760-15768.

[37] ROSENFELD Y. Relation between transport-coefficients and internal entropy of simple systems [J]. Physical Review A, 1977, 15(6): 2545-2549.

[38] DZUGUTOV M. A universal scaling law for atomic diffusion in condensed matter [J]. Nature, 1996, 381(6578): 137-139.

[39] DZUGUTOV M. Anomalous slowing down in the metastable liquid of hard spheres [J]. Physical Review E, 2002, 65(3): 032501.

[40] ROSENFELD Y. A quasi-universal scaling law for atomic transport in simple fluids [J]. Journal of Physics-Condensed Matter, 1999, 11(28): 5415-5427.

[41] MITTAL J, ERRINGTON J R, TRUSKETT T M. Relationships between self-diffusivity, packing fraction, and excess entropy in simple bulk and confined fluids [J]. Journal of Physical Chemistry B, 2007, 111(34): 10054-10063.

[42] VAZ R V, MAGALHAES A L, FERNANDES D L A, et al. Universal correlation of self-diffusion coefficients of model and real fluids based on residual entropy scaling law [J]. Chemical Engineering Science, 2012, 79: 153-162.

[43] CARMER J, GOEL G, POND M J, et al. Enhancing tracer diffusivity by tuning interparticle interactions and coordination shell structure [J]. Soft Matter, 2012, 8(15): 4083-4089.

[44] HE P, LI H Q, HOU X J. Excess-entropy scaling of dynamics for methane in various nanoporous materials [J]. Chemical Physics Letters, 2014, 593: 83-88.

[45] HE P, LIU H, ZHU J Q, et al. Tests of excess entropy scaling laws for diffusion of methane in silica nanopores [J]. Chemical Physics Letters, 2012, 535: 84-90.

[46] CHOPRA R, TRUSKETT T M, ERRINGTON J R. Excess-entropy scaling of dynamics for a confined fluid of dumbbell-shaped particles [J]. Physical Review E, 2010, 82(4): 041201.

[47] LIU Y, FU J, WU J Z. Excess-entropy scaling for gas diffusivity in nanoporous materials [J]. Langmuir, 2013, 29(42): 12997-13002.

Foundation item: supported by the National Natural Science Foundation of China (91334203, 21506051, 21176066), the Shanghai Pujiang Program (15PJ1401400), the Fundamental Research Funds for the Central Universities of China (222201414008) and the Open Project of State Key Laboratory of Chemical Engineering (SKL-ChE-15C05).

Gas adsorption and separation in metal-organic framework: application of classical density functional theory

LIU Yu1, ZHAO Shuangliang1, HU Jun2, LIU Honglai2
(1State Key Laboratory of Chemical Engineering and Department of Chemical Engineering, East China University of Science and Technology, Shanghai 200237, China;2State Key Laboratory of Chemical Engineering and Department of Chemistry, East China University of Science and Technology, Shanghai 200237, China)

Abstract:This article reviews the recent progress on predicting the adsorption properties of metal-organic framework by using classical density functional theory and focused on the application of the classical density functional theory to the high-throughput screening, which is accelerated by Fast Fourier Transform. Comparing to the conventional molecular simulations, the advantage of the accelerated classical density functional theory is the calculation speed, especially for simple small molecule systems, which makes the high-throughput screening on MOF materials feasible. However, it appears that there is a lack of efficient method to deal with the complicated molecules. How to construct a reasonable free energy functional of complicated fluid is the main challenge to state of art classical density functional theory.

Key words:metal-organic framework; classical density functional theory; adsorption; statistical thermodynamics; gas

Corresponding author:Prof. LIU Honglai, hlliu@ecust.edu.cn

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目(91334203, 21506051, 21176066)；上海市浦江人才計(jì)劃(15PJ1401400)；中央高?；究蒲袠I(yè)務(wù)費(fèi)(222201414008)；化學(xué)工程聯(lián)合國(guó)家重點(diǎn)實(shí)驗(yàn)室開放課題（SKL-ChE-15C05）。

中圖分類號(hào)：O 642.4

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：0438—1157（2016）01—0089—08

DOI:10.11949/j.issn.0438-1157.20151403