24小时热门版块排行榜    

CyRhmU.jpeg
南方科技大学公共卫生及应急管理学院2025级博士研究生招生报考通知
查看: 1684  |  回复: 6
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

janet2012222

金虫 (小有名气)

[求助] VASP直接计算弹性常数出错

大家好:我用VASP计算金刚石的弹性常数。INCAR设置如下:
  SYSTEM=diamond  
    #NPAR=4
    ISTART=0
    ICHARG=2
    ENCUT=600
    ISMEAR=0
    SIGMA=0.1
    ISIF=3
    IBRION=6
     NFREE=4
    NSW=100
    EDIFF=0.1E-04
    PREC=Accurate
    POTIM=0.05
    ISPIN=1
    LWAVE=.FALSE.
    LCHARG=.FALSE.
    LREAL=Auto
可是得到的结果却是:
TOTAL ELASTIC MODULI (kBar)
Direction    XX          YY          ZZ          XY          YZ          ZX
--------------------------------------------------------------------------------
XX       11430.2418    911.0905    599.6742      0.0000    440.4091      0.0000
YY         911.0905  11430.2418    599.6742      0.0000   -440.4091     -0.0000
ZZ         599.6742    599.6742  11741.6580      0.0000      0.0000      0.0000
XY           0.0000      0.0000      0.0000   5259.5757      0.0000    440.4091
YZ         440.4091   -440.4091      0.0000      0.0000   4948.1594      0.0000
ZX           0.0000     -0.0000      0.0000    440.4091      0.0000   4948.1594
--------------------------------------------------------------------------------
很明显不对,我不知道哪里设置得有问题。请大家帮我看看
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cindywei

木虫 (初入文坛)

ISMEAR = -5试一下吧。VASP直接计算弹性常数的原理是微小变形后静态计算的能量与形变的关系。
6楼2014-03-19 11:15:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 7 个回答

hn19870519

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
liliangfang: 金币+2, 谢谢交流 2013-06-29 14:37:39
用vasp算弹性常数,我见过多的都是根据晶格对称性,推出弹性常数表达式(应变为参数),然后再用vasp算出每个应变下的能量,再去拟合。具体参照候老师的帖子。还有jpcm的文章 ,15,5307(2003)。
2楼2013-06-27 11:19:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

janet2012222

金虫 (小有名气)

引用回帖:
2楼: Originally posted by hn19870519 at 2013-06-27 11:19:09
用vasp算弹性常数,我见过多的都是根据晶格对称性,推出弹性常数表达式(应变为参数),然后再用vasp算出每个应变下的能量,再去拟合。具体参照候老师的帖子。还有jpcm的文章 ,15,5307(2003)。

那个好像很麻烦,不是说5.2以上的版本可以直接计算出弹性常数吗?这是我用原胞计算的结果。我刚刚用金刚石晶胞计算的结果就对了。
3楼2013-06-27 11:23:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

janet2012222

金虫 (小有名气)

引用回帖:
2楼: Originally posted by hn19870519 at 2013-06-27 11:19:09
用vasp算弹性常数,我见过多的都是根据晶格对称性,推出弹性常数表达式(应变为参数),然后再用vasp算出每个应变下的能量,再去拟合。具体参照候老师的帖子。还有jpcm的文章 ,15,5307(2003)。

感觉只能计算正交的结构似的
4楼2013-06-27 11:24:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见