| 查看: 2723 | 回复: 15 | ||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||
cuicui2397新虫 (小有名气)
|
[求助]
vasp阶段能测试求助，谢谢啊！急求。
|
|
|
看的侯博士的教程，里面有测试Si的截断能，但我自己去测试老是出现问题，求解啊！谢谢。 POSCAR文件： Si-Diamond 5.430 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 2 Direct 0.0 0.0 0.0 0.25 0.25 0.25 KPOINTS文件： Automatic generation 0 Monhkorst-Pack 9 9 9 0.0 0.0 0.0 用来确定ENCUT的脚本程序为run_ecut，其内容是： #!/bin/sh rm -rf WAVECAR for i in 150 200 250 300 350 400 do cat > INCAR <<! SYSTEM = Si-Diamond ENCUT = $i ISTART = 0; ICHARG = 2 ISMEAR = -5 PREC = Accurate ! echo "ENCUT = $i eV"; mpirun -np 2 vasp E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}'' echo $i $E >>comment done 同时选择了赝势文件POTCAR，把以上四个文件POSCAR、KPOINTS、run_ecut和POTCAR放入我的计算目录下，在终端输入sh.run_ecut命令，在生成的comment里是这样的： 150 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 200 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 250 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 300 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 350 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 400 grep "TOTEN" OUTCAR |tail-1|awk {printf%12.6f\n,} 教程中第二列显示的是能量值，请问这是什么问题啊？谢谢啊！ |
回复此楼» 猜你喜欢
求助调剂，跨调
已经有20人回复
-
291 求调剂
已经有40人回复
-
300求调剂
已经有5人回复
-
急需调剂
已经有3人回复
-
22408 312求调剂
已经有21人回复
-
恳请有学校收留
已经有4人回复
-
297，工科调剂?
已经有9人回复
-
291求调剂
已经有5人回复
-
085404 22408 309分求调剂
已经有10人回复
-
296求调剂
已经有13人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- vasp运行，第一离子步结束后自动停了
已经有7人回复
- vasp中，做截断能测试时，测试收敛的标准
已经有10人回复
- 深圳超算中心vasp编译及任务提交问题
已经有25人回复
- 关于VASP一些琐事的求助（K点，吸附优化等）
已经有12人回复
- vasp并行测试时出错。
已经有7人回复
- vasp并行出错，急于求助！！！
已经有13人回复
- vasp计算：优化能收敛，但同样的设置静态计算不收敛
已经有3人回复
- VASP计算时的几个细节设置问题，希望有经验的人过来指点一下！
已经有23人回复
- VASP 截断能测试 急求
已经有8人回复
- 【求助成功】vasp测试K点收敛性的脚本问题
已经有20人回复
- 【求助】VASP怎么做晶格参数的结构优化
已经有17人回复
- 【求助】求将vasp的输出文件转化成xyz坐标的脚本
已经有11人回复
- 【求助】VASP计算电子态密度
已经有34人回复
- 【求助】vasp做表面计算时真空层的厚度
已经有13人回复
- 【求助】VASP能对体系加电场计算吗？
已经有18人回复
- 【讨论】vasp结构优化之前的测试
已经有11人回复
cuicui2397
新虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 1611.8
- 红花: 2
- 帖子: 211
- 在线: 130.4小时
- 虫号: 2198831
- 注册: 2012-12-21
- 专业: 凝聚态物性 II ：电子结构
6楼2013-11-15 10:29:02
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
cuicui2397(fzx2008代发): 金币+3 2013-11-14 21:25:06
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:10
cuicui2397: 金币+5, ★★★★★最佳答案 2013-11-23 10:10:40
感谢参与，应助指数 +1
cuicui2397(fzx2008代发): 金币+3 2013-11-14 21:25:06
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:10
cuicui2397: 金币+5, ★★★★★最佳答案 2013-11-23 10:10:40
|
那个，有可能是符号有问题，那个E='grep "TOTEN" OUTCAR |tail-1|awk '{printf"%12.6f \n",$5}'' 是`（ESC下面那个键，数字1左面那个键），而不是'(单引号)，你试试看 |
» 本帖已获得的红花（最新10朵）
2楼2013-11-14 20:30:12
cuicui2397
新虫 (小有名气)
- 应助: 5 (幼儿园)
- 金币: 1611.8
- 红花: 2
- 帖子: 211
- 在线: 130.4小时
- 虫号: 2198831
- 注册: 2012-12-21
- 专业: 凝聚态物性 II ：电子结构
送红花一朵 |
你好，首先谢谢你，我把符号改过之后，第二列出现了能量值。不过和教程上差别很大，你觉得是什么原因造成的啊？我的结果是这样的： 150 -10.744031 200 -10.795431 250 -10.822115 300 -10.836639 350 -10.840896 400 -10.842717 而教程中是这样写的： 150 -11.900655 200 -11.938864 250 -11.944599 300 -11.945248 350 -11.945503 400 -11.945622 总能变化在0.001eV左右就足够了，在这个例子中，因此，选择ENCUT=250eV。 我这样算出的没什么参考价值啊，求解，谢谢！ |
3楼2013-11-14 21:18:24
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:17
cuicui2397: 金币+12, ★★★很有帮助 2013-11-23 10:11:24
sunyang1988: 金币+1, 谢谢交流 2013-11-15 18:40:17
cuicui2397: 金币+12, ★★★很有帮助 2013-11-23 10:11:24
| 赝势是不是不一样，你有没有选正确的赝势？侯老师的那个教程，有没有说选什么赝势？ |
» 本帖已获得的红花（最新10朵）
4楼2013-11-14 23:06:04
