计算化学公社

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

[CP2K] Ni金属和复合金属催化中,vasp是否要比cp2k结构优化更快收敛

[复制链接 Copy URL]

9

帖子

0

威望

283

eV
积分
292

Level 3 能力者

如题,本人刚学了一个月cp2k,研究Ni吸附丙烷分子,切了一个Ni111表面。进行结构优化,发现,Multiwfn默认cp2k设置,SCF收敛极难(或者干脆不收敛),然后把同样晶胞放入vasp当中,只用了不到300s优化结束。当把VASP结构优化的CONTAR在CP2K当中进行单点能计算,以期望得到molden文件,进行波函数分析。但即便使用优化后的结果,单点能收敛也极为困难。所以问题:
1、是否在纯金属当中,因为金属赝势问题,导致vasp的纯平面波要优于cp2k的基函数结合平面波?
2、是否可以使用vasp结构优化后的结果放入cp2k进行单点能进行波函数分析?
3、波函数分析是否非常依赖K点,仅靠单点能的gamma点可以吗?

Ps:之前跑的结构优化cp2k程序被我删掉了,我稍微跑几步演示下,有明显的震荡。
      我通过缩小MIXING的ALPHA和增大NBROYDEN,降低TRUST_RADIUS已经能让单点能收敛。但结构优化无法进行。

202411271822254909..png (39.31 KB, 下载次数 Times of downloads: 2)

切的Ni111表面

切的Ni111表面

捕获.JPG (592.57 KB, 下载次数 Times of downloads: 2)

结构优化的震荡

结构优化的震荡

Ni111.inp

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

INCAR

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

POSCAR

3.85 KB, 下载次数 Times of downloads: 2

OUTCAR

532.6 KB, 下载次数 Times of downloads: 2

Ni111_880710.log

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

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

2#
发表于 Post on 2024-11-27 23:33:21 | 只看该作者 Only view this author
1 CP2K算Ni表面就是SCF难收敛,可以用QE等平面波程序

2 完全语义不明。CP2K又不是Multiwfn那样的波函数分析程序

3 考虑k点时CP2K根本没法产生molden文件
北京科音自然科学研究中心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

9

帖子

0

威望

283

eV
积分
292

Level 3 能力者

3#
 楼主 Author| 发表于 Post on 2024-11-28 08:58:48 | 只看该作者 Only view this author
sobereva 发表于 2024-11-27 23:33
1 CP2K算Ni表面就是SCF难收敛,可以用QE等平面波程序

2 完全语义不明。CP2K又不是Multiwfn那样的波函数 ...

谢谢社长。有些语意不明
问题2:不是使用cp2k做波函数分析。是而周期性体系里,cp2k结合Multiwfn太好用了,所以想使用丰富的研究方法,我得转回来用cp2k的molden文件。这时如果用vasp优化后的结构,是否合理?
问题3:是的,使用K点时无法产生molden文件。所以在这种情况下,我为了使用molden进行波函数分析,只有gamma点,这种情况下的分析是否合理?因为要满足30埃的条件,我应当使用331的k点进行能量计算

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

4#
发表于 Post on 2024-11-28 11:40:56 | 只看该作者 Only view this author
葡萄皮 发表于 2024-11-28 08:58
谢谢社长。有些语意不明
问题2:不是使用cp2k做波函数分析。是而周期性体系里,cp2k结合Multiwfn太好用 ...

2 VASP优化->CP2K单点得到molden文件->Multiwfn波函数分析是可以的

3 考虑k点做能量计算,扩胞成只需要gamma点就够用的情况产生molden文件用Multiwfn做波函数分析,是可以的
北京科音自然科学研究中心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

9

帖子

0

威望

283

eV
积分
292

Level 3 能力者

5#
 楼主 Author| 发表于 Post on 2024-11-28 15:25:55 | 只看该作者 Only view this author
sobereva 发表于 2024-11-28 11:40
2 VASP优化->CP2K单点得到molden文件->Multiwfn波函数分析是可以的

3 考虑k点做能量计算,扩胞成只需 ...

谢谢社长。收获颇丰

17

帖子

0

威望

376

eV
积分
393

Level 3 能力者

6#
发表于 Post on 2024-11-29 00:24:22 | 只看该作者 Only view this author
老早就在用VASP算一些CP2K不太好收敛的体系,然后结合Multiwfn进行波函数分析,今天见到社长回复心里更踏实了

123

帖子

0

威望

1550

eV
积分
1673

Level 5 (御坂)

7#
发表于 Post on 2024-11-29 09:34:51 | 只看该作者 Only view this author
sobereva 发表于 2024-11-28 11:40
2 VASP优化->CP2K单点得到molden文件->Multiwfn波函数分析是可以的

3 考虑k点做能量计算,扩胞成只需 ...

第二步这个方法里面vasp和CP2K的赝势不同,可以这样操作吗?总感觉不踏实

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

8#
发表于 Post on 2024-12-2 05:06:12 | 只看该作者 Only view this author
1138711019 发表于 2024-11-29 09:34
第二步这个方法里面vasp和CP2K的赝势不同,可以这样操作吗?总感觉不踏实

两个程序用于不同阶段,数据不直接横向对比,没有问题
北京科音自然科学研究中心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

183

帖子

0

威望

1028

eV
积分
1211

Level 4 (黑子)

9#
发表于 Post on 2025-4-16 22:01:01 | 只看该作者 Only view this author
sobereva 发表于 2024-11-27 23:33
1 CP2K算Ni表面就是SCF难收敛,可以用QE等平面波程序

2 完全语义不明。CP2K又不是Multiwfn那样的波函数 ...

老师,请问CP2K算Ni的氧化物的表面还存在SCF难收敛吗,谢谢
科音学员

6万

帖子

99

威望

5万

eV
积分
120132

管理员

公社社长

10#
发表于 Post on 2025-4-17 14:51:40 | 只看该作者 Only view this author
一条君 发表于 2025-4-16 22:01
老师,请问CP2K算Ni的氧化物的表面还存在SCF难收敛吗,谢谢

比Ni单质表面好收敛得多,但也不算很容易
北京科音自然科学研究中心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

本版积分规则 Credits rule

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

GMT+8, 2025-8-15 01:59 , Processed in 0.156676 second(s), 23 queries , Gzip On.

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