计算化学公社

标题: cp2k使用BEEFVDW泛函时报错 CPASSERT failed [打印本页]

作者
Author: 喝酸奶的巧乐兹 时间: 2025-12-28 03:09
标题: cp2k使用BEEFVDW泛函时报错 CPASSERT failed
本帖最后由喝酸奶的巧乐兹于 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。不知如何修改
输入输出文件已附上，还请各位老师不吝赐教，谢谢！

作者
Author: sobereva 时间: 2025-12-28 07:14
尝试RPBE，如果能算，用RPBE就完了，BEEF-VDW比它算化学吸附强不了多少，还更难收敛、更贵、流行度低得多

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-12-28 10:46
可以上传一下“使用PBE+DZVP-MOLOPT-SR-GTH+D3(BJ)+DFT+U完成几何结构优化”这一步能SCF收敛的输入输出文件，以及讲讲DFT+U的值是选自哪篇文献吗？

作者
Author: 喝酸奶的巧乐兹 时间: 2026-1-7 22:12

Uus/pMeC6H4-/キ发表于 2025-12-28 10:46
可以上传一下“使用PBE+DZVP-MOLOPT-SR-GTH+D3(BJ)+DFT+U完成几何结构优化”这一步能SCF收敛的输入输出文件 ...

这个输入文件

作者
Author: Uus/pMeC6H4-/キ 时间: 2026-1-7 22:30

喝酸奶的巧乐兹发表于 2026-1-7 22:12
这个输入文件

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

作者
Author: 喝酸奶的巧乐兹 时间: 2026-1-7 22:37

Uus/pMeC6H4-/キ发表于 2026-1-7 22:30
谢了，可惜不是我想看的对角化带smearing和高k点算的，即使抛开CP2K里DFT+U暂时不兼容高k点这个因素，拿O ...

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

欢迎光临计算化学公社 (http://bbs.keinsci.com/)