计算化学公社

标题: 只是改变了基组，方法没变，出错的原因 [打印本页]

作者
Author: xujian 时间: 2016-12-21 14:32
标题: 只是改变了基组，方法没变，出错的原因
本帖最后由 xujian 于 2016-12-21 15:34 编辑

是一个关于PhOO自由基的计算，根据文献SI中的坐标，用文献中的方法和基组优化（B3LYP/6-31g*），没有问题，如附件一。还用B3LYP，基组为6-311++g(3df,3pd)之后，却出现opt不收敛的错误，附件二，请问这是为什么，应该如何解决，谢谢。ps：按照sob老师关于opt不收敛的帖子中的方法试过，但没有解决。 (, 下载次数 Times of downloads: 5) (, 下载次数 Times of downloads: 5)

作者
Author: Accelerator 时间: 2016-12-21 14:54
为什么要用这么大的基组优化呢。。

作者
Author: xujian 时间: 2016-12-21 15:31

Accelerator 发表于 2016-12-21 14:54
为什么要用这么大的基组优化呢。。

想着分子比较小，机时比较富裕就用了

作者
Author: I10140317 时间: 2016-12-21 16:56
本帖最后由 I10140317 于 2016-12-21 17:01 编辑

是因为优化过渡态出现的问题，用opt=(ts,noeigentest,calcfc)可能可以避免这个报错，但不确定是不是可以优化到正确的过渡态。
(, 下载次数 Times of downloads: 62)

作者
Author: sobereva 时间: 2016-12-21 17:05

xujian 发表于 2016-12-21 15:31
想着分子比较小，机时比较富裕就用了

用了也没意义，6-311G**足矣，反倒几何优化用那么大会被懂行的人笑话
那个报错加noeigen就完了，不要以为出现l9999都是几何优化不收敛

作者
Author: xujian 时间: 2016-12-22 08:51
本帖最后由 xujian 于 2016-12-22 09:06 编辑

sobereva 发表于 2016-12-21 17:05
用了也没意义，6-311G**足矣，反倒几何优化用那么大会被懂行的人笑话
那个报错加noeigen就完了，不要 ...

好的，谢谢sob老师。高斯说明书上说“有时虽然没有通过曲率测试，但过渡态结构优化计算还是收敛，而NoEigenTest 仅用于大型的预算。”难道没有通过曲率测试对计算结果的正确性没有影响吗？

作者
Author: xujian 时间: 2016-12-22 08:51
本帖最后由 xujian 于 2016-12-22 09:06 编辑

I10140317 发表于 2016-12-21 16:56
是因为优化过渡态出现的问题，用opt=(ts,noeigentest,calcfc)可能可以避免这个报错，但不确定是不是可以优 ...

谢谢您,高斯说明书上说“有时虽然没有通过曲率测试，但过渡态结构优化计算还是收敛，而NoEigenTest 仅用于大型的预算。”难道没有通过曲率测试对计算结果的正确性没有影响吗？

作者
Author: I10140317 时间: 2016-12-22 10:19
本帖最后由 I10140317 于 2016-12-22 10:21 编辑

xujian 发表于 2016-12-22 08:51
谢谢您,高斯说明书上说“有时虽然没有通过曲率测试，但过渡态结构优化计算还是收敛，而NoEigenTest 仅用 ...

具体的我不是很理解，但从自己和实验室的经验上，如果不加noeigentest的话，这个问题非常容易遇到，加上这个关键词在许多时候是可以得到正确的结构的，前提是结构合理的话。
https://www.researchgate.net/pos ... ng_TS_with_Gaussian
这个网址是关于这个问题的一些讨论，你可以看看？我也不是主要做这个方向的，只是之前涉及过一些。我认为，这个问题主要只是优化TS算法的问题，最终判断是否是过渡态，还是要通过频率分析的。

作者
Author: sobereva 时间: 2016-12-22 12:24

xujian 发表于 2016-12-22 08:51
谢谢您,高斯说明书上说“有时虽然没有通过曲率测试，但过渡态结构优化计算还是收敛，而NoEigenTest 仅用 ...

甭管手册里的说法
最传统的牛顿法要求初猜结构（或者优化到过渡态过程中的每个结构）都必须只有一个负曲率，否则没法收敛到过渡态。但高斯用的方法思路也来自牛顿法，但经过改良，没这个限制，即便结构负曲率的数目不为1，只要初猜结构离实际过渡态偏离不很远，还是能收敛到过渡态去

更详细的看此文
过渡态、反应路径的计算方法及相关问题
http://sobereva.com/44

作者
Author: xujian 时间: 2016-12-22 14:02

I10140317 发表于 2016-12-22 10:19
具体的我不是很理解，但从自己和实验室的经验上，如果不加noeigentest的话，这个问题非常容易遇到，加上 ...

感谢！

作者
Author: xujian 时间: 2016-12-22 14:02

sobereva 发表于 2016-12-22 12:24
甭管手册里的说法
最传统的牛顿法要求初猜结构（或者优化到过渡态过程中的每个结构）都必须只有一个负 ...

thx too much

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