高斯常见错误及注意地方.docVIP

优化1在优化时采用Scf(tight)的选项，增加收敛的标准。再去计算频率。如果还有虚频，参见下一步。2. 对称性的影响，很多情况下的虚频是由于分子本身的对称性造成的。这样，在优化时，如果必要，要将对称性降低，还有，输入文件有时是用内坐标。建议如果有虚频的话，将内坐标改成直角坐标优化。3.如果上述方法还有虚频，看一下虚频，找到强度较大的，将在频率中产生的原子的振动坐标加到相应的输入文件中。这样，重新计算。直到虚频没有。4 实际上，如果分子柔性较大，很难找到最低点，这是电子结构计算的问题，这种情况下，需要动力学的 东西，用构象搜寻的办法解决。如：模拟退火，最陡下降法，淬火法等。将得到的能量最低的构象做一般的电子结构计算，这样，应当没有问题。不要讲你还没有得到最稳定的结构，那么，是你的分子有问题，要么计算错了，要么就是游离出现代计算的范畴 1。低频率振动模式时，力常数很小，必须使用opt＝tight以确保适当的收敛和随后任务步骤中频率计算的可靠度。比如我曾经算过一个体系，第一个频率振动的力常数只有10～20，很不放心，加入tight后，优出来的结构却与前面大相径庭，而且还出现了一个虚频。2。做IRC计算的时候，如果步长很小，必须要用verytight。就曾有人步长设为1，而不加verytight，算出来的反应路径是v字形倒海鸥状的。3。使用int指定用于数值积分的积分网格来用于消虚频的时候，必须注意，在比较能量的时候，对所有的计算要使用相同的积分网格。 检查错误 1． 检查是否有初始文件错误在命令行中加入 %kJob L301 or %kJob L302如果通过则一般初始文件ok。常见初级错误：a. 自旋多重度错误b. 变量赋值为整数c. 变量没有赋值或多重赋值d. 键角小于等于0度，大于等于180度e. 分子描述后面没有空行f. 二面角判断错误，造成两个原子距离过近g. 分子描述一行内两次参考同一原子,或参考原子共线2． SCF(自洽场)不收敛则一般是L502错误省却情况做64个cycle迭代(G03缺省128 cycles)a. 修改坐标，使之合理b. 改变初始猜 ?Guess=Huckel 或其他的，看Guess关键词。c. 增加叠代次数 SCFCYC=N （对小分子作计算时最好不要增加，很可能结构不合理）d. iop（5/13=1）这样忽略不收敛，继续往下做。3． 分子对称性改变a. 修改坐标，强制高对称性或放松对称性b. 给出精确的、对称性确定的角度和二面角。 如CH4的角度给到109.47122c. 放松对称性判据 ?Symm=loose? d. 不做对称性检查 iop(2/16=1) （最好加这个选项）? ? ?iop(2/16=2) 则保持新的对称性来计算4． Opt时收敛的问题a. 修改坐标，使之合理b. 增加叠代次数optcyc=N5． 优化过渡态，若势能面太平缓，则不好找到。iop(1/8=10) 默认30（下一个结构和该结构的差别0.3?），可改成10。如果每一步都要用到小的步长，应该加opt(notrustupdate)6． 在CI（组态）方法中如QCISD(T),CCSD(T),CID方法中，省却最大循环50，若出错（L913错误）解决方法：#P QCISD(maxcyc=N) ? ?注：N≤5127． 优化过渡态opt=TS (给出过渡态) ?opt=qst2 (给出反应物和产物) ?opt=qst3 (给出反应物和产物和过渡态)a. ?用G03时的出错 ? opt=ts 必须加FC (force constant)写法：opt=(TS, calcFc)or ? ?opt=(TS,calchffc)计算HF力常数，对QCISD，CCSD等方法用；or ? ?opt=(TS,modRedundant) ?(最好写这个)b. 如果计算采用QCISD计算（不好计算FC）则写为QCISD opt=(TS, calcHFFC) (用HF计算FC)8. 无法写大的Scratch文件RWFa. 劈裂RWF文件 %rwf=loc1,size1,loc2,size2,……..,locN,-1? ?b. 改变计算方法 MP2=Direct可以少占硬盘空间? ?c. 限制最大硬盘 maxdisk=N GB，****MB,有些系统写2GB会出错，可以写2000MB9. FOPT出错 原因是变量数与分子自由度数不相等。 可用POPT 或直接用OPT10. 优化过渡态只能做一个STEP 原因是负本征数目不对 添加 iop（1/11）=1或者noeigentest过度态优化 1.首先