樊群超孫衛(wèi)國?李會東馮灝
1)(西華大學(xué)物理與化學(xué)學(xué)院,先進(jìn)計算研究中心,成都610039)
2)(四川大學(xué)原子和分子物理研究所,成都610065)
(2009年11月23日收到;2010年1月30日收到修改稿)
Li2部分電子態(tài)的完全振動能譜與離解能的精確研究*
樊群超1)孫衛(wèi)國1)2)?李會東2)馮灝1)
1)(西華大學(xué)物理與化學(xué)學(xué)院,先進(jìn)計算研究中心,成都610039)
2)(四川大學(xué)原子和分子物理研究所,成都610065)
(2009年11月23日收到;2010年1月30日收到修改稿)
本文將孫衛(wèi)國等建立的精確計算雙原子分子離解能的解析表達(dá)式作為分子振動能級正確收斂的重要物理判據(jù),應(yīng)用代數(shù)方法,進(jìn)一步研究了重要的堿金屬Li2的5個電子態(tài)的完全振動能譜和離解能,獲得了這些電子態(tài)的精確振動光譜常數(shù),同時也得到了包括接近分子離解極限在內(nèi)的完全振動能譜和與實驗值符合得很好的理論離解能.這些計算結(jié)果為許多需要這些電子態(tài)振動數(shù)據(jù)的研究領(lǐng)域提供了重要的研究數(shù)據(jù).
代數(shù)方法,Li2,振動能級,離解能
PACC:3310,3520P
堿金屬Li2由于其簡單的電子結(jié)構(gòu),成為實驗和理論上用來研究中性同核雙原子分子之間相互作用的典型分子.幾十年來,國際上對Li2的分子結(jié)構(gòu)、勢能函數(shù)、振動-轉(zhuǎn)動光譜性質(zhì)以及各電子態(tài)的離解能進(jìn)行了大量的理論和實驗研究[1—9],并不斷取得新的進(jìn)展.理論上,人們對Li2的振轉(zhuǎn)能級以及部分電子態(tài)的研究主要基于從頭計算方法[1—4],直到1999年,Linton等[4]才首次應(yīng)用近離解極限展開技術(shù)(near dissociation energy expansion,簡稱NDE)對Li2的a3Σ+u電子態(tài)的低激發(fā)振動譜項值進(jìn)行了研究,并確定了該電子態(tài)的離解能.建立這些理論方法,往往需要考慮體系的多組態(tài)函數(shù)以及電子相關(guān)效應(yīng),根據(jù)經(jīng)驗選擇基函數(shù)的集合以及很多的變分參數(shù),然后經(jīng)過復(fù)雜的計算和合理的修正,最后才能夠得到收斂的計算結(jié)果.然而,要獲得體系正確的基函數(shù)集合和變分參數(shù)卻是非常困難.由于這些理論方法存在其自身的局限性,因此對于大多數(shù)電子態(tài),最終很難從理論上獲得高激發(fā)振動光譜和離解能的正確數(shù)據(jù).
實驗上,對Li2電子態(tài)振動能級和離解能的研究方法主要有微擾光學(xué)雙共振光譜法(perturbation facilitatedopticalopticaldoubleresonance spectroscopy,簡稱PFOODR)[10,11]和激光誘導(dǎo)熒光傅里葉變換光譜法(laser induced fluorescence fourier transform spectrometry,簡稱LIFFTS)[12].然而,對于大多數(shù)電子態(tài)而言,當(dāng)分子處于高激發(fā)振動態(tài)尤其是處于分子離解極限附近的時候,由于分子振動能級的間隔非常小而且密,從而對實驗儀器的精度要求很高,使得很難獲得此刻的正確振動能級.因此在實驗上利用這些方法一般只能觀測到分子處于低激發(fā)態(tài)的一些振動能級,而所獲得的分子離解能大多是通過擬合勢能曲線、外推法來得到的數(shù)據(jù)[7].目前,已發(fā)表的文獻(xiàn)中對于有關(guān)Li2長程分子作用勢、離解能和精確的高激發(fā)態(tài)振動能級的數(shù)據(jù)還比較缺乏.因此,我們?nèi)匀缓苡斜匾獙i2各電子態(tài)的完全振動能譜及離解能進(jìn)行進(jìn)一步的研究.
本文應(yīng)用孫衛(wèi)國等[13—16]建立的精確計算分子完全振動能譜代數(shù)方法(AM),從精確的實驗?zāi)芗壔蛄孔永碚撜駝幽芗壋霭l(fā),通過嚴(yán)格求解能級與光譜常數(shù)的正確代數(shù)關(guān)系式,避開了繁瑣的勢能函數(shù)和薛定諤方程求解過程,獲得了7Li2的a3Σu+,33Σg+,13Δg,b3Πu和6Li2的11Πg等5個電子態(tài)的精確振動光譜常數(shù)和完全振動能譜,同時應(yīng)用孫衛(wèi)國等[17,18]最近建立的計算精確的分子離解能的新解析表達(dá)式,計算了以上5個電子態(tài)的理論離解能DAeM,獲得了十分精確的結(jié)果.
在前一階段的工作中,我們總結(jié)了LeRoy和Bernstein[19,20]研究分子離解能的理論工作,在他們的工作基礎(chǔ)上,對于一個穩(wěn)定的雙原子分子或離子電子態(tài),我們獲得了計算這些電子態(tài)離解能的新公式[17,18]
若(2)式不滿足,則說明用于計算離解能Dceal的三個最高振動能級含有不可忽略的誤差.
為了獲得雙原子分子精確的完全振動能譜,孫衛(wèi)國等[13—16]用二階微擾理論得到了非相對論核運(yùn)動振動能級的解析展開表達(dá)式,即
AM[13—16]的基本出發(fā)點是將(3)式改寫成矩陣形式
其中振動光譜常數(shù)的向量矩陣X和振動能量矩陣E分別表示為
AM方法基于這樣的物理事實:對于絕大多數(shù)雙原子分子的穩(wěn)定電子態(tài),現(xiàn)代實驗技術(shù)和量子理論往往難以獲得高振動激發(fā)態(tài)尤其是接近離解極限區(qū)域的振動能級,總可以獲得各電子態(tài)的一部分量子態(tài)不太高的振動能級的精確數(shù)據(jù)子集合[Eν].而這些精確的實驗數(shù)據(jù)基本上包含了所有重要分子振動信息和所有重要的微觀量子效應(yīng),如電子-電子相關(guān)、自旋-軌道耦合、電磁效應(yīng)等等.因此,用AM方法在數(shù)學(xué)和物理上不加近似地從這些精確實驗?zāi)芗壸蛹陷腿〕鲆唤M包含了所有重要的高階非諧性效應(yīng)和分子振動信息的振動光譜常數(shù)X,再將這組常數(shù)代入振動能級的正確表達(dá)式(3)式,即可得到某電子態(tài)的真實振動完全能譜的一個正確表象{Eν}.所以,基于有限的精確實驗數(shù)據(jù)和嚴(yán)格的微擾理論表達(dá)式(3)所獲得的這個振動能級的完全集合{Eν}不僅能精確地重復(fù)已知實驗?zāi)芗壸蛹希跡ν],而且能正確地獲得該電子態(tài)所有高激發(fā)振動量子態(tài)的能級.而這些能級往往是實驗上很難得到的.應(yīng)用AM方法于雙原子分子的電子態(tài)時,可以從由m個已知實驗?zāi)芗壗M成的能級子集合[Eν]中選取N個小能級組,每組8個能級,然后解方程(3)N次,并獲得N組振動光譜常數(shù)X′s.這N組常數(shù)中,總有一組常數(shù)X能最好地滿足下列要求:
那么這組振動光譜常數(shù)就是該分子體系真實振動光譜常數(shù)集合的最佳物理表象之一,從而由此計算的振動能譜就是包含該體系所有真實振動能級的一組完全振動能譜{Eν}.(11)式即計算離解能的新公式(1)式,由于該式要求使用三個精確的最高振動能級,所以(1)式又可用作雙原子體系的完全振動能譜是否很好收斂的一個重要物理判據(jù).
本文應(yīng)用AM方法分別計算了7Li2的a3Σu+,33Σg+,13Δg,b3Πu和6Li2的11Πg等5個電子態(tài)的精確振動光譜常數(shù)和完全振動能譜,再根據(jù)這些電子態(tài)的正確物理能級,利用(1)式分別獲得了它們的理論離解能Dcale.表1列出了上述5個電子態(tài)的分別從文獻(xiàn)中獲得的和AM方法得到的振動光譜常數(shù)(ω0,ωe0,ωe,ωexe,ωeye,ωeze,ωete,ωese,ωere).這些電子態(tài)的實驗振動能譜和AM振動能譜以及離解能的數(shù)據(jù)分別列于表2中.
表1很清晰地顯示出,由AM方法得到的諧振常數(shù)ωe與相應(yīng)電子態(tài)的文獻(xiàn)值符合得非常好.就多數(shù)電子態(tài)而言,由AM方法得到的非諧振常數(shù)ωexe與文獻(xiàn)值符合得也比較好;而文獻(xiàn)中給出的實驗振動光譜常數(shù)還比較缺乏,尤其是高階振動光譜常數(shù)的數(shù)據(jù),其主要原因之一是文獻(xiàn)在使用類似方程(3)的振動能級展開表達(dá)式中截取的項比較少.例如,實驗上在對7Li2的33Σ+[23]電子態(tài)進(jìn)行研究時只取到了3階,所以其誤差較大,因此如果用僅有的這些低階振動光譜常數(shù)去計算這些電子態(tài)的全部振動能量,則得到的完全振動能譜尤其是高激發(fā)振動態(tài)的能譜數(shù)據(jù)將含有很大的誤差,而由AM振動光譜常數(shù)計算得到的最高振動能量與實驗離解能Dexptge非常接近,這表明由AM方法得到的振動光譜常數(shù)比文獻(xiàn)值精確得多.
表1 用AM方法和從文獻(xiàn)中分別得到的Li2部分電子態(tài)的振動光譜常數(shù)(單位:cm-1)
將表2所列的各電子態(tài)的AM能譜與文獻(xiàn)給出的實驗?zāi)芗壪啾容^可知,AM能譜的精確度完全取決于實驗?zāi)芗壍木_度,只要實驗?zāi)芗壘_可靠,則AM將產(chǎn)生一組十分準(zhǔn)確的光譜常數(shù),由此計算而得到的AM振動能譜將必然很好地再現(xiàn)實驗?zāi)茏V并正確地產(chǎn)生包括最高振動能級在內(nèi)的全部高階振動能級,從而獲得離解能的準(zhǔn)確度也應(yīng)該是相當(dāng)高的.在表2中7Li2的a3Σ+u電子態(tài),在實驗上Linton[4]小組獲得了11個振動能級,最高振動能級E(νmax=10)=333.2690cm-1,這個值和實驗值Dexpt=333.690cm-1非常接近,說明實驗數(shù)據(jù)是正e確可靠的,因此用AM方法同樣也得到了11個振動能譜,我們可以斷定該電子態(tài)的完全振動能譜只有11個,分別將最高三個精確的實驗?zāi)芗壓虯M振動能級代入新公式(1)計算得到的離解能分別為DAexMpt=333.6557cm-1和DAeM=333.2904cm-1,這些結(jié)果和Deexpt相比,誤差都非常小.
表2給出了所列各電子態(tài)的AM完全振動能譜
{Eν},對這些電子態(tài)而言,所有AM振動能譜都滿足判據(jù)條件(6)—(11)式,因此用AM方法獲得的振動能譜是體系正確的完全振動能譜.這些振動能譜都是使用AM方法,利用各電子態(tài)的已知實驗?zāi)芗壸蛹希跡ν]求解方程(3)而不用任何數(shù)學(xué)近似和物理模型而得到的(所選8個振動能級均以黑體標(biāo)出).因此這些電子態(tài)的振動光譜常數(shù)、振動能譜不僅能很好地重復(fù)已知實驗?zāi)芗壔騌KR(Rydberg-Klein-Rees)數(shù)據(jù),而且能夠合理地產(chǎn)生實驗難以得到的所有高激發(fā)態(tài)振動能級,從而獲得非常接近于離解極限的最高振動能級和最大振動量子數(shù)νmax.比較典型的例子是表2中7Li2的電子激發(fā)態(tài)b3Πu,Linton課題組[24]于1990年從實驗中對該電子態(tài)僅獲得28個振動能級,即最高振動能級E(νmax=27) =7984.1370cm-1,這個值遠(yuǎn)小于離解能的實驗值=12178.9050cm-1,而用AM方法得到該電子態(tài)的完全振動能譜{Eν}的能級數(shù)為62個,最高振動能級EAM(νmax=61)=12178.6680cm-1,將最高三個AM振動能譜代入新公式(1)計算得到的離解能為DAeM=12178.8763cm-1,此值和Deexpt相比,誤差僅為0.0287cm-1,可見獲得的AM數(shù)據(jù)是正確可靠的.
表2 Li2部分電子態(tài)的AM振動能譜和文獻(xiàn)發(fā)表的振動能量值(單位:cm-1)
續(xù)表2
對堿金屬Li2的很多電子態(tài)而言,其完全振動光譜特別是高振動激發(fā)能級還非常缺乏,而且無論是實驗上還是利用理論方法也很難獲得這些重要的振動信息.本文從實驗上獲得的各電子態(tài)的部分已知精確實驗?zāi)芗壋霭l(fā),應(yīng)用精確研究雙原子分子離子的完全振動能譜的代數(shù)方法(AM),結(jié)合基于LeRoy和Bernstein[19,20]研究分子離解能的理論基礎(chǔ)而建立的計算分子離解能的新公式,通過嚴(yán)格求解代數(shù)方程,獲得了7Li2的a3Σu+,33Σg+,13Δg,b3Πu和6Li2的11Πg等5個電子態(tài)的精確振動光譜常數(shù)、完全振動能譜和理論離解能,從而為那些需要Li2高振動激發(fā)態(tài)的精確能級的研究工作提供了重要的研究數(shù)據(jù).
[1]Konowalow D D,F(xiàn)ish J L 1983 Chem.Phys.77 435
[2]Schmidt-Mink I,Müller W,Meyer W 1985 Chem.Phys.92 263
[3]Poteau R,Spiegelmann F 1995 J.Mol.Spectrosc.171 299
[4]Linton C,MartinF,RossAJ,RussierI,CrozetP,Yiannopoulou A,Li L,Lyyra A M 1999 J.Mol.Spectrosc.196 20
[5]Cacciani P,Kokoouline V,Bouloufa N,Masnou-Seeuws F,Vetter R 2003 Phys.Rev.A 68 042506
[6]Peng Q,Lazarov G,Lyyra A M,Jeung G H 2008 J.Mol. Spectrosc.247 184
[7]Li D,Xie F,Li L,Lazoudis A,Lyyra A M 2007 J.Mol. Spectrosc.246 180
[8]Coxon J A,Melville T C 2006 J.Mol.Spectrosc.246 235
[9]Yan T M,Han Y C,Yuan K J,Cong S L 2008 Chem.Phys. 348 39
[10]Russier I,Yiannopoulou A,Crozet P,Ross A J,Martin F,Linton C 1997 J.Mol.Spectrosc.184 129
[11]Lazarov G,Lyyra A M,Li L 1997 J.Mol.Spectrosc.205 73
[12]Barakat B,Bacis R,Carrot F,Churassy S,Crozet P,Martin F 1986 J.Chem.Phys.102 215
[13]Sun W G,Hou S L,F(xiàn)eng H 2002 J.Mol.Spectrosc.215 93
[14]Sun W G,Ren W Y,Hou S L 2005 Mol.Phys.103 2335
[15]Sun W G,Liu X Y,Wang Y J,Zhan Y,F(xiàn)an Q C 2007 Prog. Phys.27 151(in Chinese)[孫衛(wèi)國、劉秀英、王宇杰、詹妍、樊群超2007物理學(xué)進(jìn)展27 151]
[16]Ren W Y,Sun W G 2005 Acta Phys.Sin.54 594(in Chinese)[任維義、孫衛(wèi)國2005物理學(xué)報54 594]
[17]Sun W G,F(xiàn)an Q C,Ren W Y 2007 Sci.Chin.Ser.G 50 611
[18]Fan Q C,Sun W G 2009 Spectrochi.Acta Part A 72 298
[19]LeRoy R J,Bernstein R B 1970 J.Chem.Phys.52 3869
[20]LeRoy R J,Bernstein R B 1971 J.Mol.Spectrosc.37 109
[21]Linton C,Murphy T L,Martin F,Bacis R,Verges J 1989 J. Chem.Phys.91 6036
[22]Linton C,Martin F,Bacis R,Verges J 1990 J.Mol.Spectrosc. 142 340
[23]Yiannopoulou A,Urbanski K,Lyyra A M,Li L,Ji B,Bahns J T,Stwalley W C 1995 J.Chem.Phys.102 3024
[24]Russier I,Yiannopoulou A,Crozet P,Ross A J,Martin F,Linton C 1997 J.Mol.Spectrosc.184 129
[25]Xie X,F(xiàn)ield R W 1986 J.Mol.Spectrosc.117 228
[26]Kaldor U 1990 Chem.Phys.140 1
PACC:3310,3520P
*Project supported by the National Natural Science Foundation of China(Grant No.10774105),the Youth Fund of the Education Department of Sichuan Province,China(Grant No.09ZB083),and the Research Fund of the Key Discipline of Atomic and Molecular Physics,Xihua University,China.
?Corresponding author.E-mail:weiguosun@x263.net
The full vibrational spectra and dissociation energies of Li2molecule*
Fan Qun-Chao1)Sun Wei-Guo1)2)?Li Hui-Dong2)Feng Hao1)
1)(Research Center for Advanced Computation,School of Physics and Chemistry,Xihua University,Chengdu610039,China)
2)(Institute of Atomic and Molecular Physics,Sichuan University,Chengdu610065,China)
(Received 23 November 2009;revised manuscript received 30 January 2010)
Based on the algebraic method(AM)proposed to calculate accurate full vibrational energies for diatomic molecules in Sun’s previous work,the parameter-free analytical formula for dissociation energy of diatomic system is used to study the dissociation energies and the full vibrational spectra of 5 electronic states of Li2molecule in this work.The results show that the AM vibrational spectra and dissociation energies are in excellent agreement with experimental results.
algebraic method,Li2,vibrational energy,dissociation energy
book=353,ebook=353
*國家自然科學(xué)基金(批準(zhǔn)號:10774105)、四川省教育廳青年基金(批準(zhǔn)號:09ZB083)和西華大學(xué)原子與分子物理重點學(xué)科基金資助的課題.
?通訊聯(lián)系人.E-mail:weiguosun@x263.net