Ⅰ 如何阅读高斯的计算输出文件
如何阅读高斯的计算输出文件
高斯输出文件的关键是看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
最后是整个计算结果的一个总结,各小节之间用\分开
单点计算看明白之后看几何优化的就比较清楚了