Ne原子與HF分子碰撞振轉激發分波截面的研究*

徐 梅 王曉璐 令狐榮鋒 楊向東

1)(貴州師范大學物理與電子科學學院，貴陽 550001)

2)(貴州師范學院物理與電子科學學院，貴陽 550018)

3)(四川大學原子與分子物理研究所，成都 610065)

(2012年3月18日收到;2012年11月8日收到修改稿)

1 引言

原子與分子的碰撞過程和相互作用，是原子與分子物理學中非常重要的研究內容.原子分子碰撞提供的各種原子分子的科學數據，在能源、環境、醫學、材料和生命科學以及國防研究中都發揮著重要的作用，在開拓高新技術產業、推動科技發展等方面有著不可忽視的重要地位.原子與分子相互作用勢在理論計算中非常重要，因為一旦確定了精確勢能面，原則上所有可觀測的物理量均可計算.因此許多實驗和理論工作者都致力于原子與分子系統相互作用勢能面的研究，這方面的工作一直在不斷地創新和發展[1-10].例如，HF分子泛頻激光具有較好的大氣傳輸性能和高亮度特征[11,12]，因而成為強激光體系的重要候選對象;而稀有氣體對HF分子具有弛豫作用[13-15]，在改善HF激光的大氣傳輸性能方面有著廣泛用途;另外用HF的水溶液氫氟酸生產的氟化鹽被廣泛地應用于食品保護、特種金屬冶煉、皮革和紡織品處理、標本保存以及核工業等行業.但是HF是一種不好打交道的化學試劑，它的強腐蝕性對實驗設備損害非常嚴重，因此從理論上研究稀有氣體Ne-HF的振轉非彈性碰撞激發過程有著重要的意義.

本文采用單雙取代并加入三重激發項校正的二次組態相互作用QCISD(T)方法，在aug-cc-pVTZ基組水平上對Ne與HF分子間的相互作用勢進行了計算.并考慮了Boys和Bernardi提出的均衡法(counterpoise method，CM)在計算的基礎上消除基組重疊誤差(BSSE).計算得到了Ne原子與HF分子碰撞系統的相互作用勢在 0°，20°，40°，60°，80°，90°，100°，120°，140°，160°，180°共 11 個方向上，步長取0.01，在2.0—10.0范圍內的891個勢能點數據.接著在Legendre多項式里展開勢能表面，進行矩陣變換，使用Huxley勢函數對數據點進行非線性最小二乘法擬合，計算了各向異性勢的徑向系數 V0，V1，V2，V3等，其函數形式能夠很好地描述 Ne原子與HF分子碰撞系統的勢能面.在此基礎上，采用密耦近似計算得到了Ne原子與HF分子碰撞系統不同能量下的總截面、彈性分波截面、非彈性分波截面.

2 研究的主要理論和方法

式中，E表示碰撞系統的總能量，R表示入射原子Ne與耙分子HF質心的相對距離，r表示分子的核間距離，(+)表示外向邊界條件，α表示一組完整的入射通道量子數.系統的總哈密頓為

如圖 1所示，式中，cosθ=R·r，μNe,HF表示系統的總約化質量，μHF表示HF分子的約化質量，V(R,r,cosθ)是碰撞系統的相互作用勢.

圖1 Ne原子與HF分子碰撞的幾何圖形

式中α為入射通道，γ為γ通道相應的量子數，knαJα可以簡寫為kα，φ(r)為雙子分子振轉波函數的徑向部分，g(R)為原子相對分子質心運動的徑向波函數，Y(R,r)為原子圍繞分子轉動和靶分子中原子轉動的總角函數，n為振動量子數，J為系統的總角動量，lγ表示γ通道中原子繞分子轉動的量子數，Jγ為γ通道中靶分子轉動的量子數.

3 Ne與HF分子間的相互作用勢研究

3.1 HF分子的振轉能級

根據波恩近似，可將分子的幾種主要運動狀態同一定的能量相互聯系起來，分子的總能量則由對應的質心平動能E，電子能Ee，振動能Ev和轉動能Er共同求和得到，即

式中，質心平動能E的改變不會引起分子內部狀態的變化，而另外三項Ee，Ev，Er的能量大小有：Ee～102Ev，Ev～102Er，所以在考慮的碰撞能量下，通常不會涉及到電子狀態的變化情況，僅可能產生振轉激發.根據經過微擾理論驗證的Herzberg的振轉能級規律表達式[16]，雙原子分子的振動光譜項為

式中，ωe表示諧振頻率，ωeχe表示非諧振頻率，Be表示剛性轉動常數，De表示轉動離心常數，αe表示非剛性轉動常數，v表示振動量子數，J表示轉動量子數.

雙原子分子HF的光譜常數如表1所示[17].

表1 HF分子的光譜常數(單位：cm-1)

根據表1和方程(10)式，即可計算得到HF分子的振轉能級的差值，如表2所示.

表2 HF分子振轉能級值(單位：eV)

3.2 HF的勢能函數

采用Murrell-Sorbie[18]函數，對Rydberg[19]提出的勢能函數

經過修正，得到函數：

其中參量a1，a2和a3可由下列公式確定：

式中μ為雙原子分子的約化質量，c是光速，ωe，ωeχe，Be，αe，De為光譜常數.

因此，根據表1和(13)到(18)式，可以計算得到基態HF分子的Murrell-Sorbie勢能函數的參量和力常數，如表3所示.

表3 HF分子勢能函數參量及力常數

表3的計算結果與文獻[20]完全一致.

3.3 Ne原子與HF分子碰撞系統的相互作用勢

采用QCISD(T)方法，在aug-cc-pVTZ基組水平上，分別對Ne與HF分子間相互作用勢進行系統研究，r取優化值0.916808.

運用Boys和Bernardi提出的CM法[21]，為了提高計算的精確性，在計算的基礎上消除基組重疊誤差(BSSE).

在 θ =0°，20°，40°，60°，80°，90°，100°，120°，140°，160°，180°共 11 個方向上的勢能曲線如圖 2所示.

研究原子與雙原子分子的碰撞，我們需要得出體系的球對稱勢和各項異性勢的表達式，一般情況下，碰撞體系的相互作用勢V(R,r,θ)可寫成：

其中Vλ(R,r)為徑向勢函數，Pλ(cosθ)為 Legendre函數，θ為原子相對于雙原子分子的方向角.

圖2 Ne-HF的勢能曲線

V0(R,r)稱為球對稱勢，V1(R,r)，V2(R,r)，V3(R,r)，V4(R,r)等也稱為各向異性勢徑向系數，表征各向異性勢的特征.V(R,r,θ)用矩陣可表示為

將上式左乘 P-1，即得：Vλ(R,r)=P-1V(R,r,θi).

使用Huxley勢函數對Ne原子與HF分子碰撞系統的徑向系數點進行非線性最小二乘法擬合，擬合公式為

式中，擬合參數為：A，Re，a1，a2，a3，a4，C6，C8，C10.前七個徑向系數 Vn(R,r)(n=0，1，2，3，4，5，6)隨 R的變化關系如圖3所示.

圖3 Ne-HF擬合勢的徑向系數

3.4 Ne原子與HF分子碰撞截面的計算

由于HF的高腐蝕性和實驗的難度，Ne原子與HF分子碰撞的分波截面實驗數據一直未見報道，也未見其他理論數據報道.因此，本文旨在為這方面的工作做一點探索，提供一定的參考.

圖4 Ne-HF碰撞系統分波截面(E=50 meV)

圖4和圖5分別是入射能量為50 meV和100 meV時，Ne原子與HF分子碰撞系統的分波截面隨量子數J的變化曲線.其中包括碰撞系統的TCS總截面，彈性00-00分波截面，以及00-01，00-02，00-03等相應系統的非彈性碰撞激發分波截面.

圖5 Ne-HF碰撞系統分波截面(E=100 meV)

圖4和圖5右上角的小圖為非彈性分波截面00-02，00-03等的坐標縮放圖.從圖中可以看出Ne-HF碰撞系統在50 meV與100 meV時，總截面和彈性00-00截面出現一個極大值后，尾部出現很強的尾部效應，然后迅速收斂.這是由于在最靠近R的遠距離處，勢的尾部有很強的吸引支產生的，彈性00-00截面對總截面的貢獻是最大的.50 meV的低轉動躍遷00→01和100 meV的00→01，00→02有兩個峰值，說明低激發態產生尾部效應，高激發態不產生尾部效應.由于非彈性截面也有比較強的尾部效應，所以總截面的尾部效應較高.另外，隨著能量的升高，分波截面的振蕩次數增多，極大值變小，振蕩極大值的位置也向J增大的方向發生移動.尾部效應相對變弱.產生非彈性激發躍遷的幾率將增大，因此，彈性分波截面極大值變小，并且極大值位置向J增加的方向移動，分波截面的振蕩次數增加，收斂速度變慢.

4 結論

本文采用單雙取代并加入三重激發項校正的二次組態相互作用(QCISD(T))方法，以及aug-ccpVTZ基組對He-HF分子間相互作用勢進行了計算，為了提高計算得到的相互作用勢離散點數據的精確性，我們考慮Boys和Bernardi提出的均衡法在計算的基礎上消除基組重疊誤差.計算得到了Ne與HF分子碰撞系統的相互作用勢在0°，20°，40°，60°，80°，90°，100°，120°，140°，160°，180°共 11個方向上，步長取0.01，在2.0—10.0范圍內，每個系統的891個勢能點數據.接著在Legendre多項式里展開勢能表面，進行矩陣變換，得到各向異性勢的徑向系數 V0，V1，V2，···.最終確定使用擬合效果較好，物理意義也比較明確的Huxley勢函數，對11個方向的徑向系數點進行非線性最小二乘法擬合，其函數形式能夠很好地描述Ne原子與HF分子碰撞系統的勢能面.

經過與可查詢的散射截面的實驗值進行了細致的對比，檢驗了擬合勢的可靠性和散射截面密偶計算結果的可靠性.結果顯示：我們構造的勢模型不僅表達形式較為簡潔，而且能夠很好地描述Ne-HF分子碰撞系統的散射截面.曾計算的微分散射截面的衍射振蕩位置和振蕩幅度與實驗光譜測量值具有相同的變化規律，而且符合得很好.因此，我們的擬合勢能夠很好地反映Ne-HF分子碰撞系統的相互作用特征，計算的碰撞激發分波截面的結果是可靠的.

[1]Scole G 1988 Atomic and Molecular Beam Methds(Oxford：Oxford University)

[2]Yu C R，Zhang J，Jiang G S 2009 Acta Phys.Sin.58 2376(in Chinese)(in Chinese)[余春日，張杰，江貴生2009物理學報58 2376]

[3]Huang D H，Wang F H 2009 Acta Phys.Sin.58 6091(in Chinese)[黃多輝，王藩侯2009物理學報58 6091]

[6]Shen G X，Linghu R F，Wang R K，Yang X D 2007 Chin.Phys.16 3352

[7]Wang R K，Shen G X，Song X S，Linghu R F，Yang X D 2008 Acta Phys.Sin.57 4138(in Chinese)[汪榮凱，沈光先，宋曉書，令狐榮鋒，楊向東2008物理學報57 4138]

[8]Wang X L，Linghu R F，Yang J H，LB，GaoT，Yang X D 2012 Acta Phys.Sin.61 093101(in Chinese)[王曉璐，令狐榮鋒，楊建會，呂兵，高濤，楊向東2012物理學報61 093101]

[9]Shen G X，Wang R K，Linghu R F，Yang X D 2011 Acta Phys.Sin.60 013101(in Chinese)[沈光先，汪榮凱，令狐榮鋒，楊向東2011物理學報60 013101]

[10]Shen G X，Linghu R F，Wang R K，Yang X D 2007 Chin.Phys.16 3352

[11]Miller J 1987 AIP Con.Proc.160 10

[12]Smith J H，Robinson D W 1981 J.Chem.Phys.74 5111

[13]Taatjes C A，Leone S R 1988 J.Chem.Phys.89 302

[14]Cuellar E，Parker J H，PrimatelG C 1974 J.Chem.Phys.61 422423

[15]Jiang G，Xie H P，Zhu Z H，Cheng H D，Yang W D，Liu D Y 1998 Chin.J.Lasers 25 817(in Chinese)[蔣剛，謝洪平，朱正和，陳涵德，楊衛東，劉冬燕1998中國激光25 817821]

[16]Herzberg G(Translated by Wang D C)1983 Molecular Spectrum and Molecular Structure(Beijing：Science Press)(in Chinese)[Herzberg G著，王鼎昌譯1983分子光譜與分子結構(北京：科學出版社)]

[17]Huber K P，Herzberg G 1979 Molecular Spectrum and Molecular Structure(IV)Constants of Diatomic Molecules(New York：Van Nostrand Reinhold Company)

[18]Huxley P，MurrellJ N，Sorbie K S 1983 Chem.Soc.Faraday Trans.II 79 323

[19]Rydberg Z R 1931 Physica 73 376

[20]Zhu Z H，Yu H G 1997 Molecular Structure and PotentialEnergy Function(Beijing：Science Press)(in Chinese)[朱正和，俞華根 1997分子結構與分子勢能函數(北京：科學出版社)]

[21]Boys S F，BernadiF 1970 Mol.Phys.19 533