如何阅读高斯的计算输出文件

如何阅读高斯的计算输出文件

高斯输出文件的关键是看Link,每个Link对应不同的部分,以单点计算为例:单点计算的流程。

输入行在#后加入P,不然的话就只会显示Link 1和0,输出文件一般很长。对于单点计算，只有300多行。

第一部分：

基本运行信息

输入输出文件名，初始命令

进入Link 1，显示进程号

首先是版权说明,

然后是作者,利用guassian计算所发的文章，参考文献上必须列的

进入

Link 101

读入输入文件，

将控制命令转化为程序能够看得懂的IOP占位段的选项和应用设置

Leave Link 101

进入

Link 202：

判断体系对称性，并决定实际计算中对称性应该如何利用

判断体系点群

例如

StoichiometryCH2O

Framework group C2V[C2(CO),SGV(H2)]

Deg. of freedom3

Full point groupC2VNOp 4

Largest Abelian subgroup C2VNOp 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

最后是整个计算结果的一个总结,各小节之间用\分开

单点计算看明白之后看几何优化的就比较清楚了

同时打开你的输入文件和输出文件，比较一下，看你的输出文件是否是你所想要的，如果是，你可以改动一下输出文件的键长，键角，或者是二面角，然后将输出文件重新保存为输入文件。如果输出文件不是你想要的，那就只能重新构建输入文件！

scf优化可以认为是对指定结构的波函数不断优化的过程，是为了找到某个指定结构下能量最低的波函数，而几何构型优化是对结构的优化的过程，是为了找到势能面上的某个能量极小结构

一次构型优化，需要多次scf优化才能最终完成

---