计算化学公社

 找回密码 Forget password
 注册 Register
Views: 9385|回复 Reply: 8
打印 Print 上一主题 Last thread 下一主题 Next thread

[ORCA] ORCA双杂化泛函计算单点能报错RIFLAG[0]!=1 (Cholesky) is not allowed here

[复制链接 Copy URL]

46

帖子

0

威望

405

eV
积分
451

Level 3 能力者

本帖最后由 魔卡少女晗 于 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


928

帖子

1

威望

8240

eV
积分
9188

Level 6 (一方通行)

2#
发表于 Post on 2021-7-9 12:17:55 | 只看该作者 Only view this author
这句不是报错,报错是下面的RIFLAG[0]!=1 (Cholesky) is not allowed here.

46

帖子

0

威望

405

eV
积分
451

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2021-7-9 13:55:34 | 只看该作者 Only view this author
hebrewsnabla 发表于 2021-7-9 12:17
这句不是报错,报错是下面的RIFLAG[0]!=1 (Cholesky) is not allowed here.

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

928

帖子

1

威望

8240

eV
积分
9188

Level 6 (一方通行)

4#
发表于 Post on 2021-7-9 14:20:47 | 只看该作者 Only view this author
本帖最后由 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,这样更方便

46

帖子

0

威望

405

eV
积分
451

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2021-7-9 16:32:24 | 只看该作者 Only view this author
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
然后算单点能想要考虑的是两个分子间弱相互作用,所以加的弥散。

928

帖子

1

威望

8240

eV
积分
9188

Level 6 (一方通行)

6#
发表于 Post on 2021-7-9 16:47:55 | 只看该作者 Only view this author
本帖最后由 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已经足够弥散了,不加弥散也行。而且有现成的辅助基组,更可靠。

46

帖子

0

威望

405

eV
积分
451

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2021-7-9 16:53:17 | 只看该作者 Only view this author
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相同的泛函基组,但是像您说的辅助基组不一样,会影响想要的结果吗?

928

帖子

1

威望

8240

eV
积分
9188

Level 6 (一方通行)

8#
发表于 Post on 2021-7-9 16:56:21 | 只看该作者 Only view this author
魔卡少女晗 发表于 2021-7-9 16:53
嗯嗯好的老师我回去试试。还有一个问题,就是我想算A、B之间的结合能,所以算了A、B、AB三个文件,如果三 ...

辅助基组也要一样的

46

帖子

0

威望

405

eV
积分
451

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2021-7-9 18:51:47 | 只看该作者 Only view this author
hebrewsnabla 发表于 2021-7-9 16:56
辅助基组也要一样的

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

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-1-25 17:58 , Processed in 3.242551 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list