用ORCA计算CCSD(T)/CC-pvqz出错

adsdfq · 发表于 Post on 2017-10-23 15:12:01

求助贴，大家好！上次参加sob老量子化学培训班，sob老师推荐用orca计算能量比较快，所有我回来后准备采用orca算能量选取的 CCSD(T)/CC-PVQZ和CCSD(T)/CC-PVTZ方法.对于CCSD(T)/CC-PVTZ, 我所有的计算都能正常结束不报错. 但是用 CCSD(T)/CC-PVQZ方法计算的时候会出现问题，具体如下：

这是我的输入文件
! cc-pvqz CCSD(T) TightSCF
%pal nprocs 16 end
%maxcore 2500

* xyz 0 2
H -2.31345100 0.44361700 0.92904500
H -2.31351500 0.44358900 -0.92893700
C -0.69392900 -0.76268300 0.00001700
H -0.59099000 -1.38949300 -0.88579600
H -0.59094600 -1.38949400 0.88582400
N 0.60459400 0.06523500 -0.00001400
O 1.63100600 -0.57160700 -0.00001600
O 0.51664800 1.26879000 -0.00003500
C -1.90681800 0.07229500 0.00004600
*
用这样的输入文件，采用CCSD(T)/CC-PVQZ计算小分子的时候也( 例如C2H4,C2H3, NO2)也都能正常结束,但是计算大分子（比如 C2H4NO2, C2H4ONO2, CH2OCH2NO 和 NO2 + C2H4反应的过渡态)就报错. 报错信息如下：
Number of amplitudes to be optimized ... 43430812
Iter E(tot) E(Corr) Delta-E Residual Time <S|S>**1/2
0 --------------------------------------------------------------------------
mpirun noticed that process rank 4 with PID 0 on node localhost exited on signal 9 (Killed).
--------------------------------------------------------------------------
ORCA finished by error termination in MDCI
Calling Command: mpirun -np 16 --allow-run-as-root /home/orca/orca_4_0_0_linux_x86-64/orca_mdci_mpi C2H4NO2-TS-CH3CHNO2.mdciinp.tmp
[file orca_tools/qcmsg.cpp, line 432]:
.... aborting the run

请问要怎么解决这个问题？？？

谢谢

sobereva · 发表于 Post on 2017-10-23 15:25:26

maxcore设小点再试，或者适当降低并行核数。
注意此级别计算超级昂贵。
用到QZ级别时建议用F12。写CCSD(T)-F12 cc-pvdz-F12 cc-pvdz-F12-CABS，精度和CCSD(T)/cc-pVQZ基本一样，耗时能省很多。

adsdfq · 发表于 Post on 2017-10-23 15:27:22

sobereva 发表于 2017-10-23 15:25
maxcore设小点再试，或者适当降低并行核数。
注意此级别计算超级昂贵。
用到QZ级别时建议用F12。写CCSD(T ...

谢谢老师

圣人惠 · 发表于 Post on 2017-11-14 08:29:02

sobereva 发表于 2017-10-23 15:25
maxcore设小点再试，或者适当降低并行核数。
注意此级别计算超级昂贵。
用到QZ级别时建议用F12。写CCSD(T ...

老师，那如果涉及到有赝势的，F12还能用吗

尚艳磊 · 发表于 Post on 2017-12-26 17:28:37

sobereva 发表于 2017-10-23 15:25
maxcore设小点再试，或者适当降低并行核数。
注意此级别计算超级昂贵。
用到QZ级别时建议用F12。写CCSD(T ...

sob老师，之前有建议用DLPNO-CCSD(T)来做，这样的方法和F12是一样的吗？
今天刚安装好ORCA，还有很多没有明白的问题。希望老师答疑解惑。

谢谢

sobereva · 发表于 Post on 2017-12-26 17:35:47

尚艳磊发表于 2017-12-26 17:28
sob老师，之前有建议用DLPNO-CCSD(T)来做，这样的方法和F12是一样的吗？
今天刚安装好ORCA，还有很多没 ...

DLPNO和F12是两码事
前者是利用定域轨道加速较大体系耦合簇计算的方法，原子数目越多效果越明显；后者是显式相关计算，使得用比较小的F12基组就可以达到较大基组的效果（比如cc-pVDZ-F12的精度和cc-pVQZ基本持平），从而节约耗时。

尚艳磊 · 发表于 Post on 2017-12-26 21:15:21

sobereva 发表于 2017-12-26 17:35
DLPNO和F12是两码事
前者是利用定域轨道加速较大体系耦合簇计算的方法，原子数目越多效果越明显；后者 ...

sob老师，我试过这个简单的测试，可以运行。
! CCSD(T) cc-pvdz
# Test for NO2

*xyz 0 2
N                0.00000000 0.00000000 0.31500700
O                0.00000000 1.09057300 -0.13781600
O                0.00000000 -1.09057300 -0.13781600
*
但是，试了一下您说的方法
! CCSD(T)-F12 cc-pvdz-F12 cc-pvdz-F12-CABS
# Test for NO2

*xyz 0 2
N                0.00000000 0.00000000 0.31500700
O                0.00000000 1.09057300 -0.13781600
O                0.00000000 -1.09057300 -0.13781600
*
就会出错了
Error (ORCA_MAIN): ... aborting the run
终止的原因是什么啊？

sobereva · 发表于 Post on 2017-12-26 23:30:44

尚艳磊发表于 2017-12-26 21:15
sob老师，我试过这个简单的测试，可以运行。
! CCSD(T) cc-pvdz
# Test for NO2

用! CCSD(T)-F12/RI cc-pVDZ-F12 cc-pVDZ-F12-CABS cc-pVTZ/C
不带RI的时候不支持UHF参考态。开RI耗时也能降低不少

尚艳磊 · 发表于 Post on 2018-1-2 11:47:32

本帖最后由尚艳磊于 2018-1-2 13:50 编辑

谢谢sob老师~

尚艳磊 · 发表于 Post on 2018-1-4 21:00:15

sobereva 发表于 2017-10-23 15:25
maxcore设小点再试，或者适当降低并行核数。
注意此级别计算超级昂贵。
用到QZ级别时建议用F12。写CCSD(T ...

sob老师，我在集群上试了简单的测试。4，8，16的核都可以算完，44核一个节点出错了。
但是所有的文件输出里面都有下面的情况：
Failed to register memory region (MR):

Hostname: n0105
Address: 60001000
Length: 4194304
Error: Cannot allocate memory
--------------------------------------------------------------------------
--------------------------------------------------------------------------
Open MPI has detected that there are UD-capable Verbs devices on your
system, but none of them were able to be setup properly. This may
indicate a problem on this system.

You job will continue, but Open MPI will ignore the "ud" oob component
in this run.

Hostname: n0105

请问是怎么回事啊

sobereva · 发表于 Post on 2018-1-5 03:40:44

尚艳磊发表于 2018-1-4 21:00
sob老师，我在集群上试了简单的测试。4，8，16的核都可以算完，44核一个节点出错了。
但是所有的文件输 ...

内存不够
maxcore设的是每个并行进程用的内存量，而不是总量。如果设了N核并行，机子里空余物理内存为M，则maxcore必须明显小于M/N才能确保运行

尚艳磊 · 发表于 Post on 2018-1-5 15:09:37

sobereva 发表于 2018-1-5 03:40
内存不够
maxcore设的是每个并行进程用的内存量，而不是总量。如果设了N核并行，机子里空余物理内存为M ...

每个节点是64GB的物理内存，我设置了maxcore=1400，还是出现了如上的提示。sob老师，是不是openmpi没装好啊？
感觉ORCA占用内存太大了，用绝对路径调用算10个原子的CCSD/cc-pVTZ都提示内存不足。

sobereva · 发表于 Post on 2018-1-6 07:10:39

尚艳磊发表于 2018-1-5 15:09
每个节点是64GB的物理内存，我设置了maxcore=1400，还是出现了如上的提示。sob老师，是不是openmpi没装好 ...

44核才64GB内存太小了
1.4GB*44=61.6GB，本身ORCA实际对内存使用量往往会高于maxcore一些，再加上系统和后台任务本身也占内存，明显超了。不是ORCA太占内存，而是计算设定不合理、本身机子内存也太小

尚艳磊 · 发表于 Post on 2018-1-7 18:40:04

sobereva 发表于 2018-1-6 07:10
44核才64GB内存太小了
1.4GB*44=61.6GB，本身ORCA实际对内存使用量往往会高于maxcore一些，再加上系统和 ...

sob老师，RI加速只支持RHF体系吗？我选择了RHF，但还是出现了问题。我的输入文件如下：
|  1> ! RHF RI-CCSD(T) cc-pVTZ cc-pVTZ/C
|  2> # Test for NO2
|  3>
|  4> *xyz 0 2
|  5>  N                0.00000000 0.00000000 0.31500700
|  6>  O                0.00000000 1.09057300 -0.13781600
|  7>  O                0.00000000 -1.09057300 -0.13781600
|  8> **
而且它提示的也是这样：WARNING: your system is open-shell and RHF/RKS was chosen
  ===> : WILL SWITCH to UHF/UKS

但是结果却说：[file orca_mdci/mdci_uhf_sigma.cpp, line 3388]: Error (ORCA_MDCI): K(C)-RI not implemented in UHF branch
请问老师，这该咋么解决啊？谢谢啦

sobereva · 发表于 Post on 2018-1-7 18:54:25

尚艳磊发表于 2018-1-7 18:40
sob老师，RI加速只支持RHF体系吗？我选择了RHF，但还是出现了问题。我的输入文件如下：
| 1> ! RHF RI- ...

CCSD(T)别写RI-，速度不会有提升
只有耦合簇结合F12的时候才有利用RI的意义

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] 用ORCA计算CCSD(T)/CC-pvqz出错

浏览过的版块