orca算单点能，SCF收敛后，MP2部分计算停止

魔卡少女晗

输入文件是这样的：
! PWPB95 D3 ma-def2-QZVPP autoaux RIJCOSX noautostart miniprint nopop tightscf
%maxcore  30000
%pal nprocs   4 end
* xyz   0   1
C                 -2.10736400    1.82412400    0.79815400
C                 -0.84999400    2.41788200    0.89723900
C                  0.16170300    2.05981400    0.00149400
C                 -0.07708300    1.09587700   -0.97652200
C                 -1.34619900    0.51973400   -1.06638600
C                 -2.37757600    0.86203000   -0.18609500
H                 -2.87921500    2.11722400    1.51042300
输出文件结尾是这样的：
1961   0.0000     198.557194      5403.0159
1962   0.0000     254.462691      6924.2818
1963   0.0000     254.773402      6932.7367
Total SCF time: 0 days 4 hours 25 min 56 sec


           ************************************************************
           *        Program running with 4 parallel MPI-processes     *
           *              working on a common directory               *
           ************************************************************

------------------------------------------------------------------------------
                                ORCA  MP2
------------------------------------------------------------------------------

Freezing NCore=44 chemical core electrons
[file orca_mp2/mp2rigrad.cpp, line 790, Process 0]: ERROR ORCA_MP2: RIFLAG[0]!=1 (Cholesky) is not allowed here.

[file orca_mp2/mp2rigrad.cpp, line 790, Process 2]: ERROR ORCA_MP2: RIFLAG[0]!=1 (Cholesky) is not allowed here.

[file orca_mp2/mp2rigrad.cpp, line 790, Process 3]: ERROR ORCA_MP2: RIFLAG[0]!=1 (Cholesky) is not allowed here.

[file orca_mp2/mp2rigrad.cpp, line 790, Process 1]: ERROR ORCA_MP2: RIFLAG[0]!=1 (Cholesky) is not allowed here.


ORCA finished by error termination in MP2
Calling Command: mpirun -np 4  /root/orca/orca_mp2_mpi ILwithoutF-R.mp2inp.tmp ILwithoutF-R
[file orca_tools/qcmsg.cpp, line 458]:
  .... aborting the run
这是为什么呢？服务器剩余内存是250G

wzkchem5

魔卡少女晗 发表于 2021-11-29 07:52
4.2版本，老师，我算的是反应物产物和过渡态，所以原子数目种类是相同的，只有反应物结构可以用这个关键 ...

先换成5.0.1试试再说

魔卡少女晗

PoorChaos 发表于 2021-11-23 16:33
我替换了你的关键词，换成了正常计算的关键词，其他不变，计算也错误

我的体系，有的算的出来有的不行，不知道为啥

魔卡少女晗

LittlePupil 发表于 2021-11-23 23:43
这个自旋多重度真的没问题吗...

GAUSSIAN可以算出来的

魔卡少女晗

wzkchem5 发表于 2021-11-23 16:28
用的是什么版本的orca？这个看起来不像内存问题，但是我如果前面关键词按你写的，分子坐标用水分子的话，重 ...

4.2版本，老师，我算的是反应物产物和过渡态，所以原子数目种类是相同的，只有反应物结构可以用这个关键词算出来

LittlePupil

这个自旋多重度真的没问题吗...

PoorChaos

魔卡少女晗 发表于 2021-11-23 16:20
应该是合理的，是gaussian优化好的结构

我替换了你的关键词，换成了正常计算的关键词，其他不变，计算也错误

wzkchem5

用的是什么版本的orca？这个看起来不像内存问题，但是我如果前面关键词按你写的，分子坐标用水分子的话，重复不出来这个问题。在你的机子上如果同样的关键词算小一些的分子，也会有这个问题吗？

魔卡少女晗

PoorChaos 发表于 2021-11-23 15:12
结构合理吗

应该是合理的，是gaussian优化好的结构

PoorChaos

结构合理吗