24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

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.

回复此楼