孫夢(mèng)遙 陳迪明

關(guān)鍵詞： 計(jì)算模擬軟件 晶態(tài)多孔材料 材料工作室 構(gòu)效-關(guān)系研究 氣體吸附分離

儲(chǔ)存和分離是化工生產(chǎn)過(guò)程中的重要操作過(guò)程。隨著現(xiàn)代工業(yè)的不斷發(fā)展，由CO2氣體排放所導(dǎo)致的溫室氣體效應(yīng)越來(lái)越受到人們的關(guān)注，因此CO2氣體的捕獲與封裝成為了科學(xué)界研究的熱點(diǎn)領(lǐng)域之一[1]。另一方面，乙炔、乙烯氣體被廣泛用作工業(yè)和日常生活中的能源物質(zhì)，然而其生產(chǎn)過(guò)程中會(huì)不可避免地混入其他氣體雜質(zhì)。開(kāi)發(fā)經(jīng)濟(jì)節(jié)能的多孔材料用于儲(chǔ)存和分離氣體，不僅從科學(xué)和工業(yè)角度上具有重要意義，并且還是能源化學(xué)與材料研究中的一個(gè)重要方面[2]。金屬-有機(jī)框架（Metal-Organic Frameworks， MOFs）作為一種新型的晶態(tài)多孔材料，具有可調(diào)節(jié)的孔徑與孔環(huán)境、大的比表面積及易于功能化等方面特點(diǎn)，使其在溫室氣體的捕獲、能源氣體的儲(chǔ)存與分離方面有著重要的應(yīng)用前景。MOFs 的上述結(jié)構(gòu)特點(diǎn)，不僅能夠使其作為吸附劑來(lái)實(shí)現(xiàn)對(duì)溫室氣體或能源氣體的高密度的儲(chǔ)存，并且能針對(duì)不同氣體分子產(chǎn)生具有差異性的主客體相互作用和/或篩分效應(yīng)來(lái)實(shí)現(xiàn)對(duì)氣體經(jīng)濟(jì)節(jié)能的分離目的[3]。

計(jì)算模擬軟件主要是依靠在計(jì)算機(jī)輸入結(jié)構(gòu)模型的基礎(chǔ)上，利用分子動(dòng)力學(xué)原理對(duì)材料可能表現(xiàn)出來(lái)的性質(zhì)進(jìn)行預(yù)測(cè)，與實(shí)驗(yàn)結(jié)果相比，使材料表現(xiàn)出來(lái)的性質(zhì)從原子層面上進(jìn)行解釋。相較于傳統(tǒng)的正交實(shí)驗(yàn)型研究，分子模擬具有成本低、實(shí)驗(yàn)風(fēng)險(xiǎn)小、實(shí)驗(yàn)結(jié)果可視化程度高等優(yōu)點(diǎn)。

隨著計(jì)算機(jī)硬件處理速度的提高，計(jì)算模擬軟件在多孔晶態(tài)材料氣體的吸附與分離機(jī)制方面展現(xiàn)出了重要的作用，具體表現(xiàn)在氣體分子初始載荷吸附位點(diǎn)的確定、氣體分子初始載荷吸附焓的計(jì)算、氣體分子在特定壓力下的吸附量及分布密度圖的計(jì)算以及特定溫度下氣體吸附等溫線(xiàn)的計(jì)算[4]。通過(guò)計(jì)算模擬軟件，可以有效預(yù)測(cè)材料的吸附分離性能，分析構(gòu)效關(guān)系以便更好地從原子層面理解材料與吸附分子之間的相互作用特性。

目前，比較常用的計(jì)算模擬軟件主要有美國(guó)西北大學(xué)Snurr 教授編寫(xiě)的Music 軟件（最新版的已更名為RASPA2）和BIOVIA 公司開(kāi)發(fā)的Materials studio（MS）軟件。前者由于只能于Unix 和Linux 操作系統(tǒng)上運(yùn)行，對(duì)于沒(méi)有相關(guān)操作系統(tǒng)基礎(chǔ)的研究人員有一定難度。而后者M(jìn)S 軟件為Windows 平臺(tái)的軟件，所有操作均可以由鼠標(biāo)點(diǎn)擊完成，并不需要輸入代碼來(lái)執(zhí)行命令，其友好的界面具有容易上手等優(yōu)點(diǎn)，并使其在晶態(tài)多孔材料氣體吸附與分離的研究工作方面得到了廣泛應(yīng)用。此外，MS 軟件還集成了分子動(dòng)力學(xué)（Forcite）、蒙卡特羅（Sorption）以及量子力學(xué)（DMol3）等多種功能模塊，使得幾乎所有計(jì)算工作可以一站式完成。在該研究中，筆者從模型構(gòu)建、吸附位點(diǎn)及吸附能計(jì)算、吸附等溫線(xiàn)計(jì)算對(duì)MS 軟件在晶態(tài)多孔材料吸與分離方面的應(yīng)用進(jìn)行介紹。

1 模型的構(gòu)建

1. 1 確定初始結(jié)構(gòu)模型

鑒于MOFs 材料結(jié)構(gòu)具有長(zhǎng)程有序的特點(diǎn)，使其絕大部分的結(jié)構(gòu)可以由X-射線(xiàn)單晶衍射分析獲得。值得注意的是，由于MOFs 材料在進(jìn)行氣體吸附與分離之前一般需要脫除客體分子，所以需要比對(duì)活化之后的樣品粉末衍射峰與合成樣品粉末衍射峰，二者區(qū)分不大的情況下可以直接用測(cè)得的樣品結(jié)構(gòu)數(shù)據(jù)作為MS 軟件計(jì)算的輸入文件（結(jié)晶的溶劑也會(huì)對(duì)粉末衍射峰產(chǎn)生貢獻(xiàn)，因此在脫除溶劑分子后MOFs 的粉末衍射峰會(huì)發(fā)生一些微小的變化）。如果MOFs 活化后粉末發(fā)生了較大的變化，例如：出現(xiàn)了整體衍射峰的偏移或部分舊的衍射峰的消失或者新的衍射峰的出現(xiàn)，則不能以合成MOFs 的結(jié)構(gòu)作為輸入文件。對(duì)得不到單晶結(jié)構(gòu)的MOFs 材料，可以通過(guò)粉末精修的方法得到MOFs 的結(jié)構(gòu)信息。ZHOU D D 等人利用Pawley 精修方法得到了一例銅基MOF 的結(jié)構(gòu)信息，并用作計(jì)算模擬的初始模型[5]。

1. 2 確定抗衡離子的位置

對(duì)于一些具有較高對(duì)稱(chēng)性或衍射點(diǎn)較弱的離子型MOFs 材料，其框架內(nèi)部的抗衡離子會(huì)遇到高度無(wú)序或衍射峰很多的情況，不能通過(guò)單晶解析確定其準(zhǔn)確位置，這時(shí)就需要利用MS 里的DMol3 或Sorption 模塊，通過(guò)幾何優(yōu)化（Geometry Optimization）或位置確定（Locate）的方法確定抗衡離子的位置。然后將得到結(jié)構(gòu)模型利用Sorption 模塊計(jì)算氮?dú)馕降葴鼐€(xiàn)與BET 表征的結(jié)果相對(duì)比來(lái)考察結(jié)構(gòu)模型的準(zhǔn)確性。CHEN D M 等利用MS 中的DMol3 模塊確定了兩個(gè)具有類(lèi)型結(jié)構(gòu)但帶有相反電荷的MOFs 材料中的抗衡離子位置，并利用Sorption 模塊研究了其吸附行為的差異[6]。

1. 3 確定MOFs 原子的電荷

確定MOFs 的原子電荷是決定計(jì)算模擬準(zhǔn)確性最為關(guān)鍵的一步。由于MOFs 材料與氣體分子之間主要依靠氫鍵作用、偶極-偶極相互作用或π-π作用等方式產(chǎn)生作用力，而這些作用力的強(qiáng)弱恰好與MOFs 材料中每個(gè)原子所攜帶的電荷密切相關(guān)。ZHOU D D 等發(fā)現(xiàn)利用DFT 計(jì)算得到的Mulliken 電荷并不能準(zhǔn)確地反映不飽和銅金屬位點(diǎn)與CO2的相互作用[5]。目前，文獻(xiàn)上報(bào)道的MOFs 材料原子電荷的計(jì)算方法主要有以下幾種。

（1）電荷平衡法（charge equilibration，Qeq）。該方法是目前文獻(xiàn)上采用較多的一種計(jì)算方法，其具有計(jì)算時(shí)間短、計(jì)算方法簡(jiǎn)單、適用范圍較廣等特點(diǎn)。CUI P 等人利用Qeq法賦予一例富氮MOF 材料的框架電荷并研究其對(duì)CO2 的選擇性吸附行為[7]；在最近報(bào)道的一篇文獻(xiàn)中，REN J 等人利用Qeq 法賦予MOF 框架電荷，并利用MS 軟件研究其對(duì)SF6 的吸附行為[8]。（2）簇基模型法。該方法首先將三維周期性的MOFs結(jié)構(gòu)轉(zhuǎn)化為若干個(gè)簇基結(jié)構(gòu)重復(fù)單元，然后再利用高斯軟件或MS 里的DMol3 模塊分別計(jì)算每個(gè)單元中原子的電荷。LIU K 等人將一例富氮MOF 材料分為3 個(gè)獨(dú)立的片段，并利用高斯軟件計(jì)算了每個(gè)片段上原子電荷來(lái)賦予整個(gè)MOF 框架電荷[9]。（3）原子連接貢獻(xiàn)法（Connectivity-Based Atom Contribution，CBAC）。該方法目前已被證實(shí)可以準(zhǔn)確地依據(jù)原子的連接類(lèi)型來(lái)賦予其電荷[10]。

2 吸附位點(diǎn)及吸附能計(jì)算

初始載荷的吸附位點(diǎn)反映的是MOFs 框架內(nèi)有利于束縛氣體分子的位置，一般通過(guò)氫鍵作用、偶極-偶極相互作用或π-π作用等方式表現(xiàn)。相互作用力越強(qiáng)，MOFs 對(duì)氣體的吸附能力也就越強(qiáng)，對(duì)應(yīng)的吸附能也就越高。吸附位點(diǎn)的計(jì)算主要是利用Sorption 模塊的固定載荷（Fixed loading）方法。吸附能的計(jì)算主要采用量子力學(xué)計(jì)算模塊DMol3 中的Energy。計(jì)算方法為用體系（通過(guò)Fixed loading 方法計(jì)算得到的最優(yōu)吸附模型）的總能量減去MOF 框架的能量再減去氣體分子的能量就得到了吸附能。CHEN D M 等人利用MS 軟件中Sorption 模塊中的Fixed loading 得到了乙炔分子的最優(yōu)吸附位置。計(jì)算結(jié)果表明，MOF 框架內(nèi)的F 原子可以與乙炔上的酸性H 原子之間有效地形成較強(qiáng)的氫鍵，其氫鍵的距離甚至小于F 原子的半徑距離，而乙烯分子的最優(yōu)吸附位點(diǎn)遠(yuǎn)離氟原子分布的位置。通過(guò)計(jì)算吸附能發(fā)現(xiàn)乙炔的吸附能是乙烯吸附能的兩倍多，這也解釋了該MOF 材料選擇性吸附乙炔的原因[11]。

當(dāng)載荷量逐漸提高時(shí)，吸附位點(diǎn)可以通過(guò)吸附密度圖的形式表現(xiàn)。例如：CHEN D M 等人通過(guò)Sorption模塊的固定壓力（Fixed pressure）方法計(jì)算了兩個(gè)具有MIL-88 結(jié)構(gòu)并修飾有不同功能基元的MOFs 材料的乙炔吸附密度圖[12]。計(jì)算結(jié)果表明，乙炔分子在MOF1中主要分布在雙核金屬位點(diǎn)周?chē)c三角雙錐的籠內(nèi)，而MOF2 中乙炔主要分布在富氮的分隔孔與三角雙錐的籠內(nèi)。考慮到雙核金屬位點(diǎn)之間的距離不能容納下兩個(gè)乙炔分子，因此其對(duì)乙炔的吸附量貢獻(xiàn)很低，這也解釋了MOF2 比MOF1 具有更高乙炔吸附量的原因。FAN W 等人利用Sorption 模塊計(jì)算了MOF 框架中丙烯的吸附密度圖，計(jì)算結(jié)果顯示，在較低壓力的時(shí)候丙烯主要吸附在苯環(huán)與吡嗪環(huán)周?chē)?，因此框架與氣體分子間的π-π堆積作用在較低壓力時(shí)對(duì)丙烯的吸附做出了主要貢獻(xiàn)[13]。

當(dāng)壓力繼續(xù)升高時(shí)，丙烯分子開(kāi)始填充其他的區(qū)域。而乙烯分子則需要在更高的壓力才能填充滿(mǎn)富含苯環(huán)與吡嗪環(huán)的區(qū)域，這就解釋了丙烯的吸附等溫線(xiàn)在低壓區(qū)域內(nèi)具有比乙烯吸附等溫線(xiàn)更陡的上升趨勢(shì)。除了計(jì)算單組分氣體之外，MS 軟件還可以用來(lái)計(jì)算雙組分混合氣體在特定條件下的分布密度圖。例如：YANG L等人計(jì)算了一例多孔Mn-MOF 在298 K、1 bar 條件下混合組分乙烯/乙烷的吸附瞬時(shí)分布圖，根據(jù)計(jì)算結(jié)果顯示，該MOF 可以吸附更多的乙烷氣體，并且乙烷氣體主要分布在吡啶氮與酰亞胺氧的周?chē)?，可以與框架之間形成較強(qiáng)的氫鍵相互作用，而乙烯氣體的周?chē)鷦t沒(méi)有能與之形成作用力的極性基團(tuán)[14]。

3 吸附等溫線(xiàn)計(jì)算

MS 軟件中Sorption 模塊的Adsorption isothem 任務(wù)可以用于MOFs 材料各種分子吸附等溫線(xiàn)的計(jì)算，例如氣體分子氮?dú)?、乙炔、乙烯，液體分子乙醇甲醇和一些常見(jiàn)的固體藥物分子5-氟尿嘧啶等。需要注意的是，計(jì)算曲線(xiàn)的平滑度與所設(shè)置的平衡步數(shù)以及生產(chǎn)步數(shù)有很大關(guān)聯(lián)，上述兩個(gè)參數(shù)設(shè)置的越大，曲線(xiàn)走勢(shì)越平滑，否則容易出現(xiàn)吸附點(diǎn)忽高忽低的情況。對(duì)于氣體分子的吸附等溫線(xiàn)，一般低壓區(qū)（P <0.1 bar）上升的趨勢(shì)要比高壓區(qū)（P >0.1 bar）更快，因此需要將logarithmic 的選項(xiàng)勾選，計(jì)算出來(lái)的吸附等溫線(xiàn)的線(xiàn)形會(huì)更貼近于實(shí)際測(cè)得的結(jié)果。需要值得注意的是，Adsorption isothem 任務(wù)中計(jì)算得到的吸附等溫線(xiàn)的橫坐標(biāo)為逸度，需要用Peng-Robinson 方程進(jìn)行壓力和逸度之間的轉(zhuǎn)換。利用MS 軟件計(jì)算77 K 下的氮?dú)馕降葴鼐€(xiàn)可以用來(lái)核實(shí)MOFs 材料是否已經(jīng)完全活化。例如：CHEN D M 等人利用MS 軟件模擬得到MOF 在77 K 下的吸附等溫線(xiàn)，根據(jù)計(jì)算結(jié)果顯示，其飽和吸附量為532 cm3/g，而實(shí)驗(yàn)測(cè)得的飽和吸附量為510 cm3/g，說(shuō)明該MOF 材料已被完全活化[15]。

4 結(jié)論與展望

利用計(jì)算模擬軟件MS 來(lái)研究晶態(tài)多孔材料的氣體吸附與分離行為，不僅可以對(duì)所合成得到的材料的吸附性能有初步的判斷，而且可以為材料的性能進(jìn)一步改善提供構(gòu)效關(guān)系研究基礎(chǔ)，為新材料的開(kāi)發(fā)設(shè)計(jì)提供理論依據(jù)。盡管MS 軟件在晶態(tài)多孔材料的吸附與分離行為研究上有著廣泛應(yīng)用，但該軟件還存在以下缺點(diǎn)：（1）MS 中Sorption 模塊不支持多核CPU 并行計(jì)算，導(dǎo)致其在設(shè)置較大平衡或生產(chǎn)步數(shù)參數(shù)后計(jì)算較慢，有時(shí)需要幾天甚至一周多來(lái)完成計(jì)算；（2）不支持主界面力場(chǎng)的直接修改，需要手動(dòng)修改力場(chǎng)文件中的參數(shù)來(lái)修正計(jì)算結(jié)果；（3）對(duì)于一些具有動(dòng)態(tài)吸附行為的MOFs 材料，MS 軟件目前不能很好地進(jìn)行適配。相信在不久的將來(lái)，隨著科學(xué)的發(fā)展與軟件技術(shù)的進(jìn)步，上述缺點(diǎn)有望得到進(jìn)一步解決。