应变-能量关系计算超胞弹性模量不准

gch123

请教一下 我现在计算hcp Mg的弹性模型，用原胞计算的结果比较合理，但用超胞后（3*3*2，36个原子），计算出来的弹性模量相对原胞有较大的差距，请问一下，这应该怎么修正呢？ 用超胞算，是因为后面要掺杂

vaspkit

检查K点密度，截断能是否保持一致？另外，如果掺杂，对称性会发生改变，可以通过SYMMETRY.in来强制指定空间群。