计算化学公社

标题: 求助计算有机分子激发态的激发。 [打印本页]

作者
Author: chiweijie 时间: 2019-1-19 10:16
标题: 求助计算有机分子激发态的激发。
各位老师好，采用高斯和Multiwfn计算一有机分子激发态的激发，高斯计算的关键词为： # td(nstates=10) scrf=(smd,solvent=ethanol) def2svp m062x IOp(9/40=3)。按照Sob老师的博文介绍（http://sobereva.com/227），采用Multiwfn得到了激发态之间的跃迁偶极矩，如下：

“Transition dipole moment between excited states (a.u.):
   i    j       X                Y                      Z          Diff.(eV)    Oscil.str
   1    1    6.6930520 -0.0004925 -0.0001133    0.00000    0.00000
   1    2    0.0000789 -0.0078447    0.0164850    1.35460    0.00001
   1    3 -0.0000026 -0.0163784    0.0146710    2.10020    0.00002
   1    4 -0.0194240    0.0000487 -0.0000036    2.39550    0.00002
   1    5    0.0015585 -0.0000940    0.0000097    3.50470    0.00000
   1    6    0.4339689 -0.0000021    0.0000081    4.46160    0.02059
   1    7 -0.0002046 -1.7645267    0.0057565    4.95930    0.37830
   1    8 -0.0000960 -0.8302505 -0.0264709    5.11630    0.08649
   1    9    0.3903106    0.0000047 -0.0000100    5.16090    0.01926
   1 10 -0.0000140 -0.0726010 -0.0235884    5.22840    0.00075
   2    2    2.0134347 -0.0001701 -0.0002913    0.00000    0.00000
   2    3 -0.3366234    0.0000023 -0.0000250    0.74560    0.00207
   2    4 -0.0000572 -0.8057725 -0.0007305    1.04090    0.01656
   2    5    0.0002725    2.9643996    0.0010382    2.15010    0.46290
   2    6 -0.0000244    0.0033175    0.0196771    3.10700    0.00003
   2    7    0.0027156 -0.0000281    0.0000012    3.60470    0.00000
   2    8    0.0079419    0.0000489 -0.0000086    3.76170    0.00001
   2    9 -0.0000071 -0.0002489    0.0007227    3.80630    0.00000
   2 10    0.0279485 -0.0000071    0.0000068    3.87380    0.00007
   3    3    2.2574521 -0.0003694 -0.0002300    0.00000    0.00000
   3    4 -0.0003548 -3.9642306 -0.0021864    0.29530    0.11369
   3    5    0.0000156    0.1775511    0.0002441    1.40450    0.00108
   3    6    0.0000770    0.0011891    0.1323648    2.36140    0.00101
   3    7 -0.0006429 -0.0000036    0.0000003    2.85910    0.00000
   3    8 -0.0122153 -0.0000035 -0.0000005    3.01610    0.00001
   3    9    0.0000537    0.0030904 -0.0771828    3.06070    0.00045
   3 10    0.0072779    0.0000073 -0.0000031    3.12820    0.00000
   4    4    2.0926180 -0.0000191 -0.0002313    0.00000    0.00000
   4    5 -0.0370742 -0.0000082    0.0000318    1.10920    0.00004
   4    6 -0.0020525 -0.0001266 -0.0000106    2.06610    0.00000
   4    7 -0.0000066    0.0035254 -0.0339418    2.56380    0.00007
   4    8    0.0000261    0.0053784    0.1124084    2.72080    0.00084
   4    9 -0.0015696 -0.0000288 -0.0000142    2.76540    0.00000
   4 10    0.0000050    0.0232386 -0.0009327    2.83290    0.00004
   5    5    2.8643443 -0.0002454 -0.0001948    0.00000    0.00000
   5    6 -0.0014088    0.0000798    0.0000009    0.95690    0.00000
   5    7 -0.0000095 -0.0034194 -0.0439076    1.45460    0.00007
   5    8 -0.0000282 -0.0074292    0.1174744    1.61160    0.00055
   5    9    0.0039813    0.0000747 -0.0000030    1.65620    0.00000
   5 10    0.0000250    0.0823461 -0.0028848    1.72370    0.00029
   6    6    1.5714494 -0.0001261 -0.0002194    0.00000    0.00000
   6    7    0.0000017    0.2819698    0.0037061    0.49770    0.00097
   6    8 -0.0000182 -1.0209036 -0.0061980    0.65470    0.01672
   6    9    1.4159085 -0.0001438 -0.0000562    0.69930    0.03435
   6 10    0.0000003 -0.0032866 -0.0103742    0.76680    0.00000
   7    7    7.3475714 -0.0005463 -0.0000250    0.00000    0.00000
   7    8    1.1662207 -0.0000488    0.0000650    0.15700    0.00523
   7    9    0.0000455    0.0457304    0.0037148    0.20160    0.00001
   7 10    0.2007910 -0.0000033 -0.0001013    0.26910    0.00027
   8    8    4.1269464 -0.0004414 -0.0002536    0.00000    0.00000
   8    9 -0.0001623 -0.5222836 -0.0169114    0.04460    0.00030
   8 10    0.0830873 -0.0000115 -0.0000341    0.11210    0.00002
   9    9    1.9564536 -0.0002454 -0.0001962    0.00000    0.00000
   9 10 -0.0000043 -0.0022603    0.0014208    0.06750    0.00000
10 10    2.2095256    0.0000166 -0.0001821    0.00000    0.00000”

接下来该如何处理才能得到激发态的激发的能量（波长）呢？

作者
Author: 让你变成回忆 时间: 2019-1-19 13:02
Gaussian的输出文件中不就有激发态的能量和激发波长吗？

作者
Author: chiweijie 时间: 2019-1-19 13:12

让你变成回忆发表于 2019-1-19 13:02
Gaussian的输出文件中不就有激发态的能量和激发波长吗？

谢谢回复，我需要的激发态的激发，以激发态为参考态，再一次激发。

作者
Author: 让你变成回忆 时间: 2019-1-19 13:18
哦哦。这样啊！

作者
Author: sobereva 时间: 2019-1-19 20:15

chiweijie 发表于 2019-1-19 13:12
谢谢回复，我需要的激发态的激发，以激发态为参考态，再一次激发。

一般就是拿相应激发态的激发能求差来得到激发态之间能级差作为这个量

作者
Author: chiweijie 时间: 2019-1-20 14:42

sobereva 发表于 2019-1-19 20:15
一般就是拿相应激发态的激发能求差来得到激发态之间能级差作为这个量

多谢sob老师的回复，我计算的各个态的激发能如下：
Exc.state#    Exc.energy(eV)    Multi. MO pairs Normalization
   1          1.64790          1       137       0.499971
   2          3.00250          1       915       0.499974
   3          3.74810          1       811       0.499999
   4          4.04340          1       801       0.499972
   5          5.15260          1       1143       0.499987
   6          6.10950          1       383       0.499936
   7          6.60720          1       216       0.499929
   8          6.76420          1       448       0.499928
   9          6.80880          1       415       0.499926
   10          6.87630          1       218       0.499908
如果我想需要S1-S2的激发跃迁能，是否就是3.00250 -1.64790=1.3546 eV?

作者
Author: sobereva 时间: 2019-1-20 15:21

chiweijie 发表于 2019-1-20 14:42
多谢sob老师的回复，我计算的各个态的激发能如下：
Exc.state# Exc.energy(eV) Multi. MO pa ...

是

实际上直接取这个值就行了
1 2 0.0000789 -0.0078447 0.0164850 1.35460 0.00001

作者
Author: chiweijie 时间: 2019-1-20 17:00

sobereva 发表于 2019-1-20 15:21
是

实际上直接取这个值就行了

好的，多谢sob老师。

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