Ⅰ 如何閱讀高斯的計算輸出文件
如何閱讀高斯的計算輸出文件
高斯輸出文件的關鍵是看Link,每個Link對應不同的部分,以單點計算為例:單點計算的流程。
輸入行在#後加入P,不然的話就只會顯示Link 1和0,輸出文件一般很長。對於單點計算，只有300多行。
第一部分：
基本運行信息
輸入輸出文件名，初始命令
進入Link 1，顯示進程號
首先是版權說明,
然後是作者,利用guassian計算所發的文章，參考文獻上必須列的
進入
Link 101
讀入輸入文件，
將控制命令轉化為程序能夠看得懂的IOP佔位段的選項和應用設置
Leave Link 101
進入
Link 202：
判斷體系對稱性，並決定實際計算中對稱性應該如何利用
判斷體系點群
例如
Stoichiometry CH2O
Framework group C2V[C2(CO),SGV(H2)]
Deg. of freedom 3
Full point group C2V NOp 4
Largest Abelian subgroup C2V NOp 4
Largest concise Abelian subgroup C2 NOp 2
再將輸入的分子坐標轉換為標准內坐標, 就是
Standard orientation:
這是程序默認，充分利用分子的對稱性來達到方便計算的目的。就是將體系的質量中心放在坐標軸的主軸或者原點上。可以用nosym來禁止這個操作
Leave Link 202
進入
Link 301 產生基組信息
基組函數對稱性
計算出核排斥能
Leave Link 301
進入
Link 302，303 ,這部分是計算積分的具體演算法,一般不會列出詳細過程,一句話帶過
leave Link302,303
進入
Link 401 在實際計算之前，必須有初始猜測，這部分的功能就是產生初始猜測，可以從chk文件讀入，也可以程序自動產生。這部分可以用guess關鍵字來指定。
我這個體系,程序默認是:
Projected INDO Guess(請仔細看手冊,不同的體系Guess是不一樣的)
Leave Link 401
進入
Link 502 （主角登場！）自恰迭代求解SCF方程
設定收斂標准：
除了對包含比Ar重的原子的分子全電子（非ECP）計算外，積分只計算到10-6的精度。l
SCF計算的能量和電子密度都收斂到10-4，或只有能量收斂到10-5為止，無論哪一個先達到。l
給出每一步迭代的信息（？）(這只有在程序輸入中加入#p才會顯示)
如果成功計算完畢
SCF Done,接下來就是列印計算出來的能量,自旋,維力值和收斂判據了.
Leave Link502(一般計算經常會在這一步掛掉,哈哈)
如果你的單點計算只寫了這些,輸出結束,如果你加入了population的內容,那麼,基於scf計算的結果,布局分析開始
進入Link601
Link 601 利用優化後的波函數，
進行Mulliken population,得到
軌道對稱性
電子態
各個軌道的本徵值
然後是一些電荷,自旋密度的分析
leav Link601
最後是整個計算結果的一個總結,各小節之間用\分開
單點計算看明白之後看幾何優化的就比較清楚了