Gaussian结合Multiwfn计算电荷转移量qCT相关问题

小螃蟹

各位大神，我使用Gaussian09软件包计算电荷转移量（qCT）,具体的输入如下所示：计算基态：
%chk=1-3-M062X-S0.chk                                                   
%nproc=16                                                               
%mem=500MW                                                              
#p M062X/6-31g(d,p) out=wfn nosymm scrf=(PCM,solvent=tetrahydrofuran)   

Title Card Required                                                     

0 1                                                                     
  6           0        6.955901    3.558118   -1.217909                 
  6           0        5.812564    3.109408   -1.884939                 
  6           0        5.218109    1.899656   -1.517672                 
......
/public/home/hengpp/DSSC5/1-3-M062X-S0.wfn
计算激发态：   
%chk=1-3-M062X-S1.chk                                                                           
%nproc=16                                                                                       
%mem=500MW                                                                                       
#p M062X/6-31g(d,p) out=wfn nosymm scrf=(PCM,solvent=tetrahydrofuran) td(root=1,singlet) density

Title Card Required                                                                              

0 1                                                                                             
  6           0        6.955901    3.558118   -1.217909                                          
  6           0        5.812564    3.109408   -1.884939                                          
  6           0        5.218109    1.899656   -1.517672
......
/public/home/hengpp/DSSC5/1-3-M062X-S1.wfn   
计算完成后，将生成的.wfn文件在Multiwfn中运行可得出qCT具体数值以及可用于做电荷密度差图的.cub文件。
但是，现在我有一些疑问如下，还请各位大神指点
1，我在计算基态与激发态时都是用的基态下的分子构型坐标（这是不是正确的 ） ，还是说需要对激发态优化用优化好的激发态构型再去计算，
2，在计算电荷转移量时，是否是默认的从基态到第一激发态的电荷转移（即S0-S1）。
谢谢各位啦                        

千域小少

1，都用基态构型没问题
2，root=1，你指定的就是第一激发态的密度
参见 http://sobereva.com/113

sobereva

两个wfn文件的坐标必须严格对应，否则基于密度差的分析没有意义。考察垂直吸收过程就用基态下极小点结构的坐标，考察垂直发射过程就用相应激发态极小点的坐标

小螃蟹

千域小少 发表于 2019-6-4 16:51
1，都用基态构型没问题
2，root=1，你指定的就是第一激发态的密度
参见 http://sobereva.com/113

好的 谢谢你，我会好好看看链接内容的。

小螃蟹

sobereva 发表于 2019-6-4 21:50
两个wfn文件的坐标必须严格对应，否则基于密度差的分析没有意义。考察垂直吸收过程就用基态下极小点结构的 ...

嗯~  谢谢sobereva大神的指点

quanjiezhong

你好，请问用Gaussian分别计算了基态和激发态能量，怎么通过Multiwfn计算qCT的具体数值。

sobereva

quanjiezhong 发表于 2019-6-17 23:22
你好，请问用Gaussian分别计算了基态和激发态能量，怎么通过Multiwfn计算qCT的具体数值。

光有能量显然没用

qCT的例子看Multiwfn手册4.18.3节

更多信息看
Multiwfn支持的电子激发分析方法一览
http://sobereva.com/437（http://bbs.keinsci.com/thread-11406-1-1.html）

weixy

你好,我想请问下啊你这个具体是咋做的?如何得到你说的电荷转移量qCT

YOUTH-ZTR

请问 为何我用s0结构下的td和opt文件作差，得到的qCT为0，但电荷转移距离不为零。

wzkchem5

YOUTH-ZTR 发表于 2023-4-22 04:41
请问 为何我用s0结构下的td和opt文件作差，得到的qCT为0，但电荷转移距离不为零。

提问格式不符合论坛要求，按http://sobereva.com/620的要求仔细修改、反复检查确保符合格式，我们再回答

YOUTH-ZTR

wzkchem5 发表于 2023-4-22 14:52
提问格式不符合论坛要求，按http://sobereva.com/620的要求仔细修改、反复检查确保符合格式，我们再回答

1、我利用上图两个算法基组先后计算得到了分子的基态构型和基态构型下的垂直跃迁能，并用计算得到的fchk文件通过Multiwfn的100-2-5功能转换为wfn文件格式。
2、随后在Multiwfn中拖入td的wfn文件，并按照手册中4.18.3中的步骤进行操作，希望计算得到基态构型下s0-s1的qCT。
3、但最终结果显示qCT为0，但dCT不为0，请问是我在哪一步计算中出现了问题。感谢。