如何使用MOMAP软件计算S1→S0的荧光辐射速率

描述

在使用量化软件BDF完成Azulene的基态和激发态的结构优化、频率计算以及非绝热耦合的计算,并计算完MOMAP输入文件中需要的参数后,接下来将使用MOMAP软件对azulene的S1→S0的辐射速率和内转换速率进行计算,通过对比辐射速率和内转换速率来解释azulene的S1→S0的荧光难以被观测到的原因。 首先计算S1→S0的荧光辐射速率,第一步为做电子振动耦合(electron-vibration coupling, EVC)计算,该计算基于量化计算输出的分子振动频率、力常数矩阵,同时在内坐标以及直角坐标系下,计算分子跃迁发生初末态间的模式位移、黄昆因子、重整能以及Duschinsky转动矩阵。 将bdf的S0优化频率计算结果改为azulene-s0.out,将S1的计算结果改为azulene-s1.out,同时放在EVC计算文件夹中。

EVC计算的输入文件momap.inp为:

do_evc = 1

&evc

ffreq(1)      = "azulene-s0.out"

ffreq(2)      = "azulene-s1.out"

proj_reorg = .t.

/

提交脚本文件momap.slurm,任务运行结束后,生成如下文件

计算

其中 evc.cart.dat 为利用直角坐标系计算得到的模式位移、黄昆因子、重整能以及 Duschinsky 转动矩阵的结果;而evc.dint.dat 为内坐标计算模式位移、黄昆因子、重整能,直角坐标计算Duschinsky转动矩阵的结果。

打开evc.cart.dat文件,查看总重整能的数值:

计算

并与evc.dint.dat文件中的数值进行比较:

计算

比较 evc.cart.dat 以及 evc.dint.dat 文件内的重整能,若重整能相差不大(< 1000 cm-1),使用evc.cart.dat文件进行后续计算,若重整能相差较大(一般情况evc.cart.dat大于evc.dint.dat),使用evc.dint.dat文件进行后续计算。这里两者较为接近,且数值较为合理(<10000 cm-1),可使用evc.cart.dat进行下一步S1→S0的荧光辐射速率计算。 此外,还可以根据EVC计算的结果文件做更多的后处理。 evc.dx.x.com和evc.dx.x.xyz为两个电子态分子叠加图,其中evc.dx.x.com可用GaussView打开,在View-Display Format-Molecule中选择Tube类型,显示如下:

计算

evc.dx.x.xyz可用Jmol软件打开,显示如下:

计算

evc.dx.v.xyz也是分子态叠加图,用Jmol软件打开,选择显示-矢量-0.1A,显示如下:

计算

evc.cart.abs为Duschinsky矩阵文件,可用来画Duschinsky矩阵二维图。可以在Device Studio中,选择Simulator-momap-analysis,打开evc.cart.abs文件,显示如下:

计算

同样的方法打开evc.dint.abs文件,可在ColorType下拉选框选择显示颜色,显示如下:

计算

Duschinsky矩阵都是用直角坐标计算的,二者选一即可。 将evc.vib1.xyz和evc.vib2.xyz文件与evc.cart.dat文件放在同一路径下,在Device Studio中,选择Simulator-momap-analysis,打开evc.cart.dat文件,出现基态和激发态下各个振动模式对应的重整能和黄昆因子图,显示如下:

计算

计算

计算

计算

  点击Choose Color可以改变线的颜色,改变Set Width中的数值可改变线的粗细。点击图中的线,右侧结构将展示相应的振动模式。振动的快慢可通过Animation Frequency调整,振动的幅度可根据Displacement Amplitude显示。   

      审核编辑:彭静
打开APP阅读更多精彩内容
声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉

全部0条评论

快来发表一下你的评论吧 !

×
20
完善资料,
赚取积分