请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 15551|回复: 16

[Gaussian/gview] IRC 出错:Maximum number of corrector steps exceded

[复制链接]

29

帖子

0

威望

84

eV
积分
113

Level 2 能力者

发表于 2014-12-21 01:34:40 | 显示全部楼层 |阅读模式
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)        。  请教大家,该如何解决这个问题?

2万

帖子

25

威望

2万

eV
积分
46931

管理员

公社社长

发表于 2014-12-21 04:07:15 | 显示全部楼层
尝试其它的IRC生成算法,比如G03默认的GS2,或者LQA算法,见手册IRC部分的关键词介绍。或者减小IRC步长再试
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。约6000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

29

帖子

0

威望

84

eV
积分
113

Level 2 能力者

 楼主| 发表于 2014-12-21 16:01:56 | 显示全部楼层
用GS2算法正常结束了,  多谢sob 老师。

3

帖子

0

威望

44

eV
积分
47

Level 2 能力者

发表于 2015-6-30 21:33:54 | 显示全部楼层
顶一个,我也遇到这样的问题了!

250

帖子

0

威望

866

eV
积分
1116

Level 4 (黑子)

发表于 2015-7-8 08:48:25 | 显示全部楼层
mei 发表于 2014-12-21 16:01
用GS2算法正常结束了,  多谢sob 老师。

请问如何用GS2算法,可以发输入文件的关键字看看吗?

2万

帖子

25

威望

2万

eV
积分
46931

管理员

公社社长

发表于 2015-7-8 18:23:45 | 显示全部楼层
lei234 发表于 2015-7-8 08:48
请问如何用GS2算法,可以发输入文件的关键字看看吗?

IRC里写上GS2即可
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。约6000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

250

帖子

0

威望

866

eV
积分
1116

Level 4 (黑子)

发表于 2015-7-8 19:35:40 | 显示全部楼层
sobereva 发表于 2015-7-8 18:23
IRC里写上GS2即可

谢谢~

145

帖子

0

威望

910

eV
积分
1055

Level 4 (黑子)

发表于 2016-1-12 00:00:05 | 显示全部楼层
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谢谢

2万

帖子

25

威望

2万

eV
积分
46931

管理员

公社社长

发表于 2016-1-12 00:14:00 | 显示全部楼层
灰飞的旋律 发表于 2016-1-12 00:00
Sob君们和楼主,你们好。
我现在也碰到这个问题了呢。我是用TS找到的过渡态,虚频-500多且振动方向也是 ...


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

不要用什么IOp(1/7),毫无意义。maxcycle设这种东西几乎一定没用。gradientonly也根本不是解决这个的。你被误导了。
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。约6000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

145

帖子

0

威望

910

eV
积分
1055

Level 4 (黑子)

发表于 2016-1-12 00:26:02 | 显示全部楼层
sobereva 发表于 2016-1-12 00:14
最简单的做法,IRC里写上LQA。如果之后发现IRC不够平滑,用calcall或者减小步长

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

(⊙v⊙)嗯。非常谢谢您这么晚还指导我,我去试试。还有个问题也一直没解决,就是经常最后出现can't be archived而且还normal结束,IRC中有时一步不走或者走一两步,这种问题里一般虚频都较小-100多。这个问题我试过缩小步长,calcfc,convergence,然都并卵。还有其他解决方法吗。

2万

帖子

25

威望

2万

eV
积分
46931

管理员

公社社长

发表于 2016-1-12 02:05:55 | 显示全部楼层
灰飞的旋律 发表于 2016-1-12 00:26
(⊙v⊙)嗯。非常谢谢您这么晚还指导我,我去试试。还有个问题也一直没解决,就是经常最后出现can't be ar ...

这说明过渡态定位不够准确。
找过渡态时用opt=tight int=ultrafine(之后IRC时也得用ultrafine)
如果还不行,有两种做法:(1)找过渡态时用calcall,治本,但是耗时高 (2)加大IRC的步长,治标
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社QQ群,1号:18616395,2号:466017436。约6000人,专门交流理论、计算化学。两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

145

帖子

0

威望

910

eV
积分
1055

Level 4 (黑子)

发表于 2016-1-12 10:31:17 | 显示全部楼层
sobereva 发表于 2016-1-12 02:05
这说明过渡态定位不够准确。
找过渡态时用opt=tight int=ultrafine(之后IRC时也得用ultrafine)
如果 ...

恩,我也觉得是过渡态没找到理想的构型导致的。再次谢谢您啦

23

帖子

0

威望

222

eV
积分
245

Level 3 能力者

发表于 2016-1-14 23:01:54 | 显示全部楼层
我看google 上有提到过这个问题,有人说 use=l115 ?  如果 过渡态 曲线 过于 平坦。。

15

帖子

1

威望

200

eV
积分
235

Level 3 能力者

发表于 2016-6-17 13:52:25 | 显示全部楼层
use=l115应该和使用GS2是一个作用,都是03默认的

5

帖子

0

威望

42

eV
积分
47

Level 2 能力者

发表于 2019-7-12 07:59:17 | 显示全部楼层
sobereva 发表于 2014-12-21 04:07
尝试其它的IRC生成算法,比如G03默认的GS2,或者LQA算法,见手册IRC部分的关键词介绍。或者减小IRC步长再试

请问sob老师,出现上面错误可能是由什么造成的?一般如何避免上述错误呢?
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2019-10-21 09:13 , Processed in 0.189504 second(s), 23 queries .

快速回复 返回顶部 返回列表