[转载]vasp计算弹性常数
(2016-06-13 19:53:24)
标签:
转载 |
分类: vasp |
原文地址:vasp计算弹性常数作者:小龙
vasp5.2计算可以直接计算弹性常数。
关键参数:
弹性矩阵会直接写在OUTCAR里。
弹性常数[elastic constant]是经常会用到的参数。
在VASP5.2之前,计算弹性常数是非常艰巨的一项工作,不信你可以阅读候柱峰博士写的文档:Use VASP to Calculate Elastic Constants [PDF下载] 。
但现在有了5.2之后,算弹性常数就变成非常容易的一件事了,具体介绍如下。
INCAR中主要是配置3个参数:将IBRION=6;ISIF=3;NFREE=4 ,INCAR源文件如下:
01
02
03
04
05
06
07
08
09
10
11
12
13
14
|
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里,如下。
01
02
03
04
05
06
07
08
09
|
TOTAL ELASTIC MODULI (kBar) Direction -------------------------------------------------------------------------------- XX YY ZZ XY YZ ZX |
上面的数据,XX, YY, ZZ, XY, YZ, ZX 分别对应弹性常数中的1-6。比如你找C44,就是5633.1588kBar,合563.31588GPa。