| 查看: 588 | 回复: 0 | |||
qphll金虫 (正式写手)
|
[求助]
收敛测试, LaFeO3, 能量变化异常
|
|
尝试用vasp计算LaFeO3, 在收敛测试的时候发现随着encut的增加, 总能反而升高. (1) INCAR脚本 m WAVECAR for i in 200 250 300 350 400 450 500 550 600 do cat > INCAR << ! SYSTEM = LF orthorhombic ENCUT = $i ISMEAR = 0 ! Gaussian smearing SIGMA = 0.1 LREAL= Auto PREC= Accurate ! echo "ENCUT= $i eV" ; time; ./vasp5211 ; E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'` echo $i $E >> compare done (2) KPOINTS 文件 Automatic generation 0 Monhkorst-Pack 4 4 4 0.0 0.0 0.0 (3) POSCAR 文件 rthorhombic 1.0 5.5349998474 0.0000000000 0.0000000000 0.0000000000 7.8880000114 0.0000000000 0.0000000000 0.0000000000 5.5989999771 12 4 4 ! O Fe La Direct 0.270900011 0.034699999 0.229399994 0.729099989 0.965300024 0.770600021 0.229099989 0.965300024 0.729399979 0.770900011 0.034699999 0.270600021 0.729099989 0.534699976 0.770600021 0.270900011 0.465299994 0.229399994 0.770900011 0.465299994 0.270600021 0.229099989 0.534699976 0.729399979 0.493999988 0.250000000 0.564400017 0.506000042 0.750000000 0.435599983 0.006000012 0.750000000 0.064400017 0.993999958 0.250000000 0.935599983 0.000000000 0.000000000 0.000000000 0.500000000 0.000000000 0.500000000 0.000000000 0.500000000 0.000000000 0.500000000 0.500000000 0.500000000 0.518400013 0.250000000 0.008900000 0.481599987 0.750000000 0.991100013 0.981599987 0.750000000 0.508899987 0.018400013 0.250000000 0.491100013 最后得到的图是这样的: http://edu.muchong.com/attachment/16/b3/89654_1333900508_628.jpg 同样, 如果修改脚本, 对sigma测试, 最后得到的图中能量是线性的: SIGMA 脚本 rm WAVECAR for i in 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28 0.30 do cat > INCAR << ! SYSTEM = LF orthorhombic ENCUT = 500 ISMEAR = 0 ! Gaussian smearing SIGMA = $i LREAL= Auto PREC= Accurate ISTART = 0 ICHARGE = 2 ! echo "SIGMA= $i eV" ; time; ./vasp5211 ; E=`grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'` echo $i $E >> sigma.compare done sigma测试能量图: http://edu.muchong.com/attachment/87/87/89654_1333900718_282.jpg 这是怎么回事呢? 可能的问题: a. INCAR中的参数设置?? b.赝势的选择?? PAW-PBE的O, La, Fe有几个, 我倒没有进行选择 c. 先优化结构, 再单点能计算?? POSCAR直接从文献报导的CIF得到. 我的打算是做完测试以后, 再做晶体的结构优化. 大家有什么建议? 多谢.  ENCUT.test  SIGMA.test [ Last edited by qphll on 2012-4-9 at 01:03 ] |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有10人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有7人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有8人回复
-
有多少人是今天查系统知道结果的？
已经有18人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有11人回复
找到一些相关的精华帖子，希望有用哦~
求助关于荧光共振能量转移的问题？
已经有8人回复
- 用Dmol优化收敛的结构，用vasp不收敛？
已经有9人回复
- 3.22文献解读：三元材料在高能量密度锂离子电池中的应用
已经有236人回复
- 干涉光能量哪里去了？
已经有5人回复
- 求助：castep计算参数设置没有收敛
已经有9人回复
- 用CHEM3D软件PM3理论水平计算分子轨道却无法显示轨道能量
已经有5人回复
- 关于dmol3的收敛性测试
已经有12人回复
- 天然气燃烧的能量利用效率应该怎么计算？
已经有4人回复
- 请教扣电测试异常
已经有4人回复
- 两个原子靠近成键的过程中，系统总能量会不会发生变化？
已经有3人回复
- XPS中能量如何校准
已经有7人回复
- 计算表面时：Vacuum真空层和slab层厚度收敛性测试，保持超晶胞体积不变吗？？
已经有10人回复
- 【其他】关于收敛测试的问题
已经有22人回复
- 【求助】收敛测试怎么做？？？
已经有17人回复
- 【求助】收敛性测试的问题
已经有40人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
10