過渡態計算

計算方法:b3lyp/6-31+g(d,p) opt(calcfc,ts,noeigentest) freq


IRC 反應路徑計算

沿著反映路徑一步步搜尋在IRC上的分子結構及能量 , 看結構是否逐漸收斂到反應物與產物

 

 

計算方法:b3lyp/6-31+g(d,p) irc(maxpoints=60,calcfc,stepsize=3,forward)

代表每0.03bohr算一次結構 , 總共算60個點

計算方法:b3lyp/6-31+g(d,p) irc(maxpoints=60,calcfc,stepsize=3,reverse)



電子激發態

CIS : 可以看成是激發態計算的HF方法
CIS(D) : 可以看成是激發態計算的MP2方法

TD-DFT
EOM-CCSD
此兩者誤差約在0.2eV

計算方法:# TD(NStates=20,Root=1) b3lyp/6-31+g(d,p) opt freq

 

TS計算方法
TD(NStates=5,root=1) B3LYP/6-31+g(d,p) opt(TS,noeigentest,calcFC) freq 顯示錯誤,說明為此方法沒有二次微分

改善方法有二:
1.TD(NStates=5,root=1) B3LYP/6-31+g(d,p) opt(TS,noeigentest,ef) freq
 
注意:須使用z-matrix , 且不得超過50項
 z-matrix取得方法 : gaussview開啟input→右鍵→save→類型改成mol2→開啟新儲存的mol2檔案→儲存成.com

2. 先計算Cis(NStates=5,root=1)/6-31+g(d,p) opt(TS,noeigentest,calcfc) freq
 再計算TD(NStates=5,root=1) B3LYP/6-31+g(d,p) opt(TS,noeigentest,readFC) geom=check guess=check freq

 

Hatree normal form TS tautomer form
S0 (ground) -379.888048 -379.783693 -379.868226
S1 (excited) -379.734235
-379.65843
-379.770481

kcal/mol normal form TS tautomer form
S0 (ground)
-238383.359
-238317.8752
-238370.9205
S1 (excited)
-238286.84
-238239.2716
-238309.5844
差距
96.51897
78.60366
61.33611