计算化学公社

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

[新手求助] OPT结果扭曲

[复制链接 Copy URL]

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

本帖最后由 小书童 于 2016-1-12 18:20 编辑

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



11.png (62.19 KB, 下载次数 Times of downloads: 29)

最终向得到的结果

最终向得到的结果

12.png (61.34 KB, 下载次数 Times of downloads: 34)

优化一定步数结果,也是想要的结构

优化一定步数结果,也是想要的结构

13.png (67.29 KB, 下载次数 Times of downloads: 42)

最终计算结果,扭曲了

最终计算结果,扭曲了

HCN-CC-NO-OPT.gjf

2.05 KB, 阅读权限: 1, 下载次数 Times of downloads: 9

输入

HCN-CC-NO-OPT.log

2.31 MB, 阅读权限: 1, 下载次数 Times of downloads: 4

输出

207

帖子

1

威望

3855

eV
积分
4082

Level 6 (一方通行)

2#
发表于 Post on 2016-1-12 19:31:05 | 只看该作者 Only view this author
你这个分子电荷是0吗?因为去了一个H。

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2016-1-12 20:13:10 | 只看该作者 Only view this author
bomsaude 发表于 2016-1-12 19:31
你这个分子电荷是0吗?因为去了一个H。

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

6万

帖子

99

威望

5万

eV
积分
120134

管理员

公社社长

4#
发表于 Post on 2016-1-12 21:10:38 | 只看该作者 Only view this author
没有虚频,说明在当前计算级别下这就是合理的极小点结构。你也可以再试试其它级别,比如M062X/TZVP,看是否都是这样。
这种结构极小点结构不是平面的是正常的事,因为有五元环和六元杂环存在,由于一些张力因素可以导致体系稍微变形成碗状。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

5#
 楼主 Author| 发表于 Post on 2016-1-12 21:42:27 | 只看该作者 Only view this author
sobereva 发表于 2016-1-12 21:10
没有虚频,说明在当前计算级别下这就是合理的极小点结构。你也可以再试试其它级别,比如M062X/TZVP,看是否 ...

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

6万

帖子

99

威望

5万

eV
积分
120134

管理员

公社社长

6#
发表于 Post on 2016-1-12 22:07:58 | 只看该作者 Only view this author
小书童 发表于 2016-1-12 21:42
那么您的意思是可能这个结构就是变形的。还有,结构优化一直平坦的走很久的这种我记得您说过怎么可以快点 ...



用gdiis,可以适当改用loose收敛限
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

7#
 楼主 Author| 发表于 Post on 2016-1-13 11:27:06 | 只看该作者 Only view this author
sobereva 发表于 2016-1-12 22:07

用gdiis,可以适当改用loose收敛限

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

6万

帖子

99

威望

5万

eV
积分
120134

管理员

公社社长

8#
发表于 Post on 2016-1-13 17:16:10 | 只看该作者 Only view this author
小书童 发表于 2016-1-13 11:27
老师按照您的方法,从之前优化的步数中选择一个比较合适的点(没扭曲的),然后用gdiis方式进行OPT,但是 ...

所以没优化到极小点。
如果得到的是平面结构,有虚频那是必然的,因为当前级别下极小点是扭曲的。
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

407

帖子

0

威望

856

eV
积分
1263

Level 4 (黑子)

9#
 楼主 Author| 发表于 Post on 2016-1-13 20:04:39 | 只看该作者 Only view this author
sobereva 发表于 2016-1-13 17:16
所以没优化到极小点。
如果得到的是平面结构,有虚频那是必然的,因为当前级别下极小点是扭曲的。

恩恩,了然了,谢谢老师!

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 08:06 , Processed in 0.219598 second(s), 23 queries , Gzip On.

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