sobereva 发表于 2025-8-23 12:54
本论坛里别随便提AI，试图拿AI创建计算化学的输入文件纯属搞笑，完全是让一个文科生干理科生的事。认真看下 ...

好的，谢谢老师

Uus/pMeC6H4-/キ 发表于 2025-8-23 11:31
我还真没意识到irc的stepsize参数只能用整数，似乎程序设计上L123接受IOp(1/39)时无法识别小数点、1D-1、1E ...

是stepsize填错了，谢谢讲解

Solitude198 发表于 2025-8-23 03:32
显然stepsize不能设置为小数
http://sobereva.com/400

解决问题了，感谢！

本论坛里别随便提AI，试图拿AI创建计算化学的输入文件纯属搞笑，完全是让一个文科生干理科生的事。认真看下文及其回帖里的讨论
计算新手不要用AI
http://bbs.keinsci.com/thread-52382-1-1.html

好好正经学一遍Gaussian比什么都强，别指望靠AI能帮你创建输入文件，纯属瞎走弯路
谈谈学量子化学如何正确地入门
http://sobereva.com/355（http://bbs.keinsci.com/thread-4447-1-1.html）
北京科音初级量子化学培训班（http://www.keinsci.com/KEQC）

Gaussian 16对每行的字符数完全没有70的限制

我还真没意识到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的步长已经挺精确了。

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

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

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