水液相下依達拉奉與超氧陰離子自由基反應(yīng)的密度泛函理論

張雪嬌,楊 應(yīng),楊文富,張 勇,姜春旭,王佐成,董雷剛

(1.白城師范學(xué)院 理論計算中心,吉林 白城137000;2.白城師范學(xué)院 傳媒學(xué)院,吉林 白城137000)

1 計算方法

采用M06-2X雜化泛函[15-16]方法,在三重分裂重原子加彌散基組6-311+G(d,p)[17-18]水平全優(yōu)化反應(yīng)勢能面上駐點物種,反應(yīng)在生命體內(nèi)環(huán)境310.15 K和1個標準大氣壓的水液相下進行.通過對反應(yīng)過渡態(tài)[19]虛頻振動模式的分析(紅外振動有且只有一個虛頻)和內(nèi)稟反應(yīng)坐標(IRC)計算[20],確認過渡態(tài)與相關(guān)穩(wěn)定點的關(guān)聯(lián)性.為得到相對精確的高水平反應(yīng)過程勢能面,采用能給出體系較精確能量的更高精度雜化泛函MN15[21]方法,在6-311++G(2df,pd) 基組下計算單點能.駐點全參數(shù)優(yōu)化及單點能的水溶劑效應(yīng)均采用SMD(solvation model density)極化連續(xù)介質(zhì)模型[22]方法處理,其中水是具有強極性的連續(xù)介質(zhì).

抽H和加成反應(yīng)相關(guān)物種的總Gibbs自由能是自由能熱校正和單點能之和,勢能面上反應(yīng)物為零勢點,其他物種的能量取其與反應(yīng)物的相對能量.利用自然鍵軌道理論(NBO)計算主要物種的自然電荷分析(NPA)電荷,利用分子中的原子理論(AIM)對駐點物種的電子密度進行拓撲分析.

電子轉(zhuǎn)移過程無過渡態(tài),自由能壘采用Marcus理論[23-24]計算,計算公式為

(1)

式中λ為重組能,計算公式為

(2)

采用Gaussian16程序[25]進行駐點結(jié)構(gòu)優(yōu)化、NPA電荷及單點能計算,采用AIM 2000程序[26]計算鍵臨界點(BCP)和環(huán)臨界點(RCP)的電荷密度ρBCP和ρRCP及其Laplace值2ρ.

2 結(jié)果與討論

圖1 Eda(C10H10N2O)的幾何構(gòu)型Fig.1 Geometric conformation of Eda (C10H10N2O)

2.1 抽H過程

圖3 超氧陰離子自由基抽H10和超氧陰離子水分子簇抽H22的反應(yīng)過程Fig.3 Reaction process of superoxide anion free radicals extracting H10 and superoxide anion water molecule cluster extracting H22

圖4 超氧陰離子自由基(水分子簇)抽H反應(yīng)的勢能面Fig.4 Potential energy surface of superoxide anion free radicals (water molecule clusters) extracting H reaction

2.2 加成過程

圖5 超氧陰離子自由基加成反應(yīng)過程Fig.5 Process of superoxide anion free radicals addition reaction

圖6 超氧陰離子自由基加成到C12和C13的反應(yīng)過程Fig.6 Reaction process of superoxide anion free radicals addition to C12 and C13

圖7 超氧陰離子水分子簇加成反應(yīng)過程Fig.7 Process of superoxide anion water molecule cluster addition reaction

圖8 超氧陰離子自由基(水分子簇)加成反應(yīng)的勢能面Fig.8 Potential energy surface of superoxide anion free radicals (water molecule clusters) addition reaction

2.3 單電子轉(zhuǎn)移過程

圖9 單電子從Eda向轉(zhuǎn)移的過程Fig.9 Process of single electron transferring from Eda to

3 結(jié) 論