计算化学公社

标题: 求助Potentially linear dependencies in the auxiliary basis (\J) [打印本页]

作者
Author: 天道啊啊 时间: 2020-2-2 19:01
标题: 求助Potentially linear dependencies in the auxiliary basis (\J)
各位老师，您好！
我在使用ORCA 4.2.0计算弱相互作用，目前在两个分子结合的：构型（构型由B97-3c AutoAux RIJCOSX优化得到）的单点能计算遇到了orca 的报错。
输入文件如下：
! wB97M-V ma-def2-TZVP AutoAux GCP(DFT/TZ) TIGHTSCF grid4 gridx4 RIJCOSX noautostart miniprint nopop %cpcm smd true SMDsolvent "water" end %maxcore 2500
%pal nprocs 5 end  %scf  SOSCFmaxit 12 directresetfreq 1 end  往下是130个原子的坐标。

输出文件中提示：
WARNING! Potentially linear dependencies in the auxiliary basis (\J).
===>> Running eigenvalue decompsition instead (RIFLAG[1]<0).
      This is not supported with all the options/modules.
      If the calculation crashes, please re-try with a different auxiliary basis.

如果无视这个警告，会有两个结果：
1. 内存莫名报错 [file orca_tools/qcmem.cpp, line 931, Process 2]:  OUT OF MEMORY ERROR! 我的电脑32G内存，按照5X2500MB 分配核心和内存会报错。就算我改成1核X1000MB仍然报错
2.又产生一条警告WARNING (ORCA_TOOLS): LoewdinSM12()has at least 5 eigenvalues smaller than 1.000e-07 ( 0.000e+00  0.000e+00  0.000e+00  0.000e+00  0.000e+00 )
  然后SCF 能量振荡，无法收敛最终报错。按照http://bbs.keinsci.com/thread-4570-1-1.html的帖子试了试各种解决SCF不收敛的方法都不行。

我自己做了以下尝试：
1.更换泛函  将wB97M-V换为PW6B95-D3BJ，M06-2X-D3，B3LYP-D3BJ。依然报错
2.更换基组  将ma-def2-TZVP  换为def2-TZVPD。仍然报错
3.不用AutoAux，使用不太匹配的辅助基组。ma-def2-TZVP或def2-TZVPD搭配def2-TZVP/c都没有问题。
4.使用更昂贵的基组  aug-cc-pVnZ + aug-cc-pVnZ/C没有问题。

我觉得问题在AutoAux上。我自己思考了一下解决方案
1. 构型问题？B97-3c AutoAux RIJCOSX 优化得到的构型有问题？换B3LYP-gCP-D3/6-31G* 作构型优化，再用wB97M-V ma-def2-TZVP AutoAux GCP(DFT/TZ)计算单点能试试？
2. 选用不太合适的辅助基组做计算？ ma-def2-TZVP或def2-TZVPD 搭配def2-TZVP/c？这样会不会误差比较大？
3.使用aug-cc-pVTZ搭配aug-cc-pVTZ/C？
4. 强攻SCF 能量振荡的问题？
请各位老师指正！感谢万分！

作者
Author: wuzhiyi 时间: 2020-2-3 02:26
我觉得这问题你发orca论坛更有用
我觉得是autoaux的问题告诉开发者说不定可以在下个版本修复

作者
Author: sobereva 时间: 2020-2-3 04:39
都3-zeta加弥散了，根本就不需要再用gCP。必须去掉。

改成def2-QZVP，不仅有现成的辅助基组（比用autoaux构建的划算得多），线性依赖问题还没有ma-def-TZVP强。怎么就不想着用def2-QZVP。下文里明确体现def2-QZVP在ORCA里对于DFT算中等体系毫无压力
大体系弱相互作用计算的解决之道
http://sobereva.com/214

作者
Author: 天道啊啊 时间: 2020-2-3 10:06
本帖最后由天道啊啊于 2020-2-3 10:33 编辑

sobereva 发表于 2020-2-3 04:39
都3-zeta加弥散了，根本就不需要再用gCP。必须去掉。

改成def2-QZVP，不仅有现成的辅助基组（比用autoau ...

谢谢老师，gcp我会删除。弥散函数应该还是要加的把？您的意思是用DEF2-QZVPD？可是程序里没有和它比较搭配的辅助基组啊。Def2-TZVP，Def2-TZVPP/C和Def2-QZVPP的辅助基组比较完备，我手动输入个Def2-QZVP/C或Def2-QZVPD/C，程序提示无法识别关键词。是用DEF2-QZVPD搭配Def2-QZVPPD/C?

作者
Author: 天道啊啊 时间: 2020-2-3 10:11

wuzhiyi 发表于 2020-2-3 02:26
我觉得这问题你发orca论坛更有用
我觉得是autoaux的问题告诉开发者说不定可以在下个版本修复

谢谢您的建议

作者
Author: mutron 时间: 2020-2-3 21:30
本帖最后由 mutron 于 2020-2-3 21:51 编辑

天道啊啊发表于 2020-2-3 10:06
谢谢老师，gcp我会删除。弥散函数应该还是要加的把？您的意思是用DEF2-QZVPD？可是程序里没有和它比较搭 ...

大体系弱相互作用计算的解决之道
http://sobereva.com/214
仔细读它几次，尤其是中间用def2-QZVPP（新版没有了阉割版）取代aug-cc-pVTZ的理由

RI and auxiliary basis sets
https://sites.google.com/site/orcainputlibrary/basis-sets/ri-and-auxiliary-basis-sets
仔细阅读最下面关于autoaux那段话

仔细看看ORCA manual里面关于/C辅助基组可以有哪些选择，不要自己写程序没有的关键词

如果你有不可抗拒的理由非要用弥散函数、又不想用autoaux，可以尝试aug-cc-pVnZ/JK，当然最好用noRI做一次单点对比（noRI读取RIJCOSX的gbw节省时间）。以前做过一些小分子用这个辅助基组带来的误差很小

作者
Author: sobereva 时间: 2020-2-4 10:07

天道啊啊发表于 2020-2-3 10:06
谢谢老师，gcp我会删除。弥散函数应该还是要加的把？您的意思是用DEF2-QZVPD？可是程序里没有和它比较搭 ...

先把我博文仔细读了再问

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