利用态指定方式考虑溶剂的TDDFT计算求助

让你变成回忆

利用态指定方式，在TDDFT级别下考虑溶剂效应，我的计算过程如下：（1）在溶剂下做基态单点计算，把溶剂信息写到chk文件里（noneq=write写在坐标后空一行的位置）；
（2）在TDDFT级别下，scrf=(externaliteration, read, solvent=XXX)  noneq=read（写在坐标后空一行的位置）；

正常计算结束以后，有如下问题：
（1）按照sob老师的帖子，激发态能应该为：E(After PCM corrections) - E(SCF Done)，前者为上面第二步的输出，后者为上面第一步的输出。
（2）按照态指定的计算流程，考虑的S1态与溶剂的迭代到自洽以后，为什么最后一步的S1态能量与上面（1）计算出来的不一致？
（3）如果按照SS方式考虑溶剂，那么要对该激发态进行分析，比如NTO、电子-空穴分析等，应该如何进行？是直接取（2）中的态吗？如果是，为什么有上面（2）的问题？

sobereva

2 看不懂。(1)是基态计算，S1态能量怎么和基态能量能一致
3 是

计算步骤的编号不要和问题的编号相同，否则容易混淆

让你变成回忆

sobereva 发表于 2019-6-21 11:45
2 看不懂。(1)是基态计算，S1态能量怎么和基态能量能一致
3 是

sob老师您好，我从新整理表述如下：
1. 在考虑SS下，S1态激发态的能量应该为：E(After PCM corrections) - E(SCF Done)，其中E(After PCM corrections)从上面的（2）步骤中读取，E(SCF Done)从上面的（1）步骤中读取，这是从sob老师的帖子中的得到的；
2. 但是在上面的步骤（2）中，我从迭代的最后一步读取S1态能量，却和1计算的不等，这是为什么？
3. 做NTO分析等等，就算上面2中读取的态进行分析。

sobereva

让你变成回忆 发表于 2019-6-21 13:36
sob老师您好，我从新整理表述如下：
1. 在考虑SS下，S1态激发态的能量应该为：E(After PCM corrections) ...

2 如果你是指直接读取的Excited State   1:      Singlet-A"     4.3060 eV  287.93 nm这种地方的数据，肯定和正规方式算的SS结果不同，这相当于当前反应场下线性响应方式给的结果
3 没法严格地在SS下做电子激发分析（除了那种基于激发态密度与基态密度差的分析）

让你变成回忆

sobereva 发表于 2019-6-21 23:16
2 如果你是指直接读取的Excited State   1:      Singlet-A"     4.3060 eV  287.93 nm这种地方的数据， ...

好的，明白了。 谢谢sob老师。