计算化学公社

标题: 求助关于旋轨耦合计算报错Spin-orbit coupling not yet implemented with ECPs [打印本页]

作者
Author: az￥ 时间: 2025-5-12 09:46
标题: 求助关于旋轨耦合计算报错Spin-orbit coupling not yet implemented with ECPs
大家好！如果计算的结构里带有I原子，我结构优化的时候没有用赝势基组来优化，整个结构用的都是pbe0/def2-TZVP，但是当我计算旋轨耦合的时候，用
B3LYP/def2-TZVP时，会有“[file orca_main/mainchk.cpp, line 2207]: ERROR (CIS/TDDFT): Spin-orbit coupling not yet implemented with ECPs”这样的报错出现，我看了sob老师发布的相关帖子，就在输入文件里加上了
%basis
NewGTO I "SARC-DKH-TZVP" end
end
这样的参数，结果是正常结束的，我想问一下，这样做可以吗？期待并感谢大家的回复！

作者
Author: Stardust0831 时间: 2025-5-12 10:10
本帖最后由 Stardust0831 于 2025-5-12 10:17 编辑

SARC-DKH-TZVP是给相对论计算用的。你要用这个基组的话把DKH2开了，同时把def2-TZVP换成DKH-def2-TZVP

一个类似的example：

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

或者参考《ORCA计算含有重原子和不含重原子的有机分子的S1-T1旋轨耦合》

作者
Author: az￥ 时间: 2025-5-13 09:30

Stardust0831 发表于 2025-5-12 10:10
SARC-DKH-TZVP是给相对论计算用的。你要用这个基组的话把DKH2开了，同时把def2-TZVP换成DKH-def2-TZVP

...

谢谢您的回复！DKH2我打开了，基组也换成DKH-def2-TZVP了，我想问一下，像我这样，优化I原子用的pbe0/def2-TZVP和计算旋轨耦合I原子用SARC-DKH-TZVP，这样计算合理吗？

作者
Author: Stardust0831 时间: 2025-5-13 09:45

az￥发表于 2025-5-13 09:30
谢谢您的回复！DKH2我打开了，基组也换成DKH-def2-TZVP了，我想问一下，像我这样，优化I原子用的pbe0/def ...

换级别是有合理理由的。带着DKH2算结构优化也很难算得动，单点又有非开DKH2不可的理由。我觉得没问题。

作者
Author: az￥ 时间: 2025-5-15 09:35

Stardust0831 发表于 2025-5-13 09:45
换级别是有合理理由的。带着DKH2算结构优化也很难算得动，单点又有非开DKH2不可的理由。我觉得没问题。

好，明白了，非常感谢您的解答

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