计算化学公社

标题: 有机发光材料的FWHM计算求助 [打印本页]

作者
Author: xuhj199508 时间: 2021-6-24 17:22
标题: 有机发光材料的FWHM计算求助
各位老师好：
最近在研究TADF客体材料半峰宽的计算。之前的计算方法是利用Gaussian通过四点法求出图中激发和发射部分的重组能λs和λs*，计算级别为PBE0/6-31G*。
我计算的λs与材料实测的FWHM的确呈正相关关系，但是计算值与实测值的拟合度较低，目前只能用作定性判断。请问大家有什么更好的方法可以更准确地预测FWHM，谢谢？

作者
Author: 心向暖阳 时间: 2021-6-24 18:05
在单分子水平的话，重组能大，振动峰强，FWHM就宽
可是实验上还涉及更多的碰撞展宽的问题，这个应该与分子间相互作用有关系

作者
Author: wzkchem5 时间: 2021-6-24 18:47
比较严格的方法是结合MOMAP或者用ORCA计算振动分辨光谱，再在振动分辨光谱上量取FWHM

作者
Author: xuhj199508 时间: 2021-6-25 18:26

wzkchem5 发表于 2021-6-24 18:47
比较严格的方法是结合MOMAP或者用ORCA计算振动分辨光谱，再在振动分辨光谱上量取FWHM

老师你好，感谢您的回答。我试过使用MOMAP 1.0和G09算光谱。
我参考了http://sobereva.com/223这篇帖子，g09算振动分辨光谱默认SpecHwHm=135 cm-1，或者自行设置。我所使用的MOMAP版本也是需要设置计算时的FWHM。
我也有使用ORCA算过SOC校正光谱http://bbs.keinsci.com/thread-12114-1-1.html，请问和您说的ORCA计算振动分辨光谱是一回事吗？

作者
Author: wzkchem5 时间: 2021-6-25 18:32

xuhj199508 发表于 2021-6-25 11:26
老师你好，感谢您的回答。我试过使用MOMAP 1.0和G09算光谱。
我参考了http://sobereva.com/223这篇帖子 ...

SOC校正光谱和振动分辨光谱完全没有关系。ORCA里面振动分辨光谱是用ESD模块做的，sob老师不一定写过这方面的博文，应该去看ORCA说明书的ESD章节。
作振动光谱时设置的FWHM和你要算的FWHM也没有关系，前者是单个振动能级峰的FWHM，后者是所有振动能级所形成的包络线的FWHM，你需要的是后者。设置前者对后者有影响，但影响极小，可能也就在1%量级。所以，你需要设FWHM和你可以从振动光谱算出FWHM不矛盾

作者
Author: 胖鱼加加 时间: 2021-8-20 12:33

wzkchem5 发表于 2021-6-25 18:32
SOC校正光谱和振动分辨光谱完全没有关系。ORCA里面振动分辨光谱是用ESD模块做的，sob老师不一定写过这方 ...

前辈，看了上文的讨论，我在这里有两个问题向您请教：
1. ORCA计算的振动分辨光谱和Gaussian计算的振动分辨的电子光谱有什么区别和联系？
2.实测PL光谱如果是固态下的测量结果，那么是不是计算用QM/MM方法模拟固态下单分子PL会比模拟气相单分子更靠谱？

作者
Author: wzkchem5 时间: 2021-8-20 16:00

胖鱼加加发表于 2021-8-20 05:33
前辈，看了上文的讨论，我在这里有两个问题向您请教：
1. ORCA计算的振动分辨光谱和Gaussian计算的振动 ...

1. ORCA提供的计算方法更多，其中最严格的那种应该和高斯等价（是否有微小区别我就不知道了），但是当激发态结构和基态差别很大、激发态有的振动频率特别小、激发态没有平衡结构等情况下，高斯的那种计算方法可能做不了，此时只能用ORCA的那些相对不太严格的，但是在这些情况下也能算的方法
2. 对

作者
Author: xuhj199508 时间: 2021-8-24 15:12

wzkchem5 发表于 2021-6-25 18:32
SOC校正光谱和振动分辨光谱完全没有关系。ORCA里面振动分辨光谱是用ESD模块做的，sob老师不一定写过这方 ...

老师你好，我算了几个分子的ESD_AH，有一个分子的拟合曲线有点问题，想请教一下老师。
计算GS，ES，ESD都正常结束。查看了ES的out文件，发现在计算振动频率时有一个虚频如下：
-----------------------
VIBRATIONAL FREQUENCIES
-----------------------

Scaling factor for frequencies =  1.000000000 (already applied!)

0:       0.00 cm**-1
1:       0.00 cm**-1
2:       0.00 cm**-1
3:       0.00 cm**-1
4:       0.00 cm**-1
5:       0.00 cm**-1
6:    -61.00 cm**-1 ***imaginary mode***
7:       9.37 cm**-1
请问这个虚频会后续ESD计算有影响吗？如果有影响，应该怎么消除？谢谢~

作者
Author: wzkchem5 时间: 2021-8-24 16:23

xuhj199508 发表于 2021-8-24 08:12
老师你好，我算了几个分子的ESD_AH，有一个分子的拟合曲线有点问题，想请教一下老师。
计算GS，ES，ESD ...

和基态消虚频的方法一样，参见http://sobereva.com/278
对于ESD计算，仅当用AHAS、VG等不基于激发态平衡结构的方法时，激发态可以有虚频，否则不行

作者
Author: 小范范1989 时间: 2021-8-24 16:33
黄昆因子和半高宽的关系呢？

作者
Author: xuhj199508 时间: 2021-8-25 14:17

wzkchem5 发表于 2021-8-24 16:23
和基态消虚频的方法一样，参见http://sobereva.com/278
对于ESD计算，仅当用AHAS、VG等不基于激发态平衡 ...

老师，按照类似博文中的方法，我在输入文件里加入tightopt def2grid关键词是否可以？
ORCA中类似于G16下opt=recalc=n的关键词是什么？

作者
Author: wzkchem5 时间: 2021-8-25 15:39

xuhj199508 发表于 2021-8-25 07:17
老师，按照类似博文中的方法，我在输入文件里加入tightopt def2grid关键词是否可以？
ORCA中类似于G16下 ...

tightopt可以尝试。但是def2grid没用，因为默认格点就是def2grid，如果怀疑是格点导致虚频，应该写defgrid3。opt=recalc=n的等价做法是在%geom里写calc_hess true recalc_hess n。
不过比较推荐的做法是用orca的ITERATIVE_OPTIMIZATION方法（从orca 5.0开始才支持，在手册的compound methods部分），这个方法先做一次普通的opt+freq，如果有虚频，自动把分子往虚频的方向扰动，然后继续优化，直至没有虚频。如果虚频不是因为格点等数值误差导致的话，这个方法可以保证找到没有虚频的结构，而且完全不需要用户干预

作者
Author: xuhj199508 时间: 2021-8-25 18:13

wzkchem5 发表于 2021-8-25 15:39
tightopt可以尝试。但是def2grid没用，因为默认格点就是def2grid，如果怀疑是格点导致虚频，应该写defgri ...

谢谢老师，我想尝试使用ITERATIVE_OPTIMIZATION方法。我看了手册，他是用! compound[protocol name]来调用ITERATIVE_OPTIMIZATION，可是我把这个加到input中后，一直无法计算。请问这个功能应该怎么正确地调用？

作者
Author: wzkchem5 时间: 2021-8-25 18:55

xuhj199508 发表于 2021-8-25 11:13
谢谢老师，我想尝试使用ITERATIVE_OPTIMIZATION方法。我看了手册，他是用! compound来调用ITERATIVE_OPTI ...

这里的问题在于ITERATIVE_OPTIMIZATION方法不能读取用户的理论级别等信息，你必须把ITERATIVE_OPTIMIZATION的compound method脚本（.cmp）拷到你输入文件的目录下，把你的理论级别和%tddft块等的信息写进那个.cmp文件，然后调用

作者
Author: xuhj199508 时间: 2021-8-27 12:49

wzkchem5 发表于 2021-8-25 18:55
这里的问题在于ITERATIVE_OPTIMIZATION方法不能读取用户的理论级别等信息，你必须把ITERATIVE_OPTIMIZATI ...

老师你好，我创建了test.inp和iterativeOptimization.cmp，在.cmp文件中写入了计算级别和核数。用orca test.inp > test.out 提交了计算，可正常运行，但是只有一个核在跑,不能识别cmp文件中的%pal nprocs信息吗？

作者
Author: wzkchem5 时间: 2021-8-27 16:19

xuhj199508 发表于 2021-8-27 05:49
老师你好，我创建了test.inp和iterativeOptimization.cmp，在.cmp文件中写入了计算级别和核数。用orca te ...

对，并行核数必须写在主输入文件里，写在cmp文件里不认

作者
Author: xuhj199508 时间: 2021-9-25 10:14

wzkchem5 发表于 2021-6-24 18:47
比较严格的方法是结合MOMAP或者用ORCA计算振动分辨光谱，再在振动分辨光谱上量取FWHM

老师，还有个问题想请教下您:
请问如何根据计算得到的UV/PL光谱分析每个尖峰是由分子结构中的哪部分区域的振动伸缩贡献得到？

作者
Author: wzkchem5 时间: 2021-9-25 16:20

xuhj199508 发表于 2021-9-25 03:14
老师，还有个问题想请教下您:
请问如何根据计算得到的UV/PL光谱分析每个尖峰是由分子结构中的哪部分区域 ...

一般不知道，因为这个谱是在时域求解，再傅里叶变换过来的，不是一个峰一个峰地求解的，所以每个峰代表什么这个信息就丢失了。
不过一般而言，不是每个峰对应一个单独的振动，很多峰对应的是同一个振动的不同的倍频。如果看到每个峰之间的距离基本是相等的，那么可以量一下峰之间的距离，换算成波数，可以知道对应的是一个多大频率的振动模式

作者
Author: zhaixiaoyi001 时间: 2023-2-14 10:23
@wzkchem5 老师，我有一个关于ORCA 5.0.3 ESD 荧光模块使用 PRINTLEVEL 4选项输出的黄昆因子的问题想要请教一下:
1. 这里输出的黄昆因子对应的是基态还是激发态的振动模式?
2. 这里输出的黄昆因子大小是否可以衡量S1态无辐射能量耗散过程的强度?

作者
Author: NGC626 时间: 2023-4-5 10:21

xuhj199508 发表于 2021-8-27 12:49
老师你好，我创建了test.inp和iterativeOptimization.cmp，在.cmp文件中写入了计算级别和核数。用orca te ...

楼主您好，我最近在给一个结构做optfreq的时候也出现了两个虚频，一直无法消除，我也想用ITERATIVE_OPTIMIZATION方法消除虚频，但是修改文件之后一直提示无法运行，请问应如何修改inp文件和cmp文件呢？这是我的输入文件夹

作者
Author: wzkchem5 时间: 2023-4-5 18:09

zhaixiaoyi001 发表于 2023-2-14 03:23
@wzkchem5 老师，我有一个关于ORCA 5.0.3 ESD 荧光模块使用 PRINTLEVEL 4选项输出的黄昆因子的问题想要请教 ...

抱歉之前没看到。最好在orca论坛上问一下或者搜一下

作者
Author: zhaixiaoyi001 时间: 2023-4-11 19:41

wzkchem5 发表于 2023-4-5 18:09
抱歉之前没看到。最好在orca论坛上问一下或者搜一下

感谢老师回答！我现在的理解是黄昆因子对应的是基态和激发态振动模式的混合，并不是对应的某一个态的振动模式，单凭黄昆因子的大小无法衡量无辐射跃迁速率，还需要考虑NAC等问题。

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