利用Multiwfn计算Gaussian输出的激发态之间的跃迁偶极矩

Phil宏爷

本人参照教程利用Multiwfn计算Gaussian输出的激发态之间的跃迁偶极矩时，遇到以下问题，希望可以请教一下大家：1.我用Gaussian输入命令：# B3LYP/6-31+G(d) TD(nstates=10) IOp(9/40=4)，请问这是对于单线态的研究还是三线态的研究。


2.怎么判断我的体系需要做多少个激发态？（教程建议是用10）
3.我的结果如下，怎么看偶极矩？S0-S1的跃迁偶极矩又该怎么计算？


Ground state dipole moment in X,Y,Z:    1.027343    0.886573    0.182031 a,u,


Transition dipole moment between excited states (a.u.):
     i     j         X             Y             Z        Diff.(eV)   Oscil.str
     1     1    -4.4257679     3.2016113    -0.1358897     0.00000     0.00000
     1     2     2.3122051    -1.6883880    -0.0260235     0.12210     0.02452
     1     3     1.1233449    -0.2004368    -0.0079681     0.25950     0.00828
     1     4     2.8638589    -0.1060974    -0.0108409     0.53800     0.10825
     1     5    -3.1581674     0.7116104    -0.1030120     0.64760     0.16645
     1     6    -0.0672630     0.1473130    -0.0022858     0.83220     0.00053
     1     7    -0.8022089     0.4011351    -0.0533246     0.88940     0.01759
     1     8    -0.2523473    -0.5977021     0.1904789     1.14140     0.01279
     1     9     1.2785860    -0.2889457     0.0783643     1.16730     0.04932
     1    10    -2.4453258    -0.6067728     0.2524216     1.25900     0.19776
     2     2     0.3214683    -0.2979580    -0.0731896     0.00000     0.00000
     2     3    -0.4005476     0.0761722    -0.0172813     0.13740     0.00056
     2     4    -1.2936675     0.0417576    -0.0208220     0.41590     0.01707
     2     5     1.0131020    -0.1508747    -0.0118191     0.52550     0.01351
     2     6     0.5045086    -0.0328673     0.0640793     0.71010     0.00452
     2     7     0.2273081    -0.1467913     0.0208487     0.76730     0.00138
     2     8     0.5647708     0.2186490     0.0454146     1.01930     0.00921
     2     9    -1.3744657     0.0680653    -0.2624825     1.04520     0.05026
     2    10     0.9201637     0.0282574    -0.1191709     1.13690     0.02400
     3     3    -1.3779219     2.2388007     0.0469366     0.00000     0.00000
     3     4    -2.3536721     0.6005828    -0.0418020     0.27850     0.04027
....
希望能有老师能解答我的疑问，谢谢

sobereva

1 如果你是在0 1下算的，基态是单重态，你这样的关键词算出来的激发态也都是单重态。
2 算多少取决于目的。比如你就算第一激发能，那你算5个态就够，如果是大体系要算整个UV-Vis光谱，那就得算几十个甚至上百个态，才能把整个UV-Vis光谱感兴趣的区域覆盖。
3 基态偶极矩就是第一行（高斯输出文件末尾也能看到）。
后面的输出是激发态之间的跃迁偶极矩。S0到S1的（以及S0到其它激发态的）跃迁偶极矩在你的高斯输出文件中直接就能看到，不需要Multiwfn。只有激发态之间的TDDFT跃迁偶极矩才是必须要用Multiwfn的。

Phil宏爷

谢谢老师。

huangyoyo

sobereva 发表于 2015-9-30 00:28
1 如果你是在0 1下算的，基态是单重态，你这样的关键词算出来的激发态也都是单重态。
2 算多少取决于目的 ...

用ORCA的输出文件（%tddft   nroots 200; tda flase; end）来产生各个态之间的跃迁偶极矩，得到的基态到各个激发态之间的振子强度与输出文件中的结果相差很大，是不是因为我没有设置Tprint导致或者没有考虑去激发的影响？

sobereva

huangyoyo 发表于 2019-5-8 09:57
用ORCA的输出文件（%tddft   nroots 200; tda flase; end）来产生各个态之间的跃迁偶极矩，得到的基态到 ...

不设tprint肯定非常不准。如果设了，结果还和直接输出的差异大，一般就是去激发没考虑的事

huangyoyo

sobereva 发表于 2019-5-9 08:55
不设tprint肯定非常不准。如果设了，结果还和直接输出的差异大，一般就是去激发没考虑的事

能用Multiwfn将跃迁偶极矩的数据画个光谱图么？

sobereva

huangyoyo 发表于 2019-5-21 16:20
能用Multiwfn将跃迁偶极矩的数据画个光谱图么？

跃迁偶极矩本身没法用来画图，画图得先转成振子强度。通过能量差和跃迁偶极矩很简单就能算出振子强度，把数据写成Multiwfn能认的文本文件形式，就可以作为输入文件画图