zjxitcc 发表于 2023-3-24 12:50 谢谢老师 |
|
常规做法：读取该稳定波函数，做TDDFT单点计算。由于基态为UHF/UDFT，此时TD(triplet)括号中的triplet无效，故而不用写，你无法控制程序仅算单重态激发态 亦或 仅算三重态激发态，只能从算出的结果里去找哪个是T1态。找到哪个是T1态后，再读取该稳定波函数，指定目标态，做TDDFT结构优化计算。 上述过程中，你可能会发现在TDDFT结果里，根据自旋多重度很难分辨哪个是T1态，这是常规TDDFT原理缺陷造成的，需要换用更严谨的方法，例如单参考方法中有Spin Flip TDDFT，多组态/多参考态方法中有CASSCF, CASPT2, NEVPT2等等。 |
|
“稳定波函数的<S**2>”是从SCF Done附近读的（也就是和正常SCF计算的读法一样），而不是从稳定性分析的输出里读的 |
|
老铁，不要随便贴个图。你贴的图里 既无 稳定波函数的电子能量，也无 稳定波函数的<S**2>。 |
bms 发表于 2023-3-23 19:30 先说说稳定的波函数 电子能量是多少，<S**2>又是多少。另外，标题中把Gaussian拼错了。 |
zjxitcc 发表于 2023-3-23 13:13 老师，我通过对优化的T1结构用了stable=opt得到了稳定的波函数，我现在需要怎么去计算T1态的的能量呢？再对稳定波函数的T1结构进行TDDFT的话得到的是不是就不是T1能量了？ |
本帖最后由 zjxitcc 于 2023-3-23 13:15 编辑 bms 发表于 2023-3-23 12:51 你把问题弄复杂了，你的计算过程在这个问题里是没有关系的，因此你不用告诉大家你是怎么算的，也不用说结构是啥样的（当然，你说出计算过程，有助于判断计算合理性，但与此贴讨论的负激发能几乎没有关系，属于另一个问题）。 stable或stable=opt 必须在出现负激发能的结构下进行（我们不管你这个结构是怎么来的）、且必须读取 负激发能的结构的波函数，你现在似乎不满足这两点，因此要先正确地检验波函数稳定性，回到本帖主题上来，以免讨论走偏。 |
wzkchem5 发表于 2023-3-22 21:31 老师，我稳定性计算的是基态结构优化后的结构，然后取得基态结构优化后的log文件去计算的稳定性，我计算S0是用优化过的基态结构做TDDFT计算出基态S0能量，得到uv-vis曲线，然后还是通过对优化过的基态结构做opt+TDDFT 0 1 single 和0 1 triplet计算S1 和T1 结构和能量，我想请问下老师我这么是对的吗？目前只有T1出现了负的能量。另外都是用的B3LYP/6-31G* 另外附上分子式 |
bms 发表于 2023-3-22 14:27 稳定性计算得读你TDDFT计算的波函数，而且还是得用你怀疑可能有不稳定问题的结构，在你的例子里也就是激发能变成负数的结构。波函数和结构两者都必须正确。否则你不是在检验TDDFT波函数的稳定性，而是重新产生了一个新的波函数然后检验它的稳定性，不能用来证明你TDDFT计算里的波函数是稳定的。 |
zjxitcc 发表于 2023-3-22 16:13 上传在楼上了 麻烦您看一下 |
| 好的老师 麻烦您看一下 |
43.56 KB, 下载次数 Times of downloads: 13
手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图
GMT+8, 2026-2-18 23:05 , Processed in 0.175664 second(s), 25 queries , Gzip On.
组图打开中，请稍候......
收藏 Add to favorites