版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (4387)
>文献求助 (326)
>虫友互识 (255)
>基金申请 (239)
>导师招生 (193)
>硕博家园 (144)
>休闲灌水 (139)
>论文投稿 (131)
>论文道贺祈福 (94)
>考博 (93)
>教师之家 (91)
>考研 (81)
>博后之家 (70)
>找工作 (49)
>公派出国 (48)
>绿色求助(高悬赏) (40)

返回列表

hbl626

铁虫 (正式写手)

应助: 0 (幼儿园)
金币: 19.8
散金: 19
红花: 2
帖子: 330
在线: 54.2小时
虫号: 546036
注册: 2008-04-15
性别: MM
专业: 金属结构材料

[交流] 用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.

回复此楼

» 猜你喜欢

微信指数没变化，科研之友没阅读已经有18人回复
基础研究怎么拉横向，学校到款任务越来越多，难以完成拉横向，都有哪些途径啊已经有8人回复
HXDI做水性聚氨酯乳液，是不是特别容易出渣已经有3人回复
小木虫看见有人已经知道结果了已经有19人回复
2026年国自然面上资助率已经有8人回复
系统今天又提示维护了，估计离放榜不远了已经有15人回复
申请2027年材料/化学类博士已经有5人回复
你们的时间戳变了吗已经有4人回复
今年的WR进展到哪一步了？已经有9人回复
产物和副产物价值比较已经有5人回复

1楼 2014-06-06 14:01:40

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

dengc09

禁虫 (正式写手)

本帖内容被屏蔽

2楼2014-06-06 14:31:43

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

675932406

铜虫 (小有名气)

应助: 1 (幼儿园)
金币: 21.6
散金: 134
帖子: 131
在线: 72.4小时
虫号: 1146053
注册: 2010-11-13
专业: 力学中的基本问题和方法

★
小木虫: 金币+0.5, 给个红包，谢谢回帖

楼主你能分享下5.3的源文件吗，

赞一下

回复此楼

3楼2014-06-06 15:03:30

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

小小静

铜虫 (小有名气)

应助: 0 (幼儿园)
金币: 1474.5
帖子: 126
在线: 285.5小时
虫号: 2090541
注册: 2012-10-27
性别: MM
专业: 理论和计算化学

★
小木虫: 金币+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 (幼儿园)
金币: 1474.5
帖子: 126
在线: 285.5小时
虫号: 2090541
注册: 2012-10-27
性别: MM
专业: 理论和计算化学

★
小木虫: 金币+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的回帖