摘 要：低共熔溶劑（DES）是一種類離子液體的綠色溶劑，具有蒸汽壓低、熱穩(wěn)定性好、制備簡(jiǎn)單等優(yōu)點(diǎn)，但高黏度特性制約了其工業(yè)應(yīng)用，同時(shí)DES 具有強(qiáng)吸濕性，因此探究體系中水的存在對(duì)其黏度的影響至關(guān)重要. 本文利用分子動(dòng)力學(xué)模擬的方法，在GAFF 力場(chǎng)下分別模擬了氯化膽堿+乙二醇/甘油/尿素等三種DES 水合物在293. 15～353. 15 K溫度范圍內(nèi)的密度和黏度，通過(guò)分析氫鍵數(shù)目和徑向分布函數(shù)，研究了含水量對(duì)DES 微觀結(jié)構(gòu)及宏觀性質(zhì)的影響. 結(jié)果表明，水的存在破壞了DES 內(nèi)膽堿基陽(yáng)離子（CHO⁺）與陰離子（Cl^?）、陰離子（Cl^?）與氫鍵供體（HBD）、氫鍵供體間形成的氫鍵；同時(shí)，水與CHO⁺、Cl^?、HBD 形成新的弱氫鍵從而降低原有DES 體系的黏度. 當(dāng)上述三種體系內(nèi)的含水質(zhì)量分?jǐn)?shù)分別達(dá)到8 wt% 、10 wt% 、10 wt% 時(shí)，仍能維持以DES 為中心的氫鍵體系，且對(duì)應(yīng)黏度分別降低了58. 3%、73. 78%、89. 41%. 因此，可以通過(guò)添加少量水來(lái)調(diào)節(jié)DES 的黏度，進(jìn)而擴(kuò)大其工業(yè)應(yīng)用范圍.

關(guān)鍵詞： 低共熔溶劑；密度； 黏度；氫鍵；分子動(dòng)力學(xué)模擬

中圖分類號(hào)： O55 文獻(xiàn)標(biāo)志碼： A DOI： 10. 19907/j. 0490-6756. 240246

1引言

2003年，Abbott 等［1］首次以1∶2的摩爾比將氯化膽堿與尿素混合加熱攪拌，冷卻后得到一種常溫下呈液態(tài)的低熔點(diǎn)溶液，稱為低共熔溶劑（DES）. DES 是由氫鍵供體（HBD）和氫鍵受體（HBA）組成的混合物，其中大量的氫鍵破壞了混合物的晶格結(jié)構(gòu)，使晶格能下降，熔點(diǎn)降低. 近年來(lái)，隨著現(xiàn)代社會(huì)對(duì)綠色能源需求的日益增長(zhǎng)，DES 作為一種類離子液體的環(huán)保溶劑，因具備合成簡(jiǎn)單、無(wú)毒可降解、導(dǎo)電性優(yōu)良和熱穩(wěn)定性好等優(yōu)勢(shì)在能源動(dòng)力領(lǐng)域得到了廣泛關(guān)注.

在吸收式制冷系統(tǒng)中，傳統(tǒng)的工質(zhì)對(duì)易結(jié)晶、易腐蝕且成本較高. DES 具有飽和蒸氣壓低、熱穩(wěn)定性好、不易結(jié)晶和腐蝕金屬部件等優(yōu)點(diǎn). 2017年，Abedin 等［2］首次提出將DES 作為吸收劑應(yīng)用在吸收式制冷系統(tǒng). 2020 年，Haghbakhsh 等［3］研究了DES+水作為工質(zhì)對(duì)在吸收式制冷系統(tǒng)中的應(yīng)用，結(jié)果表明，當(dāng)氯化膽堿和乙二醇/甘油/尿素組成的DES 作為吸收劑時(shí)，對(duì)應(yīng)的性能系數(shù)（COP）分別為0. 713、0. 716 和0. 705，驗(yàn)證了DES 作為吸收劑的可行性. 此外，在儲(chǔ)能電池的開發(fā)中，DES可作為太陽(yáng)能電池的電解液傳導(dǎo)電流［4］；在生物柴油的提純中，DES 可作為萃取劑萃取副產(chǎn)物甘油［5］；在太陽(yáng)能儲(chǔ)熱器中，DES 可作為傳熱介質(zhì)通過(guò)相變存儲(chǔ)能量［6］. 由此可見，DES 作為一種環(huán)境友好型溶劑可應(yīng)用于能量的儲(chǔ)存和傳遞過(guò)程，但它的高黏度特性也會(huì)嚴(yán)重影響傳熱傳質(zhì)性能，因此有效降低其黏度是促進(jìn)DES 廣泛應(yīng)用的關(guān)鍵.

由于DES 具有強(qiáng)吸濕性，水的存在會(huì)顯著影響其黏度等物理性質(zhì). 眾多學(xué)者在分子動(dòng)力學(xué)模擬和實(shí)驗(yàn)方面研究了DES 水合物的相關(guān)特性［7，8］，但涉及的溫度區(qū)間相對(duì)較窄. 為進(jìn)一步探究寬溫度區(qū)間下水對(duì)DES 體系微觀和宏觀特性影響的機(jī)制，本文利用分子動(dòng)力學(xué)模擬研究了溫度范圍為293. 15～353. 15 K，氯化膽堿與乙二醇/甘油/尿素形成的DES 水合物的密度和黏度，并分析了DES體系的含水量對(duì)體系微觀結(jié)構(gòu)以及宏觀熱物理性質(zhì)的影響，為DES 工業(yè)應(yīng)用提供支撐.

2模型與方法

2. 1模擬體系

目前，DES含水體系以氯化膽堿+乙二醇/甘油/尿素三種DES 水合物的熱物理性質(zhì)實(shí)驗(yàn)數(shù)據(jù)最為翔實(shí)［9，10］. 為了驗(yàn)證模擬方法的可靠性，本文選擇上述三種體系，模擬了不同溫度和含水量條件下的三種含水體系的密度與黏度，并分析了水含量對(duì)DES 結(jié)構(gòu)和性質(zhì)的影響. 所模擬體系的信息如表1 所示. 其中，DES-E-H₂O、DES-G-H₂O、DES-U-H₂O 分別代表乙二醇、甘油、尿素與氯化膽堿形成的DES水合物.

2. 2水模型的選取

在313. 15 K、 0. 1 MPa 條件下，水的摩爾分?jǐn)?shù)為20% 時(shí)，利用目前常用的OPC、SPC/E、TIP3P、TIP4P 四種水模型和DES-E 體系混合進(jìn)行了分子動(dòng)力學(xué)模擬，模擬結(jié)果如表2所示，并與實(shí)驗(yàn)值［12］進(jìn)行了對(duì)比. 結(jié)果表明：當(dāng)水模型為SPC/E 時(shí)，對(duì)DES水溶液體系的模擬最好. 因此本文模擬選擇了AMBER 力場(chǎng)下的SPC/E 水模型.

2. 3模擬方法

本文模擬采用了GROMACS 2016. 4［13］，選用的模擬力場(chǎng)為GAFF［14］力場(chǎng)，這種力場(chǎng)處理各種有機(jī)小分子具有優(yōu)勢(shì). 由于GROMACS 不自帶該力場(chǎng)，本文采用了AmberTools［15］產(chǎn)生了模擬所需要的力場(chǎng)參數(shù)，再利用ACPYPE［16］腳本文件將其轉(zhuǎn)換為可識(shí)別的形式. 模擬時(shí)，首先建立6 nm×6 nm×6 nm 大小的初始盒子并進(jìn)行填充；其次對(duì)體系進(jìn)行能量最小化，利用退火設(shè)定翻越能量勢(shì)壘，最后進(jìn)行平衡相和生產(chǎn)相的模擬，時(shí)間為10 ns. DES 是類離子液體，體系中存在游離的氯離子以及部分帶電荷的基團(tuán)，不同的電荷計(jì)算方法對(duì)模擬結(jié)果會(huì)有很大的影響［17］，因此本文利用Multiwfn 3. 7［18］程序提供的限制性靜電勢(shì)擬合（RESP，Restrained Electrostatic Potential）［19］的方法來(lái)獲得體系的電荷，獲得了較好的結(jié)果.

3結(jié)果與分析

3. 1密度

密度是流體的重要物理性質(zhì)之一，由圖2 所示，每種DES 水合物的密度隨含水量的增多非線性地降低，當(dāng)水摩爾分?jǐn)?shù)超過(guò)60% 后，密度變化顯著. 其中DES-E-H₂O、DES-G-H₂O 和DES-U-H₂O體系模擬密度與文獻(xiàn)實(shí)驗(yàn)值［9，12］的平均相對(duì)偏差分別為0. 15%、0.16% 和0. 27%，證明了模擬方法和力場(chǎng)的可靠性.

3. 2黏度

黏度是表征液體流動(dòng)阻力的熱物理參數(shù). 由圖3 所示，所有DES+H₂O 體系的黏度均隨含水量的增加明顯下降，當(dāng)水摩爾分?jǐn)?shù)僅為20% 時(shí)，DES-E-H₂O、DES-G-H₂O、DES-U-H₂O 的黏度在常溫下的變化就已非常顯著，不同于醇類體系，DES-U-H₂O 的黏度隨含水量的增加下降速度更快，這可能與1 個(gè)尿素分子含有4 個(gè)氨基H 原子可以同時(shí)與水作用有關(guān). 因此在工業(yè)應(yīng)用中如需對(duì)黏度進(jìn)行調(diào)控，可以在體系中加適量的水，以滿足特定的需求.

3. 3結(jié)構(gòu)分析

宏觀黏度和氫鍵的行為密切相關(guān). 為簡(jiǎn)潔起見，本節(jié)以DES-E-H₂O 為例研究體系內(nèi)的主要分子間相互作用力，并通過(guò)氫鍵數(shù)目和徑向分布函數(shù)，討論水的存在對(duì)結(jié)構(gòu)產(chǎn)生的影響. 在其他兩個(gè)體系中，當(dāng)含水量增加，組分之間的相互作用變化與DES-E-H₂O 呈現(xiàn)出相同的趨勢(shì).

3. 3. 1分子間相互作用力 為研究氫鍵的主要存在位點(diǎn)，當(dāng)水摩爾分?jǐn)?shù)為40% 時(shí)，在0. 1 MPa、313. 15 K 下DES-E-H₂O 的徑向分布函數(shù)，如圖4所示. 從圖4a 中可以看出，當(dāng)Cl 原子為參考原子時(shí)，可以與3 種H 原子形成氫鍵，分別是膽堿基陽(yáng)離子（CHO⁺）上的羥基H 原子、乙二醇（ETG）上的羥基H 原子以及水中的H 原子. 其中，Cl 與水的氫原子（HW1）得到的峰最高，意味著Cl 和HW1形成氫鍵的可能性最大，其次是CHO⁺上的羥基H原子，最后是ETG 上的羥基H 原子. 從距離上來(lái)看，Cl 原子與HW1 之間的最近距離為2. 24 ?，Cl離子與HC21、H8 之間的最近距離均為2. 30 ?，因此，認(rèn)為DES-E-H₂O 中Cl 離子與HW1 之間的氫鍵相互作用最強(qiáng). 值得注意的是，Cl-HW1 的徑向分布函數(shù)圖呈現(xiàn)了2 個(gè)峰. 第一個(gè)峰認(rèn)為是氫鍵相互作用，第二個(gè)峰對(duì)應(yīng)的距離為3.68 ?，大于Cl原子與H 原子的范德華半徑之和2. 86 ?，不在氫鍵作用范圍內(nèi)，認(rèn)為是水分子與Cl 原子之間的強(qiáng)范德華作用.

由圖4b 所示，以O(shè)原子為參考原子時(shí)，較明顯的峰出現(xiàn)在水中氧原子OW 與CHO⁺上的羥基氫原子HC21（距離為1. 76 ?）、乙二醇上的羥基氫原子H10上（距離為1. 78 ?），均小于O和H 的范德華半徑之和（2. 5 ?），認(rèn)為是氫鍵相互作用. 由峰值來(lái)看，羥基O 原子聚集吸引H 原子的能力不如Cl 離子.

3. 3. 2氫鍵數(shù)目 對(duì)DES-E-H₂O 中可能形成氫鍵的組分進(jìn)行了氫鍵數(shù)目的計(jì)算，如圖5 所示. 由于含水量不同時(shí)，分子數(shù)也不同，因此計(jì)算的是組分間的平均氫鍵數(shù)，即：N_Hbonds=N _總/粒子數(shù). 可以看出，CHO⁺-Cl、Cl-HBD、HBD-HBD 三種氫鍵均隨水的增加而減少. 其中Cl-HBD 類氫鍵的衰減幅度最大，即水的存在主要破壞了Cl-HBD之間的氫鍵，形成了新的氫鍵Cl-H₂O；其次破壞了CHO⁺-Cl 類氫鍵，水對(duì)HBD-HBD 氫鍵的影響相對(duì)較小. 隨著水的增多，CHO⁺、Cl、ETG 能夠和水形成的氫鍵數(shù)均有所增加. 其中Cl-H₂O 氫鍵增加最多，表明水更容易和Cl 離子形成氫鍵.

由圖5 所示，當(dāng)體系內(nèi)水的摩爾分?jǐn)?shù)升至30%（即質(zhì)量分?jǐn)?shù)w=8%）時(shí)，由于水的存在新形成的氫鍵和本就存在的氫鍵數(shù)目近似相等，因此我們認(rèn)為此時(shí)DES 水合物中的主要?dú)滏I網(wǎng)絡(luò)仍然由CHO⁺ 、Cl、ETG組成的，同時(shí)其黏度減少了58.30%. 若w（H₂O）高于8%，DES中原先存在的主要?dú)滏I網(wǎng)絡(luò)將被水破壞，氫鍵體系的中心由原先的DES 轉(zhuǎn)移為水，這可能導(dǎo)致其無(wú)法再滿足工業(yè)應(yīng)用的要求.

3. 3. 3徑向分布函數(shù) 為了深入理解水對(duì)DESE-H₂O 分子間相互作用力的影響，我們計(jì)算了水摩爾分?jǐn)?shù)為20%、40%、60%、80% 時(shí)各個(gè)組分之間的徑向分布函數(shù)，如圖6所示.

從圖6a、6b 可以看出，隨著水含量的增多，CHO⁺-Cl、Cl-ETG的峰值均明顯下降，即水的存在使得原DES體系中大部分分子間的相互作用力減弱. 由圖6c，當(dāng)水摩爾分?jǐn)?shù)分?jǐn)?shù)在0%～60% 時(shí)，ETG-ETG 峰值增強(qiáng)，我們可以認(rèn)為在此范圍內(nèi)，水的存在首先破壞了Cl 離子和乙二醇分子之間的氫鍵，隨后Cl 和水形成了新的氫鍵，因此乙二醇周圍出現(xiàn)乙二醇分子的概率變大，導(dǎo)致ETG-ETG峰值變強(qiáng)；當(dāng)水摩爾分?jǐn)?shù)超過(guò)60%，ETG-ETG 的峰值開始降低，我們認(rèn)為此時(shí)Cl-H₂O 之間的氫鍵達(dá)到飽和狀態(tài)，更多的水分子分散在乙二醇分子周圍，導(dǎo)致ETG-ETG峰值開始下降. 如圖6c 所示，在不同水摩爾分?jǐn)?shù)下，ETG-ETG 的峰位基本沒(méi)有移動(dòng)，即乙二醇分子間的短程相互作用幾乎不受水的影響. 除ETG-ETG 外，DES各個(gè)組分間首峰出現(xiàn)的位置隨水的增多逐漸右移，這表明水分子滲透進(jìn)了原DES 各組分之間，使得CHO⁺和Cl、Cl和ETG逐漸遠(yuǎn)離. 其他HBD與氯化膽堿形成的DES水合物體系均遵循相似的趨勢(shì).

4結(jié)論

本研究考察了氯化膽堿與乙二醇/甘油/尿素形成的三種DES 水合物在293. 15～353. 15 K 下的密度和黏度，并將模擬結(jié)果與已有文獻(xiàn)中的實(shí)驗(yàn)數(shù)據(jù)進(jìn)行比較，以驗(yàn)證MD模擬的可靠性.

（1） 選擇了GAFF 力場(chǎng)、RESP 電荷和SPC/E水模型再現(xiàn)DES 水合物的密度和黏度特性.

（2） 通過(guò)氫鍵、徑向分布函數(shù)可知，DES 水合物的主要?dú)滏I相互作用存在于CHO⁺-Cl、Cl-HBD、Cl-H₂O、HBD-H₂O 之間. 其中，Cl 和H₂O 上的氫原子間的相互作用最強(qiáng).

（3） 對(duì)DES 水合物的微觀結(jié)構(gòu)分析表明，水的存在擾亂了DES 內(nèi)部原有的強(qiáng)氫鍵，會(huì)與陰離子（Cl^?）和HBD 建立相對(duì)微弱的新氫鍵，因此黏度顯著下降. 由氫鍵數(shù)目分析可知，對(duì)上述三種DES水合物體系，當(dāng)體系內(nèi)含水質(zhì)量分?jǐn)?shù)達(dá)到8%、10%、10% 時(shí)，黏度已降低58. 3%、73. 78%、89. 41%，此時(shí)依舊可以維持以DES 為中心的氫鍵體系. 因此，在DES 中添加少量的水來(lái)調(diào)節(jié)其熱物理性質(zhì)，可以有效地解決DES 工業(yè)應(yīng)用中遇到的高黏度、傳熱和傳質(zhì)效率低等難題.