为什么会出现换行报错

尘星

老师们好，
我在用gaussian 16跑irc时在一行中输入了以下route line：
#p irc=(calcfc,reverse,maxpoints=200,stepsize=0.05) b3lyp/6-31g(d) scf=tight

但始终出现以下报错：
在尝试了手打输入文件、分行输入、图形化界面输入、重启后操作、ai输出后复制（与原文件相同）、ai直接输出.gjf文件运行等方法后仍然显示完全相同的报错。
此前多次出现这个现象，我注意到报错的断行都在第70个和71个字符之间，但没能找到有效的解决方法。
感谢各位老师

Solitude198

显然stepsize不能设置为小数
http://sobereva.com/400

54687

另外补充一点，从G09开始默认收敛限就是scf=tight了，写这个完全就是多余的

Uus/pMeC6H4-/キ

我还真没意识到irc的stepsize参数只能用整数，似乎程序设计上L123接受IOp(1/39)时无法识别小数点、1D-1、1E-1这样的浮点数格式。

有一种可能是你对stepsize的单位和0.05这个参数的单位理解有误：(1) stepsize为正时单位为0.01 Bohr（程序自动转换为质权坐标的amu^(1/2) Bohr处理），如果0.05的单位是Bohr那应该写成stepsize=5；(2) stepsize为负时单位为0.01 amu^(1/2) Bohr，如果0.05的单位就是这个的话那应该写成stepsize=-5（比stepsize=5的步长小）。

如果确信单位没问题，那得清楚意识到0.05 * 0.01 Bohr的步长是超级超级短的。就算没有前述的软件限制，除非势能面描述极其精确、在过渡态附近的确十分陡峭，否则很有可能走一步就碰到校正步无法收敛（用默认的HPC方法）或者因电子能变化极小而被判定为走到极小点并自动结束（用LQA方法时）。就算连这些问题都没有，那短步长意味着同样反应路径需要更多数量的点来描述，甚至可能需要成千上万个点，显然计算资源开销也是难以承受的。对于得到完整高质量的连接过渡态与极小点的IRC曲线而言，stepsize=5的步长已经挺精确了。

建议把输入文件传一下，可以顺便检查检查结构。