| 查看: 1408 | 回复: 10 | ||
董大瓜新虫 (小有名气)
|
[求助]
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 |
回复此楼» 猜你喜欢
参与限项
已经有3人回复
-
假如你的研究生提出不合理要求
已经有7人回复
-
实验室接单子
已经有4人回复
-
全日制（定向）博士
已经有4人回复
-
对氯苯硼酸纯化
已经有3人回复
-
求助:我三月中下旬出站，青基依托单位怎么办？
已经有12人回复
-
不自信的我
已经有12人回复
-
所感
已经有4人回复
-
要不要辞职读博？
已经有7人回复
-
北核录用
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 用VASP优化结构收敛不了，求指教
已经有10人回复
- <求助>VASP 结构优化收敛，静态计算不收敛。（INCAR KPOINTS POSCAR附上）
已经有4人回复
- 【专家在线答疑贴之三】杂化泛函计算FAQ
已经有28人回复
- 关于用超胞方法计算声子的一些经验
已经有21人回复
- vasp 杂化泛函HSE 它怎么就不run涅？
已经有12人回复
- vasp优化结构中出现的问题
已经有10人回复
- VASP优化结构遇到的问题。
已经有9人回复
- 【求助】vasp 5.2中用HSE06计算不能进入主循环
已经有7人回复
- 【求助】求VASP计算ZnO能带和态密度的详细输入设置
已经有15人回复
- 【求助】关于杂化泛函HSE计算
已经有20人回复
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2013-08-09 17:03:42
董大瓜
新虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 3194.1
- 红花: 6
- 帖子: 158
- 在线: 57.9小时
- 虫号: 2326901
- 注册: 2013-03-07
- 专业: 凝聚态物性 II ：电子结构
3楼2013-08-09 18:55:39
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
4楼2013-08-09 19:29:27
董大瓜
新虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 3194.1
- 红花: 6
- 帖子: 158
- 在线: 57.9小时
- 虫号: 2326901
- 注册: 2013-03-07
- 专业: 凝聚态物性 II ：电子结构
5楼2013-08-09 19:35:26
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
老板说，如果我自己出钱把组里机器升下级，取多大都行。6楼2013-08-09 21:00:09
liqizuiyang
木虫 (著名写手)
- 应助: 59 (初中生)
- 金币: 6859.5
- 散金: 261
- 红花: 94
- 帖子: 2093
- 在线: 708.9小时
- 虫号: 2366686
- 注册: 2013-03-21
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
7楼2013-08-09 21:00:50
董大瓜
新虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 3194.1
- 红花: 6
- 帖子: 158
- 在线: 57.9小时
- 虫号: 2326901
- 注册: 2013-03-07
- 专业: 凝聚态物性 II ：电子结构
没钱呀，我也还是取的小一点吧 8楼2013-08-09 21:06:07
董大瓜
新虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 3194.1
- 红花: 6
- 帖子: 158
- 在线: 57.9小时
- 虫号: 2326901
- 注册: 2013-03-07
- 专业: 凝聚态物性 II ：电子结构
9楼2013-08-10 19:33:32
董大瓜
新虫 (小有名气)
- 应助: 6 (幼儿园)
- 金币: 3194.1
- 红花: 6
- 帖子: 158
- 在线: 57.9小时
- 虫号: 2326901
- 注册: 2013-03-07
- 专业: 凝聚态物性 II ：电子结构
10楼2013-08-11 22:05:08