求助关于旋轨耦合计算报错Spin-orbit coupling not yet implemented with ECPs

az￥ · 发表于 Post on 2025-5-12 09:46:20

大家好！如果计算的结构里带有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
这样的参数，结果是正常结束的，我想问一下，这样做可以吗？期待并感谢大家的回复！

Stardust0831 · 发表于 Post on 2025-5-12 10:10:01

本帖最后由 Stardust0831 于 2025-5-12 10:17 编辑

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

一个类似的example：

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

az￥ · 发表于 Post on 2025-5-13 09:30:50

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，这样计算合理吗？

Stardust0831 · 发表于 Post on 2025-5-13 09:45:24

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

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

az￥ · 发表于 Post on 2025-5-15 09:35:00

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

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

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

[ORCA] 求助关于旋轨耦合计算报错Spin-orbit coupling not yet implemented with ECPs

浏览过的版块