计算化学公社

标题: Gaussian 16 用TDDFT优化激发态结构收敛和断点续跑问题 [打印本页]

作者
Author: 莱万少夫斯基 时间: 2023-3-21 03:53
标题: Gaussian 16 用TDDFT优化激发态结构收敛和断点续跑问题
本帖最后由莱万少夫斯基于 2023-3-22 07:40 编辑

大家好。我正在使用Gaussian 16通过TDDFT优化激发态结构，模拟UV-Vis。收敛问题。我发现在out输出文件中，似乎已经显示在某一个interation全部收敛，然而从下一个interation开始又变成全部不收敛了。我比较困惑这种情况是否正常。另外，我发现使用现在的out文件，Multiwfn已经可以模拟出UV-Vis图谱。我不知道是否可以停止优化？#P opt=restart B3LYP/gen pseudo=read TD(nstates=203,restart) guess=read。希望可以获得大佬的指点！不胜感激！

作者
Author: sobereva 时间: 2023-3-21 04:19
你犯了一个关键性错误，就是你当前优化的是第一激发态，你计算最低3~5个激发态就足够了，算203个态完全是无意义地严重白浪费时间。而且模拟UV-Vis并不需要优化激发态。算发射光谱才需要优化激发态，而此时明显又用不着算200个态。
仔细把此文看了，该说的基本知识都说了：
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314（http://bbs.keinsci.com/thread-2413-1-1.html）

作者
Author: 莱万少夫斯基 时间: 2023-3-21 05:28
本帖最后由莱万少夫斯基于 2023-3-21 06:00 编辑

sobereva 发表于 2023-3-21 04:19
你犯了一个关键性错误，就是你当前优化的是第一激发态，你计算最低3~5个激发态就足够了，算203个态完全是无 ...

谢谢社长回复。可是我是按照这个体系的相关文献来的。他们模拟了S0到S200的TDDFT去获得UV-Vis光谱。而且，在此之前我优化过类似最低的5个态，那个时候获得的光谱区间离UV-Vis相当远。当我增加到200个态的时候，就已经和文献报道的光谱区间比较吻合了。我看了那个贴文，不知道我理解的对不对。我们需要计算垂直吸收能，要保持能量基态最低结构。如果要计算到S200的结构不优化激发态的话，正确的设置是否是#P B3LYP/gen pseudo=read TD(nstates=203, root=203) guess=read？请您明示，谢谢。

作者
Author: sobereva 时间: 2023-3-21 06:19

莱万少夫斯基发表于 2023-3-21 05:28
谢谢社长回复。可是我是按照这个体系的相关文献来的。他们模拟了S0到S200的TDDFT去获得UV-Vis光谱。而且 ...

完全错误，基本概念和逻辑还没搞明白。没好好看我的文章。

作者
Author: 莱万少夫斯基 时间: 2023-3-21 07:07

sobereva 发表于 2023-3-21 06:19
完全错误，基本概念和逻辑还没搞明白。没好好看我的文章。

感谢您的回复。按照您4.1的例子，我使用基态结构
#P B3LYP/gen pseudo=read TD(nstates=203,root=200)，这样的设置是否可以正确的用于考察S0到S200的垂直吸收能，并用于UV-Vis的模拟？

作者
Author: KSeGaSn 时间: 2023-3-21 13:06
请问楼主方便提供下相关的文章吗？这吸收光谱对应到S0→S200的激发我还真是挺好奇是什么样的体系… 平常如果我算的低激发态这么红的不可接受我就换泛函了，从来没想过（也没条件）一次算203个激发能的……

作者
Author: wzkchem5 时间: 2023-3-21 16:56

莱万少夫斯基发表于 2023-3-21 00:07
感谢您的回复。按照您4.1的例子，我使用基态结构
#P B3LYP/gen pseudo=read TD(nstates=203,root=200)， ...

可以，这样会算203个激发态，其中除了最高的几个激发态可能例外以外，第n个激发态就是能量第n低的激发态。虽然不能保证前200个激发态都严格满足”第n个激发态就是能量第n低的激发态“，但是这一点引入的不确定性对于画UV-Vis图而言是可以接受的。
另外注意此时无需写root，但写了也不影响结果的可用性

作者
Author: 莱万少夫斯基 时间: 2023-3-21 19:46

wzkchem5 发表于 2023-3-21 16:56
可以，这样会算203个激发态，其中除了最高的几个激发态可能例外以外，第n个激发态就是能量第n低的激发态 ...

感谢您的回复！这样是不是还得考虑三重激发态的可能，然后把三重激发态和单重态一起对比才能说谁是能量最低的n个激发态呢？

作者
Author: wzkchem5 时间: 2023-3-21 20:16

莱万少夫斯基发表于 2023-3-21 12:46
感谢您的回复！这样是不是还得考虑三重激发态的可能，然后把三重激发态和单重态一起对比才能说谁是能量最 ...

我说得不准确，我的意思是”第n个激发态就是能量第n低的单重态激发态“，也就是你的计算即使连这个也不能保证，但是这一点对计算UV-Vis的影响可以忽略

欢迎光临计算化学公社 (http://bbs.keinsci.com/)