计算化学公社

标题: OPT结果扭曲 [打印本页]

作者
Author: 小书童 时间: 2016-1-12 16:29
标题: OPT结果扭曲
本帖最后由小书童于 2016-1-12 18:20 编辑

用高斯计算结构的OPT+freq，计算正常结束，频率分析四个YES，最终的结果扭曲了。我打开优化过程的步数看，发现优化过程中出现了我要的图像，但是不是最终结果。计算文件如附件。
问题1：优化过程很平坦，我记得老师说过怎么处理这类问题，但是忘记了，希望老师讲解下。
问题2：我觉得我想要的结构和最总的结构在能量和力上都差不多啊，为什么最终得到了个扭曲的结果呢？应该怎么处理，
问题3：我扭曲的结果对吗？下一步计算用这个扭曲的结果，还是我要重新计算，如果重新计算，应该如何处理呢？
先感谢老师的指点，

！

作者
Author: bomsaude 时间: 2016-1-12 19:31
你这个分子电荷是0吗？因为去了一个H。

作者
Author: 小书童 时间: 2016-1-12 20:13

bomsaude 发表于 2016-1-12 19:31
你这个分子电荷是0吗？因为去了一个H。

是0，这个三个活性位的碳基模型，然后HCN和NO分别吸附到碳基上，11C和10C以及25C上没有H用来模拟活性位的。

作者
Author: sobereva 时间: 2016-1-12 21:10
没有虚频，说明在当前计算级别下这就是合理的极小点结构。你也可以再试试其它级别，比如M062X/TZVP，看是否都是这样。
这种结构极小点结构不是平面的是正常的事，因为有五元环和六元杂环存在，由于一些张力因素可以导致体系稍微变形成碗状。

作者
Author: 小书童 时间: 2016-1-12 21:42

sobereva 发表于 2016-1-12 21:10
没有虚频，说明在当前计算级别下这就是合理的极小点结构。你也可以再试试其它级别，比如M062X/TZVP，看是否 ...

那么您的意思是可能这个结构就是变形的。还有，结构优化一直平坦的走很久的这种我记得您说过怎么可以快点来着，您在指教一下。

作者
Author: sobereva 时间: 2016-1-12 22:07

小书童发表于 2016-1-12 21:42
那么您的意思是可能这个结构就是变形的。还有，结构优化一直平坦的走很久的这种我记得您说过怎么可以快点 ...

是
用gdiis，可以适当改用loose收敛限

作者
Author: 小书童 时间: 2016-1-13 11:27

sobereva 发表于 2016-1-12 22:07
是
用gdiis，可以适当改用loose收敛限

老师按照您的方法，从之前优化的步数中选择一个比较合适的点（没扭曲的），然后用gdiis方式进行OPT，但是结果都有一个问题，那就是会出现一个虚频。这个时候应该如何解决呢？

作者
Author: sobereva 时间: 2016-1-13 17:16

小书童发表于 2016-1-13 11:27
老师按照您的方法，从之前优化的步数中选择一个比较合适的点（没扭曲的），然后用gdiis方式进行OPT，但是 ...

所以没优化到极小点。
如果得到的是平面结构，有虚频那是必然的，因为当前级别下极小点是扭曲的。

作者
Author: 小书童 时间: 2016-1-13 20:04

sobereva 发表于 2016-1-13 17:16
所以没优化到极小点。
如果得到的是平面结构，有虚频那是必然的，因为当前级别下极小点是扭曲的。

恩恩，了然了，谢谢老师！

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