一句话提要
VASPKIT v1.00可以批量生成施加应变以后的结构,并根据能量-应变关系求解弹性常数矩阵。
计算能带反折叠—VASPKIT v1.00新功能
计算跃迁电偶极矩—VASPKIT v1.00新功能
深入分析能带结构(二)-VASPKIT能带图计算
VASPKIT计算电荷密度差
在材料的线性形变范围内(应变较小的情况下),体系的应力与应变满足胡克定律,
其中 σ 和 ε 分别表示应力和应变,Cij是弹性刚度常数,这里 1 ≤ i ≤ 6 即应变和应力分别有6个独立分量。施加应变后体系的应变能可按应变张量 ε 进行泰勒级数展开并取二阶近似得到:
这里E(V, ε)和E(V0, 0)别表示施加应变后和基态构型的总能,V0是平衡体积。
因此,采用第一性原理计算弹性常数有两种方法,第一种方法是应力-应变方法,即通过给结构施加不同应变,分别计算出所对应的应力大小,然后利用公式(1)拟合得到一次项系数,从而得到Ci