计算化学公社

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

[CP2K] 求助cp2k的赝势和基组选择

[复制链接 Copy URL]

43

帖子

0

威望

597

eV
积分
640

Level 4 (黑子)

各位老师好,最近刚学cp2k,有以下几点求助各位老师:1、该如何基于自己的体系选择基组和赝势呢?对于不同基组的精度和介绍可以在哪看呢

2、我用DZVP-MOLOPT-SR-GTH基组和GTH-PBE赝势计算了以下Fe2SiO4的结构优化(包含晶格参数优化),发现计算得到的晶格参数比实验值小很多,随着优化步数的增加,晶格参数越来越小,是我的设置有什么问题吗?麻烦大家帮忙看一下

Fe2SiO4.inp.out.zip

483 KB, 下载次数 Times of downloads: 42

输出文件

Fe2SiO4.inp

4.46 KB, 下载次数 Times of downloads: 135

输入文件

6万

帖子

99

威望

5万

eV
积分
120169

管理员

公社社长

2#
发表于 Post on 2021-5-27 05:09:41 | 只看该作者 Only view this author
以后发帖时标题里不要带叹号,给你改了

1 一般情况优化用DZVP-MOLOPT-SR-GTH就可以(赝势自然是结合标配的GTH赝势。如果有专门对你用的泛函优化的GTH赝势就用专用的,其它情况用PBE的凑合)。能量可以用更好的TZV2P-MOLOPT-GTH(但不是所有元素都有)。不同的MOLOPT系列基组的精度对比看MOLOPT基组原文J. Chem. Phys., 127, 114105 (2007)。对于结合multiwfn做波函数分析,pob-TZVP比较适合,见比如《使用Multiwfn结合CP2K通过NCI和IGM方法图形化考察固体和表面的弱相互作用》(http://sobereva.com/588)里的说明。

2 要贴就贴带完整结构信息的输入文件。应当考虑k点,单胞跑的话至少2 4 5。当前计算也不算是夸张的越来越小,也不要指望PBE优化得到的晶胞参数能和实验特别相符。还要注意对磁性的考虑。
北京科音自然科学研究中心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

43

帖子

0

威望

597

eV
积分
640

Level 4 (黑子)

3#
 楼主 Author| 发表于 Post on 2021-5-27 10:33:40 | 只看该作者 Only view this author
sobereva 发表于 2021-5-27 05:09
以后发帖时标题里不要带叹号,给你改了

1 一般情况优化用DZVP-MOLOPT-SR-GTH就可以(赝势自然是结合标配 ...

谢谢sob老师的解答,我用的基组和赝势跟您推荐的一样,我试试设置K点计算

32

帖子

0

威望

253

eV
积分
285

Level 3 能力者

4#
发表于 Post on 2022-3-8 16:24:30 | 只看该作者 Only view this author
本帖最后由 qq870624836 于 2022-3-8 19:07 编辑
sobereva 发表于 2021-5-27 05:09
以后发帖时标题里不要带叹号,给你改了

1 一般情况优化用DZVP-MOLOPT-SR-GTH就可以(赝势自然是结合标配 ...

C:\Users\Administrator\Desktop\QQ图片20220308162033,老师好,def2tzvp这个基组应该是可以计算Pb元素的把为什么在cp2k这里会报错呀。我该怎么解决呢

QQ图片20220308162033.png (18.87 KB, 下载次数 Times of downloads: 39)

QQ图片20220308162033.png

6万

帖子

99

威望

5万

eV
积分
120169

管理员

公社社长

5#
发表于 Post on 2022-3-9 04:54:53 | 只看该作者 Only view this author
qq870624836 发表于 2022-3-8 16:24
,老师好,def2tzvp这个基组应该是可以计算Pb元素的把为什么在cp2k这里会报错呀。我该怎么解决呢

CP2K里没法用def2-TZVP赝势基组(def2系列对第五周期开始是赝势基组)
北京科音自然科学研究中心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

27

帖子

0

威望

198

eV
积分
225

Level 3 能力者

6#
发表于 Post on 12 hour ago | 只看该作者 Only view this author
sobereva 发表于 2021-5-27 05:09
以后发帖时标题里不要带叹号,给你改了

1 一般情况优化用DZVP-MOLOPT-SR-GTH就可以(赝势自然是结合标配 ...

社长,请问周期性晶体进行变胞优化的时候也可以使用DZVP-MOLOPT-SR-GTH基组和GTH赝势吗,就是Multiwfn默认设置里的组合

本版积分规则 Credits rule

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

GMT+8, 2025-8-16 23:24 , Processed in 0.168170 second(s), 23 queries , Gzip On.

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