李曉兵，崔現(xiàn)寶，馮天揚(yáng)，揭會(huì)民，徐 麗，林瑞榕

（天津大學(xué)化工學(xué)院，化學(xué)工程國(guó)家重點(diǎn)實(shí)驗(yàn)室，天津300072）

離子液體由有機(jī)陽(yáng)離子和無(wú)機(jī)或者有機(jī)陰離子組成，在100℃（大多數(shù)在常溫下）呈液態(tài)[1]。目前，常見(jiàn)的陽(yáng)離子為：N，N-二烷基咪唑離子、N-烷基吡啶離子、吡咯離子、季銨離子和季鏻陽(yáng)離子。陰離子的種類很多，常見(jiàn)的主要有：醋酸鹽、四氟硼酸鹽、六氟磷酸鹽和硫酸鹽等。離子液體具有幾乎可以忽略的飽和蒸汽壓、低熔點(diǎn)、低毒性、高密度和良好的溶解性等優(yōu)點(diǎn)，尤其是離子液體的可設(shè)計(jì)性，即通過(guò)組合不同的陰、陽(yáng)離子和支鏈，“按需定制”離子液體，吸引了很多研究者的注意[2-4]。而想要設(shè)計(jì)出具有特定性質(zhì)的離子液體，必須先要獲得離子液體的物性數(shù)據(jù)，通過(guò)結(jié)構(gòu)與其性質(zhì)的關(guān)系，進(jìn)而得到合適的結(jié)構(gòu)。

對(duì)于研究者來(lái)說(shuō)，測(cè)量如此龐大物系的物理化學(xué)性質(zhì)是一項(xiàng)耗時(shí)耗力的工作，所以開(kāi)發(fā)估算離子液體性質(zhì)的計(jì)算方法是很必要的。定量結(jié)構(gòu)-性質(zhì)關(guān)系（QSPR）是常用的定量地研究物質(zhì)結(jié)構(gòu)和其性質(zhì)關(guān)系的方法，關(guān)于離子液體的QSPR研究已經(jīng)涉及到了黏度、密度、無(wú)限稀釋活度系數(shù)和表面張力等方面，而對(duì)于光學(xué)性質(zhì)方面的研究較少[5-9]。

折射率是表征物質(zhì)光學(xué)性能的重要參數(shù)，通過(guò)研究離子液體的折射率可以了解其光學(xué)性能、確定其純度、計(jì)算混合物的濃度，而且，折射率與其它性質(zhì)，如密度、表面張力以及介電常數(shù)是相關(guān)的。因此，預(yù)測(cè)并設(shè)計(jì)出所需要的折射率的離子液體對(duì)于開(kāi)發(fā)新的光學(xué)材料等具有重要的意義。

目前，對(duì)于預(yù)測(cè)離子液體折射率（nD）的報(bào)道較少。Gardas和Coutinho[10]分析了24種咪唑類的離子液體，得到了含有10個(gè)基團(tuán)參數(shù)的基團(tuán)貢獻(xiàn)模型，其總的相對(duì)偏差僅為0.18%，但該模型沒(méi)能預(yù)測(cè)其它類型的離子液體，具有很大的局限性。Soriano等[11]和Freire等[12]使用了相同的方法，也僅僅預(yù)測(cè)了咪唑類離子液體的折射率。Ram jugernath[13]分析了97種離子液體，包含了50種陽(yáng)離子和33種陰離子，得到了含有17個(gè)基團(tuán)參數(shù)的基團(tuán)貢獻(xiàn)模型，但是基團(tuán)組合不夠靈活，不能為離子液體的設(shè)計(jì)提供方便。本論文中，搜集了115種離子液體在不同溫度下的1161個(gè)折射率數(shù)據(jù)，運(yùn)用基團(tuán)貢獻(xiàn)法，得到了39個(gè)基團(tuán)參數(shù)的模型，計(jì)算值與文獻(xiàn)值的平均相對(duì)偏差為0.452%，而且基團(tuán)之間組合靈活，為離子液體的結(jié)構(gòu)設(shè)計(jì)提供便利。

1 計(jì)算部分

1.1 實(shí)驗(yàn)數(shù)據(jù)的選取

選用的實(shí)驗(yàn)數(shù)據(jù)取自離子液體數(shù)據(jù)庫(kù)（NIST），共115種離子液體（可在本刊網(wǎng)站下載）。其中，陽(yáng)離子有5類：季銨陽(yáng)離子、季鏻陽(yáng)離子、咪唑陽(yáng)離子、吡啶陽(yáng)離子和吡咯陽(yáng)離子，陰離子有26類，詳見(jiàn)表1。

1.2 基團(tuán)的劃分

將115種離子液體的結(jié)構(gòu)按照特定的方式進(jìn)行拆分，然后歸類，共得到39個(gè)基團(tuán)，見(jiàn)表2。

如表2所示，基團(tuán)劃分考慮了陽(yáng)離子主體部分不同取代基對(duì)離子液體折射率的影響，將陽(yáng)離子劃分成主基團(tuán)、R1、R2、R3和 R45個(gè)部分，其中，主基團(tuán)是陽(yáng)離子的中心結(jié)構(gòu)，R1、R2、R3和 R4分別為支鏈結(jié)構(gòu)，依次表示第1、2、3和4個(gè)支鏈，目前搜集到的離子液體，最多含有4個(gè)支鏈。對(duì)于陽(yáng)離子含多個(gè)支鏈的離子液體，支鏈按照其復(fù)雜程度排序。

首先是含有烴基以外的官能團(tuán)的支鏈，其次是鏈長(zhǎng)較長(zhǎng)或者結(jié)構(gòu)較復(fù)雜的支鏈，例如同時(shí)含芐基與乙基2個(gè)支鏈的陽(yáng)離子，芐基比乙基復(fù)雜，因此，該陽(yáng)離子的支鏈劃分為：1個(gè) R1即 C7H7，1個(gè) R2即 C2H5，其中，R1又分為 C6H5和 CH2，R2又分為 1個(gè)CH3和1個(gè)CH2。2個(gè)部分中CH2的貢獻(xiàn)值是不同的，這樣劃分基團(tuán)雖然增加了復(fù)雜性，但是估算精度會(huì)有所提高。

為了更加清晰地說(shuō)明基團(tuán)的分類方法，以三乙基甲基四氟硼酸銨和1-芐基-3甲基咪唑氯鹽為例。前者的陰離子可以劃分為1個(gè)B和4個(gè)F，陽(yáng)離子中包含1個(gè)主基團(tuán)：N，1個(gè) R1：乙基，1個(gè) R2：乙基，1 個(gè) R3：乙基，1 個(gè) R4：甲基，其中 R1，R2，R3 又分別劃分為1個(gè)CH3和1個(gè)CH2。因此，陽(yáng)離子包含1個(gè)主基團(tuán)：N，R1中的1個(gè) CH3和1個(gè) CH2，R2中的1個(gè)CH3和1個(gè)CH2，R3中的1個(gè)CH3和1個(gè)CH2以及R4中的1個(gè) CH3。后者的陰離子為1個(gè)Cl，陽(yáng)離子中包含主基團(tuán)：咪唑（Im），1個(gè) R1：芐基，1個(gè)R2：甲基，其中R1又可以劃分為：1個(gè)C6H5和1個(gè)CH2。因此，陽(yáng)離子包含1個(gè)主基團(tuán)：Im，R1中的1個(gè) C6H5和1個(gè) CH2，R2中的1個(gè) CH3。

1.3 基團(tuán)貢獻(xiàn)模型的建立

考慮到折射率隨溫度變化的關(guān)系，首先將搜集到的115種離子液體在不同溫度下折射率的變化作圖，為了看起來(lái)更加清晰，圖1中僅給出了11種離子液體折射率隨溫度的變化，剩余的104種離子液體具有同樣的變化趨勢(shì)。

由圖1可以看出，離子液體折射率隨著溫度的變化呈線性關(guān)系，且都隨著溫度的升高有所下降。然后，根據(jù)方程1分別對(duì)115種離子液體的折射率值進(jìn)行線性擬合，得到的相關(guān)系數(shù)都大于0.99。

表1 26類陰離子Table 126 k inds of anions

其中，A是離子液體在0 K下的折射率，B是折射率隨溫度變化的系數(shù)。根據(jù) Gardas[10]提出的基團(tuán)貢獻(xiàn)方法，A和B可以表示為：

其中，k為不同基團(tuán)的總數(shù)，值為39；ni為基團(tuán)i的個(gè)數(shù)；ai，bi分別為基團(tuán)i對(duì)參數(shù)A，B的貢獻(xiàn)值。

將115種離子液體的1161個(gè)折射率值用上述模型進(jìn)行回歸，得到各個(gè)基團(tuán)的貢獻(xiàn)值見(jiàn)表2。

表2 基團(tuán)的分類及其貢獻(xiàn)值Table 2 The groups and their con tribu tion values for estim ation

2 結(jié)果與討論

2.1 標(biāo)準(zhǔn)差和相對(duì)偏差的計(jì)算

為說(shuō)明此方法的回歸效果，用式（4）計(jì)算了估算的標(biāo)準(zhǔn)差，式中n表示數(shù)據(jù)點(diǎn)數(shù)1161。結(jié)果顯示參與估算的折射率數(shù)據(jù)的標(biāo)準(zhǔn)方差為0.00897。

圖1 離子液體的折射率隨溫度的變化Fig.1 The refractive indices of ionic liquids changing with tem perature

用回歸得到的基團(tuán)的貢獻(xiàn)值，計(jì)算參與估算的1161個(gè)折射率值，用式（5）計(jì)算文獻(xiàn)值與計(jì)算值的相對(duì)偏差（RD），進(jìn)而可以求出計(jì)算的平均相對(duì)偏差（ARD）為0.452%。具體的相對(duì)偏差的分布情況見(jiàn)圖2。

圖2 離子液體折射率計(jì)算值與文獻(xiàn)值的相對(duì)偏差的分布比例Fig.2 The d istribu tion of relative deviations between calcu lated and the literature values

由圖2可以看出，97.93%的數(shù)據(jù)點(diǎn)的相對(duì)偏差小于2%，說(shuō)明該方法是比較精確的。另外有2.07%數(shù)據(jù)點(diǎn)的相對(duì)偏差大于2%，這部分?jǐn)?shù)據(jù)是由 Ziyada[14-16]測(cè) 量 的[CNC2Oim]Cl（RD=2.7346%）、[CNC2Oim]TFMS（RD=2.068%）和[CNC2Him]TFMS（RD=2.032%），測(cè)得數(shù)據(jù)的不確定度相對(duì)較小，分別是0.00057、00006和 0.00081，所以應(yīng)該不是數(shù)據(jù)測(cè)量不準(zhǔn)確，而有可能是因?yàn)殡x子液體中雜質(zhì)含量比較高：含水量分別為0.0259%、0.0179%和0.0243%，另還含有一定的Br雜質(zhì)。

2.2 基團(tuán)貢獻(xiàn)法驗(yàn)證

為了更好的檢驗(yàn)估算方法的準(zhǔn)確性，利用該方法估算未參與回歸的離子液體[EM2N（CH2）2OH][ESO4][17]和[P666（14）][Deca][18]的折射率，結(jié)果見(jiàn)表3。得到的2種具有代表性的離子液體折射率的估算值與文獻(xiàn)實(shí)驗(yàn)值的相對(duì)偏差分別為0.520%和0.284%，平均相對(duì)偏差為0.331%，說(shuō)明該方法是可靠的。但是由于所搜集的離子液體的陽(yáng)離子中很少具有3個(gè)支鏈，所以該方法對(duì)于估算陽(yáng)離子具有3個(gè)支鏈的離子液體的折射率有一定的局限性。

表3 離子液體的折射率估算結(jié)果和文獻(xiàn)實(shí)驗(yàn)測(cè)量結(jié)果對(duì)比Table 3 Com parison of estim ated resu lts and literature values of refractive indices

2.3 與文獻(xiàn)中的基團(tuán)貢獻(xiàn)法相比較

Freire等[12]之前的文獻(xiàn)中僅僅提出了預(yù)測(cè)咪唑類離子液體折射率的基團(tuán)貢獻(xiàn)法；Ram jugernath[13]的文中分析了97種離子液體，得到了17個(gè)基團(tuán)參數(shù)模型；而本論文中搜集了115種離子液體，通過(guò)新的基團(tuán)劃分方法得的39個(gè)基團(tuán)參數(shù)更能靈活的組合成新型的離子液體，而且該模型預(yù)測(cè)的離子液體折射率的平均相對(duì)偏差僅為0.331%，小于 Ramjugernath[13]提出的17參數(shù)模型的平均相對(duì)偏差為0.45%，說(shuō)明該方法更為靈活、準(zhǔn)確。

3 結(jié)論

提出一種新型的基團(tuán)分類方法，通過(guò)回歸115種離子液體的1161個(gè)折射率數(shù)據(jù)點(diǎn)，建立了1個(gè)含有39個(gè)基團(tuán)參數(shù)的基團(tuán)貢獻(xiàn)模型來(lái)預(yù)測(cè)離子液體的折射率，回歸標(biāo)準(zhǔn)差和相對(duì)偏差的值分別為0.00897和0.452%。使用該方法計(jì)算了未參與回歸的2種離子液體的5個(gè)數(shù)據(jù)點(diǎn)，得到的相對(duì)偏差分別是0.520%和0.284%，說(shuō)明該基團(tuán)貢獻(xiàn)模型和該基團(tuán)分類方法是準(zhǔn)確、可靠的。

[1] 蔣鵬忠.離子液體的應(yīng)用與發(fā)展前景[J].科教文匯（中旬刊），2008，（3）：195-196 Jiang Pengzhong.The application and development prospects of ionic liquids[J].Education Science&Culture Magazine，2008，（3）： 195-196（in Chinese）

[2] Brennecke J F，Maginn E J.Ionic liquids： Innovative fluids for chemical processing[J].AIChE Journal，2001，47（11）： 2384-2389

[3] 彭艷枚，崔現(xiàn)寶，張纓，等.甲醇-乙酸甲酯-1-丁基-3-甲基咪唑雙三氟甲磺酰亞胺鹽的等壓汽液平衡[J].化學(xué)工業(yè)與工程，2013，30（6）：27-31 Peng Yanmei，Cui Xianbao，Zhang Ying，et al.Isobaric vapor-liquid equilibrium for methanol-methyl acetate-1-butyl-3-methyl-imidazoliumbis[（trifluoromethyl） sulfonyl]im ide at 101.3 kPa[J].Chem ical Industry and Engineering，2013，30（6）： 27-31（in Chinese）

[4] 關(guān)淑霞，于海南，邰永娜.離子液體的應(yīng)用進(jìn)展[J].化學(xué)工業(yè)與工程技術(shù)，2011，32（6）：32-36 Guan Shuxia，Yu Hainan，Tai Yongna.Application progress of ionic liquids[J].Chemical Industry and Engineering，2011，32（6）： 32-36（in Chinese）

[5] 盛衛(wèi)心，戎宗明，英徐根.基團(tuán)貢獻(xiàn)法分子設(shè)計(jì)研究的進(jìn)展[J].化學(xué)工業(yè)與工程，2007，24（5）：457-465 Sheng Weixin，Rong Zongm ing，Ying Xugen.Progress in computer-aided molecular design based on group contribution method[J].Chem ical Industry and Engineering，2007，24（5）： 457-465（in Chinese）

[6] 文路.離子液體粘度的QSPR研究[D].北京：北京化工大學(xué)，2011 Wen Lu.QSPR study on the viscosity of ionic liquids[D].Beijing：Beijing University of Chemical Technology，2011（in Chinese）

[7] Qiao Y，Ma Y，Huo Y，et al.A group contribution method to estimate the densities of ionic liquids[J].The Journal of Chemical Thermodynamics，2010，42（7）： 852-855

[8] 李瑞.離子液體萃取精餾的萃取劑分子設(shè)計(jì)[D].天津：天津大學(xué)，2012 Li Rui.Molecular design of ionic liquid entrainers for extractive distillation[D].Tianjin： Tianjin University，2012（in Chinese）

[9] 馬沛生，李平，趙興民，等.用基團(tuán)貢獻(xiàn)方法估算20℃下表面張力[J].化工學(xué)報(bào)，1994，45（1）：112-116 Ma Peisheng，Li Ping，Zhao Xingmin，et al.Estimating surface tensions by a group contribution method when the temperature is less than 20℃[J].Journal of Chemical Industry and Engineering，1994，45（1）： 112-116（in Chinese）

[10] Gradas R L，Coutinho JA P.Group contribution methods for the prediction of thermophysical and transport properties of ionic liquids[J].AIChE Journal，2009，55（5）： 1274-1290

[11] Soriano A N，Doma Jr B T，Li M.Density and refractive indexmeasurements of1-ethyl-3-methylimidazoliumbased ionic liquids[J].Journal of the Taiwan Institute of Chemical Engineers，2010，41 （1）： 115-121

[12] Freire M G，Teles A R R，Rocha M A A，et al.Thermophysical characterization of ionic liquids able to dissolve biomass[J].Journal of Chem ical&Engineering Data，2011，56（12）： 4813-4822

[13] Sattari M，Kamari A，Mohammadi A H，et al.A group contribution method for estimating the refractive indices of ionic liquids[J].Journal of Molecular Liquids，2014，200，Part B：410-415

[14] Ziyada A K，W ilfred C D，Murugesan T.Densities，viscosities and refractive indices of 1-alkyl-3-propanenitrile im idazolium chloride ionic liquids[J].Physics and Chemistry of Liquids，2011，50（2）： 152-160

[15] Ziyada A K，W ilfred C D.Effect of temperature and anion on densities，viscosities，and refractive indices of 1-octyl-3-propanenitrile imidazolium-based ionic liquids[J].Journal of Chemical&Engineering Data，2014，59（5）： 1385-1390

[16] Ziyada A K，Bustam M A，W ilfred C D，et al.Densities，viscosities，and refractive indices of 1-hexyl-3-propanenitrile imidazolium ionic liquids incorporated with sulfonate-based anions[J].Journal of Chemical&Engineering Data，2011，56（5）： 2343-2348

[17] Mahrova M，Vilas M，Domínguez A，et al.Physicochemical characterization of new sulfonate and sulfate ammonium ionic liquids[J].Journal of Chemical&Engineering Data，2012，57（2）： 241-248

[18] Almeida H F D，Lopes-da-Silva JA，F(xiàn)reire M G，et al.Surface tension and refractive index of pure and watersaturated tetradecyltrihexylphosphonium-based ionic liquids[J].The Journal of Chem ical Thermodynamics，2013，57（10）： 372-379