上面表达式中σi和εj分别表示应力和应变,Cij是弹性常数,其中1≤i≥6和1≤j≥6。通过给结构施加不同应变,分别计算出所对应的应力大小,然后利用上面公式拟合得到一次项系数,从而便可得到Cij。下图给出了立方晶系的具体形式: 下面我们以金刚石结构为例讲解如何采用应力-应变函数关系计算弹性常数,详见VASPKIT/examples/...
下面我们以金刚石结构为例讲解如何采用应力-应变函数关系计算弹性常数,详见VASPKIT/examples/elastic/diamond_3D。由于金刚石具有立方晶体结构,一共有3个独立弹性常数C11、C12和C44(不明白的请看原胞转化方法以及标准原胞在计算中的重要性)。 1. 准备优化好的POSCAR文件,注意通常采用具有标准基矢形式的原胞计算弹性常数...