本帖最后由 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,结果均与上面一致。请问是何种原因导致上述问题,以及该如何解决。
|