ORCA单点任务报错 "ORCA finished by error termination in SCF"

一个用户名

在Gaussian上面做的优化，想用wB97M-V算一下单点。算几十个原子的体系没有问题，但是算百来个原子的体系就会报错：

Calculating surface distance matrix               ...        done! (  1.8s)
Performing Cholesky decomposition & store         ...        
ORCA finished by error termination in SCF
Calling Command: mpirun -np 16  /home/gaussian/orca501/orca_scf_mpi single-dsDNA-CTT-X1-CU-ORCA.gbw b single-dsDNA-CTT-X1-CU-ORCA
[file orca_tools/qcmsg.cpp, line 458]:
  .... aborting the run

Ubuntu 18.04，16核32线程，32GB；openmpi 4.1.1，ORCA 5.0.1

输入文件：

! wB97M-V def2-TZVP def2/J RIJCOSX strongSCF noautostart SMALLPRINT nopop
%maxcore 1600
%pal nprocs  16 end
%cpcm
smd true
SMDsolvent "water"
end
%basis
NewGTO Cu
"LANL2TZ(f)"
end
end
* xyz   1   2
(坐标省略)
Cu    -0.31959200   -0.62579000   -0.00793900  NewECP  "LANL2"  end
*

希望老师们不吝赐教，谢谢！

wzkchem5

检查内存、硬盘空间是否足够。注意maxcore指定的是每个核的内存（单位MB），而不是计算的总内存，maxcore乘以核数必须小于总剩余物理内存的80%

一个用户名

wzkchem5 发表于 2021-9-25 19:24
检查内存、硬盘空间是否足够。注意maxcore指定的是每个核的内存（单位MB），而不是计算的总内存，maxcore乘 ...

感谢回答。1600MB*16=25GB<32GB，应该不是内存分配过大的问题，硬盘空间也是足够的。

wzkchem5

一个用户名 发表于 2021-9-25 12:38
感谢回答。1600MB*16=25GB

抱歉我没看到你已经说了机器内存了。
这个任务如果用比较少的核跑，或者串行跑，问题仍然存在吗

一个用户名

wzkchem5 发表于 2021-9-25 20:20
抱歉我没看到你已经说了机器内存了。
这个任务如果用比较少的核跑，或者串行跑，问题仍然存在吗

感谢回答！试了一下减少到12核跑，确实可以稳定地跑起来了。而且减少核数之后，输出文件明确给出了警告：

Warning (ORCA_SCF): Not enough memory available!
                    Memory available for SCF calculation:         1600 MB
                    Memory needed (estimated)           :         1772 MB

所以我猜测之前确实是分配内存不足的问题，只是没有给出警告直接就crash了，应该算是程序设计的一个bug。

wzkchem5

一个用户名 发表于 2021-9-25 16:23
感谢回答！试了一下减少到12核跑，确实可以稳定地跑起来了。而且减少核数之后，输出文件明确给出了警告： ...

orca没有把所有内存使用都统计到，所以虽然有时可以预估到内存不足，但是不保证所有时候都会预估到，可能实际跑到那个地方才知道内存不足，但是已经晚了，被操作系统杀掉了。
另外可以检查一下同一台机子上是否跑了其他程序，可能其他程序用的内存挤占了orca的内存

七尺贱

我用wB97M-V/def2-QZVPP计算一个200原子体系的结合能任务，48核，345G，一开始给30核，5000MB都会报错，现在是15核，6000MB，会提示警告:

Warning (ORCA_SCF): Not enough memory available!
                    Memory available for SCF calculation:         6000 MB
                    Memory needed (estimated)           :         10000 MB
这意味着每个核得需要10G？这也太吃内存了吧

wzkchem5

七尺贱 发表于 2022-6-1 10:17
我用wB97M-V/def2-QZVPP计算一个200原子体系的结合能任务，48核，345G，一开始给30核，5000MB都会报错，现 ...

基组太大了，正常。可以考虑减掉f和g极化函数，或者改用ma-def2-TZVPP

七尺贱

wzkchem5 发表于 2022-6-1 17:32
基组太大了，正常。可以考虑减掉f和g极化函数，或者改用ma-def2-TZVPP

谢谢回复，我要算一个精确的弱相互作用能，所以才想用这么大的基组，而且结合了counterpoise校正

七尺贱

wzkchem5 发表于 2022-6-1 17:32
基组太大了，正常。可以考虑减掉f和g极化函数，或者改用ma-def2-TZVPP

我还有个问题就是在内存足够的情况下，这个单点任务应该不会失败了吧

wzkchem5

七尺贱 发表于 2022-6-1 10:50
我还有个问题就是在内存足够的情况下，这个单点任务应该不会失败了吧

不会因为内存问题失败，但不排除因为基组线性相关或SCF不收敛而失败

wzkchem5

七尺贱 发表于 2022-6-1 10:48
谢谢回复，我要算一个精确的弱相互作用能，所以才想用这么大的基组，而且结合了counterpoise校正

已经加了counterpoise就更不需要用quadruple zeta基组了。况且你算的是溶液相的作用能，溶剂模型的误差都明显大于ma-def2-TZVPP基组的误差

七尺贱

wzkchem5 发表于 2022-6-1 19:04
已经加了counterpoise就更不需要用quadruple zeta基组了。况且你算的是溶液相的作用能，溶剂模型的误差都 ...

好的谢谢您