计算化学公社

标题: 一个解决“Gaussian做IRC时因TS虚频过小导致走一步就停止”的方法 [打印本页]

作者
Author: wjc404 时间: 2023-12-17 20:58
标题: 一个解决“Gaussian做IRC时因TS虚频过小导致走一步就停止”的方法
在用Gaussian跑IRC时，有时会遇到在两个方向上各走一个点就停止的情况，查看输出文件发现是程序判断该点的梯度小于某个门限值，就认为到达终点了，不再往前走了。Gaussian程序的帮助里面没有说怎么把这个梯度门限值调小，不过 gaussian.com/overlay1 中介绍 IOp(1/108) 为 "Gradient magnitude for DVV stopping criteria (Crit1)."。实测这个IOp对于以HPC为IRC算法的任务亦有效，能够关闭“梯度低于某个值后停止”的逻辑。

PS: 还有一种做法是调大步长，不过这个有可能会降低IRC的精度。

作者
Author: ionexchangeC 时间: 2023-12-18 14:27
本帖最后由 ionexchangeC 于 2023-12-18 16:14 编辑

实测有效，注意用IOp(1/108=-1)时，IRC不会自行中断，会跑满设置的步数
更正：还是可以自动停止的，我触发的是Step angle

作者
Author: huanghuanhuan 时间: 2024-11-7 13:12
请问这种情况下我加了IOPirc能跑出我预期的反应物和产物，但是过渡态虚频只有-20左右实在太小，这能算正确的过渡态吗

作者
Author: magicrao 时间: 2025-3-14 13:44
感谢分享，很好的技巧！

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