计算化学公社

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

[CP2K] 求助:使用cp2k对γ-Al2O3结构优化后和文献优化结果不同

[复制链接 Copy URL]

36

帖子

0

威望

244

eV
积分
280

Level 3 能力者

各位老师好,我使用MS软件按照图1中的数据建立了γ-Al2O3的晶胞,接着在multiwfn里分别设置理论方法PBEsol,设置k点3,2,2,设置任务优化原子和坐标,其他保持默认不变,生成cp2k的输入文件进行计算,最后收敛,但是优化后的结果和别人文献中的结果有差别,想问各位老师这种情况正常吗,我优化的结果是否是正确的
文献1中优化结果为:a=5.59 Å, b=8.42 Å, c=8.09 Å , β=90.58°
文献2中优化结果为:a=5.585 Å,b=8.406 Å,c=8.079 Å,β=90.555°
我的优化结果为:a=5.54622 Å,b=8.36104 Å,c=8.07013 Å,β=90.514°

联想截图_20230912201638.png (79.08 KB, 下载次数 Times of downloads: 34)

图1 文献中γ-Al2O3的数据

图1 文献中γ-Al2O3的数据

Al2O3-3.inp

6.25 KB, 下载次数 Times of downloads: 10

cp2k输入文件

Al2O3-3-1.restart

7.07 KB, 下载次数 Times of downloads: 4

restart文件

γAl2O3.cif

1.95 KB, 下载次数 Times of downloads: 23

根据图1建立的γ-Al2O3结构

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

2#
发表于 Post on 2023-9-12 21:55:26 | 只看该作者 Only view this author
基本没相差多少
差的那一点是基组、CUTOFF、收敛限之类因素
北京科音自然科学研究中心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

36

帖子

0

威望

244

eV
积分
280

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2023-9-12 23:27:33 | 只看该作者 Only view this author
sobereva 发表于 2023-9-12 21:55
基本没相差多少
差的那一点是基组、CUTOFF、收敛限之类因素

谢谢sob老师,那我这个优化是不是也没有问题,可以继续进行后续切面等操作

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

4#
发表于 Post on 2023-9-12 23:35:14 | 只看该作者 Only view this author
卡卡西永恒 发表于 2023-9-12 23:27
谢谢sob老师,那我这个优化是不是也没有问题,可以继续进行后续切面等操作

没有明显问题
如果你较真的话,也可以把基组提升一点(TZ档次)、CUTOFF提升一点、收敛限设严一点,看看是否相符程度更高
北京科音自然科学研究中心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

36

帖子

0

威望

244

eV
积分
280

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2023-9-13 18:06:49 | 只看该作者 Only view this author
sobereva 发表于 2023-9-12 23:35
没有明显问题
如果你较真的话,也可以把基组提升一点(TZ档次)、CUTOFF提升一点、收敛限设严一点,看看 ...

sob老师,我把原来的基组改为了TZVP-MOLOPT-GTH,运算的时候就报错了,这是什么原因呢

微信图片_20230913180307.jpg (456.32 KB, 下载次数 Times of downloads: 36)

运行图

运行图

Al2O35.out

6.71 KB, 下载次数 Times of downloads: 1

out文件

Al2O35.inp

6.33 KB, 下载次数 Times of downloads: 1

输入文件

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

6#
发表于 Post on 2023-9-13 21:56:01 | 只看该作者 Only view this author
卡卡西永恒 发表于 2023-9-13 18:06
sob老师,我把原来的基组改为了TZVP-MOLOPT-GTH,运算的时候就报错了,这是什么原因呢

提示写明了,此基组在基组文件里对此元素没有定义。
BASIS_MOLOPT_UZH里的基组种类丰富得多
北京科音自然科学研究中心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

36

帖子

0

威望

244

eV
积分
280

Level 3 能力者

7#
 楼主 Author| 发表于 Post on 2023-9-14 16:05:08 | 只看该作者 Only view this author
sobereva 发表于 2023-9-13 21:56
提示写明了,此基组在基组文件里对此元素没有定义。
BASIS_MOLOPT_UZH里的基组种类丰富得多

您是指下面图中的这些基组吗,我看Multiwfn里好像没有这些基组,是要手动修改输入文件吗,请问sob老师有没有关于修改基组的推荐资料呀,对这还不太了解,我手动改了一下报错了

截图 2023-09-14 15-56-59.png (129.05 KB, 下载次数 Times of downloads: 33)

截图 2023-09-14 15-56-59.png

6万

帖子

99

威望

5万

eV
积分
120109

管理员

公社社长

8#
发表于 Post on 2023-9-14 19:00:25 | 只看该作者 Only view this author
卡卡西永恒 发表于 2023-9-14 16:05
您是指下面图中的这些基组吗,我看Multiwfn里好像没有这些基组,是要手动修改输入文件吗,请问sob老师有 ...

自己改输入文件里的基组/赝势文件名、基组/赝势名就完了

这都是CP2K最基本的使用常识知识,在北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/workshop/KFP_content.html)里面有非常全面详细的讲解
北京科音自然科学研究中心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

36

帖子

0

威望

244

eV
积分
280

Level 3 能力者

9#
 楼主 Author| 发表于 Post on 2023-9-15 10:10:54 | 只看该作者 Only view this author
sobereva 发表于 2023-9-14 19:00
自己改输入文件里的基组/赝势文件名、基组/赝势名就完了

这都是CP2K最基本的使用常识知识,在北京科音 ...

好的,谢谢sob老师

20

帖子

0

威望

140

eV
积分
160

Level 3 能力者

10#
发表于 Post on 2025-7-6 14:26:15 | 只看该作者 Only view this author
请问这是参考的那一篇文献呀

本版积分规则 Credits rule

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

GMT+8, 2025-8-14 06:46 , Processed in 0.171238 second(s), 23 queries , Gzip On.

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