计算化学公社

标题: IRC求助，未走到最小值就结束 [打印本页]

作者
Author: fresco604 时间: 2016-11-17 08:41
标题: IRC求助，未走到最小值就结束
(, 下载次数 Times of downloads: 20)

求助老师一个问题，我的IRC计算有一个方向（左边）每次都没有走到最小值就结束了，但是方向是正确的，离最小值就差几步了

我试过增加maxpoints / maxcycles 都没有用

这是我的input 文件：

# irc=(maxpoints=400,calcall,lqa,stepsize=40,maxcycle=800) b3lyp/6-31+g(d)
geom=connectivity

实际上只走了50步左右， stepsize设在40是因为其他stepsize得不到理想的结果，

右边已经基本上走到了我想要的结构左边离我想要的结构已经很接近了，但就是走不到最小值，所以斜率还很大。

求助各位老师，这种情况下该怎么办，谢谢！

作者
Author: sobereva 时间: 2016-11-17 08:58
maxcycles本来就不可能有丝毫用处，完全多余
看输出文件里什么原因导致停止。如果是提示发现已经走到了极小点而停止（程序这么认为的），尝试g03默认的GS2算法，以及不写LQA，以及小点的步长，各种情况都试
另外还要明确跑IRC目的何在，是否真有绝对必要跑到头

作者
Author: fresco604 时间: 2016-11-17 09:27

sobereva 发表于 2016-11-17 08:58
maxcycles本来就不可能有丝毫用处，完全多余
看输出文件里什么原因导致停止。如果是提示发现已经走到了极 ...

谢谢老师指教！

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