计算化学公社

标题: orca算单点能,SCF收敛后,MP2部分计算停止 [打印本页]

作者
Author:
魔卡少女晗    时间: 2021-11-23 08:40
标题: 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


作者
Author:
PoorChaos    时间: 2021-11-23 15:12
结构合理吗
作者
Author:
魔卡少女晗    时间: 2021-11-23 16:20
PoorChaos 发表于 2021-11-23 15:12
结构合理吗

应该是合理的,是gaussian优化好的结构
作者
Author:
wzkchem5    时间: 2021-11-23 16:28
用的是什么版本的orca?这个看起来不像内存问题,但是我如果前面关键词按你写的,分子坐标用水分子的话,重复不出来这个问题。在你的机子上如果同样的关键词算小一些的分子,也会有这个问题吗?
作者
Author:
PoorChaos    时间: 2021-11-23 16:33
魔卡少女晗 发表于 2021-11-23 16:20
应该是合理的,是gaussian优化好的结构

我替换了你的关键词,换成了正常计算的关键词,其他不变,计算也错误
作者
Author:
LittlePupil    时间: 2021-11-23 23:43
这个自旋多重度真的没问题吗...
作者
Author:
魔卡少女晗    时间: 2021-11-29 14:52
wzkchem5 发表于 2021-11-23 16:28
用的是什么版本的orca?这个看起来不像内存问题,但是我如果前面关键词按你写的,分子坐标用水分子的话,重 ...

4.2版本,老师,我算的是反应物产物和过渡态,所以原子数目种类是相同的,只有反应物结构可以用这个关键词算出来
作者
Author:
魔卡少女晗    时间: 2021-11-29 14:53
LittlePupil 发表于 2021-11-23 23:43
这个自旋多重度真的没问题吗...

GAUSSIAN可以算出来的
作者
Author:
魔卡少女晗    时间: 2021-11-29 14:54
PoorChaos 发表于 2021-11-23 16:33
我替换了你的关键词,换成了正常计算的关键词,其他不变,计算也错误

我的体系,有的算的出来有的不行,不知道为啥
作者
Author:
wzkchem5    时间: 2021-11-29 16:22
魔卡少女晗 发表于 2021-11-29 07:52
4.2版本,老师,我算的是反应物产物和过渡态,所以原子数目种类是相同的,只有反应物结构可以用这个关键 ...

先换成5.0.1试试再说




欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3