帮同学做ZrC (Fm-3m结构,也就是氯化钠结构。实验晶格常数为4.68埃)的EOS拟合,用的PBE-PAW赝势。选择截断能500eV,为了速度,精度选择了Normal~然后K点是13x13x13。
POSCAR的设置为: 从OUTCAR的内容来看,结构是正确读取了的: 分别使用了4.6和5.2计算,结果得到的E-V图居然是(使用a.u.单位是为了使用eos程序方便):
![]()
出现跳点的位置是4.67埃,而这个点也是直接使用结构优化得到的点……
看了一下某版热心提供的新版PBE-PAW赝势,偏偏Zr的就没有更新……由于碳的赝势已经计算过不少东西了,暂时没发现什么问题,所以现在怀疑是Zr的赝势本身的问题……毕竟就算K点不足,截断能不够,也不该就总是一个点跳吧……
计算的INCAR,KPOINTS和POSCAR都在附件中~ |