gaussian基本概念和用法.doc

《gaussian基本概念和用法.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、1计算流程上面的图应该是从L301说起的，即图中的“基组”，L301的作用是产生基组信息。（在这之前还有L1：处理计算执行路径，创建执行链接的列表，并初始化scratch文件；L101：读取标题和分子说明部分；L103：berny优化到最小值；L202：重新定位坐标，计算对称性，检查变量）“基组”之后是“分子结构”即L401，形成初始的分子轨道初猜。第一个框中有“分子轨道”和“能量”两项，我想这应该指的是L502即迭代求解SCF方程。下面红线圈的框中则是迭代求解的具体过程。每次迭代完毕都得到一个能量。即（SCFDone: E(RH

2、F)= -xx.xxxx  A.U.after xxcycles）然后是L103判断力和位移是否满足收敛条件。如果满足（即四个yes），则进行布居分析等，并完成计算。如果不满足，则继续调整分子结构（L401），再次进行迭代求解，反复这个过程，直到满足力和位移收敛条件为止。如果不是HF方法而是多体微扰，CI等后-scf方法的话，在每次迭代（L502）完成后还要多一步相关能的计算，以MP2为例，迭代完成后多出:L801(双电子积分变换的初始化),L906(半直接的MP2)和L1002(迭代求解CPHF方程)等计算。对于楼主的图，我觉得

3、力和位移判断为NO以后，箭头应该指向最头上的“基组”，因为实际过程中，每次构型调整都会重新定位坐标和产生基组信息。常见问题分析1．检查是否有初始文件错误在命令行中加入%kJobL301or%kJobL302如果通过则一般初始文件ok。常见初级错误：a.自旋多重度错误;b.变量赋值为整数;c.变量没有赋值或多重赋值;d.键角小于等于0度，大于等于180度e.分子描述后面没有空行;f.二面角判断错误，造成两个原子距离过近;g.分子描述一行内两次参考同一原子,或参考原子共线2．SCF(自洽场)不收敛则一般是L502错误,省却情况做64个

4、cycle迭代(G03缺省128cycles)a.修改坐标，使之合理;b.改变初始猜 Guess=Huckel或其他的，看Guess关键词;c.增加叠代次数SCFCYC=N1计算流程上面的图应该是从L301说起的，即图中的“基组”，L301的作用是产生基组信息。（在这之前还有L1：处理计算执行路径，创建执行链接的列表，并初始化scratch文件；L101：读取标题和分子说明部分；L103：berny优化到最小值；L202：重新定位坐标，计算对称性，检查变量）“基组”之后是“分子结构”即L401，形成初始的分子轨道初猜。第一个框中有

5、“分子轨道”和“能量”两项，我想这应该指的是L502即迭代求解SCF方程。下面红线圈的框中则是迭代求解的具体过程。每次迭代完毕都得到一个能量。即（SCFDone: E(RHF)= -xx.xxxx  A.U.after xxcycles）然后是L103判断力和位移是否满足收敛条件。如果满足（即四个yes），则进行布居分析等，并完成计算。如果不满足，则继续调整分子结构（L401），再次进行迭代求解，反复这个过程，直到满足力和位移收敛条件为止。如果不是HF方法而是多体微扰，CI等后-scf方法的话，在每次迭代（L502）完成后还要多一

6、步相关能的计算，以MP2为例，迭代完成后多出:L801(双电子积分变换的初始化),L906(半直接的MP2)和L1002(迭代求解CPHF方程)等计算。对于楼主的图，我觉得力和位移判断为NO以后，箭头应该指向最头上的“基组”，因为实际过程中，每次构型调整都会重新定位坐标和产生基组信息。常见问题分析1．检查是否有初始文件错误在命令行中加入%kJobL301or%kJobL302如果通过则一般初始文件ok。常见初级错误：a.自旋多重度错误;b.变量赋值为整数;c.变量没有赋值或多重赋值;d.键角小于等于0度，大于等于180度e.分子描

7、述后面没有空行;f.二面角判断错误，造成两个原子距离过近;g.分子描述一行内两次参考同一原子,或参考原子共线2．SCF(自洽场)不收敛则一般是L502错误,省却情况做64个cycle迭代(G03缺省128cycles)a.修改坐标，使之合理;b.改变初始猜 Guess=Huckel或其他的，看Guess关键词;c.增加叠代次数SCFCYC=N（对小分子作计算时最好不要增加，很可能结构不合理）;d.iop（5/13=1）这样忽略不收敛，继续往下做。3．分子对称性改变a.修改坐标，强制高对称性或放松对称性;b.给出精确的、对称性确定的

8、角度和二面角。如CH4的角度给到109.47122;c.放松对称性判据 Symm=loose; d.不做对称性检查iop(2/16=1),最好加这个选项）,iop(2/16=2)则保持新的对称性来计算4．Opt时收敛的问题a.修改坐标，使之合理;b