计算化学公社

标题: 求助：ORCA计算卡住，警告信息“a second arrayIndex is passed”，结果还可用吗？ [打印本页]

作者
Author: dnlx 时间: 2022-6-15 20:10
标题: 求助：ORCA计算卡住，警告信息“a second arrayIndex is passed”，结果还可用吗？
本帖最后由 dnlx 于 2022-6-15 20:34 编辑

ORCA 5.0.1，在Windows 11上，体系是乙烷在石墨烯上的吸附（比社长的石墨烯少了一圈），在Iterative Optimization过程中，出现这样信息：

Warning: Variable with the name MODES is not a 2D-Array but a second arrayIndex is passed

然后程序就卡在那里了（卡住的原因不太可能是因为我输入文件存在问题）。输出文件已经消除了虚频，并且得到了自由能。请问这样的结果（优化后的结构和计算的自由能），是否可以使用？输入输出文件在附件。这种情况遇到了两次。
多谢各位指教！
(, 下载次数 Times of downloads: 9)

作者
Author: wzkchem5 时间: 2022-6-15 20:35
虚频并没有消掉。注意即使有虚频，THERMOCHEMISTRY AT 423.15K那里是不显示虚频的，需要到VIBRATIONAL FREQUENCIES那里看。更简单的方法是在整个输出文件里搜索imag。
这个可能是程序的bug，5.0.2开始修复了compound job的一些bug，可以试试下载最新版orca（5.0.3）看看

作者
Author: dnlx 时间: 2022-6-16 06:42

wzkchem5 发表于 2022-6-15 20:35
虚频并没有消掉。注意即使有虚频，THERMOCHEMISTRY AT 423.15K那里是不显示虚频的，需要到VIBRATIONAL FREQ ...

升级到5.0.3，还是同样的问题。现在用compound job来消虚频走不通，只好手调撞运气了。

作者
Author: chrinide 时间: 2022-6-16 09:36
本帖最后由 chrinide 于 2022-6-16 09:38 编辑

其实得恭喜你找到一个过渡态

Lowest eigenvalues of augmented Hessian:
-0.000000077 0.000291379 0.002111312 0.003506927 0.005847816

FINAL ENERGY EVALUATION AT THE STATIONARY POINT

ORCA已经告诉你了，你找到了一个 STATIONARY POINT
可以把“石墨烯”固定住，做个构象搜索, 再算。全构象搜索太昂贵了

作者
Author: wzkchem5 时间: 2022-6-16 13:58

dnlx 发表于 2022-6-15 23:42
升级到5.0.3，还是同样的问题。现在用compound job来消虚频走不通，只好手调撞运气了。

另外也可以试一下linux版，我有模糊的印象好像compound job在windows上有一些问题

作者
Author: dnlx 时间: 2022-6-18 05:43
网格改成defgrid3就消掉了虚频（谢谢四楼指出有虚频的是过渡态）。出现了好多次内存不足的警告，但是好歹没崩掉。

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