24小时热门版块排行榜    

查看: 2175  |  回复: 9

Meng xf

新虫 (小有名气)


[求助] 高斯跑 irc出现123错误,有人知道怎么解决吗? 已有1人参与

高斯跑 irc出现123错误,有人知道怎么解决吗?谢谢
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : irc6-01.log
  • 2019-11-27 20:48:46, 573.61 K

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

感谢参与,应助指数 +1
是因为IRC的默认算法HPC在计算每一个点的时候都要做一个校正步。要进行迭代。当迭代到最大次数时候如果还不收敛 就会报错。这就是为什么末尾提示 corrector steps exceed的原因。比较简单的方法就是改变irc的算法 换成 LQA算法  关键词写成 irc=(calcfc,maxpoints=20,LQA) 就可以了。

» 本帖已获得的红花(最新10朵)

2楼2019-11-27 22:00:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

建建爱乖乖

金虫 (小有名气)

IRC(calcfc,maxpoints=20, lqa),还有建议优化过渡态时如果使用opt(ts, calcfc),IRC就用calcfc。如果用的是opt(ts, calcall),IRC强烈建议用calcall。
3楼2019-11-28 09:23:38
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Meng xf

新虫 (小有名气)


引用回帖:
2楼: Originally posted by paramecium86 at 2019-11-27 22:00:41
是因为IRC的默认算法HPC在计算每一个点的时候都要做一个校正步。要进行迭代。当迭代到最大次数时候如果还不收敛 就会报错。这就是为什么末尾提示 corrector steps exceed的原因。比较简单的方法就是改变irc的算法 换 ...

谢谢,解决了

发自小木虫Android客户端
4楼2019-11-30 20:13:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Meng xf

新虫 (小有名气)


引用回帖:
4楼: Originally posted by Meng xf at 2019-11-30 20:13:37
谢谢,解决了
...

要算一下键能,请问扭转能怎么算?

发自小木虫Android客户端
5楼2019-11-30 20:14:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Meng xf

新虫 (小有名气)


引用回帖:
2楼: Originally posted by paramecium86 at 2019-11-27 22:00:41
是因为IRC的默认算法HPC在计算每一个点的时候都要做一个校正步。要进行迭代。当迭代到最大次数时候如果还不收敛 就会报错。这就是为什么末尾提示 corrector steps exceed的原因。比较简单的方法就是改变irc的算法 换 ...

谢谢,就是用了这个方法,解决了

发自小木虫Android客户端
6楼2019-11-30 20:15:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Meng xf

新虫 (小有名气)


送红花一朵
引用回帖:
2楼: Originally posted by paramecium86 at 2019-11-27 22:00:41
是因为IRC的默认算法HPC在计算每一个点的时候都要做一个校正步。要进行迭代。当迭代到最大次数时候如果还不收敛 就会报错。这就是为什么末尾提示 corrector steps exceed的原因。比较简单的方法就是改变irc的算法 换 ...

请问扭转能怎么算?

发自小木虫Android客户端
7楼2019-11-30 20:16:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
Meng xf: 金币+20, ★★★★★最佳答案 2019-12-02 15:54:20
引用回帖:
7楼: Originally posted by Meng xf at 2019-11-30 07:16:12
请问扭转能怎么算?
...

计算某一个二面角的扭转能 那就沿着这个二面角做柔性扫描。 让这个二面角转动一整圈 扫描出来的势能曲线的能垒就是扭转能
8楼2019-12-01 00:15:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Meng xf

新虫 (小有名气)


9楼2019-12-02 15:54:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

LiuEve

新虫 (初入文坛)

在计算化学公社里看到的
Q:IRC任务报错啦!末尾提示Maximum number of corrector steps exceded咋办?
A:这就是前面反复说的HPC方法校正步不收敛。有以下办法:
(1)用LQA,比如IRC(calcfc,LQA),由于此时不涉及校正步了,因此100%解决问题。但这容易造成IRC不够准确、不光滑。
(2)如果你不希望用LQA在避免报错的同时牺牲IRC精度,则尝试减小步长(越小避免报错的几率越高),比如用IRC(calcfc,stepsize=5)。或者加上calcall,若嫌太昂贵就改用recalc=x(x越小避免报错的几率越高,但也越昂贵)。
(3)改用GS2算法,即IRC(calcfc,GS2),可完全避免以上报错。耗时比LQA高得多,但精度也比LQA好。
(4)IRC里加上ReCorrect=never,这使得HPC方法不做校正步,故也完全避免以上报错。此时耗时和LQA相同,但所得IRC精度不如LQA,因此强烈不建议用。
(5)IRC里加上maxcyc=N(N应大于默认的20)来加大HPC校正步迭代次数上限。笔者时常在菜鸟的IRC输入文件里看到这关键词,笔者强烈不建议用。虽然此方法不是说解决问题的可能性精确为0%,但可能性实在甚微,没有试的必要。“不收敛就直接加大循环次数上限”是菜鸟最常见的思维方式。

顺带一提,有时笔者看到有人的IRC输入文件里在IRC里用了tight,这是用来把校正步收敛限设得更严的,用这个完全是莫名其妙,也不知道从哪里学来的。本来默认的收敛限下就容易不收敛,居然还给设得更严,明显会导致出现上述报错的几率大增。
10楼2020-03-04 15:36:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Meng xf 的主题更新
信息提示
请填处理意见