计算化学公社

标题: 激发态计算疑问 [打印本页]

作者
Author: 天天向上 时间: 2018-6-23 08:53
标题: 激发态计算疑问
老师，您好。我算出来一个分子的td,得到的out文件里面显示出的激发态的能量、波长、振子强度，很多表现出两两相同的情况，不知是什么原因导致的？为什么会产生这种情况？下面是输出文件的部分信息：
%chk=SFtd.chk
%nprocshared=24
Will use up to 24 processors via shared memory.
%mem=24GB
----------------------------------------------------
#p b3lyp/6-31g(d) td(singlet,direct,nstates=20) test
----------------------------------------------------
1/38=1/1;
2/12=2,17=6,18=5,40=1/2;
3/5=1,6=6,7=1,11=9,16=1,25=1,30=1,74=-5/1,2,8,3;
4//1;
5/5=2,38=5/2;
8/6=1,10=2,108=20/1;
9/41=20,42=1,48=1,70=2/14;
6/7=2,8=2,9=2,10=2/1;
99/5=1,9=1/99;
Leave Link 1 at Sat Mar 24 15:10:06 2018, MaxMem=  3221225472 cpu:       1.9
(Enter /share/apps/gauss/d01g09/g09/l101.exe)
-----
title

Excited State 1:    Singlet-A    4.2741 eV  290.09 nm  f=0.0526  <S**2>=0.000
   82 -> 85       -0.14857
   82 -> 86       -0.10904
   83 -> 84       0.62371
   83 -> 85       0.23234
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-KS) =  -962.172829440
Copying the excited state density for this state as the 1-particle RhoCI density.

Excited State 2:    Singlet-A    4.2741 eV  290.08 nm  f=0.0526  <S**2>=0.000
   82 -> 84       -0.14861
   82 -> 87       0.10904
   83 -> 84       -0.23233
   83 -> 85       0.62370

Excited State 3:    Singlet-A    4.5681 eV  271.41 nm  f=0.0027  <S**2>=0.000
   78 -> 84       -0.11357
   81 -> 84       0.14982
   81 -> 85       -0.14005
   82 -> 84       0.11634
   82 -> 87       -0.13811
   83 -> 86       0.57977
   83 -> 87       0.23850

Excited State 4:    Singlet-A    4.5681 eV  271.41 nm  f=0.0027  <S**2>=0.000
   78 -> 85       0.11357
   81 -> 84       0.14008
   81 -> 85       0.14982
   82 -> 85       -0.11635
   82 -> 86       -0.13816
   83 -> 86       -0.23851
   83 -> 87       0.57976

Excited State 5:    Singlet-A    4.6991 eV  263.85 nm  f=0.1675  <S**2>=0.000
   81 -> 84       -0.22030
   82 -> 84       0.52108
   82 -> 85       -0.37249

Excited State 6:    Singlet-A    4.6991 eV  263.85 nm  f=0.1675  <S**2>=0.000
   81 -> 85       0.22028
   82 -> 84       0.37251
   82 -> 85       0.52108

Excited State 7:    Singlet-A    4.8866 eV  253.72 nm  f=0.0032  <S**2>=0.000
   79 -> 84       0.31367
   79 -> 85       -0.17176
   80 -> 84       -0.17180
   80 -> 85       -0.31368
   82 -> 89       -0.14297
   83 -> 88       0.45139

Excited State 8:    Singlet-A    4.8982 eV  253.12 nm  f=0.0590  <S**2>=0.000
   81 -> 85       0.14677
   82 -> 86       -0.36572
   82 -> 87       0.48502
   83 -> 85       -0.13808
   83 -> 86       0.20379
   83 -> 87       -0.15300

Excited State 9:    Singlet-A    4.8982 eV  253.12 nm  f=0.0590  <S**2>=0.000
   81 -> 84       -0.14673
   82 -> 86       0.48501
   82 -> 87       0.36573
   83 -> 84       0.13807
   83 -> 86       0.15297
   83 -> 87       0.20382

Excited State  10:    Singlet-A    4.9317 eV  251.40 nm  f=0.0000  <S**2>=0.000
   79 -> 84       0.37631
   80 -> 85       0.37630
   82 -> 88       -0.34063
   83 -> 89       0.26763

Excited State  11:    Singlet-A    5.1633 eV  240.12 nm  f=0.0000  <S**2>=0.000
   79 -> 84       -0.10254
   79 -> 85       -0.48256
   80 -> 84       0.48452
   80 -> 85       -0.10361

Excited State  12:    Singlet-A    5.1731 eV  239.67 nm  f=0.0000  <S**2>=0.000
   79 -> 84       0.23883
   79 -> 85       0.43675
   80 -> 84       0.43458
   80 -> 85       -0.23838

Excited State  13:    Singlet-A    5.2707 eV  235.23 nm  f=0.1418  <S**2>=0.000
   78 -> 84       0.11604
   78 -> 85       -0.11362
   81 -> 84       0.59938
   82 -> 86       0.21942
   83 -> 87       -0.13407

Excited State  14:    Singlet-A    5.2707 eV  235.23 nm  f=0.1418  <S**2>=0.000
   78 -> 84       -0.11365
   78 -> 85       -0.11604
   81 -> 85       0.59938
   82 -> 87       -0.21945
   83 -> 86       0.13404

Excited State  15:    Singlet-A    5.3707 eV  230.85 nm  f=0.2273  <S**2>=0.000
   79 -> 84       -0.22459
   79 -> 85       0.12304
   80 -> 84       0.12301
   80 -> 85       0.22459
   82 -> 89       0.26869
   83 -> 88       0.52072

Excited State  16:    Singlet-A    5.4880 eV  225.92 nm  f=0.0000  <S**2>=0.000
   82 -> 88       0.47291
   83 -> 89       0.50005

Excited State  17:    Singlet-A    5.5651 eV  222.79 nm  f=0.0877  <S**2>=0.000
   78 -> 84       0.17225
   81 -> 86       0.64768

Excited State  18:    Singlet-A    5.5651 eV  222.79 nm  f=0.0877  <S**2>=0.000
   78 -> 85       -0.17232
   81 -> 87       0.64766

Excited State  19:    Singlet-A    5.5728 eV  222.48 nm  f=0.0000  <S**2>=0.000
   79 -> 86       0.39632
   79 -> 87       -0.29422
   80 -> 86       0.29459
   80 -> 87       0.39651

Excited State  20:    Singlet-A    5.6060 eV  221.16 nm  f=0.0000  <S**2>=0.000

作者
Author: sobereva 时间: 2018-6-23 12:33
说明那些态简并。
不要写direct、test这些多余的关键词

作者
Author: 天天向上 时间: 2018-6-23 14:10

sobereva 发表于 2018-6-23 12:33
说明那些态简并。
不要写direct、test这些多余的关键词

谢谢老师的回复。如果基态分子td计算出的是简并的，有什么应用上的优势或者劣势吗？我是研究太阳能电池受体材料的，是有利于激子的激发吗（简并时，相同能量下的态数量多）？

作者
Author: sobereva 时间: 2018-6-24 11:51

天天向上发表于 2018-6-23 14:10
谢谢老师的回复。如果基态分子td计算出的是简并的，有什么应用上的优势或者劣势吗？我是研究太阳能电池受 ...

没听过这种说法。而且这个的重点原理上也不在于是否存在简并
出现简并大多是由于体系对称性所致的

作者
Author: 让你变成回忆 时间: 2018-6-25 08:26
经常遇到一些D-A-D化合物（两个D单元是相同的）的激发态计算就会遇到如LZ说的激发态之间存在两两简并的情况，但是也不是像LZ这样简并的两个态能量完全相同。正如sob老师所说，简并往往就是体系具有对称性所导致的。

作者
Author: 天天向上 时间: 2018-6-25 09:55
了解了，谢谢！

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