双金属协同催化过渡态计算总报错

nono114514

计算过渡态时，在计算hessian矩阵时就出现了如图的错误，然后就报错了。atom number 在180-260之间，会有一定概率出现这个情况，160以下没有遇到过。
计算水平如图，已经尝试的方法有：
1.不知道是不是def的自带赝势的问题，所以用gencep试了下ccpvdz，不太行，还是会有问题。
2.换用低水平，b3lyp/6-31g(d)/sdd，只有一个跑动了，1/6，换到def2svp guess=read 继续跑，目前正常。
3.虽然也没说convergence有问题但是我还是对scf动了心思，在思想家公社看了解决scf不收敛的文章然后试了试 scf=(novaracc,noincfock) int=acc2e=12 来搞，1/5的成功率，感觉关系不是很大..个人感觉是calcfc的时候的问题。
所有截至目前6个要算的任务只跑起来了两个..


另，由于双金属协同催化过渡态位阻比较大，我尝试过从产物出发进行优化然后拉长C-C键去opt=ts，但是我发现，上述方法没有搞定的分子，在产物优化的最后一步计算频率时也会报错！就是出现很多NaN导致out文件都无法直接用GW打开，而是要用记事本打开从第一个NaN直接删到结尾的内种..

求老师们解答！初学者，对这个错误产生的原理完全不懂，难道是G09的bug要换用G16？？

nono114514

sobereva 发表于 2022-4-26 08:02
可能是优化算法的bug，没有具体输出文件没法更进一步说

盲猜回答，可以尝试opt里加上cartesian，或者用r ...

昨天尝试了换用不同坐标系！cartesian可以解决一部分问题!多谢老师

nono114514

lao7 发表于 2022-4-25 16:16
你这优化到最后可能是结构无法改变了 试着改变初始构型  通过旋转 间距等

调整过了呜呜呜，卡在这儿1周了

sobereva

可能是优化算法的bug，没有具体输出文件没法更进一步说

盲猜回答，可以尝试opt里加上cartesian，或者用recalc等等试试
实在不行也可以尝试用ORCA优化

lao7

你这优化到最后可能是结构无法改变了 试着改变初始构型  通过旋转 间距等