從頭算分子動力學研究O-與CH3F反應的產物通道

吳琍霞 于 鋒 劉 靜 戴靜華 周曉國 劉世林

(合肥微尺度物質科學國家實驗室(籌),中國科學技術大學化學物理系,合肥 230026)

從頭算分子動力學研究O-與CH3F反應的產物通道

吳琍霞 于 鋒 劉 靜 戴靜華 周曉國*劉世林

(合肥微尺度物質科學國家實驗室(籌),中國科學技術大學化學物理系,合肥 230026)

在B3LYP/6-31+G(d,p)理論水平下采用基于波恩-奧本海默近似的從頭算分子動力學方法重新研究了O-與CH3F反應經抽氫生成OH-和生成H2O的兩條產物通道.反應軌線從反應初始過渡態開始,采用300 K時的熱取樣確定初始條件,同時為對比不同的初始碰撞平動能條件下產物通道的變化,分別限定過渡矢量上的能量為2.1、36.8及62.8 kJ·mol-1進行軌線計算,所有軌線計算的結果表明抽氫生成OH-的過程始終為主要的產物通道.我們的計算不僅進一步證實了以往實驗的結論,而且描繪了抽氫生成OH-和生成H2O這兩個產物通道在反應出口勢能面上的動態反應路徑,更為深刻地揭示了該反應的微觀機理.

從頭算分子動力學;氧負離子;氟甲烷;反應機理;過渡態

負離子與分子的反應在大氣化學、燃燒化學、有機化學等研究領域中占有重要的地位[1-4],深入了解這些反應的機理有助于我們更清楚地認識離子-分子化學反應的本質,并為探索液相化學反應過程提供重要的線索.氧負離子(O-)是一種高反應活性的負離子,它所參與的化學反應是有機含氧負離子的最直接來源,在有機化學、電離層大氣化學、質譜探測分析等諸多領域都有大量涉及,因此吸引了實驗工作的廣泛研究[5].

在眾多的相關反應中,O-和鹵代烷烴的反應作為典型代表受到了格外的關注.Futrell和Tiernan[6]利用串聯質譜研究了在較低的碰撞平動能條件下O-與CH3F的反應,由測量得到的反應速率常數推知產物離子的分布約為OH-(82%)和CHF-(18%). Tanaka等[7]在流動余輝裝置中測量得到的產物離子分支比為OH-占98%、CHF-和F-均占1%.Dawson等[8]在離子回旋共振質譜實驗中觀測得到的結果則為OH-占80%、CHF-占20%.Peverall等[9]則采用離子選擇流動管研究了300 K和66.67 Pa條件下O-和O-2與鹵代烷烴分子的反應,獲得了反應速率常數和產物通道分支比等信息,其中O-與CH3F的反應速率常數為1.4×10-9cm3·molecule-1·s-1,離子產物分別是OH-(68%)、CHF-(31%)和F-(1%).顯然,生成OH-和CHF-離子的過程是該反應的主要產物通道,而生成F-的反應通道分支比非常小.但是,不同實驗觀測的結果也存在著一些差異,而產生差異的原因被認為是由于實驗條件和負離子探測手段的不同而造成的,其中最重要的因素可能是在流動管等實驗研究過程中普遍存在的二次反應的影響[5].例如,在O-與CH3F的反應過程中就可能發生初生產物OH-與CH3F的反應,進而影響OH-和F-這兩種負離子產物的分支比.

由于質譜實驗僅能給出離子產物的部分信息,因此真實的化學反應過程顯然無法從這些有限的信息中得到準確的描述.相對而言,高水平的量子化學理論計算則能夠為澄清反應機理提供全面的信息,同時也能夠證實實驗觀測結果的合理性.這樣,開展該反應的理論研究顯然對于揭示其反應機理、深入理解分子－離子作用機制十分必要.事實上, Yamamoto等[10]就曾采用密度泛函理論(DFT)計算了O-與CH3F反應的部分勢能面特征,試圖闡明實驗觀測的離子產物生成機理.針對這個反應中可能存在的如下三個產物通道:

Yamamoto等分別在B3LYP/6-31+G(d,p)和B3LYP/ aug-cc-pVTZ理論水平下優化了反應勢能面上中間體和過渡態的分子結構,結合內稟反應坐標(IRC)理論計算給出了該反應的勢能剖面和反應機理.此外,通過比較O-與CF4、CHF3、CH2F2和CH3F的反應速率常數,他們認為離子偶極相互作用是決定O-與這類分子發生碰撞反應時相互作用的關鍵因素.然而,他們未能在此基礎上進一步研究該體系的動力學反應特征.

最近,我們分別利用IRC和基于波恩-奧本海默近似的分子動力學(Born-Oppenheimer molecular dynamics,BOMD)理論計算對比研究了該反應體系SN2產物通道的靜態和動態反應路徑,揭示了一條全新的動力學反應產物通道,即生成HF和CH2O-的過程[11].有意思的是,通過對靜態和動態反應路徑的比較,我們發現加在該反應體系原子核上的平動能大小對動力學反應路徑產生決定性的影響,因此像Yamamoto等單純靠尋找最小能量反應途徑的方法來研究反應機理是遠遠不夠的.更深入的計算還表明,SN2過渡態的過渡矢量與其余“實”振動模式之間的耦合強弱是決定生成F-+CH3O還是生成HF+CH2O-產物通道的根本因素.

相對分支比小的SN2產物通道,抽氫生成OH-和生成H2O+CHF-的過程是該反應的主要產物通道,因此開展其相關的反應動力學研究對于全面了解反應機理、判斷實驗觀測結果的合理性具有更為重要的作用.為此,我們采用從頭算分子動力學方法重新研究了這兩個產物通道,并且通過改變加在反應初始過渡態的過渡矢量上的能量來模擬在不同碰撞條件下的反應情形.最終,通過與Yamamoto等反應勢能剖面的比較能夠更為深刻地揭示了該反應的微觀反應機理,并為進一步研究O-與其他有機分子的反應機理提供了參考.

1 計算方法

本文的計算工作均在Gaussian 03程序包[12]中完成.首先采用密度泛函理論中的B3LYP方法[13-14],在6-31+G(d,p)水平下優化了O-離子進攻CH3F中氫原子的過渡態分子構型,即Yamamoto等[10]獲得的反應勢能面上對應抽氫生成OH-和生成H2O+ CHF-的過程的初始過渡態TS1.同時,為了研究不同方法和基組水平對該過渡態優化構型可能存在的影響,采用二階微擾MP2方法[15-16]及較大的基組aug-cc-pvTZ[17]重新進行構型優化,從而對比該過渡態在B3LYP/aug-cc-pvTZ,MP2/6-31+G(d,p)及MP2/ aug-cc-pvTZ水平下的優化構型.最終,從B3LYP/6-31+G(d,p)水平下獲得的TS1的優化構型開始,在相同計算水平下進行BOMD計算[18-20].具體計算過程中采用了預測-校正算法,并在計算新的解析Hessian矩陣(能量的二階偏導數矩陣)前在5個梯度點上更新Hessian[21-22].動力學軌線計算的初始條件由300 K溫度下的熱取樣決定[23-24],即針對TS1的振動取樣溫度和轉動取樣溫度均設為300 K,該溫度與前人實驗條件[9]相符.計算過程中,軌線的步長設定為0.25 amu1/2·bohr,而對于每一條軌線我們設定的最大計算步數為5000.當生成的產物很好地分離(即不同產物之間的質心距離超過13 bohr)或者計算的總步數超過5000時,軌線計算將自動終止.

除了正常的熱取樣確定初始條件外,為了模擬O-與CH3F在較低、中等和較高的初始碰撞平動能條件下產物通道的變化,我們還對加在TS1虛頻振動模式即過渡矢量(transition vector)上的能量加以限定比較.具體計算中,分別采用2.1、36.8及62.8 kJ·mol-1平動能代表反應物在較低、中等和較高碰撞能的初始條件.

一般來講,要獲得多通道化學反應的精確通道分支比至少需要計算幾百條反應軌線.如果能夠擬合得到反應的全維勢能面,那么實現大量軌線的計算是可行的.然而,對O-+CH3F的反應體系而言,構建其全維反應勢能面將非常困難.這主要是因為以下三個原因:首先,在該體系勢能面上的一些區域自恰場計算極難收斂;其次,由于六原子反應體系,自由度較多,構建高精度的全維勢能面計算量非常大;更為重要的是,由于對該類反應體系的作用機制并不明確,因此缺乏恰當的擬合函數.在這種情況下,采用BOMD理論方法可以避開全維反應勢能面的構建,大大節省了計算資源.我們前期的理論計算工作[11]也證實了其結果的可靠性.當然,即便如此這種動力學計算同樣需要非常大的計算量.考慮到計算時間和軌線計算過程中的自恰場計算收斂的問題,我們總共計算了97條軌線.其中,利用正常的熱取樣確定初始條件,不對分配在過渡矢量上的能量做任何限定條件下,總共計算了57條軌線,用以評估實際反應中的主要產物通道;當限定在過渡矢量上的能量為2.1 kJ·mol-1時,計算了10條軌線,而限定在過渡矢量上的能量為36.8和62.8 kJ·mol-1時,則分別計算了15條軌線,由此對比反應在不同的初始碰撞平動能條件下產物通道的變化.需要說明的是,盡管當前的計算中僅涉及有限數目的反應軌線,但是如下文所示,我們的結果依然和實驗結果符合得很好.

此外,軌線計算所得產物的電荷分布由Mulliken電荷布居分析[25]確定.我們以往的理論工作[11,26-28]也證實了該方法可以很好地描述離子－分子反應中產物的電荷分布特征.本文工作中,我們以軌線達到停止標準時的分子結構為研究對象,對該分子結構分別做單點能量計算和Mulliken電荷布居分析,根據獲得的負電荷分布狀況確定最終的離子產物和中性碎片.

2 結果與討論

2.1 抽氫反應通道和生成水反應通道的入口勢能面特征

在離子誘導偶極作用下,當O-靠近CH3F的氫原子端時可迅速形成離子絡合物CH3F…O-,之后經初始過渡態TS1發生抽氫過程,生成離子偶極絡合物OH-…CH2F.圖1顯示了在B3LYP/6-31+G(d,p), B3LYP/aug-cc-pvTZ,MP2/6-31+G(d,p)及MP2/augcc-pvTZ水平下優化得到的TS1分子構型.顯然,各理論水平下分子優化構型沒有明顯差異,這表明當前反應途徑中,采用B3LYP/6-31+G(d,p)方法來處理是合理的.在B3LYP/6-31+G(d,p)獲得的優化構型中,C和H以及O和H之間的距離分別為0.1365和0.1146 nm,而O、H和C之間的夾角為151.79°.

圖1 O-+CH3F反應的初始階段過渡態TS1的優化分子構型Fig.1 Optimized geometry of the initial transition state TS1 of the O-+CH3F reactionThe geometry parameters from the upper to the bottom are obtained at the B3LYP/6-31+G(d,p),B3LYP/aug-cc-pvTZ,MP2/6-31+G(d,p), and MP2/aug-cc-pvTZ levels,.respectively.

正如Yamamoto等[10]的計算所示,抽氫和生成水的反應產物均是由絡合物OH-…CH2F經碰撞解離或進一步異構解離得到.具體而言,絡合物OH-…CH2F可以經碰撞直接解離,生成OH-+CH2F(即抽氫產物通道);絡合物OH-…CH2F還可經另一過渡態異構為H2O與CHF-形成的離子偶極絡合物H2O…CHF-,繼而碰撞解離生成產物H2O和CHF-.

2.2 抽氫反應和生成水產物通道的反應軌線總結

從反應入口的初始過渡態TS1開始,在正常的熱取樣確定初始條件,并不對分配在過渡矢量上的能量做任何限定的條件下,我們總共計算了57條軌線,如表1所示.計算獲得的軌線中總能量變化在0.16-7.33 kJ·mol-1范圍內,始終低于化學精度(8.4 kJ·mol-1),而總的角動量變化則在10-9?-10-7?(?為約化普朗克常量)范圍.盡管通過減小計算步長可以改善能量守恒的情況,但是計算量會相應大大增加,因此我們只需控制反應過程中能量守恒滿足化學精度要求即可.此外,考慮到該反應是在基態雙重態勢能面上發生的,我們也考查了軌線計算過程中的自旋污染情況.其中對于該雙重態體系而言,投影前自旋角動量平方〈S2〉的本征值接近0.75,這說明自旋污染并不嚴重,當前多選取的B3LYP/6-31+G(d,p)計算水平是合適的.最終的反應產物由Mulliken電荷布居分析可知,分別對應OH-+CH2F和H2O+ CHF-產物,而不是OH+CH2F和H2O+CHF.

在這57條軌線中,35條軌線對應生成的產物是OH-+CH2F,6條軌線對應生成H2O+CHF-產物,2條軌線顯示的產物是CH2OH+F-,14條軌線未在5000步時完成產物分離,停留在離子-分子絡合物狀態.當然,隨著反應步數的進一步增加,這14條軌線中的離子-分子絡合物最終分離,分別對應OH-+ CH2F和H2O+CHF-產物.

為了對比反應在不同的初始碰撞平動能條件下產物通道的變化,我們在熱取樣確定初始條件的前提下,分別限定在過渡矢量上的能量為2.1 kJ·mol-1時,計算了10條反應軌線,而限定在過渡矢量上的能量為36.8和62.8 kJ·mol-1時各自計算了15條軌線.在這總共40條軌線中,大部分軌線的總能量變化低于4.0 kJ·mol-1,其中能量守恒最差的一條軌線能量變化達到8.0 kJ·mol-1,但仍然在化學精度的范圍內.總的角動量變化始終在10-8?左右.

表1 初始條件由300 K熱取樣確定的條件下O-與CH3F的反應軌線總結Table 1 Trajectory summary for the reaction of O-with CH3F,with the starting conditions determined by thermal sampling at 300 K

當限定在過渡矢量上的能量為2.1 kJ·mol-1時, O-與CH3F的反應近似是在熱碰撞條件下發生,相應的相對平動能非常小.在我們計算的10條軌線中,6條軌線最終生成產物OH-+CH2F,即抽氫反應通道;2條軌線最終生成產物H2O+CHF-,即生成H2O的反應通道;一條軌線在5000步時停留在H2O…CHF-中間體絡合物,在增加反應時間后最終可以解離為H2O+CHF-;此外,還有一條軌線生成了產物HF+CH2O-,這顯然是由于計算過程中反應軌線耦合到別的產物路徑上去了.

當加在過渡矢量上的能量為36.8 kJ·mol-1時,我們計算的15條反應軌線中有11條生成產物OH-+CH2F,3條生成H2O+CHF-產物,另有1條軌線由于意外的耦合生成產物F-+CH2OH.而當加在過渡矢量上的能量為62.8 kJ·mol-1時,計算的15條反應軌線中11條對應抽氫反應過程,生成OH-+CH2F產物,其余4條軌線對應生成水的產物通道,得到H2O+CHF-的產物.

當前計算的所有軌線及相應的產物歸納在表1中.顯然,在300 K熱取樣條件下,O-+CH3F反應中通過抽氫過程生成OH-+CH2F產物為主要通道,而且即便是考慮到加在初始反應過渡態TS1的過渡矢量上能量變化,抽氫生成OH-離子的產物通道始終占據主要地位.此外,從我們有限數量的軌線計算來看,加在初始反應過渡態TS1的過渡矢量上能量增加,并沒有使抽氫反應或生成水的產物通道分支比有明顯的上升.

2.3 抽氫反應通道和生成水反應通道的動態反應路徑

通常,動態反應通道與靜態反應通道在產物端會有所差異.當前計算所涉及的抽氫產物通道和生成水的產物通道也是如此,其對應的靜態反應途徑在前面已簡要說明,而動態反應過程則略有不同.

圖2列舉了我們計算得到的一條典型抽氫產物通道的反應軌線結果,及反應體系相對勢能隨時間的變化.反應經歷最初的過渡態TS1后,迅速抽氫形成OH-…CH2F中間體絡合物.之后在絡合物解離的過程中,在強烈的離子誘導偶極作用下OH-圍繞CH2F運動持續一段時間,期間甚至形成弱的氫鍵,最終生成OH-和CH2F.

圖2 生成OH-+CH2F產物的典型反應軌線中勢能隨時間變化的曲線Fig.2 Potential energy profile for a typical reaction trajectory leading to OH-+CH2F

類似地,生成H2O+CHF-產物的典型反應軌線中隨時間變化的勢能曲線見圖3.顯然,該產物通道的實現是分步進行的.起初的過程與前述的抽氫通道相似,反應在經歷最初的過渡態TS1后,迅速抽氫形成OH-…CH2F中間體絡合物.之后OH-在圍繞CH2F運動時迅速奪取CH2F上的一個質子形成另一種中間體絡合物H2O…CHF-,隨后在新的絡合物解離過程中,同樣在離子誘導偶極作用下CHF-和H2O分子相伴運動一段時間,最終生成CHF-和H2O.對比這種動態反應機理,其過程與基于勢能剖面的反應機理(靜態反應途徑)[10]是一致的.

由于當前計算中軌線數目相對較少,很難對涉及的中間體絡合物的壽命進行統計分析.然而,在分析具體的反應軌線時,我們發現當過渡矢量上的能量為62.8 kJ·mol-1時,有一條生成OH-+CH2F的軌線僅僅用了大約310 fs的時間就完成了,說明反應過程中OH-…CH2F中間體絡合物的壽命極其短暫,甚至可以認為是從過渡態TS1直接解離生成產物OH-+CH2F的.

2.4 與實驗結果的比較

圖3 生成H2O+CHF-的典型反應軌線中勢能隨時間變化的曲線Fig.3 Potential energy profile for a typical reaction trajectory leading to H2O+CHF-

仔細對比以往的實驗研究,Futrell和Tiernan[6]在串聯質譜中觀測的產物分布為OH-約占82%、CHF-約占18%,而Tanaka等[7]在流動余輝實驗中給出的產物分支比為OH-占98%、CHF-和F-均占1%.考慮到兩者的實驗條件非常類似,均是由近300 K溫度下熱碰撞引起的化學變化,因此兩實驗中O-與 CH3F反應的初始碰撞平動能均較低.然而由于負離子探測效率的不同,兩實驗測量得到的CHF-產物分支比存在較大差異.此外,Dawson等[8]測量得到的產物分布為OH-占80%、CHF-占20%;Peverall等[9]觀測的反應負離子產物主要有OH-(占68%)、CHF-(占31%)和F-(占1%).盡管這些實驗結果彼此之間存在較大的差異,而且實驗測量具有較大的誤差,如Peverall等[9]報道的實驗誤差約為20%,但是所有的實驗研究[6-9]都證實抽氫生成OH-+CH2F是最主要的產物通道,這顯然與我們當前的計算結果吻合得很好.這樣,我們的計算將為進一步研究O-與有機分子的反應機理提供了有價值的參考.

3 結論

我們在B3LYP/6-31+G(d,p)理論水平下采用BOMD方法研究了O-與CH3F反應經抽氫生成OH-和生成H2O的兩個產物通道.通過對比靜態反應途徑和動態反應通道的差別,我們詳細描繪了抽氫生成OH-和生成H2O這兩個產物通道在反應出口勢能面上的動態反應路徑,更為深刻地揭示了該反應的微觀機理,并且證實抽氫生成OH-離子是該反應的主要產物通道,這與以往實驗測量的結論完全一致.

此外,我們還通過改變加在反應初始過渡態的過渡矢量上的能量來模擬在不同碰撞條件下的反應軌線.反應軌線依然從初始過渡態TS1開始,當加在過渡矢量上的能量分別為2.1、36.8以及62.8 kJ· mol-1時,抽氫生成OH-離子的過程始終表現為主要的產物通道,這進一步證實了以往實驗的結論.對該反應體系的計算將為深入研究O-與有機分子的反應機理提供依據和參考.

1 Bowers,M.T.Gas phase ion chemistry.Vol.1.New York: Academic Press,1984:45-82

2 Wayne,R.P.Chemistryofatmosphere.Oxford:Clarendon Press, 1991:225-256

3 Fialkov,A.B.Prog.Energy Combust.Sci.,1997,23:399

4 DePuy,C.H.J.Org.Chem.,2002,67:2393

5 Lee,J.;Grabowski,J.J.Chem.Rev.,1992,92:1611,and references therein

6 Futrell,J.H.;Tiernan,T.O.Ion molecule reactions.Vol.2. Franklin,J.L.Ed.London:Butterworths,1972:485-551

7 Tanaka,K.;Mackay,G.I.;Payzant,J.D.;Bohme,D.K.Can.J. Chem.,1976,54:1643

8 Dawson,J.H.J.;Jennings,K.R.J.Chem.Soc.Faraday Trans 2, 1976,72:700

9 Peverall,R.;Kennedy,R.A.;Mayhew,C.A.;Watts,P.Int.J.Mass Spectrom.Ion Processes.,1997,171:51

10 Yamamoto,M.;Yamashita,K.;Sadakata,M.J.Mol.Struct.-Theochem,2003,634:31

11 Yu,F.;Wu,L.X.;Liu,S.L.;Zhou,X.G.J.Mol.Struct.-Theochem,2010,947:1

12 Frisch,M.J.;Trucks,G.W.;Schlegel,H.B.;et al.Gaussian 03. Reversion B.05.Pittsburgh,PA:Gaussian Inc.,2003

13 Lee,C.;Yang,W.;Parr,R.G.Phys.Rev.B,1988,37:785

14 Becke,A.D.J.Chem.Phys.,1993,98:5648

15 M?ller,C.;Plesset,M.S.Phys.Rev.,1934,46:618

16 Frisch,M.J.;Head-Gordon,M.;Pople,J.A.Chem.Phys.Lett., 1990,166:281

17 Kendall,R.A.;Dunning,T.H.;Harrison,R.J.J.Chem.Phys., 1992,96:6796

18 Thompson,D.L.Modern methods for multidimensional dynamics computations in chemistry.Singapore:World Scientific,1998: 143-189

19 Schlegel,H.B.J.Comput.Chem.,2003,24:1514

20 Schlegel,H.B.Bull.Korean Chem.Soc.,2003,24:837

21 Millam,J.M.;Bakken,V.;Chen,W.;Hase,W.L.;Schlegel,H.B. J.Chem.Phys.,1999,111:3800

22 Bakken,V.;Millam,J.M.;Schlegel,H.B.J.Chem.Phys.,1999, 111:8773

23 Schleyer,P.v.R.;Allinger,N.L.;Clark,T.;Gasteiger,J.; Kollman,P.A.;Scheafer III,H.F.;Schreiner,P.R.Encyclopedia of computational chemistry.Chichester:Wiley,1998:402-407

24 Peslherbe,G.H.;Wang,H.;Hase,W.L.Adv.Chem.Phys.,1999, 105:171

25 Mulliken,R.S.J.Chem.Phys.,1955,23:1833

26 Yu,F.;Zhao,Y.G.;Wang,Y.;Zhou,X.G.;Liu,S.L.Acta Chim. Sin.,2007,65:899 [于 鋒,趙英國,王 勇,周曉國,劉世林.化學學報,2007,65:899]

27 Wang,X.L.;Yu,F.;Xie,D.;Liu,S.L.;Zhou,X.G.Acta Chim. Sin.,2008,66:2499 [王新磊,于 鋒,謝 丹,劉世林,周曉國.化學學報,2008,66:2499]

28 Wang,J.X.;Yu,F.;Liu,J.;Liu,S.L.;Zhou,X.G.Acta Phys.-Chim.Sin.,2008,24:1393 [王俊霞,于 鋒,劉 靜,劉世林,周曉國.物理化學學報,2008,24:1393]

Ab initio Molecular Dynamics Investigation on the Production Channels for the Reaction of O-with CH3F

WU Li-Xia YU Feng LIU Jing DAI Jing-Hua ZHOU Xiao-Guo*LIU Shi-Lin
(Hefei National Laboratory for Physical Sciences at the Microscale,Department of Chemical Physics, University of Science and Technology of China,Hefei230026,P.R.China)

H-atom abstraction and H2O production channels for the reaction of O-with CH3F were reinvestigated using the ab initio molecular dynamics method at the B3LYP/6-31+G(d,p)level of theory and based on the Born-Oppenheimer approximation.The reactive trajectories were initiated at the transition state of H-atom abstraction. Thermal sampling at 300 K was chosen to determine the initial conditions.Additionally,the energies added to the transition vector of the barrier were restricted to 2.1,36.8,and 62.8 kJ·mol-1,separately,to reveal the impact of different initial collision energies on the reaction pathways.The results of all the trajectory calculations demonstrate that the H-atom abstraction channel is the dominant production channel.Therefore,our calculations are consistent with previous experimental conclusions.Furthermore,the dynamic reaction pathways for H-atom abstraction and the H2O production channels on the exit-channel potential energy surface are described based on our calculations and thus a comprehensive reaction mechanism is revealed at the microscopic level.

Ab initio molecular dynamics;Atomic oxygen radical anion;Methyl fluoride;Reaction mechanism; Transition state

O643

Received:April 12,2010;Revised:May 24,2010;Published on Web:July 9,2010.

*Corresponding author.Email:xzhou@ustc.edu.cn;Tel/Fax:+86-551-3600031.

The project was supported by the National Natural Science Foundation of China(20603033,10979042)and National Key Basic Research Program of China(973)(2007CB815204).

國家自然科學基金(20603033,10979042)和國家重點基礎研究發展規劃項目(973)(2007CB815204)資助

?Editorial office of Acta Physico-Chimica Sinica