计算化学公社

标题: IRC正常结束但两端能量明显没有走平 [打印本页]

作者
Author: naoki 时间: 2021-10-17 22:54
标题: IRC正常结束但两端能量明显没有走平
本帖最后由 naoki 于 2021-10-17 22:57 编辑

如题，走IRC时两端明显没有走平，也没有走到maxpoints设定的步数，就正常结束了，反应物和产物结构看起来都还没走到最终的样子。想请教一下各位老师，是什么原因导致出现了这种情况，是否可以直接取两端结构进行几何优化了呢，谢谢大家~
输入指令：irc=(calcfc,maxpoints=300,recalc=3,lqa) b3lyp/6-31g(d,p) em=gd3bj scrf=(smd,solvent=n-hexane)
(, 下载次数 Times of downloads: 16) (, 下载次数 Times of downloads: 24)
放大看两端是这样
(, 下载次数 Times of downloads: 24) (, 下载次数 Times of downloads: 19)

作者
Author: GoldenBaby 时间: 2021-10-17 23:08
直接取两边优化就可以了，你这个IRC走的点已经足够多了。跑IRC又不是几何优化，走不到极小点很正常吧。

作者
Author: naoki 时间: 2021-10-17 23:27

GoldenBaby 发表于 2021-10-17 23:08
直接取两边优化就可以了，你这个IRC走的点已经足够多了。跑IRC又不是几何优化，走不到极小点很正常吧。

谢谢，以前做IRC两端都走得相对平一些，这次看着总感觉很别扭

作者
Author: sobereva 时间: 2021-10-18 01:59
Gaussian判断IRC什么时候自动结束（预测继续跑下去能量会上升）的方法不是很理想，而且里面也不给用户留出可调参数，没法总是期望能跑到接近平坦（如果允许Gaussian根据IRC方向受力大小判断什么时候结束，而且阈值可调的话，会好得多）
使用不同的IRC参数、不同的算法，IRC最多能跑多长都是不同的。对于优化两侧结构的目的，IRC稍微跑跑就行，当前情况完全够了

作者
Author: naoki 时间: 2021-10-18 09:37

sobereva 发表于 2021-10-18 01:59
Gaussian判断IRC什么时候自动结束（预测继续跑下去能量会上升）的方法不是很理想，而且里面也不给用户留出 ...

原来如此，谢谢sob老师解惑

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