DTO体系分子反应动力学

第25卷第6期原子与分子物理学报2008年12月JOURNAL OF ATOMIC AND MOLECULAR PHYSICSDec.2008文章编号:10000364(2008)06-143908DTO体系分子反应动力学阮文12,谢安东1,余晓光1,罗文浪12,朱正和2,傅依备2(1.井冈山大学数理学院,吉安343009;2,四川大学原子与分子物理研究所,成都610065摘要:基于DO分子(X1A1)的氢同位素效应,得到修正的 Born-Oppenheimer(BO)理论下多体展式分析势能函数,用准经典的 Monte-Carlo轨迹法研究了D+TO(0,0)和T+DO(0,0)的分子反应动力学过程结果表明:在碰撞能量较低时(<209.2km-1),D+TO(0,0)和T+DO(0,0)反应主要生成O+DT,并且该反应是无阂能的;有少量的交换反应产物DO(TO)生成,并伴有极少量的络合物产生,碰撞能大于209.2 kI.mol-后,逐渐出现分子被完全碰散成D,T,O原子的情形,反应D+TO(0,0)→OD+T和T+DO(0,0)→OT+D是无热反应但是有阂能存在,由于D和T原子的同位素效应,置换产物轨线存在非致性关键词:DTO;势能函数;分子反应截面;轨线非对称性中图分类号:O643.1文献标识码:ATheoretical study on molecular reaction dynamics of the DTo systemRUAN Wen,2, XIE An-Dong, YU Xiao-Guang, LUO Wen-Lang, 2, ZHU Zheng-He2, FU Yi-Bei2(1. College of Mathematics and Physics, Jinggangshan University, Ji'an 343009, China:2. Institute of Atomic and Molecular Physics, Sichuan University, Chengdu 610065, China)Abstract: Based on the hydrogen isotopic effect and its nuclear motion, the analytical potential energy functionAPEF)of Dto(X' A1)is derived by many-body expansion method(MEM)under the correctional Born-Oppenheimer(B-O). Then the atomic and molecular reaction dynamics for D+ TO and T+DO have beenstudied on the potential energy function DTo(X'A1)by Monte-Carlo quasi-classical trajectory approach. Itis shown that the main reactions of D+ To--DT+O and T+Do--DT+O along with a few long-lived com-plex compounds have no threshold energy with low collision energy, which agrees with the potential surfaceWhen the collision energy is greater than 209.2 k/mol", the DTO molecules will decompose into D, T and Oatom completely, and these reactions have threshold energy. The other reactions of D+ToTO+D and T+DO-DO+ T are the reactions of non-decalescence, but they have threshold energy. However, the trajectoriesnd collision sections of reactions D+ To (0, 0)and T+ Do(0, 0) present variance distribution due to the iso-topic effect of deuterium and tritium atom, which corresponds to the potential surfaceKey words: DTO, potential function, reactive section, trajectory asymmetrical distribution收稿日期:2008-0508基金项目:国家自然科学基金(10676022)资助项目;江西省教育厅科技研究ayH中国煤化工CNMHG省科技厅工业攻关项目[200621];井冈山大学校级自然科学研究项目作者简介:阮文. E-mail: ruanwensongyan@126.cm通讯作者:朱正和.Emal:zhum@scu.edu.cn原子与分子物理学报第25卷1引言量和核振动能计算,用核振动能修正电子能量获得DTO分子的势能值.将势能函数展开成键长和键众所周知,氚对能源、环保、核材料和军事具有角的泰勒级数19,用这组修正后的势能值,通过最非常重要的意义,如何从重水中提取氚是当今的一小二乘法拟合,得到泰勒级数二次项的系数,则可大主要问题目前,从重水中提氚有多种方法,如氢获得DIO分子的谐性力常数计算结果见表1水同位素交换-3,由D2交换分离出DT,是通过势能函数采用多体项展式理论确定41.令反应D2(气)+DTO(蒸汽)→DT(气)+D2O(蒸R1,R2和R3分别表示DTO分子的DO,TO和汽).长期以来,人们对水及其同位素分子进行了大DT的键长,并设基态原子能量为零,则满足离解量的研究, Cheng等(4采用同步加速辐射测定了氢极限的多体项展式分析势能函数为同位素水分子在室温下的吸收光谱;AtlG等5V=wg(D)+VB(R1,A2∑”)采用高精度(包括核、价键、相对论、量子电动力学以及 Born-Oppenheimer修正对角化)的完全势能V(R2,A2∑)+VB(R3)+面( CVRQD)方法获得了H2O和D2O分子的平衡vBor(ri, R2, R3)几何结构;wang等6研究了HO分子的势能面及选择基态原子的能量为零,v8(D2)为激发态原OH+D→OD+H的交换反应速率常数;Pab等]子O(D2)的能量,其值为1.967eV0,R1=R2采用 CVRQD方法获得了H2O及氢同位素水分子=R,R3=Rm,两体项选用Mure! Sorbie势能的平衡几何结构、正则振动频率及H2O分子的势函数1:114能面;我们曾采用密度泛函B3lyp6-311++gV=-De(1+a1P+a2p2+ a3p)exp(-a1p)方法,使用多体项展式理论,结合同位素效应得到式中p=r-x2,r,分别为核间距和其平衡距离,了 Born-Oppenheimer(BO)近似理论下由同位素De为离解能,a1,a2,a3为势能函数参数两体项效应修正的解析势能函数8.事实上,O在微量参数见表2DT,以及D在微量OT存在下的反应过程及机理三体项vBr(R1,R2,R3)一般可写成十分复杂,有很多问题还未能给出明确回答因此,VEOr(Ri, R2, R3)=PT(1)进一步开展有关理论研究,对深入认识这种复杂反式(1)中的P为(R1,R2,R3)的多项式.T为量程应体系机理有积极意义函数,量程函数的一种可能选择为量子力学从头计算是研究分子结构及其相作T=Ⅱ[1-tanh(yS2)],用势分子微观反应动力学及化学反应过程机理的式(2)中%为参数,S1为对称内坐标,当S趋于无有效方法,理论研究结果与实验结果基本一致而穷大时,T趋于零,从而保证当任何一个原子移动用密度泛函理论方法,研究分子体系的势能函数,无穷远时,VB必为零的条件进而研究反应体系微观分子反应动力学过程,已经由于三体项需要表示为对称内坐标的函数,选为不少研究者成功应用9-18.本文利用密度泛函用C2构型的对称变换换矩阵21,以R1:=1.528理论结合同位素效应和多体展式理论方法8),得A,R2=R32=0.962A为参考点,则内坐标向对到基态DTO分子的结构与分子势能函数,基于准称内坐标的变换如下经典分子反应理论,研究反应体系微观分子反应动力学过程,进一步研究DTO微观反应过程机理S2=01/21/22同位素效应修正下的势能函数s301/2-1/2人(p3P2=R1-R(i=1,2,3)本文通过计算同位素分子核的振动能,修正则式(1)中P和T变换为S的多项式其形式为:BO近似下电子能量的方法来计算DTO的力常数8,运用B3yp/6-311++g“方法,在DTO平衡+C1S1+C2S2+C3S2+C4S3+键长、键角值两边邻域内选取若干个点,使这些点中国煤化工C7S+CsS1(s+处于近似谐性振动范围,然后逐一进行单点电子能CNMHG第6期文等:DTO体系分子反应动力学T=[1-tanh(Y1 S:/2)][1-tanh(2S2/2)][1-tanhr3 S3/2)三体项参数见表3.表1基态DTO的结构与力常数Table 1 Equilibrium Structure and Force constants of DToEquilibrium structure R/Aor/Dissociation energy D/evForce constants /auR1=0.962;R2=0.962f1=0.52567,f2=0.53084,f12=-0.0099351105,065fa=0.020371,fan=0.020539,f=0.14863表2Dr(x()和OD,OT(A2∑)的离解能及势能函数的系数Table 2 Dissociation energy, Force constant and Coefficients of 2-body term of DT(X' >')and OD, OT(A221A-1ADT(X∑)3,9233.90673.3980.7414OD(A2∑+)5.56678.5474Or(A2(+)5.64367.11051,01242.5361表3DIO分子三体项分析势能函数的参数Table 3 Coefficients of 3-body term of potential energy function for DhTOC0=-1.87335C1=1,66980C2=1.35304C3=-9.23512C4=-5.6874Cs=3.78958C6=4.14632C9=-0.42493y1=1.75外,由于D+OT和T+OD反应生成DTO基态分1:-123eV2:-11.8eV;3:-10.5eV子存在很深的势阱,所以可以生成稳定的DTO分4:-9.0eV;5-7.5ev6:-6.0ev子,并且通过对势能面逐点进行单点能计算,发现7:4.5eV:8:-3.5eV9:-2.5e该反应过程中不存在势垒因此,说明这是一个无10:-2.0eV;ll:-1.0eV;12:0.5eV阈能的反应.从等值势能图中还可以发现D+OT和T+OD两个反应并不是严格对称的(如等势线8的位置分布可见),这可能是由于氘氚原子的质量不同而引起的同位素效应所有这些合理地体现了DTO的结构特征3分子反应动力学计算方法OD/0.1 nm势能函数确定以后,由准经典的 Monte- carlo图1DIO分子的伸缩等值势能面轨线法131,2求解三粒子保守体系的 HamiltonFig 1 Stretched vibrational contour of potential energy运动方程.设三原子体系为ABC,其质量分别为function for DTOmA,mg,mc,在直角坐标系下,它们的坐标分别为DTO分子的等值势能面如图1,在图1中固定A(q1,q2,q3),B(q4,q5,q6)和C(q7,q8,q).故三粒子体系的 Hamilton运动方程有18个,即∠DOT=105.065°,OD和OT键对称伸缩振动从图中发现,在RDo=Ror=0.96201A附近,出现一个很深的势阱(-12.3eV),也即在平衡点H中国煤化工…,9)(R1=R2=0.96201A),准确地再现了基态DTO在质CNMH、离为质心运动的分子的离解能12306eV,构型是C2v的结构另 Hamilton函数How和相对运动的 Hamilton函数,前者包含六个描述质心运动的变量,后者包1442原子与分子物理学报第25卷含了其余十二个变量.这样描述粒子相对运动的的轨线数,bm为最大碰撞参数,E1为初始相对平Hamilton方程简化为十二个独立变量的微分方动能,V,J分别为振动和转动量子数在计算中N程引入广义相对坐标,经坐标变换,获得十二个相取有限值,其大小由结果所允许误差的范围来确对运动方程,即:定,本文实际计算20000条轨线,bmx值的确定方BnQ=1(=12,…,6),(6)法参照文献23ahne /aP:= Q31D(2S2)+OT(A2∑+)碰撞分子反应动力学(Q1,Q2,Q3)表示粒子C相对粒子B为坐标原点结果与分析的直角坐标,(Q4,Q5,Q6)表示粒子A相对于分对于D(2S)+O(A2∑+)体系的碰撞反应产子BC的质心为坐标原点的直角坐标,(Q7,Q8,物有5种通道Q9)表示整个三粒子体系质心的直角坐标对应于D(2S2)+TO(A2∑+)→Q的共轭动量为P,上式中 Hamilton函数HmO(D2)+DT(X∑)1通道为D(2S)+10(A2∑+)2通道Hd: P2+1 2P2+U(Q1T(2s)+OD)(A2∑+)3通道(10DTO(XA14通道其中OD2)+D(2s2)+T(2S2)5通道11+11=1当D原子的初始平动能为E1=0.4184~627.60ABC B c A, BCkJ/mol时,D(2S2)+OT(A2∑+)反应产物分布及式(7)中,U(Q1,Q2,…,Q6)为分析势能函数经反应截面,见表4从表4可见,随D原子初始相过坐标变换,可以建立内由坐标(RAB,RBC,RcA)对平动能的增加,交换反应1通道的数量逐渐降表示的势能函数的 Hamilton运动方程十二个,粒低,而交换反应3通道的数量随D原子初始相对子相对运动的 Hamilton函数(7)式变为平动能的增加逐渐增加,至E2=209.20 kJ/mol后Hrd=22P2+2u12I随D原子初始相对平动能的增加又逐渐降低;非P2+反应2通道的数量随D原子初始相对平动能的增V(RAB, RBC, RcA),(8)加而逐渐增加;反应4通道(生成络合物的通道)的通过对这十二个方程进行 Runge-Kutta-Gill法和数量很少且随D原子初始相对平动能的增加而迅Adms- Moulton法联合数值求解可以确定出ABC速减小,对此特别计算了一下:当D原子初始相对三粒子体系中碰撞运动的轨线,从而计算出化学反平动能到了E1=20.920kJ/mo时就没有络合物应的反应截面,反应阈能、反应活化能以及反应产生成了;5通道(碰散的通道)在E4>209.20kJ物的分布和各量子态粒子数分布等信息,实现从原mol后随D原子初始相对平动能增加逐渐增加,其子分子水平上研究分子反应动力学反应截面用下主要原因是由于TO的离解能为2.5361eV,只要式计算D原子的初始平动能约为244.68596kJ/mol就有碰散的可能性,因此到E1>209.20kJ/mol时5通a、(E,V,J)=2xP(E,V,,b)如b=道就有反应从整个计算结果看,在低能部分ME,如mN(E,V,0.484-41.8401/m0,D原子把O置换出来相对其他反应较为容易;而在高能部分209.20D)627.60kJ/mol,以反应碰撞和5通道(碰散的通式中N与N分别为计算的总轨线与指定生成物道)为主中国煤化工CNMHG第6期阮文等:DO体系分子反应动力学1443表4D+TO(0,0)产物分布及bm和a与E的关系Table 4 Distributions of products, bmx, nd a, for different channels versus relative translational energy e for D+ TOEra,(1)I channel 2 channel 3 channel 4 channel 5 channelAkJ.mo/nm20.41841505.57300.73480.00790.0034143334.78900.51650.00680.00120.00880.000820.92038540.014541.8401263362510000.02700.0416111900.08360.05980.000040612.98800.01830.02940.0000136344532.87900.00990.01390.0000439.32139342.86200.00950,0104502.081457145342.75500.00720.00340.0000151662.43500.0610.001000(3)0.70(4)0030.0b0000.0020.001-00010501001503003504004505005500.0图3D+TO→T+ OD and DTO反应的a-E2曲线图2D+TO→DI+O反应的a,bmx-E4曲线Fig 3 The reactive cross-section vs Et for D+ To-T+Fig.2 The reactive cross-section and the maximum colli-OD and DTosion parameter vs E for D+ TO--DT+O初始相对平动能E:与反应截面a和最大碰交换反应3通道大体是无热反应,但是有阈能,因撞参数bm的关系曲线如图2,3.在所计算的能量为截面a随平动能E先增加后减小,只不过阈能范围内D(2S)+O(A2∑+)的反应碰撞主要表值很小而已当平动能增加到334.72 kJmol时现为交换反应1通道和非反应碰撞2通道交换反由于5通道轨线数增加,导致3通道轨线数又开始应(1通道)是强烈无阈能放热反应,在0.4184-降低20.920kJ/mdl低能量范围伴随有络合物生成,由3.2r(s2)+DO(A2∑+)碰撞的分子反应动力学于生成络合物的反应不满足动力学轨道受陷条件结果与分析所以低能时虽然有极少量络合物,但是其反应截面+∩TA2+)体系的碰撞反应产远小于1通道这说明D()+Or(A2)(O物也中国煤化工(Dg)+DT(X(∑g)交换反应很容易发生.而CNMHG1444第25卷T(2S2)+DO(A2∑+)对平动能的变化,该体系的反应特点与D(2s)+O(D2)+D(X∑g)1通道OT(A2∑+)体系的反应特点类似,但各对应初始(2s2)+DO(A2∑+)2通道相对平动能下的最大碰撞参数bmx比较D(2S2)+D(2s)+OT(A2∑+)3通道(11)OmA2∑+)有所增大,1和3通道的反应截面aDTO(XAJ)4通道也较D(2s2)+Or(A2∑+)有所增大,相反4通道O(D2)+D(2S2)+T(2S)5通道的反应截面σ,有所下降,反应4通道(生成络合物当T原子的初始平动能为E=0,4184~627.60的通道)的数量也较D(2S)+Or(A2∑+)更少,但kJ/mol时,T(s2)+DO(A2+)反应产物分布及是出现络合物的初始相对平动能增加到了E=反应截面a见表5.从表5可见,随T原子初始相表5T+DO(0,0)产物分布及bm和a与E的关系Table 5 Distributions of products, bmax and o, for different channels versus relative translational energy Et for T+DOI channel 2 channel 3 channel 4 channel 5 channel0.41845.63300.75020.01270.00124.1840146644.84900.54170.01280.0006144054.60100.47910.01390.000320.920137634.26200.39280.02230.0000314540.31640,000033.83300.26710.0490.000083,6801064266203.71500.23070.0593209.2065753.31000.1132334.7223983.10200.03630.07451577376850962.99600.02220.0531439.3215430.02160.0479502,080.01830.0341627,6001111582.82500.01630.01430.自Eb0.002-0.001500图4T+DO→·DT+O的a,bm-E2曲线图5T+DO→D+TO和DTO的a1-E2曲线Fig. 4 The reactive cross-section and the maximum colli中国煤化工EfrT+DO→D+sion parameter vs Et for T+ Do-DT+OTTHCNMHG阮文等:DO体系分子反应动力学1445初始相对平动能E1对反应截面a,和最大碰(2S)+TO(A2∑+)和T(2s2)+DO(A22+)两个撞参数b。的关系曲线如图4,5在所计算的能量反应体系的反应特点相似;但是在相同碰撞能量范围内,交换反应(1通道)是强烈无阈能放热反时,最大碰撞参数以及各对应反应通道的碰撞截面应在0.4184~41.840kJ/mol低能量范围伴随有等存在一定的差异.可以设想若体系为普通水分子极少量络合物生成,由于生成络合物的反应也不满(HO)的反应碰撞,这种差异应该是不存在的,本足动力学轨道受陷条件,所以低能时虽然有极少量文出现的这种差异以及前面势能面上的轻微不对络合物的反应发生,但是其反应截面远小于1通称现象可能是由于氘氚原子的质量不同而引起的道这说明T(2Sx)+DO(A2+)→O(Dg)+DT同位素效应(x∑g)交换反应容易得多交换反应3通道也是参考文献有阈能的无热反应总体反应特点与D(2s2)+TO(A2+)的反应碰撞特点相似可以设想如果是普10uM,LoYM,Hmn,aa., Experimental study on通水H2O,两个体系反应的产物轨线分布、各对应Hydrogen isotopic Deuterium and Tritium from gas通道的反应截面σ、最大碰撞参数bm等与能量phase to liquid phase by Hydrogen-water exchange [J]的对应关系应该是完全相同的,但是通过比较两个At.Ener.Sci. Tech, 2005, 41: 425(in Chinese古梅,罗阳明,韩军,等.氢同位素D、T从气相到体系的轨线分布及各通道的反应截面,不难发现两液相的氢H水交换实验研究[J].原子能科学技术,个体系的反应动力学存在诸多不一致的情形,这应2007,41:425]该是氘氚原子的质量不同而引起的同位素效应,这[2]FuZH, LUOY m, WanghU,eal,. A pilot test on一点有待进一步研究performance evaluation of Pt-SDB catalyst for hD(g)HO(1)isotopic exchange reaction[J]. NucL. Tech4结果与讨论2004,27:448( in Chinese)[傅中华,罗阳明,王和用半经验轨线法考察了D(2s2)+TO(A义,等. Pt-SDB疏水催化剂应用于氢?水同位素交∑+)和T(S2)+DO(A2∑+)两体系的分子动力换的实验研究[J].核技术,2004,27:448]学在计算结果的基础上,比较两个反应体系的动[3] Zhong Z K, Zhang L,SuY,eral. Calculation of theseparation factors for hydr力学特点这对了解重水的分子结构以及重水提氚and water isotopiechange reaction [J]. At. Ener. Sci. Tech, 2004的研究具有一定的指导意义38:148( in chinese)[钟正坤,张莉,孙颖,等.氢水(1)低能碰撞时,D(2s3)+TO(A2∑+)和T同位素交换分离因子理论计算[J.原子能科学技(2s)+DO(A2∑+)两体系的碰撞,反应主要为第术,2004,38:1481通道为主,即:D2s)+T0(A2+)→O(Dg)+[4】 Cheng B M, Chung CY, Mohammed Bahou,aDT(X12)和T(S2)+DO(A2∑+)→O(Dg)+Quantitative spectroscopic and theoretical study of theDT(X∑g),同时伴随少量产物第3通道产物高otical absorption spectra of H2O, HOD, and D,O inthe 125 145 nm region[J].J. Chem. Phys., 2004,能碰撞时,主要为非反应碰撞和第5通道产物120:2242)计算结果表明:反应体系D(2s)+TO(A2[s] Attila g,Caor, abor czako, Tibor furtenbach,∑+)和T(2S2)+DO(A2∑+)由于均不满足动力et al. On equilibrium structures of the water molecule学轨道受陷条件,故仅在低能条件下有极少量的络[J].J.Chem.Phy,2005,122:214305合物生成并且从计算结果还可知:D(s)+Or6 Wang WL, Emanuel santos, odo branda.Thee(A2∑+)和T(2S)+DO(A2∑+)反应生成DTOcal rate coefficients for the exchange reaction OH+D\的反应截面非常小,基本不发生反应;而D与TOOD+H[J].J.Chem,Phys.,2006,124:074305或者T与DO反应生成O+DT要容易得多,这是[7] Paolo barletta, Sergei V,Shin, Nikolai Fal. CvRQd ab initio ground-state adiabatic因为OT(A2+)和OD(A2∑+)的离解能分别为2.5361eV和2.5360eV,D原子接近TO或T原子接近DO时,比较容易离解TO或DO分子生成[8]dYHale[j].J.中国煤化工CNMHG. Analytical poten-DT分子ial energy function for tritium water molecule T,O(X(3)从表4,5以及图2-5可以发现,尽管DA1)[J. Acta Phys.Sin,,2008,57(8):4833(in1446原子与分子物理学报第25卷Chinese)[罗文浪,阮文,张莉,朱正和氢同位素氚激光振转能级粒子数分布的影响[].原子分子物水T2O(XA1)的解析势能函数[J].物理学报,理学报,2000,17:1442008,57(8):4833][16] Ran M, Huang P. The effect of N2 for vibration-ro[9] Luo D L, Sun Y, Liu X Y, et al. Structure and poten-tional populations of HF overtone laser[J].J,Attial energy function investigation on UH and UH2 moM.Phw.,2002,19:487( in Chinese)[冉鸣,黄claes[J]. Acta Phys. Sin., 2001, 50: 1896(in Chi萍.N2对H泛频激光振转能级粒子数分布的影响nese)[罗德礼,孙颖,刘晓亚,等,UH和UH2分子J.原子与分子物理学报,2002,19:487]构与势能函数[J.物理学报,2001,50:1896117] Chen Z Q,XeHP. The paraliel complex reactions[10]Xue WD, Wang H Y, Zhu ZH, et aL. The structuredynamics of Si(3P4)+O(X∑”)[J.J.Aand potential energy function of CUO[J]. Acta PhysM.Phy,2006,23:579( in Chinese)[陈泽琴,谢Sin.,2002,51:2480( in Chinese)[薛卫东,王红洪平.S(P)+OO(X∑+)体系的平行络合竞争动艳,朱正和,等,CUO分子结构与势能函数[.物力学[J.原子与分子物理学报理学报,2002,51:2480][11] Li Q, Lu H. Study on Molecular Reaction Dynami[18] Chen Z Q. Study on the reactions dynamics for synfor Pu(Fg)+H2(X2*, 0, 0)[]. Acta chemthesizing SiC[J].J. At. Mol. Phys.,2007,241229( (in Chinese)陈泽琴.碳化硅合成反应动力学Sin.,2003,11:1881( in Chinese)[李权,卢红.P研究[].原子与分子物理学报,2007,24:1229](Fg)+H(X1∑,0,0)的分子反应动力学[.[19] Sorbie K s, Murrell N. Analytical potentials化学学报,2003,11:1881[12 Luo D L, Meng D Q, Zhu Z H. Potential energyatomic molecules from spectroscopic data [J1975functions and molecular reaction dynamic s for LiHLiO and LiOH[J]. Acta Phys.Sin,2003,52[20] Moore C E. Atomic Energy Levels [M]. WashingtonDC: U.S. Government Printing Office, 19712438( in Chinese)[罗德礼蒙大桥,朱正和.LH,[21] Xue WD, Zhu ZH,RanM,etal. Theoretical studyLO和LOH的分析势能函数与分子反应动力学[J.物理学报,2003,52:2438]on molecular reaction dynamics of the U+ CO system[13] Liu X Y, Zhu Z H, Jiang G, et aL. The reaction dy[JI. Acta Phys. Sin., 2002, 51: 2503( in Chinese[薛卫东,朱正和,冉鸣,等.U+OO体系的分子反namics of U and Hydrogen[J]. Sci. Chin, Series G应动力学研究[J.物理学报,2002,51:2503]2005,35:415( in Chinese)[刘晓亚,朱正和,蒋刚等.铀与氢气的反应动力学[J.中国科学:G辑,[22] Yu s Q. Micro-chemical reaction[M].Hefei: anhui2005,35:415]Science and Technology Press, 1983(in Chinese)lAr书勤.微观化学反应[M].合肥:安徽科学技术出[14] Zhu ZH, Yu H G. In Molecular Structure and po.版社,1983tential Energy Function[M]. Science Press, Beijin1997( in Chinese)[朱正和,俞华根.分子结构与分[23] Gao T, Wang H Y, Zhu ZH子势能函数[M].北京:科学出版社,199y]study on molecular reaction dynamics of Pu[15] Liu X Y, Jiang G, Wang FH, et al. The effect of HOn(xS2)System [J]. J. At. Mol. Physfor vibration-rotational populations of HF overtone1999,16:162( in Chinese)[高涛,王红艳,朱正和,laser [J]. J. At. Mol. Phys., 2000, 17: 144(in等.PLO2体系的分子反应动力学研究[J].原子与Chinese)[刘晓亚,蒋刚,王藩侯,等.H对HF泛频分子物理学报,199,16:162]中国煤化工CNMHG