求助，磷光发射能是负值

qingqing11

Stardust0831 发表于 2026-5-28 13:10
那基态大概率是自旋极化单重态了。
应该是那个闭壳层单重态的波函数不稳定，具体可以看stable=opt那一步 ...

Stability analysis using <AA,BB:AA,BB> singles matrix:
***********************************************************************
1PDM for each excited state written to RWF  633
Ground to excited state transition densities written to RWF  633

Eigenvectors of the stability matrix:

Excited state symmetry could not be determined.
Eigenvector   1:      Triplet-?Sym  Eigenvalue=-0.0695813  <S**2>=2.000
      93 -> 94         0.69520
The wavefunction has an RHF -> UHF instability.

Stability analysis using <AA,BB:AA,BB> singles matrix:
***********************************************************************
1PDM for each excited state written to RWF  633
Ground to excited state transition densities written to RWF  633

Eigenvectors of the stability matrix:

Excited state symmetry could not be determined.
Eigenvector   1:  1.122-?Sym  Eigenvalue= 0.0585066  <S**2>=0.065
     93B -> 94B        0.99129
The wavefunction is stable under the perturbations considered.
Leave Link  914 at Thu May 28 11:55:52 2026, MaxMem= 12884901888 cpu:            5785.2 elap:             241.6
(Enter /usr/local/cui/gaussian/gaussian16-a03-avx2/g16/l508.exe)
The wavefunction is already stable.
您好，上面提示已经从RKS到UKS，我再往下找就发现提示波函数稳定了，q想请问您这个可以判断这个结构的三重态结构下的基态是开壳层的吗

Stardust0831

qingqing11 发表于 2026-5-28 15:00
Stability analysis using  singles matrix:
****************************************************** ...

嗯，那S0是开壳层的。用stable=opt后得到的单重态能量当S0。

qingqing11

Stardust0831 发表于 2026-5-28 13:19
基于方法二，都带上stable=opt，重算一下试试看

没有明白您的意思，是先第一步：stable=opt 0 3,再第二步：opt 0 3,第三步：用第二步优化好的结构stable=opt 0 1,第四步：读第三步的波函数去算，最后第四步的能量减去第二步的吗

Stardust0831

qingqing11 发表于 2026-5-28 16:22
没有明白您的意思，是先第一步：stable=opt 0 3,再第二步：opt 0 3,第三步：用第二步优化好的结构stable= ...

基于T1结构，同时算S0和T1能量（也就是0 1和0 3直接两个单点），这两个单点都需要stable=opt。

qingqing11

Stardust0831 发表于 2026-5-28 16:25
基于T1结构，同时算S0和T1能量（也就是0 1和0 3直接两个单点），这两个单点都需要stable=opt。

好的好的，谢谢您，我试一下

qingqing11

Stardust0831 发表于 2026-5-28 16:25
基于T1结构，同时算S0和T1能量（也就是0 1和0 3直接两个单点），这两个单点都需要stable=opt。

0 3：-1013.686621；0 1：-1013.686514。最后相减是-0.067286825kcal/mol。您好，最后算出来还是负值

Stardust0831

qingqing11 发表于 2026-5-28 17:02
0 3：-1013.686621；0 1：-1013.686514。最后相减是-0.067286825kcal/mol。您好，最后算出来还是负值

其实0 3是基态也完全可能的。现在这个能量差其实算是非常接近了。
或许可以上一些多参考或者SF-TDDFT的方法细致看看。

qingqing11

Stardust0831 发表于 2026-5-28 18:42
其实0 3是基态也完全可能的。现在这个能量差其实算是非常接近了。
或许可以上一些多参考或者SF-TDDFT的 ...

您好，我用CAS计算了这个结构的磷光发射能，只有1kcal/mol,这个是不是不太合理，我调了好多个轨道，磷光发射能都很低，想请教一下这种情况合理吗?我应该怎么解决呢？补充：T1跟S0结构相比发生了断键，结构变化还是比较大的

Stardust0831

qingqing11 发表于 2026-6-2 16:47
您好，我用CAS计算了这个结构的磷光发射能，只有1kcal/mol,这个是不是不太合理，我调了好多个轨道，磷光 ...

正常的，如实报道即可。