计算化学公社

标题: IRC 出错：Maximum number of corrector steps exceded [打印本页]
作者
Author: mei    时间: 2014-12-21 01:34
标题: IRC 出错：Maximum number of corrector steps exceded
g09 算IRC时出错：“Maximum number of corrector steps exceded”  。log 文件的最后部分为：
“ Error in corrector energy =           0.0000000000
Magnitude of corrector gradient =     0.0001224981
Magnitude of analytic gradient =      0.0001224743
Magnitude of difference =             0.0000004525
Angle between gradients (degrees)=    0.2114
Pt  1 Step number1000 out of a maximum of1000
Modified Bulirsch-Stoer Extrapolation Cycles:
    EPS =    0.000010000000000
Maximum DWI energy   std dev =  0.000000310 at pt     2
Maximum DWI gradient std dev =  0.872506745 at pt    73
WARNING: Bulirsch-Stoer Method is not Converging
CORRECTOR INTEGRATION CONVERGENCE:
   Recorrection delta-x convergence threshold:    0.010000
   Delta-x Convergence NOT Met
Maximum number of corrector steps exceded.
Error termination via Lnk1e in /home/gauss/g09d/g09/l123.exe ”
                                                                     IRC输入为：  irc(calcfc,maxpoint=100,maxcyc=1000,forward)        。  请教大家，该如何解决这个问题？
作者
Author: sobereva    时间: 2014-12-21 04:07
尝试其它的IRC生成算法，比如G03默认的GS2，或者LQA算法，见手册IRC部分的关键词介绍。或者减小IRC步长再试
作者
Author: mei    时间: 2014-12-21 16:01
用GS2算法正常结束了，  多谢sob 老师。
作者
Author: ian    时间: 2015-6-30 21:33
顶一个，我也遇到这样的问题了！
作者
Author: lei234    时间: 2015-7-8 08:48
mei 发表于 2014-12-21 16:01
用GS2算法正常结束了，  多谢sob 老师。

请问如何用GS2算法，可以发输入文件的关键字看看吗？
作者
Author: sobereva    时间: 2015-7-8 18:23
lei234 发表于 2015-7-8 08:48
请问如何用GS2算法，可以发输入文件的关键字看看吗？

IRC里写上GS2即可
作者
Author: lei234    时间: 2015-7-8 19:35
sobereva 发表于 2015-7-8 18:23
IRC里写上GS2即可

谢谢~
作者
Author: 灰飞的旋律    时间: 2016-1-12 00:00
sobereva 发表于 2015-7-8 18:23
IRC里写上GS2即可

Sob君们和楼主，你们好。
我现在也碰到这个问题了呢。我是用TS找到的过渡态，虚频-500多且振动方向也是朝着生成产物的方向振动。我做IRC（RCFC）时，forward得到的是产物，reverse时就出现这个问题。试了好多参数都不行，首先是加上了maxcycle，又换成gradientonly，又加上了IOP(1/7=50),最后我看到您这个帖子回复，我就把gradientonly换成GS2，去掉1/7=50。正在计算中，但我在log文件中搜索L123发现，进入这个链接后输出的信息还是
(Enter /public/home/users/shzu001/g09/l123.exe)
IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC
Error in corrector energy =          -0.0000720839
Magnitude of corrector gradient =     0.0000390019
Magnitude of analytic gradient =      0.0004986517
Magnitude of difference =             0.0004956828
Angle between gradients (degrees)=   83.3962
Pt  1 Step number   6 out of a maximum of 400
Using modified Bulirsch-Stoer corrector integration.
   Using distance weighted interpolants to fit PES.
PreDWI: IReCr0= 2 IReCr1= 1 IReCr2= 1 IReCr3= 1
Modified Bulirsch-Stoer Extrapolation Cycles:
    EPS =    0.000100000000000
    Cycle     1 NS=    48
    Cycle     2 NS=    64 Truncation Error =  0.000237058
    Cycle     3 NS=    96 Truncation Error =  0.000000090
BS completed in    3 cycles and     208 integration steps.
Maximum DWI energy   std dev =  0.000114085 at pt    -1
Maximum DWI gradient std dev =  1.331170304 at pt    36
SUMMARY OF CORRECTOR INTEGRATION:
   Predictor End Point Energy = -2278.219096
   Old End Point Energy       = -2278.219078
   Corrected End Point Energy = -2278.219262
   Predictor End-Start Dist.  =     0.054017
   Old End-Start Dist.        =     0.000508
   New End-Start Dist.        =     0.045749
   New End-Old End Dist.      =     0.045337
CORRECTOR INTEGRATION CONVERGENCE:
   Recorrection delta-x convergence threshold:    0.000050
   Delta-x Convergence NOT Met
IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC
Leave Link  123 at Mon Jan 11 20:30:05 2016, MaxMem=  2621440000 cpu:         0.4
以前正常结束的IRC的log文件是
(Enter /public/home/users/shzu001/g09/l123.exe)
IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC
Pt  7 Step number   1 out of a maximum of 400
Using modified Bulirsch-Stoer corrector integration.
   Using distance weighted interpolants to fit PES.
PreDWI: IReCr0= 2 IReCr1= 1 IReCr2= 1 IReCr3= 1
Modified Bulirsch-Stoer Extrapolation Cycles:
    EPS =    0.000100000000000
    Cycle     1 NS=    48
    Cycle     2 NS=    64 Truncation Error =  0.000002098
BS completed in    2 cycles and     112 integration steps.
Maximum DWI energy   std dev =  0.000000890 at pt    47
Maximum DWI gradient std dev =  0.002961393 at pt    38
SUMMARY OF CORRECTOR INTEGRATION:
   Predictor End Point Energy = -2278.220653
   Old End Point Energy       = -2278.220653
   Corrected End Point Energy = -2278.220653
   Predictor End-Start Dist.  =     0.054017
   Old End-Start Dist.        =     0.054017
   New End-Start Dist.        =     0.054017
   New End-Old End Dist.      =     0.000274
CORRECTOR INTEGRATION CONVERGENCE:
   Recorrection delta-x convergence threshold:    0.000050
   Delta-x Convergence NOT Met
IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC-IRC
Leave Link  123 at Thu Jan  7 12:26:50 2016, MaxMem=  2621440000 cpu:         2.3
而且到现在已经经过好多次L123了，还是一个point都没找出来。。

非常期待您的回复，O(∩_∩)O谢谢
作者
Author: sobereva    时间: 2016-1-12 00:14
灰飞的旋律 发表于 2016-1-12 00:00
Sob君们和楼主，你们好。
我现在也碰到这个问题了呢。我是用TS找到的过渡态，虚频-500多且振动方向也是 ...


最简单的做法，IRC里写上LQA。如果之后发现IRC不够平滑，用calcall或者减小步长

不要用什么IOp(1/7)，毫无意义。maxcycle设这种东西几乎一定没用。gradientonly也根本不是解决这个的。你被误导了。
作者
Author: 灰飞的旋律    时间: 2016-1-12 00:26
sobereva 发表于 2016-1-12 00:14
最简单的做法，IRC里写上LQA。如果之后发现IRC不够平滑，用calcall或者减小步长

不要用什么IOp(1/7) ...

(⊙v⊙)嗯。非常谢谢您这么晚还指导我，我去试试。还有个问题也一直没解决，就是经常最后出现can't be archived而且还normal结束，IRC中有时一步不走或者走一两步，这种问题里一般虚频都较小-100多。这个问题我试过缩小步长，calcfc，convergence，然都并卵。还有其他解决方法吗。
作者
Author: sobereva    时间: 2016-1-12 02:05
灰飞的旋律 发表于 2016-1-12 00:26
(⊙v⊙)嗯。非常谢谢您这么晚还指导我，我去试试。还有个问题也一直没解决，就是经常最后出现can't be ar ...

这说明过渡态定位不够准确。
找过渡态时用opt=tight int=ultrafine（之后IRC时也得用ultrafine）
如果还不行，有两种做法：(1)找过渡态时用calcall，治本，但是耗时高 (2)加大IRC的步长，治标
作者
Author: 灰飞的旋律    时间: 2016-1-12 10:31
sobereva 发表于 2016-1-12 02:05
这说明过渡态定位不够准确。
找过渡态时用opt=tight int=ultrafine（之后IRC时也得用ultrafine）
如果 ...

恩，我也觉得是过渡态没找到理想的构型导致的。再次谢谢您啦
作者
Author: ximi1986    时间: 2016-1-14 23:01
我看google 上有提到过这个问题，有人说 use=l115 ?  如果 过渡态 曲线 过于 平坦。。
作者
Author: Merlot    时间: 2016-6-17 13:52
use=l115应该和使用GS2是一个作用，都是03默认的
作者
Author: 1316615557    时间: 2019-7-12 07:59
sobereva 发表于 2014-12-21 04:07
尝试其它的IRC生成算法，比如G03默认的GS2，或者LQA算法，见手册IRC部分的关键词介绍。或者减小IRC步长再试

请问sob老师，出现上面错误可能是由什么造成的？一般如何避免上述错误呢？
作者
Author: sobereva    时间: 2019-7-13 05:09
1316615557 发表于 2019-7-12 07:59
请问sob老师，出现上面错误可能是由什么造成的？一般如何避免上述错误呢？

在Gaussian中计算IRC的方法和常见问题
http://sobereva.com/400（http://bbs.keinsci.com/thread-7736-1-1.html）
作者
Author: 1316615557    时间: 2019-7-17 09:37
好的，谢谢sob老师
作者
Author: syuusama    时间: 2024-7-12 11:39
1316615557 发表于 2019-7-12 07:59
请问sob老师，出现上面错误可能是由什么造成的？一般如何避免上述错误呢？

出现在这个问题的原因是HPC方法校正步不收敛，http://sobereva.com/400，sob老师在里面有详细的原因和解决方法



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3