计算化学公社

标题: 求助ORCA计算H的单点能时报错 [打印本页]

作者
Author: GeminiWang 时间: 2024-3-13 21:10
标题: 求助ORCA计算H的单点能时报错
各位老师好
ORCA输入文件
inp
%maxcore  1000
%pal nprocs  4 end
* xyz 0 2
  H       0.0000000       0.0000000       0.0000000
*

方法文件
%compound "CBS_rohf.cmp"

variable FinalEnergy;

New_Step
!CCSD(T)-F12/RI cc-pVTZ-F12 cc-pVTZ-F12-CABS cc-pVQZ/C tightSCF noautostart nopop rijdx
%mdci UseQROs true end
Step_End

read FinalEnergy = MDCI_Total_Energy[1];

End

报错
Error (ORCA_MDCI): Number of processes (4) in parallel calculation exceeds number of pairs (0)

ORCA finished by error termination in MDCI
Calling Command: mpirun -np 4  /home/test/orca/orca_mdci_mpi reac2_orca_Compound_1.mdciinp.tmp reac2_orca_Compound_1
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run

[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run

恳请各位老师指教，调试一天了，没什么效果，相同的计算方法和脚本提交H2以及其他的分子都没有问题。按照其他帖子，改变内存和核心数，并没有作用。
服务器64核心，350GB

作者
Author: zjxitcc 时间: 2024-3-13 21:24
本帖最后由 zjxitcc 于 2024-3-13 21:25 编辑

氢原子只有1个电子，没有电子相关能（也可以说电子相关能为0），无法使用电子相关方法进行计算，至少2个电子才有电子相关。有的程序会判断一下电子数，输出0的相关能，有的程序直接就报错（这种一般是认为用户缺乏基本知识，需要停下来想一想；或者假设用户懂基本原理，不会去算氢原子）。此时取ROHF能量=Full CI能量，取ROHF即可。

作者
Author: hebrewsnabla 时间: 2024-3-13 21:46
参见orca官方的说明 https://orcaforum.kofo.mpg.de/viewtopic.php?f=18&t=7913

作者
Author: GeminiWang 时间: 2024-3-13 22:05

zjxitcc 发表于 2024-3-13 21:24
氢原子只有1个电子，没有电子相关能（也可以说电子相关能为0），无法使用电子相关方法进行计算，至少2个电 ...

谢谢老师的解答，问题已经解决了。我确实是门外汉，缺乏很多基础知识，甚至不太会linux基础操作。以后多多在社区里学习，再次感谢。

作者
Author: GeminiWang 时间: 2024-3-13 22:05

hebrewsnabla 发表于 2024-3-13 21:46
参见orca官方的说明 https://orcaforum.kofo.mpg.de/viewtopic.php?f=18&t=7913

感谢您的解答！

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