有机发光材料的FWHM计算求助

xuhj199508

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

心向暖阳

在单分子水平的话，重组能大，振动峰强，FWHM就宽
可是实验上还涉及更多的碰撞展宽的问题，这个应该与分子间相互作用有关系

wzkchem5

比较严格的方法是结合MOMAP或者用ORCA计算振动分辨光谱，再在振动分辨光谱上量取FWHM

xuhj199508

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计算振动分辨光谱是一回事吗？

wzkchem5

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不矛盾

胖鱼加加

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

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

wzkchem5

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

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

xuhj199508

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计算有影响吗？如果有影响，应该怎么消除？谢谢~

wzkchem5

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

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

小范范1989

黄昆因子和半高宽的关系呢？

xuhj199508

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

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

wzkchem5

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，如果有虚频，自动把分子往虚频的方向扰动，然后继续优化，直至没有虚频。如果虚频不是因为格点等数值误差导致的话，这个方法可以保证找到没有虚频的结构，而且完全不需要用户干预

xuhj199508

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

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

wzkchem5

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

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

xuhj199508

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

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