关于MOF-801体积模量的计算值与实验或文章中的理论值差异较大的问题

jqian0402

我在使用PBE-sol + DZVP-MOLOPT-SR-GTH的配置对MOF-801材料计算体积模量的过程中发现，计算得到的数值“35.78 GPa”远高于文献中的理论值“7.6 GPa”，也远高于我使用力场计算得到的数值。但我针对其他MOF材料进行的同类计算得到的结果与其他文献的结果又是相对匹配的。想问一下我的计算出了什么问题，才导致如此巨大的数值差异？

我的计算方法是：
1、对MOF-801进行结构优化。
2、计算晶胞边长分别变化1%、2%、3%时的压力。
3、根据公式：B = -P / (ΔV / V)，计算体积模量。

附件是我边长增大1%后的结构文件和配置文件。

archer

VASPKIT可以生成要计算的单点，拟合弹性模量