计算化学公社

标题: Multiwfn计算的跃迁电偶极矩和gaussian输出文件内的跃迁电偶极矩是同样的概念吗？ [打印本页]

作者
Author: knite 时间: 2023-11-25 18:03
标题: Multiwfn计算的跃迁电偶极矩和gaussian输出文件内的跃迁电偶极矩是同样的概念吗？
各位大佬，请问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

作者
Author: wzkchem5 时间: 2023-11-25 18:32
是不是分子的标准取向问题导致的？http://sobereva.com/297

作者
Author: knite 时间: 2023-11-25 18:53

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

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

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

作者
Author: knite 时间: 2023-11-26 22:03

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

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

作者
Author: sobereva 时间: 2023-11-27 03:23
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

作者
Author: knite 时间: 2023-11-27 08:47

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

好的，谢谢sob老师。

欢迎光临计算化学公社 (http://bbs.keinsci.com/)