计算化学公社

标题: 求助，ORCA 自定义Cu元素基组报错(已解决) [打印本页]

作者
Author: love_yy 时间: 2023-12-28 15:09
标题: 求助，ORCA 自定义Cu元素基组报错(已解决)
本帖最后由 love_yy 于 2024-1-3 15:22 编辑

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

(, 下载次数 Times of downloads: 27)

报错信息提示如下：

(, 下载次数 Times of downloads: 24)

作者
Author: KSeGaSn 时间: 2023-12-28 17:33
赝势 newECP "SDD" end
赝势基组 newGTO的话要写完整的定义 SDD没用

作者
Author: love_yy 时间: 2023-12-28 20:22

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，没太明白啥意思，这个关键词手册里说可以写呀

作者
Author: KSeGaSn 时间: 2023-12-29 09:01

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

作者
Author: KSeGaSn 时间: 2023-12-29 09:04

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

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

作者
Author: love_yy 时间: 2023-12-29 09:27

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

好

，谢谢

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

作者
Author: love_yy 时间: 2023-12-30 17:33

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

好嘞，谢谢老师

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