Multiwfn计算的跃迁电偶极矩和gaussian输出文件内的跃迁电偶极矩是同样的概念吗？

knite

各位大佬，请问Multiwfn计算的跃迁电偶极矩和gaussian输出文件内的跃迁电偶极矩是同样的概念吗？
两个地方得到的基态与激发态的跃迁电偶极矩数值并不相等，是因为这二者不是同一概念还是其他原因的？
Multiwfn计算得到的数值如下：
Transition electric dipole moment between singlet states (a.u.):
   i        j          X             Y             Z         Diff.(eV)
S  0 --  S  0    -0.0000123     0.0391263     0.0000911      0.0000
S  0 --  S  1    -0.0000677    -1.5239175    -0.0000810      2.2717
S  0 --  S  2     2.5761388     0.0000990     0.5263812      2.5282
S  0 --  S  3    -1.7972609    -0.0000228     0.0979431      2.6117
S  0 --  S  4    -0.0000632     0.3821283     0.0000457      2.8848
S  0 --  S  5     2.8921224    -0.0001207    -0.8714163      2.9004
S  0 --  S  6    -0.5709858    -0.0000095    -0.0746964      2.9837
S  0 --  S  7    -0.0000336    -0.1063722    -0.0000287      3.0316
S  0 --  S  8     0.0004794     0.1255494     0.0000067      3.1899
S  0 --  S  9     0.0000829    -0.3183567    -0.0001884      3.2343
S  0 --  S 10     0.0833302     0.0000744    -0.5167079      3.2604


gaussian输出文件内的电偶极矩如下（因为用了50-50，里面含有三重态的部分，对应偶极矩为0）：
Ground to excited state transition electric dipole moments (Au):
       state          X           Y           Z        Dip. S.      Osc.
         1         0.0000     -0.0000     -0.0000      0.0000      0.0000
         2        -0.0000     -0.0000      0.0000      0.0000      0.0000
         3         0.0000     -1.3386      0.0000      1.7920      0.0997
         4         0.0000      0.0000     -0.0000      0.0000      0.0000
         5        -0.0000     -0.0000     -0.0000      0.0000      0.0000
         6         2.3683      0.0000      0.3646      5.7417      0.3556
         7         0.0000     -0.0000      0.0000      0.0000      0.0000
         8        -1.6204     -0.0000      0.0685      2.6303      0.1683
         9        -0.0000      0.0000      0.0000      0.0000      0.0000
        10        -0.0000      0.0000     -0.0000      0.0000      0.0000
        11         0.0000     -0.0000     -0.0000      0.0000      0.0000
        12         0.0000     -0.0000      0.0000      0.0000      0.0000
        13         0.0000     -0.0000     -0.0000      0.0000      0.0000
        14        -0.0000      0.3313     -0.0000      0.1098      0.0078
        15         2.6928     -0.0000     -0.6117      7.6252      0.5418
        16        -0.5344      0.0000     -0.0567      0.2888      0.0211
        17        -0.0001     -0.1215     -0.0000      0.0148      0.0011
        18        -0.0001      0.1100     -0.0001      0.0121      0.0009
        19         0.0000     -0.2873      0.0000      0.0825      0.0065
        20         0.0676      0.0005     -0.3547      0.1304      0.0104

wzkchem5

是不是分子的标准取向问题导致的？http://sobereva.com/297

knite

wzkchem5 发表于 2023-11-25 18:32
是不是分子的标准取向问题导致的？http://sobereva.com/297

如果是朝向问题的话应该模一样吧，但似乎这两个偶极矩的模是不一样的

sobereva

别用50-50，掺和在一起都不好对比
别一会儿偶极矩一会儿跃迁偶极矩，完全是不同的概念
确保写了IOp(9/40=4)。如果写了，把Gaussian的fch和out文件压缩后上传，到时候我告诉你你哪弄错了

knite

sobereva 发表于 2023-11-26 02:47
别用50-50，掺和在一起都不好对比
别一会儿偶极矩一会儿跃迁偶极矩，完全是不同的概念
确保写了IOp(9/40= ...

因为文件有些大，所以我上传到网盘了，麻烦sob老师解答一下，十分感谢！！输出文件

sobereva

IOp(9/40=3)输出的组态系数不够多
我用你的结构重新做了TDDFT计算，用了IOp(9/40=4)，只算单重态激发态，关键词#p td(nstates=10) b3lyp/def2svp em=gd3bj scrf=(smd,solvent=toluene) iop(9/40=4) nosymm

Gaussian输出结果为
Ground to excited state transition electric dipole moments (Au):
       state          X           Y           Z        Dip. S.      Osc.
         1         0.0000     -1.3386      0.0000      1.7920      0.0997
         2         2.3683      0.0000      0.3645      5.7415      0.3556
         3        -1.6206     -0.0000      0.0686      2.6309      0.1683
         4        -0.0000      0.3309     -0.0000      0.1095      0.0077
         5         2.6929     -0.0000     -0.6118      7.6260      0.5419
         6        -0.5343      0.0000     -0.0566      0.2887      0.0211
         7        -0.0001     -0.1212     -0.0000      0.0147      0.0011
         8        -0.0000      0.1096     -0.0000      0.0120      0.0009
         9         0.0000     -0.2884      0.0000      0.0832      0.0066
        10         0.0679      0.0005     -0.3547      0.1304      0.0104

Multiwfn（目前官网上最新版）结果为
Transition electric dipole moment between ground state (0) and excited states (a.u.)
     i     j         X             Y             Z        Diff.(eV)   Oscil.str
     0     1     0.0002211    -1.3683808    -0.0000895     2.27170     0.10421
     0     2     2.4039406     0.0000802     0.3963247     2.52820     0.36767
     0     3    -1.6589634     0.0000607     0.0785839     2.61170     0.17649
     0     4     0.0000285     0.3452126    -0.0000619     2.88480     0.00842
     0     5     2.7178414    -0.0001899    -0.6492479     2.90040     0.55484
     0     6    -0.5460487    -0.0000142    -0.0564698     2.98370     0.02203
     0     7    -0.0000395    -0.1196411    -0.0000451     3.03160     0.00106
     0     8    -0.0001559     0.1073801     0.0001215     3.18990     0.00090
     0     9     0.0001640    -0.2937128     0.0001138     3.23430     0.00684
     0    10     0.0648073     0.0009447    -0.3878216     3.26040     0.01235

可见已经挺相符了。如果需要吻合程度更高，需要9/40=5

knite

sobereva 发表于 2023-11-27 03:23
IOp(9/40=3)输出的组态系数不够多
我用你的结构重新做了TDDFT计算，用了IOp(9/40=4)，只算单重态激发态， ...

好的，谢谢sob老师。