计算化学公社

 找回密码 Forget password
 注册 Register
Views: 32988|回复 Reply: 8
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] Gaussian计算中结构优化不收敛问题

[复制链接 Copy URL]

27

帖子

0

威望

1365

eV
积分
1392

Level 4 (黑子)

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

TS-4-Re6.gjf

3.53 KB, 下载次数 Times of downloads: 157

TS-4-Re6.out

957.05 KB, 下载次数 Times of downloads: 85

27

帖子

0

威望

1365

eV
积分
1392

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2017-10-19 09:26:02 | 只看该作者 Only view this author
ssp2493187 发表于 2017-10-17 17:16
我少贴了帖子 不过老师们都已经提出中肯建议了 我有其他别的想法说一下:
首先,做计算不能拿着东西就扔给 ...

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

97

帖子

0

威望

448

eV
积分
546

Level 4 (黑子)

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

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

27

帖子

0

威望

1365

eV
积分
1392

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2017-10-16 18:32:51 | 只看该作者 Only view this author
liyuanhe211 发表于 2017-10-16 18:25
这是SCF不收敛,不是几何优化不收敛。
先把基组换了,再出问题再说

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

3098

帖子

28

威望

1万

eV
积分
16892

Level 6 (一方通行)

6#
发表于 Post on 2017-10-16 18:25:15 | 只看该作者 Only view this author
柚子不开花 发表于 2017-10-16 18:21
非常感谢Sob老师的回复,我是看每次都在34步左右任务就终止了,并且从优化曲线上看能量和受力还在继续下 ...

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

27

帖子

0

威望

1365

eV
积分
1392

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2017-10-16 18:21:47 | 只看该作者 Only view this author
sobereva 发表于 2017-10-16 17:37
不知道你说的“基础课”是什么,如果是科音的基础量化班,每次课上都着重强调不要试图盲目用maxcyc解决问题 ...

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

5万

帖子

99

威望

5万

eV
积分
112408

管理员

公社社长

4#
发表于 Post on 2017-10-16 17:37:12 | 只看该作者 Only view this author
不知道你说的“基础课”是什么,如果是科音的基础量化班,每次课上都着重强调不要试图盲目用maxcyc解决问题,而且每次培训班里也强调绝对不要用6-31G。换个像样的基组,不行到时候再说(如果之前的人就是用的6-31G,重复这种糟糕的数据也没有任何意义)
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入“北京科音”微信公众号获取北京科音培训的最新消息、避免错过网上有价值的计算化学文章!
欢迎加入人气非常高、专业性特别强的综合性理论与计算化学交流QQ群“思想家公社QQ群”:1号:18616395,2号:466017436,3号:764390338,搜索群号能搜到哪个说明目前哪个能加,合计9000人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大的量子化学波函数分析程序)
ResearchGate:https://www.researchgate.net/profile/Tian_Lu
Money and papers are rubbish, get a real life!

27

帖子

0

威望

1365

eV
积分
1392

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2017-10-16 17:31:12 | 只看该作者 Only view this author
ssp2493187 发表于 2017-10-16 17:26
有没有尝试过用这个帖子解决 http://sobereva.com/164

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

97

帖子

0

威望

448

eV
积分
546

Level 4 (黑子)

2#
发表于 Post on 2017-10-16 17:26:51 | 只看该作者 Only view this author
有没有尝试过用这个帖子解决 http://sobereva.com/164

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2024-11-25 21:51 , Processed in 0.342925 second(s), 24 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list