求助：S1-T1之间绝热能为负值

Lapheal

目的：计算得到S1-T1之间绝热能
思路：分别优化S1和T1的结构，分别读取FINAL SINGLE POINT ENERGY的能量
方法：inp中坐标一致，只改变mult分别为1和3，进行结构优化
!opt freq def2-TZVPP def2/J def2-QZVPP/C Normalprint Printbasis PrintMOs
!moread
%moinp "a.gbw"

%casscf
    nel  10
    norb 10
    mult 1
    nroots  10
    weights[0] = 1,0,0,0,0,0,0,0,0,0
end
问题：1，S1结构下的极小点能量<T1结构下的极小点能量（都无虚频），请问这是有可能发生的正常情况吗？因为不是在同一结构下得到的能量？
2，将优化得到的S1&T1的结构分别拿来做单点能计算在NEVPT2级别下，得到的能量依旧是S1态更低。

cokie

用casscf算是因为你的材料有多参考态的特征吗？如果是的话，如果你的材料基态不是单重态而是三重态，那你的T1（这个时候就应该理解为基态了）肯定是要比S1（激发态）低的。

Lapheal

cokie 发表于 2024-7-23 15:56
用casscf算是因为你的材料有多参考态的特征吗？如果是的话，如果你的材料基态不是单重态而是三重态，那你的 ...

是有机自由基体系，基态是单重态，吸光后被激发到S1态，通过ISC到达T1态。

cokie

Lapheal 发表于 2024-7-24 03:51
是有机自由基体系，基态是单重态，吸光后被激发到S1态，通过ISC到达T1态。

基态是单重态的有机自由基，那就是双自由基咯，如果你确定T1的存在是光激发到S1之后系间窜越得到的（毕竟双自由基的T态也可以靠热激活得到），那计算结果和实验结果对不上大概率是参考态选的不一样。

hebrewsnabla

你贴的这个输入算的是 S0

Lapheal

hebrewsnabla 发表于 2024-7-24 15:23
你贴的这个输入算的是 S0

请问优化S1的输入部分，mult应该怎么设置呢，该怎么区分S0和S1呢？他们都是单重态

hebrewsnabla

Lapheal 发表于 2024-7-24 15:28
请问优化S1的输入部分，mult应该怎么设置呢，该怎么区分S0和S1呢？他们都是单重态

weights[0] = 0,1,0,...

当基态是单重态时，S0 指的是最低的单重态（基态），S1 指的是第二低的单重态。