计算化学公社

 找回密码 Forget password
 注册 Register

使用CP2K计算Fe2O3表面吸附NH3的吸附能与文献差距过大

查看数: 255 | 评论数: 4 | 收藏 Add to favorites 2
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2025-12-15 23:03

正文摘要:

之前发过使用Fe2O3吸附NH3吸附能与文献有差距使用CP2K计算Fe2O3表面吸附结合能与文献差距较大(http://bbs.keinsci.com/forum.php ... 57378&fromuid=82794),经过sob老师指点得知所得结果可以使用,随后计算Fe2O ...

回复 Reply

喝酸奶的巧乐兹 发表于 Post on 2026-1-7 22:14:53
后记:经过各种条件比较,问题解决了,是Fe2O3的表面自旋设置不恰当,表面从自旋向上改为自旋向下就合理了
sobereva 发表于 Post on 2025-12-16 00:33:00
喝酸奶的巧乐兹 发表于 2025-12-16 00:23
谢谢sob老师回答!我再试试换更高精度基组。
另外还有个问题,我使用TZVP-MOLOPT-GTH计算能量时,给Fe的 ...

+U是理论方法层面的,和基组是两码事

参考北京科音CP2K第一性原理计算培训班(http://www.keinsci.com/KFP)讲DFT+U的幻灯片:

喝酸奶的巧乐兹 发表于 Post on 2025-12-16 00:23:06
sobereva 发表于 2025-12-16 00:02
注意下文说的各种问题
谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678(http:// ...

谢谢sob老师回答!我再试试换更高精度基组。
另外还有个问题,我使用TZVP-MOLOPT-GTH计算能量时,给Fe的DFT+U哪怕只有1eV都会报warning:DFT+U energy contibution is negative!进而scf难收敛。而关闭DFT+U后则可以收敛,请问这是说明精度高的基组不需要DFT+U吗?
sobereva 发表于 Post on 2025-12-16 00:02:19
注意下文说的各种问题
谈谈重复不出来计算化学文献里的数据的可能原因
http://sobereva.com/678http://bbs.keinsci.com/thread-38911-1-1.html
诸如色散校正要和文献里用的尽量一样

建议基组改用TZV2P-MOLOPT-GTH
如果能找到CP2K算此体系的合适的+U参数,算吸附能时也要加上
注意所有计算细节,如自旋多重度、是否收敛到正确的磁性状态等

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

GMT+8, 2026-1-23 21:54 , Processed in 0.441342 second(s), 25 queries , Gzip On.

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