计算化学公社

标题: FCClasses3 计算RISC速率为负数 [打印本页]

作者
Author: lennylyp 时间: 2023-6-25 00:06
标题: FCClasses3 计算RISC速率为负数
本帖最后由 lennylyp 于 2023-6-25 00:09 编辑

想用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

作者
Author: lennylyp 时间: 2023-6-25 10:15
额外发现温度设置得低一点，可以避免出现负数。但这又是什么原因导致的呢？

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

作者
Author: 15689900879@163 时间: 2024-2-19 17:07
本帖最后由 15689900879@163 于 2024-2-21 14:33 编辑

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

作者
Author: ABQTrap 时间: 2024-2-26 20:55
本帖最后由 ABQTrap 于 2024-2-26 20:57 编辑

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

作者
Author: Eclair 时间: 2024-2-29 15:03

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

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

作者
Author: ABQTrap 时间: 2024-3-2 17:41

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

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

作者
Author: wzkchem5 时间: 2024-3-3 05:55

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

作者
Author: T9的梦 时间: 2024-4-2 19:53

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）感觉非常迷茫

作者
Author: wzkchem5 时间: 2024-4-2 21:00

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

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

作者
Author: T9的梦 时间: 2024-4-3 10:34

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

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

作者
Author: cokie 时间: 2024-5-17 21:36
楼主您好，请问您在做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吗？）

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