大气污染的日益严重和全球变暖的加速到来^[1]，原因之一是在生活和生产中产生大量尾气、废气和烟气，这些气体含有含碳、氮、硫^[2]等元素的化合物不断地污染环境. 目前已经有多种方法吸收这些污染物，但或多或少的存在着一些不足. 例如工业上被广泛应用于吸收二氧化硫的石灰石技术^[3]会耗费大量的水资源并排放出二氧化碳，且这种工艺不可逆；通过离子液体吸收废气^[4-5]的能力很强且没有副产物生成，但由于离子液体价格昂贵、黏度高、有毒性等因素制约了其发展^[6].

低共熔溶剂（deep eutectic solvent, 简称DESs）是一类新型离子液体（ionic liquid，简称IL）类似物^[7]，其具有离子液体的许多特征和优良性质，同时相较于离子液体，DES合成过程简单、成本低并且无毒性、可生物降解^[8]，所以DES被视为离子液体的绿色替代溶剂并在溶解^[9-10]、材料制备^[11]、催化^[12-13]以及电化学^[14-16]等众多领域都有广泛的应用，目前在污染物气体吸收领域也展现出广阔的前景. DES的最大优点是可以通过简单地改变氢键供体与受体的种类而设计出成千上万种不同的DESs，且不同DESs吸收某种气体的能力各有高低，因此，寻找适用于能高效吸收不同气体的DES是一个重点也是一个难点. 如果通过在实验室开展有限次实验制备DESs再进行吸收不同气体的实验，会耗费大量的人力物力，而使用计算机进行分子模拟实验可以方便地建立各种DESs的模型，对DES的结构、性质以及和气体的相互作用进行模拟计算^[17]. 目前国内外已经陆续有学者使用分子模拟技术来模拟DES对气体的吸收，这种方法有利于高效地选择DESs并预测实验结果，同时通过各种数据分析可以更加清晰地展示其中的作用机理. 本文对低共熔溶剂和分子模拟技术进行简要介绍，并对DES吸收有害气体的分子模拟研究进展进行综述，为人们设计不同的DES合成出更加高效的气体捕捉剂提供基本的理论指导.

1.   低共熔溶剂概述（ Overview of deep eutectic solvents）

在2003年Abbott等^[18]将氯化胆碱和尿素经过混合搅拌制得一种低共熔混合物，这种混合物的熔点比其任意组分的熔点显著降低（图1），通常表现出液态溶剂的性质，故将其命名为低共熔溶剂，也称为低共熔离子液.

图 1  双组分相图上DESs低熔点的示意图^[21]

Figure 1.  Schematic diagram of the low melting point of DESs on the two-component phase diagram^[21]

下载: 全尺寸图片 幻灯片

DES的制备方法非常简单，只需要将两种或多种组分混合加热便可以获得相应的低共熔混合物^[19] 这些不同组分之间一般通过氢键作用力进行结合，由此可将其定义为氢键供体（hydrogen-bonding donor,HBD）和氢键受体（hydrogen-bonding acceptor, HBA）^[20].

随着研究的不断深入，近年来人们用不同种类、摩尔比的氢键供体和氢键受体组合制备出多种多样的低共熔溶剂. Smith等^[22]提出低共熔溶剂可以由一个通式Cat⁺X⁻zY来描述，在这个公式中Cat⁺代表所有的铵、磷、锍阳离子，X表示碱，一般为卤化物阴离子，Y是一类酸分子，而z则表示其数量. 根据DESs组分的不同可以大致分为4类，如表1所示，类型Ⅰ是季铵盐和金属卤化物体系，但可以形成DESs的非水合金属卤化物较少；而水和金属卤化物与季铵盐组合而成的是第Ⅱ类DESs；类型Ⅲ是季铵盐和HBD形成的一类可以溶解过渡金属的DESs，通过调整HBD就可以改变DESs的物理性质；类型Ⅳ是金属卤化物与HBD的混合物.

表 1  DESs的分类^[22]

Table 1.  Classification of DESs^[22]

类型 Type	通式 General formula	实例 Terms
类型Ⅰ	Cat+X−zMClx	M = Zn,Sn,Fe, Al,Ga,In
类型Ⅱ	Cat+X−zMClx·yH2O	M = Cr,Co,Cu,Ni,Fe
类型Ⅲ	Cat+X−zRZ	Z = CONH2, COOH, OH
类型Ⅳ	MClx+ RZ = MClx−1+·RZ +MClx+1−	M = Al, Zn Z=CONH2, COOH, OH

 | Show Table

DownLoad: CSV

低共熔溶剂有着一系列优良的物理化学性质，这使其在溶解^[23]、萃取分离、材料制备、催化以及电化学等众多领域都有广泛的应用^[24-29]，如使用氯化胆碱类DESs对生活中常见的丝瓜络进行处理，可成功分离提取出较高含量的丝瓜络纤维素^[30]，进一步研究发现氯化胆碱/二水合草酸体系的DESs具有分离三大素和促进纤维素纳米化的双重作用，结合超声辅助处理，可制备含木质素的纳米纤维素^[31]. 也有学者以氯化胆碱/三乙醇胺DES为反应介质，用不同原料合成的具有介孔结构的新型钠离子电池的负极材料，且有着很好的电化学性能^[32]. 王淑波等^[33]利用乙二醇与氯化胆碱（物质的量比为2∶1）DES可以有效破坏乙醇/水混合溶液共沸点的特性，对乙醇/水溶液进行蒸馏脱水,制备出高浓度的燃料级乙醇，由于这种DES的沸点高，萃取过程中能量消耗少，相比于有毒的乙二醇做萃取剂的方法有着显著的优势. 还有研究者将尿素、金属氯化物、三聚氰胺合成的DESs进行热解，使金属氢氧化物均匀的负载到C₃N₄上得到了一种高效的光催化固氮复合材料^[34]. 由此可见，可以通过简单地改变氢键供体与受体的种类而设计出成千上万种具有不同功能作用的DESs，利用DESs的这种可设计性，可以根据实际需要，选取不同的原料，实现其特定结构与功能的构建.

2.   分子模拟技术（ Molecular simulation technology）

分子模拟技术（molecular simulation，MS）是近年来飞速发展的一种理论计算技术，是采用合适的计算机软件，模拟化学体系的微观结构，利用数值运算、分析统计等方法对系统的平衡热力学、统计动力学等性质进行理论预测，建立多维场景以模拟分子的静态、动态结构变化，实现实验手段难以考察的化学分子作用过程与作用位点^[35-36]，通过计算机对物质的微观结构和运动进行模拟计算可以得到有关物质的一系列性质参数以及宏观特性，并帮助预测实验结果以及阐述作用机理^[37]. 研究者只需对简单分子结构进行模拟并拟合其热力学性质就能对DESs进行快速筛选，从大量配体中选择满足需求的组分，同时还能预测各种气体在DESs中的溶解度，并分析其作用机理. 相较于传统实验，分子模拟实验所需成本较低、安全性高、结果准确度高、可以大量节省时间和资源，并且计算机模拟可以得到一些无法从实验室获取的信息^[38]. 因此，运用分子模拟技术进行模拟实验正在成为目前的研究热点之一.

2.1   分子模拟方法

根据分子模拟方法所运用的理论基础可将其分为量子力学方法和经典力学方法两大类^[39]，在经典力学方法中又可以划分出多种模拟方法，如蒙特卡罗方法（Monte Carlo method， MC）、分子动力学模拟（molecular dynamics simulation, MD）、分子力学方法（molecular mechanics, MM）等，这些分子模拟方法基于不同的原理通过计算机对原子或分子模型进行计算，进而得出现代实验方法难以测算的数据参数.

2.1.1   量子力学方法

量子力学计算方法是基于量子力学（quantum mechanics）^[40]，几乎可以计算出和分子有关的一切性质且与实际实验结果非常相近. 量子力学方法使用最多的计算方法是从头算方法（abinitio method）^[41]，这种方法通过一定的原理将电子的波函数转化为其他函数形式进行计算以研究分子中电子的行为. 从头算方法计算结果极为精确，但其可以计算的体系非常小并且计算速度极慢. 随后人们通过研究不同的波函数转化方式开发出半经验方法（semi-empirical method），这种方法相比于从头算方法提升了计算速度且扩大了计算体系. 量子力学方法一般被用于计算精确的DESs构型，保证模拟结果的准确性，同时对DESs与气体的相互作用形式进行定性、定量计算. 量子力学方法尽管计算结果精确，但其对电子的计算过程困难,故不太适用于较大体系，人们逐渐开始研究非量子计算方法并取得一系列成果.

2.1.2   蒙特卡罗方法

蒙特卡罗方法^[42]是最早用于研究大分子体系的非量子计算方法. 蒙特卡罗方法通过建立体系中分子或原子的随机运动的模型，选取多个随机数，应用统计学中的方法计算概率，从而得到体系的一系列计算结果，随机数选取越多则结果误差越小. 但这种方法只能计算出结果的平均值，并不能对体系的动态过程进行研究，故逐渐被后续发展起来的分子动力学模拟所取代.

2.1.3   分子力学方法

分子力学方法又被叫做力场方法（force field method），是目前应用较多的研究手段，主要是用来计算分子的能量并找出其平衡构型^[43]. 分子力学的基本思想是由Andrews在1930年提出^[44]：每个分子的化学键都有其特殊的键长、键角值，分子会减小内部的非键作用（van der Waals力）并使键长、键角值接近原本的数值，从而变为平衡构型. 根据玻恩-奥本海默近似（Born-Oppenheimer approximation，简称BO近似，又称绝热近似）原理，在分子力学中不会去计算电子相互作用，所以计算量远远小于从头算量子力学方法，且在某些情况下分子力学的计算结果与量子力学计算结果相差不大. 因此研究者一般将分子力学方法用于DESs团簇构型的初筛选，随后再进行量子力学计算可以大大缩短模拟时间成本. 分子力学方法很适合研究计算量庞大的大分子体系，如药物、生物大分子等.

2.1.4   分子动力学模拟

分子动力学模拟是目前使用最广泛的研究复杂微观系统的技术^[45]. 在分子动力学模拟中，一个体系内每个原子核都会受到其他所有原子核以及电子形成的势场作用而做牛顿运动，计算机求解牛顿运动方程就可以得到粒子的运动规律，从而得到体系的种种性质^[46]. 分子动力学模拟的动态特性可以快速分析出DESs配体之间的相互作用位点，并提供大量可能的作用构型并结合分子力学模拟以及量子力学模拟进行筛选. 同时分子动力学模拟考虑了更多影响体系变化的因素如温度、压强等，从而使DESs吸收气体模拟结果更加符合实际. 这种模拟方法的结果虽然不如从头算动力学方法计算精确，但其优点就在于可以计算的体系非常大并且计算速度较快. 目前已经有多款分子动力学模拟计算软件被开发并应用于物理、化学、医学等多个领域.

这些模拟方法各有所长，所应用的领域也各不相同，量子力学方法与分子动力学方法是研究分子间相互作用以及其反应性最常用的方法. 量子力学方法研究的是真空条件下分子之间由电子变化而导致的一系列性质和作用，分子动力学模拟则是考虑分子在给定的力场条件下进行牛顿运动，对其运动规律求解得到体系的动力学性质. 蒙特卡洛模拟和分子力学方法由于其计算速度快，相对精度较低，一般只会用于模拟前的大量构型搜索和初步结构优化，减轻后续计算负担.

2.2   分子模拟常用软件介绍

随着分子模拟技术的广泛应用，各大公司不断开发更新不同的模拟软件以满足人们的不同研究需求. 以下将介绍几款目前使用人数较多的分子模拟软件.

2.2.1   Gaussian

Gaussian^[47]软件是目前较为流行的一款计算功能强大、应用范围广的量子化学软件. 这款软件基于量子力学理论，通过一些简单的命令就可以计算并预测体系几乎所有的性质，如过渡态能量和结构、分子轨道、原子电荷和电势、红外光谱等. Gaussian软件可以研究气、液、固相多种物质形态，并可以很快计算出体系的最稳定结构使结果更加精确.

2.2.2   ORCA

ORCA^[48]是一款完全免费的量子化学程序，无需对其代码进行编译可以直接从其官网下载安装，随着ORCA不断对其程序进行更新发展，其功能更加完善且在计算速度、结果精确度上都有着很大的提高，越来越受到研究人员的欢迎. 这款软件最大的优势在于其充分利用了密度拟合技术结合独家的COSX方法，并在计算时默认开启RI加速，使得其在处理较大体系的计算速度以及精度上都有着惊人的提升. 但是，ORCA的短板在于其没有开发出相应的图形界面且输入文件相对较难书写，计算精度相较于Gaussian还有一定不足.

2.2.3   GROMACS

GROMACS^[49]是一个可以进行分子动力学模拟并对产生的轨迹数据进行分析的程序包. 除了分子动力学方法，GROMACS还可以应用随机动力学和路径积分方法来模拟体系中的分子. 由于程序对键合作用与非键合作用的强大计算能力，使其可以对晶体、液晶、聚合物、生物大分子溶液等多种物质进行模拟研究，并且在GROMACS中不需要再编写代码来分析轨迹数据，程序会直接将其整合输出为图表的形式.

2.2.4   Materials Studio

Materials Studio^[50]（简称MS）是一款由美国Accelrys公司开发研制的专门应用于材料科学领域的模拟计算软件，该软件可以很方便地建立模型并对其进行模拟计算得到系列相关性质. 软件涵盖了分子力学、量子力学的密度泛函理论（DFT）^[51]、介观模型、分子动力学和统计相关等多种模拟方法，根据计算所使用的不同模拟方法软件细分出多个不同的模块以方便使用，可以帮助解决目前在催化剂、聚合物、固体与表面、晶体与衍射、结晶学和材料特性等多方面的问题.

以上各种软件各有特色，且都有其适用领域，随着分子模拟技术的不断发展进步，这些软件逐渐成为研究各种材料、溶剂和生物医药等领域难以或缺的工具.

3.   DESs吸收有害气体的分子模拟研究（Molecular simulation of harmful gas absorption by DESs）

目前已经有很多研究报道了使用低共熔溶剂吸收各类有害气体如：酸性气体SO₂^[52-56]和CO₂^[57-62]、氮氧化物^[63]、NH₃^[64]等. 由于所研究的气体大多对人身体有害，通过计算机分子模拟实验可以避免直接接触，并且相对于实验室实验更加方便，还能加深人们对其中作用机理的理解.

3.1   DESs对硫化物的吸收

气态硫化物是大气污染和损害人身体的主要物质之一，例如在工业上冶炼金属和化石燃料的燃烧等所产生大量的SO₂，在一些化工过程还会产生副产品H₂S等^[65]，如果不加以控制，这些气体会慢慢毒害人们的身体，甚至形成酸雨造成极大的破坏^[66]. 传统脱硫技术主要是通过碱性脱硫剂与烟气中的硫化物反应而将其固定，但传统方法效率低、成本高还容易产生二次污染^[67]，发展新技术进行脱硫处理并改善环境问题是社会发展的必然趋势. 使用DESs进行硫化物的捕集很好的弥补了传统技术的不足，是一种极具发展前途的方法. 已报导的可以吸收二氧化硫的DESs已经有很多种，不同种类甚至不同配比的DESs对SO₂的吸收能力都不同，通过分子模拟实验找到成本低、稳定并且吸收效果好的DESs对此类研究的后续发展有着重要意义. 其中以氯化胆碱类DESs居多，如Li^[68]应用基于量子力学的多种方法对DESs-SO₂的相互作用进行研究，希望能帮助人们设计出可以高效吸收SO₂的低共熔溶剂. 通过对体系的电荷分析表明，氯化胆碱-甘油在和SO₂分子相互作用时会发生电荷的转移，从而使DESs分子带正电，SO₂分子带负电，这是DESs吸收SO₂的重要因素之一.

Korotkevich等^[69]使用LAMMPS软件对氯化胆碱-甘油DES以及和SO₂的混合物（图2）进行了从头算分子动力学（AIMD）研究，通过径向分布函数（RDF）分析体系中不同组分间的相互作用力，发现在二氧化硫加入后，由于二氧化硫分子与氯离子形成了强络合物，影响了DESs中原有的氢键网络，原本DESs中氯离子与羟基间的相互作用都变弱，这也解释了加入二氧化硫后体系流动性增强的原因.

图 2  分子模拟盒示意图以及组分的球棍模型^[69]

Figure 2.  Schematic diagram of molecular simulation box and ball and stick model of components^[69]

下载: 全尺寸图片 幻灯片

郭连晓^[70]研究了6种氯化胆碱型DESs（配体分别是乙二醇、丙三醇、硫脲、尿素、丙二酸和乳酸）对SO₂的吸收能力和机理，绝大部分DESs吸收二氧化硫的过程都是由二氧化硫分子中的硫与氯离子之间的强相互作用主导的，而二氧化硫分子中的氧与配体和阳离子上的氢形成的相互作用较弱. 在有少量水存在时，DESs吸收二氧化硫的位点并没有改变，并且吸收二氧化硫的结合能要大于无水体系，即水分子提高了对二氧化硫的吸收能力. 水分子的存在还会降低DESs的黏度并提升二氧化硫的吸收速率.

García^[71]和其研究团队建立了3种氯化胆碱基DESs与SO₂ 、CO₂和模拟的烟气成分的界面接触模型，并进行分子动力学模拟，以进一步研究DESs在吸收酸性气体时的界面作用机理. 结果表明SO₂和CO₂的吸收过程相似，都是先在DESs表面大量聚集，经过一段时间后DESs的吸附层进行重组，气体分子之间的相互作用减弱而与DESs阳离子的相互作用增强，气体分子便逐渐开始进入DESs内部，由于这几种DESs和SO₂之间更强的相互作用使SO₂的吸收过程更加迅速.

刘妍^[72]对氯化胆碱和配体（乙二醇、丙三醇、丙二酸、硫脲、尿素）的电荷分布以及结合能进行分析发现，氯化胆碱和配体的物质的量比为1∶2时所形成的DES团簇最稳定. 将5种DESs吸收SO₂的过程分别进行分子模拟，发现DESs与二氧化硫的结合方式可以分为两种情况：第一种适用于配体为醇类和酸类的DESs，二氧化硫会进入DESs内部并破坏内部氢键与氯离子形成S—Cl⁻相互作用. 另一种情况则是DESs的配体是硫脲和尿素时，二氧化硫会附着在DESs表面，不破坏DESs的构型与HBD中电负性较大的原子结合形成O—S、S—S相互作用. 在这两种不同的结合机理下合理进行预测，若DESs中HBD与氯离子的作用力越弱或者二氧化硫与HBD的作用力越强，则DESs吸收二氧化硫的能力越强.

在实际应用中，不仅要求DESs对二氧化硫的吸收能力要好，对SO₂的解吸能力一样重要，很多DESs由于和SO₂之间强烈的化学相互作用导致后续难以对SO₂进行解吸. Hou等^[73]以1-丁基-3-甲基咪唑鎓盐（BmimCl）、乙烯脲（EU）和4-甲基咪唑（4-CH₃-Im）为原料制备出一种三元DES，这种DES可以在烟气中高效的捕获低浓度的SO₂，并具有出色的解吸能力如图3. 使用Gaussian 09软件进行量子化学计算，发现在吸收过程中4-甲基咪唑中的N原子会和SO₂中的S原子逐渐靠近，作用力增强. DES中的氯离子会转移N原子上的H⁺使SO₂形成S—OH结构，而乙烯脲可能正是通过影响这一过程来调控DES对SO₂的吸收和解吸.

图 3  三元DESs捕获SO₂示意图^[73]

Figure 3.  Ternary DESs capture SO₂ schematic diagram^[73]

下载: 全尺寸图片 幻灯片

Atilhan等^[74]应用密度泛函理论和分子中原子的量子理论（QTAIM）分析了氯化胆碱、1-丁基-3-甲基咪唑氯化物（Bmim）和1-乙基-3-甲基咪唑氯化物（EMIM）作为HBA，乙酰胺、乙酰丙酸、柠檬酸、果糖、乙二醇、甘油和乳酸作为HBD的共11种不同DESs与SO₂的相互作用. 首先使用Avogadro软件对所有DESs-SO₂构型进行结构和能量优化，并选择其中能量最低的构型进行研究. 对DESs与SO₂的相互作用过程模拟发现DESs在吸收SO₂过程中自身结构并不会被破坏，理论上可以重复使用. QTAIM分析可以计算出体系中各种相互作用的类型与大小，Reduced Density Gradient（RDG）分析进一步证明了DESs和SO₂间主要为范德华型相互作用.

Liu等^[75]使用蒙特卡洛（Monte Carlo method）模拟和量子力学方法研究了氯化胆碱-尿素吸收H₂S气体的机理，这是目前为数不多的关于DESs吸收H₂S气体的分子模拟研究. 体系中的相互作用被RDF分析（图4）清晰的展现出来，吸收过程主要是通过氯离子和H₂S中氢原子之间较强的相互作用主导，因此可以通过调整DESs中组分的摩尔比来改变其吸收H₂S的能力.

图 4  DES-H₂S内部相互作用的径向分布函数^[75]

Figure 4.  Radial distribution function of DESs-H₂S internal interaction^[75]

下载: 全尺寸图片 幻灯片

综上所述，通过分子模拟研究DESs对硫化气体的吸收主要体现在对SO₂气体的吸收，而各种不同的DESs吸收SO₂的机理有所差异，总结来说可以分为以下两种：一是认为二氧化硫中的硫原子一般会与DESs中的氯离子形成强相互作用，而氧原子与DESs中某些原子形成弱相互作用. 二是认为DESs在吸收SO₂过程中可能会发生电荷转移导致两者之间相互吸引. 而DESs吸收H₂S气体的分子模拟研究不多，主要观点也是认为H₂S是通过氢原子和DESs氯原子之间的氢键作用被吸收.

3.2   DESs对CO₂的吸收

随着科技的不断发展进步，各种新型绿色可再生能源也正在不断被开发利用，但化石燃料在未来几年仍然会在能源生产中占很大比重，这不可避免地导致大量CO₂排放而加剧全球变暖和海洋酸化，因而CO₂气体的捕获自然也成为近年来的研究热点. 目前工业上主要使用烷醇胺^[76]吸收二氧化碳气体，但这种方法成本高，且吸收剂容易蒸发或降解，开发合适的二氧化碳捕获技术同样显得尤为关键.

García等^[77]首次使用密度泛函理论（DFT）在分子层面上研究了二氧化碳吸收的短程相互作用. 所选用的DESs为氯化胆碱-甘油和氯化胆碱-丙二酸，并且考虑了氯化胆碱与氢键供体在3种不同物质的量比情况下捕获二氧化碳的能力. 对DESs几何构型与能量的计算表明DESs中的主要相互作用在氯离子和羧基或者羟基中的氢原子之间产生. 而不同的DESs构型在吸收CO₂时所产生的分子间相互作用不同，CO₂分子主要和氯离子、HBD分子作用紧密，同时还会产生部分羟基与CO₂的相互作用. 由于氯化胆碱的吸湿性，这类DESs在使用过程中不可避免地会吸收水分，通过研究有水体系下DESs吸收CO₂的作用位点与相互作用能表明，H₂O分子与DESs的亲和力更强，会抢占CO₂的最佳作用位点，但同时会通过电荷的转移增强CO₂与DESs的结合强度，相互作用能的计算结果显示水的存在并不会影响DESs对CO₂的捕获能力. 为了较为真实地模拟烟气，后续计算模型中加入了一定比例的O₂、CO₂、N₂、H₂O分子，自扩散系数的计算结果表明体系中的其他分子不会对CO₂吸收造成影响，这也验证了他们之前实验中水与DESs捕获二氧化碳能力关系的结论^[71].

Ullah等^[78]通过实验与分子模拟相结合的方法，研究了氯化胆碱-乙酰丙酸（LEV）DES的密度和黏度，得出了分子模拟的结果与实验数据较为吻合的结论. 使用密度泛函理论研究了体系中的短程属性，对纯DESs系统和DES-CO₂系统进行了径向分布函数和空间分布函数计算，发现DES中的氢键相互作用主要是乙酰丙酸与氯化胆碱中的氯原子产生，在吸收CO₂时这种相互作用不会改变，同时CO₂分子与乙酰丙酸分子间形成比与离子更强的相互作用.

Salehi等^[79]使用蒙特卡罗（MC）模拟方法研究了CO₂、CO、H₂S等气体在氯化胆碱-尿素、氯化胆碱-乙二醇中的溶解度，溶解度通过亨利常数来表征，研究所采用的力场为GAFF （Generalized AMBER Force Field） . 在分析数据时发现DESs的密度以及径向分布函数（RDF）的计算结果与其他学者的结果有些许出入，这主要是由于选择了不同的模拟方法和不同力场参数所导致. 作者通过降低离子的电荷标度因子（一种力场参数），发现离子间的相互作用也随之降低，体系的密度随之降低，气体在DESs中的溶解度会升高，所以为体系选择合适的力场参数可以使模拟更加接近真实水平. 计算出的亨利常数表明H₂S、CO₂在这两种DESs中的溶解性较好，CO则相对难溶.

Altamash等^[80]使用MDynaMix v.5.2软件对氯化胆碱-苯乙酸的结构、性质进行动力学模拟并研究了其界面行为来推断对CO₂捕集的可行性. 采用密度泛函理论计算体系的相互作用能表明，氯化胆碱-苯乙酸相比于之前报道的氯化胆碱-乙酰丙酸对CO₂有更好的捕获性能. 在吸收CO₂分子过程DES的体积膨胀非常小，即DES可以在不削弱其分子间相互作用的情况下重新排列自身结构. 这个现象在Wang等^[81]的实验中也被发现并佐证，在吸收CO₂前后，DES的径向分布函数峰的位置与大小几乎相差不大，且DES体系中组分的相互作用能前后变化不超过2%. 对DES吸收烟气中二氧化碳的界面行为模拟如图5表明^[80]，在吸收CO₂时界面处会出现3个不同的吸附区域：最外层的N₂吸附层、靠近cl⁻的水分子吸附层和靠近苯乙酸的CO₂分子层. 水分子层的存在虽然不会影响CO₂的吸收但会降低被吸附的CO₂分子向溶剂内部扩散的速度.

图 5  DESs吸收CO₂时的快照^[80]

Figure 5.  Snapshot of CO₂ absorption by DESs^[80]

下载: 全尺寸图片 幻灯片

Shen等^[82]利用DESs对不同气体的溶解度不同这一特点对甲烷/二氧化碳气体进行分离，希望找到较为有效的气体分离系统. 所使用的DESs体系为氯化胆碱/乙二醇和氯化胆碱/乙酰丙酸，物质的量比均为1∶2，随后建立了带有孔隙的石墨与二氧化钛模型，让混合气体通过填充有DESs的孔隙，主要是通过比较不同条件下气体的溶解性、扩散性和渗透性来表征DESs分离二氧化碳的能力. 然而气体在DESs膜内的渗透性与吸附和扩散过程都紧密相关，不同的模型中哪种行为占据了主导地位就会导致不同的渗透性.

以上研究表明，不同的氯化胆碱基DESs对CO₂均有较好的吸收效果，一般是通过氯化胆碱中的氯离子或HBD与CO₂产生相互作用. CO₂与DESs所产生的氢键相互作用以及范德华相互作用是吸收的主导因素，因此通过选择与CO2具有强相互作用的配体制备DESs可以达到更好的吸收效果. 而在实际应用中还要考虑温度、压强、湿度等因素影响以及吸收后的解吸能力等，这些是分子模拟研究难以完成的.

3.3   DESs对含氮化合物的吸收

有毒有害的含氮化合物一般包括氮氧化物（NO）和NH₃等，其中氮氧化物会破坏臭氧层、形成酸雨造成大气污染，还会影响身体健康导致一系列疾病^[83]，但目前关于DESs对这类物质吸收的相关研究相对还较少. Hizaddin等^[84]通过COSMO-RS模拟对94种DESs脱氮能力进行了实验预测. 发现季铵盐和磷组成的DESs有较好的脱氮潜能，而这种能力的强弱与DESs中阴、阳离子的种类和HBD及其摩尔比都相关.

魏光森^[85]使用1,3-二甲基硫脲/四丁基氯化膦DESs吸收NO气体，通过核磁共振氢谱图发现在吸收NO以后DES的氢谱图出现了新的峰，分析发现是1,3-二甲基硫脲中的一个N—H键脱氢后吸收NO分子形成新的结构. 随后通过Gaussian09软件对红外吸收峰的计算也验证了氢谱图的结果，DESs中的一个N原子作为吸收位点吸收了两摩尔的NO分子形成了—N—（NO—）—N=O的结构如图6.

图 6  DESs在吸收NO后的最优构型（N1为吸收位点）^[85]

Figure 6.  Optimal configuration of DESs after absorbing NO （N1 is the absorption site） ^[85]

下载: 全尺寸图片 幻灯片

张吕鸿等^[86]通过Gaussian软件对以苯甲酸、尿素、咪唑和硫脲为氢键供体，离子液为氢键受体的4种DESs进行电荷密度计算，并将结果与NO吸收能力对比发现，DESs中核心原子电荷密度越低，则其去质子化能力越强，吸收NO的能力越强. DESs吸收NO是通过化学吸附作用，故提高DESs中可与NO反应的组分的比例，就可以提高吸收NO的能力. 在这项实验中以苯甲酸为HBD的DESs吸收NO的效率最高.

除了氮氧化物，氨气（NH₃）也是一种有害气体，同时它又是一种重要化工原料，无论是出于大气污染的治理还是原料的回收，研究对氨气的吸收都有重要意义. Deng等^[87]用NH₄SCN作为氢键受体，甘油、乙二醇、尿素、乙酰胺作为氢键供体制备的DESs吸收NH₃，发现每克NH₄SCN-甘油（物质的量比为2:3）DESs在303 K和0.10 MPa条件下可以吸收0.223 g的NH₃，通过与其他文献的数据比较发现，这个吸收能力高于所有已经报道的DESs. 通过分子动力学模拟和密度泛函理论计算发现NH₄SCN-甘油DESs与NH₃之间存在的是一种较弱的相互作用. 为了更好地实施绿色可持续发展概念，研究者们^[88]选择氯化胆碱和天然糖（木糖、核糖和果糖）等均无毒且可生物降解的原料，混合制备成真正意义上“绿色溶剂”的DESs用于吸收NH₃，通过核磁共振氢谱表明在吸收了氨气以后，DESs中羟基的H⁺峰消失了，说明这类DESs中的羟基是NH₃的吸收位点，随后进行的量子化学计算也证明了羟基与NH₃之间形成了氢键相互作用.

由于普通二元DESs在吸收氨气时，HBD的酸性强会导致DESs结构被破坏，吸收时产生固体，影响吸收能力，而HBD酸性低或使用中性HBD时氨气的溶解度很低. 因而，李煜惠^[89]利用COSMO-RS模拟研究了三元杂化DES（氯化胆碱/间苯二酚/甘油）内部氢键网络的结构与相互作用，探讨了该三元DES吸收NH₃的影响因素. 径向分布函数分析表明这种杂化DES组分之间的相互作用在吸收氨气过程中没有较大变化，故其稳定性较高. 且DESs与NH₃分子间存在多重氢键作用，这种杂化DES对氨气的溶解能力很强.

从上述研究可以发现，分子模拟作为目前常用的理论分析手段被广泛应用于有害气体治理领域，尤其是结合目前新兴的低共熔溶剂，利用分子结构模型，通过对电荷、能量、分子轨道和相互作用等进行分析，可以帮助研究者有选择性的调控低共熔溶剂组成，并根据应用需要选择合适的溶剂，这种方法已被众多学者证明其准确性(见表2). 这种理论预测方法可以节约大量溶剂筛选的时间和经济成本，机理分析结合实验验证也更加具有说服力，通过低共熔溶剂吸收有害气体的分子模拟研究是非常有潜力和意义的一种方法.

表 2  已报导的低共熔溶剂吸收有害气体的分子模拟研究

Table 2.  Reported Molecular Simulation Study of Harmful Gas Absorption by Deep Eutectic Solvents

有害气体 Harmful gases	低共熔溶剂 DESs	软件 Software	相互作用机制 Interaction mechanism	参考文献 Reference
SO2	氯化胆碱/丙三醇	LAMMPS	SO2···氯离子	[69]
SO2	氯化胆碱/乙二醇	Gaussian	S（SO2）···氯离子	[72]
SO2	氯化胆碱/丙三醇	Gaussian	S（SO2）···氯离子	[72]
SO2	氯化胆碱/丙二酸	Gaussian	S（SO2）···氯离子	[72]
SO2	氯化胆碱/尿素	Gaussian	S（SO2）···O（尿素）	[72]
SO2	氯化胆碱/硫脲	Gaussian	S（SO2）···S（硫脲）	[72]
SO2	氯化胆碱/丙三醇	GaussianGAMESS	S（SO2）与氯离子、电荷转移相互作用	[68]
CO2	氯化胆碱/丙二酸	Gaussian	CO2···氯离子、CO2··· -COOH	[77]
CO2	氯化胆碱/丙三醇	Gaussian	CO2···氯离子、CO2··· -OH	[77]
CO2	氯化胆碱/乙酰丙酸	GaussianMDynaMix	C（CO2）··· -COOH（乙酰丙酸）	[78]
CO2	氯化胆碱/苯乙酸	GaussianMDynaMix	C（CO2）···苯环（苯乙酸） 、（CO2）··· -OH（氯化胆碱）	[80]
H2S	氯化胆碱/尿素	RASPA Gaussian	H（H2S）···氯离子、H（H2S）···O（尿素）	[75]
NO	1,3-二甲基硫脲/四丁基氯化膦	Gaussian	NO···1,3-二甲基硫脲化学反应	[85]
NO	苯甲酸/氯化四丁基膦	Gaussian	NO···苯甲酸化学反应	[86]
NH3	硫氰酸铵/丙三醇	Materials studio	NH3··· -OH（丙三醇）	[87]
NH3	氯化胆碱/木糖	Gaussian	NH3···氯离子、NH3··· -OH（木糖）	[88]
CH4	氯化胆碱/乳酸	ORCA MDynaMix	H（CH4） ···O（氯化胆碱）	[90]
CH4	甜菜碱/乳酸	ORCA MDynaMix	C（CH4）··· -COOH（乳酸）	[90]
CH4	丙氨酸/乳酸	ORCA MDynaMix	C（CH4）··· -OH（乳酸）	[90]

 | Show Table

DownLoad: CSV

4.   总结与展望（Summary and Outlook）

从目前众多学者所研究的结果来看，低共熔溶剂这类溶解性强、可自由调控且制备简单的绿色溶剂的确是未来用于吸收各种有害气体的理想溶剂. 而随着计算机技术的不断发展，软件模拟的准确性越来越高，计算机模拟技术成为预测反应进程、分析反应机理的有效手段，不断有学者将其用于DESs吸收有害气体的研究.

对这些研究进行汇总比较发现，大部分气体的吸收过程是通过氯化胆碱类DESs与气体之间的分子间作用力而产生的物理相互作用，但也存在一些DESs与气体分子之间发生化学作用生成化学键. 分子模拟的结果表明在DESs内部组分之间存在较强的氢键网络，吸收有害气体时DESs组分中的某些原子会与气体分子之间产生相互作用并互相吸引以达到吸收气体的目的，这种相互作用的种类以及强弱与DESs和所吸收有害气体的种类有关，同时DESs组分的配比也会对吸收能力产生一定影响. 大部分DESs吸收前后的相互作用对比表明在吸收气体前后内部结构并没有发生变化，即这些DESs在对气体进行解吸后可以循环使用. 但目前缺少对不同体系DESs在相同条件下吸收有害气体的分子模拟对比研究，难以从理论的角度确定适合某种有害气体的最佳吸收剂.

另外在进行分子模拟时也存在一定问题，分子模拟时各种模拟参数的选择尤为重要，合适的参数会使模拟结果更加符合真实情况，例如分子动力学模拟中的分子力场和量子力学中的基组方法都是模拟结果准确与否的关键性因素. 而目前没有明确的使用关系，只能依靠研究者们的不断探索积累经验而确定更加合适的参数，这在一定程度上制约了计算机分子模拟的快速发展. 另一方面，建模以及模拟过程中的操作不同也会对结果造成影响，目前这方面的计算机模拟方法和软件较多，不同研究者使用的软件和操作等的差异会造成无法避免的误差并使实验结果之间的比较分析产生争议. 但在实际应用中，DESs的黏度、pH、极性以及所含基团等因素都会影响到气体吸收效果，同时在吸收过程中，温度、压力、时间、气体流速和水含量等也会对吸收容量及效率产生影响，这些因素很难通过分子模拟体现出来. 希望在不久的将来更多操作简单、适应性强、模拟实验数据准确的分子模拟软件可以被开发出来.