ORCA计算SOC：体系偶极矩较大是否会影响SOC矩阵元结果？

大太阳

各位前辈好， 我是计算初学者，在用 ORCA 做自旋-轨道耦合（SOC）计算时遇到了一些疑惑，想请教大家。

一、计算背景
体系：含I的有机金属化合物，原子数超过300
优化方法：使用高斯，B3LYP-D3，C、H、N、O 用 6-31G(d)，I 用 LANL2DZ
SOC计算：使用 5.0.3版本ORCA，B3LYP/G D3 6-31G(d)  I 用LANL2DZ二、具体问题
体系是高斯优化后，在ORCA计算SOC，计算无报错正常结束，但偶极矩为 300 Debye。
疑问是：
1.偶极矩较大是否会对 SOC 计算结果产生影响？计算结果正常结束，对于偶极矩过大问题是否可以忽略？

在使用6.1.1版本ORCA计算时，加入了NewECP后，计算报错找不到ECP。
输入文件格式为：
! B3LYP/G D3 6-31G(d) SARC/J RIJCOSX tightSCF miniprint
%basis
NewGTO I "LANL2DZ" end
NewECP I "LANL2DZ" end
end
%tddft
nroots 3
TDA true
dosoc true
printlevel 5
end
* xyz 0 1
疑问：
1.这段输入文件的写法是否有错误？


望各位大佬老师们解惑！

sobereva

跟偶极矩大小毫无关系

当前完全是胡算，认真把下文看了，正确使用相对论哈密顿、正确选择全电子基组
使用ORCA在TDDFT下计算旋轨耦合矩阵元和绘制旋轨耦合校正的UV-Vis光谱
http://sobereva.com/462（http://bbs.keinsci.com/thread-12114-1-1.html）

大太阳

谢谢大佬！

大太阳

老师，我仔细看了文章，看了输入文件，我在优化以及算ΔE_s1-t1使用的是B3LYP-D3 6-31G(d)，I用LANL2DZ，后续为匹配 SOC 自旋轨道耦合计算、评估分子系间窜越（ISC）性能，ORCA 计算也沿用了完全一致的泛函与基组。此前 TDA 设为 true 是程序报错后，参考相关技术帖调整的参数。请问是对于I、Br等重原子是否不适合搭配 LANL2DZ，计算 SOC 时必须更换配套自旋轨道哈密顿对应的赝势基组？？麻烦老师指点，谢谢老师！🙏🌹

wal

大太阳 发表于 2026-8-1 21:35
老师，我仔细看了文章，看了输入文件，我在优化以及算ΔE_s1-t1使用的是B3LYP-D3 6-31G(d)，I用LANL2DZ，后 ...

旋轨耦合计算无法使用赝势，必须使用全电子基组计算
你的输入文件根本没用赝势，一个lanl2dz空壳摆在那，所以才能正常（个锤子）计算旋轨耦合，但也导致偶极矩见了鬼一样高

大太阳

wal 发表于 2026-8-2 10:47
旋轨耦合计算无法使用赝势，必须使用全电子基组计算
你的输入文件根本没用赝势，一个lanl2dz空壳摆在那 ...

好的好的，感谢大佬！❀