如何阅读高斯的计算输出文件
高斯输出文件的关键是看Li
k每个Li
k对应不同的部分以单点计算为例单点计算的流程见附件输入行在后加入P不然的话就只会显示Li
k1和0输出文件一般很长。对于单点计算，只有300多行。第一部分：基本运行信息输入输出文件名，初始命令进入Li
k1，显示进程号首先是版权说明然后是作者利用guassia
计算所发的文章，参考文献上必须列的进入Li
k101读入输入文件，将控制命令转化为程序能够看得懂的IOP占位段的选项和应用设置LeaveLi
k101进入Li
k202：判断体系对称性，并决定实际计算中对称性应该如何利用判断体系点群例如StoichiometryCH2OFrameworkgroupC2VC2COSGVH2Degoffreedom3Fullpoi
tgroupC2VNOp4LargestAbelia
subgroupC2VNOp4Largestco
ciseAbelia
subgroupC2NOp2再将输入的分子坐标转换为标准内坐标就是Sta
dardorie
tatio
这是程序默认，充分利用分子的对称性来达到方便计算的目的。就是将体系的质量中心放在坐标轴的主轴或者原点上。可以用
osym来禁止这个操作LeaveLi
k202进入Li
k301产生基组信息基组函数对称性
f计算出核排斥能LeaveLi
k301进入Li
k302，303这部分是计算积分的具体算法一般不会列出详细过程一句话带过leaveLi
k302303进入Li
k401在实际计算之前，必须有初始猜测，这部分的功能就是产生初始猜测，可以从chk文件读入，也可以程序自动产生。这部分可以用guess关键字来指定。我这个体系程序默认是ProjectedINDOGuess请仔细看手册不同的体系Guess是不一样的LeaveLi
k401进入Li
k502（主角登场！）自恰迭代求解SCF方程设定收敛标准：除了对包含比Ar重的原子的分子全电子（非ECP）计算外，积分只计算到106的精度。λSCF计算的能量和电子密度都收敛到104，或只有能量收敛到105为止，无论哪一个先达到。λ给出每一步迭代的信息（？）这只有在程序输入中加入p才会显示如果成功计算完毕SCFDo
e接下来就是打印计算出来的能量自旋维力值和收敛判据了LeaveLi
k502一般计算经常会在这一步挂掉哈哈如果你的单点计算只写了这些输出结束如果你加入了populatio
的内容那么基于scf计算的结果布局分析开始进入Li
k601Li
k601利用优化后的波函数，进行Mullike
populatio
得到轨道对称性电子态各个轨道的本征值然后是一些电荷自旋密度的分析leavLi
k601最后是整个计算结果的一个总结各小节之间用分开
f单点计算看明白之后看几何优化的就比较清楚了优化计算的第一步必须和以后各步骤分开处理，因为有几个操作只需做一次。例如r