计算化学公社

标题: 求助通过Multiwfn计算q值(第2个q值) [打印本页]

作者
Author:
yezhonghua    时间: 2016-4-29 00:06
标题: 求助通过Multiwfn计算q值(第2个q值)
非常感谢sober老师在q值计算过程中给与的帮助和指导。
根据sober老师的指导,第一个q值(用来计算EVA),其和文献完全一致,如图(1)所示。即1-21(donor)对于HOMO、LUMO轨道贡献的差值。
但是在计算第2个q值(EVE,甲苯溶液)时遇到了些问题。第2个q(图1描红)按照文献的说法(如图2),其计算是通过S1 opt得到,采用的泛函为EVA-HF%中的optimal HF,即为MPWLYP1M(5%)。根据第一个q值的计算方法,计算结果为0.89和文献的0.186相差巨大。
不知道是什么原因引起,是否是我的第2个q的计算方法理解有误。
谢谢大家的指导。

作者
Author:
sobereva    时间: 2016-4-29 03:30
看一下S1结构下S1-S0跃迁对应的各对轨道贡献是多少。
从当前fch文件看,HOMO和LUMO空间分布相差巨大,所以拿它算出来的q显然很接近于1。
用HOMO、LUMO来算q的基本前提是S1-S0电子态跃迁可以基本近似用HOMO-LUMO跃迁来表示,如果这个近似不成立,那就得用相应的其它的主导的轨道对来计算q。

作者
Author:
yezhonghua    时间: 2016-4-29 16:56
sobereva 发表于 2016-4-29 03:30
看一下S1结构下S1-S0跃迁对应的各对轨道贡献是多少。
从当前fch文件看,HOMO和LUMO空间分布相差巨大,所以 ...

sober老师,感谢你的解答。
           S0-S1的跃迁是64-66,就是HOMO-LUMO+1;S1 opt (S1-S0)的跃迁轨道为从开始的64-66优化到最后的64-65(最后一个S1)。
也就是S1-S0的跃迁确实是homo-lumo。
无论是64-65,还是64-66的q值都不是0.186。
谢谢!

作者
Author:
sobereva    时间: 2016-4-30 04:50
你先别自己优化,就直接用那篇文章补充材料里直接提供的PhCz的S1结构算个TDDFT,然后看看
作者
Author:
yezhonghua    时间: 2016-5-3 17:53
sobereva 发表于 2016-4-30 04:50
你先别自己优化,就直接用那篇文章补充材料里直接提供的PhCz的S1结构算个TDDFT,然后看看

恩恩,谢谢sober老师!




欢迎光临 计算化学公社 (http://ccc.keinsci.com/) Powered by Discuz! X3.3