铀与氢气的反应动力学

中国科学G辑物理学力学天文学2005,35(4):415-423415铀与氢气的反应动力学刘晓亚朱正和蒋刚陈涵德(①中国工程物理研究院核物理与化学研究所,绵阳621900:②四川大学原子分子物理研究所,成都610065)摘要运用相对论有效原子实模型RECP和B3LYP方法计算了UH和UH2的结构.并在此基础上,利用多体项展式理论方法,导出三原子分子UH2的解析势能函数.根据UH2势能函数,采用准经典的 Monte- Carlo轨线法研究了U+H2体系的碰撞过程.依据反应截面与能量之间的关系,指出U+H2(v=j=0是无阈能反应,主要生成UH2,而反应H+UH(v=j=0)生成UH2的反应截面相对来说要小得多,可以忽略.如此证明铀与氢气反应生成UH3的中间产物是UH2,而不是UH这些结果为金属铀的抗氢化腐蚀研究提供了理论基础关键词UH2势能函数分子反应动力学铀的氢化腐蚀自1789年由 Klaproth从沥青铀矿中发现铀以来,铀及其化合物越来越受到人们重视.它具有良好的金属特性和独特的核性能,不仅在国防上,而且在人类开发和利用核能中起着极其重要的作用.但是,铀的化学性质非常活泼,容易和环境气氛如O2,H2和HO2汽发生反应.铀长期处在H2气氛中表面会发生严重的腐蚀,其核性能会大大受到影响.目前,国外已在U-H2(ν)体系的表面氧化反应进行了较为深入的研究4.长期以来人们虽然在氢化反应动力学和反应机理研究方面进行了大量的研究⑤8.但过去的研究工作主要着眼于铀金属块(bulkuranium)与H2在高温0~500℃)、高压(102-10Pa)条件下的氢化反应.对铀金属在H气氛中的表面化学反应尚缺乏系统研究,用原子分子微观反应动力学研究铀金属与H2的反应机制更是少有报道.国内在上述领域的研究工作目前刚刚开展.1UH2基态的分子结构与势能函数1.1UH2分子平衡结构与振动频率U原子采用14个价电子(6s27s36p6d5f的相对论有效原子实势及(5413d4)2004-08-03收稿,200505-16收修改稿中国煤化工国家表面物理和化学重点实验室基金资助项目(批准号联系人,E-mai:zhuxi@scu.edu.cn;peonyzhu@sina.comHCNMHGSCIENCE IN CHINA Ser. G Physics, Mechanics Astronomy416中国科学G辑物理学力学天文学第35卷3S3p2d2收缩价基集合,H原子采用非相对论6-311+G*全电子基集合,在U原子(RECP近似下,用B3LYP( Becke3个参数交换函数与Lee- Yang-Parr相关函数组成的杂化密度泛函理论(DFT方法)计算了体系基态的平衡几何、谐性力常数和谐振频率υ1,υ2及υ3.计算表明UH2分子最稳态为C2构型,具有B2对称性的电子态,平衡核间距RuH=0.21325m,∠HUH=1084146°,如表1.表1UH2(XB,)的结构与性质参数平衡结构RuH=0.21325nm,∠HUH=108.4146°谐振频率U1=308.1526cm,u2=1308.5678cm-,u3=1329.806cm离解能5.1635eV力常数f(HU)=0656318,fR(UH)=0.065313.fR=0.18152×102fa=-0.372984×102,fR=-0.3743939×10-2,fa(∠HUH)=0.0209712UH2分子多体项展式的势能函数UH体系的多体项展式理论方法的势能函数应写为v(R1,R2,R3)=Vu(2(R1)+vu2(R2)+Vm2(R3)+Vm3(R1R2R3其中R1=R2=RuH,R3=Rm,(1)式中的两体项采用 Murrell- Sorbie势能函数UH分子的两体项势能函数数据见表2,H2(Xx)的势能函数引用文献[0的值.(1)式中的vum23(R1,R2,R3)为三体项,有10个线性系数(C,C1-C和3个非线性系数(Y,γ2,γ),对势能表面进行非线性优化,确定出3个非线性系数,而10个线性系数由10个已知条件确定,根据表1,UH2分析势能函数的参数可由此而确定,列于表3表2UH(Xm)分子势能曲线参数R/nmUH(XTI)1.77240.227321980123.38369.4表3UH2(⑧B2)分析势能函数的三体项参数7501494×102.4138476×106.7268183×10-31614000×10-2.125886×102.1222724×11.270089X102946270×10C9-8.69857X10-2y=0.9p=0.63=3.5图1和2是根据UH2的分析势能函数绘制中国煤化工地再现了UH2的结构特征.图1是U-H键与UH键的伸HCNMHG们可以精SCIENCE IN CHINA Ser. G Physics, Mechanics Astronomy第4期刘晓亚等:铀与氢气的反应动力学4171为-5153eV;2为-49cV3为-45ev4为-3eV5为-2eV6为-1.5eV;7为-ev0.20图1UH2的伸缩振动势能图0.20.12为-46eV3为-35eV4为-2.5eV0.05为-1.72ev6为-1.00eV7为-4.6ev0.10.2.2-0.10.00.10.2X/nm图2UH2的旋转势能图确地找到UH2(ⅹB)的平衡结构(RuH=0.21325mm)和它的离解能De=5,16356cV而且在两个等价的通道上U+HH→UH2没有鞍点存在,这说明U原子与H2的反应是没有阈能的放热反应图2是把RH=0.34569mm固定在ⅹ轴上,U原子绕H-H旋转的等值势能图.图中仍然再现了UH2(文B2)的平衡结构的特征,有两个等价的C2V角型极小(Y=±0.128nm,X=0.0nm),在这两个极小之间有一个线性鞍点(X=Y=0.0nm,E≡-3.0cV)存在于H-H键的中间和两个等价的线性鞍点(X=±0.3945nm,Y=00nm,E=-1.68eV)在H2的两边.这说明:如果U原子要穿过H2的中间进行内迁移需要克服2.16356εⅴ的能垒;如果穿过H的两边进行內迁移则需翻越3,.484eV的能垒,以此可以看出U原子能比较容易拉长HH键而生成UH.所有这些合理地再现了UH2的C2结构特征,为UH2体系THE中国煤化工了数据CNMHGhina co418中国科学G辑物理学力学天文学第35卷2UH2体系的分子反应动力学过程选用前面计算得到的UH2分子的全空间解析势能函数,建立12个联立的微分方程组,采用 Monte-Carlo半经验轨线法( General Trajectory Program)求解得到原子间距随时间的变化关系(即轨线),通过对这些轨迹的统计分析可以得到反应产物的分布与初始平动能的关系21U(L)+H2(X)碰撞的分子反应动力学过程对于UCL)+H2(X∑t)体系的碰撞反应产物有5种通道H('S)+UH by通道UCLu+HH(o*2通道U+H(Xo:→{H2S2)+UH(t3通道UHH( B2)4通道,UL)+H(2S2)+H(S2)通道当U原子的初始平动能为E=0.4184-8368k/mo时,U+HxX2)反应产物的结果见表4.从表4可知:当U原子的初始平动能在E1=04184-271.96kJ/mol时,U+HxXΣ反应生成络合物和非弹性散射产物.当U原子的初始平动能在E27196-836.8kJ/mo时,U+H2(X2)反应生成UH和非弹性散射产物当U原子的初始平动能E1≥6276k/mol时,U+H2(XΣ反应产物U+H+H开始出现.这可以由UH2的解析势能函数分析得出初步的结论.从UH2的旋转图可知,由于在U向H2分子接近生成UH2时,无论从哪个方向都没有能垒出现,也就是说U+H2→表4U+H2(X∑反应产物的分布E/kJmo1通道2通道3通道5通道26664580251.04292.884874355.6447794858103836.8YH出中国煤化工204NMHSCIENCE IN CHINA Ser. G Physics, Mechanics Astronomy第4期刘晓亚等:铀与氢气的反应动力学419UH2的反应是无阈能的放热反应.所以在较低碰撞能量下,主要以络合物和非弹性碰撞产物为主,碰撞能量稍高一点后,U原子就有可能把H2分子的键打断,而与H原子重组导致交换反应的发生;在较高碰撞能量下,U原子与H原子还没来得及重组就很快完全分解,则导致第5通道的产物U+H+H产生为了确定反应UXL)+HH(v=j=0)→UH+H以及UH2的碰撞截面,计算了大量不同初始状态的碰撞轨线,使其具有统计分布的规律.在确定反应截面σr时,取N为5000条表5和6以及图3和4记录了U+H2UH2和U+H2→UH+H的反应截面表5U+H2→UH2的反应截面E/kJ·mol104.6251.040.9210.7780.6850.6170.2350.2050.1390.16513330.8230.2820.0890003000040.00021.9×10-3表6U+H2→UH+H的反应截面E/kJ·mol27196334.72355,64418,4627,6836.80.2000.1660.1500.1362.9413.1654.3974.8532.4573.3211.7661.25一最大磁撞参数/nm磁撞截面/nm0.500.000.25020406080100120140160180200220240260280E/kJ. mol图3U+H2→UH2的反应截面在所计算的能量范围内,U(L)+H(X)的反应碰撞主要表现为非反应碰撞,在0.4184-271.96kJmo能量范围有络合物生成,在271,96-836.8kJ/mol能量范围有交换反应产生.在低能范围(0.4184-271.96kJ/mol)反应U+H2→UH2是强烈的无阈能放热反应,在碰撞能量到271.96 kJ/mol时开始有UH+H产生(交换反应),但是交换反应的截面远小于生成络合中国煤化工低能时络合物很容易生成,而交换反应的发生相对要困HCNMHGwww.scichina.com420中国科学G辑物理学力学天文学第35卷最大磁撞参磁撞截面E/kJ. mol图4U+H2→UH+H的反应截面22H(S)+UH(Xm碰撞的分子反应动力学过程对于H(S)+UHX)体系的碰撞反应产物有5种通道H(S)+UH(X"II)通道H(S)+UH(X D)2通道H(Se)+UH'(X"ID)U()+H,(x'=g3通道UH2(X B2)4通道UCLu)+H(S)+H(SE)5通道当U原子的初始平动能为E=04184-8368kJ/mol时,U+H(xX2)反应产物的结果见表7表7H2Sn)+UH(X1m反应产物的分布E/kJ·m1通道2通道3通道4通道5通道0.41841503346l4.184339l1373197198l41.840.0.685167.360.54141560.5323514561836.84902中国煤化工CNMHGSCIENCE IN CHINA Ser. G Physics, Mechanics Astronomy第4期刘晓亚等:铀与氢气的反应动力学421从表7可见,随H原子初始相对平动能的增加,交换反应1通道的数量逐渐增加,在E1>41.84kJ/mol以后又逐渐降低,而交换反应3通道的数量随H原子初始相对平动能的增加逐渐減少;非反应2通道的数量随H原子初始相对平动能的增加而逐渐增加;反应4通道(生成络合物的通道)的数量随H原子初始相对平动能的增加而增加,稍后又减小,对此特别计算了一下:当H原子初始相对平动能到了E=34.3088kJ/mol时就没有络合物生成了;第5通道(碰散的通道)在E>167.36 kJ/mol后,逐渐增加而后开始减少,其主要原因是由于UH的离解能为1.7724eV,只要H原子的初始平动能大于171.0kJ/mo就有碰散的可能性,因此到E>167.36kJ/mol时第5通道就有反应了.从整个计算结果看,在低能部分0.4184-41.84kJ/mol,H原子把U置换出来相对其他反应较为容易;而在高能部分41.84~836.8kJ/mol,非反应碰撞占主要地位,即H+UH碰撞主要产物为U+H2和非反应碰撞.H(S)+UH(v=j=0体系的反应截面σ计算结果列于表8中,初始相对平动能E对反应截面σ和最大碰撞参数bmnx的关系曲线如图5表8H(S2)+UH(v=j=0体系碰撞反应的E和o通道3通道4通道0.41840.39321.5034.9121.5030.072412553.16380.097001060.910.508l0.6850.20970.64840.1012167.360.54140.54130.06440.53190.53180.3911627.60.28730.0015836.80.2790.00020.30430.0001在所计算的能量范围内,H(S2+UH(Xm)的反应碰撞主要表现为交换反应H(S2)+UH(XI)→UXL)+H(XΣ和非反应碰撞,在0.4184-41:84kJ/mo低能量范围伴随有HS)+UH(Xm)→UH(X1+H(S2及络合物生成,在4848368kJ/mol能量范围主要为非反应碰撞.H(S2)+UH(XI)UXL)+H2(Xt交换反应3通道反应)是强烈无阔能的放热反应,H(S2)+UH(XTm→UH(XTm+H(S2)反应及生成络合物的反应,虽然在低能时有少量的发生,但是其截面远小于3通道的反应截面.这说明H(S2)+UHXm)→U(XL)+H2xX2交换反应容易得多,而其他反应相对要困难一些中国煤化工CNMHGhina co42中国科学G辑物理学力学天文学第35卷一最大磁撞参数1通道3通道4通道0100200300400500600700800900E/k. mol图50-E1与bmxE1关系曲线3结果与讨论用半经验轨线法考察了UH2体系的分子动力学过程.在计算结果的基础上提岀了铀的氢化腐蚀的机理.这对铀的抗腐蚀研究具有较为重要的意义.计算了U(L)+H2(X∑,H(S)+UH(XT两体系的分子动力学过程.对于U(L)+H2(X碰撞反应,U(L+H1(X的反应主要生成UH2(x3B2),基本不生成H(S2)+UH(XT).对于HS)+UH(XT碰撞,反应主要为第3通道的反应即H(S2)+UHXT)→U(L)+H2(X,以及非反应碰撞.此计算结果说明:铀材料在氢气氛中容易生成UH2这种中间产物,多数的文献山报道,U(La)+H2(X2)反应的最终产物是UH,但我们可以推断:U与H2反应生成UH3的过程为两步:(i)U(L)+HX→H2(X3B2),(i)UH+H2→UH2,而不是U(L)+H(X∑→H(2S2)+UH(Xm,UH+H2→UH因为从计算结果可知:U(L)+H2xXΣ)反应生成UH的反应截面非常小,基本不发生反应,U与H2反应生成UH2要容易得多.这一论断从前面的UH2势能函数分析也可证明这一点,另一方面,由于H2的离解能为4.746eV,U原子接近H2时,它较难于离解H2生成UH,而拉长H一H键重组生成UH2要相对容易得多,所以U+H2UH2→UH2+H2UH3是合理的过程,而不是U+H2→UH+H→UH+H2→UH参考料I Colmenares C A. Oxidation mechanisms and catalytic properties of the actinides. Prog Solid State Chem1984,15:257~3642 Allen g c. stevens jCH The behavior of uranium metaI Cham Soc, FaradayTrans I,1988,84:165-174中国煤化工3 Alirer m. mueller b A. Peterson C L. et al. reactionsCNMHGChem PhysSCIENCE IN CHINA Ser. G Physics, Mechanics Astronomy第4期刘晓亚等:铀与氢气的反应动力学421970,52:37-464 Condon J B, Larson E A Kinetics of the uranium-hydrogen system J Chem Phys, 1973, 59: 855-8655 Stakebake J L. Kinetics for the reaction of hydrogen with uranium powder. J Electronhem Soc, 1979, 1261596~16016 Condon JB, Calculation of oxidation rate of zirconium Corrosion 1981, 37: 721Bloch J, Mintz M H. Type of hydride phase development in bulk uranium and holmium. J Nucl Mater982,110:251~258 Bloch J, Cimca F, Kroup M, et al. The initial kinetcs of uranium hydride formats studied by a hot stagemicroscope technique. J Less-Common Met, 1984, 103: 165-171Hay P J, Richard L M. Theoretical studies of the structures and vibrational frequencies of actinompounds using relativistic effective core potentials with Hartree-Fock and density functional methodUF6 and PuF6. J Chem Phys, 1998, 109: 387510朱正和,俞华根.分子结构与分子势能函数.北京:科学出版社,199711 Libowitz GG. The Actinide Hydrides. In: Mueller W, Blackedege J P, Libowitz GG, eds Metal HydridesNew York: Academic Press, 196812 Cordfunke E H P In The Chemistry of Uranium. Elsevier Pub Co Amsterdam, 196913 Katz J J, Rabinowitch E, In The Chemistry of Uranium. New York: Dover Pub, Inc, 19514 Bennett MJ, Price JB. The oxidation behavior of uranium in air at 348-765. J Nucl Mater, 1981. 1015 Stakebake J L Kinetics for the reaction of hydrogen with uranium powder. J Electerchem Soc, 1979, 126495-49616 Toteneier T C. A review of corrosion and pyrophoricity behavior of uranium and plutonium, ArgonneNational Laboratory Technical Report ANW/ED/95-2, 19957 MulfordR NR, Ellinger F H. A new form of uranium hydride. J Amer Chem Soc, 1954, 76: 297-298中国煤化工CNMHGhina. com