化工学报, 2022, 73(4): 1606-1614 doi: 10.11949/0438-1157.20211853

分离工程

离子液体萃取硝酸中Ce（Ⅳ）的动力学研究

李春晖,1, 何辉1,2, 何明键1, 张萌1, 高杨1, 矫彩山,1

1.哈尔滨工程大学核科学与技术学院，黑龙江 哈尔滨 150001

2.中国原子能科学研究院放射化学研究所，北京 102413

Extraction kinetics of Ce(Ⅳ) from nitric acid solutions using ionic liquid

LI Chunhui,1, HE Hui1,2, HE Mingjian1, ZHANG Meng1, GAO Yang1, JIAO Caishan,1

1.College of Nuclear Science and Technology, Harbin Engineering University, Harbin 150001, Heilongjiang, China

2.Department of Radiochemistry, China Institute of Atomic Energy, Beijing 102413, China

通讯作者: 矫彩山（1964—），男，硕士，教授，jiaocaishan@163.com

收稿日期: 2021-12-30   修回日期: 2022-01-21  

基金资助: 国家自然科学基金项目.  U1967219

Received: 2021-12-30   Revised: 2022-01-21  

作者简介 About authors

李春晖（1991—），男，博士研究生，lichunhui@hrbeu.edu.cn , E-mail：lichunhui@hrbeu.edu.cn

摘要

通过恒界面池法研究离子液体1-丁基-3-甲基咪唑双三氟甲磺酰亚胺盐（C4mimNTf2）萃取硝酸溶液中四价铈的动力学，并结合分子动力学模拟分析C4mimNTf2在萃取相界面附近的分布及结构特征。萃取动力学实验通过考察搅拌速率、相界面面积、温度以及两相不同组分浓度对Ce(Ⅳ) 萃取速率的影响，表明萃取过程为水相扩散-界面反应控制模式，且磷酸三丁酯（TBP）和C4mimNTf2间存在动力学的反协同萃取效应。分子动力学模拟分析不同组分在两相中的分布规律，结果验证了C4mimNT2在相界面处存在明显吸附现象。研究表明界面过程对离子液体萃取体系的动力学行为有重要影响，可以为离子液体萃取体系的萃取过程和机理的研究提供一定参考。

关键词： 溶剂萃取 ; 动力学 ; 离子液体 ; 分子模拟 ; 反协同效应 ; 界面

Abstract

The kinetics of extraction of tetravalent cerium from nitric acid solution by ionic liquid 1-butyl-3-methylimidazolium bis(trifluoromethanesulfonyl)amide (C4mimNTf2) in 1,2-dichloroethane was studied by the constant interface cell method. The molecular dynamics (MD) simulations were performed to reveal the effect of the interface on the extraction kinetics. In the experimental section, the extraction rate regime was analyzed according to the dependences of extraction rates on the stirring speed, interfacial area, and temperature. The effect of solution compositions on the extraction rate was also investigated. The experimental results show that a mixed diffusional-kinetic regime occurs, where the diffusion resistance is in the aqueous phase, and chemical reactions occur at the interphase. The aqueous nitric acid concentration significantly influences the extraction rate, because it affects the coordination behaviors of tetravalent ceric ions. The apparent rate constant under the experimental condition was determined as 0.094 L·cm·mol-1·min-1 by fitting the experimental data. It also indicates that the additional tri-n-butyl phosphate (TBP) in the organic phase exerts the antagonistic effect on the kinetic with C4mimNTf2. The MD simulations employed the general AMBER force field (GAFF) combined with the restrained electrostatic potential (RESP) atom charge. The scaled charge scheme with a scaled factor of 0.8 was adopted for ionic liquids to avoid overestimating the interaction energy between cations and anions. A mixed dissociation model of nitric acid was employed to consider the state of nitric acid molecules in the aqueous phase. In simulations, the number and charge density distribution along the axial, the radial distribution functions (RDF), and diffusion coefficients were calculated. The simulation results show that C4mimNTf2 in the organic bulk tends to adsorb at the interface with the polar groups toward the aqueous phase. The center-of-mass RDF between C4mim+ and NTf2- depicts that C4mimNTf2 exists in the form of ion pairs or clusters. With the addition of TBP in the organic phase, C4mimNTf2 can bond with TBP by the hydrogen bonds between the imidazolium ring and the terminal oxygen atom of the P̿    O group in TBP. Meanwhile, the simulated diffusion coefficients of C4mimNTf2 and TBP decrease with the increased organic phase TBP concentration. The hydrogen bond interactions and the variation of diffusion coefficients lead to the antagonistic effect. This study demonstrates that the interface behavior significantly affects the extraction of Ce(Ⅳ) from nitric acid solutions, which provides some reference for the extraction process and mechanism of the ionic liquid-based extraction system

Keywords： solvent extraction ; kinetics ; ionic liquid ; molecular simulation ; antagonistic effect ; interface

PDF (1577KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

李春晖, 何辉, 何明键, 张萌, 高杨, 矫彩山. 离子液体萃取硝酸中Ce（Ⅳ）的动力学研究[J]. 化工学报, 2022, 73(4): 1606-1614. doi:10.11949/0438-1157.20211853

LI Chunhui, HE Hui, HE Mingjian, ZHANG Meng, GAO Yang, JIAO Caishan. Extraction kinetics of Ce(Ⅳ) from nitric acid solutions using ionic liquid[J]. CIESC Journal, 2022, 73(4): 1606-1614. doi:10.11949/0438-1157.20211853

引言

离子液体是由阴阳离子构成，在室温下呈液态的有机盐类化合物。以离子液体为溶剂的萃取体系对核电站乏燃料水法后处理中关键核素通常表现出更高的萃取效率和选择性，但其萃取行为和机理通常更为复杂[1-4]。同时，由于离子液体的阴离子通常具有一定亲水性，因此离子液体在溶剂萃取过程中不仅作为溶剂，也可以通过阴离子交换萃取水相中以配阴离子或含氧酸根形式存在的金属离子，如Pu(Ⅳ)[4]、Ce(Ⅳ)[5-6]、Tc(Ⅶ)[7]等。因此在离子液体体系萃取水相阴离子组分的研究中，需考虑离子液体自身的萃取行为，否则实验结果的解释可能偏离实际。Rout等[8-9]提出对离子液体/有机分子溶剂体系的研究，有利于理解离子液体体系的萃取行为及萃取机理；该思路也已被应用于研究离子液体和萃取剂间的协同萃取行为[10-11]。Li等[6]将1-丁基-3甲基咪唑双三氟甲磺酰亚胺盐（C4mimNTf2）溶解于1,2-二氯乙烷中，考察了其萃取硝酸中Ce(Ⅳ)的行为和机理，结果表明C4mimNTf2更倾向于通过阴离子交换萃取单价组分HCe(NO3)6-，且磷酸三丁酯（TBP）和C4mimNTf2之间对Ce(Ⅳ)的萃取存在热力学的反协同效应。尽管对离子液体萃取体系热力学已开展了丰富的工作，但是相关动力学研究却少有报道。恒界面池法是研究溶剂萃取动力学的有效方法之一[12-13]，目前在离子液体体系的萃取动力学研究中应用较多，研究结果表明由于离子液体的高黏度，有机相内扩散及化学反应过程通常为离子液体萃取体系的控制步骤之一[14-17]。同时，离子液体萃取体系相界面的分子动力学研究表明，由于离子液体的离子性，溶质离子在相界面附近可以被离子液体阳（或阴）离子溶剂化，进而转移至有机相中[18-19]。萃取动力学以及相界面行为研究对于揭示离子液体体系的萃取机理具有重要的意义。

由于具有相近的离子半径和配位结构，Ce(Ⅳ)通常作为Pu(Ⅳ)的非放射性替代元素研究其萃取行为[20]。同时，由于离子液体C4mimNTf2具有较高的化学和辐照稳定性，以及相对较低的黏度，在离子液体萃取体系中应用较为广泛。因此，本文通过恒界面池法，以1,2-二氯乙烷为有机相稀释剂，研究C4mimNTf2对硝酸介质中Ce(Ⅳ)的萃取动力学，并考察有机相中萃取剂TBP对体系萃取动力学的影响；同时结合分子动力学模拟研究萃取体系相界面的结构及特征。

1 研究方法

1.1 实验试剂

若无额外说明，实验药品均为分析纯试剂且均未经任何预处理。其中，C4mimNTf2购买于林州科能有限公司，> 99.0%；1,2-二氯乙烷购于麦克林公司，为色谱纯试剂，> 99.7%；TBP购于麦克林公司，为分析纯试剂；硝酸钠购于天津福晨化学试剂有限公司，为分析纯试剂；水相Ce(Ⅳ)溶液由硝酸铈铵晶体（麦克林，分析纯）溶于硝酸溶液（大茂，分析纯）制成。实验用水均为实验室自制一级去离子水。

1.2 萃取动力学实验及数据处理

恒界面池装置由中国原子能研究院设计并加工[21]，恒界面池结构及参数在图1中给出。实验前，两相料液分别水浴加热至预定温度。实验中，按一定时间间隔取水相料液测定Ce(Ⅳ)浓度。水相 Ce(Ⅳ)浓度通过分光光度法测定[22]。若无说明，所有萃取动力学实验均在25℃下进行，水相Ce(Ⅳ)浓度为7.0 mmol·L-1，水相硝酸浓度为3.0 mol·L-1，有机相C4mimNTf2浓度为0.5 mol·L-1，两相体积均为15 ml，相界面面积为10.4 cm2。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   恒界面池装置

Fig.1   Diagram of constant interfacial area cell


萃取动力学数据通过微分法处理。初始萃取速率为：

r0=-dc[Ce(Ⅳ)]/dtt = 0=k(A/V)c0[Ce(Ⅳ)]ac0(C4mimNTf2)b

式中，k为表观速率常数；A为相界面面积；V为水相体积；a、b为反应级数；下角标0表示初始时刻。实验数据通过五次函数拟合，r0为t = 0时刻的一阶导数值。

1.3 分子动力学模拟

通过分子动力学模拟研究C4mimNTf2 + 1,2-二氯乙烷 / H2O + HNO3体系的相界面性质。为考虑水相硝酸形态的影响，模拟设定50%硝酸以HNO3分子形式存在，50%硝酸电离成H3O+和NO3-[23-24]。其中，水分子选择TIP3P模型[25]，H3O+力场参数及原子电荷取自文献[26]。C4mim+、NTf2-、TBP、1,2-C2H4Cl2、HNO3、NO3-的力场参数均取自GAFF（general AMBER force field）力场[27]，其原子电荷由量子化学计算获得，流程如下：通过ORCA 4.2.1程序[28]，在B3LYP [29]-D3BJ[30] / def2-TZVP[31]级别下对体系分子进行结构优化，并通过频率分析表明不存在虚频；通过Multiwfn 3.8程序[32]内置高性能代码[33]计算溶剂化条件下的限制性拟合静电势（RESP）电荷，作为各组分的原子电荷，其中离子液体阴阳离子原子电荷乘以系数0.8，以避免过度预测离子间的相互作用。

通过Packmol程序[34]构建初始体系结构，如图2所示：水相中水分子个数为8000，硝酸分子、H3O+和NO3-个数均为216；有机相中1,2-二氯乙烷分子个数为2000，C4mim+和NTf2-个数为200。分子动力学模拟通过GROMACS[35]实现，对于静电相互作用，原子间距小于1.2 nm时，显式计算静电作用；当大于1.2 nm时，通过SPME（Smooth Particle-Mesh Ewald）方法[36]计算；范德华作用计算使用截断方法，截断距离为1.2 nm；对能量和压力均做范德华长程校正；模拟步长2 fs；使用三维周期性边界条件。模拟流程如下：首先通过共轭梯度法消除初始构型中的不合理接触；随后通过等温等压系综（NPT）在298.15 K和0.1 MPa下模拟10 ns，使体系达到平衡，该过程使用V-scaled热浴[37]和Berendsen压浴[38]；最后在相同条件下通过NPT模拟80 ns产生用于统计体系性质的产生相，该过程使用V-scaled热浴和Parrinello-Rahman压浴[39]。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   萃取两相体系示意图

Fig.2   Snapshot of the two phases system


2 结果与讨论

2.1 萃取动力学模式的判定

在恒界面池法的萃取动力学研究中，一般通过综合分析搅拌转速n、比界面面积α（α = A/V）及温度T对萃取速率r0的影响判定萃取动力学模式。实验分别测定水相及有机相中搅拌转速对Ce(Ⅳ)萃取速率的影响，如图3(a)所示。实验表明，萃取速率随水相中搅拌速率的增加而缓慢增加，但由于1,2-二氯乙烷黏度与水相近，当水相搅拌速率超过约95 r·min-1时，两相界面出现扰动，萃取速率突增；而在有机相中所考察的搅拌速率范围内，萃取速率基本维持恒定。由此可初步认为萃取过程为扩散模式或者动力学模式，主要阻力均存在于水相中。在后续实验中，两相内搅拌速率均设定为80 r·min-1。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   萃取模式判定

Fig.3   Determination of extraction control mode


图3(b)表明Ce(Ⅳ)的萃取速率r0随两相界面面积A增加而增加，二者间呈正比关系且经过原点，因此该萃取过程受界面过程影响，即萃取过程为界面化学反应控制过程[40]。

为进一步说明体系的萃取动力学模式，实验测定了20~40℃范围内萃取速率随温度的变化。图3(c)给出了lnr0关于1/T的斜率。根据Arrhenius方程拟合反应活化能：

lnr0 = -(Ea/R)T -1 + C

式中，Ea为活化能；R为摩尔气体常数；T为温度；C为常数。通常而言，化学反应相较于扩散对温度变化更为敏感，化学反应活化能通常大于42 kJ·mol-1，而扩散过程一般不超过21 kJ·mol-1[41]。拟合得到的活化能为29.9 kJ·mol-1，位于二者之间。

综上所述，可以初步判断C4mimNTf2萃取硝酸介质中的Ce(Ⅳ)是由水相扩散以及界面反应共同控制的“混合模式”。

2.2 两相组成对萃取速率的影响

2.2.1 硝酸浓度影响

在硝酸溶液中，Ce4+与硝酸根通过逐级配位生成系列配合物，Ce(Ⅳ)的物种分布与硝酸浓度密切相关。同时由于阴离子交换萃取机理，C4mimNTf2仅萃取水相中的阴离子组分，因此，水相硝酸浓度对Ce(Ⅳ)的萃取有明显影响。不同硝酸浓度下Ce(Ⅳ)的萃取速率在图4中给出。硝酸浓度低于2 mol·L-1时，Ce(Ⅳ)萃取速率较低且几乎不随硝酸浓度增加而改变；随硝酸浓度继续增加，Ce(Ⅳ)萃取速率增加；当硝酸浓度高于5 mol·L-1时，Ce(Ⅳ)萃取速率几乎不变。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   水相硝酸浓度对萃取速率的影响

Fig.4   Dependence of extraction rate on initial aqueous nitric acid concentration


分配比实验表明，C4mimNTf2 + 1,2-二氯乙烷体系倾向于萃取硝酸溶液中HCe(NO3)6-组分[6]。由 Ce(Ⅳ)在硝酸溶液中的水溶液化学研究可知，硝酸浓度小于2 mol·L-1时，水相中Ce(Ⅳ)主要以阳离子型配合物存在[42]，Ce(Ⅳ)需要通过逐级配位取代转化为阴离子型可萃组分，水相中逐级配位取代反应为控制步骤，因此萃取速率较低；当硝酸浓度高于 2 mol·L-1时，水相中Ce4+和NO3-间形成的配合物主要以Ce(NO3)62-形式存在[42]，并进一步与H+结合生成可萃的HCe(NO3)6-，水相中的配位取代反应对萃取速率的影响减弱，水相内可萃组分扩散过程逐渐成为水相内的控制步骤。

为进一步说明硝酸浓度影响的本质，分别固定水相NO3-及H+浓度为3.0 mol·L-1，测定H+及NO3-浓度对Ce(Ⅳ)萃取速率的影响，结果分别在图5和图6中给出。实验表明，不同于传统萃取体系（如TBP），水相H+浓度的增加不会导致硝酸的竞争效应[43]，在研究浓度区间内，H+浓度的增加会提高萃取效率，是由于H+的盐析效应使被萃组分在溶液中的活度增加引起。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   水相氢离子浓度对萃取速率的影响（水相硝酸根浓度c(NO3-) = 3.0 mol·L-1）

Fig.5   Dependence of extraction rate on aqueous hydrogen ion concentration with 3.0 mol·L-1 aqueous nitrate ion


图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   水相硝酸根浓度对萃取速率的影响（水相氢离子浓度c(H+) = 3.0 mol·L-1）

Fig.6   Dependence of extraction rate on aqueous nitrate ion concentration with 3.0 mol·L-1 aqueous hydrogen ion


而固定水相氢离子浓度为3.0 mol·L-1，Ce(Ⅳ)萃取速率随水相NO3-浓度增加而降低。这是由于NO3-与Ce(Ⅳ)阴离子配合物组分间的竞争作用引起。

2.2.2 Ce（Ⅳ）浓度影响

水相初始Ce(Ⅳ)浓度对萃取速率的影响在图7中给出。如图所示，Ce(Ⅳ)的萃取速率随Ce(Ⅳ)浓度增加而增加，且lgc[Ce(Ⅳ)]和lgr0呈良好的线性关系，斜率为1.18，因此可认为 Ce(Ⅳ)的萃取速率与水相Ce(Ⅳ)浓度为一级关系。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   水相Ce(Ⅳ)浓度对萃取速率的影响

Fig.7   Dependence of extraction rate on aqueous Ce(Ⅳ) concentration


2.2.3 离子液体浓度影响

图8为恒定水相组成的条件下，Ce(Ⅳ)萃取速率与有机相C4mimNTf2浓度的变化关系。萃取速率随C4mimNTf2浓度增加而增加，且lgc(C4mimNTf2)和lgr0呈直线关系，斜率为0.98，因此可认为Ce(Ⅳ)的萃取速率与C4mimNTf2浓度为一级关系。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8   有机相C4mimNTf2浓度对萃取速率的影响

Fig.8   Dependence of extraction rate on organic C4mimNTf2 concentration


综上所述，在25℃，水相硝酸浓度为3.0 mol·L-1时，C4mimNTf2萃取Ce(Ⅳ)的初始萃取速率方程可写为：

r0= k(A/V)c0[Ce(Ⅳ)]c0(IL)

通过拟合实验数据，表观速率常数k为0.094 L·cm·mol-1·min-1。

2.2.4 水相NTf2-浓度的影响

由于萃取体系的阴离子交换机理，NTf2-在萃取过程中将转移至水相，因此，实验测定了水相NTf2-浓度对萃取速率的影响，如图9所示，Ce(Ⅳ)萃取速率随水相NTf2-浓度增加而降低。离子液体萃取体系中离子交换过程会引起离子液体的组分在水相中的积累，除可能引起离子液体损耗以及设备腐蚀问题外，也会导致体系萃取动力学性能降低。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9   水相NTf2-浓度对萃取速率的影响

Fig.9   Dependence of extraction rate on aqueous NTf2- concentration


2.2.5 有机相TBP的影响

在实际中，离子液体萃取体系的有机相通常为萃取剂和离子液体混合溶液，因此，萃取剂和离子液体混合物对水相组分萃取行为的研究更具应用价值，同时对理解离子液体萃取体系的机理具有重要意义。研究表明，硝酸溶液中Ce(Ⅳ)的动力学过程在TBP体系中主要为扩散控制[43]，与C4mimNTf2体系间存在一定差异。

萃取热力学实验表明，由于有机相中TBP和C4mimNTf2通过PO基团和咪唑环上质子间的氢键作用形成复合物，TBP和C4mimNTf2的有效浓度降低，二者间存在反协同效应[6]。为进一步了解TBP对C4mimNTf2体系萃取动力学的影响，实验测定了TBP和C4mimNTf2摩尔比对Ce(Ⅳ)萃取速率的影响，其中TBP摩尔分数x(TBP)定义如下：

x(TBP)= n(TBP) / [n(TBP) + n(C4mimNTf2)]

式中，n为溶液中对应组分物质的量。实验中固定有机相中TBP和C4mimNTf2总浓度为0.5 mol·L-1，改变二者摩尔比，结果在图10中给出。当有机相中TBP摩尔分数为1时，即有机相中不含C4mimNTf2，Ce(Ⅳ)的萃取速率为0.80 mmol·L-1·min-1，其显著大于相同浓度C4mimNTf2萃取Ce(Ⅳ)的速率0.30 mmol·L-1·min-1。

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10   TBP摩尔分数对萃取速率的影响（TBP和C4mimNTf2总浓度为0.5 mol·L-1）

Fig.10   Dependence of extraction rate on the molar fraction of TBP with total concentration of TBP and C4mimNTf2 at 0.5 mol·L-1


Ce(Ⅳ)萃取速率与TBP摩尔分数间的下凹曲线关系，说明TBP和C4mimNTf2间存在动力学的反协同效应。其主要原因同样是由于有机相中TBP和咪唑阳离子间通过氢键形成复合物，复合物需要解离生成自由TBP，通过PO基团上的O原子与Ce4+配位，才能进一步萃取水相中的Ce(Ⅳ)。

2.3 萃取体系界面的分子动力学模拟

C4mimNTf2萃取硝酸溶液中Ce(Ⅳ)的动力学实验研究表明，体系相界面行为是影响萃取动力学的重要因素之一。因此通过分子动力学模拟相平衡状态下的两相溶液，从微观层面理解萃取体系的相界面行为和性质。

2.3.1 C4mimNTf2 + 1，2-二氯乙烷体系

通过分子动力学模拟研究体系中各组分沿z轴方向在两相中的分布。对于水相-有机相界面，使用“90-90”标准定义界面区域[44]，即为水和1,2-二氯乙烷在相主体浓度均为90%之间的区域，相界面区域在图11(a)中以阴影部分表示。如图11(a)所示，C4mim+和NTf2-在相界面附近存在显著的吸附现象；同时，图11(b)的电荷密度分布进一步表明，C4mim+和NTf2-中的极性区域（带电荷基团）倾向于朝水相方向分布，而非极性区域则倾向于朝有机相方向分布。

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11   C4mim+和NTf2-数密度及电荷密度在z轴方向分布

Fig.11   Number and charge density profile of C4mim+ and NTf2- along z axis


如图12所示，C4mim+和NTf2-质心间径向分布函g(r)仅在0.52 nm处存在显著的尖峰，表明二者在萃取体系中倾向于以离子对或团簇的形式存在，与C4mimNTf2在水及醇等溶液中存在形态的研究结论相一致[45]。

图12

新窗口打开| 下载原图ZIP| 生成PPT

图12   C4mim+和NTf2-质心间径向分布函数

Fig.12   Center of mass radial distribution function between C4mim+ and NTf2-


综上所述，结合萃取动力学实验及分子动力学模拟，可以认为C4mimNTf2萃取硝酸介质中Ce(Ⅳ)主要包括以下步骤：（1）Ce(Ⅳ)在硝酸溶液水相中生成可萃的阴离子HCe(NO3)6-；（2）HCe(NO3)6-从水相主体扩散至相界面；（3）C4mimNTf2从有机相主体扩散至相界面；（4）HCe(NO3)6-在界面附近与C4mimNTf2通过离子交换形成萃合物C4mimHCe(NO3)6；（5）C4mimHCe(NO3)6从界面扩散至有机相主体；（6）NTf2-从界面扩散至水相主体。萃取过程为水相扩散和界面反应控制的“混合模式”，控制步骤为水相扩散过程（2），以及界面反应（4）。

2.3.2 C4mimNTf2 + TBP + 1，2-二氯乙烷体系

当有机相中加入TBP后，如图12所示，C4mim+和NTf2-间径向分布函数峰高随体系中TBP分子个数 N(TBP)增加而增加，但峰形保持不变，表明体系中TBP不会破坏C4mimNTf2的离子对或团簇结构。同时，TBP中PO基团上的氧原子与C4mim+中咪唑环上的氢原子间的径向分布函数在图13中给出。如图所示，在0.26 nm处存在强烈的尖峰且峰高随TBP浓度增加而降低，而峰形保持不变。径向分布函数验证了TBP与C4mim+间通过氢键作用在体系中形成复合物，进而使体系中二者的有效浓度降低，表现为反协同萃取效应。

图13

新窗口打开| 下载原图ZIP| 生成PPT

图13   PO中氧原子和咪唑环上氢原子间径向分布函数

Fig.13   Radial distribution function between the oxygen atom in PO and hydrogen atom in imidazolium ring


图14给出体系中C4mim+、NTf2-及TBP的扩散系数D与TBP含量的关系，可见随体系中TBP含量增加各组分扩散系数均呈现下降趋势，即组分在有机相内的扩散速率降低。因此，有机相主体中TBP和C4mimNTf2扩散系数降低是导致动力学反协同萃取效应的重要原因之一。

图14

新窗口打开| 下载原图ZIP| 生成PPT

图14   C4mim+、NTf2-和TBP扩散系数

Fig.14   Diffusion coefficient of C4mim+, NTf2- and TBP


3 结 论

结合恒界面池法和分子动力学模拟研究了C4mimNTf2萃取硝酸介质中Ce(Ⅳ)的动力学，结果表明：

（1）C4mimNTf2存在明显的界面吸附行为，在实验条件下，萃取过程为扩散-界面化学反应共同控制的动力学模式，扩散阻力主要存在于水相中，萃取过程的表观速率常数为0.094 L·cm·mol-1·min-1；

（2）萃取体系中，C4mimNTf2存在明显的界面吸附行为，Ce(Ⅳ)配阴离子在界面附近与C4mimNTf2发生阴离子交换，并进一步转移至有机相；

（3）有机相中引入TBP后，TBP与C4mimNTf2间通过氢键作用结合，同时溶液中各组分扩散系数降低，使TBP和C4mimNTf2间表现为动力学的反协同萃取效应。

参考文献

刘海望, 杨涛, 陈庆德, 等.

离子液体体系的萃取行为及其在乏燃料后处理中的应用前景

[J]. 核化学与放射化学, 2015, 37(5): 286-309.

[本文引用: 1]

Liu H W, Yang T, Chen Q D, et al.

Extraction behaviors of ionic liquid systems and application perspectives in reprocessing of spent nuclear fuel

[J]. Journal of Nuclear and Radiochemistry, 2015, 37(5): 286-309.

[本文引用: 1]

袁立永, 石伟群, 蓝建慧, 等.

离子液体在核燃料后处理中的应用

[J]. 科学通报, 2012, 57(8): 581-593.

Yuan L Y, Shi W Q, Lan J H, et al.

Application of ionic liquids in nuclear fuel reprocessing

[J]. Chinese Science Bulletin, 2012, 57(8): 581-593.

黄焜, 李晓佩, 徐怡庄, 等.

萃取分离体系分子间弱相互作用的研究进展

[J]. 化工学报, 2016, 67(1): 152-164.

Huang K, Li X P, Xu Y Z, et al.

Research progress on intermolecular weak interaction in extraction and separation system

[J]. CIESC Journal, 2016, 67(1): 152-164.

Sivaramakrishna M, Raut D R, Nayak S, et al.

Unusual selective extraction of Pu4+ by some novel diamide ligands in a room temperature ionic liquid

[J]. Separation and Purification Technology, 2017, 181: 69-75.

[本文引用: 2]

Zuo Y, Liu Y, Chen J, et al.

The separation of cerium(Ⅳ) from nitric acid solutions containing thorium(Ⅳ) and lanthanides(Ⅲ) using pure [C8mim]PF6 as extracting phase

[J]. Industrial & Engineering Chemistry Research, 2008, 47(7): 2349-2355.

[本文引用: 1]

Li C H, He H, He M J, et al.

An experimental study on the extraction mechanisms of Ce(Ⅳ) from HNO3 solutions using C4mimNTf2 as extractant

[J]. Journal of Radioanalytical and Nuclear Chemistry, 2021: 1-9.

[本文引用: 4]

Chotkowski M, Połomski D.

Extraction of pertechnetates from HNO3 solutions into ionic liquids

[J]. Journal of Radioanalytical and Nuclear Chemistry, 2017, 314(1): 87-92.

[本文引用: 1]

Rout A, Venkatesan K A, Srinivasan T G, et al.

Tuning the extractive properties of purex solvent using room temperature ionic liquid

[J]. Separation Science and Technology, 2013, 48(17): 2576-2581.

[本文引用: 1]

Rout A, Venkatesan K A, Srinivasan T G, et al.

Ionic liquid extractants in molecular diluents: extraction behavior of europium (Ⅲ) in quarternary ammonium-based ionic liquids

[J]. Separation and Purification Technology, 2012, 95: 26-31.

[本文引用: 1]

Turanov A N, Karandashev V K, Khvostikov V A.

Synergistic extraction of lanthanides (Ⅲ) with mixtures of TODGA and hydrophobic ionic liquid into molecular diluent

[J]. Solvent Extraction and Ion Exchange, 2017, 35(7): 461-479.

[本文引用: 1]

Turanov A N, Karandashev V K, Boltoeva M, et al.

Synergistic extraction of uranium(Ⅵ) with TODGA and hydrophobic ionic liquid mixtures into molecular diluent

[J]. Separation and Purification Technology, 2016, 164: 97-106.

[本文引用: 1]

李洲.

HS-LIX65N从氯化物介质中萃铜的萃取动力学研究

[J]. 化工学报, 1988, 39(3): 299-308.

[本文引用: 1]

Li Z.

Studies on the Kinetics of extracting copper with HS-LIX65N from chloride mediums

[J]. Journal of Chemical Industry and Engineering(China), 1988, 39(3): 299-308.

[本文引用: 1]

蒋盼, 张盼良, 曾乐林, 等.

反应萃取分离2-苯基丁酸对映体的动力学

[J]. 化工学报, 2017, 68(1): 163-169.

[本文引用: 1]

Jiang P, Zhang P L, Zeng L L, et al.

Kinetics of reaction extraction separation of 2-phenylbutyric acid enantiomers

[J]. CIESC Journal, 2017, 68(1): 163-169.

[本文引用: 1]

Wang Y L, Wang Y B, Zhou H Y, et al.

Extraction kinetics of mixed rare earth elements with bifunctional ionic liquid using a constant interfacial area cell

[J]. RSC Advances, 2017, 7(63): 39556-39563.

[本文引用: 1]

Yang H L, Chen J, Zhang D L, et al.

Kinetics of cerium(Ⅳ) and fluoride extraction from sulfuric solutions using bifunctional ionic liquid extractant (Bif-ILE) [A336][P204]

[J]. Transactions of Nonferrous Metals Society of China, 2014, 24(6): 1937-1945.

Yang H L, Chen J, Wang W, et al.

Extraction kinetics of lanthanum in chloride medium by bifunctional ionic liquid [A336][CA-12] using a constant interfacial cell with laminar flow

[J]. Chinese Journal of Chemical Engineering, 2014, 22(10): 1174-1177.

Sypula M, Ouadi A, Gaillard C, et al.

Kinetics of metal extraction in ionic liquids: Eu3+/HNO3// TODGA/[C1C4im][Tf2N] as a case study

[J]. RSC Advances, 2013, 3(27): 10736-10744.

[本文引用: 1]

Sieffert N, Wipff G.

The

[BMI][Tf 2N] ionic liquid/water binary system: a molecular dynamics study of phase separation and of the liquid–liquid interface[J]. The Journal of Physical Chemistry B, 2006, 110(26): 13076-13085.

[本文引用: 1]

Chevrot G, Schurhammer R, Wipff G.

Molecular dynamics simulations of the aqueous interface with the

[BMI][PF 6] ionic liquid: comparison of different solvent models[J]. Physical Chemistry Chemical Physics, 2006, 8(36): 4166-4174.

[本文引用: 1]

邰德荣, 仝继红, 李守成.

30%TBP(煤油)/UO2(NO3)2-HNO3体系中铈的某些化学性质

[J]. 核科学与工程, 1996, 16(2): 154-159.

[本文引用: 1]

Tai D R, Tong J H, Li S C.

Studies on some chemical properties of cerium in the system of 30% TBP(OK)/UO2(NO3)2-HNO3

[J]. Chinese Journal of Nuclear Science and Engineering, 1996, 16(2): 154-159.

[本文引用: 1]

王辉, 贾永芬, 魏艳, 等.

TBP-煤油从硝酸介质中萃取Tc(Ⅶ)的动力学研究

[J]. 核化学与放射化学, 2010, 32(4): 193-198.

[本文引用: 1]

Wang H, Jia Y F, Wei Y, et al.

Extraction kinetics of technetium(Ⅶ) between TBP-kerosine and nitric acid medium

[J]. Journal of Nuclear and Radiochemistry, 2010, 32(4): 193-198.

[本文引用: 1]

吴瑞林, 朱利亚.

某些稀土元素共存下银合金中铈的光度测定

[J]. 冶金分析, 1990, 10(5): 17-19.

[本文引用: 1]

Wu R L, Zhu L Y.

Photometric determination of cerium coexistent with some other rare earth elements in silver alloys

[J]. Metallurgical Analysis, 1990, 10(5): 17-19.

[本文引用: 1]

Singh M B, Patil S R, Lohi A A, et al.

Insight into nitric acid extraction and aggregation of N,N,N′,N′-tetraoctyl diglycolamide (TODGA) in organic solutions by molecular dynamics simulation

[J]. Separation Science and Technology, 2018, 53(9): 1361-1371.

[本文引用: 1]

Gaillard C, Mazan V, Georg S, et al.

Acid extraction to a hydrophobic ionic liquid: the role of added tributylphosphate investigated by experiments and simulations

[J]. Physical Chemistry Chemical Physics, 2012, 14(15): 5187-5199.

[本文引用: 1]

Jorgensen W L, Chandrasekhar J, Madura J D, et al.

Comparison of simple potential functions for simulating liquid water

[J]. The Journal of Chemical Physics, 1983, 79(2): 926-935.

[本文引用: 1]

Noroozi J, Smith W.

Force-field-based computational study of the thermodynamics of a large set of aqueous alkanolamine solvents for post-combustion CO2 capture

[J]. Journal of Chemical Information and Modeling, 2021, 61(9): 4497-4513.

[本文引用: 1]

Wang J M, Wolf R M, Caldwell J W, et al.

Development and testing of a general amber force field

[J]. Journal of Computational Chemistry, 2004, 25(9): 1157-1174.

[本文引用: 1]

Neese F.

The ORCA program system

[J]. WIREs Computational Molecular Science, 2012, 2(1): 73-78.

[本文引用: 1]

Stephens P J, Devlin F J, Chabalowski C F, et al.

Ab initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields

[J]. The Journal of Physical Chemistry, 1994, 98(45): 11623-11627.

Grimme S, Ehrlich S, Goerigk L.

Effect of the damping function in dispersion corrected density functional theory

[J]. Journal of Computational Chemistry, 2011, 32(7): 1456-1465.

[本文引用: 1]

Schäfer A, Huber C, Ahlrichs R.

Fully optimized contracted Gaussian basis sets of triple zeta valence quality for atoms Li to Kr

[J]. The Journal of Chemical Physics, 1994, 100(8): 5829-5835.

[本文引用: 1]

Lu T, Chen F W.

Multiwfn: a multifunctional wavefunction analyzer

[J]. Journal of Computational Chemistry, 2012, 33(5): 580-592.

[本文引用: 1]

Zhang J, Lu T.

Efficient evaluation of electrostatic potential with computerized optimized code

[J]. Physical Chemistry Chemical Physics, 2021, 23(36): 20323-20328.

[本文引用: 1]

Martínez L, Andrade R, Birgin E G, et al.

PACKMOL: a package for building initial configurations for molecular dynamics simulations

[J]. Journal of Computational Chemistry, 2009, 30(13): 2157-2164.

[本文引用: 1]

Van Der Spoel D, Lindahl E, Hess B, et al.

GROMACS: fast, flexible, and free

[J]. Journal of Computational Chemistry, 2005, 26(16): 1701-1718.

[本文引用: 1]

Essmann U, Perera L, Berkowitz M L, et al.

A smooth particle mesh Ewald method

[J]. The Journal of Chemical Physics, 1995, 103(19): 8577-8593.

[本文引用: 1]

Bussi G, Donadio D, Parrinello M.

Canonical sampling through velocity rescaling

[J]. The Journal of Chemical Physics, 2007, 126(1): 014101.

[本文引用: 1]

Berendsen H J C, Postma J P M, van Gunsteren W F, et al.

Molecular dynamics with coupling to an external bath

[J]. The Journal of Chemical Physics, 1984, 81(8): 3684-3690.

[本文引用: 1]

Parrinello M, Rahman A.

Polymorphic transitions in single crystals: a new molecular dynamics method

[J]. Journal of Applied Physics, 1981, 52(12): 7182-7190.

[本文引用: 1]

Xiao C X, Huang K, Liu H Z.

A novel constant interfacial area cell for determining the extraction kinetics of Er(Ⅲ) from chloride medium

[J]. Chinese Journal of Chemical Engineering, 2018, 26(6): 1435-1441.

[本文引用: 1]

於静芬, 刘涤非, 吴雄武, 等.

界面化学反应控制的萃取动力学研究——新协萃体系D2EHPA-MEHPA萃取Al3+动力学及反应机理

[J]. 铀矿冶, 1989, 8(2): 18-28.

[本文引用: 1]

Yu J F, Liu D F, Wu X W, et al.

Investigation on kinetics and mechanism of aluminum extraction by a new synergistic system D2EHPA-MEHPA

[J]. Uranium Mining and Metallurgy, 1989, 8(2): 18-28.

[本文引用: 1]

Levin V I, Korpusov G V, Man'ko N M, et al.

Extraction of tetravalent cerium by organic solvents

[J]. Soviet Atomic Energy, 1964, 15(2): 828-835.

[本文引用: 2]

秦启宗, 周祖铭, 杨永乐, 等.

溶剂萃取动力学的研究(Ⅰ): 磷酸三丁酯萃取硝酸铈(Ⅳ)的速率与机理

[J]. 复旦学报(自然科学版), 1981, 20(2): 148-154.

[本文引用: 2]

Qin Q Z, Zhou Z M, Yang Y L, et al.

Kinetic studies on the solvent extraction (Ⅰ):extraction rate and mechanism of cerium (Ⅳ) nitrate in the TBP-HNO3 system

[J]. Journal of Fudan University (Natural Science), 1981, 20(2): 148-154.

[本文引用: 2]

刘梅堂, 浦敏锋, 马鸿文.

阴离子Gemini表面活性剂在油/水界面行为的分子动力学模拟

[J]. 高等学校化学学报, 2012, 33(6): 1319-1325.

[本文引用: 1]

Liu M T, Pu M F, Ma H W.

Molecular dynamics simulation on adsorbing behavior of anionic gemini surfactants at decane/water interface

[J]. Chemical Journal of Chinese Universities, 2012, 33(6): 1319-1325.

[本文引用: 1]

Matsumoto R, Thompson M W, Cummings P T.

Ion pairing controls physical properties of ionic liquid-solvent mixtures

[J]. The Journal of Physical Chemistry B, 2019, 123(46): 9944-9955.

[本文引用: 1]

/

〈 〉

版权所有 © 2019 《化工学报》编辑部
北京市东城区青年湖南街13号（100011）
电话：010-64519489，010-64519362，010-64519485，010-64519490，010-64519451
E-Mail：hgxb@cip.com.cn
京ICP备12046843号-7
京公网安备 11010102001995号
本系统由北京玛格泰克科技发展有限公司设计开发