ORCA 中 计算 绝热△EST 问题 求教

waker

近期在做关于绝热△Est的计算，计算方法如下：
（1） 用TD-DFT优化分子的S1态结构，然后计算激发态单点能；
（2） 用UKS 法优化分子的T1态结构，然后用TD-DFT法计算激发态单点能；
（3）用上述（1）值-（2）值
结果计算出来的△Est 是 负数  为  - 0.065 eV（S1能量 < T1 能量）。T1 态结构下计算出来的△Est = 0.056 eV


各位大佬，想问一下是哪里出了问题吗？

waker

wzkchem5 发表于 2023-5-17 16:47
应该用TDDFT优化T1结构，否则结构优化的理论级别不统一

知道了，谢谢老师

wzkchem5

应该用TDDFT优化T1结构，否则结构优化的理论级别不统一

waker

我又补算了个T2，感觉还是不对的样子。! PBE0 D3 def2-TZVP(-f) def2/J RIJCOSX tightSCF noautostart miniprint nopop
%maxcore  8000
%pal nprocs   18 end
%tddft
nroots 5
IROOT 2
IROOTMULT Triplet
end
* xyzFILE   0   1 T2.xyz

结果：
CIS/TD-DFT TOTAL ENERGY
-----------------------

    E(SCF)  =  -2129.896007770 Eh
    DE(CIS) =      0.102586415 Eh (Triplet root  2)
    ----------------------------- ---------
    E(tot)  =  -2129.793421355 Eh



-------------------------------------------------------------------------------
                          DFT DISPERSION CORRECTION                           

                              DFTD3 V3.1  Rev 1                                
                          USING Becke-Johnson damping                          
-------------------------------------------------------------------------------
-------------------------   ----------------
Dispersion correction           -0.124759586
-------------------------   ----------------


-------------------------   --------------------
FINAL SINGLE POINT ENERGY     -2129.918180941001
-------------------------   --------------------

waker

S1，T1的单点能计算输出文件如附件