24小时热门版块排行榜    

查看: 2339  |  回复: 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的回帖

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的回帖

herohou

金虫 (小有名气)

楼主的问题搞定了么?直接算弹性常数可行么?
5楼2014-03-17 21:35:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cindywei

木虫 (初入文坛)

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

JJ20

金虫 (小有名气)

你把nsw 去掉,POTIM设到0.1试试,应该是可以的
选择放弃只为保留尊严
7楼2014-03-20 09:26:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 janet2012222 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 11:14 by xSrQNSSjcgoX
[考研] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 09:35 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 SX5mJmPWDfOO 2026-11-13 4/200 2026-11-15 07:34 by xSrQNSSjcgoX
[基金申请] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-13 3/150 2026-11-15 04:16 by xSrQNSSjcgoX
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 6/300 2026-11-15 02:59 by xSrQNSSjcgoX
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 7/350 2026-11-15 02:48 by xSrQNSSjcgoX
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 1Nb0IVlmV6v0 2026-11-11 4/200 2026-11-13 11:18 by Equinoxhua
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 FChKX8meu1ds 2026-11-09 10/500 2026-11-13 11:08 by Equinoxhua
[考博] 27博士申请 +5 白白白ly 2026-11-09 5/250 2026-11-13 05:01 by wlf1982
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 4/200 2026-11-12 12:51 by ECLTXnJnRSXf
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 TgvkbZZl9tqB 2026-11-11 6/300 2026-11-12 12:42 by ECLTXnJnRSXf
[硕博家园] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +4 cdqGD6ftHmbl 2026-11-11 4/200 2026-11-12 12:26 by ECLTXnJnRSXf
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:33 by HajghKOMynmz
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 Z7qfCsdIiD1T 2026-11-10 7/350 2026-11-12 11:31 by HajghKOMynmz
[考研] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +7 czAiolzR0IJC 2026-11-10 9/450 2026-11-12 11:00 by HajghKOMynmz
[硕博家园] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 C6CRn66RYYXS 2026-11-10 6/300 2026-11-12 10:16 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 FChKX8meu1ds 2026-11-09 8/400 2026-11-12 09:46 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 12/600 2026-11-12 09:21 by HajghKOMynmz
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 TgvkbZZl9tqB 2026-11-11 3/150 2026-11-12 08:26 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +4 FChKX8meu1ds 2026-11-09 7/350 2026-11-12 05:19 by 9fOTqj074tIg
信息提示
请填处理意见