ORCA 5 ,Yes

zjxitcc

wzkchem5 发表于 2021-7-2 15:01
就是基组线性相关问题吧。。。
之前这边有人有计划调整AutoAux产生的辅助基组，来缓解基组线性相关问题 ...

你是指辅助基的线性相关问题？不是原基组的线性相关问题？

Satoru

zjxitcc 发表于 2021-7-2 14:40
给ma-def2系列基组用AutoAux翻车 能否介绍一下是什么情况？

ma-def2-QZVP或者QZVPP加autoaux有时会提示辅助基组线性相关问题，双杂化基本到MP2部分必报错，wB97M-V在SCF开始前随机报错，看了半天手册没搞定怎么改...

一次wB97M-V/ma-def2-QZVP/smd RI+RIJCOSX算一个中性有机Pd(II)的过渡态单点，初猜是某斯PBE0/def2-SVP几何优化的波函数，算出来比底物高大概540 kcal/mol，怀疑是没收敛到基态；改用slowconv以后收敛到了比底物低90 kcal/mol左右。底物往前几步相对能量都很正常，感觉SCF没问题。
后来改用了maug-cc-pvqz和-PP来算，主族用预设的/JK基组，Pd用autoaux，算出来就很正常。两个方法算底物绝对能量大概差0.03 Ha，所以感觉是之前autoaux收敛到了奇怪的能量过低的解。

hdhxx123

Satoru 发表于 2021-7-2 15:13
ma-def2-QZVP或者QZVPP加autoaux有时会提示辅助基组线性相关问题，双杂化基本到MP2部分必报错，wB97M-V在 ...

感觉这是操作的问题，如果要读波函数的话一般需要是同zeta的，基组SVP和QZVP是2zeta到4zeta，波函数质量差得有点太多了，按照你的这个orca计算级别，某斯应该是pbe0/def2tzvp起步

Satoru

hdhxx123 发表于 2021-7-2 15:16
感觉这是操作的问题，如果要读波函数的话一般需要是同zeta的，基组SVP和QZVP是2zeta到4zeta，波函数质量 ...

比较懒x之前还没翻过车
针对这个例子，某斯wB97X-D/def2-TZVP波函数初猜我也尝试过，不记得是收敛到哪个解了，反正是不对...

zjxitcc

Satoru 发表于 2021-7-2 15:19
比较懒x之前还没翻过车
针对这个例子，某斯wB97X-D/def2-TZVP波函数初猜我也尝试过，不记得是收敛到哪个 ...

辅助基线性相关问题不好解决。SCF收敛性问题、收敛到正确的解上，这十分好解决。

wzkchem5

zjxitcc 发表于 2021-7-2 08:08
你是指辅助基的线性相关问题？不是原基组的线性相关问题？

对，这个是autoaux的主要问题，容易有线性相关，导致SCF出现variational instability。
我不知道ORCA 5有没有把他们说的改进加进去，如果加了的话，ORCA 5的autoaux会有改善

Satoru

zjxitcc 发表于 2021-7-2 15:24
辅助基线性相关问题不好解决。SCF收敛性问题、收敛到正确的解上，这十分好解决。

qwq所以收敛到比（理想的）正确解低90 kcal/mol左右也可以救回来嘛

wzkchem5

Satoru 发表于 2021-7-2 08:32
qwq所以收敛到比（理想的）正确解低90 kcal/mol左右也可以救回来嘛

这个属于线性相关问题。。。他应该是说除了基组线性相关以外的SCF收敛问题好解决

hdhxx123

wzkchem5 发表于 2021-7-2 15:41
这个属于线性相关问题。。。他应该是说除了基组线性相关以外的SCF收敛问题好解决

nearly all done

Satoru

wzkchem5 发表于 2021-7-2 15:41
这个属于线性相关问题。。。他应该是说除了基组线性相关以外的SCF收敛问题好解决

感谢orz这两天去试试拿5.0重算下

biogon

wzkchem5 发表于 2021-7-2 01:50
提醒一下大家：
虽然ORCA对学术用户是免费的，但是原则上讲，只有从官网（https://orcaforum.kofo.mpg.d ...

我算DLPNO-STEOM出了如下的错误是bug吗？
------------------------------------------------------------------------------
                           ORCA POPULATION ANALYSIS
------------------------------------------------------------------------------
Input electron density              ... 1.mdcip
BaseName (.gbw .S,...)              ... 1
[pc-AS-2023US-TR4:200440] *** Process received signal ***
[pc-AS-2023US-TR4:200440] Signal: Floating point exception (8)
[pc-AS-2023US-TR4:200440] Signal code: Integer divide-by-zero (1)
[pc-AS-2023US-TR4:200440] Failing at address: 0xba0de2
[pc-AS-2023US-TR4:200438] *** Process received signal ***
[pc-AS-2023US-TR4:200438] Signal: Floating point exception (8)
[pc-AS-2023US-TR4:200438] Signal code: Integer divide-by-zero (1)
[pc-AS-2023US-TR4:200438] Failing at address: 0xba0de2
[pc-AS-2023US-TR4:200438] [ 0] /lib/x86_64-linux-gnu/libc.so.6(+0x3efd0)[0x155112408fd0]
[pc-AS-2023US-TR4:200438] [ 1] [pc-AS-2023US-TR4:200441] *** Process received signal ***

wzkchem5

biogon 发表于 2021-7-2 09:28
我算DLPNO-STEOM出了如下的错误是bug吗？
---------------------------------------------------------- ...

建议发到orca论坛的bugs版面
顺便你可以测一下，在你的机子上其他类似的方法会不会在population那里报错，比如EOM-CCSD，这样方便我们锁定问题所在

biogon

wzkchem5 发表于 2021-7-2 16:37
建议发到orca论坛的bugs版面
顺便你可以测一下，在你的机子上其他类似的方法会不会在population那里报错 ...

这个体系四十多原子，没DLPNO的恐怕搞不动吧

biogon

wzkchem5 发表于 2021-7-2 16:37
建议发到orca论坛的bugs版面
顺便你可以测一下，在你的机子上其他类似的方法会不会在population那里报错 ...

我刚才试了跑别的体系，给64核scf迭代的时候就崩了，32核就没事

wzkchem5

biogon 发表于 2021-7-2 10:09
我刚才试了跑别的体系，给64核scf迭代的时候就崩了，32核就没事

你的服务器有多少个核？有跨节点跑吗？
比较新的OpenMPI好像要求任务用的核数不能超过服务器的物理核数。比较老的版本可以超过物理核数，甚至超过逻辑核数可能都行