计算化学公社

标题: 求助过渡金属反应，为什么irc步数设置不同，分别单独进行优化后得到的结构也不同 [打印本页]

作者
Author: zh7729 时间: 2023-6-22 19:53
标题: 求助过渡金属反应，为什么irc步数设置不同，分别单独进行优化后得到的结构也不同

各位老师好，我现在在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）我的计算过程哪一步出错了，需要怎么进行调整

谢谢各位老师

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

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

作者
Author: zh7729 时间: 2023-6-22 20:33
本帖最后由 zh7729 于 2023-6-22 20:42 编辑

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

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

作者
Author: wzkchem5 时间: 2023-6-22 20:44

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

http://gaussian.com/opt/

作者
Author: zh7729 时间: 2023-6-22 20:46

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

好的，谢谢老师！

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

作者
Author: zh7729 时间: 2023-6-23 09:06

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

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

作者
Author: 神仙难救得自救 时间: 2024-9-14 15:49
一般设置多少步啊

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