| 查看: 1714 | 回复: 3 | ||
[求助]
TiO2 晶胞参数优化问题
|
|
最近在这里卡了好久 还请各位前辈指点一下看看是哪里出了问题。 （用的原胞如图所示） INCAR TiO2-unit-cell ENCUT=450 ISTART=0;ICHARG=2 ISMEAR=0;SIGMA=0.2 NSW=60;IBRION=2 ISIF=3 POTIM=0.2 EDIFF=1E-6;EDIFFG=-1E-3 PREC=Accurate NPAR=4 KPOINTS Automatic generation 0 Monhkorst-Pack 9 9 9 0.0 0.0 0.0 脚本 #!/bin/sh rm WAVECAR for i in 9.436 9.446 9.456 9.466 9.476 9.486 9.496 9.506 9.516 9.526 9.536 do cat > POSCAR << ! SYSTEM=TiO2-unit $i 0.3980602994 0.0000000000 0.0000000000 0.0000000000 0.3980602994 0.0000000000 0.0000000000 0.0000000000 1.0000000000 4 8 Direct 0.0000000000 0.0000000000 0.0000000000 0.5000000000 0.5000000000 0.5000000000 0.0000000000 0.5000000000 0.2500000000 0.5000000000 0.0000000000 0.7500000000 0.0000000000 0.0000000000 0.2080000000 0.5000000000 0.5000000000 0.7080000000 0.0000000000 0.5000000000 0.0420000000 0.5000000000 0.0000000000 0.5420000000 0.0000000000 0.0000000000 -0.2080000000 0.5000000000 0.5000000000 0.2920000000 0.0000000000 0.5000000000 0.4580000000 0.5000000000 0.0000000000 0.9580000000 ! echo "a=$i angstrom" ; time mpiexec -np 16 vasp E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f\n",$5}'` V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.4f\n",$5}'` echo $E $V >> EtVo.dat echo $i : $V $E >> SUMMARY cp OUTCAR OUTCAR.$i cp CONTCAR CONTCAR.$i done 计算完成后得到的不是二次曲线 而是锯齿状的。如图所示。不知道这个问题出在哪里。 还有就是用二步法优化的话出来的结果也不是二次曲线，而是斜率为负的直线。也请各位前辈帮忙看看问题出在哪里 INCAR.relax TiO2-unit ENCUT=450 ISTART=0;ICHARG=2 ISMEAR=0;SIGMA=0.02 NSW=60;IBRION=2 ISIF=4 POTIM=0.2 EDIFF=1E-5;EDIFFG=-1E-3 PREC=Accurate NPAR=4 INCAR.static TiO2-unit ENCUT=450 ISTART=0;ICHARG=2 ISMEAR=-5 NPAR=4 PREC=Accurate KPOINTS: Automatic generation 0 Monhkorst-Pack 9 9 9 0.0 0.0 0.0 脚本 #!/bin/sh rm WAVECAR for i in 循环值 do cat > POSCAR << ! SYSTEM=TiO2-unit $i 0.3980602994 0.0000000000 0.0000000000 0.0000000000 0.3980602994 0.0000000000 0.0000000000 0.0000000000 1.0000000000 4 8 Direct 0.0000000000 0.0000000000 0.0000000000 0.5000000000 0.5000000000 0.5000000000 0.0000000000 0.5000000000 0.2500000000 0.5000000000 0.0000000000 0.7500000000 0.0000000000 0.0000000000 0.2080000000 0.5000000000 0.5000000000 0.7080000000 0.0000000000 0.5000000000 0.0420000000 0.5000000000 0.0000000000 0.5420000000 0.0000000000 0.0000000000 -0.208000000 0.5000000000 0.5000000000 0.2920000000 0.0000000000 0.5000000000 0.4580000000 0.5000000000 0.0000000000 0.9580000000 ! cp INCAR.relax INCAR echo "relax a=$i angstrom" ; time mpiexec -np 16 vasp >/dev/null cp CONTCAR POSCAR cp INCAR.static INCAR echo "static a=$i angstrom" ; time mpiexec -np 16 vasp >/dev/null E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f\n",$5}'` V=`grep "volume" OUTCAR | tail -1 | awk '{printf "%12.4f\n",$5}'` echo $i : $V $E >> SUMMARY cp OUTCAR OUTCAR.$i done  optonce.JPG  originalcell.JPG |
回复此楼» 猜你喜欢
324求调剂
已经有12人回复
-
311求调剂
已经有7人回复
-
085404，285分求调剂
已经有3人回复
-
275 求调剂
已经有8人回复
-
调剂
已经有3人回复
-
生物与医药273求调剂
已经有16人回复
-
求调剂
已经有13人回复
-
327求调剂
已经有11人回复
-
274求调剂求调剂
已经有12人回复
-
0703化学
已经有26人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- MS建无定形晶胞
已经有21人回复
- 关于vasp超胞中LREAL如何设置
已经有7人回复
- 请教如何用VASP计算氧化物的晶格常数？以alpha-Al2O3为例
已经有7人回复
- vasp 计算单分子/原子 时晶胞大小的选取
已经有6人回复
- MS 结构优化的结果不能理解 求高人指点
已经有10人回复
- ZnO原胞原子数
已经有8人回复
- CASTEP计算出错，好像是并行出了问题，请大家帮我看一下
已经有7人回复
- 两个晶格常数值相差多少（百分比）可以认为相等或不等？
已经有7人回复
- 关于晶格常数和DOS
已经有8人回复
- VASP并行可执行软件包，可对晶胞参数进行部分优化
已经有17人回复
- 请问：SCLAE S=2，对于三维体系，意思是建一个2*2*2的超晶胞吗？
已经有9人回复
- TiO2（101）晶面的优化
已经有13人回复
- VASP一次性优化锐钛矿TiO2晶格常数时出错！！
已经有6人回复
- 【原创】EMASS是CPMD最重要的参数
已经有20人回复
- 【求助成功】关于ZnO理想晶胞建模的问题
已经有20人回复
- 【求助】勾选与不选optimize cell 的差异
已经有15人回复
- 【求助】LDA+U优化问题
已经有10人回复
- 【求助】在计算红外时k－point的选择及晶胞是否优化？
已经有5人回复
- 【求助】如何解决自旋极化计算时的难以收敛问题？
已经有16人回复
saikun
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 2483.2
- 散金: 891
- 红花: 5
- 帖子: 424
- 在线: 89.1小时
- 虫号: 197920
- 注册: 2006-02-25
- 性别: GG
- 专业: 古陶瓷与传统陶瓷
2楼2013-08-29 14:02:32
3楼2013-08-30 10:46:07
saikun
木虫 (正式写手)
- 应助: 4 (幼儿园)
- 金币: 2483.2
- 散金: 891
- 红花: 5
- 帖子: 424
- 在线: 89.1小时
- 虫号: 197920
- 注册: 2006-02-25
- 性别: GG
- 专业: 古陶瓷与传统陶瓷
4楼2013-08-30 13:06:41
100