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

wjc404

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

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

ionexchangeC

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

huanghuanhuan

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

magicrao

感谢分享，很好的技巧！