关于四点法算激发态重组能的问题

wILlD

根据下图，我的理解激发态重组能应该是E2-E3+E4-E1.E1是opt基态的能量；E2是基态的结构加td算出来的能量；E4的能量是基态结构td+opt后的体系的能量，E3是E4的结构再加td得出的体系能量。请问这个想法是正确的吗？

peptide

cokie 发表于 2024-12-25 17:10
以采用TD-DFT方法在gaussian中的计算结果为例
你可以在做TD OPT任务的输出文件中，同时读到E3和E4，分别 ...

豁然开朗的感觉，非常感谢老师指点！！！

cokie

peptide 发表于 2024-12-25 17:01
谢谢老师答疑，您回答中最后一句“基于优化后的激发态几何结构算单点能得到能量E3”，我还是有点不太懂， ...

以采用TD-DFT方法在gaussian中的计算结果为例
你可以在做TD OPT任务的输出文件中，同时读到E3和E4，分别位于输出文件的末尾，HF = xxxx（E3），以及优化最后一步所输出的“Total Energy, E(TD-HF/TD-DFT)”所在行的电子能量（E4）。
E3是基于优化后的激发态的几何结构，计算出的基态的电子能量（尚未弛豫到基态极小点or势能面全局最小点）。

peptide

cokie 发表于 2024-12-25 16:46
涉及激发态的重组能。
以S0-S1为例，优化基态几何结构并得到能量E1，基于优化后的基态结构算TD单点能得 ...

谢谢老师答疑，您回答中最后一句“基于优化后的激发态几何结构算单点能得到能量E3”，我还是有点不太懂，优化后的激发态结构算单点能，不还是激发态最小点E4吗，从激发态那个势能面垂直发射回到基态势能面的E3，再不需要TD了吗，我把我目前理解的关键词列在下面请老师帮看看哈：
E1: opt freq计算后，读取能量值；
E2: 以E1优化的结果，进行td计算，读取第一激发态能量值；
E4: 同样以E1优化的结果，进行td opt计算，即优化到S1极小点，读取能量值；
E3：以E4优化的结果，进行td计算，读取哪个值（这里就是我糊涂的地方，难道也是类似吸收过程，读取第一激发态能量值？如果按照老师说的，以E4优化的结果，算单点能，我糊涂的认为还是E4能量），
也可能是我对单点能计算概念有误解，是不是单点能计算过程中关键词可以设置是哪个态，麻烦老师能不能给我提示一下这个E3单点能关键词是啥，非常感谢！！！

cokie

peptide 发表于 2024-12-25 15:48
学习完本帖老师的讨论，再跟帖请教一下老师们一个问题，按照http://sobereva.com/330里的第一个例子，计算 ...

涉及激发态的重组能。
以S0-S1为例，优化基态几何结构并得到能量E1，基于优化后的基态结构算TD单点能得到能量E2，基于基态几何结构做激发态几何结构优化并得到能量E4，基于优化后的激发态几何结构算单点能得到能量E3

peptide

学习完本帖老师的讨论，再跟帖请教一下老师们一个问题，按照http://sobereva.com/330里的第一个例子，计算重组能是需要计算空穴重组能和电子重组能，即两个势能面分别是中性态和带电荷状态，计算输入文件不涉及TD计算，只涉及电荷和自旋的设置；那第二个例子是计算S0-T1重组能，则需要涉及TD算激发态，此时再不需要调整电荷了，请问老师我这样理解正确吗？
如果我想计算一个分子的重组能，用于与其发光情况建立定性关系，那我应该算涉及激发态的重组能，还是算涉及电荷改变状态的重组能呢，问题可能比较幼稚，还请老师不吝解惑，谢谢啦!

身在盘丝洞

cokie 发表于 2024-7-17 18:54
先不管你说的E3,E4和楼主说的一不一样，咱先统一按博文330里E3和E4的指代来进行后续对话：

330博文里E ...

E4是这个# opt b3lyp/6-311g(d,p) em=gd3bj   E3是# td=(singlets,nstates=10,root=1) b3lyp/6-311g(d,p) em=gd3bj他们理论级别是一致的吗？算出来像你说的那样“E3的单点能和TD+OPT得到的E4输出文件中最后的“HF=”后面的能量一致”

E1和E2适用与这个吗？

身在盘丝洞

cokie 发表于 2024-7-17 18:54
先不管你说的E3,E4和楼主说的一不一样，咱先统一按博文330里E3和E4的指代来进行后续对话：

330博文里E ...

此时如果计算级别严格统一的话，E3的单点能和TD+OPT得到的E4输出文件中最后的“HF=”后面的能量一致），这个是为什么？
我给了关键词了不就是为了说和楼主不一样吗
D3BJ这个是漏了
激发态结构呀，log文件最主要的不就是关键词吗

cokie

身在盘丝洞 发表于 2024-7-17 15:48
你不觉得楼主搞反了吗我说的E3和E4全是按照http://sobereva.com/330博文里来的。我给了关键词已经证明我 ...

先不管你说的E3,E4和楼主说的一不一样，咱先统一按博文330里E3和E4的指代来进行后续对话：

330博文里E4和E3都是基于S1的结构来算的，具体顺序是先TD+OPT得到E4，然后基于该结构DFT算单点得到E3（此时如果计算级别严格统一的话，E3的单点能和TD+OPT得到的E4输出文件中最后的“HF=”后面的能量一致），如果你是按照这个顺序算的，在保持计算级别/程序完全一致的前提下，无论如何也得不到E4能量低于E3的结果。

最后，你提问中上来就说E3和E4，本贴中一共就只有一张标了E3和E4是什么的图，不管楼主是否有指代错误，你这么问，回答你的人不肯定都按主贴图中的E3和E4来理解吗？你说的E3,E4是别的地方的，你好歹附上链接或者自己贴图；自己在另一个贴中问了同样的问题，但那个贴你的TD-B3LYP是加了D3BJ的，这个贴又不加；问问题的时候，也不说自己E3和E4分别是基于什么结构做的计算，也不放log文件，你这么个问法，怎么可能让人给你一击必中的回答到点上？他人的指代不清不是你指代不清的理由。

身在盘丝洞

cokie 发表于 2024-7-17 15:36
按楼主的E1-E4的定义，E4本来不就应该比E3低吗？

另外，在确保你优化结构用B3LYP不需要D3BJ的结果确定 ...

你不觉得楼主搞反了吗我说的E3和E4全是按照http://sobereva.com/330博文里来的。我给了关键词已经证明我和他的不一样了。楼主说的E4的能量是基态结构td+opt后的体系的能量，E3是E4的结构再加td得出的体系能量和他给的图都不一样，这句话本来就是错的，td完了之后再td是要干什么

cokie

身在盘丝洞 发表于 2024-7-17 15:25
老师，我E3和E4分别用的# b3lyp/6-311g(d,p)和# opt td=(singlets,nstates=10,root=1) b3lyp/6-311g(d,p) ...

按楼主的E1-E4的定义，E4本来不就应该比E3低吗？

另外，在确保你优化结构用B3LYP不需要D3BJ的结果确定可以用的前提下，算重组能时候你用TDDFT做完优化以后，直接在优化最后一步“Total Energy, E(TD-HF/TD-DFT) =”里就读出了E3的能量，然后在最后的“HF=”后面就读出E4能量了

身在盘丝洞

snljty2 发表于 2024-7-17 14:42
一般情况下可以，但是四点法计算重组能的时候不行，原因见为什么重组能计算的结构优化和单点能计算必须在 ...

老师，我E3和E4分别用的# b3lyp/6-311g(d,p)和# opt td=(singlets,nstates=10,root=1) b3lyp/6-311g(d,p)但是算完E4要比E3还小是什么原因

snljty2

身在盘丝洞 发表于 2024-7-17 14:23
老师，优化结构和计算能量用不同的泛函吗？

一般情况下可以，但是四点法计算重组能的时候不行，原因见为什么重组能计算的结构优化和单点能计算必须在同一级别？

身在盘丝洞

小范范1989 发表于 2021-11-1 14:45
看到你的问题，我感觉需要澄清两个点。
首先，优化S1的文件，S1的能量从哪里读取？可以是你说的“excite ...

老师，优化结构和计算能量用不同的泛函吗？

wILlD

小范范1989 发表于 2021-11-1 20:57
“E3是S1优化后的结构加td按照E2那么读取呢，还是直接就可以从s1优化的结构中读state 1的total energy.”
...

哈哈哈哈哈，确实。