计算化学公社

标题: orca5.0.4关键词def2-QZVPP(-g,-f)不存在 [打印本页]
作者
Author: mzf    时间: 2023-7-31 12:46
标题: orca5.0.4关键词def2-QZVPP(-g,-f)不存在
大家好，
我最近阅读了sob的帖子大体系弱相互作用计算的解决之道（http://sobereva.com/214），其中4.8节提到了对于几十个原子的体系建议使用def2-QZVPP(-g,-f)。我尝试直接复制原文关键词（RI-PWPB95 D3 def2-QZVPP(-g,-f) def2-QZVPP/J def2-QZVPP/C RIJCOSX grid4 tightscf）提示报错，但是如果把def2-QZVPP(-g,-f)改为def2-QZVPP就不会报错了。由于sob提到“def2-QZVPP(-g,-f)提升到def2-QZVPP给总耗时带来了近一半的增加”，我还是想尝试一下这个基组。因此这里想问一下在orca中如何实现def2-QZVPP(-g,-f)？我使用的版本是orca5.0.4，尝试搜索了官方例子里面的文件也都是def2-QZVPP而不是def2-QZVPP(-g,-f)。感谢各位解惑！
作者
Author: Deepast    时间: 2023-7-31 13:07
2018-Apr-30补充：注意，从ORCA4.0开始，辅助基组定义方式相较于ORCA3.0有了极大的变化。上文的def2-SVP/J应改为def2/J，而def2-QZVPP/JK应改为def2/JK。下文提到的def2-QZVPP(-g,-f)基组在ORCA4.0中也没法直接用相应关键词使用了。——大体系弱相互作用计算的解决之道（http://sobereva.com/214）

你实在想用可以自定义
作者
Author: mzf    时间: 2023-7-31 13:13
Deepast 发表于 2023-7-31 13:07
你实在想用可以自定义

谢谢回复！请问就是把def2-QZVPP 自定义，然后把里面关于g和f的部分删除吗？我是小白，不太了解感谢
作者
Author: Deepast    时间: 2023-7-31 14:39
mzf 发表于 2023-7-31 13:13
谢谢回复！请问就是把def2-QZVPP 自定义，然后把里面关于g和f的部分删除吗？我是小白，不太了解:handshak ...

不是全部删
去基组数据库找找
下个ORCA4用orca_exportbasis输出基组可能也行
作者
Author: zjxitcc    时间: 2023-7-31 20:28
本帖最后由 zjxitcc 于 2023-7-31 20:31 编辑

博文里有说“def2-QZVPP(-g,-f)代表从重原子和氢原子上分别去掉最高角动量函数g和f”。实际操作可以有很多种，比如
（1）楼上说的借用旧版ORCA功能。
（2）到BSE基组网站下载ORCA格式的def2-QZVPP基组数据下来自己删减，在ORCA输入文件里用自定义基组。
（3）到BSE基组网站下载Gaussian格式的def2-QZVPP基组数据，写个gjf文件，以水分子为例，
  1. %chk=h2o.chk
  2. %mem=1GB
  3. %nprocshared=1
  4. #p RHF gen nosymm int=nobasistransform guess(only,save)

  6. title

  8. 0 1
  9. O   -0.65151515    0.25757575    0.0
  10. H    0.30848485    0.25757575    0.0
  11. H   -0.97196973    1.16251158    0.0

  13. H     0
  14. S    4   1.00
  15.     190.6916900              0.70815167D-03
  16.      28.6055320              0.54678827D-02
  17.       6.5095943              0.27966605D-01
  18.       1.8412455              0.10764538
  19. S    1   1.00
  20.       0.59853725             1.0000000
  21. S    1   1.00
  22.       0.21397624             1.0000000
  23. S    1   1.00
  24.       0.80316286D-01         1.0000000
  25. P    1   1.00
  26.       2.29200000             1.0000000
  27. P    1   1.00
  28.       0.83800000             1.0000000
  29. P    1   1.00
  30.       0.29200000             1.0000000
  31. D    1   1.00
  32.       2.06200000             1.0000000
  33. D    1   1.00
  34.       0.66200000             1.0000000
  35. <b>F    1   1.00
  36.       1.39700000             1.0000000</b>
  37. ****
  38. O     0
  39. S    8   1.00
  40. 116506.4690800              0.40383857939D-04
  41.   17504.3497240              0.31255139004D-03
  42.    3993.4513230              0.16341473495D-02
  43.    1133.0063186              0.68283224757D-02
  44.     369.99569594             0.24124410221D-01
  45.     133.62074349             0.72730206154D-01
  46.      52.035643649            0.17934429892
  47.      21.461939313            0.33059588895
  48. S    2   1.00
  49.      89.835051252            0.96468652996D-01
  50.      26.428010844            0.94117481120
  51. S    1   1.00
  52.       9.2822824649           1.0000000
  53. S    1   1.00
  54.       4.0947728533           1.0000000
  55. S    1   1.00
  56.       1.3255349078           1.0000000
  57. S    1   1.00
  58.       0.51877230787          1.0000000
  59. S    1   1.00
  60.       0.19772676454          1.0000000
  61. P    5   1.00
  62.     191.15255810             0.25115697705D-02
  63.      45.233356739            0.20039240864D-01
  64.      14.353465922            0.93609064762D-01
  65.       5.2422371832           0.30618127124
  66.       2.0792418599           0.67810501439
  67. P    1   1.00
  68.       0.84282371424          1.0000000
  69. P    1   1.00
  70.       0.33617694891          1.0000000
  71. P    1   1.00
  72.       0.12863997974          1.0000000
  73. D    1   1.00
  74.       3.77500000             1.0000000
  75. D    1   1.00
  76.       1.30000000             1.0000000
  77. D    1   1.00
  78.       0.44400000             1.0000000
  79. F    1   1.00
  80.       2.66600000             1.0000000
  81. F    1   1.00
  82.       0.85900000             1.0000000
  83. <b>G    1   1.00
  84.       1.84600000             1.0000000</b>
  85. ****

复制代码
<b></b>里的数据表示(-f,-g)，也就是你要删除的部分（由于论坛代码没法显示粗体，我只好用<b></b>表示）。提交Gaussian计算，1s完成，因为only表示只产生初猜，不做SCF计算。然后执行
  1. fch2mkl h2o.chk
复制代码
会生成h2o_o.inp和h2o_o.mkl文件，此时这两个文件里都是def2-QZVPP(-f,-g)基组，还有电荷、自旋多重度、坐标这些基本信息都写好了，在h2o_o.inp文件里添加关键词做你要想要的计算就行。fch2mkl小程序可以免费下载。
作者
Author: mzf    时间: 2023-7-31 22:34
Deepast 发表于 2023-7-31 14:39
不是全部删
去基组数据库找找
下个ORCA4用orca_exportbasis输出基组可能也行

非常感谢！我去尝试一下
作者
Author: mzf    时间: 2023-7-31 22:35
zjxitcc 发表于 2023-7-31 20:28
博文里有说“def2-QZVPP(-g,-f)代表从重原子和氢原子上分别去掉最高角动量函数g和f”。实际操作可以有很多 ...

感谢您如此具体的解答！我会去尝试一下。另想问一下因为关键词是-g，-f，是否F相关的数据也应该一并删除呢
作者
Author: zjxitcc    时间: 2023-8-1 09:13
mzf 发表于 2023-7-31 22:35
感谢您如此具体的解答！我会去尝试一下。另想问一下因为关键词是-g，-f，是否F相关的数据也应该一并删除 ...

建议你再阅读几遍我的回答，我写的很详细，你问的我都回答了，你只是重复问了个问题。
作者
Author: mzf    时间: 2023-8-1 12:45
zjxitcc 发表于 2023-8-1 09:13
建议你再阅读几遍我的回答，我写的很详细，你问的我都回答了，你只是重复问了个问题。

你好，之前沒有看到段落末尾的<\b>，看到了之後就理解了，并且也成功運行程序了，再次感謝！



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3