用Gaussian算angle constant

Gaussian套裝軟體是利用量子力學原理計算原子、分子的各式性質。因此，下面要介紹如何求水分子的平衡鍵角以及鍵角常數。

用Gaussian view畫出一團水分子

為了加快運算，我選擇Hartree Fock來計算potential energy surface (PES)，
Basis Sets選6-31G，角度從60掃到160，下面是我的script檔:

%nprocshared=4
%mem=200MW
%chk=C:\Users\Lynn\Desktop\gaussian.chk
# opt=modredundant hf/6-31g geom=connectivity

Title Card Required

0 1
O 0.00000000 0.00000000 0.17474200
H 0.00000000 0.50443600 -0.69896700
H 0.00000000 -0.50443600 -0.69896700

1 2 1.0 3 1.0
2
3

A 2 1 3 S 100 1.000000

跑完之後Results->Scan可以得到下面angle potential vs. theta的趨勢圖

結果平衡角度是111.5與理論值109.5大約差2度。
從上圖迴歸可以得到Angle constant，公式如下;
> y=K(theta-theta0)^2
> logy=logK+2*log(theta-theta0)
(1 Hartree=627.509 kcal mol^-1)

不過，負值的harmonic angle potential似乎是錯的...............................
未完待續