ORCA DLPNO-CCSD(T)单点报错求助

量子极限

我用的版本是ORCA 4.2.1，ORCA输入文件为参考文献提供的orca输入文件修改的，我只修改了%maxcore 6000 和%pal nprocs  4 其他没有改动。服务器是28核 56G内存.
这是我的输入文件：
! DLPNO-CCSD(T) verytightscf def2-tzvp def2/JK def2-tzvp/C RIJCOSX gridx7 LED
%maxcore  6000
%pal nprocs   4 end
%basis aux "def2/jk" end
%mdci
TCutPairs 1e-5
PrintLevel 3
End
* xyz   0   1
坐标
*

报错信息如下：
--------------------------------
LOCAL RI TRANSFORMATION (VABPAO)
--------------------------------

Number of PAOs:    3246
Basis functions:   3246 (1290 shells)
Aux. functions:    8024 (2512 shells)

Processing maps (0.1 sec)
Average map sizes:
   Aux shells -> PAOs           2730.0
   PAOs       -> AO shells      1215.9

Calculating integrals Error in AddMatrix: Failed to add 1176 0 to Cat_MA-TS_Major_degree-thred-10_1_TS_CCSD-8.VABPAO.tmp.proc0

ORCA finished by error termination in MDCI
Calling Command: mpirun -np 4  -machinefile Cat_MA-TS_Major_degree-thred-10_1_TS_CCSD-8.nodes /public/sourcecode/orca_4_2_1_linux_x86-64_openmpi314/orca_mdci_mpi Cat_MA-TS_Major_degree-thred-10_1_TS_CCSD-8.mdciinp.tmp Cat_MA-TS_Major_degree-thred-10_1_TS_CCSD-8
[file orca_tools/qcmsg.cpp, line 458]:
  .... aborting the run



想请教一下大家，这可能是哪里出现问题，应该怎么处理
谢谢大家

菓瀛

zjxitcc 发表于 2025-5-17 10:38
默认不检验。然而，你这么提问似乎在暗示 你不知道如何在cc-pVTZ基组下检验波函数稳定性，建议阅读
http ...

好的老师 我明白啦 感谢您的耐心解答

zjxitcc

菓瀛 发表于 2025-5-16 22:29
好呢 谢谢老师

我使用cc-pVTZ的波函数去计算cc-pVQZ下的单点能 关键词这样写您看可以不

默认不检验。然而，你这么提问似乎在暗示 你不知道如何在cc-pVTZ基组下检验波函数稳定性，建议阅读
http://bbs.keinsci.com/forum.php ... 02&fromuid=2632

菓瀛

zjxitcc 发表于 2025-5-16 18:28
可以的，cc-pVTZ基组下的稳定波函数，在经过cc-pVQZ读取后，极大概率仍然是稳定的，因此无需在cc-pVQZ下 ...

好呢 谢谢老师

我使用cc-pVTZ的波函数去计算cc-pVQZ下的单点能 关键词这样写您看可以不


还有请问一下用哪个关键词可以规避在cc-pVQZ下检验波函数的稳定性呢 还是系统会默认不去检验呢 <跪谢大佬指点>

zjxitcc

菓瀛 发表于 2025-5-16 15:46
老师 我做的是非限制性开壳层计算 使用cc-pvtz可以成功算出来 cc-pvqz算不出来 (原因是机子内存不足)
按 ...

可以的，cc-pVTZ基组下的稳定波函数，在经过cc-pVQZ读取后，极大概率仍然是稳定的，因此无需在cc-pVQZ下做波函数稳定性检验。

菓瀛

zjxitcc 发表于 2025-5-10 13:23
（1）体系尺寸较小，直接用DLPNO近似，无需CIM近似。减小核数，增加每核内存。
（2）由于你做的是 非限 ...

老师 我做的是非限制性开壳层计算 使用cc-pvtz可以成功算出来 cc-pvqz算不出来 (原因是机子内存不足)
按您说的(3) 我可以使用cc-pvtz计算单点能 去检验波函数的稳定性 然后使用cc-pvqz基组读取UHF波函数进行单点能计算吗

菓瀛

zjxitcc 发表于 2025-5-10 13:23
（1）体系尺寸较小，直接用DLPNO近似，无需CIM近似。减小核数，增加每核内存。
（2）由于你做的是 非限 ...

谢谢老师的解答 我做的就是非限制性开壳层计算 我去试试您说的方法

菓瀛

wzkchem5 发表于 2025-5-10 12:23
报内存不足既可能是maxcore设得太小，也可能是maxcore设得太大。不要一味增大maxcore。
这个体系用一般 ...

谢谢老师的解答

zjxitcc

菓瀛 发表于 2025-5-10 10:29
老师 不好意思打扰您 我的体系有30个原子 计算的是含能分子HMX和NO2的氢提取反应 采用CIM-DLPNO-CCSD(T) ...

（1）体系尺寸较小，直接用DLPNO近似，无需CIM近似。减小核数，增加每核内存。
（2）由于你做的是 非限制性开壳层 计算，计算量比相应的闭壳层任务 大得多，因此建议你基组降为cc-pVTZ，如果能够成功算出，再考虑增加精度不迟。路还是一步一步走，不宜一口气吃成胖子。如果连cc-pVTZ都算不出来，说明对于目标计算而言 当前机器配置过低。
（3）由于是 非限制性开壳层 计算，建议分两步算：第一步做UHF单点计算，并检验波函数稳定性，确保获得了稳定的UHF波函数；第二步读取稳定的UHF波函数进行! UHF DLPNO-CCSD(T)计算。

wzkchem5

菓瀛 发表于 2025-5-10 11:25
谢谢老师 我之前用的是DLPNO-CCSD(T)方法计算的单点能 但是会报错内存不足的 所以我想试一下CIM 以为它不 ...

报内存不足既可能是maxcore设得太小，也可能是maxcore设得太大。不要一味增大maxcore。
这个体系用一般的机子不应该算不动，如果真算不动，可以做DLPNO-CCSD(T)-F12/cc-pVTZ-F12

菓瀛

wzkchem5 发表于 2025-5-10 11:17
这么小的分子没必要用CIM，只用DLPNO就行

谢谢老师 我之前用的是DLPNO-CCSD(T)方法计算的单点能 但是会报错内存不足的 所以我想试一下CIM 以为它不会那么吃内存

wzkchem5

菓瀛 发表于 2025-5-10 10:29
老师 不好意思打扰您 我的体系有30个原子 计算的是含能分子HMX和NO2的氢提取反应 采用CIM-DLPNO-CCSD(T) ...

这么小的分子没必要用CIM，只用DLPNO就行

菓瀛

zjxitcc 发表于 2021-7-12 10:36
如果你的体系本身比较大，适合用CIM-DLPNO-CCSD(T)方法，这是两种线性标度方法CIM与DLPNO二者的结合，其 ...

老师 不好意思打扰您 我的体系有30个原子 计算的是含能分子HMX和NO2的氢提取反应 采用CIM-DLPNO-CCSD(T)方法计算单点能 图1是输入文件 C:\Users\Administrator\Desktop\捕获2.PNG 图2是报错信息C:\Users\Administrator\Desktop\捕获.PNG 文件是输出文件 1.out (114.34 KB, 下载次数 Times of downloads: 0)

2025-5-10 10:31 上传 Uploaded

点击下载
Click to download

想请教您一下 这可能是哪里出现问题 要如何解决呢

413

shalene 发表于 2021-6-25 20:45
估计你的硬件完不成这个任务。给你个参考:2088 basis functions, maxcore设为12000MB, normalPNO。 实际内 ...

你好，请问你们当时用的机器总内存多少？
你当时跑那个作业的时候用了多少个核心数还能记得吗？

量子极限

zjxitcc 发表于 2021-7-12 10:36
如果你的体系本身比较大，适合用CIM-DLPNO-CCSD(T)方法，这是两种线性标度方法CIM与DLPNO二者的结合，其 ...

感谢两位老师的解答，目前是已经找到了合适的超算（2TB或3TB的内存，硬盘几TB大小）能算这个较大的体系。主要是参考一篇比较相关的文献，学习并按照他的方法来做的benchmark，所以方法也不太方便改，如果最后还是不行的话可能试一下CIM-DLPNO-CCSD(T)加校正的方法。