24小时热门版块排行榜

返回列表

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

董大瓜

新虫 (小有名气)

应助: 6 (幼儿园)
金币: 3194.1
红花: 6
帖子: 158
在线: 57.9小时
虫号: 2326901
注册: 2013-03-07
专业: 凝聚态物性 II ：电子结构

[求助] VASP计算ZnO的性质，为什么ENCUT的测试曲线是这个样子的？

我初学VASP，好多问题不懂。现在想计算ZnO的性质，顺便学习以下，因此决定先测试ENCUT。下面是我的输入文件和运行脚本：
---------------------------------------------------------------------------------------------------
POSCAR（MS直接导出的原胞）
SYSTEM = ZnO
1
2.813950363754660    -1.624635000000000    0.000000000000000
0.000000000000000    3.249270000000000    0.000000000000000
0.000000000000000    0.000000000000000    5.205440000000000
O Zn
2 2
Slective Dynamics
Direct
0.3333333333333330 0.6666666666666670 0.3825999999999999 F F F
-0.3333333333333330  -0.6666666666666670 0.8825999999999999 F F F
0.3333333333333330 0.6666666666666670 0.0000000000000000 F F F
-0.3333333333333330  -0.6666666666666670 0.5000000000000000 F F F
-------------------------------------------------------------------------------------------------
KPONITS
------------
Automatic Generation
0
Monhkorst-Pack
13 13 9
0 0 0
--------------------------------------------------------------------------------------------------
POTCAR：
-------------
用的是paw_GGA
grep "ENMAX" POTCAR 后，显示如下：
ENMAX  =  400.000; ENMIN  =  300.000 eV
ENMAX  =  276.749; ENMIN  =  207.562 eV
--------------------------------------------------------------------------------------------------
运行脚本（run.ENCUT.sh）
---------------------------------
#!/bin/sh
rm Encut_Energy.dat
for i in 400 450 500 550 600 650 700 750 800 850 900 950 1000
do
rm WAVECAR
cat > INCAR <<!
SYSTEM = ZnO
ENCUT = $i
ISTART = 0
ICHARG = 2
ISMEAR = -5
PREC = Accurate
!
echo "ENCUT = $i eV" ; time mpirun -np 8 vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f\n",$5}'`
echo $i $E >>Encut_Energy.dat
done
cat Encut_Energy.dat
-------------------------------------------------------------------------------------------------
最后计算出来的结果是这个样子的：
400 -18.078393
450 -18.047978
500 -18.039628
550 -18.040343
600 -18.043653
650 -18.046893
700 -18.049904
750 -18.051407
800 -18.052394
850 -18.052619
900 -18.052538
950 -18.052588
1000 -18.052498
----------------------
MATLAB做图后是这个样子的：
见附件。
--------------------------------------------------------------------------------------------------
为什么我的测试曲线会先上升，然后才趋于平缓呢？400eV是赝势库里的值，而按照这个测试的结果应该取800eV，但是这是的能量值比400eV时的要高很多，在接下来的计算中我该怎么选ENCUT 的值呢，400 or 800eV？
或者说我的计算本身就有问题？
求帮助

ENCUT_Energy.jpg

回复此楼