24小时热门版块排行榜    

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

herohou

金虫 (小有名气)

楼主的问题搞定了么?直接算弹性常数可行么?
5楼2014-03-17 21:35:28
已阅   回复此楼   关注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的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 11:25 by xSrQNSSjcgoX
[博后之家] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 TZqyNTBa07NP 2026-11-13 3/150 2026-11-15 10:19 by xSrQNSSjcgoX
[考博] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 TZqyNTBa07NP 2026-11-13 4/200 2026-11-15 09:57 by xSrQNSSjcgoX
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 SX5mJmPWDfOO 2026-11-13 3/150 2026-11-15 08:07 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 05:00 by xSrQNSSjcgoX
[硕博家园] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 uKOSCocCb6Fm 2026-11-12 6/300 2026-11-15 02:59 by xSrQNSSjcgoX
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +7 0TEMgxdT1eIw 2026-11-09 10/500 2026-11-15 00:39 by Equinoxhua
[公派出国] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +4 G2rhSvqvhNzl 2026-11-10 5/250 2026-11-13 11:37 by Equinoxhua
[教师之家] 售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,我:8O.551.O54,科目全,可十急 +8 zSbvVX9JFl9J 2026-11-09 9/450 2026-11-13 11:11 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
[考博] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 Z7qfCsdIiD1T 2026-11-10 5/250 2026-11-12 11:25 by HajghKOMynmz
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 czAiolzR0IJC 2026-11-10 7/350 2026-11-12 11:09 by HajghKOMynmz
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +6 YHiw1gou1Eda 2026-11-10 6/300 2026-11-12 10:08 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文章,我:8O.5.5.1O.54,科目全,可十急 +7 0TEMgxdT1eIw 2026-11-09 9/450 2026-11-12 09:24 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文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 C6CRn66RYYXS 2026-11-10 5/250 2026-11-12 06:03 by 9fOTqj074tIg
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +5 zSbvVX9JFl9J 2026-11-09 8/400 2026-11-12 04:57 by 9fOTqj074tIg
信息提示
请填处理意见