求助三重态(T2)一直未能收敛

Lym

使用关键字如下: #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

输出图像如下：

输入输出文件如下：

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

Lym

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

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

paramecium86

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

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

Lym

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

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

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

Lym

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呢？还是说我应该更换隐形溶剂模型，缩小步长先去得到一个优化成功的结构呢？
希望老师能指点我一下，麻烦您了

sobereva

上传压缩后的输出文件才好判断可能原因

但凡关键词涉及到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）

paramecium86

从那个Gradent出现尖峰的地方开始后面T2态能量都已经是负数了。这样继续优化没有什么用了。需要考虑势能面交叉的问题。