求助过渡金属反应，为什么irc步数设置不同，分别单独进行优化后得到的结构也不同

zh7729 · 发表于 Post on 2023-6-22 19:53:46

各位老师好，我现在在B3LYP-(D3)/6-31G**(Ru使用SDD基组和SDD赝势)水平下寻找一个Ru配合物参与的异构化反应，在找到过渡态以后跑irc，步数设为50步，得到下图(1)所示的曲线，其中左端为反应物，取左端的结构在相同计算水平下单独进行优化，结果如图(2)所示，能量反而上升，但从irc图的趋势上看，能量应该是继续下降的，我一开始怀疑是因为irc跑的步数太少而导致取一端结构进行优化后跑到了错误的结构；所以我接着单独又让irc只向反应物一侧跑并将步数设为100步，得到如图(3)所示的irc曲线，同样取反应物一端的结构在相同计算水平下单独进行优化，如图(4)，能量仍然上升，并且图(4)最终得到的结构与图(2)最终得到的结构不同，同时也测试了图(2)和图(4)这两个反应物的波函数，均提示稳定，所以想请教各位老师

（1）为什么irc步数设定不一样，单独优化后得到的结果也不一样，图(2)和图(4)得到的结构哪个比较准确

（2）取一端的结构单独进行优化的过程中能量为什么会上升

（3）我的计算过程哪一步出错了，需要怎么进行调整

谢谢各位老师

wzkchem5 · 发表于 Post on 2023-6-22 20:00:32

优化的时候能量上升说明优化算法对于这个体系而言不太好。可以考虑在优化的时候算一下初始Hessian

zh7729 · 发表于 Post on 2023-6-22 20:33:03

本帖最后由 zh7729 于 2023-6-22 20:42 编辑

wzkchem5 发表于 2023-6-22 20:00
优化的时候能量上升说明优化算法对于这个体系而言不太好。可以考虑在优化的时候算一下初始Hessian

老师可以请问一下怎么计算初始Hessian吗，是在结构优化过程中使用opt(calcall)吗？

wzkchem5 · 发表于 Post on 2023-6-22 20:44:03

zh7729 发表于 2023-6-22 13:33
老师可以请问一下怎么计算初始Hessian吗，是在结构优化过程中使用opt(calcall)吗？

http://gaussian.com/opt/

zh7729 · 发表于 Post on 2023-6-22 20:46:32

wzkchem5 发表于 2023-6-22 20:44
http://gaussian.com/opt/

好的，谢谢老师！

sobereva · 发表于 Post on 2023-6-23 01:25:08

opt里加上calcfc，并且可以尝试较小的步长上限
至于是否有什么其它硬伤（比如误带了ts关键词）没输入文件没法判断

zh7729 · 发表于 Post on 2023-6-23 09:06:28

sobereva 发表于 2023-6-23 01:25
opt里加上calcfc，并且可以尝试较小的步长上限
至于是否有什么其它硬伤（比如误带了ts关键词）没输入文件 ...

好的我试一下，谢谢老师！

神仙难救得自救 · 发表于 Post on 2024-9-14 15:49:21

一般设置多少步啊

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

[Gaussian/gview] 求助过渡金属反应，为什么irc步数设置不同，分别单独进行优化后得到的结构也不同