应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 用VASP5.3计算bcc-Fe的弹性系数步骤
201/2返回列表12下一页
查看: 3844  |  回复: 19
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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.
回复此楼

» 猜你喜欢
1楼 2014-06-06 14:01:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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

dengc09

禁虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽
5楼2014-06-06 15:54:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小静

铜虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
5楼: Originally posted by dengc09 at 2014-06-06 15:54:17
一般来说，貌似energy-strain方法算出来的准一些，因为这样做直接一些，但stress-strain算出来也不会相差太大吧，你看看计算弹性常数之前结构是不是弛豫好了...

哦  好的 我再精确优化一下试试看
赞一下 回复此楼
让短暂的生命绽放最绚烂的光彩。
6楼2014-06-09 15:46:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dengc09

禁虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
本帖内容被屏蔽
7楼2014-06-09 15:49:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小静

铜虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
7楼: Originally posted by dengc09 at 2014-06-09 15:49:26
看一下我上传的那个文件吧，我给算了一下，就在下面...

请问哪里？没有看到啊
赞一下 回复此楼
让短暂的生命绽放最绚烂的光彩。
8楼2014-06-09 16:37:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dengc09

禁虫 (正式写手)
本帖内容被屏蔽
9楼2014-06-09 16:38:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小小静

铜虫 (小有名气)
引用回帖:
9楼: Originally posted by dengc09 at 2014-06-09 16:38:58
是资源帖

哦 谢谢  我去看看
回复此楼
让短暂的生命绽放最绚烂的光彩。
10楼2014-06-09 16:40:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 hbl626 的主题更新
201/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程085601，270求调剂 +35 @ASDF1234 2026-04-08 37/1850 2026-04-13 01:18 by 幸免 ..
[考研] 一志愿浙大生物325分求调剂 +9 zysheng 2026-04-12 9/450 2026-04-12 22:31 by yuyin1233
[考研] 求助调剂，跨调 +12 X十甫寸Y 2026-04-11 12/600 2026-04-12 10:20 by zhouyuwinner
[考研] 290求调剂 +13 luoziheng 2026-04-10 13/650 2026-04-11 23:15 by labixiaoqiao
[考研] 269求调剂 +11 啊啊我我 2026-04-07 11/550 2026-04-11 16:45 by vgtyfty
[考研] 085600材料与化工329分求调剂 +16 叶zilin 2026-04-10 16/800 2026-04-11 11:04 by may_新宇
[考研] 300分求调剂 （085501机械专硕，本科扬大） +8 xu@841019 2026-04-11 8/400 2026-04-11 10:46 by qingpingzhu
[考研] 083200 305分 求二轮调剂 不接受跨专业 +9 Claireyyyy 2026-04-09 10/500 2026-04-10 21:21 by Claireyyyy
[考研] 0858求调剂 5+5 Gky09300550， 2026-04-10 8/400 2026-04-10 19:13 by chemisry
[论文投稿] mdpi小修rvr时间四五天了 20+3 哈哈high 2026-04-08 5/250 2026-04-10 16:02 by 北京莱茵润色
[考研] 一志愿西交机械专硕求调剂 +8 求上岸的小王 2026-04-10 8/400 2026-04-10 15:09 by hemengdong
[考研] 一志愿中国科学院上海有机所，有机化学356分找调剂 +11 Nadiums 2026-04-09 11/550 2026-04-09 18:04 by lijunpoly
[考研] 283电子信息求调剂 +4 三石WL 2026-04-08 4/200 2026-04-09 10:21 by wp06
[考研] 招收有机化学、化工，药学，食品灯专业学生 +3 yrfhjgdj 2026-04-08 3/150 2026-04-09 10:15 by QYQX_123
[考研] 考研求调剂 +4 雯??? 2026-04-08 4/200 2026-04-08 21:44 by 土木硕士招生
[考研] 22408 266求调剂 +11 masss11222 2026-04-07 14/700 2026-04-08 11:06 by yulian1987
[考研] 313求调剂 +3 十六拾陆 2026-04-07 3/150 2026-04-07 23:20 by lbsjt
[考研] 一志愿南科大生物学297分，求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 325 调剂 +6 QQ小虾 2026-04-07 6/300 2026-04-07 15:17 by Ccclqqq
[考研] 信工所11408 340分 本科西安交大自动化 +3 moontrek 2026-04-06 3/150 2026-04-07 09:56 by chongya
信息提示
关闭
请填处理意见
关闭