第17卷第5期化學物理學報Vo.17,No.52004年10月CHINESE JOURNAL OF CHEMICAL PHYSICS20041003-7713/2004/05-543-4堿金屬團簇碰撞動力學研究徐秀蓮*,陳紅霞,周朋霞",曾祥華a.揚州大學物理科學與技術(shù)學院,復(fù)雜性科學中心,揚州225002b,鹽城師范學院物理系,鹽城224002)摘要:采用距離相關(guān)緊密束縛的分子動力學模型,研究了鈉原子團簇的碰撞動力學行為.研究得出,對Na3D)十Nas系統(tǒng),碰撞參數(shù)為b=Ωa。時,對應(yīng)于深度非彈性碰撞過程(DC),在b=13a。時,則對應(yīng)于準彈性碰撞過程.另外,得到碰撞過程中存在轉(zhuǎn)動,其能量大小與碰撞能量E-/n以及碰撞參數(shù)b均有關(guān)系.隨著E/n的增大,相對轉(zhuǎn)動減小;隨著b的變化,轉(zhuǎn)動能量先增大后減小.在b=7a。時轉(zhuǎn)動最大,最大轉(zhuǎn)動能量可達總能量的十分之關(guān)鍵詞:鈉原子團簇;緊束縛分子動力學模型;相對轉(zhuǎn)動;深度非彈性碰撞中圖分類號:O456文獻標識碼:AStudy on the dynamics of Sodium Cluster CollisionsXu Xiulian", Chen Hongxia.b, Zhou Pengxia, Zeng Xianghua"(a. College of Physics Science and Technology, Com plexity Science Center of Yangzhou University, Yangzhou 225002b. De partment of Physics, Yancheng Teachers College, Yancheng 224002)Abstract Within the framework of distance dependent tight-binding molecular dynamics DDtB-MD), the collision dynamics of sodium cluster Na, has been studied systematically. Some phenomena havebeen observed at different impact parameters b, such as fusion reaction, deep inelastic collision(DIC)andquasi-elastic collision, which are similar to the nuclear heavy-ion collisions (HIC). For the system of nap(3D)+Nas, the reaction mechanism at b=9ao is DIC, but when b is equal to 13ap. it corresponds to quasi-elastic collision. Furthermore the rotation processes during the collisions, are related to the collision energy and parameter. The larger collision energy is, the earlier relative rotation will occur, and the relaxationtime becomes shorter and the relative rotation energy is much smaller. There exists maximal relative rotation energy which corresponding to b=7ao. When b is smaller than 7ao, the rotation energy increases withb increasing, otherwise the energy decreases. And the maximal relative rotation energy is corresponding toDIC process. The maximal rotation energy can reach one tenth of total energy, which is much less thanthKeywords Nan, DDTB-MD, Relative rotation, DIC1引言由中子和質(zhì)子組成的原子核體系與由價電子和離子組成的堿金屬團簇體系很類似,它們都是典型的多金屬團簇因為具有許多新奇特異的物理、化學體費體但由工閉能休系中電磁相互作用的性質(zhì),一直是人們關(guān)注的熱點,最近,人們發(fā)現(xiàn)長程中國煤化為有別于原子和原子堿金屬團簇表現(xiàn)出典型的集體運動模式,這是因為核體、CNMHG金屬團簇碰撞表現(xiàn)出Supported by the National Natural Science Foundation of China (60072032)and Foundation of the Education Commission of芳數(shù)據(jù) of china(031B140152* Corresponding author, E-mail: xuxiulian@ sina, com Received 19 June 2003; in final form 5 December 2003544化學物理學報第17卷來了與重離子碰撞不同的現(xiàn)象.原子團簇碰撞的和湮滅算符;矩陣元h為通過用3s軌道做基并包研究為人們認識原子團簇的動力學性質(zhì)提供了可括3p軌道微擾的距離相關(guān)的量;p(R)、t(R)和能2.目前,包括了原子核自由度的典型理論方法有。(R)分別對應(yīng)原子間距離為R時的離子離子斥兩種:一種是量子化學類型的計算,如采用局域密力、ss和spn的轉(zhuǎn)換積分,所用參數(shù)見文獻[12]度近似( Local Density Approximation)LDA的Car在 DDTB- MD計算中,對初態(tài)團簇中的每一個Parinello技術(shù),但只能用于處理很小數(shù)目的原子原子,坐標空間的位置和其具有最低能量的基態(tài)位團簇;第二種如考慮了原子核自由度的 Jellium模置相同初始時,入射彈團簇和碰撞靶團簇相距大約型和經(jīng)驗的 Nilsson勢模型等,往往需預(yù)先考慮40a0(此處a0=0.529A,是Bohr半徑),使兩者之些對稱性近年來,分子束、團簇束技術(shù)的發(fā)展已經(jīng)間無相互作用在動量空間,初始時給團簇一個的熱能夠進行原子團簇之間的碰撞,于是迫切需要理論激發(fā),原子隨時間演化滿足牛頓方程:方法對團簇碰撞中岀現(xiàn)的現(xiàn)象給予解釋。.在d 2Born- Oppenheimer近似下,人們發(fā)展了一個距離相VER關(guān)的緊束縛分子動力學模型( Tight-binding Molec-其中E,為 Born-Oppenheimer勢能面;m;和r分ular Dynamics) DDTB-MD6.它是一個簡化的別為第i原子的質(zhì)量和坐標.作用到每個原子上的Born-Oppenheimer勢能面上的從頭計算法(abin力由 Hellmann-Feynman定理給出,即tio)MD.采用 DDTB-MD系統(tǒng)研究鈉原子團簇Nan的結(jié)構(gòu)、熱力學和碰撞動力學性質(zhì)0-151,計算結(jié)果在碰撞過程中,時間步長為1fs,熱化時間取能夠很好地符合 ab initio的計算結(jié)果和實驗結(jié)果為2ps體系的總能量E、勢能E和動能E關(guān)系在以往的團簇研究中0:,我們得到大分子團簇的形成與初始團簇的幾何構(gòu)型有關(guān),與入射道參數(shù)(能量、不同的體系以及不同的碰撞參數(shù)等)相關(guān)EL=ER+E, Ek并且在碰撞過程中也觀察到了與重離子碰撞相似的n1、n2分別為彈、靶團簇的原子數(shù)目動能E4由相反應(yīng)機制.本工作將繼續(xù)研究團簇碰撞過程的動力對動能E和熱能E1組成:學的行為,通過跟蹤各體系碰撞過程中的空間結(jié)構(gòu),Er= Ere+e討論碰撞過程中轉(zhuǎn)動的行為,以及由轉(zhuǎn)動引起的能erl=1量耗散現(xiàn)象.這方面深入的討論并不多見.對NaPi)存在平面結(jié)構(gòu)的Dh(2D),勢能為-3.311eV,和空定義團簇之間相對轉(zhuǎn)動能量為間結(jié)構(gòu)的C(3D),勢能為-3.349eV.我們討論的是結(jié)構(gòu)較為穩(wěn)定的Na6(3D)、Na兩團簇之間的碰Em=Em彈+E吧=21+2l撞這里J彈、J粑、Ⅰ彈、l粑分別是彈、靶團簇的角動量和轉(zhuǎn)動慣量,它們分別定義為2模型J×J靶在計算團簇的動力學過程中,要確定團簇的Born-Oppenheimer勢能面,它是通過距離相關(guān)的緊mr:, I束縛的單電子哈密頓量本征值的和來確定的,即3結(jié)果與分析3.1反應(yīng)機制對碰撞參數(shù)的依賴關(guān)系hi= h o)+h(2)中國煤化工原子碰撞能量En/n0(Ra)-t2(R)CNMH(果表明,團簇間發(fā)生相0+h2互作用的時間為2.7ps,相互作用使彈靶團簇體系-t,(R,))t(r)t(RikRR發(fā)生形變,該形變能使團簇發(fā)生轉(zhuǎn)動,體系激發(fā)能增ERr大.在轉(zhuǎn)動過程中不同的原子所獲得的能量不同,部其中,a和種教癆別為s軌道上i基和j基的產(chǎn)生分轉(zhuǎn)動能量高的原子處于較高的激發(fā)態(tài)很容易被發(fā)第5期徐秀蓮等:堿金屬團簇碰撞動力學研究545射,使體系的能量降低并被束縛在團簇的位壘中,成隨時間的變化關(guān)系圖.在此,我們用E。/E來描述為一個穩(wěn)定的復(fù)合體系.這種將初始動能轉(zhuǎn)變成復(fù)相對轉(zhuǎn)動.圖2b描述的是該體系在Epo/n=1eV合體系的內(nèi)部激發(fā)能,高激發(fā)的復(fù)合體系通過發(fā)射b分別為0a。、9a、13a0時轉(zhuǎn)動能量E隨時間的單聚體(或二聚體)而衰變的能量耗散過程,與重離變化.從圖2可以看出轉(zhuǎn)動主要發(fā)生在碰撞階段,碰子碰撞中的DC反應(yīng)機制相似.圖1a是相對動能、撞前及平衡后轉(zhuǎn)動均很小幾乎可以忽略.且Es/n動能以及熱能隨時間的變化關(guān)系.可以看岀初始的越大轉(zhuǎn)動發(fā)生得越早,但持續(xù)的時間以及相對轉(zhuǎn)動集體能量大部分轉(zhuǎn)化為系統(tǒng)的內(nèi)部激發(fā)能,這種過卻越小.由此,我們可以知道轉(zhuǎn)動主要是由彈、靶團程說明了DlC的存在.當b=13a。時,由于碰撞參數(shù)簇間的相互作用引起的,碰撞過程中彈、靶團簇靠得增大,彈、靶團簇中的原子幾乎沒有參與直接相互作較近,相互作用較大,團簇的各原子所受作用力不均用,而且由于初始的能量較大,這種相互作用只改變等,且彈、靶團簇所受作用力不在一直線上,從而產(chǎn)了團簇的形狀,隨著時間的弛豫各自成為獨立的團簇,碰撞為準彈性碰撞.如圖1b所示,熱激發(fā)能改變rc=2 ev/n很小,團簇主要以初始的集體速度運動pro =1.5 eV02至23000Time/fs4000Time/fs0250200.150.1000060006o008000Time/fs圖2a.Em/Epm隨時間的變化;隨時間的變化圖1動能、相對動能、熱能隨時間的變化Fig. 1 Changes of Et, Erel, Eth with the timeb. Changes of Erot with the time生了力矩,最終引起了轉(zhuǎn)動.碰撞前、碰撞后,彈、靶3.2團簇碰撞過程中的轉(zhuǎn)動團簇相距較遠,勢場作用很小,幾乎不會引起轉(zhuǎn)動.在碰撞的初始時刻,我們給靶團簇隨機地扭轉(zhuǎn)個角度以保證其碰撞的隨機性,同時去除彈、靶團E大提地早可以解釋為Emn/n越簇的轉(zhuǎn)動能量,這樣保證彈、靶團簇在沒有轉(zhuǎn)動的情大時彈、靶團簇相互靠近得越?jīng)r下碰撞,但是隨著碰撞的進一步深入,彈、靶團簇CNMHG時相對運動速度越快,彈、靶團簇靠近過后背離得也越快,所以持續(xù)的時間間發(fā)生了相對轉(zhuǎn)動,該轉(zhuǎn)動在深度非彈性碰撞中尤反而短至于E/n越大相對轉(zhuǎn)動卻越小,這是因為突出.圖2a是碰撞體系Na(3D)+Nas在b=為,E/n越大,越過勢壘的能力越強,勢場的相對9a,E萬5、1.0cV時相對轉(zhuǎn)動能影響就越小,勢場引起的相對轉(zhuǎn)動也就越小546化學物理學報第17卷從圖2b我們看出,b對轉(zhuǎn)動能量E的大小也統(tǒng),碰撞參數(shù)b為Ωaσ對應(yīng)DIC(深度非彈性碰撞),有很大的影響為進一步討論E對b的依賴關(guān)系,為13a。則對應(yīng)準彈性碰撞;在碰撞過程會產(chǎn)生了較圖3給出了碰撞體系Na(3D)十Na8在不同碰撞參大的轉(zhuǎn)動,該轉(zhuǎn)動與碰撞能量E/n以及碰撞參數(shù)數(shù)b下最大轉(zhuǎn)動能量的結(jié)果.中心碰撞時,彈、靶體b有關(guān).隨著E。/n的增大,相對轉(zhuǎn)動發(fā)生得越早但系的相互作用較強,初始運動的能量大部分轉(zhuǎn)化為是持續(xù)時間以及相對轉(zhuǎn)動能量越來越小.隨著b的熱能,轉(zhuǎn)動很小.隨著b的增大,參與直接相互作用變化,轉(zhuǎn)動能量先增大后又減小,在b=7a。的DIC的Na原子數(shù)減小,彈、靶團簇中原子受力不均勻,中轉(zhuǎn)動能量達到最大值,稱b=7a。為極值碰撞參參與作用的區(qū)域產(chǎn)生很大的集體流,為保持體系的數(shù).此時,轉(zhuǎn)動能量大約為總能量的十分之一,這要平衡,出現(xiàn)較大的轉(zhuǎn)動.在b=7a。時,轉(zhuǎn)動達最大,比HIC過程中的轉(zhuǎn)動小得多最大轉(zhuǎn)動能可達總動能的十分之一.而隨著的繼續(xù)參考文獻增大,兩者之間的重疊部分逐漸變小,參與作用的粒子數(shù)減少,相互作用時間變短且作用強度減弱,集體1] Zhang Wenxian(張文獻), Liu lei(劉磊), Li Yufen(李流也減小,彈、靶主要以各自的初始速度運動,導致郁芬),etal.Chin.J.Atom.Mol.Phys.(原子分子物理轉(zhuǎn)動越來越弱.由此得到,轉(zhuǎn)動集體流從近中心到周學報),1999,4:592邊碰撞先増大后減小,最后在非常周邊碰撞處消失[2] Zhang Fengshou(張豐收), Ge lingxiao(葛凌霄也就是說,中心碰撞時相對轉(zhuǎn)動能量并不是最大;存Nuclear Physics Review, 1998, 15: 38在一個極值碰撞參數(shù)(bn=7a);小于bn時,隨著b3] Zhang Fengshou(張豐收), Ge lingxiao(葛凌霄).Nuclear Multi-fragmentation(原子核多重碎裂), Beijing增加,相對轉(zhuǎn)動變強,轉(zhuǎn)動能量也增大;大于b時,(北京): Science press(科學出版社),1998.10轉(zhuǎn)動能量隨b的增加迅速地減小,到最后幾乎沒有41 Hortig g, Muller M.z.Phy、,1969,21:11相對轉(zhuǎn)動,因為在b很大時,彈、靶團簇間的相互作5 Li Quan(李權(quán)), Lu hong(盧紅), Wang hongyan(王用已經(jīng)很微弱了紅艷),etal.Chin.J.Chem.Phys.(化學物理學報)2003,16:996 Kroto H W, Heath JR, O'Brien SC, Curl R F, Small04[7] Iijima S. Nature,1991,354:56[ 8]Tenne R, Margulis L, Genut M, Hodes G. Nature0.061992,360:444[9] Margulis L, Salitra G, Tenne R, Talianker M. Nature1993,365:113[10] Xu Xiulian(徐秀蓮), Zeng xianghua(曾祥華),WangFeng(王鋒),etal. Acta Phys.Sin.(物理學報)24681012142002,51:31[1] Zeng xianghua(曾祥華), Xu Xiulian(徐秀蓮),WangFeng(王鋒),etal. Acta Phys.-Chim.Sin.(物理化學圖3(a)、(b)圖分別是最大轉(zhuǎn)動能量以及最大轉(zhuǎn)動能量學報),2002,18:2與對應(yīng)時刻總動能的比值隨b的變化[12 Zhang FS, Suraud E, Spiegelman F, et al. Z. Phys.Fig 3 The changes with the collision parameter b1995,D35:131(a)Maximal rotational energy; (b) The ratio of[13] Zhang F S, Spiegelman F, Suraud E, et al. Physmaximal rotational energy to total kinetic energy中國煤化工[14]F. et al. ch4結(jié)論CNMHG[15] Zhang Fengshou(張豐收), Ge lingxiao(葛凌霄)采用距離相關(guān)緊密束縛的分子動力學模型,研Chin.J.Aom.Mol.Phys.(原子分子物理學報),究鈉原子團簇的碰撞動力學一些性質(zhì).我們發(fā)現(xiàn):碰1997,14:583撞反應(yīng)機隨參數(shù)有關(guān)對Na(3D)+Na系