激发态计算疑问

天天向上

老师，您好。我算出来一个分子的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

sobereva

说明那些态简并。
不要写direct、test这些多余的关键词

天天向上

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

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

sobereva

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

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

让你变成回忆

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

天天向上

了解了，谢谢！