计算化学公社

标题: 求助multistep gaussian计算，第一步QST3可以算出来，第二部IRC算不出来 [打印本页]

作者
Author: fresco604 时间: 2016-9-18 23:30
标题: 求助multistep gaussian计算，第一步QST3可以算出来，第二部IRC算不出来
新手求助，我要寻找两个生物分子之间的过渡态和pathway，我使用multi-step Gaussian计算，第一步QST3计算没有任何问题，可以生成过渡态，但第二部从过渡态出发进行IRC计算始终不能成功，比如得到这样的错误信息：

CORRECTOR INTEGRATION CONVERGENCE:
Recorrection delta-x convergence threshold: 0.010000
Delta-x Convergence NOT Met
Maximum number of corrector steps exceded.

这是我的input文件，省去了分子结构：

%chk= 文件名
%mem=8000MB
# hf/6-31+g(d) opt=qst3 freq scf=(maxconventionalcycles=700,xqc)

分子1

0 1
分子结构1

分子2

0 1
分子结构2

过渡态

0 1
过渡态结构

--Link1--
%chk=
%mem=8000MB
# irc=(maxpoints=50,rcfc) hf/6-31+g(d) geom=checkpoint Guess=Read

过渡态

0 1

请诸位大神帮忙看看，是不是需要在IRC计算里加上别的什么关键词？我总是不理解为什么QST3可以计算出来IRC算不出来，他们不是相同的pathway吗？

小弟在此谢过！

作者
Author: liyuanhe211 时间: 2016-9-19 01:59
本帖最后由 liyuanhe211 于 2016-9-20 21:59 编辑

irc=lqa
两步独立，不一回事遇到问题先简化、排除一些无关变量是科研人员应有的基本素质和技能。例如把第一步的坐标拿到第二步算就知道这个问题和多步计算没关系。

作者
Author: fresco604 时间: 2016-9-19 03:48

liyuanhe211 发表于 2016-9-19 01:59
irc=lqa
两步独立，不一回事遇到问题先简化、排除一些无关产量是科研人员应有的基本素质和技能。例如把第 ...

多谢帮助，如果我在第一步中计算力常数，加上opt=calcfc, 然后第二部当中用 irc=rcfc 是不是就是读取第一步中算出来的力常数？

作者
Author: liyuanhe211 时间: 2016-9-19 07:28

fresco604 发表于 2016-9-19 03:48
多谢帮助，如果我在第一步中计算力常数，加上opt=calcfc, 然后第二部当中用 irc=rcfc 是不是就是读取第一 ...

读freq

作者
Author: fresco604 时间: 2016-9-19 09:03
本帖最后由 fresco604 于 2016-9-19 09:08 编辑

liyuanhe211 发表于 2016-9-19 07:28
读freq

哦对对，我忘了，多谢指教

我再重新试一下irc=lqa, 刚才把checkpoint file写错了

作者
Author: fresco604 时间: 2016-9-20 21:59

liyuanhe211 发表于 2016-9-19 01:59
irc=lqa
两步独立，不一回事遇到问题先简化、排除一些无关产量是科研人员应有的基本素质和技能。例如把第 ...

加上irc=lqa还真算出来了，尽管算了600多步....谢谢老师！

作者
Author: liyuanhe211 时间: 2016-9-20 22:00

fresco604 发表于 2016-9-20 21:59
加上irc=lqa还真算出来了，尽管算了600多步....谢谢老师！

IRC跑多少个点明明是你自己指定的

作者
Author: fresco604 时间: 2016-9-20 22:28

liyuanhe211 发表于 2016-9-20 22:00
IRC跑多少个点明明是你自己指定的

是的，但是跑少了就跑不到我的local minimum, 所以我就随机给了一个很大的值，比如maxpoints=400,maxcycle=800，最后跑了593个点

作者
Author: liyuanhe211 时间: 2016-9-20 22:51

fresco604 发表于 2016-9-20 22:28
是的，但是跑少了就跑不到我的local minimum, 所以我就随机给了一个很大的值，比如maxpoints=400,maxcycl ...

跑一点IRC做优化就可以了。非要靠IRC找到可以把步长搞大点儿

作者
Author: fresco604 时间: 2016-9-21 00:33

liyuanhe211 发表于 2016-9-20 22:51
跑一点IRC做优化就可以了。非要靠IRC找到可以把步长搞大点儿

好的，那我改一改stepsize

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