计算化学公社

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

[CP2K] cp2k使用BEEFVDW泛函时报错 CPASSERT failed

[复制链接 Copy URL]

25

帖子

0

威望

260

eV
积分
285

Level 3 能力者

跳转到指定楼层 Go to specific reply
楼主
本帖最后由 喝酸奶的巧乐兹 于 2025-12-28 03:09 编辑

各位老师好!
我想用cp2k计算Fe2O3 001表面吸附NO的吸附能,使用PBE+DZVP-MOLOPT-SR-GTH+D3(BJ)+DFT+U完成几何结构优化,使用的是二维周期性MT
随后使用优化完成的结构进行能量计算,使用multiwfn生成输入文件,设置为BEEFVDW+TZV2P-MOLOPT-SR-GTH+DFT+U,使用的是二维周期性MT,运行后报错 CPASSERT failed。不知如何修改
输入输出文件已附上,还请各位老师不吝赐教,谢谢!

Fe2O3-NO-energyshuiping-before-2.inp

14.97 KB, 下载次数 Times of downloads: 3

输入文件

Fe2O3-NO-energyshuiping-before-2.out

149.38 KB, 下载次数 Times of downloads: 0

输出文件

6万

帖子

99

威望

5万

eV
积分
124662

管理员

公社社长

2#
发表于 Post on 2025-12-28 07:14:32 | 只看该作者 Only view this author
尝试RPBE,如果能算,用RPBE就完了,BEEF-VDW比它算化学吸附强不了多少,还更难收敛、更贵、流行度低得多
北京科音自然科学研究中心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

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

3#
发表于 Post on 2025-12-28 10:46:58 | 只看该作者 Only view this author
可以上传一下“使用PBE+DZVP-MOLOPT-SR-GTH+D3(BJ)+DFT+U完成几何结构优化”这一步能SCF收敛的输入输出文件,以及讲讲DFT+U的值是选自哪篇文献吗?
√546=23.36664289109

25

帖子

0

威望

260

eV
积分
285

Level 3 能力者

4#
 楼主 Author| 发表于 Post on 2026-1-7 22:12:35 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2025-12-28 10:46
可以上传一下“使用PBE+DZVP-MOLOPT-SR-GTH+D3(BJ)+DFT+U完成几何结构优化”这一步能SCF收敛的输入输出文件 ...

这个输入文件

geo.inp

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

1208

帖子

6

威望

2430

eV
积分
3758

Level 5 (御坂)

傻傻的木瓜

5#
发表于 Post on 2026-1-7 22:30:27 | 只看该作者 Only view this author

谢了,可惜不是我想看的对角化带smearing和高k点算的,即使抛开CP2K里DFT+U暂时不兼容高k点这个因素,拿OT算过渡金属氧化物表面能SCF收敛也心里不踏实(参考http://bbs.keinsci.com/thread-27793-1-1.html的意见)。
√546=23.36664289109

25

帖子

0

威望

260

eV
积分
285

Level 3 能力者

6#
 楼主 Author| 发表于 Post on 2026-1-7 22:37:12 | 只看该作者 Only view this author
Uus/pMeC6H4-/キ 发表于 2026-1-7 22:30
谢了,可惜不是我想看的对角化带smearing和高k点算的,即使抛开CP2K里DFT+U暂时不兼容高k点这个因素,拿O ...

是的,Fe2O3表面吸附小分子后,即便使用OT收敛后的波函数,使用gamma+对角化计算能量都不一定能收敛

本版积分规则 Credits rule

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

GMT+8, 2026-1-23 13:15 , Processed in 0.169847 second(s), 23 queries , Gzip On.

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