求助计算有机分子激发态的激发。

chiweijie · 发表于 Post on 2019-1-19 10:16:32

各位老师好，采用高斯和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”

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

让你变成回忆 · 发表于 Post on 2019-1-19 13:02:08

Gaussian的输出文件中不就有激发态的能量和激发波长吗？

chiweijie · 发表于 Post on 2019-1-19 13:12:05

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

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

让你变成回忆 · 发表于 Post on 2019-1-19 13:18:20

哦哦。这样啊！

sobereva · 发表于 Post on 2019-1-19 20:15:20

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

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

chiweijie · 发表于 Post on 2019-1-20 14:42:49

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?

sobereva · 发表于 Post on 2019-1-20 15:21:52

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

chiweijie · 发表于 Post on 2019-1-20 17:00:16

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

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

好的，多谢sob老师。

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[综合交流] 求助计算有机分子激发态的激发。

浏览过的版块