| 查看: 3620 | 回复: 19 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
用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. |
回复此楼» 猜你喜欢
想换工作。大多数高校都是 评职称时 认可5年内在原单位取得的成果吗？
已经有7人回复
-
需要合成515-64-0，50g，能接单的留言
已经有4人回复
-
自荐读博
已经有4人回复
-
上海工程技术大学【激光智能制造】课题组招收硕士
已经有4人回复
-
写了一篇“相变储能技术在冷库中应用”的论文，论文内容以实验为主，投什么期刊合适？
已经有6人回复
-
带资进组求博导收留
已经有10人回复
-
最近几年招的学生写论文不引自己组发的文章
已经有11人回复
-
中科院杭州医学所招收博士生一名（生物分析化学、药物递送）
已经有3人回复
-
临港实验室与上科大联培博士招生1名
已经有8人回复
16楼2014-06-10 09:01:01
|
本帖内容被屏蔽 |
2楼2014-06-06 14:31:43
675932406
铜虫 (小有名气)
- 应助: 1 (幼儿园)
- 金币: 21.6
- 散金: 134
- 帖子: 131
- 在线: 72.4小时
- 虫号: 1146053
- 注册: 2010-11-13
- 专业: 力学中的基本问题和方法
3楼2014-06-06 15:03:30
4楼2014-06-06 15:22:07