计算化学公社

标题: ORCA NEVPT2不明原因报错 [打印本页]

作者
Author: ZetaFunction 时间: 2025-6-16 05:06
标题: ORCA NEVPT2不明原因报错
最近在尝试使用ORCA6.0.0做NEVPT2计算，对77个原子的双镍离子配位体系取CAS(14o, 24e)的活性空间，CASSCF成功收敛到了期望的轨道空间，然而后续的NEVPT2计算却出了问题。在使用def2-SVP基组时，默认的FIC_NEVPT2方法出现以下报错：

Class TUAB ... Skipped = 0 of 71267364 LowestEV= 1.266251e+00 ... done in 0.9 sec
Class IJTA ... [file orca_casscf/cas_ficpt.cpp, line 1884, Process 2]: ERROR (ORCA_NEVPT2): V1_ija MatrixContainer IKJA
[file orca_casscf/cas_ficpt.cpp, line 1884, Process 10]: ERROR (ORCA_NEVPT2): V1_ija MatrixContainer IKJA
[file orca_casscf/cas_ficpt.cpp, line 1884, Process 11]: ERROR (ORCA_NEVPT2): V1_ija MatrixContainer IKJA
[file orca_casscf/cas_ficpt.cpp, line 1884, Process 12]: ERROR (ORCA_NEVPT2): V1_ija MatrixContainer IKJA
...
ORCA finished by error termination in CASSCF
Calling Command: mpirun -np 128 /dssg/home/acct-clsst/clsst/orca_6_0_0/orca_casscf_mpi --energy IM1-CAS14-SVP-FIC_NEVPT2.casinp.tmp IM1-CAS14-SVP-FIC_NEVPT2
[file orca_tools/qcmsg.cpp, line 394]:
.... aborting the run

复制代码

换用SC_NEVPT2方法后，计算成功，耗时不到半小时。然而再尝试将基组加大到def2-TZVP后，无论是使用FIC_NEVPT2还是SC_NEVPT2都出现了同样的报错：

--------------------
MULT 1, ROOT 0
--------------------
--- Canonicalize Internal Space
--- Canonicalize External Space
3 total processes killed (some possibly by mpirun during cleanup)
ORCA finished by error termination in CASSCF
Calling Command: mpirun -np 64 /dssg/home/acct-clsst/clsst/orca_6_0_0/orca_casscf_mpi --energy IM1-CAS14-TZVP-SC_NEVPT2.casinp.tmp IM1-CAS14-TZVP-SC_NEVPT2
[file orca_tools/qcmsg.cpp, line 394]:
.... aborting the run

复制代码

一开始考虑可能是内存问题，学校的服务器每个节点有64个核心，每个核心8G内存，总共512G，输入文件中每个进程给了6G内存。我尝试过将内存改为每个进程4G，同时使用两个节点共128G，理应不会出现内存不足的情况，但是依然发送同样的报错。请问各位，这是什么原因导致的？是程序存在BUG吗？

后来又尝试了DLPNO_NEVPT2方法并同时开启RI加速，def2-TZVP基组下计算成功了，但是用时15个小时，对比def2-SVP下的SC_NEVPT2用时似乎增加的太多了，不知道这个方法是否稳健。

我也尝试使用PySCF做NEVPT2计算，但是PySCF似乎会在计算途中产生特别大的临时文件，大约2T左右，超过了服务器上的剩余可用存储空间。

作者
Author: wzkchem5 时间: 2025-6-16 08:27
应该是硬盘不足。MatrixContainer是把矩阵存到硬盘上的一种数据结构，MatrixContainer报错估计是硬盘空间不够

作者
Author: ZetaFunction 时间: 2025-6-16 15:11

wzkchem5 发表于 2025-6-16 08:27
应该是硬盘不足。MatrixContainer是把矩阵存到硬盘上的一种数据结构，MatrixContainer报错估计是硬盘空间不 ...

感谢王老师回复！

服务器上的剩余空间有2T左右，对于77个原子的NEVPT2(14o,24e)来说也不够吗？只有在def2-SVP基组下的FIC_NEVPT2计算才会遇到MatrixContainer报错，而def2-SVP基组下的SC_NEVPT2计算和def2-TZVP基组下的DLPNO_NEVPT2计算却成功了，是因为FIC_NEVPT2特别吃存储空间吗？

现在只有DLPNO_NEVPT2可以成功完成计算，请问这个方法的精度对比普通的NEVPT2来说怎么样？

作者
Author: wzkchem5 时间: 2025-6-16 15:13

ZetaFunction 发表于 2025-6-16 15:11
感谢王老师回复！

服务器上的剩余空间有2T左右，对于77个原子的NEVPT2(14o,24e)来说也不够吗？只有在d ...

FIC-NEVPT2是这3个方法里最精确的，也是计算资源消耗最多的。开RI有助于加速FIC-NEVPT2，不过当体系比较大时，开RI不如开DLPNO（DLPNO同时也会开RI）。
关于DLPNO-NEVPT2的精度，参见https://doi.org/10.1063/1.4942769。可以认为非常接近FIC-NEVPT2

作者
Author: ZetaFunction 时间: 2025-6-16 15:45

wzkchem5 发表于 2025-6-16 15:13
FIC-NEVPT2是这3个方法里最精确的，也是计算资源消耗最多的。开RI有助于加速FIC-NEVPT2，不过当体系比较 ...

所以说即使是def2-SVP基组下77个原子的FIC_NEVPT2(14,24)也需要好几T的空间用来储存临时文件吗？如果是这样那我还是用DLPNO了。

作者
Author: wzkchem5 时间: 2025-6-16 15:48

ZetaFunction 发表于 2025-6-16 15:45
所以说即使是def2-SVP基组下77个原子的FIC_NEVPT2(14,24)也需要好几T的空间用来储存临时文件吗？如果是这 ...

我不确定，但是这个很好验证，去看一下临时文件夹有没有被撑爆就知道了

作者
Author: ZetaFunction 时间: 2025-6-16 15:55

wzkchem5 发表于 2025-6-16 15:48
我不确定，但是这个很好验证，去看一下临时文件夹有没有被撑爆就知道了

临时文件是那些.tmp.文件吗？我在计算途中检查过，文件夹内的文件总量没有超过100G，不过也可能是最大的文件因为空间不足还没产生任务就结束了。我尝试过使用PySCF的NEVPT2，在报错的输出中确实发现了程序试图生成一个超过2T的临时文件。

作者
Author: ZetaFunction 时间: 2025-6-18 19:15

wzkchem5 发表于 2025-6-16 08:27
应该是硬盘不足。MatrixContainer是把矩阵存到硬盘上的一种数据结构，MatrixContainer报错估计是硬盘空间不 ...

我又尝试了一个很小的体系，只有12个原子的CAS(2,2)，在def2-SVP基组下做FIC_NEVPT2依然有同样的报错，这次总不可能是存储空间不足的问题吧。

作者
Author: wzkchem5 时间: 2025-6-18 20:15

ZetaFunction 发表于 2025-6-18 19:15
我又尝试了一个很小的体系，只有12个原子的CAS(2,2)，在def2-SVP基组下做FIC_NEVPT2依然有同样的报错，这 ...

如果这个问题用ORCA 6.1.0也会出现，可以在ORCA论坛上报一下bug

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