计算弹性常数

准备文件：INCAR POSCAR POTCAR KPOINTS
INCAR设置：#==================以内的内容
#====================================================
ENCUT = 600
EDIFF = 1e-6
ISMEAR = 0
SIGMA = 0.05
POTIM = 0.1
LCHARG = FALSE
LWAVE = FALSE
IBRION=6                #弹性常数关键的三个参数
ISIF=3                  #
NFREE=4                 #
EDIFFG = -0.001

LMAXMIX = 4             !d电子为4，f电子为6
LDAU=T                  !这两项不要改
LDAUTYPE=2              !
LDAUL= 2 -1 -1
LDAUU= 4.4 0.0 0.0      !加u的那个值对应的不为0
LDAUJ= 0.0 0.0 0.0
LASPH=T

ISYM=2
SYMPREC=1E-7
#====================================================
弹性矩阵会直接写在OUTCAR里。
在outcar中搜索：TOTAL ELASTIC
计算得到的elastic moduli 是直角坐标下的，只要在同一直角坐标下，不同取lattice，结果相同。测试过
不能设置NPAR， 不然会报错: VASP internal routines  have requested a change of the k-point set.