计算化学公社»论坛首页 › 理论与计算化学 (Theoretical and Computational Chemistry) › 量子化学 (Quantum Chemistry) › IRC用HPC第一步就校正失败

IRC用HPC第一步就校正失败

查看数: 3172 | 评论数: 14 | 收藏 Add to favorites 0

关灯 | 提示：支持键盘翻页<-左右->

帖子模式

mfdsrax2

发布时间: 2022-5-5 22:28

正文摘要:

本帖最后由 mfdsrax2 于 2022-5-5 22:38 编辑用M062X/6-311G(2D,P)优化了过渡态，过渡态有53i的虚频，对应氯离子的离去，与反应方向一致。然后在相同级别上计算IRC，用默认的HPC方法，第一步就校正失败，仔细 ...

回复 Reply

sobereva 发表于 Post on 2022-5-7 15:48:32

mfdsrax2 发表于 2022-5-7 10:46
校正步不收敛，默认的20次

甭用HPC就完了
又不想用LQA的话就用GS2

mfdsrax2 发表于 Post on 2022-5-7 10:46:13

sobereva 发表于 2022-5-7 10:09
具体是什么不收敛？

校正步不收敛，默认的20次

sobereva 发表于 Post on 2022-5-7 10:09:03

mfdsrax2 发表于 2022-5-7 09:38
卢老师，我觉得我尽力了，先用OPT=TIGHT INT=SUPERFINE重新优化了过渡态，然后在相同级别上，用IRC=(RCFC ...

具体是什么不收敛？

mfdsrax2 发表于 Post on 2022-5-7 09:38:05

sobereva 发表于 2022-5-6 18:59
想更平滑就calcall
用OPT=TIGHT INT=SUPERFINE没什么意义

卢老师，我觉得我尽力了，先用OPT=TIGHT INT=SUPERFINE重新优化了过渡态，然后在相同级别上，用IRC=(RCFC,STEPSIZE=1,RECALC=1) INT=SUPERFINE再次计算IRC，仍然第一个点就不收敛，还有比这更强的条件吗？STEPSIZE设置成负数的话，步长会更小吗？

sobereva 发表于 Post on 2022-5-6 18:59:24

mfdsrax2 发表于 2022-5-6 14:28
减小步长我试过STEPSIZE=15，也不行，正逆都是1个点就停了，只有20才行。而且20也必须是LQA方法，用GS2也 ...

想更平滑就calcall
用OPT=TIGHT INT=SUPERFINE没什么意义

RAL 发表于 Post on 2022-5-6 16:15:41

mfdsrax2 发表于 2022-5-6 15:44
这是一个项目，前期交的计算材料里面已经包含IRC的能量和梯度图，后面的需要和前面保持一致

前面的图 ...

这坑挖得挺深

可以考虑用AIMD把这一步反应跑通来代替TS和IRC

RAL 发表于 Post on 2022-5-6 14:49:06

mfdsrax2 发表于 2022-5-6 14:28
减小步长我试过STEPSIZE=15，也不行，正逆都是1个点就停了，只有20才行。而且20也必须是LQA方法，用GS2也 ...

是什么目的要求梯度也要像能量曲线那样平滑呢

sobereva 发表于 Post on 2022-5-6 13:37:01

想更光滑就RECALC=3

用HPC时，提供更好的Hessian和更小的步长都能减小出现校正步不收敛的概率

Freeman 发表于 Post on 2022-5-6 12:37:21

本帖最后由 Freeman 于 2022-5-6 12:41 编辑

比较平的过渡态就是这样的，正常现象。如果你只要过渡态旁边的极小点，可以从过渡态沿着振动方向拉一下分子，然后用calcall或recalcfc优化；如果你要的是整个路径，可以用orca的neb功能。

mfdsrax2 发表于 Post on 2022-5-6 10:07:01

一颗赛艇发表于 2022-5-6 08:46
我怀疑你的过渡态找得就不对

波数是不大，但是振动方向和预期一样，LQA跑出来的IRC也和反应方向一致

一颗赛艇 发表于 Post on 2022-5-6 08:46:48

我怀疑你的过渡态找得就不对

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register