计算化学公社

标题: gaussian 16跑irc时为什么会出现换行报错 [打印本页]

作者
Author: 尘星 时间: 2025-8-23 03:08
标题: gaussian 16跑irc时为什么会出现换行报错
老师们好，
我在用gaussian 16跑irc时在一行中输入了以下route line：
#p irc=(calcfc,reverse,maxpoints=200,stepsize=0.05) b3lyp/6-31g(d) scf=tight

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

作者
Author: Solitude198 时间: 2025-8-23 03:32
显然stepsize不能设置为小数
http://sobereva.com/400

作者
Author: 54687 时间: 2025-8-23 10:50
另外补充一点，从G09开始默认收敛限就是scf=tight了，写这个完全就是多余的

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-8-23 11:31
我还真没意识到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的步长已经挺精确了。

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

作者
Author: sobereva 时间: 2025-8-23 12:54
本论坛里别随便提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的限制

作者
Author: 尘星 时间: 2025-8-23 14:46

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

解决问题了，感谢！

作者
Author: 尘星 时间: 2025-8-23 14:47

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

是stepsize填错了，谢谢讲解

作者
Author: 尘星 时间: 2025-8-23 14:48

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

好的，谢谢老师

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