计算化学公社

标题: Gaussian计算中结构优化不收敛问题 [打印本页]

作者
Author: 柚子不开花 时间: 2017-10-16 17:22
标题: Gaussian计算中结构优化不收敛问题
老师，您好！我最近在用高斯重复前辈文献的工作，寻找P45O和甲苯反应的中间体时（体系有55个原子，包含C、H、O、N、S、Fe），当体系的自旋多重度为二重态时，成功收敛，虚频唯一，过渡态基本正确。但改变多重度为四重态后，关键词不变，可体系一直无法收敛，显示L502错误，观察优化曲线，有收敛趋势比较平滑，但每次都在34步左右时就会报错。尝试了基础课上学的，增加步数上限，调整步长上限，因为含有过渡金属还尝试了能级移动的方法，但依然报相同的错误，在34步左右任务就死掉了。记得老师基础课上讲过，对于金属配合物，自旋多重度设定的不合理会导致SCF不收敛。但我是重复前辈的工作，他可以成功找到，应该不会是这个原因，对吗？附件中给出了输入和输出文件，问题还没能解决，望老师百忙之中予以指点，不胜感激。

作者
Author: ssp2493187 时间: 2017-10-16 17:26
有没有尝试过用这个帖子解决 http://sobereva.com/164

作者
Author: 柚子不开花 时间: 2017-10-16 17:31

ssp2493187 发表于 2017-10-16 17:26
有没有尝试过用这个帖子解决 http://sobereva.com/164

谢谢您的回复，我仔细阅读国这个帖子，基组和方法是文献给出的不能更换，其他的相关的方法也尝试了的，还是同样的错误。还有我有些不细心帖子里有个错误，是在寻找过渡态，而不是中间体。

作者
Author: sobereva 时间: 2017-10-16 17:37
不知道你说的“基础课”是什么，如果是科音的基础量化班，每次课上都着重强调不要试图盲目用maxcyc解决问题，而且每次培训班里也强调绝对不要用6-31G。换个像样的基组，不行到时候再说（如果之前的人就是用的6-31G，重复这种糟糕的数据也没有任何意义）

作者
Author: 柚子不开花 时间: 2017-10-16 18:21

sobereva 发表于 2017-10-16 17:37
不知道你说的“基础课”是什么，如果是科音的基础量化班，每次课上都着重强调不要试图盲目用maxcyc解决问题 ...

非常感谢Sob老师的回复，我是看每次都在34步左右任务就终止了，并且从优化曲线上看能量和受力还在继续下降，所以试了试增加步长上限。老师我还想问一下，什么原因会导致能量和受力还在继续下降但任务就终止了呢？可以尝试怎么解决呀？

作者
Author: liyuanhe211 时间: 2017-10-16 18:25

柚子不开花发表于 2017-10-16 18:21
非常感谢Sob老师的回复，我是看每次都在34步左右任务就终止了，并且从优化曲线上看能量和受力还在继续下 ...

这是SCF不收敛，不是几何优化不收敛。
先把基组换了，再出问题再说

作者
Author: 柚子不开花 时间: 2017-10-16 18:32

liyuanhe211 发表于 2017-10-16 18:25
这是SCF不收敛，不是几何优化不收敛。
先把基组换了，再出问题再说

好的，我试试。谢谢老师。

作者
Author: ssp2493187 时间: 2017-10-17 17:16
我少贴了帖子不过老师们都已经提出中肯建议了我有其他别的想法说一下：
首先，做计算不能拿着东西就扔给服务器就去算，要考虑是否有意义，是否能够达到预期的目标。
其次，如果基础不太好，一定要看书做练习，如果可以的话去sob老师的培训班，就不会那么多关键词了和错误关键词了。
最后，提问一定要给出自己尝试过的解决的方式，为什么没有成功，报错是什么，尽量避免一些关键词写错了，内存不足或者某些设置不合理的情况。

以上都不是针对楼主的，是对所有新接触量子化学计算的新手包括我自己。大家一起跟着老师们进步吧！

作者
Author: 柚子不开花 时间: 2017-10-19 09:26

ssp2493187 发表于 2017-10-17 17:16
我少贴了帖子不过老师们都已经提出中肯建议了我有其他别的想法说一下：
首先，做计算不能拿着东西就扔给 ...

很多都是我的问题，非常感谢您的建议。

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