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

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看能量已经维持在一个固定值优化多次没有变化，遇到这样的问题我应该怎么处理？谢谢老师！

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

好的，谢谢老师

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

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方向得到了平滑的曲线。我接下来应该怎么办？是不是说明我的过渡态存在问题？谢谢老师

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

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

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

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.（我看您的之前的帖子说这不是报错）正常结束。但是我在查看路径时发现只有两个点，这样可以吗？我应该怎么办？谢谢您！

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

好的，谢谢您！

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

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

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

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

用stepsize减小步长（步长越小越平滑），或者用recalc=N设定每N步重新精确计算一次Hessian（N越小越平滑但也越耗时）