daihui 发表于 2025-3-2 18:09
这样的话，跟E3算出来有区别吗

在S1结构下，计算TD单点能，那么你可以得到两个能量，其一是当前构型下的激发态能量（E3）和当前构型下的激发能（E3-E4）。

举个例子

Excited State   1:      Singlet-A      2.3163 eV  535.27 nm  f=0.5838  <S**2>=0.000
     142 ->143         0.70556
This state for optimization and/or second-order correction.
Total Energy, E(TD-HF/TD-DFT) =  -3361.64685560   
Copying the excited state density for this state as the 1-particle RhoCI density.

其中E3可以在“Total Energy, E(TD-HF/TD-DFT) =  -3361.64685560”此行读取
E4则为“(-3361.64685560) - (2.3163/27.2114)”

当前结构下的基态能量也可以在最后的HF = xxx读到，对一对就会发现，数值和“(-3361.64685560) - (2.3163/27.2114)”一致。

了解一下势能面的横纵坐标分别代表什么，就很好理解了。

daihui 发表于 2025-3-2 18:09
这样的话，跟E3算出来有区别吗

我想了一下。你算单点能用DFT时，不管是否是S1的结构，你只能得到基态的单点能，所以用b3lyp/6-31g(d) em=gd3bj算出来的应该不是E3
S1的结构下，用b3lyp/6-31g(d) em=gd3bj计算单点能应该算的是E4。
E3应该是用S1结构，在b3lyp/6-31g(d,p) td=(nstates=5,root=1) em=gd3bj下算单点能。

ljh123 发表于 2025-3-2 16:25
应该是在S1的结构下，用b3lyp/6-31g(d) em=gd3bj计算单点能

这样的话，跟E3算出来有区别吗

应该是在S1的结构下，用b3lyp/6-31g(d) em=gd3bj计算单点能