计算化学公社

标题: 量子化学qm9数据库里面的dft泛函和基组的问题 [打印本页]

作者
Author: hanshan 时间: 2023-4-13 16:49
标题: 量子化学qm9数据库里面的dft泛函和基组的问题
在量子化学qm9数据库里，http://www.doi.org/10.1038/sdata.2014.22
提到dft计算采用的是B3LYP/6-31G(2df,p)，Gaussian 09

如果同样参数用orca5来计算，是不是要改成
B3LYP/G 6-31G(2df,2p)
？
因为orca5的manual里面，只有
6-31G(2df)
6-31G(2df,2p)
似乎没有6-31G(2df,p)，程序会对6-31G(2df,p)报错的

谢谢！

作者
Author: wzkchem5 时间: 2023-4-13 17:04
可以考虑，但即便orca支持6-31G(2df,p)，结果和高斯也不严格可比，因为高斯默认用笛卡尔d函数而orca只能用球谐d函数

作者
Author: sobereva 时间: 2023-4-13 17:30
值得一提的是，6-31G(2df,p)对于DFT来说本身不是什么好的选择。把基组从2-zeta增大到3-zeta比加昂贵的f极化通常重要得多。而且重原子用到2df这么多极化的时候只给氢用p也显得严重不均衡

作者
Author: hanshan 时间: 2023-4-13 17:43

wzkchem5 发表于 2023-4-13 17:04
可以考虑，但即便orca支持6-31G(2df,p)，结果和高斯也不严格可比，因为高斯默认用笛卡尔d函数而orca只能用 ...

原来您是orca开发团队成员，谢谢您的答复！
其实我还有一个疑问，以下报错可能是什么情况呢？

第一次报错
Error (ORCA_SCFHESS): Available memory of 1000 MB is insufficient!
   ------->       Increase MaxCore to 2271 MB for OnDisk storage
                  - or -             13406 MB for InCore storage

修改MaxCore后第二次报错，
----------------------------------------------
Forming right-hand sides of CP-SCF equations    ...
----------------------------------------------
One electron integral derivatives  (SHARK)    ... done (    12.4 sec)
Transforming the overlap derivative matrices    ... done (    19.0 sec)
Making the Q(x) pseudodensities                ... done (    21.3 sec)
Adding the E*S(x)*S(y) terms to the Hessian    ... done (    12.3 sec)
Calculating energy weighted overlap derivatives  ... done (    1.2 sec)

Two electron integral derivatives (SHARK-RI)    ...
IBatch  1 (of  2):
ORCA finished by error termination in SCF Hessian
Calling Command: mpirun -np 32  /public/home/bin/orca_5_0_3_linux_x86-64_shared_openmpi411/orca_scfhess_mpi 1.scfhess.inp 1
[file orca_tools/qcmsg.cpp, line 465]:
  .... aborting the run

作者
Author: hanshan 时间: 2023-4-13 17:51

sobereva 发表于 2023-4-13 17:30
值得一提的是，6-31G(2df,p)对于DFT来说本身不是什么好的选择。把基组从2-zeta增大到3-zeta比加昂贵的f极化 ...

qm9数据库论文是02 April 2014投稿的，可能也有时代的原因
您是建议用def2-TZVP，甚至def2-QZVP吧

作者
Author: sobereva 时间: 2023-4-13 18:36

hanshan 发表于 2023-4-13 17:51
qm9数据库论文是02 April 2014投稿的，可能也有时代的原因
您是建议用def2-TZVP，甚至def2-QZVP吧

那个基组在14年也是明显不合时宜的。def2早就有了
是。def2-TZVP算大批量分子偏贵，比6-31G(2df,p)更便宜的def-TZVP都比6-31G(2df,p)强

作者
Author: wzkchem5 时间: 2023-4-13 18:41

hanshan 发表于 2023-4-13 10:43
原来您是orca开发团队成员，谢谢您的答复！
其实我还有一个疑问，以下报错可能是什么情况呢？

maxcore必须大于计算所需的每核最小内存（输出文件给的就是每核需要的内存），但是小于当前可用物理内存除以所用核数。可能你设的maxcore满足了前一个条件，但没有满足后一个条件

作者
Author: hanshan 时间: 2023-4-13 21:13

wzkchem5 发表于 2023-4-13 18:41
maxcore必须大于计算所需的每核最小内存（输出文件给的就是每核需要的内存），但是小于当前可用物理内存 ...

确实，128G内存的服务器，32核心并行，maxcore设置成3000就报错了
跳到2500，暂时还没有问题

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