| 查看: 574 | 回复: 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 ] |
回复此楼» 猜你喜欢
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
-
filecode与中标关系的预测
已经有3人回复
-
filecode变化情况
已经有22人回复
-
十年后又回来了，论文投稿求助
已经有4人回复
找到一些相关的精华帖子，希望有用哦~
求助关于荧光共振能量转移的问题？
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
100