计算化学公社

标题: cp2k单原子能量计算问题 [打印本页]

作者
Author: Penson 时间: 2023-12-14 14:32
标题: cp2k单原子能量计算问题
用cp2k算了一些单原子的单点能计算，在固定的的盒子尺寸内均放置一个原子（带相应电荷，当作离子），发现了一些结果比较怪异的地方：Al3+的绝对能量远高于Mg2+，这个是正常现象吗，具体的输入文件关键参数我都贴在下方，请各位大佬指教一下，是否有其合理性，万分感谢！
(, 下载次数 Times of downloads: 20) (, 下载次数 Times of downloads: 20) (, 下载次数 Times of downloads: 16) (, 下载次数 Times of downloads: 17)

(, 下载次数 Times of downloads: 17) (, 下载次数 Times of downloads: 20) (, 下载次数 Times of downloads: 17) (, 下载次数 Times of downloads: 18)

作者
Author: Assiduitas 时间: 2025-4-14 10:18
请问楼主最后单原子能量计算的问题解决了吗？听说VASP存在单原子计算能量不准确的问题，cp2k不确定是否也存在这个问题。

作者
Author: Penson 时间: 2025-4-20 11:08

Assiduitas 发表于 2025-4-14 10:18
请问楼主最后单原子能量计算的问题解决了吗？听说VASP存在单原子计算能量不准确的问题，cp2k不确定是否也存 ...

并没有呢，不确定问题出在哪里，感觉Al3+体系的结果问题比较大

作者
Author: Uus/pMeC6H4-/キ 时间: 2025-4-20 11:35
本帖最后由 Uus/pMeC6H4-/キ于 2025-4-20 11:36 编辑

比较不同体系的绝对能量一般没啥意义。如果一定要这样讨论，至少先考虑这三点：(1)统一赝势基组/赝势保留的非赝化的价电子范围，具体来说就是把Mg从默认的q10版换成q2版以与Al的q3版对应（Al没有q11版）；(2)统一原子在盒子内坐标，最好也对盒子尺寸做一个收敛性测试（虽然看当前正交盒子感觉是为了与某个基底吸附离子的体系一致而设置的，但不妨换成方便讨论的立方盒子）；(3)尝试按http://bbs.keinsci.com/thread-52560-1-1.html二楼这样的写法设置合理的初猜电子状态。

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