计算化学公社

标题: 过渡态TS优化报错l9999求助，往上找具体错误原因时并未找到opt stopped [打印本页]

作者
Author: swyrzsyf 时间: 2023-12-18 17:38
标题: 过渡态TS优化报错l9999求助，往上找具体错误原因时并未找到opt stopped
(, 下载次数 Times of downloads: 9)
图为初猜过渡态结构（左）和产物（右）
刚入门，在优化过渡态的时候出现了很多种类的错误，都在公社中一一找到并解决了，但每次都有新的问题出现。
经常优化的时候我发现有些键会莫名其妙断掉。
这次的问题就是out文件报错l9999，并且向上搜索opt stopped并没有搜索到，只找到一个opt aborted，也看了优化曲线，没有震荡的情况。

我想知道是我初猜的有问题吗？我这个反应就是想计算氧自由基进攻碳碳双键形成环氧的一个过程。
所用的基组为#p opt=(calcfc,ts,maxcyc=500,noeigen,cartesian) freq ub3lyp/6-31g(d,p)，电荷和自旋多重度为0 3。
gjf文件和报错的out文件如下，麻烦老师帮我看看，谢谢。
(, 下载次数 Times of downloads: 22)

作者
Author: Yaocy 时间: 2023-12-18 20:18

Item Value Threshold Converged?
Maximum Force 0.024277 0.000450 NO
RMS Force 0.007693 0.000300 NO
Maximum Displacement 0.000000 0.001800 YES
RMS Displacement 0.000000 0.001200 YES
Predicted change in Energy= 0.000000D+00
Optimization aborted.
-- No acceptable step.

复制代码

输出文件中有这一段，看来几何优化没有收敛，而且对应的 -DE/DX 也有很多星号

作者
Author: sobereva 时间: 2023-12-18 21:03
过渡态初猜太糟糕
而且应当当单重态来算
没必要写cartesian

作者
Author: swyrzsyf 时间: 2023-12-19 08:46

sobereva 发表于 2023-12-18 21:03
过渡态初猜太糟糕
而且应当当单重态来算
没必要写cartesian

那请问初猜应该是怎样的呢？
三重态是因为我觉得有个氧自由基所以用三重态算的
好的，我把cartesian删掉试试

作者
Author: swyrzsyf 时间: 2023-12-19 08:47

Yaocy 发表于 2023-12-18 20:18
输出文件中有这一段，看来几何优化没有收敛，而且对应的 -DE/DX 也有很多星号

那您知道这种情况应该怎样去尝试解决吗？或者说朝哪个方向去修改试错呢？

作者
Author: sobereva 时间: 2023-12-19 12:30

swyrzsyf 发表于 2023-12-19 08:46
那请问初猜应该是怎样的呢？
三重态是因为我觉得有个氧自由基所以用三重态算的
好的，我把cartesian删 ...

不牵扯反应的区域，显然不该有奇怪的键长、诡异的结构扭曲。结合基本结构化学常识判断
并且先把产物优化出来，基于产物的结构去调整结构得到过渡态初猜，这样容易得多
单个氧是三重态，绝对不代表过渡态时候是三重态。单个氧是三重态是因为2p原子轨道简并，然而在过渡态时氧已经部分成键了，简并性显然会破坏

作者
Author: swyrzsyf 时间: 2023-12-19 14:11

sobereva 发表于 2023-12-19 12:30
不牵扯反应的区域，显然不该有奇怪的键长、诡异的结构扭曲。结合基本结构化学常识判断
并且先把产物优化 ...

好的，我明白了，我去再试试，谢谢老师。

作者
Author: swyrzsyf 时间: 2023-12-20 09:36
本帖最后由 swyrzsyf 于 2023-12-20 09:42 编辑

swyrzsyf 发表于 2023-12-19 14:11
好的，我明白了，我去再试试，谢谢老师。

成功计算出来了

作者
Author: swyrzsyf 时间: 2023-12-20 09:38

sobereva 发表于 2023-12-19 12:30
不牵扯反应的区域，显然不该有奇怪的键长、诡异的结构扭曲。结合基本结构化学常识判断
并且先把产物优化 ...

接过渡态优化问题，上次经过初猜修正后重新进行TS计算，这次能够成功计算出输出文件了，在这里谢谢老师的指导，点一下就把我自个研究一个周的困惑解决了。
但是现在出现了一个新的问题。
最终得到的结果只有一个虚频，但是该虚频振动并不对应我的反应方向，而是甲基自己在旋转。
如下图左为我初猜结构，图右是过渡态优化得到的最终结构（这对应上了我的目标生成物......）

(, 下载次数 Times of downloads: 10)
想请教下虚频不是我要的虚频这种情况是什么原因产生的呢？应该怎么解决呢？
输入文件和输出文件如下
(, 下载次数 Times of downloads: 4)

作者
Author: sobereva 时间: 2023-12-21 03:13

swyrzsyf 发表于 2023-12-20 09:38
接过渡态优化问题，上次经过初猜修正后重新进行TS计算，这次能够成功计算出输出文件了，在这里谢谢老师的 ...

初猜结构不当，优化到了无意义的过渡态

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