求助：MDCI错误，显示无法打开GBW文件，Cannot open GBW file

路飞07

本人初入量化，目前遇到一个不知道是系统bug的问题还是本人操作有误，特来求助
在算一个非常简单的碳原子的sp，有的价态以及电子数下的就能算出来，但有的任务就会失败，重复提交也还是一样的失败，

错误代码：
（之前都正常，直到这一步）

Localizing the valence orbitals
[file orca_mdci/mdci_util.cpp, line 960]: Error (ORCA_MDCI): Cannot open GBW file: //w0/tmp/slurm_br466165.43818652/43818652/example15.loc


ORCA finished by error termination in MDCI
Calling Command: /cvmfs/restricted.hpc.rwth.de/Linux/RH8/x86_64/ISV/ORCA/5.0.4-gompi-2022a/bin/orca_mdci //w0/tmp/slurm_br466165.43818652/43818652/example15.mdciinp.tmp
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run




详细原文件请看附件
已参考过本站一些求助帖，例如改为提交单核、单任务、调整内存、allow RHF、SCF不收敛的情况下改为SOSCF等灯，均不奏效，希望大佬能帮忙解决一下，万分感谢！

以及，是否存在一些不稳定的或者很奇怪的东西，他本身就不能算sp或者不能进行一些计算？？因此无论什么情况 不管怎么调惨，都不会得出结果？

感谢大佬指点！！

zjxitcc

这与MDCI、gbw都没有关系。你算的是+3价的孤立C原子，体系只有3个电子。然后用的是post-HF方法，默认冻结1s^2芯电子，所以体系仅剩1个价电子可以用于电子相关计算，然而至少需要2个电子才能有电子相关，因此此时post-HF电子相关能为0，不需要算（算了程序报错，也是可以理解的）。

要想算这个，你有两种办法：
（1）不做post-HF计算，就用HF能量。
（2）做post-HF计算，这需要添加关键词指定 不冻结芯电子，这样体系就有3个电子参与post-HF计算。要注意的是，如果你有其他体系与这个体系进行比较，计算相对能量，需要保持一致的冻芯设置。

MrMr浩

zjxitcc 发表于 2024-2-26 13:06
这与MDCI、gbw都没有关系。你算的是+3价的孤立C原子，体系只有3个电子。然后用的是post-HF方法，默认冻结1s ...

老师您好，我也遇到了类似的问题，我想算的是一个+3价的Al原子，也设置了不冻结芯电子，但是还是出现了一样的报错，附件中是输出文件，想向您请教这是怎么回事呢，谢谢！！

zjxitcc

MrMr浩 发表于 2024-8-15 14:51
老师您好，我也遇到了类似的问题，我想算的是一个+3价的Al原子，也设置了不冻结芯电子，但是还是出现了一 ...

除了冻芯问题外，你俩报错看似报错有点像，实际原因不同。你已经设置不冻芯了，你的报错估计是ORCA的bug，可以去ORCA论坛反映一下。

针对你的问题，目前临时解决办法是在inp文件里写上UHF（就是! UHF NoFrozenCore 其他关键词），你这体系是闭壳层的，原本用默认的RHF参考波函数就可以，更省时间；但由于ORCA可能处理这个简单体系有点问题，你暂时用UHF也可以得到相同的结果，这可以在输出文件中看到<S**2>=0，仍然是闭壳层的。等以后ORCA新版本解决了这个问题，也会算出一样的结果。另外，不用写TSCALEPNO_CORE 0.01。

一般来说太小的体系直接用CCSD(T)计算即可，不需要DLPNO近似；但有的人计算复合物用了DLPNO近似，再算其中的片段时为了保持计算级别一致，不得不使用DLPNO近似。

MrMr浩

zjxitcc 发表于 2024-8-22 19:38
除了冻芯问题外，你俩报错看似报错有点像，实际原因不同。你已经设置不冻芯了，你的报错估计是ORCA的bug ...

收到，多谢老师耐心解答！！