X F2(X=B，N)分子基態的結構與勢能函數*

肖夏杰韓曉琴 劉玉芳

X F2(X=B，N)分子基態的結構與勢能函數*

肖夏杰1)2)韓曉琴3)劉玉芳1)

1)(河南師范大學物理與信息工程學院，新鄉453007)
2)(河南質量工程職業學院，平頂山467000)
3)(商丘師范學院物理與信息工程系，商丘476000)
(2010年8月9日收到;2010年9月10日收到修改稿)

基于Gaussian03計算軟件利用QCISD方法，選用不同基組對X F2(X=B，N)分子基態結構進行了幾何優化，在此基礎上選出最優基組D95(df，pd)和D95+(df，pd)分別對BF2和NF2分子的諧振頻率、力常數等進行了計算．推導出X F2(X=B，N)分子基態的多體展式勢能函數，同時根據勢能函數繪制了X F2(X=B，N)分子基態的3種等值勢能圖，從不同角度驗證了該勢能面符合3原子分子幾何構型，清晰準確的再現了X F2(X=B，N)分子基態的結構特征．

BF2，NF2，結構，勢能函數

PACS:31．50．Bc，33．15．－e，33．15．Fm

1.引言

氟與第Ⅲ主族元素形成的化合物一直是無機與有機化學研究的熱點［1］．BF3是有機合成和石油化工廣泛應用的一種重要酸性催化劑，在很多有機化學反應過程中都有應用［2］，同時BF3是制備鹵化硼、元素硼、硼烷、硼氫化鈉等的主要原料，對BF3理論和實驗研究一直較多［3］．NF3氣體在高能化學激光、電子工業以及太陽能光電產業等方面具有非常廣泛的應用，同時NF3具有在等離子條件下釋放氟的能力，是一種熱力學穩定的氧化劑［4］;隨著科技的飛速發展，對NF3理論和實驗的研究受到人們高度的重視［5，6］．研究BF，BF2，NF，NF2分子的結構與性質是得到BF3，NF3分子勢能函數與相關動力學參數的基礎．1972年Colin等［7］用從頭計算的多組態自洽場理論(MCSCF)研究了BF2，NF2的結構和原子化能，1981年Kurtz小組［8］用相對論Hartree Fock(RHF)理論研究了BF的基態分子結構，1993年Papakondylis等［9］研究了NF2的分子結構，1999年Charles等［3］用B3LYP/6-311+G(2df)方法研究了BF，BF2，BF+， BF+2的結構、諧振頻率、電離能和原子化能，2003年彭謙等［10］用密度泛函方法計算了BF的電子結構與光譜數據，2005年謝安東、朱正和等［2］利用SAC/ SAC方法，對BF分子的勢能函數進行了研究．2008年Jiri Czernek等［11］研究了NF2的分子結構和離解能．盡管還有對BF，BF2，NF，NF2基態結構的實驗與理論研究結果［12—22］，但均未涉及BF2，NF2力常數和勢能函數的研究，用QCISD方法對BF2，NF2分子基態勢能函數的研究，未見報道．

本文基于Gaussian03計算軟件利用QCISD方法，分別選用D95(df，pd)和D95+(df，pd)基組對BF2，NF2分子基態的平衡幾何、諧振頻率、力常數等進行了計算，并在此基礎上推導出了BF2，NF2分子的多體展式勢能函數，同時根據勢能函數討論了勢能面的靜態特征．

2.理論計算與結果討論

2.1.X F2(X=B，N)分子基態的平衡結構、振動頻率與力常數

基于Gaussian03計算軟件利用QCISD方法，選常用基組分別對BF2，NF2分子進行結構優化，得到其鍵長RF—X、鍵角α、離解能De見表1．同時得到它們的基態構型都為C2v，電子態分別為X2A1［3，7，21］和X2B1

［7，9，11］，與文獻一致．

表1 X F2(X=B，N)基態幾何優化計算結果

表1中不同基組計算得到的鍵長RX—F、鍵角α與實驗結果比較，可以看出BF2分子采用D95(df，pd)基組和NF2分子采用D95+(df，pd)基組得到的結果最接近實驗值．分別使用這兩個基組對BF2，NF2分子基態進一步計算，得到振動頻率ν(對稱伸縮振動頻率ν1(α1)、彎曲振動頻率ν2(α1)、反對稱伸縮振動頻率ν3(b1))和力常數f見表2．

表2 X F2(X=B，N)的平衡結構與諧振頻率

2.2.X F2(X=B，N)分子基態的多體項展式勢能函數

根據原子分子反應靜力學原理［13］，推導出BF2(X2A1)、NF2(X2B1)可能的離解極限分別為

對于基態X F2(X=B，N)分子，以R1，R2，R3分別表示X—F，F—X，F—F的間距，選擇基態原子能量為零，則滿足上述離解極限的多體項展式為

其中R1=R2=RXF，R3=RFF，V(X2F)(R1)，V(F2X)(R2)， V(F2F)(R3)分別為兩體項BF(X1Σ+)或NF(X3Σ－)， F2(X1Σ+)的勢能函數，V(3)XF2(R1，R2，R3)為三體項BF2(X2A1)，NF2(X2B1)的勢能函數，兩體項與三體項勢能函數的計算結果如下:

2.2.1.兩體項BF，NF的勢能函數

為了得到基態X F2(X=B，N)分子的多體項展式勢能函數，本文對雙原子分子采用Murrell-Sorbie勢能函數形式:

其中De為雙原子分子的離解能，ρ=R－Re，而R，Re分別為雙原子分子的核間距和平衡核間距，a1，a2，a3為擬合參數．

分別采用QCISD，CCSD，B3 LYP，B3 P86等方法，選常用基組D95+(df，pd)，D95+(2 df，pd)，6-311++g(df，pd)，6-311++g(2 df，pd)，6-311 ++g(3df，2 pd)，6-311++g(3 df，3 pd)等對BF (X1Σ+)，NF(X3Σ－)分子進行幾何優化和頻率計算．通過與文獻值、實驗值的對比，發現對BF (X1Σ+)用CCSD/6-311++g(3 df，3 pd)方法計算得到的結果最優，對NF(X3Σ－)用B3 LYP/6－311 ++g(df，pd)方法計算得到的結果最優．在此基礎上分別對BF(X1Σ+)，NF(X3Σ－)進行單點能掃描，進一步將掃描結果擬合為Murrell-Sorbie勢能函數．BF(X1Σ+)，NF(X3Σ－)，F2(X1Σ+)基態的Murrell-Sorbie勢能函數的各項參數見表3．由表3可以看出擬合得到的BF(X1Σ+)，NF(X3Σ－)基態的Murrell-Sorbie勢能函數的各項參數與實驗值相比達到了較高的精度．

表3 BF(X1Σ+)，NF(X3Σ－)，F2(X1Σ+)的Murrel-Sorbie勢能函數參數

圖1和圖2分別表示出了BF(X1Σ+)，NF(X3Σ－)分子基態的勢能曲線．其中，實線為擬合得到的Murrel-Sorbie勢能函數曲線，圓圈線為單點能掃描曲線．圖1，圖2中，擬合函數曲線與單點能掃描曲線比較，相應位置處的點、線之間符合得較好，說明擬合得到的Murrell-Sorbie勢能函數能夠正確地反映BF(X1Σ+)，NF(X3Σ－)分子基態．由Murrell-Sorbie勢能函數與力常數的關系以及力常數與光譜數據的關系，可以分別得到BF(X1Σ+)，NF(X3Σ－)，F2(X1Σ+)的光譜數據，結果見表4．由表4可知，計算得到的BF(X1Σ+)，NF(X3Σ－)的光譜數據與實驗值相比達到了較高的精度．

圖1 BF基態的勢能曲線

圖2 NF基態的勢能曲線

表4 BF(X1Σ+)，NF(X3Σ－)，F2(X1Σ+)的光譜常數

2.2.2.三體項BF2，NF2的勢能函數

三體項BF2(X2A1)，NF2(X2B1)的勢能函數，采用文獻［22—25］中的形式:

其中，P為對稱內坐標Si的多體項，T為量程函數，它們的形式為

基態BF2(X2A1)，NF2(X2B1)的平衡構型都為C2v對稱結構，為了方便研究其勢能函數，根據勢能面的結構特征，采用優化內坐標．取X F2(X=B，N)的兩個平衡鍵長為參考結構，計算中所使用的內坐標ρi經以下變換成為優化內坐標Si:

式中ρi=Ri－R0i(i=1，2，3)，其中S2對R1和R2的交換是反對稱的，為了滿足這一物理意義，S2只能含偶次項．

由以上可知，要得到三體項BF2(X2A1)，NF2(X2B1)的勢能函數V(3)XFx2(R1，R2，R3)就必需確定7個線性系數(C0，C1，C2，C3，C4，C5，C6)和兩個非線性系數λ1，λ2．對勢能面進行非線性優化，可以確定出兩個非線性系數λ1，λ2，根據表2列出的BF2(X2A1)，NF2(X2B1)的平衡結構常數、離解能及力常數可以確定出(6)式中的7個線性系數(C0，C1，C2，C3，C4，C5，C6)．結果見表5．

表5 BF2(X2A1)，NF2(X2B1)的解析勢能函數的三體項參數

根據基態X F2(X=B，N)分子的多體項展式勢能函數，進一步可繪制出其等值勢能圖，如圖3—8所示．圖像清晰地反映了X F2(X=B，N)分子基態的平衡結構特征．

圖3是根據BF2(X2A1)的解析勢能函數，在固定∠FBF=122.928°的情況下，得到F—B鍵對稱伸縮振動的等值勢能圖．圖中很好地反映了BF2分子基態的離解能及平衡結構特征．在平衡點RBF=RFB=0.131 nm位置處有一勢阱出現，阱深約為3.9374 eV，由此說明在該處形成穩定的BF2分子，這與結構優化結果及得到的離解能計算結果相符．在BF+ F→BF2兩個等價的通道上存在兩個對稱的鞍點，位于(0.14789，0.24263)和(0.24263，0.14789)處，活化能大約為30.875 kJ/mol，并且在RBF=RFB= 0.227 nm位置處有一勢壘，勢壘高度約為0.7 eV，說明反應只能通過兩個等價的通道越過勢壘才能進行．

圖3 BF2分子基態的伸縮振動勢能圖(等值線的間隔為0.38 eV)

圖4 和圖5分別為F原子繞B—F鍵與B原子繞F—F鍵旋轉時形成的等值勢能圖，它們都能表明在RBF=0.227 nm，∠FBF=122.928°處，BF2分子的離解能最低，約為3.9374 eV．由此可以說明在該處形成穩定的BF2分子，與結構優化結果及得到的離解能計算結果相符．

圖6是根據NF2(X2B1)的解析勢能函數，在固定∠FNF=103.35°的情況下，得到F—N鍵對稱伸縮振動的等值勢能圖，很好地反映了NF2分子基態的離解能及平衡結構特征．由圖6可知，NF2分子基態的平衡結構RN—F=RF—N=0.135 nm位置處有一勢阱出現，阱深約為4.7351 eV，由此可以說明在該處形成穩定的NF2分子，準確再現了分子基態的C2v結構，這與結構優化結果及得到的離解能計算結果相符，并且在NF+F→NF2反應中不存在鞍點，該反應為無閾能反應．

圖7是F原子繞N—F鍵旋轉時形成的等值勢能圖，即固定N—F鍵在X軸上，以N—F鍵的中點為原點建立坐標系．由該圖可以看出，當F原子旋轉到距離N原子約0.135 nm處且∠FNF=103.35°時，NF2分子的離解能最低，約為4.7351 eV．

圖4 F原子繞B—F鍵旋轉時的等值勢能圖(等值線的間隔為0.5 eV)

圖5 B原子繞F—F鍵旋轉時的等值勢能圖(等值線的間隔為0.4 eV)

圖6 NF2分子基態的伸縮振動勢能圖(等值線的間隔為1 eV)

圖7 F原子繞N—F鍵旋轉時的等值勢能圖(等值線的間隔為1 eV)

圖8 是N原子繞F—F鍵旋轉時形成的等值勢能圖，即固定F—F鍵在X軸上，以F—F鍵的中點為原點建立坐標系．該圖反映了當N原子旋轉到(0.000 nm，0.135 nm)位置時，NF2分子的能量最低，約為4.7351 eV．由此可以說明在該處形成穩定的NF2分子，準確地再現了NF2分子基態的C2v結構，這與結構優化結果及得到的離解能計算結果相符．

圖8 N原子繞F—F鍵旋轉時的等值勢能圖(等值線的間隔為1 eV)

BF2和NF2分子都是C2v結構且鍵長、鍵角相差無幾，但其形成分子的機理完全不同，BF+F→BF2反應為越過勢壘才能進行的有閾能反應，NF+F→NF2為容易進行的無閾能反應．

3.結論

基于Gaussian03計算軟件利用QCISD方法，分X F(X=B，N)分子基態的平別選用不同的基組對2衡幾何、諧振頻率、力常數等進行了計算，進一步計算得到了X F2(X=B，N)分子基態的離解能，并在此基礎上推導出了X F2(X=B，N)分子基態的多體展式勢能函數，同時根據勢能函數繪制了X F2(X=B，N)分子基態的等值勢能圖進而討論了勢能面的靜態特征．X F2(X=B，N)分子基態都為C2v結構，BF2分子的基電子態為X2A1，平衡核間距為0.13138 nm，鍵角為122.928°，NF2分子的基電子態為X2B1，平衡核間距為0.13532 nm，鍵角為103.35°．根據勢能函數繪制的X F2(X=B，N)分子基態的3種等值勢能圖從不同角度驗證了該勢能面符合三原子分子的幾何構型，清晰準確地再現了X F2(X=B，N)分子基態的結構特征．同時等值勢能圖還反映出X F2(X=B，N)形成分子的機理完全不同．這些不但為研究X F+F(X=B，N)體系的分子反應動力學提供了理論依據，同時也為X F3(X=B，N)的理論和實驗研究奠定了理論基礎．

［1］Yu H T，Huang X R，Chi Y J 2001 Acta Chim．Sin．59 1055 (in Chinese)［于海濤、黃旭日、池玉娟2001化學學報59 1055］

［2］Xie A D，Zhu Z H 2005 Acta Chim．Sin．63 2126(in Chinese)［謝安東、朱正和2005化學學報63 2126］

［3］Charles W，Bauschlicher Jr，Alessandra Ricca 1999 J．Phys．Chem．A 103 313

［4］Peng L P，Wang S B 2007 Chem．Ind．Engng．24 86(in Chinese))［彭立培、王少波2007化學工業與工程24 86］

［5］Zhang H B，Wei G H，Zhou S R 2001 J．Chin．Mass Sp．Soc．22 51(in Chinese)［張洪彬、韋桂歡、周升如2001質譜學報22 51］

［6］Du W H 2009 Low Temp．Spec．Gas．27 14(in Chinese)［杜偉華2009低溫與特氣27 14］

［7］Colin T，Douglas A B 1972 Chem．Phys．Lett．16 573

［8］Kurtz H A，Jordan K D 1981 Chem．Phys．Lett．81 104

［9］Papakondylis A，Mavridis A 1993 Chem．Phys．Lett．216 167

［10］Peng Q，Zhai G H，Wang Y B 2003 Acta Chim．Sin．61 1375 (in Chinese)［彭謙、翟高紅、王育彬2003化學學報61 1375］

［11］Jiri Czernek，Oldrich Zivny 2008 Chem．Phys．Lett．457 54

［12］Moeller M B，Stuart J 1973 Chem．Phys．Lett．19 78

［13］Zhu Z H，Yu H G 1997 Molecular Structuer and Potential Energy Function(Bering:Science Press)p102(in Chinese)［朱正和、俞華根1997分子結構與分子勢能函數(北京:科學出版社)第102頁］

［14］Huber，K P，Herzberg G 1979 Molecular Spectrum and Molecular Structure(Vol．4)(New York:van Nostrand Reinhold Company)250

［15］Yu H T，Huang X R，Chi Y J 2002 Chim J Chin．Univ．23 888 (in Chinese)［于海濤、黃旭日、池玉娟2002高等化學學報23 888］

［16］Muller H S P，Loblein K 2008 J．Mol．Spectrosc．251 185

［17］Karthikeyan B，Bagare S P，Rajamanickam N，Raja V 2009 Astropart．Phys．31 6

［18］Rode M F，Sadlej J 2002 Chem．Phys．Lett．358 237

［19］Brown R D，Burden F R，Godfrey P D 1974 J．Mol．Spectrosc．52 301

［20］Liu Y C，Jiang G，Zhu Z H 2003 Acta Phys．Chim．Sin．19 858 (in Chinese)［劉幼成、蔣剛、朱正和2003物理化學學報19 858］

［21］Yoshimine M，McLean A D M 1967 Int．J．Quantum Chem．1 313

［22］Jiang L J，Liu Y F，Liu Z Z，Hang X Q 2009 Acta Phys．Sin．58 201(in Chinese)［蔣麗娟、劉玉芳、劉振中、韓曉琴2009物理學報58 201］

［23］Shi D H，Zhang J P，Sun J F，Liu Y F，Zhu Z L 2009 Acta Phys．Sin．58 5329(in Chinese)［施德恒、張金平、孫金峰、劉玉芳、朱遵略2009物理學報58 5329］

［24］Zhang J P，Shi D H，Sun J F 2008 J．Atom．Mol．Phys．25 1404(in Chinese)［張金平、施德恒、孫金峰2008原子與分子物理學報25 1404］

［25］Han X Q，Jiang L J，Liu Y F 2010 Acta Phys．Sin．59 1000(in Chinese)［韓曉琴、蔣麗娟、劉玉芳2010物理學報59 1000］

PACS:31．50．Bc，33．15．－e，33．15．Fm

*Project supported by the National Natural Science Foundation of China(Grant No．60977063)，the Innovation Scientists and Technicians Troop Construction Projects of Henan Province，China(Grant No．084100510011)，and the Innovation Talents Program of Institution of Higher Education of Henan Province，China(Grant No．2006 KYCX002)．

Corresponding author．E-mail:yf-liu@henannu．edu．cn

Structure and potential energy function of X F2(X=B，N)molecular ground state*

Xiao Xia-Jie1)2)Han Xiao-Qin3)Liu Yu-Fang1)
1)(Collegel of Physics and Information Engineering，Henan Normal University，Xinxiang 453007 China)
2)(Henan Quality Polytechnics，Pingdingshan 467000 China)
3)(Department of Physics and Information Engineering，Shangqiu Normal College，Shangqiu 476000 China)
(Received 9 August 2010;revised manuscript received 10 September 2010)

Based on the Gaussian03 calculation software，QCISD method is used to optimize the possible ground-state structures of X F2(X=B，N)molecules with the different basis sets．On this basis，the resonant frequencies and force constants are calculated with the D95(df，pd)and D95+(df，pd)basis sets respectively．The potential energy functions of X F2(X= B，N)are derived from the many-body expansion theory．At the same time，according to the potential energy functions，three equivalent potential energy diagrams for the X F2(X=B，N)ground states are drawn．The potential energy surface is verified to be consistent with the three atomic molecular geometry．The structures and the features of X F2(X=B，N)can correctly reappear on the potential surface．

BF2，NF2，structure，potential energy function

*國家自然科學基金(批準號:60977063)、河南省創新型科技人才隊伍建設工程(批準號:084100510011)、河南省高等學校杰出科研人才創新工程(批準號:2006 KYCX002)資助的課題．

．E-mail:yf-liu@henannu．edu．cn