劉雅茹，劉寶友

（1.河北科技大學(xué)環(huán)境科學(xué)與工程學(xué)院，河北石家莊 050018；2.河北省污染防治生物技術(shù)實(shí)驗(yàn)室，河北石家莊 050018）

離子液體（ILs）作為一種新型綠色溶劑，已在能源、環(huán)境、材料等領(lǐng)域中展現(xiàn)出廣闊的應(yīng)用前景［1］。但I(xiàn)Ls的種類繁多，性質(zhì)多種多樣，而工業(yè)應(yīng)用過程中需要完整的數(shù)據(jù)資料，完全依賴實(shí)驗(yàn)測定性質(zhì)既不可能也沒必要。因此，開發(fā)合適的預(yù)測模型估算ILs的熱物理性質(zhì)顯得極為迫切。基團(tuán)貢獻(xiàn)法是預(yù)測化合物物化性質(zhì)的有效方法，被廣泛應(yīng)用于純物質(zhì)及混合體系的物性估算。該方法可以通過少量的實(shí)驗(yàn)數(shù)據(jù)來預(yù)測大量未知化合物的物性數(shù)值。相比于其他方法，基團(tuán)貢獻(xiàn)法具有預(yù)測ILs較為全面；預(yù)測過程較為簡單；能夠反向設(shè)計ILs；相對偏差較小等特點(diǎn)。本文結(jié)合近幾年的文獻(xiàn)報道，按ILs性質(zhì)及吸收酸性氣體進(jìn)行綜述，以期為新型ILs的設(shè)計及實(shí)現(xiàn)工業(yè)化應(yīng)用提供指導(dǎo)。

1 基團(tuán)貢獻(xiàn)法預(yù)測離子液體的物性

有關(guān)預(yù)測ILs物性的方法有很多，可分為實(shí)驗(yàn)?zāi)Ｐ秃屠碚撃Ｐ蛢深?。其中?shí)驗(yàn)類有基團(tuán)貢獻(xiàn)法、定量結(jié)構(gòu)性質(zhì)（QSPR）法等，理論類包括分子設(shè)計法、狀態(tài)方程（EoS）法、統(tǒng)計關(guān)聯(lián)式（Corr）法等［2］。QSPR 法很難為具體要求提供需要ILs的最佳性質(zhì)。分子設(shè)計法是復(fù)雜而漫長的計算過程［3－4］，因此，對研究人員來說不是一種合適的預(yù)測模型。EoS法僅適合于理想氣體的模型，所以，對實(shí)際氣體只能近似地應(yīng)用于低壓力范圍［5］。基團(tuán)貢獻(xiàn)法已應(yīng)用于各類ILs的研究，具有通用性強(qiáng)、預(yù)測過程簡單等特點(diǎn)。

1.1 密度

近幾年，國內(nèi)外文獻(xiàn)中報道了多種有關(guān)預(yù)測ILs密度的基團(tuán)貢獻(xiàn)法［6－14］。其預(yù)測方式由需要多個溫度點(diǎn)下的數(shù)據(jù)進(jìn)行預(yù)測發(fā)展到有一個溫度點(diǎn)就可以預(yù)測其他溫度下的ILs密度。

2010年，QIAO［11］等報道了一種預(yù)測ILs密度的基團(tuán)貢獻(xiàn)法，將離子液體按陰陽離子劃分。其中陽離子均含有一個中心和取代基，將取代基的烷基鏈切斷成單個烷基。他們共預(yù)測了123種純ILs和13種混合ILs（包含咪唑類、吡啶類、季銨鹽類等ILs），其平均偏差分別為0.88%和1.27%。再挑選188個數(shù)據(jù)點(diǎn)利用最小二乘法，計算了3 種典型的ILs的密度，其相對平均偏差為0.27%。其數(shù)學(xué)模型見式（1）：

式中：ρ表示密度（g／cm3）；P表示壓強(qiáng)（MPa）；T表示溫度（K），A，B，C是3 個參數(shù)（各個基團(tuán)貢獻(xiàn)值的總和）。與人工神經(jīng)網(wǎng)絡(luò)（ANN）相比，該模型的平均相對偏差較小、應(yīng)用范圍較廣。其不僅能預(yù)測純ILs的密度，也能預(yù)測二元混合ILs的密度。

2011年，SHAHBAZ［13］等提出了一種改進(jìn)的Lydersen－Joback－Reid基團(tuán)貢獻(xiàn)模型，該模型綜合了Lydersen模型對臨界壓力、體積及Joback－Reid模型對臨界溫度和沸點(diǎn)溫度估算精度較高的優(yōu)勢。用該模型預(yù)測了9 種低共熔離子液體在不同溫度（298.15～368.15 K）下的密度，其相對偏差為1.9%。該模型的優(yōu)點(diǎn)在于僅需測定一個溫度點(diǎn)下的密度，就可以預(yù)測其他各個溫度下的ILs密度，且計算精度較高。2014年，EVANGELISTA［14］等研究了一種新的基團(tuán)貢獻(xiàn)法（GCVOL）來預(yù)測純ILs的密度。相比與之前的研究，該研究更為廣泛。他們在溫度為251.62～473.15 K，壓力為0.1～300 MPa的條件下，選取了864種ILs共21 845個數(shù)據(jù)點(diǎn)進(jìn)行預(yù)測，其平均相對偏差為0.83%。該模型的預(yù)測精度較好，可用來預(yù)測ILs密度。

1.2 黏度

黏度的測量是實(shí)現(xiàn)安全生產(chǎn)、節(jié)約和開發(fā)能源的重要手段。其在理論研究中對ILs其他性質(zhì)（電導(dǎo)率等）也有影響。工業(yè)的發(fā)展和大量設(shè)備的使用，潤滑問題越來越受到關(guān)注，液體黏度研究受到人們的廣泛重視。已經(jīng)有文章［15－20］報道了用基團(tuán)貢獻(xiàn)法預(yù)測ILs的黏度。

2008年，GARDAS［18］等在溫度為293～393K和黏度為4～21 000Pa·S范圍內(nèi)，用基團(tuán)貢獻(xiàn)法預(yù)測了29種ILs（包含咪唑類、吡啶類、吡咯類）的黏度，對大約500個數(shù)據(jù)點(diǎn)進(jìn)行研究，結(jié)果表明平均相對偏差為7.78%。使用的Orrick－Erbar估算模型見式（2）：

式中：ρ和η分別表示ILs的密度和黏度；M和T分別代表相對分子質(zhì)量和溫度；N是阿伏伽德羅常數(shù)，V是分子體積；a，b，c是各個基團(tuán)的貢獻(xiàn)值；A，B是a，b貢獻(xiàn)值的總和。該模型的基團(tuán)劃分：將ILs中的陽離子中心與陰離子劃分為2 個基團(tuán)，剩下的烷基取代基切成甲基或亞甲基。

2012年，GHARAGHEIZI［19］等發(fā)現(xiàn)了一種在常壓下預(yù)測ILs 黏度的基團(tuán)貢獻(xiàn)法，比GARDAS［18］等的模型更有通用性。他們研究了443種離子液體的1 672 個具有黏度貢獻(xiàn)值的基團(tuán)數(shù)據(jù)點(diǎn)。為了區(qū)別陰、陽離子結(jié)構(gòu)對黏度的影響，故將陰陽離子劃分開。預(yù)測結(jié)果發(fā)現(xiàn)，該模型的平均相對偏差低于6.4%，優(yōu)于MOUSAVISAFAVI［20］等在2013年提出的QSPR 模型，其預(yù)測結(jié)果的平均相對偏差為9%。

2014年，BAJIC［21］等研究了預(yù)測ILs黏度的基團(tuán)貢獻(xiàn)法。與GARDAS［18］等和GHARAGHEIZI［19］等的研究相比，主要預(yù)測ILs與有機(jī)相混合體系的黏度。運(yùn)用兩種不同的基團(tuán)貢獻(xiàn)法UNIFAC－VISCO和ASOG－VISCO 對10種混合體系進(jìn)行預(yù)測，加入新影響參數(shù)的UNIFAC－VISCO 模型相對偏差較高（0.27%～11.76%），而ASOG－VISCO 模型預(yù)測結(jié)果比UNIFAC－VISCO 模型好。該模型可用于預(yù)測某些ILs的二元體系，對于其他體系的相對偏差在11%以上。

1.3 電導(dǎo)率

電導(dǎo)率是ILs的重要性質(zhì)之一，其可以幫助研究人員了解離子的結(jié)構(gòu)，以便控制設(shè)計工業(yè)進(jìn)程。在電化學(xué)領(lǐng)域方面，高電導(dǎo)率的ILs可以用來制作高性能的電池。隨著ILs的快速發(fā)展，其在許多科學(xué)領(lǐng)域得到廣泛應(yīng)用。因此，預(yù)測電導(dǎo)率對科學(xué)地設(shè)計ILs非常重要［21－24］。

2007年，MATSUDA［22］等提出了預(yù)測ILs電導(dǎo)率的基團(tuán)貢獻(xiàn)模型。這次研究在溫度為－30～70 ℃范圍內(nèi)，共選取206個數(shù)據(jù)點(diǎn)作為參數(shù)對ILs進(jìn)行預(yù)測，發(fā)現(xiàn)其實(shí)驗(yàn)值與計算值之間的相關(guān)系數(shù)R2為0.909 4。該模型主要依靠陽離子的類型、側(cè)鏈長度及陰離子類型，即可預(yù)測ILs的性質(zhì)。他們還利用該模型對ILs進(jìn)行反向設(shè)計，可根據(jù)需要，設(shè)計出具有理想電導(dǎo)率的ILs。

2009年，GARDAS［23］等應(yīng)用基團(tuán)貢獻(xiàn)法預(yù)測ILs的電導(dǎo)率。他們在溫度為258.10～433.15K，電導(dǎo)率為0.01～12.68S／m 進(jìn)行研究。利用300個數(shù)據(jù)點(diǎn)預(yù)測了15種ILs（包括咪唑類、吡啶類、吡咯類、季銨鹽類）的電導(dǎo)率。其數(shù)學(xué)模型（VTF 方程）見式（3）：

式中：λ是ILs的電導(dǎo)率；ni為基團(tuán)i出現(xiàn)的次數(shù)；k為基團(tuán)i出現(xiàn)的總數(shù)；ai，λ，bi，λ分別為對應(yīng)基團(tuán)的貢獻(xiàn)值；T0λ＝165.06 K。該模型與MARSUDA 提出的模型均適用于預(yù)測黏度。兩種模型相比，VTF 方程更為簡單，且精度高。2014年，GHARAGHEIZI［24］等將最小二乘向量機(jī)與基團(tuán)貢獻(xiàn)法聯(lián)用，選取1 077個實(shí)驗(yàn)數(shù)據(jù)點(diǎn)，對54種ILs的電導(dǎo)率進(jìn)行預(yù)測，其平均相對偏差低于3.3%。

表1列出了基團(tuán)貢獻(xiàn)法對ILs密度、黏度、電導(dǎo)率預(yù)測的一些參數(shù)及相對偏差?？梢钥闯?，基團(tuán)貢獻(xiàn)法對比于其他方法有明顯優(yōu)勢：預(yù)測范圍廣，適用條件要求低，相對偏差較?。搭A(yù)測結(jié)果準(zhǔn)確）。

表1 不同模型對ILs物性預(yù)測的典型對比Tab.1 A typical comparison of different models for the prediction of the physical properties of ILs

2 基團(tuán)貢獻(xiàn)法預(yù)測ILs對煙道氣體的吸收

煙道氣的主要成分為CO2，CO，NO，NO2，SO2，H2S等，對環(huán)境有較嚴(yán)重的破壞。ILs是一種環(huán)境友好的綠色材料，它對氣體有較好的分離吸收功能。因此，對煙道氣的治理提供了一種綠色環(huán)保的替代方案。預(yù)測煙道氣的主要成分在ILs中的溶解度，設(shè)計理想的吸收過程顯得尤為重要。截至目前，對于采用基團(tuán)貢獻(xiàn)法對H2S，NO2，CO，NO 的溶解性還未見有文獻(xiàn)報道，以下僅對CO2和SO2的預(yù)測進(jìn)行綜述。

2.1 對CO2 氣體吸收

CO2作為重要的溫室氣體，與其他氣體的分離在工業(yè)應(yīng)用上尤為重要。實(shí)驗(yàn)表明，ILs對CO2有較好的吸收分離效果。

2011年，伍艷輝［30］等利用基團(tuán)貢獻(xiàn)法預(yù)測了常壓下CO2在含極性基團(tuán)的咪唑類離子液體中的溶解度。通過計算得到亨利系數(shù)，進(jìn)而得到氣體在該類離子液體中的溶解選擇性。預(yù)測結(jié)果表明，極性基團(tuán)的引入提高了ILs對CO2的吸收，計算值與實(shí)驗(yàn)結(jié)果一致。其數(shù)學(xué)模型見式（4）和式（5）：

其中，A是取決于氣體和離子液體的種類，不隨體系溫度和壓力變化的參數(shù)，B是與氣體和離子液體的種類以及溫度有關(guān)的參數(shù)。在給定的溫度下，A，B均可以視為常數(shù)。δIL，δi分別為氣體和離子液體的溶解度參數(shù)。Hi是亨利常數(shù)，aij是氣體在離子液體中的溶解度，F(xiàn)j是基團(tuán)j的物質(zhì)的量吸引常數(shù)。該模型可以用于預(yù)測多種氣體在一種離子液體中的溶解選擇性，可根據(jù)基團(tuán)貢獻(xiàn)來設(shè)計離子液體，從而使其更好地應(yīng)用在工業(yè)中。

2012年，ASHRAFMANSOURI［31］等報道了一種估算氣體在ILs中的溶解度的基團(tuán)貢獻(xiàn)模型，該模型是基于統(tǒng)計締合流體理論（SAFT－γ）。在溫度為313.15～353.15K，壓力達(dá)10MPa時，該模型能較好地描述ILs與氣體的相行為。他們通過對參數(shù)的優(yōu)化設(shè)計預(yù)測了一些ILs與CO2系統(tǒng)，取得了較好的預(yù)測結(jié)果，其平均相對偏差為1.309 7%。該模型能夠很好地預(yù)測這種系統(tǒng)的相行為，可用于未來的ILs研究中。

2.2 對SO2 氣體吸收

SO2是大氣中數(shù)量最大的有害成分，80%的SO2來源于煙道氣排放。為了減少或消除其危害，開發(fā)了各種治理方法。其中將ILs吸收SO2應(yīng)用到工業(yè)中，實(shí)現(xiàn)清潔生產(chǎn)，具有廣泛的應(yīng)用前景［32－37］。

2010年，MARTIN［32］等應(yīng)用一種基團(tuán)貢獻(xiàn)狀態(tài)方程對二元體系（SO2＋I(xiàn)Ls）的熱力學(xué)模型進(jìn)行優(yōu)化。他們認(rèn)為氣體溶解度是吸引力相和自由體積相的加和。其中自由體積相由臨界直徑dc和標(biāo)準(zhǔn)范德華體積關(guān)聯(lián)計算，吸引力相由純ILs基團(tuán)常數(shù)、純基團(tuán)能參數(shù)和基團(tuán)相互影響參數(shù)決定。這些參數(shù)可從文獻(xiàn)［33］中查到，為了計算簡單，假設(shè)基團(tuán)間的相互影響參數(shù)為1或0。在溫度為283～323K 范圍內(nèi)，選取11個數(shù)據(jù)點(diǎn)對（SO2＋［bmpyrr］［Tf2N］）進(jìn)行預(yù)測，其平均相對壓力偏差為6.2%。同年，CARVALHO［34］等用Flory－Huggins模型對SO2在ILs中的溶解度進(jìn)行預(yù)測，預(yù)測結(jié)果出現(xiàn)正負(fù)2種偏差，在－0.042～0.059之間。兩種模型相比，基團(tuán)貢獻(xiàn)法實(shí)驗(yàn)次數(shù)少，更為簡單，且準(zhǔn)確度較高。

2013年，王衛(wèi)［35］擬合了ILs－SO2體系的UNFIAC模型。在298.15 K 下，對5 種離子液體與SO2用UNFIAC 模型進(jìn)行模擬，并與COSMO－RS模型以及COSMO－thermX 模型進(jìn)行對比，發(fā)現(xiàn)只有UNIFAC 模型有較好的實(shí)驗(yàn)結(jié)果，［EMIM］［BF4］的平均相對偏差為3.25%。在298.15K 下，不同壓力下不同ILs對SO2的吸收差別不大。其計算氣體溶解度模型見式（6）：

x1和x2分別代表氣體在液相和氣相中的物質(zhì)的量分?jǐn)?shù)；φ（T，P，x1）表示氣相中氣體的逸度系數(shù)；P是系統(tǒng)的壓力；Ps是氣體的飽和蒸汽壓，γ是氣體在液相中的活度系數(shù)，可以由基團(tuán)貢獻(xiàn)UNIFAC模型求得。

3 結(jié) 語

目前，已有大量的文獻(xiàn)報道采用基團(tuán)貢獻(xiàn)法對ILs性質(zhì)進(jìn)行預(yù)測，并取得了良好的預(yù)期效果。與其他方法相比，基團(tuán)貢獻(xiàn)法具有應(yīng)用范圍廣、與計算機(jī)聯(lián)用預(yù)測過程更簡單、通用性強(qiáng)等優(yōu)點(diǎn)?；鶊F(tuán)貢獻(xiàn)法的核心應(yīng)是基團(tuán)的劃分方式，因此，對基團(tuán)劃分的策略統(tǒng)一是未來研究的重點(diǎn)。對于同一性質(zhì)的預(yù)測模型較多，因此，建立在統(tǒng)一尺度上的模型及參數(shù)的選取，顯得尤為重要。模擬預(yù)測與實(shí)驗(yàn)驗(yàn)證的結(jié)合，才能更好地指導(dǎo)進(jìn)行ILs結(jié)構(gòu)設(shè)計，實(shí)現(xiàn)由“經(jīng)驗(yàn)研究”到“定向優(yōu)化”的研究飛躍。對煙道氣在ILs溶解度的預(yù)測研究還處于早期階段，研究設(shè)計基于基團(tuán)貢獻(xiàn)法的通用的氣體吸收模型，必將有利地推動ILs應(yīng)用于煙道氣治理的工業(yè)化進(jìn)程。

／References：

［1］ PENA－PEREIRA F，NAMIESNIK J.Ionic liquids and deep eutectic mixtures：Sustainable solvents for extraction processes［J］.Chem Sus Chem，2014，7（7）：1784－1800.

［2］ COUTINHO J A P，CARVALHO P J，OLIVEIRA N M C.Predictive methods for the estimation of thermophysical properties of ionic liquids［J］.RSC Advances，2012，2（19）：7322－7346.

［3］ MAGINN E J，ELLIOTT J R.Historical perspective and current outlook for molecular dynamics as a chemical engineering tool［J］.Industrial ＆ Engineering Chemistry Research，2010，49（7）：3059－3078.

［4］ RAI N，MAGINN E J.Vapor－liquid coexistence and critical behavior of ionic liquids via molecular simulations［J］.Journal of Physical Chemistry Letters，2011，2（12）：1439－1443.

［5］ VEGA L F，VILASECA O，LLOVELL F，et al.Modeling ionic liquids and the solubility of gases in them：Recent advances and perspectives［J］.Fluid Phase Equilibria，2010，294（1／2）：15－30.

［6］ LARRIBA M，GARCIA S，NAVARRO P，et al.Physical properties ofn－butylpyridinium tetrafluoroborate andn－butyl pyridinium bis（trifluoromethylsulfonyl）imide binary ionic liquid mixtures［J］.Journal of Chemical ＆ Engineering Data，2012，57（4）：1318－1325.

［7］ PADUSZYNSKI K，DOMANSKA U.A new group contribution method for prediction of density of pure ionic liquids over a wide range of temperature and pressure［J］.Industrial ＆Engineering Chemistry Research，2012，51（1）：591－604.

［8］ BHATTACHARJEE A，VARANDA C，F(xiàn)REIRE M G，et al.Density and viscosity data for binary mixtures of 1－alkyl－3－methylimidazolium alkylsulfates＋water［J］.Journal of Chemical ＆Engineering Data，2012，57（12）：3473－3482.

［9］ LAZZUS J A.A group contribution method to predictρ－T－Pof ionic liquids［J］.Chemical Engineering，2010，197（7）：974－1015.

［10］ 喬燕.離子液體／二氧化碳體系密度黏度和相平衡研究［D］.天津：天津大學(xué)，2010.QIAO Yan.Study on the Densities，Viscosities and Phase Equilibria of Ionic Liquids／Carbon Dioxide System［D］.Tianjin：Tianjin University，2010.

［11］ QIAO Yan，MA Youguang，HUO Yan，et al.A group contribution method to estimate the densities of ionic liquids［J］.Journal of Chemical Thermodynamics，2010，42（7）：852－855.

［12］ 沈沖.甲醇－乙醇－水－離子液體體系汽液平衡的測定及熱力學(xué)模型研究［D］.北京：北京化工大學(xué)，2011.SHEN Chong.Measurement and Thermodynamic Modeling of Vapor－Liquid Equilibria for Systems Containing Methanol，Ethanol，Water and Ionic Liquids［D］.Beijing：Beijing University of Chemical Technology，2011.

［13］ SHAHBAZ K，MJALLI F S，HASHIM M A，et al.Prediction of deep eutectic solvents densities at different temperatures［J］.Thermochimica Acta，2011，515（1／2）：67－72.

［14］ EVANGELISTA N S，do CARMO F R，de SANTIAGOAGUIAR R S，et al.Development of a new group contribution method based on GCVOL model for the estimation of pure ionic liquid density over a wide range of temperature and pressure［J］.Industrial ＆ Engineering Chemistry Research，2014，53（22）：9506－9512.

［15］ ADAMOVA G，GARDAS R L，NIEUWENHUYZEN M，et al.Alkyltributylphosphonium chloride ionic liquids：Synthesis，physicochemical properties and crystal structure［J］.Dalton Transactions，2012，41（27）：8316－8332.

［16］ ADAMOVA G，GARDAS R L，REBELO L P N，et al.Alkyltrioctylphosphonium chloride ionic liquids：Synthesis and physicochemical properties［J］.Dalton Transactions，2011，40（47）：12750－12764.

［17］ COSTA A J L，ESPERANCA J M S S，MARRUCHO I M，et al.Densities and viscosities of 1－ethyl－3－methylimidazoliumnalkyl sulfates［J］.Journal of Chemical ＆ Engineering Data，2011，56（8）：3433－3441.

［18］ GARDAS R L，COUTINHO J A P.A group contribution method for viscosity estimation of ionic liquids［J］.Fluid Phase Equilibria，2008，266（1／2）：195－201.

［19］ GHARAGHEIZI F，IIANI－KASHKOULI P，MOHAMMADI A H，et al.Development of a group contribution method for determination of viscosity of ionic liquids at atmospheric pressure［J］.Chemical Engineering Science，2012，80：326－333.

［20］ MOUSAVISAFAVI S，MIRKHANI S，GHARAGHEIZI F，et al.A predictive quantitative structure－property relationship for glass transition temperature of 1，3－dialkyl imidazolium ionic liquids（Part 1）：The linear approach［J］.Journal of Thermal Analysis and Calorimetry，2013，111（1）：235－246.

［21］ BAJIC D M，SERBANOVIC S P，ZIVKOVIC E M，et al.Prediction and correlation of viscosity of binary mixtures of ionic liquids with organic solvents［J］.Journal of Molecular Liquids，2014，197：1－6.

［22］ MATSUDA H，YAMAMOTO H，KURIHARA K，et al.Computer－aided reverse design for ionic liquids by QSPR usingdescriptors of group contribution type for ionic conductivities and viscosities［J］.Fluid Phase Equilibria，2007，261（1／2）：434－443.

［23］ GARDAS R L，COUTINHO J A P.Group contribution methods for the prediction of thermophysical and transport properties of ionic liquids［J］.AIChE Journal，2009，55（5）：1274－1290.

［24］ GHARAGHEIZI F，IIANI－KASHKOULI P，SATTARI M，et al.Development of a LSSVM－GC model for estimating the electrical conductivity of ionic liquids［J］.Chemical Engineering Research and Design，2014，92（1）：66－79.

［25］ JACQUEMIN J，GE R，NANCARROW P，et al.Prediction of ionic liquid properties（I）：Volumetric properties as a function of temperature at 0.1 MPa［J］.Journal of Chemical ＆Engineering Data，2008，53（3）：716－726.

［26］ VALDERRAMA J O，ZARRICUETA K.A simple and generalized model for predicting the density of ionic liquids［J］.Fluid Phase Equilibria，2009，275（2）：145－151.

［27］ LAZZUS J A.ρ（T，p）model for ionic liquids based on quantitative structure－property relationship calculations［J］.Journal of Physical Organic Chemistry，2009，22（12）：1193－1197.

［28］ KOLLER T，RAUSCH M H，SCHULZ P S，et al.Viscosity，interfacial tension，self－diffusion coefficient，density，and refractive index of the ionic liquid 1－ethyl－3－methylimidazolium tetracyanoborate as a function of temperature at atmospheric pressure［J］.Journal of Chemical ＆ Engineering Data，2012，57（3）：828－835.

［29］ ZHAO He，LIANG Zhongcheng，LI Fang.An improved model for the conductivity of room－temperature ionic liquids based on hole theory［J］.Journal of Molecular Liquids，2009，149（3）：55－59.

［30］ 伍艷輝，于世昆，段永超，等.基團(tuán)貢獻(xiàn)法預(yù)測含極性基團(tuán)離子液體對氣體的溶解選擇性［J］.化工學(xué)報，2011，62（10）：2684－2690.WU Yanhui，YU Shikun，DUAN Yongchao，et al.Determination of gas selectivity of solution in RTILs with polar group by group contribution method［J］.CIESC Journal，2011，62（10）：2684－2690.

［31］ ASHRAFMANSOURI S S，RAEISSI S.Modeling gas solubility in ionic liquids with the SAFT－γ group contribution method［J］.Journal of Supercritical Fluids，2012，63：81－91.

［32］ MARTIN A，MENDEZ D，BERMEJO M D.Application of a group contribution equation of state for the thermodynamic modeling of binary systems（gas＋ionic liquids）with bis［（trifluoromethyl）sulfonyl］imide anion［J］.Journal of Chemical Thermodynamics，2010，42（4）：524－529.

［33］ FORNARIA T.Revision and summary of the group contribution equation of state parameter table：Application to edible oil constituents［J］.Fluid Phase Equilibria，2007，262（1／2）：187－209.

［34］ CARVALHO P J，COUTINHO J A P.Non－ideality of solutions of NH3，SO2，and H2S in ionic liquids and the prediction of their solubilities using the Flory －Huggins model［J］.Energy Fuels，2010，24（12）：6662－6666.

［35］ 王衛(wèi).UNIFAC模型用于離子液體－氣體體系的研究［D］.北京：北京化工大學(xué)，2013.WANG Wei.UNIFAC Model for Ionic Liquid－Gas Systems［D］.Beijing：Beijing University of Chemical Technology，2013.

［36］ 楊會龍，劉寶友，王園園.氨基功能化離子液體表征及吸收SO2的實(shí)驗(yàn)研究［J］.河北科技大學(xué)學(xué)報，2011，32（3）：220－224.YANG Huilong，LIU Baoyou，WANG Yuanyuan.Characterization of amino－functionalized ionic liquids and their application in SO2absorption［J］.Journal of Hebei University of Science and Technology，2011，32（3）：220－224.

［37］ 章淼淼，任愛玲，關(guān)亞楠，等.離子液體乙醇溶液吸收SO2的物性研究［J］.河北科技大學(xué)學(xué)報，2011，32（5）：507－512.ZHANG Miaomiao，REN Ailing，GUAN Ya＇nan，et al.Study on physical properties of ionic liquid and ethanol solution for absorbing SO2［J］.Journal of Hebei University of Science and Technology，2011，32（5）：507－512.