| 查看: 493 | 回复: 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 ] |
回复此楼» 猜你喜欢
本人考085602 化学工程 专硕
已经有23人回复
-
求调剂院校信息
已经有4人回复
-
085600材料与化工306
已经有4人回复
-
286求调剂
已经有10人回复
-
328求调剂，英语六级551，有科研经历
已经有9人回复
-
一志愿北京化工大学070300 学硕336求调剂
已经有4人回复
-
286分人工智能专业请求调剂愿意跨考！
已经有8人回复
-
资源与环境 调剂申请(333分)
已经有5人回复
-
280求调剂
已经有12人回复
-
269专硕求调剂
已经有5人回复
找到一些相关的精华帖子，希望有用哦~
求助关于荧光共振能量转移的问题？
已经有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人回复
点击这里搜索更多相关资源科研从小木虫开始，人人为我，我为人人
40