24小时热门版块排行榜    

Znn3bq.jpeg
查看: 3857  |  回复: 19
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

hbl626

铁虫 (正式写手)

[交流] 用VASP5.3计算bcc-Fe的弹性系数步骤 已有6人参与

最近一直在探索用VASP5.3计算弹性系数,以bcc-Fe为例,发现VASP 手册上面说的很笼统,仅仅用优化好的晶格常数修改IBRION =  6,ISIF   =      3,NFREE = 4后计算得到的弹性系数并不准确。然后在初步优化好的晶格常数的基础上采用下面的INCAR 设置再提高精度进一步优化:
SYSTEM = graphite(100)

Start parameter for this Run:
#   ISTART =      1
  # NPAR   =     2
#  LREAL  = .TRUE.
   ENCUT  =  400
   ISPIN  =   2
   MAGMOM =  2*3
Electronic Relaxation 1
   ISMEAR =  1
   SIGMA  =  0.2
   NELMIN =     2
   NELM   =    80
   EDIFF  =   1E-05   stopping-criterion for ELM
Ionic Relaxation 1
   NSW    =   100
   IBRION =  1
   EDIFFG =    -0.001   stopping-criterion for ELM
  ISIF   =      3
DOS related values:
  LORBIT = 10
  #NFREE = 4
# ISYM=0
  POTIM=0.015
PREC=Accurate
这时候优化好的CONTCAR 为:
Fe bulk
   1.00000000000000
     2.8306224245125096    0.0000000000000000    0.0000000000000000
     0.0000000000000000    2.8306224245125096    0.0000000000000000
     0.0000000000000000    0.0000000000000000    2.8306224245125096
   Fe
     2
Selective dynamics
Direct
  0.0000000000000000  0.0000000000000000  0.0000000000000000   T   T   T
  0.5000000000000000  0.5000000000000000  0.5000000000000000   T   T   T
这时候将CONTCAR 拷贝到POSCAR中修改INCAR中参数
   NSW    =   1
   IBRION =  6
   去掉NFREE前的#
重新提交(删去WAVCAR)
最后在OUTCAR里面可以得到弹性系数:
TOTAL ELASTIC MODULI (kBar)
Direction    XX          YY          ZZ          XY          YZ          ZX
--------------------------------------------------------------------------------
XX        2666.2622   1478.7279   1478.7279      0.0000      0.0000      0.0000
YY        1478.7279   2666.2622   1478.7279      0.0000      0.0000      0.0000
ZZ        1478.7279   1478.7279   2666.2622      0.0000      0.0000      0.0000
XY           0.0000      0.0000      0.0000    985.3189      0.0000      0.0000
YZ           0.0000      0.0000      0.0000      0.0000    985.3189      0.0000
ZX           0.0000      0.0000      0.0000      0.0000      0.0000    985.3189
--------------------------------------------------------------------------------
C_11=267 GPa  C_21=148 GPa  C_44=99 GPa.
回复此楼

» 猜你喜欢

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

dengc09

禁虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本帖内容被屏蔽

7楼2014-06-09 15:49:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 20 个回答

dengc09

禁虫 (正式写手)

本帖内容被屏蔽

2楼2014-06-06 14:31:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

675932406

铜虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主你能分享下5.3的源文件吗,
3楼2014-06-06 15:03:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小静

铜虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
这样算出来的就比较准确吗?我之前也是用IBRION =  6,ISIF = 3,NFREE = 4来计算弹性常数,计算的结果和MS计算出来的比较吻合,但是和实验室相差甚远啊,不知道这是为什么?请问楼主你觉得这会是什么原因呢?
让短暂的生命绽放最绚烂的光彩。
4楼2014-06-06 15:22:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 297,工科调剂? +4 河南农业大学-能 2026-04-14 4/200 2026-04-16 22:52 by wulijun2012
[考研] 291求调剂 +9 关忆北. 2026-04-14 9/450 2026-04-16 22:49 by cfdbai
[考研] 22专硕求调剂 +10 haoyun上岸 2026-04-11 12/600 2026-04-16 22:21 by 猪会飞
[考研] 材料相关专业344求调剂双非工科学校或课题组 +23 hualkop 2026-04-12 25/1250 2026-04-16 22:12 by SUSE_CL
[考研] 22408 312求调剂 +23 门路摸摸 2026-04-14 25/1250 2026-04-16 21:21 by Art1977
[考研] 294求调剂 +14 淡然654321 2026-04-15 14/700 2026-04-16 21:01 by lpl364211
[考研] 297,工科调剂?河南农业大学本科 +14 河南农业大学-能 2026-04-14 14/700 2026-04-16 14:41 by dingyanbo1
[考研] 290调剂生物0860 +38 哇哈哈,。 2026-04-11 44/2200 2026-04-16 09:52 by cuisz
[考研] 297,工科调剂? +10 河南农业大学-能 2026-04-14 10/500 2026-04-15 21:50 by noqvsozv
[考研] 366求调剂 +11 不知名的小卅 2026-04-11 11/550 2026-04-14 15:50 by zs92450
[考研] 105500药学求调剂 +4 x_skys 2026-04-12 4/200 2026-04-14 13:37 by rndfc
[考研] 食品与营养(0955)271求调剂 +15 升格阿达 2026-04-12 16/800 2026-04-14 13:18 by 浮若_安生
[考研] 085408光电信息工程专硕355一志愿长春光机所调剂 +6 王ymaa 2026-04-13 13/650 2026-04-14 11:33 by 王ymaa
[考研] 085600材料与化工349分求调剂 +16 李木子啊哈哈 2026-04-12 17/850 2026-04-14 09:11 by fenglj492
[考研] 考研英一数一338分 +9 长江大学东校区 2026-04-13 10/500 2026-04-14 00:41 by 王珺璞
[考研] 297工科,求调剂? +13 河南农业大学-能 2026-04-12 13/650 2026-04-13 14:12 by dingyanbo1
[考研] 290求调剂 +18 柯淮然 2026-04-12 20/1000 2026-04-13 12:56 by cyh—315
[考研] 一志愿085802 323分求调剂 +13 drizzle_9 2026-04-12 14/700 2026-04-13 10:26 by Faiz5552
[考研] 0831一轮调剂失败求助 +10 小熊睿睿_s 2026-04-11 10/500 2026-04-12 22:43 by 长弓傲
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
信息提示
请填处理意见