sobereva 发表于 2015-5-5 12:11
用此关键词，看看最后有无虚频、最后一次输出的四个收敛判断标准各是多少。
# M062X/6-311+G(d,p) int=ult ...

好的，我试试。谢谢Sob老师。

用此关键词，看看最后有无虚频、最后一次输出的四个收敛判断标准各是多少。
# M062X/6-311+G(d,p) int=ultrafine opt(tight,calcall) freq

补充一下哈，初猜结构是先用B3LYP优化过一次的，不过当时也是计算正常结束，结果中后两个判据是NO但是受力已小于判据100倍而被认为收敛了。这次我用的方法是# M062X/6-311+G(d,p) opt=calcall，结果还是后两个判据是NO但是受力已小于判据100倍而被认为收敛了，同时出现了虚频。我优化的是丙醇二聚体。谢谢sob老师。

1 这么说是错误的。可能只是你用的初猜或者优化方式不合理导致优化没完成，这是数值上的问题，不能从化学意义上说明什么。
并且要弄清楚到底是什么情况，是最后根本没收敛，还是说后两个判据是NO但是受力已小于判据100倍而被认为收敛了。

2 可能高斯根据calcall最后一步的Hessian自动进行了振动分析（因为没有额外耗时），此时的振动频率是没问题的。