震荡严重，优化不出来结构，是因为初始结构不合理吗？如果直接用这个结构跑freq，虚频比较多，但是可以得到想要虚频，该怎么解决？求解答

从你这个电子能量来看，分子也不是很大。
那过渡态这个样子大概率是初始结构不够好导致的。要不可能就是该加溶剂模型的没加或者不该加溶剂模型的加了。之类的问题

科研小宝 发表于 2020-5-11 14:52
从你这个电子能量来看，分子也不是很大。
那过渡态这个样子大概率是初始结构不够好导致的。要不可能就是该 ...

您好！是加了溶剂模型，溶剂会影响优化和过渡态结果吗？

疯不觉 发表于 2020-5-11 14:54
您好！是加了溶剂模型，溶剂会影响优化和过渡态结果吗？

溶剂模型肯定会有影响。不过90%以上的过渡态不收敛都是因为初始结构不够好导致的。 可以尝试调整初始结构。或者你先把这个过渡态两侧的中间体找到。根据中间体结构再去猜测过渡态结构。

科研小宝 发表于 2020-5-11 14:55
溶剂模型肯定会有影响。不过90%以上的过渡态不收敛都是因为初始结构不够好导致的。 可以尝试调整初始结构 ...

老师，不是过渡态振荡，而是前一步优化就出现了振荡，气相中优化正常，加了溶剂模型，优化不出来。

疯不觉 发表于 2020-5-11 14:59
老师，不是过渡态振荡，而是前一步优化就出现了振荡，气相中优化正常，加了溶剂模型，优化不出来。
...

这个状况也不少见。可以试试别的溶剂模型看看。

科研小宝 发表于 2020-5-11 14:59
这个状况也不少见。可以试试别的溶剂模型看看。

老师，之前的计算都是用的这个模型！如果就想用这个模型优化，还有没有其他办法？出现这种情况还是因为结果不合理吗？

疯不觉 发表于 2020-5-11 15:00
老师，之前的计算都是用的这个模型！如果就想用这个模型优化，还有没有其他办法？出现这种情况还是因为结 ...

对付几何优化不收敛。而你不想换溶剂模型那就只能用常规手段来一套。

1        OPT=GDIIS (或者OPT=RFO)

高斯的默认优化算法是OPT=GEDIIS ，在优化比较“soft” 比较巨大的体系时候可能不好用。没准换成别的优化算法可以帮助收敛。
2        OPT=CALCALL (或者 OPT=RECALC=3)
OPT=CACALL这个关键词帮助高斯在优化的步骤中进行严格的HESSIAN矩阵计算。这样会帮助收敛。但是慎用，对于原子很多的体系每一步都计算HESSIAN其实就是每一步都算一次频率一样。会让计算总时长翻十倍都不止可能。只有在非常难收敛的时候再考虑使用。OPT=RECALC=3 意思是每三步计算一次HESSIAN，会比CALCALL快很多。G16中支持。
3        OPT=MAXCYC=200
这个就是增加优化步数的上限。这个仅仅对于原子数小于10个的体系使用才有效。因为原子数目比较少的时候，高斯的默认步数太小了 （最大步数=2倍的冗余内坐标数目）
4        OPT(MAXSTEP=10, NOTRUST)   (默认值是30 * 0.01Bohr)
这个关键词是缩短优化过程每一步的最大步长。对于优化中反复出现震荡的体系可能会有效。（查看优化进程是否出现振荡，查看本文GaussView使用部分3.1）NOTRUST的意思是让优化中每一步的步长上限都是10.不写这个的话从第二个点开始步长又回到30.
5        OPT=LOOSE
对于有的很大的体系，就是很难收敛这时候可以释放放宽收敛限。但是如果要进行频率计算。则这个做法无效。因为频率需要很准的结构，一般至少是OPT=TIGHT 甚至要OPT=VERYTIGHT才行
6    以上提到的方法可以各种组合，帮助收敛。比如
        OPT(MAXSTEP=5, NOTRUST, RECALC=3,GDIIS)
  而且如果分子又高对称性也最好打破对称性去算。要不然也容易振荡。