返回首页

VASP5.2直接计算弹性常数

时间:2012-11-09 00:38来源:知行网www.zhixing123.cn 编辑:麦田守望者

弹性常数[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。

------分隔线----------------------------
标签(Tag):MATLAB Matlab技巧 Matlab实例教程 matlab源代码 matlab基础教程
------分隔线----------------------------
推荐内容
猜你感兴趣