呂洪彬，左　臣，晏太紅

中國(guó)原子能科學(xué)研究院 放射化學(xué)研究所，北京　102413

?

磷酸三丁酯和磷酸二丁酯與Pu(Ⅳ)配合物的密度泛函研究

呂洪彬，左臣，晏太紅

中國(guó)原子能科學(xué)研究院 放射化學(xué)研究所，北京102413

采用B3LYP密度泛函方法，Pu(Ⅳ)和其他原子分別采用相對(duì)論有效原子實(shí)勢(shì)(RECP)和6-311 g(d,p)全電子基組，研究了Pu(NO3)4TBP2和Pu(NO3)2DBP2的幾何結(jié)構(gòu)、熱力學(xué)參數(shù)及紫外可見(jiàn)光譜。幾何結(jié)構(gòu)計(jì)算結(jié)果表明，由于配體空間位阻的影響，Pu與DBP的鍵長(zhǎng)更短，自然布居分析(NPA)結(jié)果也表明，Pu與DBP之間的電荷轉(zhuǎn)移更多，因此二者之間相互作用更強(qiáng)。熱力學(xué)分析顯示，Pu(NO3)4TBP2與HDBP發(fā)生配體交換的反應(yīng)自由能為-178.9 kJ/mol，說(shuō)明DBP與Pu之間的相互作用更強(qiáng)。Pu(NO3)4TBP2的紫外吸收光譜主要來(lái)源于硝酸根到Pu的f軌道的躍遷，TBP貢獻(xiàn)不大；而Pu(NO3)2DBP2的躍遷同時(shí)來(lái)源于硝酸根和DBP到Pu的f軌道的配體-金屬電荷轉(zhuǎn)移躍遷(LMCT)。

钚；磷酸三丁酯(TBP)；磷酸二丁酯(HDBP)；密度泛函理論

利用磷酸三丁酯(TBP)作為萃取劑的普雷克斯(PUREX)流程是核燃料后處理中最為成熟的技術(shù)。在萃取過(guò)程中TBP和稀釋劑都會(huì)在輻射條件下發(fā)生輻解，產(chǎn)生包括磷酸二丁酯(HDBP)在內(nèi)的輻解產(chǎn)物[1-2]。在實(shí)驗(yàn)和工廠運(yùn)行中存在Pu(Ⅳ)不能被充分反萃的現(xiàn)象，而有機(jī)相中TBP的輻解產(chǎn)物(HDBP)與Pu(Ⅳ)的結(jié)合被認(rèn)為是使其無(wú)法被充分反萃的重要因素之一[3]。對(duì)TBP和HDBP與Pu(Ⅳ)的配合物結(jié)構(gòu)和性質(zhì)的研究將會(huì)有助于萃取過(guò)程的機(jī)理研究和Pu(Ⅳ)回收率的提高。已有文獻(xiàn)[4-6]在實(shí)驗(yàn)上對(duì)Pu(Ⅳ)和TBP、HDBP形成的配合物組成、平衡常數(shù)等性質(zhì)進(jìn)行了研究，但由于Pu的毒性與放射性，對(duì)其結(jié)構(gòu)的研究比較少。對(duì)钚配合物的理論研究既是相對(duì)匱乏的實(shí)驗(yàn)數(shù)據(jù)的補(bǔ)充，還能夠加深對(duì)钚配合物的認(rèn)識(shí)。近年來(lái)，采用B3LYP、PBE0等泛函的密度泛函理論(DFT)方法被應(yīng)用于錒系元素配合物的結(jié)構(gòu)、性質(zhì)研究及核燃料后處理流程相關(guān)的研究中[7-10]。相對(duì)論效應(yīng)、溶劑效應(yīng)等因素的引入使得理論計(jì)算的結(jié)果精確度明顯提高。本工作擬在已有實(shí)驗(yàn)基礎(chǔ)上，采取DFT方法對(duì)Pu(Ⅳ)和TBP、HDBP形成的配合物結(jié)構(gòu)進(jìn)行研究，并對(duì)其熱力學(xué)參數(shù)進(jìn)行計(jì)算，同時(shí)對(duì)其紫外光譜進(jìn)行考察和分析，這些結(jié)果將有助于對(duì)有機(jī)相中TBP和HDBP的Pu(Ⅳ)配位性質(zhì)的理解，同時(shí)為工藝改進(jìn)提供理論基礎(chǔ)。

1　方法和基組

所有的計(jì)算均采用密度泛函理論，利用Gaussian 09軟件包[11]，選取B3LYP雜化泛函進(jìn)行不考慮自旋-軌道耦合的標(biāo)量相對(duì)論計(jì)算。對(duì)钚原子的計(jì)算考慮相對(duì)論效應(yīng)，對(duì)內(nèi)層電子采用相對(duì)論有效原子實(shí)勢(shì)(RECP)，對(duì)Pu(Ⅳ)配合物采用五重態(tài)自旋態(tài)，對(duì)U(Ⅳ)配合物采取單重態(tài)自旋態(tài)。對(duì)钚原子采用相對(duì)論有效原子實(shí)勢(shì)基組，在進(jìn)行結(jié)構(gòu)計(jì)算時(shí)，對(duì)钚原子選取德國(guó)斯圖加特大學(xué)ECP78MWB_SEG 基組進(jìn)行計(jì)算[12]，這個(gè)基組包含78個(gè)中心電子和16個(gè)價(jià)電子，在進(jìn)行自由能及紫外光譜計(jì)算時(shí)，采取德國(guó)斯圖加特大學(xué)ECP60MWB_SEG 基組，這個(gè)基組包含60個(gè)中心電子和34個(gè)價(jià)電子。本次計(jì)算中沒(méi)有考慮自旋-軌道耦合效應(yīng)。對(duì)其他原子采取6-311g(d,p)全電子基組[13]，有研究結(jié)果表明，采用B3LYP/RECP/6-311g(d,p)的理論計(jì)算方法對(duì)錒系元素及其與有機(jī)配體形成的配合物進(jìn)行計(jì)算得到的結(jié)果是可靠的[9]。所有的結(jié)構(gòu)優(yōu)化都是在氣相環(huán)境中，采用B3LYP/RECP/6-311g(d，p)的理論計(jì)算方法進(jìn)行的。在優(yōu)化結(jié)構(gòu)的基礎(chǔ)上，采用相同的方法對(duì)振動(dòng)頻率進(jìn)行計(jì)算，在氣象環(huán)境中得到配合物的自由能以及零點(diǎn)能。采用Multiwfn程序?qū)ψ匀辉与姾珊蚆ayer軌道鍵級(jí)進(jìn)行分析。在氣相環(huán)境優(yōu)化結(jié)構(gòu)的基礎(chǔ)上，利用似導(dǎo)體屏蔽模型得到有機(jī)溶劑環(huán)境下的自由能以及零點(diǎn)能。利用含時(shí)-密度泛函方法(TD-DFT)對(duì)其紫外光譜進(jìn)行模擬同時(shí)對(duì)紫外峰對(duì)應(yīng)的躍遷進(jìn)行了分析。為節(jié)省計(jì)算量，將TBP及HDBP中的丁基替換為甲基。去質(zhì)子化的HDBP用DBP表示。

2　計(jì)算結(jié)果

2.1Pu(NO3)4TBP2和Pu(NO3)2DBP2的幾何構(gòu)型

依據(jù)文獻(xiàn)[4]中钚的TBP和DBP配合物組成，優(yōu)化得到了氣相環(huán)境中Pu(NO3)4TBP2和Pu(NO3)2DBP2兩種化合物的結(jié)構(gòu)，示于圖1。

Pu(NO3)4TBP2和Pu(NO3)2DBP2的結(jié)構(gòu)參數(shù)與自然布居分析(NPA)結(jié)果列入表1。NPA分析顯示，Pu(NO3)4TBP2中Pu電荷為1.863，Pu(NO3)2DBP2中Pu電荷為1.601，說(shuō)明DBP與Pu原子之間電荷轉(zhuǎn)移更多，這導(dǎo)致了DBP與Pu的結(jié)合能比TBP與Pu的結(jié)合能更大。計(jì)算表明，去離子化的HDBP中的配位氧原子電荷(-1.116)比TBP中的配位氧原子電荷(-1.074)更多，說(shuō)明DBP中氧原子的配位能力比TBP中氧原子的配位能力更強(qiáng)。

淺藍(lán)色代表Pu，紅色代表O，藍(lán)色代表N，橙色代表P，灰色代表C，白色代表H圖1　Pu(NO3)4TBP2(a)和Pu(NO3)2DBP2(b)的結(jié)構(gòu)Fig.1　Optimized structure of Pu(NO3)4TBP2(a) and Pu(NO3)2DBP2(b)

表1　Pu(NO3)4TBP2和Pu(NO3)2DBP2的結(jié)構(gòu)參數(shù)與NPA分析結(jié)果Table 1　Structure parameters and NPA analyze results of Pu(NO3)4TBP2 and Pu(NO3)2DBP2

對(duì)配合物中Pu-O的Mayer鍵級(jí)進(jìn)行了分析[15]，Mayer鍵級(jí)越高，化學(xué)鍵的共價(jià)成分就越多。在Pu(NO3)4TBP2中Pu-OTBP的Mayer鍵級(jí)為0.28左右，Pu-ONO3的鍵級(jí)為0.27～0.43，在Pu(NO3)2DBP2中Pu-ODBP的Mayer鍵級(jí)為0.52～0.55，Pu-ONO3的鍵級(jí)為0.40～0.51，這些結(jié)果說(shuō)明Pu與O之間有一定的共價(jià)成分。Pu-OTBP的Mayer鍵級(jí)比Pu-OTBP的鍵級(jí)高，說(shuō)明二者之間的共價(jià)成分更多，因此，Pu與DBP的相互作用更強(qiáng)。Pu(NO3)2DBP2中Pu-ONO3的鍵級(jí)比Pu(NO3)4TBP2中Pu-ONO3的鍵級(jí)更高，說(shuō)明Pu(NO3)2DBP2中Pu-ONO3的作用更強(qiáng)，這主要是由于空阻原因?qū)е翽u(NO3)2DBP2中Pu-ONO3的距離更短，使得二者相互作用更強(qiáng)。與自由Pu(NO3)4中Pu-ONO3的鍵級(jí)(0.49～0.54)相比，TBP和DBP與Pu配位后，Pu-ONO3的鍵級(jí)都有減少，Pu(NO3)4TBP2中減少的更加明顯。

2.2Pu(NO3)4TBP2和Pu(NO3)2DBP2的熱力學(xué)性質(zhì)研究

實(shí)驗(yàn)結(jié)果表明，在有機(jī)相環(huán)境中，DBP與Pu(Ⅳ)配合物的穩(wěn)定常數(shù)比TBP與Pu(Ⅳ)的穩(wěn)定常數(shù)大，說(shuō)明DBP與Pu(Ⅳ)的相互作用更強(qiáng)。為了研究TBP和DBP對(duì)Pu(Ⅳ)的熱力學(xué)選擇性，首先在氣相條件下對(duì)反應(yīng)(1)的熱力學(xué)參數(shù)進(jìn)行了計(jì)算，結(jié)果列入表2。由表2可知：ΔH=20.36 kJ/mol；ΔS=44.40 J/K，說(shuō)明該反應(yīng)是熵驅(qū)動(dòng)的反應(yīng)。配體交換過(guò)程是配體解離，新配體再重新配位的過(guò)程。配體解離導(dǎo)致熵增，配體配位導(dǎo)致熵減，反應(yīng)過(guò)程中伴隨著TBP和DBP交換的過(guò)程有兩個(gè)硝酸根從配合物中解離下來(lái)，因此導(dǎo)致整個(gè)反應(yīng)的熵增加。雖然反應(yīng)熵值為正值，但是由于焓變的正值較大，所以計(jì)算得到的反應(yīng)自由能ΔG=7.13 kJ/mol，因此在這種條件下，配位反應(yīng)難以發(fā)生。

Pu(NO3)2DBP2+2TBPHNO3

(1)

ΔG=ΔG(Pu (NO3)2DBP2)+

2ΔG(TBPHNO3)-ΔG(Pu(NO3)4TBP2)-

2ΔG(HDBP)

(2)

雖然氣相計(jì)算結(jié)果表明，TBP與Pu(Ⅳ)的結(jié)合能力更強(qiáng)，為使計(jì)算結(jié)果與實(shí)際狀態(tài)更加接近，因此在考慮溶劑情況下進(jìn)行了計(jì)算。因?yàn)镻u萃取實(shí)驗(yàn)的有機(jī)相通常為煤油或正十二烷，因此采用正己烷作為模擬溶劑模型，計(jì)算結(jié)果列入表2。由表2可知：與在氣相條件中計(jì)算結(jié)果不同，在溶劑條件下計(jì)算得到反應(yīng)(1)的ΔH=-154.97kJ/mol，在溶劑條件下反應(yīng)焓變轉(zhuǎn)為負(fù)；ΔS=80.27J/K，在溶劑條件下熵變?cè)龃?，更有利于反?yīng)的進(jìn)行，由此計(jì)算得到的反應(yīng)自由能ΔG=-178.89kJ/mol，在液相中反應(yīng)轉(zhuǎn)變?yōu)闊崃W(xué)有利的反應(yīng)。熱力學(xué)計(jì)算結(jié)果說(shuō)明，DBP能夠和Pu(NO3)4TBP2配位，說(shuō)明與TBP相比，DBP與Pu(Ⅳ)的配位能力更強(qiáng)。

表2　利用B3LYP/6-311G(d,p) 計(jì)算得到的反應(yīng)(1)的熱力學(xué)參數(shù)Table 2　Calculated thermodynamic parameters of reaction (1) by B3LYP/6-311G(d,p) methods

2.3Pu(NO3)4TBP2和Pu(NO3)2DBP2吸收光譜的理論研究

利用TD-DFT方法對(duì)氣相中Pu(NO3)4TBP2和Pu(NO3)2DBP2的紫外光譜進(jìn)行了模擬，分子軌道示于圖2、3。Pu(NO3)4TBP2的吸收峰主要為377、407、581、607 nm，分別對(duì)應(yīng)HOMO-2到LUMO+2、HOMO-3到LUMO、HOMO-2到LUMO和HOMO-1到LUMO的躍遷。由圖2可知，Pu(NO3)4TBP2的HOMO-3、HOMO-2、HOMO-1和HOMO軌道其主要由硝酸根中的N、O原子的p軌道構(gòu)成，LUMO軌道主要由Pu原子f軌道組成，因此其紫外光譜主要對(duì)應(yīng)于硝酸根中的N、O原子的p軌道到Pu原子f軌道的配體-金屬電荷轉(zhuǎn)移躍遷(LMCT)，TBP分子的軌道對(duì)紫外光譜躍遷貢獻(xiàn)不大。

圖2　氣相中Pu(NO3)4TBP2的分子軌道Fig.2　Molecular orbitals of Pu(NO3)4TBP2 in gas phase

Pu(NO3)2DBP2的吸收峰主要在405、420、432、487、513、612、630 nm，分別對(duì)應(yīng)于HOMO-4到LUMO+2、HOMO-2/HOMO-1到LUMO+2、HOMO-6到LUMO/HOMO-3到LUMO+1、HOMO-2/HOMO-4到LUMO+1、HOMO到LUMO+1、HOMO到LUMO的躍遷。由圖3可知，HOMO-6、HOMO-4、HOMO-3、HOMO-1主要由DBP中原子軌道組成，HOMO-2和HOMO主要由硝酸根中原子軌道組成，LUMO、LUMO+1和LUMO+2主要由Pu的f軌道組成，因此405、420、32 nm處吸收峰主要對(duì)應(yīng)于DBP到Pu的配體-金屬電荷轉(zhuǎn)移躍遷(LMCT)，487 nm吸收峰由DBP到Pu的LMCT和硝酸根到钚的LMCT混合而成，513、612、630 nm則主要是硝酸根到钚的LMCT。

圖3　氣相中Pu(NO3)2DBP2的分子軌道Fig.3　Molecular orbitals of Pu(NO3)2DBP2 in gas phase

3　結(jié)　論

利用B3LYP密度泛函理論方法，采用Gaussian 09程序和相對(duì)論贗勢(shì)基組，對(duì)Pu(NO3)4TBP2和Pu(NO3)2DBP2兩種化合物進(jìn)行了結(jié)構(gòu)優(yōu)化。結(jié)果表明：Pu(NO3)4TBP2的優(yōu)化結(jié)構(gòu)與實(shí)驗(yàn)結(jié)果相吻合，Pu(NO3)2DBP2的優(yōu)化結(jié)構(gòu)中，Pu與DBP和硝酸根的鍵長(zhǎng)都比Pu(NO3)4TBP2中Pu與TBP和硝酸根的鍵長(zhǎng)要短，相互作用更強(qiáng)。在優(yōu)化結(jié)構(gòu)基礎(chǔ)上進(jìn)行NPA分析，發(fā)現(xiàn)在Pu(NO3)2DBP2中，DBP到Pu的電荷轉(zhuǎn)移比Pu(NO3)4TBP2中TBP與Pu的電荷轉(zhuǎn)移更多，且DBP與Pu之間的Mayer鍵級(jí)比TBP與Pu的鍵級(jí)更大。利用TD-DFT方法對(duì)紫外可見(jiàn)光譜進(jìn)行了模擬并對(duì)吸收峰所對(duì)應(yīng)的軌道躍遷進(jìn)行了分析。對(duì)紫外光譜的比較可以發(fā)現(xiàn)，Pu(NO3)2DBP2在400～430 nm處的紫外吸收峰更強(qiáng)，這一部分主要對(duì)應(yīng)于DBP到Pu的LMCT躍遷，而在Pu(NO3)4TBP2中，TBP分子則對(duì)紫外吸收貢獻(xiàn)較少。

[1]Mincher B J, Elias G, Martin L R, et al. Radiation chemistry and the nuclear fuel cycle[J]. J Radioanal Nucl Chem, 2009, 282(2): 645-649.

[2]Bellido A V, Rubenich M N. Influence of the diluents on the radiolytic degradation of TBP in TBP, 30% (V/V)-diluent-HNO3systems[J]. Radiochim Acta, 1984, 36: 61-64.

[3]Mincher B J, Mezyk S P. Radiation chemical effects on radiochemistry: a review of examples important to nuclear power[J]. Radiochim Acta, 2009, 97(9): 519-534.

[4]May I, Taylor R J, Wallwork A L, et al. The influence of dibutylphosphate on actinide extraction by 30% tributylphosphate[J]. Radiochim Acta, 2000, 88(5/2000): 283-290.

[5]Sokhina L P, Solovkin A S, Teterin E G, et al. Dibutylphosphates of tetravalent plutonium, zirconium and thorium, formed in solutions of tri-n-butyl phosphate (TBP)[J]. Radiokhimiva, 1978, 20(1): 28-31.

[6]Sulka M, Cantrel L, Vallet V. Theoretical study of plutonium(Ⅳ) complexes formed within the PUREX process: a proposal of a plutonium surrogate in fire conditions[J]. J Phys Chem A, 2014, 118(43): 10073-10080.

[7]王東琪,van Gunsteren W F.錒系計(jì)算化學(xué)進(jìn)展[J].化學(xué)進(jìn)展,2011,23(7)：1566-1581.

[8]Steele H, Taylor R J. A theoretical study of the inner-sphere disproportionation reaction mechanism of the pentavalent actinyl ions[J]. Inorg Chem, 2007, 46(16): 6311-6318.

[9]Xiao C L, Wu Q Y, Wang C Z, et al. Quantum chemistry study of uranium(Ⅵ), neptunium(Ⅴ), and plutonium(Ⅳ,Ⅵ) complexes with preorganized tetradentate phenanthrolineamide ligands[J]. Inorg Chem, 2014, 53(20): 10846-10853.

[11]Frisch M J，Trucks G W，Schlegel H B，et al. Gaussian 09, Revision A.1[R]. U S A: Gaussian, Inc. Wallingford, CT, 2009.

[12]Kuchile W, Dolg M, Stoll H, et al. Energy-adjusted pseudopotentials for the actinides. parameter sets and test calculations for thorium and thorium molecules[J]. J Chem Phys, 1994, 100: 7535-7542.

[13]Dithchfield R, Hehre W J, Pople J A. Self-consistent molecular orbital methods extended gaussian-type basis for molecular-orbital studies of organic molecules[J]. J Chem Phys, 1971, 54: 724-728.

[14]Den Auwer C, Revel R, Charbonnel M C, et al. Actinide coordination sphere in various U, Np and Pu nitrato coordination complexes[J]. J Synchrotron Radiat, 1999, 6: 101-104.

[15]Lu T, Chen F W. Multiwfn: a multifunctional wavefunction analyzer[J]. J Comput Chem, 2012, 33(5): 580-592.

Density Functional Theoretical Analysis of Complexes of Tributylphosphate and Dibutylphosphate of Pu(Ⅳ)

LU Hong-bin, ZUO Chen, YAN Tai-hong

China Institute of Atomic Energy, P. O. Box 275(26), Beijing 102413, China

The geometry, thermodynamic selectivity and UV-vis spectra of Pu(NO3)4TBP2and Pu(NO3)2DBP2in the gas phase were investigated by using B3LYP density functional theory method with the RECP basis for Pu and 6-311 g(d,p) basis for other atoms. The structure optimization results show the bond length between Pu and DBP is shorter than the length between Pu and TBP for the steric hinder of the ligands, and NPA charge analysis shows more electron trasnfer between DBP and Pu ion indicating the stronger coordination ability of DBP. The free energy of ligand exchange between Pu(NO3)4TBP2and HDBP was -178.9 kJ/mol, indicating the higher coordination affinity between DBP and Pu. The UV-vis absorption spectra of Pu(NO3)4TBP2in the gas phase were pmainly from nitrate to f orbital of Pu, TBP ligand made few contribution, yet the absorption spectrum of Pu(NO3)2DBP2were both from nitrate to f orbital of Pu and the ligand-to-metal charge transfer between DBP and Pu.

plutonium; TBP; HDBP; density functional theory

2016-03-03；

2016-09-07

O634.12

A

0253-9950(2016)05-0282-06

10.7538/hhx.2016.38.05.0282