利用Gaussian求得分子單點能量、Frequency

以C2H5Br+Cl-→C2H5Cl+Br-為例

1.建立input file

利用Gaussview建立input file,複製分子座標
C2H5Br 與 Cl- 距離無窮遠
C2H5Br 與 Cl- (尚未形成過渡態)
Transition state
C2H5Cl 與 Br- (剛形成C2H5Cl)
C2H5Cl 與 Br- 距離無窮遠

point group
C2H5Br與C2H5Cl : Cs
Transition state : C1

2.建立檔案,輸入計算方法&基底函數,將複製的座標貼至檔案內

C2H5Br 與 Cl- 距離無窮遠
C2H5Br 與 Cl- (尚未形成過渡態)
Transition state
C2H5Cl 與 Br- (剛形成C2H5Cl)
C2H5Cl 與 Br- 距離無窮遠

3.離開vi,執行檔案,檢視output file

可以看到以下資訊
Energy(以hartree為單位, 1 hartree=627.51kcal/mol)
如何可以找到 Energy

選用HF/6-31+G(d),計算C2H5Br to transition barrier


參考值barrier來源 :
http://pubs.acs.org/doi/pdf/10.1021/ja960565o




4.將output file 利用gaussview檢視結構、過渡態振動方向

Transition state 振動方向
Freq : -354.80
Infrared : 1015.4041
Raman : 20.6049

補充

以C2H5Br 與 Cl- 距離無窮遠,利用不同計算方法&基底函數所得到的barrier(單位 : kcal/mol)

計算方法/基底函數 C2H5Br 與 Cl- Transition state C2H5Cl 與 Br- C2H5Cl 與 Br- 距離無窮遠 barrier
HF/6-31G -12.75 -2.87 -12.32 2.86 -2.87
HF/6-31G(d) -11.64 2.79 -14.87 -2.40 2.79
HF/6-31+G(d) -11.76 3.65 -16.34 -6.10 3.65
HF/aug-cc-pVDZ -10.50 1.68 -18.82 -10.19 1.68
參考值barrier -0.7

計算方法/基底函數 C2H5Br 與 Cl- Transition state C2H5Cl 與 Br- C2H5Cl 與 Br- 距離無窮遠 barrier
HF/6-31+G(d) -11.76 3.65 -16.34 -6.10 3.65
MP2/6-31+G(d) -11.42 5.23 -16.00 -6.09 5.23
B3LYP/6-31+G(d) -12.32 -1.58 -14.79 -3.75 -1.58
參考值barrier         -0.7




Linux

HF/6-31G


HF/6-31G(d)


HF/6-31+G(d)


HF/aug-cc-pVDZ


MP2/6-31+G(d)


B3LYP/6-31+G(d)

回每週進度