使用ORCA做DLPNO-NEVPT2计算报错不知如何解决

阿甘

各位老师好，学生最近在使用ORCA 5.0.3做DLPNO-NEVPT2计算有机物（体系共90原子，全是CHNO原子）的singlet-triplet gap的时候报错了，搜索了本论坛和ORCA论坛但未找到解决办法，不知如何解决。烦请各位老师帮我看看，不胜感激。以下是输出文件报错附近的几行：
------------------------------------------
DLPNO-NEVPT2 CALCULATION FOR Viab  CLASS
------------------------------------------

DOI based multipole  prescreening for Viab     ... used
Estimated Energy For SKIPPED Pairs             ...       0.000000000000 Eh
DLPNO Viab  Class PreScr                       ... done (3443.7 sec)
The Iteration Cycle OF #    0 Eigen            ... start
LMP2-Iter=   0:  ELMP2= -0.192111762050  DE=  -9.2e-03 RMAX=   1.6e-03
LMP2-Iter=   1:  ELMP2= -0.196089743518  DE=  -4.0e-03 RMAX=   6.7e-04
LMP2-Iter=   2:  ELMP2= -0.196979052061  DE=  -8.9e-04 RMAX=   2.3e-04
LMP2-Iter=   3:  ELMP2= -0.197063763148  DE=  -8.5e-05 RMAX=   6.0e-05
LMP2-Iter=   4:  ELMP2= -0.197067132572  DE=  -3.4e-06 RMAX=   1.4e-05
LMP2-Iter=   5:  ELMP2= -0.197067420017  DE=  -2.9e-07 RMAX=   3.4e-06
LMP2-Iter=   6:  ELMP2= -0.197067482790  DE=  -6.3e-08 RMAX=   8.2e-07
LMP2-Iter=   7:  ELMP2= -0.197067485155  DE=  -2.4e-09 RMAX=   1.8e-07 => CONVERGED
The Iteration Cycle OF #    1 Eigen            ... start
LMP2-Iter=   0:  ELMP2= -0.038247508322  DE=  -1.7e-03 RMAX=   7.9e-04
LMP2-Iter=   1:  ELMP2= -0.038970521498  DE=  -7.2e-04 RMAX=   3.2e-04
LMP2-Iter=   2:  ELMP2= -0.039121455697  DE=  -1.5e-04 RMAX=   1.0e-04
LMP2-Iter=   3:  ELMP2= -0.039133719465  DE=  -1.2e-05 RMAX=   2.2e-05
LMP2-Iter=   4:  ELMP2= -0.039134081312  DE=  -3.6e-07 RMAX=   4.8e-06
LMP2-Iter=   5:  ELMP2= -0.039134133870  DE=  -5.3e-08 RMAX=   1.1e-06
LMP2-Iter=   6:  ELMP2= -0.039134147873  DE=  -1.4e-08 RMAX=   2.5e-07 => CONVERGED
The Iteration Cycle OF #    2 Eigen            ... start
LMP2-Iter=   0:  ELMP2=  0.000000000000  DE=   1.4e-07 RMAX=   0.0e+00
Error in LINEQ_Gauss_Jordan:
The solution could not be computed.

ORCA finished by error termination in CASSCF
Calling Command: mpirun -np 24  /public1/home/sch4760/apps/orca_5_0_3_linux_x86-64_shared_openmpi411/orca_casscf_mpi DLPNO-NEVPT2.casinp.tmp DLPNO-NEVPT2
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run



任务已经在CASSCF上收敛了，但是不知这个报错信息应该如何处理。以下是输入文件的命令：
! def2-TZVPD def2-TZVP/C RIJCOSX def2/J
! moread
%moinp "CASSCF-RI-BS11.gbw"
%scf
thresh 1e-12
tcut 1e-14
sthresh 1e-5
end
%casscf nel    2
        norb   4
        nroots 1
        mult   1
        OrbStep SuperCI_PT
        SwitchStep DIIS
        SwitchConv 0.003
        SwitchIter 30
        shiftup 2
        shiftdn 2
        TrafoStep RI
        PTMethod DLPNO_NEVPT2
         PTSettings
          TCutPNO 1e-8
          MaxIter 20
          MaxDIIS 7
         end
        end
%maxcore  10000
%pal nprocs  24 end
* xyz   0   1

shalene

我也不知道问题怎么解决。
不过你可以换个小体系算算，验证下所有的关键词

另外，def2-TZVPD def2-TZVP/C ?   应该是def2-tzvpd/c，没d的话，数值误差是不是太大了？

阿甘

shalene 发表于 2023-1-6 11:44
我也不知道问题怎么解决。
不过你可以换个小体系算算，验证下所有的关键词

感谢老师的建议和提醒，我原来测试时用的ma加弥散，换了基组之后忘记改辅助基组了，马上把它换了。

yuyangzap

八成错误都是内存不足，减小并行核心，加大内存。

阿甘

yuyangzap 发表于 2023-1-6 14:18
八成错误都是内存不足，减小并行核心，加大内存。

感谢老师的建议，我试试看能不能解决这个报错。