计算化学公社

标题: 走IRC出错 No lower point found -- run aborted. [打印本页]

作者
Author: 嘻嘻哈哈哈 时间: 2020-7-23 11:12
标题: 走IRC出错 No lower point found -- run aborted.
走IRC用#p IRC=(rcfc,reverse,noeigen,maxpoints=300,stepsize=8,LQA) ub3lyp/genecp em=gd3bj SCF=QC guess=read，走了74步后出现Search did not lower the energy significantly. No lower point found -- run aborted.

后来看有类似情况，说去掉noeigen就可以了，结果去掉之后还出现了这个问题。
问题：
1：出现这个错误该怎么处理？是过渡态找的不对吗？过渡态用的#p opt=(calcfc,ts,noeigen,cartesian,gdiis) freq ub3lyp/genecp scf=qc
em=gd3bj 混合基组使用的是def2svp 和sdd,有一个虚频-301，振动方向看着也符合。
2：之前走的IRC也添加了noeigen，没有出错，添加noeigen对于已经走成功的IRC来说，结果差别大吗？需要去掉noeigen再走一遍吗？
3：走forward方向时IRC开始曲线也不光滑，最后走出来的结构可以用作接下来的计算吗？

作者
Author: biogon 时间: 2020-7-23 15:01
可能是过渡态定位精度不够

作者
Author: 嘻嘻哈哈哈 时间: 2020-7-23 20:00

biogon 发表于 2020-7-23 15:01
可能是过渡态定位精度不够

好的，那已经走出来IRC的别的过渡态还需要重新提高精度吗？还是只把出错的这一个过程的过渡态增加精度就可以？

作者
Author: mutron 时间: 2020-7-23 21:33
本帖最后由 mutron 于 2020-7-23 21:35 编辑

IRC里用noeigen没意义，那是寻找过渡态时候不对每一步做Hessian矩阵本征值数目的检测
简谈Gaussian里找过渡态的关键词opt=TS和QST2、3
http://sobereva.com/460

你这个报错信息应该是scf（qc方法）不收敛。如果你的体系不是十分难收敛，不建议一上来就用QC，建议scf=xqc，甚至可以设置MaxConventionalCycles=64先用效率高的DIIS方法迭代64次才切换成二次收敛方法
解决SCF不收敛问题的方法
http://sobereva.com/61

如果你的过渡态是一直用qc方法找出来的，注意检查一下最后得到的波函数稳定性

如果嫌IRC不光滑，把步长再缩小到5，并考虑加上recalc=n

作者
Author: sobereva 时间: 2020-7-24 07:14
仔细看此文，该注意的这里都说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

绝对不要轻易用scf=QC

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