计算化学公社

标题: 求助：CASSCF计算的SOC校正吸收光谱中振子强度与跃迁电偶极矩的关系 [打印本页]

作者
Author: warfarin 时间: 2023-6-28 18:55
标题: 求助：CASSCF计算的SOC校正吸收光谱中振子强度与跃迁电偶极矩的关系
本帖最后由 warfarin 于 2023-6-28 18:58 编辑

各位老师好，我在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之后振子强度的计算需要加入修正项吗，还是我的计算方式有问题？
(, 下载次数 Times of downloads: 23)
2）这个计算得到的振子强度都太小了，所以我希望振子强度的小数保留位数能够多一些，请问ORCA是否有命令能够调整小数保留位数？

作者
Author: sobereva 时间: 2023-6-29 07:03
别乱用“常数”这个词。本来振子强度就不是什么常数，而是一个依赖于各种因素的标量值。都不知道你算出来了什么东西

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

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