vasp 计算杨氏弹性_基于能量应变关系计算弹性常数—VASPKIT v1.00新功能

一句话提要

VASPKIT v1.00可以批量生成施加应变以后的结构，并根据能量-应变关系求解弹性常数矩阵。

计算能带反折叠—VASPKIT v1.00新功能

计算跃迁电偶极矩—VASPKIT v1.00新功能

深入分析能带结构(二)-VASPKIT能带图计算

VASPKIT计算电荷密度差

在材料的线性形变范围内(应变较小的情况下)，体系的应力与应变满足胡克定律，

其中 σ 和 ε 分别表示应力和应变，Cij是弹性刚度常数，这里 1 ≤ i ≤ 6 即应变和应力分别有6个独立分量。施加应变后体系的应变能可按应变张量 ε 进行泰勒级数展开并取二阶近似得到：

这里E(V, ε)和E(V0, 0)别表示施加应变后和基态构型的总能，V0是平衡体积。

因此，采用第一性原理计算弹性常数有两种方法，第一种方法是应力-应变方法，即通过给结构施加不同应变，分别计算出所对应的应力大小，然后利用公式(1)拟合得到一次项系数，从而得到Ci