黃恩玲 裴詩恩 王杰 等

摘 要：在密度泛函DFT/B3LYP理論水平優(yōu)化結(jié)構(gòu)，結(jié)合Aug?CC?pVDZ基組計(jì)算了17種取代乙烷交叉式與重疊式的電子能量。結(jié)果發(fā)現(xiàn)，相關(guān)能差值（ΔEc）與總能量差值（ΔE）的相關(guān)回歸系數(shù)最高（R2=0.95），交換能差值（ΔEx）與 ΔE相關(guān)回歸系數(shù)次之（R2=0.93）。選擇 ΔEc與ΔEx對ΔE進(jìn)行多元線性擬合，得到其相關(guān)回歸系數(shù)（R2=0.95），再用擬合公式計(jì)算ΔE（實(shí)際）， 并與計(jì)算得到的ΔE（理論）進(jìn)行線性擬合， 得到其相關(guān)回歸系數(shù)R2=0.95。我們認(rèn)為， 交叉式取代乙烷的構(gòu)象穩(wěn)定性主要起源于其分子內(nèi)的相關(guān)能（Ec）和交換能（Ex）作用。

關(guān) 鍵 詞：取代乙烷； 交叉式； 重疊式； 穩(wěn)定性； 密度泛函

中圖分類號：O 644.32 文獻(xiàn)標(biāo)識碼： A 文章編號： 1671-0460（2016）03-0460-02

Study on the Stability of Substituted Ethane With Eclipsed

Form and Staggered Form by Density Functional Theory

HUANG En-ling， PEI Shi-en， WANG Jie， WANG Sai-long， CHEN Hao， ZHONG Ai-guo

（Taizhou University， Zhejiang Taizhou 318000，China）

Abstract： Through optimizing the structure in the density functional theory DFT/B3LYP level， combined with Aug-CC-pVDZ basis set， electron energy of 17 kinds of substituted ethane with eclipsed form and staggered form was calculated. The results show that the correlation regression coefficient of correlation energy difference （ΔEc） and total energy difference （ΔE） is the highest （R2 = 0.95）， the correlation regression coefficient of exchange energy difference （ΔEx） and total energy difference （ΔE） is the second （R2 = 0.93）.ΔEc and ΔEx were used to carry out multivariate linear regression on ΔE， correlation regression coefficient was obtained （R2 = 0.95）， and then fitting formula was used to calculate ΔE （actual）， and the linear fitting of ΔE （actual） and calculated ΔE （theoretical） was carried out to obtain its correlation regression coefficient R2 = 0.95. Its pointed out that the conformational stability of substituted ethane with staggered form depends on correlation energy （Ec） and exchange energy （Ex）.

Key words： Substituted ethane； Staggered form； Eclipsed form； Stability； DFT

在藥物分子設(shè)計(jì)和酶促反應(yīng)里，化學(xué)鍵的旋轉(zhuǎn)及其受阻是一重要事件[1]。分子之間或者分子內(nèi)的立體效應(yīng)和超共軛效應(yīng)常常被用來判斷鍵的旋轉(zhuǎn)位壘的起源[2，3]。然而，即便是分子內(nèi)存在較顯著的立體效應(yīng)和超共軛效應(yīng)的最簡單的乙烷分子，人們至今還是沒有取得一致的看法[4，5]。究其原因，一方面是這些效應(yīng)不能使用公式定量表達(dá)出來進(jìn)行比較；另一方面，可能是沒有建立起統(tǒng)一的理論去區(qū)別與明晰這些效應(yīng)及其物理意義。這里我們將嘗試使用密度泛函理論能量分解方法，對最簡單的乙烷及其17種取代衍生物，探討其C-C鍵旋轉(zhuǎn)位壘的關(guān)鍵因子， 明晰交叉式取代乙烷構(gòu)象額外的穩(wěn)定性起源。

1 計(jì)算方法

在G09W軟件的密度泛函DFT理論中， 分子的總能量使用（1）式分解為六項(xiàng)。

E=Ee-e+En-n+En-e+Ek+Ex+Ec （1）

式中： Eee — 電子排斥能；

Enn — 核間排斥能；

Ene — 單電子勢能；

Ek — 單電子動(dòng)能；

Ex — 交換能；

Ec — 相關(guān)能。

選擇乙烷作為母體（見圖1），分別用吸電子和推電子取代基（?F、?Cl、?Br、?OH、?CN、?NO2、?CH3、?NH3）作其一元取代和二元取代，得到17種乙烷的衍生物。在微機(jī)上使用密度泛函DFT/B3LYP理論水平優(yōu)化結(jié)構(gòu)，選擇Aug?CC?pVDZ基組計(jì)算其重疊式與交叉式分子的總能量E，利用公式（2）， 計(jì)算其重疊式與交叉式總能量差ΔE及其分項(xiàng)能量差ΔEi。

E=ΔEee+ΔEnn+ΔEne+ΔEk+ΔEx+ΔEc （2）

（a） 交叉式 （b） 重疊式

2 結(jié)果和討論

表1 顯示了在C-C鍵的旋轉(zhuǎn)過程中，一些在其重疊式與交叉式總能量差ΔE及其分項(xiàng)能量差ΔEi之間的單因子相關(guān)關(guān)系。計(jì)算結(jié)果如表1所示。從表1中可以發(fā)現(xiàn)：在研究的17種取代乙烷中，ΔE理論值總是大于零，意味著乙烷及其取代乙烷（交叉式）總是比其重疊式的分子總能量低，亦即取代乙烷。

為了從中找出能較大程度地反映其能量的變量，分別將ΔEee、ΔEnn、ΔEne、ΔEk、ΔEx、ΔEc對ΔE作單因數(shù)線性回歸分析見圖2（a）-（f）。

ΔEee對ΔE的線性回歸結(jié)果如圖2（a），ΔEee與ΔE為正相關(guān)，但其相關(guān)系數(shù)較小，線性回歸系數(shù)R2=0.68。ΔEnn對ΔE的線性回歸結(jié)果如圖2（b），ΔEn-n與ΔE為正相關(guān)，其相關(guān)系數(shù)較小，線性回歸系數(shù)R2=0.68。ΔEne對ΔE的線性回歸結(jié)果如圖2（c），ΔEne與ΔE為負(fù)相關(guān)，其相關(guān)系數(shù)較小，線性回歸系數(shù)R2=0.68。ΔEk對ΔE的線性回歸結(jié)果如圖2（d），ΔEk與ΔE為正相關(guān)，其相關(guān)系數(shù)較小，線性回歸系數(shù)R2=0.82。

ΔEx對ΔE的線性回歸結(jié)果如圖2（e），ΔEx與ΔE為負(fù)相關(guān)，其相關(guān)系數(shù)很大，線性回歸系數(shù)R2=0.93。ΔEc對ΔE的線性回歸結(jié)果如圖2（f），ΔEc與ΔE為負(fù)相關(guān)，其相關(guān)系數(shù)很大，線性回歸系數(shù)R2=0.93。

3 結(jié) 論

（1）由此認(rèn)為ΔE與ΔEne、ΔEx和ΔEc有關(guān)，但是ΔEne的貢獻(xiàn)不是很大，主要取決于ΔEc和ΔEx。為了能定量描述其間的關(guān)系，我們將ΔEc和ΔEx對ΔE進(jìn)行二元線性擬合，其結(jié)果如圖2所示。