FCClasses3 计算RISC速率为负数

lennylyp

想用FCClasses3算一下TADF分子的逆向系间窜越速率，结果得到的是负数……除了旋轨耦合和ΔEst之外其他参数都一样，有些分子是正的，有些是负的，求问各位前辈是哪里出了问题……

输入文件参考了这个网页：https://www.homechemer.com/computchem/9
构型优化和频率分析用Gaussian在B3lyp/6-31g(d)下做的，SOC用ORCA算的，ΔEst用SCS-CC2算的

$$$
PROPERTY     =   NR0  ; OPA/EMI/ECD/CPL/RR/TPA/MCD/IC
MODEL        =   AH   ; Adiabatic Hessian; AS/ASF/AH/VG/VGF/VH
DE           =   0.181107455 ； Energy difference in ev
NR0_COUPL    =   6.47E-07 ; Spin-orbit coupling integral in a.u.
TEMP         =   273.15 ; Temperature in K
BROADFUN     =   GAU  ; GAU/LOR/VOI
HWHM         =   0.06 ; Broadening width in eV
METHOD       =   TD   ; Time-dependent; TI/TD
ROT          =   1    ; Rotate state 1
;VIBRATIONAL ANALYSIS
NORMALMODES  =   COMPUTE   ; COMPUTE/READ/IMPLICIT
COORDS       =   CARTESIAN ; CARTESIAN/INTERNAL
;INPUT DATA FILES
STATE1_FILE  =   m-s1.fcc
STATE2_FILE  =   m-t1.fcc

lennylyp

额外发现温度设置得低一点，可以避免出现负数。但这又是什么原因导致的呢？

Novice

我也用FCClasses3算过，不知道是因为不熟悉导致创建的输入文件不对还是什么，得到的结果也不合理（数量级严重偏差），但是运行示例文件却没有问题。我当时也有把温度改为300k，不知道是不是有bug?

15689900879@163

试试不同的展宽函数和系数，设置rotate=-1，internal内坐标。

ABQTrap

你好，请问计算ISC时，输入文件中的旋轨耦合数值是哪个态的SOC？初态？末态？还是MECP处的？谢谢！

Eclair

ABQTrap 发表于 2024-2-26 20:55
你好，请问计算ISC时，输入文件中的旋轨耦合数值是哪个态的SOC？初态？末态？还是MECP处的？谢谢！

没有哪个态的soc这一个说法，soc是两个态之间的相互作用，

ABQTrap

Eclair 发表于 2024-2-29 15:03
没有哪个态的soc这一个说法，soc是两个态之间的相互作用，

比如说S1和T1间的ISC，soc是在S1优化的结构基础上进行SOC计算，还是优化后的T1结构，还是优化的MECP处的结构进行计算SOC？谢谢！

wzkchem5

ABQTrap 发表于 2024-3-2 10:41
比如说S1和T1间的ISC，soc是在S1优化的结构基础上进行SOC计算，还是优化后的T1结构，还是优化的MECP处的 ...

这个问题我在这个论坛上回答了无数次，随便一搜就知道了
http://bbs.keinsci.com/forum.php ... 6&fromuid=21811
http://bbs.keinsci.com/forum.php ... 4&fromuid=21811

T9的梦

15689900879@163 发表于 2024-2-19 17:07
试试不同的展宽函数和系数，设置rotate=-1，internal内坐标。

大佬您好，我试着调了各种展宽函数和系数，发现VOI函数和HWHM=40 eV时能基本得到想要的速率，但HWHM默认值是0.01eV，偏差了太多。请问这样合理吗？
（我不清楚HWHM是否是实验/计算光谱的半峰宽，比如1240/40 =31 nm，而默认值0.01算出来的半峰宽有12400 nm）感觉非常迷茫

wzkchem5

T9的梦 发表于 2024-4-2 12:53
大佬您好，我试着调了各种展宽函数和系数，发现VOI函数和HWHM=40 eV时能基本得到想要的速率，但HWHM默认 ...

HWHM指的是你把谱的横坐标换算成eV以后，量出来的半峰宽的一半。而不是你先把半峰宽量出来再换算成eV。这个计算次序不能错。
这样你就能看出来0.01eV起码对于气相分子而言是很合理的（一般UV-Vis峰的宽度在0.1~0.3eV量级的比较多，但FCClasses的振动关联函数已经考虑了振动展宽，所以这里只需要考虑自然展宽和溶剂展宽，自然展宽一般远小于振动展宽，气相分子没有溶剂展宽），而40eV根本是不可能的

T9的梦

wzkchem5 发表于 2024-4-2 21:00
HWHM指的是你把谱的横坐标换算成eV以后，量出来的半峰宽的一半。而不是你先把半峰宽量出来再换算成eV。这 ...

谢谢您的详细解答！原来HWHM还有计算次序这样的细节问题，感觉豁然开朗！

cokie

楼主您好，请问您在做RISC计算时，SOC（in a.u.）是怎么得到的呢？是用比如ORCA计算出<T|HSO|S>值，除以219474.6363，转换成原子单位的吗？（比如ORCA算出的S1和T1之间的SOC是0.1 cm^-1，那么写在FCclasses里的SOC值就是 0.1/219474.6363 = 4.55E-7吗？）