TDDFT没有轨道越迁信息

一颗赛艇

大家好！最近再算一个氙团簇的TDDFT，用的M06-2X，从单重态到三重态的激发，但是算完了发现前四个越迁没有轨道信息。请问这是哪里出错了？谢谢

1. ***********************************************************************
   
2. Excited states from <AA,BB:AA,BB> singles matrix:
   
3. ***********************************************************************
   
4. 
   
5. Ground to excited state transition densities written to RWF  633
   
6. Ground to excited state transition electric dipole moments (Au):
   
7.        state          X           Y           Z        Dip. S.      Osc.
   
8.          1        -0.0000     -0.0000     -0.0000      0.0000      0.0000
   
9.          2        -0.0000     -0.0000     -0.0000      0.0000      0.0000
   
10.          3        -0.0000     -0.0000     -0.0000      0.0000      0.0000
    
11.          4        -0.0000     -0.0000     -0.0000      0.0000      0.0000
    
12.          5        -0.0000     -0.0000     -0.0000      0.0000      0.0000
    
13. Ground to excited state transition velocity dipole moments (Au):
    
14.        state          X           Y           Z        Dip. S.      Osc.
    
15.          1         0.0000      0.0000      0.0000      0.0000      0.0000
    
16.          2         0.0000      0.0000      0.0000      0.0000      0.0000
    
17.          3         0.0000      0.0000      0.0000      0.0000      0.0000
    
18.          4         0.0000      0.0000      0.0000      0.0000      0.0000
    
19.          5         0.0000      0.0000      0.0000      0.0000      0.0000
    
20. Ground to excited state transition magnetic dipole moments (Au):
    
21.        state          X           Y           Z
    
22.          1         0.0000      0.0000      0.0000
    
23.          2         0.0000      0.0000      0.0000
    
24.          3         0.0000      0.0000      0.0000
    
25.          4         0.0000      0.0000      0.0000
    
26.          5         0.0000      0.0000      0.0000
    
27. Ground to excited state transition velocity quadrupole moments (Au):
    
28.        state          XX          YY          ZZ          XY          XZ          YZ
    
29.          1         0.0000      0.0000      0.0000      0.0000      0.0000      0.0000
    
30.          2         0.0000      0.0000      0.0000      0.0000      0.0000      0.0000
    
31.          3         0.0000      0.0000      0.0000      0.0000      0.0000      0.0000
    
32.          4         0.0000      0.0000      0.0000      0.0000      0.0000      0.0000
    
33.          5         0.0000      0.0000      0.0000      0.0000      0.0000      0.0000
    
34. <0|del|b> * <b|rxdel|0> + <0|del|b> * <b|delr+rdel|0>
    
35. Rotatory Strengths (R) in cgs (10**-40 erg-esu-cm/Gauss)
    
36.        state          XX          YY          ZZ    R(velocity)    E-M Angle
    
37.          1         0.0000      0.0000      0.0000      0.0000       90.00
    
38.          2         0.0000      0.0000      0.0000      0.0000       90.00
    
39.          3         0.0000      0.0000      0.0000      0.0000       90.00
    
40.          4         0.0000      0.0000      0.0000      0.0000       90.00
    
41.          5         0.0000      0.0000      0.0000      0.0000       90.00
    
42. 1/2[<0|r|b>*<b|rxdel|0> + (<0|rxdel|b>*<b|r|0>)*]
    
43. Rotatory Strengths (R) in cgs (10**-40 erg-esu-cm/Gauss)
    
44.        state          XX          YY          ZZ     R(length)
    
45.          1         0.0000      0.0000      0.0000      0.0000
    
46.          2         0.0000      0.0000      0.0000      0.0000
    
47.          3         0.0000      0.0000      0.0000      0.0000
    
48.          4         0.0000      0.0000      0.0000      0.0000
    
49.          5         0.0000      0.0000      0.0000      0.0000
    
50.   1/2[<0|del|b>*<b|r|0> + (<0|r|b>*<b|del|0>)*] (Au)
    
51.        state          X           Y           Z        Dip. S.   Osc.(frdel)
    
52.          1        -0.0000     -0.0000     -0.0000      0.0000      0.0000
    
53.          2        -0.0000     -0.0000     -0.0000      0.0000      0.0000
    
54.          3        -0.0000     -0.0000     -0.0000      0.0000      0.0000
    
55.          4        -0.0000     -0.0000     -0.0000      0.0000      0.0000
    
56.          5        -0.0000     -0.0000     -0.0000      0.0000      0.0000
    
57. 
    
58. Excitation energies and oscillator strengths:
    
59. 
    
60. Excited State   1:      Triplet-A      8.9180 eV  139.03 nm  f=0.0000  <S**2>=2.000
    
61. This state for optimization and/or second-order correction.
    
62. Total Energy, E(TD-HF/TD-KS) =  -48406.4518915   
    
63. Copying the excited state density for this state as the 1-particle RhoCI density.
    
64. 
    
65. Excited State   2:      Triplet-A      9.0818 eV  136.52 nm  f=0.0000  <S**2>=2.000
    
66. 
    
67. Excited State   3:      Triplet-A      9.0940 eV  136.34 nm  f=0.0000  <S**2>=2.000
    
68. 
    
69. Excited State   4:      Triplet-A      9.0944 eV  136.33 nm  f=0.0000  <S**2>=2.000
    
70. 
    
71. Excited State   5:      Triplet-A      9.1354 eV  135.72 nm  f=0.0000  <S**2>=2.000
    
72.     1910 ->1912        0.33531
    
73.     1910 ->1921        0.10411
    
74.     1910 ->1944        0.13779
    
75.     1910 ->1996        0.10031
    
76. SavETr:  write IOETrn=   770 NScale= 10 NData=  16 NLR=1 NState=    5 LETran=     100.
    
77. Hyperfine terms turned off by default for NAtoms > 100.
    

复制代码

丁越

zjxitcc 发表于 2021-12-18 10:06
S0->T1是有跃迁信息的。是跃迁偶极矩为零，导致振子强度为零，但跃迁信息还是有的。你看第5个态也是三重 ...

呃呃，懂了。谢谢老师。是我看错楼主的表述了，还以为是为啥没有跃迁偶极矩信息

zjxitcc

丁越 发表于 2021-12-18 09:06
这是由于你算的是单重态到三重态的激发，在不考虑旋轨耦合的情况下是跃迁禁阻的，即跃迁偶极矩为零。你只能 ...

S0->T1是有跃迁信息的。是跃迁偶极矩为零，导致振子强度为零，但跃迁信息还是有的。你看第5个态也是三重态，就有跃迁信息。
估计是2L说的输出阈值的问题。

丁越

这是由于你算的是单重态到三重态的激发，在不考虑旋轨耦合的情况下是跃迁禁阻的，即跃迁偶极矩为零。你只能得到单重态向三重态激发的激发能，而得不到跃迁信息。

Novice

可能你前几个态任何一对轨道的贡献都小于Gaussian默认的0.1，所以没有输出。你在关键字一栏加上IOP(9/40=4)，这样输出文件就能把贡献大于0.0001的都显示出来了。