cokie 发表于 2024-12-25 17:10 豁然开朗的感觉，非常感谢老师指点！！！ |
本帖最后由 cokie 于 2024-12-25 17:12 编辑 peptide 发表于 2024-12-25 17:01 以采用TD-DFT方法在gaussian中的计算结果为例 你可以在做TD OPT任务的输出文件中，同时读到E3和E4，分别位于输出文件的末尾，HF = xxxx（E3），以及优化最后一步所输出的“Total Energy, E(TD-HF/TD-DFT)”所在行的电子能量（E4）。 E3是基于优化后的激发态的几何结构，计算出的基态的电子能量（尚未弛豫到基态极小点or势能面全局最小点）。 |
| 参与人数Participants 1 | eV +4 | 收起 理由Reason |
|---|---|---|
|
| + 4 | 谢谢 |
本帖最后由 peptide 于 2024-12-25 17:03 编辑 cokie 发表于 2024-12-25 16:46 谢谢老师答疑，您回答中最后一句“基于优化后的激发态几何结构算单点能得到能量E3”，我还是有点不太懂，优化后的激发态结构算单点能，不还是激发态最小点E4吗，从激发态那个势能面垂直发射回到基态势能面的E3，再不需要TD了吗，我把我目前理解的关键词列在下面请老师帮看看哈： E1: opt freq计算后，读取能量值； E2: 以E1优化的结果，进行td计算，读取第一激发态能量值； E4: 同样以E1优化的结果，进行td opt计算，即优化到S1极小点，读取能量值； E3：以E4优化的结果，进行td计算，读取哪个值（这里就是我糊涂的地方，难道也是类似吸收过程，读取第一激发态能量值？如果按照老师说的，以E4优化的结果，算单点能，我糊涂的认为还是E4能量）， 也可能是我对单点能计算概念有误解，是不是单点能计算过程中关键词可以设置是哪个态，麻烦老师能不能给我提示一下这个E3单点能关键词是啥，非常感谢！！！ |
peptide 发表于 2024-12-25 15:48 涉及激发态的重组能。 以S0-S1为例，优化基态几何结构并得到能量E1，基于优化后的基态结构算TD单点能得到能量E2，基于基态几何结构做激发态几何结构优化并得到能量E4，基于优化后的激发态几何结构算单点能得到能量E3 |
|
本帖最后由 peptide 于 2024-12-25 15:52 编辑 学习完本帖老师的讨论，再跟帖请教一下老师们一个问题，按照http://sobereva.com/330里的第一个例子，计算重组能是需要计算空穴重组能和电子重组能，即两个势能面分别是中性态和带电荷状态，计算输入文件不涉及TD计算，只涉及电荷和自旋的设置；那第二个例子是计算S0-T1重组能，则需要涉及TD算激发态，此时再不需要调整电荷了，请问老师我这样理解正确吗？ 如果我想计算一个分子的重组能，用于与其发光情况建立定性关系，那我应该算涉及激发态的重组能，还是算涉及电荷改变状态的重组能呢，问题可能比较幼稚，还请老师不吝解惑，谢谢啦! 
|
本帖最后由 身在盘丝洞 于 2024-7-17 22:20 编辑 cokie 发表于 2024-7-17 18:54 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适用与这个吗？ |
本帖最后由 身在盘丝洞 于 2024-7-17 22:04 编辑 cokie 发表于 2024-7-17 18:54 此时如果计算级别严格统一的话，E3的单点能和TD+OPT得到的E4输出文件中最后的“HF=”后面的能量一致），这个是为什么？ 我给了关键词了不就是为了说和楼主不一样吗 D3BJ这个是漏了 激发态结构呀，log文件最主要的不就是关键词吗 |
本帖最后由 cokie 于 2024-7-17 18:59 编辑 身在盘丝洞 发表于 2024-7-17 15:48 先不管你说的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文件，你这么个问法，怎么可能让人给你一击必中的回答到点上？他人的指代不清不是你指代不清的理由。 |
本帖最后由 身在盘丝洞 于 2024-7-17 15:50 编辑 cokie 发表于 2024-7-17 15:36 你不觉得楼主搞反了吗我说的E3和E4全是按照http://sobereva.com/330博文里来的。我给了关键词已经证明我和他的不一样了。楼主说的E4的能量是基态结构td+opt后的体系的能量，E3是E4的结构再加td得出的体系能量和他给的图都不一样，这句话本来就是错的，td完了之后再td是要干什么 |
本帖最后由 cokie 于 2024-7-17 15:39 编辑 身在盘丝洞 发表于 2024-7-17 15:25 按楼主的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还小是什么原因 |
身在盘丝洞 发表于 2024-7-17 14:23 一般情况下可以，但是四点法计算重组能的时候不行，原因见为什么重组能计算的结构优化和单点能计算必须在同一级别？ |
小范范1989 发表于 2021-11-1 14:45 老师，优化结构和计算能量用不同的泛函吗？ |
小范范1989 发表于 2021-11-1 20:57 哈哈哈哈哈，确实。 |
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2025-8-12 21:13 , Processed in 0.258470 second(s), 26 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites