用orca算SOC矩阵元，其对几何结构高度敏感？

皖酒

我是打算用orca算S1-T1之间的ISC，看了一些帖子以及orca手册，需要先优化S1和T1态、生成.hess文件，并结合SOC矩阵元计算。我之前已经用Gaussian优化过了S1和T1（然后用这个S1优化后的结构，用orca算了一次SOC,与文献相对符合），然后将他们作为初猜，用ORCA重新优化并计算频率，取优化后的S1结构，用orca再次做SOC计算，结果相差非常大，望老师们指点。
SOC_1输入：
! M06 6-31G* miniprint tightSCF pal4 SMD(Acetonitrile)
%tddft
nroots 5
dosoc true
tda false
printlevel 3
end

SOC_1输出如下图：


S1的优化输入：
! 6-31G* tightopt SMD(Acetonitrile) OPT FREQ miniprint nopop RIJCOSX
%method
functional M06
end
%maxcore 2000
%pal nprocs 20 end
%tddft
nroots 5
tda false
IROOT 1
end
SOC_2输入：
! M06 6-31G* miniprint tightSCF pal4 SMD(Acetonitrile)
%tddft
nroots 5
dosoc true
tda false
printlevel 3
end
SOC_2输出如下图：

wal

由于结构变化，透热态的能量顺序变了，而soc很大程度上是跟着透热态走的。

sobereva

帖子里的描述和标题十分不符，电子结构差异和几何结构差异是两码事，先弄清楚差异的源头主要在哪。算SOC之前先明确ORCA对Gaussian优化完的激发态进一步优化后，和Gaussian原本给出的结构有什么差异，结构没什么差异再说别的；结构若都有显著差异，先从影响结构优化的因素上弄清楚是怎么回事（如RI的影响、收敛限的影响、俩程序默认的D基函数的球谐型和笛卡尔型的差异等等）

皖酒

wal 发表于 2026-7-2 00:54
由于结构变化，透热态的能量顺序变了，而soc很大程度上是跟着透热态走的。

好的，谢谢老师，我去看看

皖酒

sobereva 发表于 2026-7-2 02:15
帖子里的描述和标题十分不符，电子结构差异和几何结构差异是两码事，先弄清楚差异的源头主要在哪。算SOC之 ...

已修改，谢谢老师