求助通过Multiwfn计算q值（第2个q值）

yezhonghua

非常感谢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的计算方法理解有误。
谢谢大家的指导。

sobereva

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

yezhonghua

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。
谢谢！

sobereva

你先别自己优化，就直接用那篇文章补充材料里直接提供的PhCz的S1结构算个TDDFT，然后看看

yezhonghua

sobereva 发表于 2016-4-30 04:50
你先别自己优化，就直接用那篇文章补充材料里直接提供的PhCz的S1结构算个TDDFT，然后看看

恩恩，谢谢sober老师！