VASP5.2直接计算弹性常数
弹性常数[elastic constant]是经常会用到的参数。
在VASP5.2之前,计算弹性常数是非常艰巨的一项工作,不信你可以阅读候柱峰博士写的文档:Use VASP to Calculate Elastic Constants [PDF下载] 。
但现在有了5.2之后,算弹性常数就变成非常容易的一件事了,具体介绍如下。
INCAR中主要是配置3个参数:将IBRION=6;ISIF=3;NFREE=4 ,INCAR源文件如下:
SYSTEM = diamond PREC = Accurate ENCUT = 600.0 EDIFF = 1e-6 ISMEAR = 0 SIGMA = 0.05 POTIM = 0.100 LCHARG = FALSE LWAVE = FALSE IBRION=6 ISIF=3 NFREE=4 EDIFFG = 1e-5 EDIFFG = -0.001
弹性常数就在OUTCAR里,如下。
TOTAL ELASTIC MODULI (kBar) Direction XX YY ZZ XY YZ ZX -------------------------------------------------------------------------------- XX 10542.0731 1253.9340 1253.8826 -33.2003 0.0000 0.0000 YY 1253.9340 10542.0731 1253.8826 -33.2003 0.0000 0.0000 ZZ 1253.8826 1253.8826 10542.0889 -4.3655 0.0000 0.0000 XY -33.2003 -33.2003 -4.3655 5633.1588 0.0000 0.0000 YZ 0.0000 0.0000 0.0000 0.0000 5633.9957 -11.1251 ZX 0.0000 0.0000 0.0000 0.0000 -11.1251 5633.9957
上面的数据,XX, YY, ZZ, XY, YZ, ZX 分别对应弹性常数中的1-6。比如你找C44,就是5633.1588kBar,合563.31588GPa。
参考:小木虫 http://emuch.net/html/201205/4533931.html
仅有一条评论
admin好,
这个方法试了一下,ISMEAR=0计算NaCl结构的话,c44会出现负值;设置ISMEAR=-5的话,对NaCl的计算结果,会比较准确。