计算化学公社

标题: IRC的两个方向跑的步数不一致，只走几个点 [打印本页]

作者
Author: 茶新味 时间: 2017-9-13 17:30
标题: IRC的两个方向跑的步数不一致，只走几个点
本帖最后由茶新味于 2017-9-13 17:32 编辑

Sob老师，您好。我找了过渡态，所用的关键字是# mp2/aug-cc-pvdz opt(ts,tight,calcall) int=ultrafine freq，然后分别往两边走IRC，
分别设置如下关键字
#mp2/aug-cc-pvdz irc=(forward,calcall,maxpoints=100,stepsize=10,LQA) int=ultrafine
# mp2/aug-cc-pvdz irc=(reverse,calcall,maxpoints=100,stepsize=10,LQA) int=ultrafine
其中stepsize分别用10,20,30代替，可是当往reverse方向走的，都只走几步，如图reverse
往前走最多也就是走49个点，如图forward，图中红色字体代表步长。不同步数的文件都能正常结束。
我想用精确的IRC证明过渡态，请问老师有什么办法跑完给定的步数呢？

作者
Author: sobereva 时间: 2017-9-13 17:59
int=ultrafine只影响DFT结果，对MP2设它无用

用了calcall通常就没有必要再用LQA。用LQA走步不如默认的算法准，曲线容易波动。除非calcall依然有提示重校正步不收敛的情况才需要加LQA

reverse设10的时候不是已经走了好多步了么。

用了LQA，或者步长太大，会导致IRC走的不准，可能导致提前判断为极小点而自发终止任务。

总之，建议先去了LQA，步长就用默认的。

作者
Author: 茶新味 时间: 2017-9-14 10:18
Sob老师，您好。我按照你的指导走IRC关键字为：# mp2/aug-cc-pvdz irc=(reverse,calcfc,maxpoints=100)，可是只走了两个点，并且这两个点还垂直。如图一，out文件出现了图二错误。
查了李博士的文章Gaussian FQA文章，提到成因：Gaussian的默认IRC算法HPC为了使曲线平滑需做重校正，但常有校正步不收敛的问题。换用GS2算法时，其限制性优化有时也不收敛。解决：增加IRC=LQA选项，例如IRC=(calcfc,LQA)。
于是又改变关键字从新尝试：# mp2/aug-cc-pvdz irc=(reverse,calcfc,maxpoints=100,LQA) 。正常结束如图三，但是只走了5步。Sob老师，是我的方法基组不合理吗？我研究的弱相互作用。

作者
Author: sobereva 时间: 2017-9-14 10:48

茶新味发表于 2017-9-14 10:18
Sob老师，您好。我按照你的指导走IRC关键字为：# mp2/aug-cc-pvdz irc=(reverse,calcfc,maxpoints=10 ...

我没让你用 mp2/aug-cc-pvdz irc=(reverse,calcfc,maxpoints=100)，怎么把calcall给去了？
跟基组完全无关

作者
Author: 茶新味 时间: 2017-9-14 11:27
Sob老师，去掉LQA，用# mp2/aug-cc-pvdz irc=(reverse,calcall,maxpoints=100)，也只走两步，out文件，出现下列错误。

作者
Author: sobereva 时间: 2017-9-14 12:19

茶新味发表于 2017-9-14 11:27
Sob老师，去掉LQA，用# mp2/aug-cc-pvdz irc=(reverse,calcall,maxpoints=100)，也只走两步，out文件，出 ...

在有calcall的前提下才去掉LQA。你没有正确理解我的文字。
没有calcall时不用LQA经常重校正步不收敛，而有calcall时出现这种问题的概率较低，没必要再用LQA，此时再用LQA只会使得IRC精度变差，甚至判断错是否已经到了极小点而应当自发终止

作者
Author: keyanxu666 时间: 2020-8-13 15:25
请问你这个问题解决了吗，我也遇到了这种问题

作者
Author: sobereva 时间: 2020-8-17 09:13

keyanxu666 发表于 2020-8-13 15:25
请问你这个问题解决了吗，我也遇到了这种问题

该说的这里都说了
在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）

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