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

2022-04-26李春晖何辉何明键张萌高杨矫彩山

化工学报 2022年4期

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

（1 哈尔滨工程大学核科学与技术学院，黑龙江哈尔滨 150001； 2 中国原子能科学研究院放射化学研究所，北京 102413）

引言

离子液体是由阴阳离子构成，在室温下呈液态的有机盐类化合物。以离子液体为溶剂的萃取体系对核电站乏燃料水法后处理中关键核素通常表现出更高的萃取效率和选择性，但其萃取行为和机理通常更为复杂[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 恒界面池装置Fig.1 Diagram of constant interfacial area cell

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

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

1.3 分子动力学模拟

通过分子动力学模拟研究C4mimNTf2+ 1,2−二氯乙烷/H2O+HNO3体系的相界面性质。为考虑水相硝酸形态的影响，模拟设定50%硝酸以HNO3分子形式存在，50%硝酸电离成H3O+和NO−3[23−24]。其中，水分子选择TIP3P 模型[25]，H3O+力场参数及原子电荷取自文献[26]。C4mim+、NTf−2、TBP、1,2−C2H4Cl2、HNO3、NO−3的力场参数均取自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+和NO−3个数均为216；有机相中1,2−二氯乙烷分子个数为2000，C4mim+和NTf−2个数为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 萃取两相体系示意图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(b)表明Ce(Ⅳ)的萃取速率r0随两相界面面积A增加而增加，二者间呈正比关系且经过原点，因此该萃取过程受界面过程影响，即萃取过程为界面化学反应控制过程[40]。

图3 萃取模式判定Fig.3 Determination of extraction control mode

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

式中，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 水相硝酸浓度对萃取速率的影响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+和NO−3间形成的配合物主要以Ce(NO3)2−6形式存在[42]，并进一步与H+结合生成可萃的HCe(NO3)−6，水相中的配位取代反应对萃取速率的影响减弱，水相内可萃组分扩散过程逐渐成为水相内的控制步骤。

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

图5 水相氢离子浓度对萃取速率的影响（水相硝酸根浓度c(NO−3)=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 水相硝酸根浓度对萃取速率的影响（水相氢离子浓度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(Ⅳ)萃取速率随水相NO−3浓度增加而降低。这是由于NO−3与Ce(Ⅳ)阴离子配合物组分间的竞争作用引起。

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

图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 有机相C4mimNTf2浓度对萃取速率的影响Fig.8 Dependence of extraction rate on organic C4mimNTf2 concentration

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

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

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

图9 水相NTf−2 浓度对萃取速率的影响Fig.9 Dependence of extraction rate on aqueous NTf−2 concentration

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

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

式中，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 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，通过P O 基团上的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+和NTf−2在相界面附近存在显著的吸附现象；同时，图11(b)的电荷密度分布进一步表明，C4mim+和NTf−2中的极性区域（带电荷基团）倾向于朝水相方向分布，而非极性区域则倾向于朝有机相方向分布。

图11 C4mim+和NTf−2 数密度及电荷密度在z轴方向分布Fig.11 Number and charge density profile of C4mim+and NTf−2 along z axis

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

图12 C4mim+和NTf−2 质心间径向分布函数Fig.12 Center of mass radial distribution function between C4mim+and NTf−2

综上所述，结合萃取动力学实验及分子动力学模拟，可以认为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）NTf−2从界面扩散至水相主体。萃取过程为水相扩散和界面反应控制的“混合模式”，控制步骤为水相扩散过程（2），以及界面反应（4）。

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