環丙酮分子及離子的解離的研究

崔金玉

（綏化學院電氣工程學院，綏化152061；吉林大學原子與分子物理研究所，長春130012）

1 引 言

酮類也是抗菌藥物的主要成分，如第三代喹諾酮類藥物常用的有氧氟沙星、環丙沙星等曾被廣泛使用，1988 年以來研究出第四代喹諾酮類藥物，如莫西沙星、格帕沙星、加替沙星等抗菌作用更強.研究環丙酮的性質，有利于提煉新的藥物成分.環丙酮是極不穩定的化合物，國內外關于環丙酮的研究報道不是很多，早在1966 年，Hammond［1］和DeBoer等［2］用乙烯酮和重氮甲烷在氟三氯甲烷中成功制備了環丙酮.2007年艾月潔等人［3］結合實驗運用量子化學計算方法CASSCF，B3LYP 和MP2研究在292～365nm 波長的光的激發下，從鍵長鍵角的角度研究環丙酮被激發至激發態S1，然后αC-C鍵斷裂.并且找到在αCC 鍵斷裂途徑上，存在基態和第一激發勢能面的交叉點TS1，然后發生另一個C-C 鍵的斷裂，生成基態產物一氧化碳和乙烯.我們用密度泛函理論B3LYP和TS方法計算環丙酮的解離過程［4，5］.我們的工作一是把分子的解離過程更清晰的描述出來，另外還對離子的解離過程做了計算.主要從計算過渡態角度，用IRC 方法查看過渡態連接的反應物和生成物，確定反應的方向.對環丙酮分子穩定的基態計算，找到兩個過渡態，一個過渡態是C-C單鍵斷裂，欲成開環結構，再由這樣的開環結構，計算得另一C-C鍵斷裂，形成CO 和C2H4（這一過程相當于把環丙酮分子看成等腰三角形，兩腰斷裂）.另一過渡態是環丙酮的另一個C-C鍵（相當于等腰三角形的底邊）斷裂，形成異構體，由異構體解離成C2H2O 和CH2.對環丙酮離子構型的計算，離子構型是不穩定構型，也是一種過渡態，找到其開環結構的穩定離子構型.

2 基礎理論和計算方法

量子力學理論是適用于微觀領域的基本理論，微觀領域的一切現象都可以用量子力學的知識來解釋.所以，我們利用現代量子化學的方法和原子分子理論來研究分子問題［6］.用密度泛函理論［7］UB3LYP方法解決這一問題.密度泛函理論是計算較大分子體系結構的強有力的計算工具，它比Hartree Fock（簡稱HF）和Moller Plesset（簡稱MP）方法節省時間.同樣的方法計算頻率值.在過渡態transition state（簡稱TS）附近開始用opt=TS優化，初始的TS我是利用固定反應坐標進行優化得到的，用freq 進行TS 判據.同時，IRC計算確認過渡態是反應物和產物之間的過渡態.利用計算機軟件Gaussian03，采用基組6-31+G（d，p）首先opt優化環丙酮分子的幾何構型，得到穩定的幾何構型，計算頻率，無虛頻，是穩定構型［10-12］.結構見圖1中.

3 計算結果

3.1 對分子的計算

環丙酮分子式為C3H4O，用Gaussion03軟件計算其幾何構型，結構見圖2 中所示，能量-191.8966Ha，無虛頻.其振動頻率及相對強度如表1所示.從表1中可見頻率為1930.79cm-1時，其振動幅值最大，達358.371KM／Mole.其次是頻率964.151cm-1，強度為108.244 KM／Mole.另外，在頻率為1018.02cm-1時，強度為27KM／Mole.在頻率為1072.33cm-1時，強度為15KM／Mole.其它情況很弱，不去分析，僅對這幾個頻率作分析如下：

表1 環丙酮分子的基態頻率和相對強度Table 1 Frequency and Intensity of Cyclopropanone molecule

（1）在頻率freq=1930.79cm-1時，強度最大358.371KM／Mole.此時1號O 原子和2號C 原子振動最強烈，可能它們與分子其它部分斷開，形成CO 和C2H4，也可能O 與C之間斷開，形成O 和C3H4.但是，這兩種情況都沒有找到過渡態［8］，可能O 與C之間雙鍵結構很難斷開，另外，2C-3C和2C-4C雙鍵同時斷開也很難，也沒有找到雙鍵同時斷開的過渡態，所以很可能先斷開單鍵，成開環結構，然后再斷開另一個C-O 鍵，見下面（2）的分析情況.

（2）在頻率freq=964.15cm-1時，強度為108.242KM／Mole，振動幅度也比較大，分析其振動情況［9］，2C左右擺動幅度很大，可能會使分子的2C-4C或2C-3C之間的鍵斷裂，成開環結構.因為分子對稱，兩種斷裂情況的幾率相同，斷裂后的分子結構、振動情況都相同，所以只分析一種情況即2C-4C 鍵斷裂.尋找這樣的過渡態［13］，得到TS1，見圖1 中所示，其能量為-191.8350 Ha，有一個虛頻值為-642.68cm-1.分析此時的振動情況［14］，2C向外振動劇烈，很可能2C-3C之間鍵斷裂，并且用IRC 優化TS2，所得兩個方向的結果，一個方向可以生成反應物，把另一個方向的結果繼續優化，得產物CO 和C2H4的共同體，不是穩定構型，能量為-191.9171 Ha，進一步分解成產物CO 和C2H4.這一過程見圖1中所示.

圖1 環丙酮分子的解離過程（Ⅰ）Fig.1 Dissociation of Cyclopropanone molecule（Ⅰ）

（3）在頻率為-1018.02cm-1時，強度為27 KM／Mole，也進行了分析，3C-4C 各向分子外振動，尋找3C-4C鍵斷裂的過渡態TS2，結構見圖2中所示，虛頻率為-525.1544cm-1，能量為-191.8287Ha，IRC計算表示，它是反應物和異構體INT1的中間過渡態，異構體INT1見圖2中，是穩定構型，能量為-191.8799 Ha.由INT1的振動模式分析，2C與3C各向分子兩側振動劇烈，欲使2C與3C之間鍵斷開，找到這種情況的過渡態TS3，圖2中所示，能量為-191.6539 Ha，虛頻值為-46.53cm-1，由此生成產物C2H2O 和CH2共同體，無虛頻是穩定構型［15］.但這一TS3勢壘較高，欲解離成產物C2H2O 和CH2遠不如（2）的情況容易.

圖2 環丙酮分子的解離過程（Ⅱ）Fig.2 Dissociation of cyclopropanone molecule（Ⅱ）

（4）頻率為1072.33cm-1時的情況，強度為15KM／Mole.3C-4C 振動方向仍相反，但與（3）情況不同，不是向分子外振動，而是上下振動，分子欲發生扭轉，可能2C-3C 或2C-4C 之間斷裂，情況如2情況相似；也可能3C-4C 鍵斷裂［16］，情況同（3）情況相似.

3.2 對離子的異構化分析計算

將環丙酮分子的穩定結構的簡正坐標作為輸入，對離子結構進行優化［4，5，17］，最后得到過渡態（TS4），能量為-191.5612 Ha，有一個虛頻為-285.23cm-1，說明環丙酮離子極不穩定，觀察虛頻下的振動模式，O 和2C一起向離子左右兩邊振動，且振動強烈，可能2C-3C鍵斷開，也可能2C-4C鍵斷開，兩種情況相同，只分析一種，得離子的開環結構，即穩定的離子的異構體INT2（2C-4C鍵斷開）.見圖3中所示.能量為-191.5798 Ha，無虛頻穩定.分析INT2振動最強烈的情況，1O和2C向離子外振動強烈，計算過渡態，得到過渡態TS5，頻率為-164.15 cm-1，能量值為-191.5525Ha.IRC 顯示TS5是INT2和INT3（2C-3C鍵斷開）之間的過渡態.最后得離子的穩定構型就是開環結構INT2和INT3.

圖3 環丙酮離子的異構化Fig.3 Isomerism of Cyclopropanone ion

表2 環丙酮分子到解離產物及其異構體的各狀態能量值（Hartree）及頻率（cm-1）情況（1Ha=27.2114eV） Table 2 Energy（Hartree）and frequency（cm-1）of from Cyclopropanone molecule to product and isomer（1Ha=27.2114eV）

4 結 語

圖4 環丙酮分子解離過程的勢能面Fig.4 The energies surface of Dissociation of Cyclopropanone molecule and Ion

最后得出環丙酮分子及離子解離過程各狀態能量高低見圖4所示，能量大小見表2中.

本文用密度泛函DFT（B3LYP）計算環丙酮分子及離子穩定結構，并用TS計算過渡態.從計算結果看出，環丙酮分子的解離產物CO 和C2H2，這與艾月潔等人［3］計算結果相同.但是，計算結果表明還有產物C2H2O 和CH2.但從能量上看，前者容易實現（TS1=-191.8350eV），后者過渡態能量大（TS3=-191.6539eV），所以實現起來有些困難，但在能量足夠大的時候，會有產物C2H2O 和CH2生成.另外，實現前者時，單鍵容易斷開，雙鍵同時斷開很難，沒有找到雙鍵同時斷開的過渡態，但找到單鍵斷開的過渡態及解離到產物的過程.同樣的方法計算環丙酮離子C3H4O+，它是不穩構型，有一個虛頻值，而它的開環異構體穩定，我們找到它的異構化過程.所以從以上計算表明，環丙酮離子異構體穩定，分子解離產物是一氧化碳和乙烯，及存在異構體，在能量足夠高時，解離成亞甲基CH2和C2H2O 共同體.

［1］ Turro N J，Hammond W B.Tetramethylcyclopropanone，Ⅱ.Mechanism of the Favorskii Rearrangment［J］.J.Am.Chem.Soc.，1965，87：3258.

［2］ Schaafsma S E，Steinberg H，Deboer T J.The synthesis of cyclopropanone［J］.Rec.Trav.Chim.，1966，85，1170.

［3］ Ai Y J，Lin L，Fang W H.Ab inito study mechanisms of photodissociation and thermal isomerization of cyclopropanone［J］.Acta Chim.Sin.，2007，65（2）：129（in Chinese）［艾玥潔，林 玲，方維.環丙酮光解離和熱異構機理的從頭算研究［J］.化學學報，2007，65（2）：129］

［4］ Cui J Y，Jin MX，Wang Z G，et al.Dissociation of 2-methyl-cyclopentanone Ion［J］.At.Mol.Phys.，2009，26（5）：787（in Chinese）［崔金 玉，金明星，王志剛，等.2-甲基環戊酮離子的解離［J］.原子與分子物理學報，2009，26（5）：787］

［5］ Cui J Y.Dissociation of the first six excited states of 2-methyl-cyclopentanone ion excited state［J］.J.Mol.Sci.，2011，27（3）：199（in Chinese）［崔金玉.2-甲基環戊酮離子前六個激發態的解離［J］.分子科學學報，2011，27（3）：199］

［6］ Tang A Q，Yang Z Z，Li Q S.Quantum chemistry［M］.Beijing：Science Press，1982.（in Chinese）［唐敖慶，楊忠志，李前樹.量子化學［M］.北京：科學出版社，1982］

［7］ Jiang F L.Principles of quantum chemistry ［M］.Shanghai：Fudan University Press，1990（in Chi-nese）［江逢林.量子化學原理［M］.上海：復旦大學出版社，1990］

［8］ Roothan C C J，Bagus P S.Methods of molecular quantum mechanics［M］.Academic Press，1964.

［9］ Gunnarsson O，Jones R O.Density Functional Calculations for Atoms，Molecules and Clusters［J］，Phys.Scr.，1980，21：394.

［10］ Becke A D.Density-functional thermochemistry.Ⅲ.the role of exact exchange［J］.Chem.Phys.，1993，98：56.

［11］ Lee C，Yang W，Parr R G.Development of the Colle-Salvetti correlation-energy formula into a functional of the electron-density［J］.Phys.Rev.B，1988，37：785.

［12］ Simons J P.The photochemical and spectroscopic approach to molecular reaction dynamics［J］.Comments on At.Mol.Phys.，1985，16（3）：157.

［13］ Leone S R.Laser probing of chemical reaction dynamics［J］.Science，1985，227.

［14］ Bernstein R B.Chemical dynamics via molecular beam and laser techniques［M］.New York：Oxford University Press，1982.

［15］ Smith I V M.Kinetics and dynamics of elementary gas reaction［M］.London：Butter Worth，1980.

［16］ Wang Z G，Pan S F，Yao M G.The two hydrogen transfer dissociation channels of nicotine ［J］.J.Mol.Struct.：Theochem.，2006，767：99.

［17］ Jin M X，Cui Q L，Emad Mukhtar，et al.Secondharmonic-generation measurements on ZnSe under high pressure ［J］.J.Phys.：Condens.Matt.，2002，14：11037.