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

gch123 · 发表于 2025-10-14 09:53:13

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

vaspkit · 发表于 2025-10-14 14:20:42

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

gch123 · 发表于 2025-10-15 14:50:10

原胞的晶格常数为 3.21*3.21*5.21，K点为18*18*12，扩胞后为9.63*9.63*10.42，K点为6*6*6，计算出来的杨氏模量为44Gpa，当K点变为7*7*7时，计算出来的杨氏模量为61.889Gpa，想请教一下是K点还不够吗

		自动登录	找回密码
密码			立即注册