计算化学公社

标题: 求助：cp2k计算Na的内聚能与实验值相差较大 [打印本页]

作者
Author: Momooooooooooo 时间: 2024-4-23 20:57
标题: 求助：cp2k计算Na的内聚能与实验值相差较大
老师们好
本人是cp2k新手，最近想要计算钠的内聚能，泛函用的是PBE，基组用的是TZV2P-MOLOPT-PBE-GTH-q9，POTENTIAL用GTH-PBE，截断能用的是1000 rel，cutoff是60

1、钠晶体是从materials project下载的cif进行2x2x1扩胞后的结构（含8个原子），k点取的是6x6x6，我在用multiwfn生成输入文件时不知道选择进行什么结构优化，所以3（固定晶胞的结构优化）和4（晶胞+结构优化）都尝试了，再用shermo进行振动分析
2、钠原子的电子能量计算的输入文件也是用multiwfn生成的，设置周期性为none，自选多重度为2，计算级别和钠晶体的opt相同
具体结果如下：
      固定晶胞的结构优化晶胞+结构优化
U(s)    -1282.376                -1283.144             eV/atom
U(g)    -1281.934                -1281.934          eV/atom
ΔU       0.441                         0.441             eV/atom
实验值：1.113 eV/atom （Charles Kittel的Introduction to Solid State Physics）
具体的输入文件也都在下方了请各位老师们帮忙看看我的输入文件是不是哪里出错了非常感谢

作者
Author: sobereva 时间: 2024-4-24 06:55
直接拿Na的原胞做优化多好，干嘛扩胞，很莫名其妙，只会额外浪费时间。
这么小的晶胞做振动分析算热力学量没有意义，如果用Shermo算热力学量需要对足够大的超胞做振动分析。有对称性的小晶胞应当结合用phonopy算热力学量，因为phonopy能考虑晶体对称性，靠有限差分计算受力的次数少得多得多。北京科音CP2K第一性原理计算培训班（http://www.keinsci.com/workshop/KFP_content.html）
里的“热力学量与能量相关属性的计算”我给了具体的CP2K+phonopy算热力学量的例子。

另外，PBE计算Na的内聚能这事过于普通，因此大概率能从文献或书籍里找到其他人计算的值。先跟别人的PBE的结果比，之后再说和实验比的事。而且最好先去比基于电子能量算的内聚能，再说牵扯到热力学校正量的情况。

作者
Author: Momooooooooooo 时间: 2024-4-24 09:24
谢谢sob老师把我的问题解释的这么清楚感谢

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