求助，ORCA 自定义Cu元素基组报错(已解决)

love_yy

输入信息是按照Multiwfn 生成的，如下：





报错信息提示如下：

KSeGaSn

赝势 newECP "SDD" end
赝势基组 newGTO的话要写完整的定义 SDD没用

love_yy

KSeGaSn 发表于 2023-12-28 17:33
赝势 newECP "SDD" end
赝势基组 newGTO的话要写完整的定义 SDD没用

好嘞谢谢，我参考了这篇帖子（http://bbs.keinsci.com/thread-10553-1-1.html）和有些资料，现在写的是这种输入方式，还是报错“Error: Pure ECPs not allowed in GTO definition” （我的版本是5.0.4）

%basis
NewGTO Cu "SDD" end
NewECP  "SDD" end
NewAuxjGTO Cu "autoaux" end
end

您说的写完整定义SDD，没太明白啥意思，这个关键词手册里说可以写呀

KSeGaSn

love_yy 发表于 2023-12-28 20:22
好嘞谢谢，我参考了这篇帖子（http://bbs.keinsci.com/thread-10553-1-1.html）和有些资料，现在写的是这 ...

因为ORCA只有内置的SDD赝势而没有对应的（价）电子基组 显然你这相当于只有内层电子的信息没有价电子/价轨道的信息啊 所以需要在newGTO的部分写出电子基组的信息 手册说的很明白“the NewECP keyword……only affect the ECP and not the valence basis set”。所以正确的输入方法是

%basis
NewGTO Cu
S   3
1        27.6963200              0.2311320
2        13.5053500             -0.6568110
3         8.8153550             -0.5458750
S   1
1         2.3808050              1.0000000
S   1
1         0.9526160              1.0000000
S   1
1         0.1126620              1.0000000
S   1
1         0.0404860              1.0000000
S   1
1         0.0100000              1.0000000
P   2
1        93.5043270              0.0228290
2        16.2854640             -1.0095130
P   2
1         5.9942360              0.2464500
2         2.5368750              0.7920240
P   1
1         0.8979340              1.0000000
P   1
1         0.1317290              1.0000000
P   1
1         0.0308780              1.0000000
D   4
1        41.2250060              0.0446940
2        12.3432500              0.2121060
3         4.2019200              0.4534230
4         1.3798250              0.5334650
D   1
1         0.3834530              1.0000000
D   1
1         0.1000000              1.0000000
end
NewAuxJGTO Cu "autoaux" end
NewECP Cu "SDD" end
end

KSeGaSn

love_yy 发表于 2023-12-28 20:22
好嘞谢谢，我参考了这篇帖子（http://bbs.keinsci.com/thread-10553-1-1.html）和有些资料，现在写的是这 ...

这段newGTO的定义是从basisexchange上下载下来的，Stuttgart RSC 1997，其中描述就说明了这包含两个部分ECP + Valence basis。

love_yy

KSeGaSn 发表于 2023-12-29 09:04
这段newGTO的定义是从basisexchange上下载下来的，Stuttgart RSC 1997，其中描述就说明了这包含两个部分E ...

好，谢谢

zjxitcc

顺便提醒一下，你这个体系有明显、不少负电荷，建议给涉及明显负电荷的元素使用ma-TZVPP或def2TZVPPD，给其他元素使用def2-TZVPP。讲道理没有必要给Cu单独用SDD，毕竟def2系列基组对Cu有定义（不过是全电子基组，而非像SDD那样带赝势）。

love_yy

zjxitcc 发表于 2023-12-29 11:37
顺便提醒一下，你这个体系有明显、不少负电荷，建议给涉及明显负电荷的元素使用ma-TZVPP或def2TZVPPD，给其 ...

好嘞，谢谢老师