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

标签: vasp

相关文章推荐

添加新评论 (无需注册,可直接评论)

仅有一条评论

  1. lammps_vasp

    admin好,
    这个方法试了一下,ISMEAR=0计算NaCl结构的话,c44会出现负值;设置ISMEAR=-5的话,对NaCl的计算结果,会比较准确。