orca5.0.4关键词def2-QZVPP(-g,-f)不存在

mzf · 发表于 Post on 2023-7-31 12:46:55

大家好，
我最近阅读了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)。感谢各位解惑！

Deepast · 发表于 Post on 2023-7-31 13:07:57

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）

你实在想用可以自定义

mzf · 发表于 Post on 2023-7-31 13:13:30

Deepast 发表于 2023-7-31 13:07
你实在想用可以自定义

谢谢回复！请问就是把def2-QZVPP 自定义，然后把里面关于g和f的部分删除吗？我是小白，不太了解

感谢

Deepast · 发表于 Post on 2023-7-31 14:39:49

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

不是全部删
去基组数据库找找
下个ORCA4用orca_exportbasis输出基组可能也行

zjxitcc · 发表于 Post on 2023-7-31 20:28:17

本帖最后由 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文件，以水分子为例，

%chk=h2o.chk
%mem=1GB
%nprocshared=1
#p RHF gen nosymm int=nobasistransform guess(only,save)
title
0 1
O -0.65151515 0.25757575 0.0
H 0.30848485 0.25757575 0.0
H -0.97196973 1.16251158 0.0
H 0
S 4 1.00
190.6916900 0.70815167D-03
28.6055320 0.54678827D-02
6.5095943 0.27966605D-01
1.8412455 0.10764538
S 1 1.00
0.59853725 1.0000000
S 1 1.00
0.21397624 1.0000000
S 1 1.00
0.80316286D-01 1.0000000
P 1 1.00
2.29200000 1.0000000
P 1 1.00
0.83800000 1.0000000
P 1 1.00
0.29200000 1.0000000
D 1 1.00
2.06200000 1.0000000
D 1 1.00
0.66200000 1.0000000
F 1 1.00
1.39700000 1.0000000
****
O 0
S 8 1.00
116506.4690800 0.40383857939D-04
17504.3497240 0.31255139004D-03
3993.4513230 0.16341473495D-02
1133.0063186 0.68283224757D-02
369.99569594 0.24124410221D-01
133.62074349 0.72730206154D-01
52.035643649 0.17934429892
21.461939313 0.33059588895
S 2 1.00
89.835051252 0.96468652996D-01
26.428010844 0.94117481120
S 1 1.00
9.2822824649 1.0000000
S 1 1.00
4.0947728533 1.0000000
S 1 1.00
1.3255349078 1.0000000
S 1 1.00
0.51877230787 1.0000000
S 1 1.00
0.19772676454 1.0000000
P 5 1.00
191.15255810 0.25115697705D-02
45.233356739 0.20039240864D-01
14.353465922 0.93609064762D-01
5.2422371832 0.30618127124
2.0792418599 0.67810501439
P 1 1.00
0.84282371424 1.0000000
P 1 1.00
0.33617694891 1.0000000
P 1 1.00
0.12863997974 1.0000000
D 1 1.00
3.77500000 1.0000000
D 1 1.00
1.30000000 1.0000000
D 1 1.00
0.44400000 1.0000000
F 1 1.00
2.66600000 1.0000000
F 1 1.00
0.85900000 1.0000000
G 1 1.00
1.84600000 1.0000000
****

复制代码

里的数据表示(-f,-g)，也就是你要删除的部分（由于论坛代码没法显示粗体，我只好用表示）。提交Gaussian计算，1s完成，因为only表示只产生初猜，不做SCF计算。然后执行

fch2mkl h2o.chk

复制代码

会生成h2o_o.inp和h2o_o.mkl文件，此时这两个文件里都是def2-QZVPP(-f,-g)基组，还有电荷、自旋多重度、坐标这些基本信息都写好了，在h2o_o.inp文件里添加关键词做你要想要的计算就行。fch2mkl小程序可以免费下载。

mzf · 发表于 Post on 2023-7-31 22:34:15

Deepast 发表于 2023-7-31 14:39
不是全部删
去基组数据库找找
下个ORCA4用orca_exportbasis输出基组可能也行

非常感谢！我去尝试一下

mzf · 发表于 Post on 2023-7-31 22:35:11

zjxitcc 发表于 2023-7-31 20:28
博文里有说“def2-QZVPP(-g,-f)代表从重原子和氢原子上分别去掉最高角动量函数g和f”。实际操作可以有很多 ...

感谢您如此具体的解答！我会去尝试一下。另想问一下因为关键词是-g，-f，是否F相关的数据也应该一并删除呢

zjxitcc · 发表于 Post on 2023-8-1 09:13:32

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

建议你再阅读几遍我的回答，我写的很详细，你问的我都回答了，你只是重复问了个问题。

mzf · 发表于 Post on 2023-8-1 12:45:20

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

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

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

[ORCA] orca5.0.4关键词def2-QZVPP(-g,-f)不存在

评分 Rate

浏览过的版块