计算化学公社

标题: 求助三重态(T2)一直未能收敛 [打印本页]

作者
Author: Lym 时间: 2022-7-8 09:05
标题: 求助三重态(T2)一直未能收敛
本帖最后由 Lym 于 2022-7-8 14:23 编辑

使用关键字如下: #p opt(MaxStep=5,NoTrust,GDIIS) freq td=(root=2,triplet)/6-311g(d,p) s crf=(smd,solvent=generic,read) em=gd3 pbe1pbe
但是目前一直是，后两项不能收敛
Maximum Force          0.000158    0.000450    YES
RMS    Force          0.000023    0.000300    YES
Maximum Displacement    0.029325    0.001800    NO
RMS    Displacement    0.005810    0.001200    NO

输出图像如下：
(, 下载次数 Times of downloads: 11) (, 下载次数 Times of downloads: 6)
输入输出文件如下：

希望能得到老师们的指点，麻烦啦

作者
Author: paramecium86 时间: 2022-7-8 10:35
从那个Gradent出现尖峰的地方开始后面T2态能量都已经是负数了。这样继续优化没有什么用了。需要考虑势能面交叉的问题。

作者
Author: sobereva 时间: 2022-7-9 04:07
上传压缩后的输出文件才好判断可能原因

但凡关键词涉及到read坐标后的内容的时候，后面的内容必须在提问时一并交代出来。

不建议优化时候用SMD，下文说了
谈谈隐式溶剂模型下溶解自由能和体系自由能的计算
http://sobereva.com/327（http://bbs.keinsci.com/thread-3345-1-1.html）

另注意看此文
谈谈势能面交叉对激发态优化的影响
http://sobereva.com/468（http://bbs.keinsci.com/thread-12318-1-1.html）

作者
Author: Lym 时间: 2022-7-9 18:15

paramecium86 发表于 2022-7-8 10:35
从那个Gradent出现尖峰的地方开始后面T2态能量都已经是负数了。这样继续优化没有什么用了。需要考虑势能面 ...

嗯嗯，谢谢老师的回复。我当前想要解释的问题是S1与T2发生系间穿越，导致无辐射反应的发生。我已经分别优化了此分子的S0,S1,T1的结构，但是在进一步优化T2的时候就发现一直未能得到其稳定的结构。我已经看过sob老师推荐的文章（关于势能面交叉对于激发态优化的影响）,但是我仍然还是有一些不太理解
1. 透热态绝热态的区别是想要的那个激发态和最后得到能量的顺序如果一直是透热态否则是绝热态嘛？也就是我说我当前优化T2，在交叉点右侧能量excited energy 2<excited energy 1则表明分子走到了绝热态了嘛？
2. 老师所说的继续优化已经没有什么用，是因为势能面发生交叉后，原本想得到T2的结构，现在已经跑到了T1上？我本来是想要得到T2的结构后，去讨论S1->T2之间的ISC，但是我现在有点迷茫不知道下一步应该是优化交叉点？用S1和交叉点去讨论ISC呢？还是说我应该更换隐形溶剂模型，缩小步长先去得到一个优化成功的结构呢？
希望老师能指点我一下，麻烦您了

作者
Author: Lym 时间: 2022-7-9 22:32
本帖最后由 Lym 于 2022-8-2 20:43 编辑

sobereva 发表于 2022-7-9 04:07
上传压缩后的输出文件才好判断可能原因

但凡关键词涉及到read坐标后的内容的时候，后面的内容必须在提问 ...

sob老师，这是我的输入输出文件
read后面是输入的混合溶剂参数，您今天上一天的课应该会比较累了，等老师不忙的时候麻烦再帮我看一眼了。谢谢老师。

作者
Author: paramecium86 时间: 2022-7-9 23:07

Lym 发表于 2022-7-9 18:15
嗯嗯，谢谢老师的回复。我当前想要解释的问题是S1与T2发生系间穿越，导致无辐射反应的发生。我已经分别优 ...

不仅仅是T1 T2势能面交叉的问题。而是如果激发能是负的说明你用来当做参考态的“基态”已经不再是基态。那这样，TD算法本身是失效的。

作者
Author: Lym 时间: 2022-7-10 09:40

paramecium86 发表于 2022-7-9 23:07
不仅仅是T1 T2势能面交叉的问题。而是如果激发能是负的说明你用来当做参考态的“基态”已经不再是基态。 ...

嗯嗯，老师那我去检查一下波函数的稳定性，然后在稳定波函数上继续进行opt td(root=2,triplet，nastates=5) 改为iefpcm溶剂在去优化一下，去得到T2的结构，之后如果得到T2的结构之后再去改为SMD溶剂模型对嘛？

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