一个解决“Gaussian做IRC时因TS虚频过小导致走一步就停止”的方法

wjc404 · 发表于 Post on 2023-12-17 20:58:30

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

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

ionexchangeC · 发表于 Post on 2023-12-18 14:27:04

本帖最后由 ionexchangeC 于 2023-12-18 16:14 编辑

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

huanghuanhuan · 发表于 Post on 2024-11-7 13:12:28

请问这种情况下我加了IOPirc能跑出我预期的反应物和产物，但是过渡态虚频只有-20左右实在太小，这能算正确的过渡态吗

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

[Gaussian/gview] 一个解决“Gaussian做IRC时因TS虚频过小导致走一步就停止”的方法

评分 Rate