吳　陽　岳麗麗　劉巧珍　王　霞

(遼寧大學(xué)化學(xué)院，沈陽110036)

?

咪唑醋酸離子液體熱力學(xué)性質(zhì)的理論研究

吳陽*岳麗麗劉巧珍王霞

(遼寧大學(xué)化學(xué)院，沈陽110036)

咪唑醋酸離子液體在催化、電化學(xué)、萃取等領(lǐng)域具有潛在的應(yīng)用價值，對其熱力學(xué)性質(zhì)的深入研究將為其應(yīng)用提供理論依據(jù)。本文采用密度泛函理論(DFT)方法和Born-Fajans-Haber(BFH)循環(huán)法對咪唑醋酸離子液體[Cnmim][OAc](n=1－6)進行熱力學(xué)性質(zhì)的理論研究。計算其相變過程中的解離焓、汽化焓、熔化焓、晶格焓、溶解焓等，并與已有實驗值進行比較。利用Gaussian 09程序在B3LYP/6-311+G(d,p)和M062X/TZVP兩種水平下計算解離焓值，同時通過計算得到分子體積和總氣相能的焓修正值，借助Matlab計算軟件擬合得到汽化焓值，取得與已有實驗值很好的一致性。使用Jenkins公式求得晶格能，計算得到晶格焓，最后利用BFH循環(huán)計算得到溶解焓。

密度泛函理論；功能化離子液體；BFH循環(huán)；解離焓；汽化焓；晶格焓；溶解焓

www.whxb.pku.edu.cn

1　引言

離子液體又稱室溫離子液體，是由特定陽離子和陰離子構(gòu)成的，離子液體的陽離子大多為季銨類離子，包括咪唑型、吡啶型、季銨型等。其中咪唑型陽離子構(gòu)成的離子液體因其易合成、易分離、功能性強等優(yōu)點得到了化學(xué)家們的廣泛關(guān)注，如[C2mim][PF6](1-乙基-3-甲基-咪唑六氟磷酸鹽離子液體)，[C2mim][OAc](1-乙基-3-甲基-咪唑醋酸鹽離子液體)和[C4mim][HSO4](1-丁基-3-甲基-咪唑硫酸氫鹽離子液體)等1,2。常用的陰離子有鹵素、氨基酸、四氟硼酸根等，如構(gòu)成[C2mim]Cl，[C2mim][Gly](1-乙基-3-甲基-咪唑甘氨酸鹽離子液體)，[C4mim][BF4](1-丁基-3-甲基-咪唑四氟硼酸離子液體)等3。在室溫下，離子液體可認為是一種有機液體物質(zhì)，具有蒸汽壓低、不易揮發(fā)、沸點高、導(dǎo)電性好等優(yōu)點，可循環(huán)使用，有更寬的電化學(xué)窗口，還有優(yōu)良的催化性能以及助催化性能等獨特性質(zhì)4－6。這使其在醫(yī)學(xué)、國防等眾多領(lǐng)域中都有良好的應(yīng)用前景。

2000年以后，離子液體發(fā)展的主流是新型功能性離子液體，研究者們聚焦于某一個行業(yè)或者某一個科研項目進行專一性離子液體的研發(fā)。最近MacFarlane等7將離子液體分為酸性、堿性和中性。堿性功能離子液體的陰離子基團為路易斯酸，例如CH3COO－、PhCOO－

等。其中具有堿性功能基團CH3COO－的咪唑醋酸系列離子液體，集合成簡易、蒸汽壓低、熱穩(wěn)定性好、溶解能力強等優(yōu)點于一身，隨著研究的深入，咪唑醋酸類堿性離子液體在更廣的領(lǐng)域得到了更多的應(yīng)用。相關(guān)報道中，劉曄等8于2007年報道的復(fù)合離子液體[C4mim][OAc]和[C4mim][TPPMS](1-丁基-3-甲基咪唑三苯基磷單磺酸鹽)，其中堿性離子液體[C4mim] [OAc]避免了反應(yīng)體系中粘稠固體物的形成，有利于鈀催化劑的穩(wěn)定性，且[C4mim][OAc]作為縛酸劑和配體[C4mim][TPPMS]綜合影響了鈀催化劑在Heck偶聯(lián)反應(yīng)中的活性。Chakraborti等9于2009年報導(dǎo)了[C4mim][OAc]離子液體在Boc2O(二碳酸二叔丁酯)與ArOH和Boc2O與H2O兩個相似的反應(yīng)體系中的催化作用機理，研究表明了咪唑醋酸離子液體陰陽離子各自發(fā)揮作用使Boc2O成為了親電親核雙活化試劑，并對[C4mim][HSO4]、[C4mim] [OAc]和[C4mim][BF4]這三種離子液體進行了對比，證明了堿性離子液體[C4mim][OAc]具有更好的催化作用。Wang等10于2010年報道了K3PO4、K2HPO4和K2CO3無機鹽水溶液分別與系列咪唑醋酸離子液體[C2mim][OAc]、[C4mim][OAc]和[C6mim] [OAc]混合，并將其應(yīng)用在L-色氨酸的提取與分離上，并通過實驗證明醋酸離子液體用于氨基酸的提取與分離具有明顯的優(yōu)勢。在高分子和生物化學(xué)領(lǐng)域，Rodriguez等11于2010年報道了[C2mim] [OAc]與特定的聚合物混溶能表現(xiàn)出更好的物理化學(xué)性質(zhì)。

物質(zhì)的性質(zhì)決定其功能，近年來離子液體性質(zhì)的預(yù)測方法發(fā)展備受關(guān)注12－14，從理論上建立體系的熱力學(xué)模型和預(yù)測體系的熱力學(xué)性質(zhì)數(shù)據(jù)對于指導(dǎo)科學(xué)研究和生產(chǎn)實踐過程發(fā)揮著重要的作用。本文選取[Cnmim][OAc](n=1－6)為研究體系，為了解其熱力學(xué)性質(zhì)，我們將在298.15 K的基準溫度下計算咪唑醋酸離子液體[Cnmim][OAc] (n=1－6)的解離焓、汽化焓、晶格焓、熔化焓和溶解焓。并與實驗及其它理論計算值進行比較。深入探討陰、陽離子對的熱力學(xué)性質(zhì)與宏觀性質(zhì)之間的關(guān)系，從而為離子液體的功能化設(shè)計提供理論指導(dǎo)，使其更好地應(yīng)用于環(huán)境、電化學(xué)、吸附等領(lǐng)域。

2　計算細節(jié)

本文采用密度泛函理論方法研究系列咪唑醋酸離子液體的離子對，由于B3LYP方法計算速度快，能節(jié)省大量的計算機資源，目前已廣泛應(yīng)用到離子液體體系的研究中，我們利用B3LYP/6-311+G(d,p)方法對咪唑醋酸離子液體[Cnmim] [OAc](n=1－6)氣態(tài)陰、陽離子及其離子對結(jié)構(gòu)進行幾何構(gòu)型的優(yōu)化，在優(yōu)化構(gòu)型基礎(chǔ)上進行振動頻率分析，結(jié)果表明均無虛頻存在。由于B3LYP方法在計算弱相互作用時也存在缺陷15，M062X方法能有效修正小分子雜化過程中的色散作用16，進而本文在B3LYP/6-311+G(d,p)水平優(yōu)化基礎(chǔ)上利用M062X/TZVP方法進行了單點能計算，給出結(jié)構(gòu)單元間的相互作用能量，并與B3LYP/6-311+G (d,p)方法所得值進行比較。所有的電子結(jié)構(gòu)計算工作利用Gaussian 0917程序完成。

由于功能性離子液體性質(zhì)的特殊性，研究者們對其熱力學(xué)性質(zhì)的研究與日俱增18。BFH循環(huán)廣泛應(yīng)用于熱力學(xué)性質(zhì)的許多領(lǐng)域，如化學(xué)反應(yīng)的原料和產(chǎn)物，或者是同一物質(zhì)的不同狀態(tài)的熱力學(xué)研究19－21。對于離子液體，一個合適的BFH循環(huán)可以使整個相變過程關(guān)聯(lián)起來，同時一些難以通過實驗測得的某些焓值可以利用已有的其他熱力學(xué)焓值得到，并對一些物理性質(zhì)做出預(yù)測。

式中任意三項已知，即可以求出第四項。這也體現(xiàn)了通過BFH循環(huán)的構(gòu)建，利用已有數(shù)據(jù)可以估算很多重要而又不易測得的數(shù)據(jù)。

3　結(jié)果和討論

3.1幾何構(gòu)型

Dupont22曾報道咪唑類離子液體中氫鍵的作用對于離子液體的結(jié)構(gòu)至關(guān)重要，進而影響離子液體的物化性質(zhì)，由此本文首先探討了咪唑醋酸離子液體[Cnmim][OAc](n=1－6)的氫鍵相互作用。我們在對咪唑氨基酸離子液體的氫鍵相互作用的研究中發(fā)現(xiàn)，[C2mim]+咪唑環(huán)周圍有S1、S2、S3、S4和S5五個區(qū)域可以和氨基酸陰離子形成氫鍵，其中S1區(qū)(C2―H與咪唑上甲基之間的區(qū)域) 和S2區(qū)(C2―H與咪唑上較長烷基側(cè)鏈之間的區(qū)域)與陰離子形成的氫鍵最穩(wěn)定23。這樣我們利用Gaussian 09程序，在B3LYP/6-311+G(d,p)水平下對醋酸系列離子液體的氣態(tài)陰陽離子對進行了優(yōu)化，在其勢能面上的最穩(wěn)定幾何結(jié)構(gòu)也是在相應(yīng)的S1、S2區(qū)域形成了穩(wěn)定的氫鍵(圖2)，這不但與氨基酸離子液體性質(zhì)相似，也與文獻報道的[C2mim]Cl體系計算的結(jié)構(gòu)相一致24。從幾何構(gòu)型中的氫鍵鍵長參數(shù)可以看出，陰、陽離子對之間都形成了2個氫鍵，其中一個氫鍵由C2―H與[OAc]－

圖1　離子液體的BFH循環(huán)Fig.1　Born-Fajans-Haber(BFH)cycle of ionic liquids

中的羰基氧形成，鍵長在0.1625－0.1645 nm之間。另一個氫鍵由烷基側(cè)鏈氫與[OAc]－中的羰基氧形成，鍵長在0.2024－0.2209 nm之間。相比較之下，C2―H形成的氫鍵相互作用更大，也更穩(wěn)定，毫無疑問，這是由于C2―H具有很強的酸性24。從計算結(jié)果可以看出，隨著咪唑烷基支鏈的增長，兩個氫鍵的鍵長基本上呈逐漸增大趨勢，并都小于O原子和H原子的van der Waals半徑之和(0.272 nm)。氨基酸離子液體[C2mim][Gly]最穩(wěn)定構(gòu)型的兩個氫鍵分別為0.1645和0.2039 nm23，傳統(tǒng)離子液體[C2mim][PF6]的兩個氫鍵分別為0.225和0.245 nm25，相比較而言，咪唑類陽離子與醋酸陰離子之間的氫鍵較短，這與醋酸陰離子中羰基氧帶有較大的負電荷有關(guān)，所以咪唑醋酸類離子液體應(yīng)該具有較高的熔點、較大的粘度及較小的電導(dǎo)率等性質(zhì)。

3.2解離焓

解離焓(ΔdisH)是氣相陰陽離子對[AB]g斷裂形成孤立氣相陰陽離子B－g和A+g過程中的能量變化。測定離子液體解離焓熱力學(xué)數(shù)據(jù)的實驗方法主要包括自由基動力學(xué)方法26和液相條件下測定的方法等27,28。但由于自由基壽命較短、溶劑化作用，以及實驗方法本身的缺陷等問題，使得實驗測定解離焓的方法仍然存在諸多困難。目前量子化學(xué)計算已被證明是獲得解離焓的最有效方法29－31。

利用Gaussian 09程序，在B3LYP/6-311+G(d, p)水平下對[Cnmim][OAc](n=1－6)氣態(tài)陰陽離子及其離子對結(jié)構(gòu)進行幾何構(gòu)型的優(yōu)化，得到相應(yīng)的焓值，再利用如下(2)式計算解離焓。式中，ΔdisH為咪唑醋酸離子液體的解離焓，H[A]+和H[B]－分別為通過Gaussian 09程序計算出的氣態(tài)陽離子和陰離子的焓值，H[AB]為Gaussian 09程序計算出的氣態(tài)陰陽離子對的焓值?？紤]到離子液體之間弱相互作用力較強，在B3LYP方法的基礎(chǔ)上又用M062X/TZVP方法計算了解離焓。

圖2　B3LYP/6-311+G(d,p)水平上優(yōu)化得到的[Cnmim][OAc](n=1－6)離子對的幾何結(jié)構(gòu)Fig.2　Geometric structures of ion pairs[Cnmim][OAc](n=1－6)optimized at B3LYP/6-311+G(d,p)level

表1給出了B3LYP和M062X兩種方法計算得到的結(jié)果，由此可見，兩種方法所得數(shù)據(jù)的最大差值約30 kJ·mol－1，如[C6mim][OAc]的兩種解離焓值分別為411.62和442.34 kJ·mol－1。此外，Preiss 等32估算的[Cnmim][OAc](n=2,3)解離焓值均為426 kJ·mol－1，與我們的計算結(jié)果具有很好的一致性。傳統(tǒng)咪唑離子液體的解離焓32，如[Cnmim] [BF4](n=2,4)的解離焓均為372 kJ·mol－1，[Cnmim] [NTf2](n=2,4,6)(1-烷基-3-甲基咪唑雙三氟甲基磺酰胺型離子液體)的解離焓分別為344、351、358 kJ·mol－1，[Cnmim]Cl(n=2,3,4)的解離焓分別為394、395、394 kJ·mol－1，由此可見，咪唑醋酸離子液體的解離焓更大，毫無疑問，這也與醋酸陰離子中羰基氧帶有較大的負電荷密切相關(guān)，也進一步證明咪唑類陽離子與醋酸陰離子之間存在較大的氫鍵相互作用。同時，具有相同陰離子的離子液體的解離焓與陽離子烷基鏈長度有關(guān)聯(lián)33，從表1中可以看出，當陽離子烷基側(cè)鏈中碳原子數(shù)少時，隨著烷基鏈的增長，解離焓數(shù)值呈現(xiàn)減小趨勢，但當烷基鏈中碳原子數(shù)增加到5以上時，解離焓值明顯增大，M062X/TZVP方法的計算結(jié)果尤為明顯，這與體系中形成了穩(wěn)定的納微結(jié)構(gòu)有關(guān)34。

表1　B3LYP/6-311+G(d,p)和M062X/TZVP水平下得到的[Cnmim][OAc](n=1－6)解離焓值Table 1　Dissociation enthalpies of[Cnmim][OAc](n= 1－6)at B3LYP/6-311+G(d,p)and M062X/TZVPlevels

3.3汽化焓

汽化焓(ΔvapH)是液體化合物汽化所必需的焓。由于汽化焓不改變物質(zhì)的溫度，只改變物質(zhì)的相，所以又稱為汽化潛熱。根據(jù)物質(zhì)分子運動論的觀點，物質(zhì)分子之間存在較強的吸引力，物質(zhì)從液相轉(zhuǎn)變?yōu)闅庀啵仨毧朔肿娱g的引力而做內(nèi)功，而物質(zhì)從液相變?yōu)闅庀鄷r，又必須反抗大氣壓力而做外功，液體蒸發(fā)或沸騰時需要的能量就是液體有汽化焓的原因。自從離子液體在非常高的溫度和低的壓強下的蒸發(fā)現(xiàn)象被發(fā)現(xiàn)后，離子液體的汽化焓就備受關(guān)注。

ΔvapH可采用不同的實驗方法來確定，通過分子動力學(xué)模擬估算35,36。Deyko等37將汽化焓考慮為陽離子和陰離子之間庫侖力和范德華相互作用力的總和。Kabo等38通過實驗測定表面張力與分子體積進而計算汽化熱，然而這種方法只適用于個別離子液體。在本文中，我們利用Krossing經(jīng)驗方法模型32對咪唑醋酸離子液體的汽化焓進行計算。Krossing經(jīng)驗方法考慮了離子液體在汽化過程中影響熱力學(xué)作用的基本因素：汽化焓與分子表面積之間的相關(guān)性。另外，由于汽化過程中物質(zhì)在氣態(tài)和液態(tài)時的焓值發(fā)生了變化，必須要考慮額外的自由度。計算公式如下：

式中，Vm為氣相離子對體積，Hg*為總氣相能的焓修正值。我們首先采用Gaussian 09程序計算[Cnmim][OAc](n=1－6)咪唑醋酸離子液體的Vm和Hg*(表2)，再利用已報道過的[Cnmim][OAc](n=2, 4,6)咪唑醋酸離子液體的汽化焓實驗數(shù)據(jù)，借助Matlab軟件擬合，得到的適用于咪唑醋酸離子液體的參數(shù)，a=－886.4 kJ·mol－1·nm－2，b=0.2745，c=305.5 kJ·mol－1。最后計算得到[Cnmim][OAc] (n=1－6)離子液體的汽化焓(表2)。Gaussian 09程序計算得到的[Cnmim][OAc](n=1－6)分子體積和汽化焓值與實驗值具有高度一致性39,40。與[Cnmim] [NTf2](n=2,4,6)的實驗汽化焓(137、136、141 kJ·mol－1)32相比較，咪唑醋酸系列的汽化焓與之相近，如[Cnmim][OAc](n=2,4,6)的汽化焓計算結(jié)果分別為134.41、134.90、140.79 kJ·mol－1；與[Cnmim] [Lact](n=2,3,5)(1-烷基-3-甲基咪唑乳酸離子液體)的實驗汽化焓(138.8、140.8、145.4 kJ·mol－1)41相比較，咪唑醋酸系列離子液體的汽化焓較小，這可能是由于離子液體[Cnmim][OAc](n=1－6)的極性相對較小。

3.4熔化焓

熔化焓(ΔfusH)是物質(zhì)從各向異性的晶體相轉(zhuǎn)變?yōu)楦飨蛲缘囊后w相所需要的焓。由于離子液體的特殊性質(zhì)，咪唑陽離子型離子液體的熔化焓通常較小32，在9.5－30.7 kJ·mol－1范圍內(nèi)，顯而易見，如此小的數(shù)據(jù)使其在實驗中很難測定，另外水含量對熔化焓的影響也會導(dǎo)致較大的實驗誤差42。由此，對于不易測得的熔化焓，研究者們通常取用平均值21.4 kJ·mol－1。

3.5晶格焓

晶格焓(ΔlatH)是離子液體從理想有序的各向異性晶體相向離子無限分離的氣體相進行轉(zhuǎn)變所需要的焓。這個過程中短程或長程有序排列被徹底破壞。因此，晶格焓可被視為一個晶格總熱力學(xué)穩(wěn)定性的量度。晶格焓可以通過修正的價鍵理論(VBT)公式計算43，也可以利用生成焓通過實驗方法確定44或使用晶體結(jié)構(gòu)作為起始點進行模擬等45。在眾多方法中，Jenkins公式被認為是獲得離子液體晶格焓較好的方法45,46。計算公式如下：其中公式(4a)計算晶格能ΔlatU，a=234.6 nm·kJ· mol－1，b=103.8 kJ·mol－1，是經(jīng)驗參數(shù)，Vm為分子體積。公式(4b)計算晶格焓，其中單原子的n+= n－ =3，線性離子的n+=n－=5，非線性離子的n+= n－ =6，R為氣體普適常數(shù)，T為氣體的熱力學(xué)溫度。表3中列出了利用Jenkins公式計算得到的[Cnmim][OAc](n=1－6)的晶格焓數(shù)值，由此可見，咪唑醋酸系列離子液體的晶格焓與陽離子烷基鏈長度密切關(guān)聯(lián)，隨著陽離子烷基鏈的增長，晶格焓數(shù)值減小，這說明隨著陽離子體積增大，晶體結(jié)構(gòu)坍塌所需的能量逐漸減小。文獻32中利用BHF循環(huán)計算給出[C2mim][OAc]和[C3mim][OAc]的晶格焓值為493和480 kJ·mol－1，與我們的計算結(jié)果一致性較好。

3.6溶解焓

溶解焓(ΔsolH)是離子液體從氣相的兩個無限分離的離子變化為均一的液相離子對過程中的焓變，可以定義為分離的氣態(tài)單離子在本體液相中的溶劑化作用。它不能被直接測量，可使用COSMO、IPCM、SCI-PCM等溶劑模型47來進行計算，但是，這些模型計算離子液體溶解焓的精確度目前尚不明確48。我們采用BFH循環(huán)計算溶解焓，公式如下：

表2　[Cnmim][OAc](n=1－6)的Vm,Hg*和?vapH以及實驗值Table 2　Calculated Vm,Hg*,and?vapH of[Cnmim][OAc](n=1－6)and corresponding experimental values

表3　[Cnmim][OAc](n=1－6)的晶格焓和溶解焓值Table 3　Calculated lattice and solvation enthalpies of [Cnmim][OAc](n=1－6)

利用B3LYP/6-311+G(d,p)和M062X/TZVP兩種方法計算得到了解離焓，因此也相應(yīng)得到了兩種結(jié)果的溶解焓值，列于表3中。利用BFH循環(huán)計算得到的溶解焓與COSMO-RS方法32得到的[C2mim][OAc]和[C3mim][OAc]的溶解焓值(－559和－560 kJ·mol－1)相比，相差較小，這說明可以利用BFH循環(huán)計算一些難以用實驗方法測得的溶解焓。另外，從表3中可以看出，隨著咪唑烷基側(cè)鏈的增長，溶解焓的數(shù)值相差不大，但[C6mim][OAc]具有最大的溶解焓(絕對值)，這表明隨著碳鏈增長，碳鏈之間的van der Waals力增大，較難溶解。

4　結(jié)論

由于離子液體數(shù)目巨大、種類繁多，在應(yīng)用中經(jīng)常會遇到缺乏離子液體熱力學(xué)數(shù)據(jù)的問題。據(jù)此本文采用密度泛函理論方法和Born-Fajans-Haber循環(huán)綜合計算了咪唑醋酸離子液體[Cnmim] [OAc](n=1－6)相變過程中的解離焓、汽化焓、熔化焓、晶格焓、溶解焓等熱力學(xué)性質(zhì)。在B3LYP/ 6-311+G(d,p)水平下，優(yōu)化得到[Cnmim][OAc](n= 1－6)離子液體氣態(tài)陰、陽離子對的穩(wěn)定結(jié)構(gòu)。由于醋酸陰離子中羰基氧電負性較強，優(yōu)化得到的咪唑醋酸離子液體氫鍵較短，解離焓數(shù)值較大，這導(dǎo)致咪唑醋酸類離子液體具有較高的熔點、較大的粘度及較小的電導(dǎo)率等性質(zhì)，計算表明解離焓值與陽離子烷基鏈長相關(guān)。利用Krossing經(jīng)驗方法和Matlab擬合得到咪唑醋酸離子液體的汽化焓，與實驗值具有較好的一致性。同時，本文使用Jenkins公式計算得到咪唑醋酸系列離子液體的晶格焓，且晶格焓與陽離子烷基鏈長度密切關(guān)聯(lián)。隨著陽離子體積增大，晶體結(jié)構(gòu)不規(guī)則性增大，晶格焓數(shù)值逐漸減小。最后，利用BFH循環(huán)計算溶解焓，計算表明隨著碳鏈增長，碳鏈之間的van der Waals力增大，使得[C6mim][OAc]較難溶解，具有最大的溶解焓(絕對值)。本文系統(tǒng)計算了咪唑醋酸離子液體[Cnmim][OAc](n=1－6)的解離焓、汽化焓、晶格焓、溶解焓值等，這些熱力學(xué)數(shù)據(jù)將為離子液體的功能化設(shè)計提供更多的理論指導(dǎo)，使其更好地應(yīng)用于生物化工、環(huán)境等更多的領(lǐng)域。

References

(1)Zhang,B.J.;Chen,L.;Xie,F.W.;Li,X.X.;Truss,R.W.;

Halley,P.J.;Shamshina,J.L.;Rogers,R.D.;McNally,T.Phys.

Chem.Chem.Phys.2015,17(21),13860.doi:10.1039/

C5CP01176K

(2)Zhang,Q.B.;Hua,Y.X.Acta Phys.-Chim.Sin.2011,27(3), 655.[張啟波,華一新.物理化學(xué)學(xué)報,2011,27(3),655.] doi:10.3866/PKU.WHXB20110339

(3)Dong,L.L.;He,L.;Tao,G.H.;Huang,M.;Hu,C.W.J.Chem. Eng.Data 2013,58(5),1176.doi:10.1021/je301285c

(4)Giammanco,C.H.;Kramer,P.L.;Yamada,S.A.;Nishida,J.; Tamimi,A.;Fayer,M.D.J.Chem.Phys.2016,144(10), 104506/1.doi:10.1063/1.4943390

(5)Beigbaghlou,S.S.;Marjani,K.;Habibi,A.;Atghia,S.V.RSC Adv.2016,6(24),20306.doi:10.1039/C5RA25297K

(6)Endres,F.;Zein ElAbedin,S.Phys.Chem.Chem.Phys.2006,8 (18),2101.doi:10.1039/b600519p

(7)MacFarlane,D.R.;Pringle,J.M.;Johansson,K.M.;Forsyth,S. A.;Forsyth,M.Chem.Commun.2006,No.18,1905. doi:10.1039/b516961p

(8)Liu,Y.;Li,M.;Lu,Y.;Wu,H.H.;Gao,G.H.Chem.J.Chin. Univ.2007,28(4),723.[劉曄,李敏,路勇,吳海虹,高國華.高等學(xué)?；瘜W(xué)學(xué)報,2007,28(4),723.]

(9)Chakraborti,A.K.;Roy,S.R.J.Am.Chem.Soc.2009,131 (20),6902.doi:10.1021/ja900076a

(10)Li,Z.Y.;Pei,Y.C.;Liu,L.;Wang,J.J.J.Chem.Thermodyn. 2010,42(7),932.doi:10.1016/j.jct.2010.03.010

(11)Rodriguez,H.;Rogers,R.D.Fluid Phase Equilib.2010,294 (1－2),7.doi:10.1016/j.fluid.2009.12.036

(12)Huang,Y.;Zhang,X.P.;Zhao,Y.S.;Zeng,S.J.;Dong,H.F.; Zhang,S.J.Phys.Chem.Chem.Phys.2015,17(40),26918. doi:10.1039/C5CP03446A

(13)Hu,Y.F.;Zhang,X.M.;Qi,J.G.;Yin,L.Y.J.Chem.Phys. 2015,143(20),204501/1.doi:10.1063/1.4936282

劉鐵頭帶來了豬頭肉和花生米，還提了兩瓶酒。我也沒客氣，簡單拾掇了一下就和劉鐵頭對坐下了。劉鐵頭不會找我喝閑酒，他來應(yīng)該有事的，擱平時我肯定要多問幾句。不過今天沒有，我覺得我今天需要喝點酒，劉鐵頭的酒來得很及時。

(14)Sprenger,K.G.;Jaeger,V.W.;Pfaendtner,J.J.Phys.Chem.B2015,119(18),5882.doi:10.1021/acs.jpcb.5b00689

(15)Gao,H.Y.;Zhang,Y.;Wang,H.J.;Liu,J.H.;Chen,J.M. J.Phys.Chem.A 2010,114(37),10243.doi:10.1021/jp1047 75z

(16)Zhao,Y.;Truhlar,D.G.Theor.Chem.Acc.2008,120(1－3), 215.doi:10.1007/s00214-007-0310-x

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

(18)Cervinka,C.;Padua,A.A.H.;Fulem,M.J.Phys.Chem.B 2016,120(9),2362.doi:10.1021/acs.jpcb.5b11070

(19)Born,M.Verh.Dtsch.Phys.Ges.1919,21,679.

(20)Fajans,K.Verh.Dtsch.Phys.Ges.1919,21,714.

(21)Haber,F.Verh.Dtsch.Phys.Ges.1919,21,750.

(22)Dupont,J.J.Braz.Chem.Soc.2004,15(3),341.doi:10.1590/ S0103-50532004000300002

(23)Wu,Y.;Zhang,T.T.J.Phys.Chem.A 2009,113(46),12995. doi:10.1021/jp906465h

(24)Mou,Z.X.;Li,P.;Bu,Y.X.;Wang,W.H.;Shi,J.Y.;Song,R. J.Phys.Chem.B 2008,112(16),5088.doi:10.1021/jp711358r

(25)Matsumoto,K.;Hagiwara,R.;Mazej,Z.;Benkic,P.;Zemva,B. Solid State Sciences 2006,8(10),1250.doi:10.1016/j. solidstatesciences.2005.12.018

(26)Berkowitz,J.;Ellison,G.B.;Gutman,D.J.Phys.Chem.1994, 98(11),2744.doi:10.1021/j100062a009

(27)Bordwell,F.G.;Zhang,X.M.Accounts Chem.Res.1993,26 (9),510.doi:10.1021/ar00033a009

(28)Long,Z.;Zhou,X.B.;Shen,X.D.;Li,D.L.;Liang,D.Q.Ind. Eng.Chem.Res.2015,54(46),11701.doi:10.1021/acs. iecr.5b03480

(29)Chan,B.;Radom,L.J.Phys.Chem.A 2013,117(17),3666. doi:10.1021/jp401248r

(30)Akhgarnusch,A.;Hoeckendorf,R.F.;Beyer,M.K.J.Phys. Chem.A 2015,119(39),9978.doi:10.1021/acs.jpca.5b06 975

(31)Santos,R.C.;Agapito,F.;Goncalves,E.M.;Martinho,S.;Jose, A.;Borges dos Santos,R.M.J.Chem.Thermodyn.2013,61, 83.doi:10.1016/j.jct.2013.01.028

(32)Preiss,U.;Verevkin,S.P.;Koslowski,T.;Krossing,I.Chem. Eur.J.2011,17(23),6508.doi:10.1002/chem.201003150

(33)Izgorodina,E.I.;Bernard,U.L.;MacFarlane,D.R.J.Phys. Chem.A 2009,113(25),7064.doi:10.1021/jp8107649

(34)Merunka,D.;Peric,M.;Peric,M.J.Phys.Chem.B 2015,119 (7),3185.doi:10.1021/jp512487y

(35)Esperanca,J.M.S.S.;Canongia Lopes,J.N.;Tariq,M.;Santos, L.M.N.B.F.;Magee,J.W.;Rebelo,L.P.N.J.Chem.Eng. Data 2010,55(1),3.doi:10.1021/je900458w

(36)Koddermann,T.;Paschek,D.;Ludwig,R.Chem.Phys.Chem. 2008,9(4),549.doi:10.1002/cphc.200700814

(37)Deyko,A.;Lovelock,K.R.J.;Corfield,J.A.;Taylor,A.W.; Gooden,P.N.;Villar-Garcia,I.J.;Licence,P.;Jones,R.G.; Krasovskiy,V.G.;Chernikova,E.A.;Kustov,L.M.Phys. Chem.Chem.Phys.2009,11(38),8544.doi:10.1039/b90 8209c

(38)Zaitsau,D.H.;Kabo,G.J.;Strechan,A.A.;Paulechka,Y.U.; Tschersich,A.;Verevkin,S.P.;Heintz,A.J.Phys.Chem.A 2006,110(22),7303.doi:10.1021/jp060896f

(39)Ma,X.X.;Wei,J.;Zhang,Q.B.;Tian,F.;Feng,Y.Y.;Guan,W. Ind.Eng.Chem.Res.2013,52(27),9490.doi:10.1021/ ie401130d

(40)Wei,J.;Bu,X.X.;Guan,W.;Xing,N.N.;Fang,D.W.;Wu,Y. RSC Adv.2015,5(86),70333.doi:10.1039/C5RA126 26F

(41)Wei,J.;Li,Z.;Gu,C.;Pan,Y.;Xing,N.N.;Tong,J.;Guan,W. J.Therm.Anal.Calorim.2016,125(1),547.doi:10.1007/ s10973-016-5318-9

(42)Troncoso,J.;Cerdeirina,C.A.;Sanmamed,Y.A.;Romani,L.; Rebelo,L.P.N.J.Chem.Eng.Data 2006,51(5),1856. doi:10.1021/je060222y

(43)Mallouk,T.E.;Rosenthal,G.L.;Muller,G.;Brusasco,R.; Bartlett,N.Inorg.Chem.1984,23(20),3167.doi:10.1021/ic 00188a028

(44)Souza,A.G.;Melo,A.M.S.;Santos,L.C.R.;Espinola,J.G. P.;Oliveira,S.F.;Airoldi,C.Thermochim.Acta 1998,313(2), 175.doi:10.1016/S0040-6031(98)00244-5

(45)Jenkins,H.D.B.;Roobottom,H.K.;Passmore,J.;Glasser,L. Inorg.Chem.1999,38(16),3609.doi:10.1021/ic9812 961

(46)Jenkins,H.D.B.;Pratt,K.F.Proc.R.Soc.A 1977,356(1684), 115.doi:10.1098/rspa.1977.0124

(47)Foresman,J.B.;Keith,T.A.;Wiberg,K.B.;Snoonian,J.; Frisch,M.J.J.Phys.Chem.1996,100(40),16098. doi:10.1021/JP960488J

(48)Klamt,A.;Schueuermann,G.J.Chem.Soc.Perkin Trans.2 1993,No.5,799.doi:10.1039/P29930000799

A Theoretical Study of the Thermodynamic Properties of Imidazolium Acetate Ionic Liquids

WU Yang*YUE Li-LiLIU Qiao-ZhenWANG Xia
(College of Chemistry,Liaoning University,Shenyang 110036,P.R.China)

Imidazolium acetate ionic liquids have a wide range of applications in catalysis,electrochemistry and extraction,and an improved understanding of their thermodynamic properties will provide a theoretical basis for future usages in related fields.In this work,density functional theory(DFT)and the Born-Fajans-Haber(BFH) cycle were used to study the thermodynamic properties of imidazolium acetate ionic liquids[Cnmim][OAc](n= 1－6).The dissociation,vaporization,latticeand solvation phase transition enthalpies of these liquids were calculated and compared with the corresponding experimental values.The dissociation enthalpies were determined using DFT at the B3LYP/6-311+G(d,p)and M062X/TZVP levels,and the results were compared with literature values.In addition,enthalpy of vaporization values were obtained from molecular volumesin conjunction with enthalpic corrections to the total gas-phase energy and Matlab fitting parameters,and the results were in good agreement with experimental data.Lastly,using the Jenkins formula,the lattice energies and the corresponding lattice enthalpies were determined,while the solvation enthalpies were calculated by employing the BFH cycle.

Density functional theory;Functionalized ionic liquid;Born-Fajans-Haber cycle;Dissociation enthalpy;Vaporization enthalpy;Lattice enthalpy;Solvation enthalpy

March 29,2016;Revised:May 11,2016;Published on Web:May 12,2016.

O641

10.3866/PKU.WHXB201605123

*Corresponding author.Email:wuyang@lnu.deu.cn;Tel:+86-24-62207802.

The project was supported by the National Natural Science Foundation of China(21373104)and Program for Liaoning Excellent Talents in University,China(LJQ2015042).

國家自然科學(xué)基金(21373104)和遼寧省高等學(xué)校優(yōu)秀人才支持計劃(LJQ2015042)資助項目

?Editorial office ofActa Physico-Chimica Sinica

［Article］