计算化学公社

标题: ORCA双杂化泛函计算单点能报错RIFLAG[0]!=1 (Cholesky) is not allowed here [打印本页]

作者
Author: 魔卡少女晗 时间: 2021-7-9 11:02
标题: ORCA双杂化泛函计算单点能报错RIFLAG[0]!=1 (Cholesky) is not allowed here
本帖最后由魔卡少女晗于 2021-7-9 13:37 编辑

SCF可以收敛，体系里面含有Ag但是计算单点能的时候提示RIFLAG[0]!=1 (Cholesky) is not allowed here，各位老师能帮我看看是怎么回事嘛输入文件是这个
! PWPB95 D3 ma-def2-QZVPP autoaux RIJCOSX grid6 gridx6 miniprint nopop
%maxcore  3000
%pal nprocs 8 end
* xyz 0 1
Ag -1.51111127 0.82222227 0.00000000
Cl -3.77721144 2.07700137 0.87785907
Cl -3.91889545 5.09636760 -0.87960407
Cl -6.80328068 3.72533149 0.63550005
Cl -5.47606556 2.28275339 -2.19309517
Al -4.97543453 3.30993646 -0.38546903
Ag    3.14286109 3.96055651 -0.66956405
Cl    5.40896128 2.70577741 -1.54742312
Cl    5.55064526 -0.31358881 0.21004002
Cl    8.43503047 1.05744729 -1.30506410
Cl    7.10781541 2.50002540 1.52353111
Al    6.60718436 1.47284232 -0.28409502
Cl    2.45275904 4.19040253 -3.31227025

，然后输出文件里面是这样报错的
Total SCF time: 0 days 2 hours 7 min 53 sec

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

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

Freezing NCore=176 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 4]: 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.

[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 5]: ERROR ORCA_MP2: RIFLAG[0]!=1 (Cholesky) is not allowed here.

[file orca_mp2/mp2rigrad.cpp, line 790, Process 6]: 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 7]: ERROR ORCA_MP2: RIFLAG[0]!=1 (Cholesky) is not allowed here.

ORCA finished by error termination in MP2
Calling Command: mpirun -np 8  /root/orca/orca_mp2_mpi px.mp2inp.tmp px
[file orca_tools/qcmsg.cpp, line 458]:
  .... aborting the run

作者
Author: hebrewsnabla 时间: 2021-7-9 12:17
这句不是报错，报错是下面的RIFLAG[0]!=1 (Cholesky) is not allowed here.

作者
Author: 魔卡少女晗 时间: 2021-7-9 13:55

hebrewsnabla 发表于 2021-7-9 12:17
这句不是报错，报错是下面的RIFLAG[0]!=1 (Cholesky) is not allowed here.

我改好了，您知道报错这个是因为啥吗

作者
Author: hebrewsnabla 时间: 2021-7-9 14:20
本帖最后由 hebrewsnabla 于 2021-7-9 14:24 编辑

魔卡少女晗发表于 2021-7-9 13:55
我改好了，您知道报错这个是因为啥吗

orca 论坛有人问过了 https://orcaforum.kofo.mpg.de/vi ... hilit=riflag#p31139

不要用AutoAUx，用def2/J aug-cc-pVQZ/C。ma-def2-QZVPP没有太合适的辅助基，只能凑合一下。

看起来这个体系不需要弥散啊，如果用def2-QZVPP基组的话，辅助基用def2/J def2-QZVPP/C，这样更方便

作者
Author: 魔卡少女晗 时间: 2021-7-9 16:32

hebrewsnabla 发表于 2021-7-9 14:20
orca 论坛有人问过了 https://orcaforum.kofo.mpg.de/vi ... hilit=riflag#p31139

不要用AutoAUx，用 ...

老师，加上您给的这个辅助基组之后，提示The basis set was either not assigned or not available for this element - Aborting the run

然后算单点能想要考虑的是两个分子间弱相互作用，所以加的弥散。

作者
Author: hebrewsnabla 时间: 2021-7-9 16:47
本帖最后由 hebrewsnabla 于 2021-7-9 16:51 编辑

魔卡少女晗发表于 2021-7-9 16:32
老师，加上您给的这个辅助基组之后，提示The basis set was either not assigned or not available for t ...

给Au用aug-cc-pvqz-pp/C，Cl用aug-cc-pvqz/C，或者统一用def2-qzvppd/c。但是这些都是借用辅助基组的做法，我不能保证靠谱。

或者你把基组也换掉，给Au用aug-cc-pvqz-pp，Cl用aug-cc-pvqz，或者统一用def2-qzvppd，这样就配套了。

def2-qzvpp已经足够弥散了，不加弥散也行。而且有现成的辅助基组，更可靠。

作者
Author: 魔卡少女晗 时间: 2021-7-9 16:53

hebrewsnabla 发表于 2021-7-9 16:47
给Au用aug-cc-pvqz-pp/C，Cl用aug-cc-pvqz/C，或者统一用def2-qzvppd/c。但是这些都是借用辅助基组的做法 ...

嗯嗯好的老师我回去试试。还有一个问题，就是我想算A、B之间的结合能，所以算了A、B、AB三个文件，如果三个文件用ORCA相同的泛函基组，但是像您说的辅助基组不一样，会影响想要的结果吗？

作者
Author: hebrewsnabla 时间: 2021-7-9 16:56

魔卡少女晗发表于 2021-7-9 16:53
嗯嗯好的老师我回去试试。还有一个问题，就是我想算A、B之间的结合能，所以算了A、B、AB三个文件，如果三 ...

辅助基组也要一样的

作者
Author: 魔卡少女晗 时间: 2021-7-9 18:51

hebrewsnabla 发表于 2021-7-9 16:56
辅助基组也要一样的

谢谢老师！更换了辅助基组后可以继续算了！

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