加载中…
个人资料
  • 博客等级:
  • 博客积分:
  • 博客访问:
  • 关注人气:
  • 获赠金笔:0支
  • 赠出金笔:0支
  • 荣誉徽章:
正文 字体大小:

[转载]vasp计算弹性常数

(2016-06-13 19:53:24)
标签:

转载

分类: vasp
原文地址:vasp计算弹性常数作者:小龙

vasp5.2计算可以直接计算弹性常数。

关键参数:  将IBRION=6;ISIF=3;NFREE=4;

弹性矩阵会直接写在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       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

0

  

新浪BLOG意见反馈留言板 欢迎批评指正

新浪简介 | About Sina | 广告服务 | 联系我们 | 招聘信息 | 网站律师 | SINA English | 产品答疑

新浪公司 版权所有