计算化学公社

标题: 关于IRC同时用lqa后的结果是否正确求助各位老师 [打印本页]

作者
Author: chrishuo 时间: 2017-8-24 15:45
标题: 关于IRC同时用lqa后的结果是否正确求助各位老师
我找到了一个过渡态，有一个虚频，振动强度为-142，振动模式是我想要的结果，但是在做IRC分析时总是出现一个突变点，结果报错无法收敛（Delta-x Convergence NOT Met），除了此突变点reverse和forward趋势都是正确的。（forward、reverse两个方向分别计算）
后来我在论坛上看到可以加上lqa的命令，重新计算正常结束，找到了反应物和产物，但是路径分析的图中出现很多拐点，不平滑。
(, 下载次数 Times of downloads: 53) 这样可以认为正确吗？
(, 下载次数 Times of downloads: 22) 我又用同样的方法计算了另外一个我得到的过渡态，有的拐点甚至高于初始点，是不是我找的这个过渡态存在问题？谢谢大家！

作者
Author: sobereva 时间: 2017-8-24 15:49
用stepsize减小步长（步长越小越平滑），或者用recalc=N设定每N步重新精确计算一次Hessian（N越小越平滑但也越耗时）

作者
Author: chrishuo 时间: 2017-8-24 15:53

sobereva 发表于 2017-8-24 15:49
用stepsize减小步长（步长越小越平滑），或者用recalc=N设定每N步重新精确计算一次Hessian（N越小越平滑但 ...

谢谢sob老师，我尝试一下减小步长。老师，我得到的第二个图那样的结果算不算是过渡态存在问题呀？我需要重新找吗？

作者
Author: liyuanhe211 时间: 2017-8-24 16:19

chrishuo 发表于 2017-8-24 15:53
谢谢sob老师，我尝试一下减小步长。老师，我得到的第二个图那样的结果算不算是过渡态存在问题呀？我需要 ...

这曲线看着像数值问题，先做平滑了再说。

作者
Author: chrishuo 时间: 2017-8-24 16:24

liyuanhe211 发表于 2017-8-24 16:19
这曲线看着像数值问题，先做平滑了再说。

好的，谢谢您！

作者
Author: chrishuo 时间: 2017-8-24 19:50

liyuanhe211 发表于 2017-8-24 16:19
这曲线看着像数值问题，先做平滑了再说。

我针对第二个图减小了步长（之前的步长是30），#p irc=(forward,calcfc,LQA,stepsize=5) b3lyp/6-31g(d)任务正常结束，输出文件中出现了This type of calculation cannot be archived.（我看您的之前的帖子说这不是报错）正常结束。但是我在查看路径时发现只有两个点，这样可以吗？我应该怎么办？谢谢您！

作者
Author: sobereva 时间: 2017-8-24 20:22

chrishuo 发表于 2017-8-24 19:50
我针对第二个图减小了步长（之前的步长是30），#p irc=(forward,calcfc,LQA,stepsize=5) b3lyp/6-31g(d) ...

(, 下载次数 Times of downloads: 38)

作者
Author: chrishuo 时间: 2017-8-25 20:35

sobereva 发表于 2017-8-24 20:22

好的，我再试试看，谢谢老师！

作者
Author: chrishuo 时间: 2017-8-26 20:04

sobereva 发表于 2017-8-24 20:22

sob老师，我用这样的命令#p opt=(ts,calcall,noeigentest,tight) int=ultrafine freq b3lyp/6-31g(d)重新计算了过渡态，对得到的过渡态进行IRC分析时，forward（给的命令是#p irc=(forward,calcfc,recalc=3,LQA,stepsize=10) int=ultrafine b3lyp/6-31g(d)）方向仍旧只有两个点，reverse方向得到了平滑的曲线。我接下来应该怎么办？是不是说明我的过渡态存在问题？谢谢老师

作者
Author: sobereva 时间: 2017-8-26 20:14

chrishuo 发表于 2017-8-26 20:04
sob老师，我用这样的命令#p opt=(ts,calcall,noeigentest,tight) int=ultrafine freq b3lyp/6-31g(d)重新 ...

stepsize设到15或20，如果能走下去但是还不平滑，把LQA去掉试试（由于用了recalc=3了，重校正步不收敛的错误几率已经降低不少）看是否已经足够平滑，如果不行，再把recalc设2乃至1（等价于calcall）。
或者尝试用GS2算法产生IRC

作者
Author: chrishuo 时间: 2017-8-26 20:38

sobereva 发表于 2017-8-26 20:14
stepsize设到15或20，如果能走下去但是还不平滑，把LQA去掉试试（由于用了recalc=3了，重校正步不收敛 ...

好的，谢谢老师

作者
Author: chrishuo 时间: 2017-8-29 09:20

sobereva 发表于 2017-8-26 20:14
stepsize设到15或20，如果能走下去但是还不平滑，把LQA去掉试试（由于用了recalc=3了，重校正步不收敛 ...

sob老师，我在用#p opt=(ts,calcall,noeigentest,tight) int=ultrafine freq b3lyp/6-31g(d)对已经找到的过渡态重新计算时，计算了将近4天还没有结束，我用chemcraft看能量已经维持在一个固定值优化多次没有变化，遇到这样的问题我应该怎么处理？谢谢老师！

作者
Author: sobereva 时间: 2017-8-29 09:48

chrishuo 发表于 2017-8-29 09:20
sob老师，我在用#p opt=(ts,calcall,noeigentest,tight) int=ultrafine freq b3lyp/6-31g(d)对已经找到的 ...

用gview看，甭用chemcraft，直接就能给你能量和受力变化，看有无收敛的迹象。
calcall显然慢。用了calcall就没有太大必要用tight，这导致难收敛，而且也没必要用ultrafine

作者
Author: chrishuo 时间: 2017-8-29 09:58
本帖最后由 chrishuo 于 2017-8-29 10:10 编辑

sobereva 发表于 2017-8-29 09:48
用gview看，甭用chemcraft，直接就能给你能量和受力变化，看有无收敛的迹象。
calcall显然慢。用了calca ...

老师，我用gaussview看了一下，发现能量也不再变化了 (, 下载次数 Times of downloads: 28)

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