计算化学公社

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

[新手求助] Zn3O3优化无法收敛该怎么办?

[复制链接 Copy URL]

21

帖子

0

威望

67

eV
积分
88

Level 2 能力者

各位老师及大佬长辈,我做的Zn3O3优化,使用了opt=calcall,opt=noeigentest,SCF=QC都无法收敛,报错了,这个该怎么办呢?请大佬指点下。

捕获.PNG (107.04 KB, 下载次数 Times of downloads: 6)

捕获.PNG

捕获2.PNG (148.81 KB, 下载次数 Times of downloads: 5)

捕获2.PNG

捕获3.PNG (173.49 KB, 下载次数 Times of downloads: 4)

捕获3.PNG

21

帖子

0

威望

67

eV
积分
88

Level 2 能力者

2#
 楼主 Author| 发表于 Post on 2023-10-13 13:41:25 | 只看该作者 Only view this author
我做优化的常规写法为
%nprocshared=20
%mem=25GB
#p opt freq b3lyp/genecp maxdisk=50GB empiricaldispersion=gd3bj
temperature=673.15

Title Card Required

0 1
Zn                -1.30585400   -1.10076000    0.04221500
O                  1.01511000    1.41459100   -0.14275800
Zn                -0.77796200    1.47388400   -0.01388400
O                 -2.28277300    0.44284400    0.10072300
Zn                 1.40139700   -0.48070800   -0.07246900
O                  0.29493900   -1.95342800   -0.02436700

O 0
6-311g(d,p)
****
Zn 0
SDD
****

Zn 0
SDD

这个有什么问题吗?该怎么改才能收敛呢?




205

帖子

3

威望

1894

eV
积分
2159

Level 5 (御坂)

3#
发表于 Post on 2023-10-13 13:59:24 | 只看该作者 Only view this author
光看这个结构感觉很诡异,可能是无法收敛的主要原因,建议确定模型构建是否合理;

此外不建议用b3lyp做过渡金属的计算

21

帖子

0

威望

67

eV
积分
88

Level 2 能力者

4#
 楼主 Author| 发表于 Post on 2023-10-13 14:06:35 | 只看该作者 Only view this author
我的整个体系都是b3lyp做的,这个结构和别人发的文章结构也是一样的,但是我的就是无法收敛呢?

5万

帖子

99

威望

5万

eV
积分
112451

管理员

公社社长

5#
发表于 Post on 2023-10-13 15:51:40 | 只看该作者 Only view this author
没事甭用SCF=QC,好好看下文
解决SCF不收敛问题的方法
http://sobereva.com/61
而且这是针对SCF不收敛的,跟你几何优化不收敛没丝毫关系

甭写noeigentest,本来优化极小点根本就不会做本征值的测试

确保结构已经做了对称化

北京科音自然科学研究中心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!

21

帖子

0

威望

67

eV
积分
88

Level 2 能力者

6#
 楼主 Author| 发表于 Post on 2023-10-13 17:01:34 | 只看该作者 Only view this author
老师,您说的确保结构做了对称性是什么意思呢?我就只是将以上结构做优化,并没有做什么对称性,所以是这个问题吗?

111

帖子

0

威望

2082

eV
积分
2193

Level 5 (御坂)

7#
发表于 Post on 2023-10-14 11:33:07 | 只看该作者 Only view this author
这能量看起来不是已经快收敛了吗?加大optcyc试试

5万

帖子

99

威望

5万

eV
积分
112451

管理员

公社社长

8#
发表于 Post on 2023-10-14 17:54:40 | 只看该作者 Only view this author
wurui 发表于 2023-10-13 17:01
老师,您说的确保结构做了对称性是什么意思呢?我就只是将以上结构做优化,并没有做什么对称性,所以是这个 ...

GaussView里用点群工具做对称化再保存输入文件
有高对称性的体系如果结构轻微偏离对称性,很容易在优化时震荡。而先做对称化后,优化得又快又准又容易收敛
北京科音自然科学研究中心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!

本版积分规则 Credits rule

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

GMT+8, 2024-11-26 10:02 , Processed in 1.079418 second(s), 23 queries , Gzip On.

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