離子體系解析勢能函數、光譜常數和振動能級研究

陳治礎, 岳現房, 孟凡華, 滿忠曉

(1. 曲阜師范大學 物理工程學院, 曲阜273165; 2. 濟寧學院 物理科學與智能工程學院, 曲阜273155)

1 引 言

2 理論方法

2.1 單點能量的計算

2.2 CBS方法

為減少由于基組的不完備性而產生的BSSE，我們使用Halkier等人[22]和Feller[26]應用過的外推方案進行外推：

EN=E∞+A·exp(-αN)

(1)

在式(1)中N代表基組的基數；E∞，系數A和α可以由能量點擬合得到.

2.3 Murrell-Sorbie函數

對于中性雙原子分子和單電荷離子，勢能曲線只有一個最小值，可以用Murrell-Sorbie(M-S)函數[27]進行擬合. M-S函數可由下式描述：

(2)

在式(2)中，ρ=R-Re，R是兩個原子的核間距離，Re是平衡核間距. 基于一系列計算得到的單點能量，平衡核間距Re和擬合參數ai可以通過擬合得到.

3 結果與討論

3.1 解析勢能函數

圖從頭算勢能曲線Fig. 1 Ab initio potential energy curves of

圖2 的M-S勢能曲線和從頭算能量點Fig. 2 M-S potential energy curve and ab initio points of

圖3 的M-S勢能曲線和從頭算能量點Fig. 3 M-S potential energy curve and ab initio points of

(3)

(4)

(5)

(6)

式中De為離解能，9個擬合參數和二至四階力常數已在表1中列出.

表的擬合參數和力常數

3.2 光譜常數

(7)

(8)

(9)

(10)

3.3 振動能級

(11)

式中V(r)為擬合得到的APEF，v和j分別為振動量子數和轉動量子數. 在給定振動能級的基礎上，轉動能級表達式為：

Ev，j=E(v)+Bv[j(j+1)]-

Dv[(j(j+1)]2+Hv[j(j+1)]3+

Lv[j(j+1)]4+Mv[j(j+1)]5+

Nv[j(j+1)]6+Ov[j(j+1)]7

(12)

式中E(v)為體系的振動能級，Bv為體系的慣性轉動常數；Dv，Hv，Lv，Mv，Nv，Ov分別是6個離心畸變常數.

表的光譜常數

表 3 離子在j=0時前10個振動能級E(v)和慣性轉動常數Bv

圖的振動能級Fig. 4 Vibrational levels of

3.4 離心畸變常數

表 4 離子在j=0時前10個振動態的6個離心畸變常數Dv，Hv，Lv，Mv，Nv，Ov

表 5 離子在j=0時前10個振動態的6個離心畸變常數Dv，Hv，Lv，Mv，Nv，Ov

圖的振動能級Fig. 5 Vibrational levels of

4 結 論