计算化学公社

标题: TS计算正常，但IRC只跑产物一边，这是什么原因呢？ [打印本页]

作者
Author: jun_UuU 时间: 2021-4-2 16:03
标题: TS计算正常，但IRC只跑产物一边，这是什么原因呢？
本帖最后由 jun_UuU 于 2021-4-2 16:22 编辑

已经认真看过IRC的常见问题博文了，但还是不明白为什么，只能来求助了，拜托！ http://bbs.keinsci.com/thread-7736-1-1.html
用linux的g16跑的，TS结果收敛，虚频正确，但IRC重复多次，每次都只能跑产物的一边，反应物只跑出来两次，一次是ReCorrect=never来试试，另一次跑出来了但任务仍然显示错误。
总结：
TS的关键词：# Opt=(CalcFC,TS,Noeigen) Freq B3LYP/6-31+G(d) SCRF=(SMD,Solvent=Generic,Read) EM=GD3BJ
IRC_1的关键词：# IRC=(Calcfc,maxpoints=300,stepsize=5) B3LYP/6-31+G(d) SCRF=(SMD,Solvent=Generic,Read) EM=GD3BJ
IRC_2是# IRC=(Calcfc,recorrect=never) （后面省略，没有改动）
IRC_3是# IRC=(Calcfc,maxpoints=200,stepsize=6)
IRC_4是# IRC=(Calcfc,maxpoints=200,stepsize=3)
IRC_5是# IRC=(Calcfc,maxpoints=200,stepsize=3,GS2)
IRC_6是# IRC=(Calcfc,recalc=3,maxpoints=200,stepsize=3)
结果：
1、3、4、5报l123错误，6报l502错误，2是没有矫正步，正确但精度太低了。附件中有1、2、6的IRC结果图。

随后对TS重新tight了一下，得到的结果如下：
IRC_new使用关键词：# IRC=(Calcfc,maxpoints=200,stepsize=3) B3LYP/6-31+G(d) SCRF=(SMD,Solvent=Generic,Read) EM=GD3BJ
这个图应该是可以的，看了一下.out文件中有PES minimum detected。.out文件太大了传不了，虽然方向是反的，但是仍然报错l123。请问是什么原因呢？
与我设置的较大maxpoint有关系吗？

溶剂模型是环氧丙烷，数据如下：
SolventName=PO
EpsInf=1.8670
Eps=20.1046

作者
Author: jun_UuU 时间: 2021-4-2 16:12
本帖最后由 jun_UuU 于 2021-4-2 16:14 编辑

同样的在另一个IRC计算中，课题组同学的关键词设为# IRC=(Calcfc,maxpoints=120,stepsize=2) 跑出来的图非常漂亮，
但是我用同样的坐标，方法，只是改变maxpoint和stepsize就是跑不出来，只有产物的一边，
IRC_1是# IRC=(Calcfc,maxpoints=200,stepsize=3)
IRC_2是# IRC=(Calcfc,maxpoints=200,stepsize=4)
IRC_3是# IRC=(Calcfc,maxpoints=200,stepsize=2) （这个没跑出来真的搞不懂，产物的IRC能量已经是挺平的一条直线了，但还是错误。）
报的错误都是l123， Maximum number of corrector steps exceded.
接下来是不是要降低maxpoint来算呀？但这其中是为什么呢？

作者
Author: 喵星大佬 时间: 2021-4-3 00:24

jun_UuU 发表于 2021-4-2 16:12
同样的在另一个IRC计算中，课题组同学的关键词设为# IRC=(Calcfc,maxpoints=120,stepsize=2) 跑出来的图非 ...

写了是校正步不收敛啊，你可以设为never或者用LQA方法

作者
Author: sobereva 时间: 2021-4-3 16:07
博文都已经明确说了，不建议用ReCorrect=never，要用就用LQA，明显没仔细follow博文

跑IRC不成功的时候，如果用了SMD，把opt=TS和IRC任务都改成默认的IEFPCM，博文里也已经明确说了

用SMD溶剂模型时，也可能个别时候由于数值噪音问题出现IRC走几步就停了的现象。可将优化和走IRC用的溶剂模型都改为IEFPCM再试，说不定能解决。

再认真一个字一个字阅读博文，不要试没意义的关键词而忽略了该考虑的解决办法

作者
Author: jun_UuU 时间: 2021-4-8 20:31
好的，谢谢因为之前说LQA精度不够，就先尝试的其他几个方法，不过现在已经用LQA算好啦，结果还不错。

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