计算化学公社

标题: 关于OfakeG（1.3.3）转换过渡态计算文件的一个小建议 [打印本页]

作者
Author: ABetaCarw 时间: yesterday 21:22
标题: 关于OfakeG（1.3.3）转换过渡态计算文件的一个小建议
我发现当用ORCA计算过渡态时候，比如用类似如下关键词
! B3LYP D4 def2-SVP def2/J optTS freq tightSCF noautostart miniprint
%maxcore    4000
%pal nprocs  8 end
%geom Calc_Hess true end

然后计算时候会在一开始做一次hess
当这个计算最终未收敛的时候，比如报如下错误
                        ERROR !!!
The optimization did not converge but reached the maximum
number of optimization cycles.
As a subsequent Frequencies calculation has been requested
ORCA will abort at this point of the run.
Please restart the calculation with the lowest energy geometry and/or
a larger maxiter for the geometry optimization.

此时用OfakeG转换输出文件
会发现在打开伪高斯输出文件的时候，在最后一帧的结构上可以正常看到频率分析结果。但是这个频率结果其实是过渡态计算一开始做的那次hess

这个很可能会让没有检查原始ORCA输出文件的人错误判断计算状态。

我的建议是，对于没有收敛的ORCA过渡态计算输出文件在转换时候给予报错提示，或者直接屏蔽掉这种状态的频率信息的转换，当用户发现没有转换出频率结果时候也会自发去查看原始数据（不排除直接伸手发帖问的）

作者
Author: Hugo_314cat 时间: 21 min ago
这个问题我也发现了，感谢您在此提出。虽然个人觉得点开看一下输出文件状态是基本素养，但是如果能够提供一个报错提示（例如通过fake方式在生成文件后用GaussView打开也能弹出一个GView的报错框），我觉得能有不少人因此而获益。

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