求助：cp2k计算Na的内聚能与实验值相差较大

Momooooooooooo

老师们好
本人是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）
具体的输入文件也都在下方了 请各位老师们帮忙看看我的输入文件是不是哪里出错了 非常感谢

Momooooooooooo

谢谢sob老师把我的问题解释的这么清楚 感谢

sobereva

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

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