计算化学公社

标题: 高斯结构优化遇到L502报错 [打印本页]

作者
Author: 枕上卢生 时间: 2020-7-25 16:38
标题: 高斯结构优化遇到L502报错
老师：我使用mp2/aug-cc-pvtz方法优化，第一优化l502 出错,使用您介绍办法进行调整，可还是L502 错误，太苦恼了，不知该如何处理，希望老师帮忙解答？

作者
Author: ldatea 时间: 2020-7-25 17:12
本帖最后由 ldatea 于 2020-7-27 18:10 编辑

（第二段有错误，已修改，Dunning相关一致性基组还是可以用的，may-cc-pvtz就足够理想了）
以下大部分内容来源都是sob老师的博文，建议阅读的帖子6楼已经给出。
先回答你提出的问题，
你这是SCF超过最大收敛圈数了（默认情况下128圈就报错结束）用SCF=maxcyc=300在此处基本没用，你看一下趋势，就会发现是震荡，Gaussian默认情况下收敛判据是RMSDP<10-8，MaxDP<10-6，energy变化<10-6。你这里RMSDP和MaxDP经常出现在10-7，10-5，而且没有逐渐减小的趋势。再设大也基本没有用。
你根本没好好看sob老师的博文http://sobereva.com/61（《解决SCF不收敛问题的方法》），谁要是一上来就是增大收敛圈数，一定是菜鸟。推荐你这么干的也一定是菜鸟。
你如果想收敛，先尝试SCF(noincfock)，如果不行，用SCF(noincfock) int=acc2e=12。

再说你的方法问题，
如果不是因为重现结果，不要用MP2这种吃力不讨好的方法。大部分情况下MP2的表现还不如选得合适的DFT杂化泛函。中性有机分子几何优化对基组的要求比较低，def2，def基组也够用了，研究弱相互作用再加一层最低限度的弥散函数也够了。更严重的是，用aug-，加这么多弥散函数就是找虐。耗时急剧增长，很容易因为不收敛报错。

如果没有3原子处于同一条直线或者是含很多分子的团簇，用opt=cartesian是吃力不讨好的，非常愚蠢。

我只能给很粗糙的建议，
如果离收敛还很远（你目前的情况里收敛就很远），先用便宜的方法和小一点的基组比如B3LYP/6-31+g(d) em=gd3bj粗优化，接着用更好的方法和基组，根据你所需的精度。我没有足够经验，接下来的处理我不能给合理的建议。

作者
Author: thanhtam 时间: 2020-7-25 17:20
不要用mp2，用dft就行了。

作者
Author: qinzhong605 时间: 2020-7-25 17:26

这么大的组合在做优化，想啥呢

作者
Author: 枕上卢生 时间: 2020-7-25 17:34
各位boss,那这个问题如何解决呢？
更换什么基组最为合适，我们想研究他们间的弱作用

作者
Author: hebrewsnabla 时间: 2020-7-25 17:54
本帖最后由 hebrewsnabla 于 2020-7-25 18:01 编辑

枕上卢生发表于 2020-7-25 17:34
各位boss,那这个问题如何解决呢？
更换什么基组最为合适，我们想研究他们间的弱作用

二楼已经说了。

基组选择看 http://sobereva.com/336 ， http://sobereva.com/387
泛函选择看 http://sobereva.com/272

另外，既然是弱相互作用，还需要看
http://sobereva.com/210
http://bbs.keinsci.com/forum.php ... 72&fromuid=5840

作者
Author: sobereva 时间: 2020-7-26 01:23
如置顶的新社员必读贴和论坛首页的公告栏所示，求助帖必须在帖子标题明确体现出此帖内容是求助或提问，并反映出帖子具体内容，仔细看http://bbs.keinsci.com/thread-9348-1-1.html。我已把你的不恰当标题“高斯结构优化”改了，以后务必注意

mp2/aug-cc-pvtz真是超级费力不讨好，没常识的人才这么做

作者
Author: 枕上卢生 时间: 2020-7-29 14:57
好的，谢谢Boss

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