计算化学公社

标题: 双金属协同催化过渡态计算总报错 [打印本页]

作者
Author: nono114514 时间: 2022-4-25 10:43
标题: 双金属协同催化过渡态计算总报错
计算过渡态时，在计算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？？

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

作者
Author: sobereva 时间: 2022-4-26 08:02
可能是优化算法的bug，没有具体输出文件没法更进一步说

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

作者
Author: nono114514 时间: 2022-4-26 14:13

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

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

作者
Author: nono114514 时间: 2022-4-26 14:14

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

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

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)