计算化学公社

标题: IRC问题求助，只跑出两个点任务结束 [打印本页]

作者
Author: leo1992 时间: 2021-3-23 21:49
标题: IRC问题求助，只跑出两个点任务结束
本帖最后由 leo1992 于 2021-3-23 22:56 编辑

（已经认真阅读“在Gaussian中计算IRC的方法和常见问题http://sobereva.com/400”帖子并尝试解决,无果，特来求助）。在对过渡态跑IRC过程中遇到问题描述如下：过渡态结构如附件TS文件夹中4，采用B3LYP/6-311G** opt=(calcall,ts,noeigen,tight) freq int=superfine 寻找过渡态和振动分析。过渡态虚频振动方向与反应物产物趋势一致，是构象互变过渡态计算。（之前还采用过B3LYP/6-31G* opt=(calcall,ts,noeigen) freq进行过渡态计算，见文件夹TS中1和2信息。但IRC出现同样下述问题，遂改用当前关键词计算过渡态）

在对过渡态进行IRC。

第一次计算：文件夹IRC中meso22，关键词B3LYP/6-311G** IRC(calcall,maxpoint=20) int=superfine guess=read，报错信息Maximumnumber of corrector steps exceded.根据帖子内容改用GS2或LQA

第二次计算：文件夹IRC中meso23，关键词B3LYP/6-311G** IRC(calcall,maxpoint=20,GS2) int=superfine guess=read，输出文件最后Thistype of calculation cannot be archived.IRC路径上只有2个点，且两个结构有重合，见meso23.log。

考虑到可能因为过渡态与反应物产物能量相差较小，计算过渡态与反应物和产物电子能相差10kJ/mol，势能面变化较小，遂改变各种步长进行计算。

第345次计算：尝试使用步长分别为stepsize=20，30，60，也尝试过3，5）见IRC文件家中文件meso24-26。计算结果与第二次结果一致，均只有两个点。

同时也尝试过改用IRC中采用LQA计算，以及采用M062X寻找过渡态TS，跑IRC，结果均与上面一致。请问是何种原因导致上述问题，以及该如何解决。

作者
Author: 星斗如盘 时间: 2021-3-24 22:03
使用LQA算法一般没问题，把stepsize设为20或者40就可以了

作者
Author: leo1992 时间: 2021-3-25 09:18

星斗如盘发表于 2021-3-24 22:03
使用LQA算法一般没问题，把stepsize设为20或者40就可以了

问题已解决，感谢

作者
Author: 丛炜熹 时间: 2023-3-23 19:34

leo1992 发表于 2021-3-25 09:18
问题已解决，感谢

也遇到了相同问题，想请教一下前辈是怎么解决的？

作者
Author: ThualishuT 时间: 2023-3-24 00:46

丛炜熹发表于 2023-3-23 19:34
也遇到了相同问题，想请教一下前辈是怎么解决的？

一要确保过渡态定位足够精确。二如果过渡态虚频较小，优先尝试IRC步长调大，很可能是势能面过于平缓导致的。

作者
Author: 丛炜熹 时间: 2023-3-25 17:23

ThualishuT 发表于 2023-3-24 00:46
一要确保过渡态定位足够精确。二如果过渡态虚频较小，优先尝试IRC步长调大，很可能是势能面过于平缓导致 ...

好的，非常感谢前辈的回复，我按前辈的方法再试一下。麻烦前辈了

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