求助：CASSCF计算的SOC校正吸收光谱中振子强度与跃迁电偶极矩的关系

warfarin

各位老师好，我在ORCA中通过CASSCF计算得到了Er离子团簇考虑SOC的吸收光谱，输入文件如下：
!DKH2 DKH-def2-SVP Moread tightSCF miniprint nopop
%moinp "HF.gbw"
%maxcore 8000
%pal nprocs  48 end
%basis
NEWGTO Er "SARC2-DKH-QZVP" end
end
%SCF
MaxIter 300
DirectResetFreq=1
CNVShift true
LShift 0.5
ShiftErr 0.1
end
%casscf
nel 11
norb 7
mult 4,2
nroots 35,112

rel
   dosoc true
   NInitstates 16
end
end

* xyz   0   4
[geometry]
*
我对其中振子强度的计算方式有一些疑问，想请教一下各位老师。
1）我在sob老师的博文中曾经看到过振子强度计算公式f=2/3*(Ej-Ei)*|<i|-r|j>|^2。按照此公式，无论用什么单位，f/(Ej-Ei)/|<i|-r|j>|^2都应该会是一个常数，但我使用输出文件中不同跃迁给出的振子强度、能级差和跃迁电偶极矩计算得到的数值不是常数，请问这是因为考虑SOC之后振子强度的计算需要加入修正项吗，还是我的计算方式有问题？

2）这个计算得到的振子强度都太小了，所以我希望振子强度的小数保留位数能够多一些，请问ORCA是否有命令能够调整小数保留位数？

sobereva

别乱用“常数”这个词。本来振子强度就不是什么常数，而是一个依赖于各种因素的标量值。都不知道你算出来了什么东西

振子强度非常小就当成0就完了，保留更多没有任何意义