如何给特定原子添加ma-方式的弥散函数

cokie · 发表于 Post on 2024-1-22 15:18:40

各位老师好，我根据博文《给def2以ma-方式加弥散函数的Gaussian格式的基组定义文件（含所有def2支持的元素）》（http://sobereva.com/509）来调用“ma-def2-TZVP”。不过我想做的是给特定几个原子加弥散（水解的羧基），其他原子用def2-TZVP。遇到两个问题。我的输入文件如下：

#p PBE1PBE/gen em=gd3bj TD(singlets,NSTATES=15,ROOT=1) IOp(9/40=4)

0 1
……
空行
1 2 5 6 29 30 0
@C:\orcajob\HP\C3T-PBE0-diffuse-def2tzvp\ma-TZVP.txt
****
3-4 7-28 31-110 0
def2tzvp
****
空行
空行

我才用了sob老师提供的“file”里的.txt文件，但是提交任务后报错显示：定义的基组错误，从“-H”开始这一行就无法读取。之后我选择钯.txt里的“-C还有-O”的负号去了（加弥散的原子），可以正确读取并计算，但程序显示：

****
Warning:  center 9 has no basis functions!
Warning:  center  10 has no basis functions!
Warning:  center  23 has no basis functions!
Warning:  center  24 has no basis functions!
Warning:  center  35 has no basis functions!
Warning:  center  36 has no basis functions!
Warning:  center  37 has no basis functions!
Warning:  center  38 has no basis functions!
Warning:  center  42 has no basis functions!
Warning:  center  43 has no basis functions!
Warning:  center  45 has no basis functions!
Warning:  center  46 has no basis functions!
Warning:  center  48 has no basis functions!
Warning:  center  49 has no basis functions!
Warning:  center  51 has no basis functions!
Warning:  center  52 has no basis functions!
Warning:  center  54 has no basis functions!
Warning:  center  55 has no basis functions!
Warning:  center  57 has no basis functions!
Warning:  center  58 has no basis functions!
Warning:  center  60 has no basis functions!
Warning:  center  61 has no basis functions!
Warning:  center  63 has no basis functions!
Warning:  center  64 has no basis functions!

…………………………………………………………

是我不能用这种方式来给局部加弥散函数吗？正确的方式应该是什么呢？恳请老师们赐教！谢谢！

wzkchem5 · 发表于 Post on 2024-1-22 17:35:37

报错的这些center对应的原子分别是什么元素？

cokie · 发表于 Post on 2024-1-22 17:47:07

本帖最后由 cokie 于 2024-1-22 17:48 编辑

wzkchem5 发表于 2024-1-22 17:35
报错的这些center对应的原子分别是什么元素？

C，H，O，N元素都有。谢谢老师，我目前已经能正确运行了（虽然不知道我解决的方案对不对）
我把自定义基组的地方改写作了这样：

……………………
空行
3-4 7-28 31-110 0
def2tzvp
****
1 2 5 6 0
S 6 1.00
27032.3826310 0.21726302465E-03
4052.3871392 0.16838662199E-02
922.32722710 0.87395616265E-02
261.24070989 0.35239968808E-01
85.354641351 0.11153519115
31.035035245 0.25588953961
S 2 1.00
12.260860728 0.39768730901
4.9987076005 0.24627849430
S 1 1.00
1.1703108158 1.0000000
S 1 1.00
0.46474740994 1.0000000
S 1 1.00
0.18504536357 1.0000000
P 4 1.00
63.274954801 0.60685103418E-02
14.627049379 0.41912575824E-01
4.4501223456 0.16153841088
1.5275799647 0.35706951311
P 1 1.00
0.52935117943 1.0000000
P 1 1.00
0.17478421270 1.0000000
D 1 1.00
2.31400000 1.0000000
D 1 1.00
0.64500000 1.0000000
F 1 1.00
1.42800000 1.0000000
S 1 1.00
0.616817878567E-01 1.000
P 1 1.00
0.582614042333E-01 1.000
****
29 30 0
S 6 1.00
13575.3496820 0.22245814352E-03
2035.2333680 0.17232738252E-02
463.22562359 0.89255715314E-02
131.20019598 0.35727984502E-01
42.853015891 0.11076259931
15.584185766 0.24295627626
S 2 1.00
6.2067138508 0.41440263448
2.5764896527 0.23744968655
S 1 1.00
0.57696339419 1.0000000
S 1 1.00
0.22972831358 1.0000000
S 1 1.00
0.95164440028E-01 1.0000000
P 4 1.00
34.697232244 0.53333657805E-02
7.9582622826 0.35864109092E-01
2.3780826883 0.14215873329
0.81433208183 0.34270471845
P 1 1.00
0.28887547253 1.0000000
P 1 1.00
0.10056823671 1.0000000
D 1 1.00
1.09700000 1.0000000
D 1 1.00
0.31800000 1.0000000
F 1 1.00
0.76100000 1.0000000
S 1 1.00
0.317214800093E-01 1.000
P 1 1.00
0.335227455700E-01 1.000
****
空行
空行

复制代码

就可以了，我猜测可能是因为sob老师的.txt文件是针对元素定义的，但是我加弥散的六个原子是C和O元素，和sob老师写的.txt里的基组定义冲突了。所以我就分别把里面对O和C的定义copy出来分开加的，现在就能正常跑了。（应该这种处理方式是对的吧）

sobereva · 发表于 Post on 2024-1-22 23:00:16

cokie 发表于 2024-1-22 17:47
C，H，O，N元素都有。谢谢老师，我目前已经能正确运行了（虽然不知道我解决的方案对不对）
我把自定义基 ...

对

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

[新手求助] 如何给特定原子添加ma-方式的弥散函数