关于弹性常数的单位换算

之前我们讲解过如何用VASP计算弹性常数：

计算弹性常数 – Qian Xia

最终得到的弹性常数矩阵的单位是KBar，但是有时候我们在文献里面(PHYSICAL REVIEW B 105, 235303 (2022))看到的结果是这样的：

可见，文献中的单位是N/m。

具体的，我们可以这样进行换算：
1 bar=100 KPa=100000 Pa,
1 Pa=1 N/m^2,
因此1bar=1*10^5 N/m^2
假设晶格沿着c方向上的厚度为a，最终1 bar=1*10^5*a N/m

我们可以把上面的文献重复一下，首先优化好结构，并设置真空层为20埃（20*10^-10 m），计算得到的弹性常数矩阵如下：

可以看到c11=467.8245 kBar，
转化一下：c11=467.8245 kBar=467.8245*10^8 N/m^2=467.8245*10^8*20*10^-10 N/m=93.56 N/m，和文献中的98.3 N/m很接近，还有一点误差可能是结构弛豫的原因或者U值的选取等等。