请教对于激发态能量的计算

lmyqwq

各位老师好，学生正在尝试绘制如图1所示的能级图，用于分析分子在不同激发态下的能量。参考了Sob老师《Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法》（http://sobereva.com/314），以下是学生的计算步骤，恳请各位老师指点是否正确。
1.基态优化:优化得到S0的极小点结构，读取SCF Done后的能量，记录其能量为 E1。将E1定义为相对能量零点（0 eV）。
2.垂直吸收能：在S0最优构型下进行TDDFT单点计算，从输出文件中读取第一激发态（S1）的垂直激发能，即图1中S0构型处S0与S1的垂直能量差。
3.激发态的优化（S1）：对S1进行结构优化，得到S1的极小点结构。从优化输出文件中读取Total Energy后的能量，作为S1态的总能量，记录为E3，绝热激发能=E3-E1。
4.垂直发射能：在S1最优构型下进行单点能计算，读取SCF Done后面的能量，记录为E4，垂直发射能=E3-E4,对应图中的S1最优构型下的S1与S0的垂直能量差。
能级图如图2展示的这样。学生的问题是：（1）以上对S0/S1的处理方法是否正确？（2）若正确，该方法是否可以推广用于计算T1（三重态）作为最优构型下的相关能量？
期待各位老师的批评指正，非常感谢！

Kevin_Chen

看上去是没什么问题，但是图1中的内容不叫Electron-hole analysis，这张图是不同几何结构下的激发能级的示意图。此外我看你左上角的那个正负电荷分离好像是有电荷转移，电荷转移不要用B3LYP算，用CAM-B3LYP之类的范围分离泛函会好很多，因为电荷转移激发需要长程高HF成分保持正确的渐近行为，B3LYP 20%的HF成分对于长程太低了

sobereva

“SOC老师” 是什么鬼
提到博文时要主动贴链接，给你补上了

“对优化得到的结果进行激发” 描述得完全语义不明。如果你要说对S0极小点结构用TDDFT方法计算激发态，就要明确这么写清楚。“S0作为最优构型”这个描述也很莫名其妙，大抵你是想说“使用S0极小点结构作为当前TDDFT计算用的结构”。对S1的描述也是非常迷。描述问题时务必要用专业术语准确、没有歧义地表达。
当前你的描述过于混乱，没法确切回答做法是否正确。

lmyqwq

Kevin_Chen 发表于 2026-8-23 00:15
看上去是没什么问题，但是图1中的内容不叫Electron-hole analysis，这张图是不同几何结构下的激发能级的示 ...

嗯嗯感谢老师，我后续计算电荷转移现象比较明显的分子是使用您提到的这个泛函，图1下面的那一段是我截图没截好没描述清楚，抱歉。

lmyqwq

sobereva 发表于 2026-8-23 02:30
“SOC老师” 是什么鬼
提到博文时要主动贴链接，给你补上了

非常抱歉老师，昨天晚上发帖的时候马虎了。我重新编辑了一下帖子，修改了一下之前的描述内容，以为我再发帖的时候仔细检查一下，给您添麻烦了。

fatpig

看起来是正确。但第二问：若正确，该方法是否可以推广用于计算T1（三重态）作为最优构型下的相关能量...

有没有一种可能，是先有T1的计算，再推广到S1的计算...